title: Bipartitioning of directed and mixed random graphs

abstract: We show that an intricate relation of cluster properties and optimal bipartitions, which takes place in undirected random graphs, extends to directed and mixed random graphs. In particular, the satisfability threshold coincides with the relative size of the giant OUT component reaching~{1/2}. Moreover, when counting undirected links as two directed ones, the partition cost, and cluster properties, as well as location of the replica symmetry breaking transition for these random graphs depend primarily on the total number of directed links and not on their specific distribution.


title: A computational method for estimating Burr XII parameters with complete and multiple censored data

abstract: Flexibility in shape and scale of Burr XII distribution can make close approximation of numerous well-known probability density functions. Due to these capabilities, the usages of Burr XII distribution are applied in risk analysis, lifetime data analysis and process capability estimation. In this paper the Cross-Entropy (CE) method is further developed in terms of Maximum Likelihood Estimation (MLE) to estimate the parameters of Burr XII distribution for the complete data or in the presence of multiple censoring. A simulation study is conducted to evaluate the performance of the MLE by means of CE method for different parameter settings and sample sizes. The results are compared to other existing methods in both uncensored and censored situations.


title: Orthogonal Time Frequency Space (OTFS) Modulation Based Radar System

abstract: Orthogonal time frequency space (OTFS) modulation was proposed to tackle the destructive Doppler effects in wireless communications, with potential applications to many other areas. In this paper, we investigate its application to radar systems, and propose a novel efficient OTFS-based matched filter algorithm for target range and velocity estimation. The proposed algorithm not only exhibits the inherent advantages due to multi-carrier modulation of the existing orthogonal frequency division multiplexing (OFDM-) based radar algorithms but also provides additional benefits to improve radar capability. Similar to OFDM, OTFS spreads the transmitted signal in the entire time--frequency resources to exploit the full diversity gains for radar processing. However, OTFS requires less cyclic prefix, and hence, shorter transmission duration than OFDM, allowing longer range radar and/or faster target tracking rate. Additionally, unlike OFDM, OTFS is inter-carrier interference-free, enabling larger Doppler frequency estimation. We demonstrate the performance of the proposed algorithm using numerical results under different system settings.


title: Codimension two cycles in Iwasawa theory and tensor product of Hida families

abstract: The purpose of this paper is to build on results in {\it{higher codimension Iwasawa theory}}. The setting of our results involves Galois representations arising as cyclotomic twist deformations associated to (i) the tensor product of two cuspidal Hida families $F$ and $G$, and (ii) the tensor product of three cuspidal Hida families $F$, $G$ and $H$. On the analytic side, we consider (i) a pair of $3$-variable Rankin-Selberg $p$-adic $L$-functions constructed by Hida and (ii) a balanced $4$-variable $p$-adic $L$-function (due to Hsieh and Yamana) and an unbalanced $4$-variable $p$-adic $L$-function (whose existence is currently conjectural). In each of these setups, when the two $p$-adic $L$-functions generate a height two ideal in the corresponding deformation ring, we use codimension two cycles of that ring to relate them to a pair of pseudo-null modules.


title: Open clusters in APOGEE and GALAH: Combining Gaia and ground-based spectroscopic surveys

abstract: Context: Open clusters are ideal laboratories to investigate a variety of astrophysical topics, from the properties of the Galactic disk to stellar evolutionary models. Knowing their metallicity and possibly detailed chemical abundances is therefore important. However, the number of systems with chemical abundances determined from high resolution spectroscopy is still small.\\ Aims: To increase the number of open clusters with radial velocities and chemical abundances determined from high resolution spectroscopy we used publicly available catalogues of surveys in combination with Gaia data.\\ Methods: Open cluster stars have been identified in the APOGEE and GALAH spectroscopic surveys by cross-matching their latest data releases with stars for which high-probability astrometric membership has been derived in many clusters on the basis of the Gaia second data release.\\ Results: Radial velocities have been determined for 131 and 14 clusters from APOGEE and GALAH data, respectively. This is the first radial velocity determination from high resolution spectra for 16 systems. Iron abundances have been obtained for 90 and 14 systems from APOGEE and GALAH samples, respectively. To our knowledge 66 of these clusters (57 in APOGEE and 9 in GALAH) do not have previous determinations in the literature. For 90 and 7 clusters in the APOGEE and GALAH samples, respectively, we have also determined average abundances for Na, Mg, Al, Si, Ca, Cr, Mn, and Ni.


title: Asymptotics of maximum likelihood estimation for stable law with continuous parameterization

abstract: Asymptotics of maximum likelihood estimation for $\alpha$-stable law are analytically investigated with a continuous parameterization. The consistency and asymptotic normality are shown on the interior of the whole parameter space. Although these asymptotics have been provided with Zolotarev's $(B)$ parameterization, there are several gaps between. Especially in the latter, the density, so that scores and their derivatives are discontinuous at $\alpha=1$ for $\beta\neq 0$ and usual asymptotics are impossible. This is considerable inconvenience for applications. By showing that these quantities are smooth in the continuous form, we fill gaps between and provide a convenient theory. We numerically approximate the Fisher information matrix around the Cauchy law $(\alpha,\beta)=(1,0)$. The results exhibit continuity at $\alpha=1,\,\beta\neq 0$ and this secures the accuracy of our calculations.


title: Energy loss of heavy and light quarks in holographic magnetized background

abstract: We systematically study holographic effects on the magnetic field dependence of the drag force, diffusion coefficient, jet quenching parameter of heavy quarks and the shooting string energy loss of light quarks in the RHIC and LHC energy regions by using the AdS/CFT correspondence in this paper. This study is motivated by the phenomena of strong magnetic field and jet quenching, which have been found in relativistic heavy ion collisions. The probe's direction of motion is perpendicular and parallel to the direction of magnetic field $B$. The effects of magnetic field on energy loss when moving perpendicular to the magnetic field direction are larger than moving parallel to the magnetic field direction, which implies that the magnetic field tends to suppress more quarks and jets when moving in the transverse direction than in the parallel direction. It is found that the diffusion coefficient decreases with the magnetic field in the transverse direction, but increases with the magnetic field in the parallel direction, which indicates that the quark may diffuse farther when moving parallel to the magnetic field direction. We also find that the magnetic field will enhance the energy loss of the light quarks when moving in the transverse direction than in the parallel direction.


title: Ripple patterns in in-plane velocities of OB stars from LAMOST and Gaia

abstract: With about 12 000 OB type stars selected from the LAMOST and Gaia survey, we study their 3 dimensional velocity distribution over the range of galactocentric radius from 6 to 15 kpc in the Galactic disk plane. A clear ripple pattern in the radial velocity ($V_R$) map is shown. The median $V_R$ reaches $-8$ km s$^{-1}$ at $R\sim9$ kpc, then increases to $\sim0$ km s$^{-1}$ at $R\sim12$ kpc, and later declines to below $-10$ km s$^{-1}$ beyond $R\sim13$ kpc. The median azimuthal velocity ($V_\phi$) map shows a similar pattern but has roughly $1/4$ phase difference with the radial velocity. Although the ripple of negative $V_R$ at $\sim9$ kpc extends to about 40$^\circ$ in the azimuth angle, it does not align with either the Local or the Perseus spiral arm. Moreover, the farther ripple beyond 13 kpc does not match the Outer spiral arm either. This indicates that the non-axisymmetric kinematic features are not induced by perturbations of known spiral structures. The central rotating bar can not lead to such patterns in the outer disk either. External perturbation of a dwarf galaxy or a dark matter sub-halo can induce such patterns but requires more evidence from both observations and simulations. The $V_\phi$ map in the $Z$--$V_Z$ plane of the OB stars is also investigated. Despite asymmetry to some degree, no spiral pattern is found. This is reasonable since most of the OB stars have ages much younger than 100 Myrs, which is smaller than one orbital period around the Galactic center.


title: Strain-Induced Room-Temperature Ferromagnetic Semiconductors with Large Anomalous Hall Conductivity in Two-Dimensional Cr2Ge2Se6

abstract: By density functional theory calculations, we predict a stable two-dimensional (2D) ferromagnetic semiconductor Cr$_2$Ge$_2$Se$_6$, where the Curie temperature $T$$_c$ can be dramatically enhanced beyond room temperature by applying a few percent strain. In addition, the anomalous Hall conductivity in 2D Cr$_2$Ge$_2$Se$_6$ and Cr$_2$Ge$_2$Te$_6$ is predicted to be comparable to that in ferromagnetic metals of Fe and Ni, and is an order of magnitude larger than that in diluted magnetic semiconductor Ga(Mn,As). Based on superexchange interactions, the enhanced $T$$_c$ in 2D Cr$_2$Ge$_2$Se$_6$ by strain can be understood by the decreased energy difference between 3$d$ orbitals of Cr and 4$p$ orbitals of Se. Our finding highlights the microscopic mechanism to obtain the room temperature ferromagnetic semiconductors by strain.


title: HetMEC: Latency-optimal Task Assignment and Resource Allocation for Heterogeneous Mobile Edge Computing

abstract: Driven by great demands on low-latency services of the edge devices (EDs), mobile edge computing (MEC) has been proposed to enable the computing capacities at the edge of the radio access network. However, conventional MEC servers suffer disadvantages such as limited computing capacity, preventing the computation-intensive tasks to be processed in time. To relief this issue, we propose the heterogeneous MEC (HetMEC) where the data that cannot be timely processed at the edge are allowed be offloaded to the upper-layer MEC servers, and finally to the cloud center (CC) with more powerful computing capacity. We design the latency minimization algorithm by jointly coordinating the task assignment, computing and transmission resources among the EDs, multi-layer MEC servers, and the CC. Simulation results indicate that our proposed algorithm can achieve a lower latency and higher processing rate than the conventional MEC scheme.


title: Energy-Efficient Resource Allocation for Secure UAV Communication Systems

abstract: In this paper, we study the resource allocation and trajectory design for energy-efficient secure unmanned aerial vehicle (UAV) communication systems where a UAV base station serves multiple legitimate ground users in the existence of a potential eavesdropper. We aim to maximize the energy efficiency of the UAV by jointly optimizing its transmit power, user scheduling, trajectory, and velocity. The design is formulated as a non-convex optimization problem taking into account the maximum tolerable signal-to-noise ratio (SNR) leakage, the minimum data rate requirement of each user, and the location uncertainty of the eavesdropper. An iterative algorithm is proposed to obtain an efficient suboptimal solution. Simulation results demonstrate that the proposed algorithm can achieve a significant improvement of the system energy efficiency while satisfying communication security constraint, compared to some simple scheme adopting straight flight trajectory with a constant speed.


title: High-dimensional statistical arbitrage with factor models and stochastic control

abstract: The present paper provides a study of high-dimensional statistical arbitrage that combines factor models with the tools from stochastic control, obtaining closed-form optimal strategies which are both interpretable and computationally implementable in a high-dimensional setting. Our setup is based on a general statistically-constructed factor model with mean-reverting residuals, in which we show how to construct analytically market-neutral portfolios and we analyze the problem of investing optimally in continuous time and finite horizon under exponential and mean-variance utilities. We also extend our model to incorporate constraints on the investor's portfolio like dollar-neutrality and market frictions in the form of temporary quadratic transaction costs, provide extensive Monte Carlo simulations of the previous strategies with 100 assets, and describe further possible extensions of our work.


title: Absence of dissipationless transport in clean 2D superconductors

abstract: Dissipationless charge transport is one of the defining properties of superconductors (SC). The interplay between dimensionality and disorder in determining the onset of dissipation in SCs remains an open theoretical and experimental problem. In this work, we present measurements of the dissipation phase diagrams of SCs in the two dimensional (2D) limit, layer by layer, down to a monolayer in the presence of temperature (T), magnetic field (B), and current (I) in 2H-NbSe2. Our results show that the phase-diagram strongly depends on the SC thickness even in the 2D limit. At four layers we can define a finite region in the I-B phase diagram where dissipationless transport exists at T=0. At even smaller thicknesses, this region shrinks in area. In a monolayer, we find that the region of dissipationless transport shrinks towards a single point, defined by T=B=I=0. In applied field, we show that time-dependent-Ginzburg-Landau (TDGL) simulations that describe dissipation by vortex motion, qualitatively reproduce our experimental I-B phase diagram. Last, we show that by using non-local transport and TDGL calculations that we can engineer charge flow and create phase boundaries between dissipative and dissipationless transport regions in a single sample, demonstrating control over non-equilibrium states of matter.


title: Q-learning with UCB Exploration is Sample Efficient for Infinite-Horizon MDP

abstract: A fundamental question in reinforcement learning is whether model-free algorithms are sample efficient. Recently, Jin et al. \cite{jin2018q} proposed a Q-learning algorithm with UCB exploration policy, and proved it has nearly optimal regret bound for finite-horizon episodic MDP. In this paper, we adapt Q-learning with UCB-exploration bonus to infinite-horizon MDP with discounted rewards \emph{without} accessing a generative model. We show that the \textit{sample complexity of exploration} of our algorithm is bounded by $\tilde{O}({\frac{SA}{\epsilon^2(1-\gamma)^7}})$. This improves the previously best known result of $\tilde{O}({\frac{SA}{\epsilon^4(1-\gamma)^8}})$ in this setting achieved by delayed Q-learning \cite{strehl2006pac}, and matches the lower bound in terms of $\epsilon$ as well as $S$ and $A$ except for logarithmic factors.


title: Robust Dynamic Resource Allocation via Probabilistic Task Pruning in Heterogeneous Computing Systems

abstract: In heterogeneous distributed computing (HC) systems, diversity can exist in both computational resources and arriving tasks. In an inconsistently heterogeneous computing system, task types have different execution times on heterogeneous machines. A method is required to map arriving tasks to machines based on machine availability and performance, maximizing the number of tasks meeting deadlines (defined as robustness). For tasks with hard deadlines (eg those in live video streaming), tasks that miss their deadlines are dropped. The problem investigated in this research is maximizing the robustness of an oversubscribed HC system. A way to maximize this robustness is to prune (ie defer or drop) tasks with low probability of meeting their deadlines to increase the probability of other tasks meeting their deadlines. In this paper, we first provide a mathematical model to estimate a task's probability of meeting its deadline in the presence of task dropping. We then investigate methods for engaging probabilistic dropping and we find thresholds for dropping and deferring. Next, we develop a pruning-aware mapping heuristic and extend it to engender fairness across various task types. We show the cost benefit of using probabilistic pruning in an HC system. Simulation results, harnessing a selection of mapping heuristics, show efficacy of the pruning mechanism in improving robustness (on average by 25%) and cost in an oversubscribed HC system by up to 40%.


title: Trigonometric Parallaxes of Star Forming Regions Beyond the Tangent Point of the Sagittarius Spiral Arm

abstract: As part of the BeSSeL Survey, we report trigonometric parallaxes and proper motions of molecular maser sources associated with 13 distant high mass star forming regions in the Sagittarius spiral arm of the Milky Way. In particular, we obtain improved parallax distance estimates for three well studied regions: 1.9 +0.1/-0.1 kpc for M17, 5.3 +1.3/-0.9 kpc for W51, and 7.9 +0.9/-0.7 kpc for GAL 045.5+00.0. Peculiar motions for all but one source are less than 20 km/s. We fit a log-periodic spiral to the locations and estimate an average pitch angle of 7.2+-1.9 deg. We find that the section of the arm beyond the tangent point in the first quadrant of the Milky Way appears 15 pc below the IAU-defined Galactic plane.


title: On Symmetric Losses for Learning from Corrupted Labels

abstract: This paper aims to provide a better understanding of a symmetric loss. First, we emphasize that using a symmetric loss is advantageous in the balanced error rate (BER) minimization and area under the receiver operating characteristic curve (AUC) maximization from corrupted labels. Second, we prove general theoretical properties of symmetric losses, including a classification-calibration condition, excess risk bound, conditional risk minimizer, and AUC-consistency condition. Third, since all nonnegative symmetric losses are non-convex, we propose a convex barrier hinge loss that benefits significantly from the symmetric condition, although it is not symmetric everywhere. Finally, we conduct experiments to validate the relevance of the symmetric condition.


title: Asynchronous Early Output Block Carry Lookahead Adder with Improved Quality of Results

abstract: A new asynchronous early output block carry lookahead adder (BCLA) incorporating redundant carries is proposed. Compared to the best of existing semi-custom asynchronous carry lookahead adders (CLAs) employing delay-insensitive data encoding and following a 4-phase handshaking, the proposed BCLA with redundant carries achieves 13% reduction in forward latency and 14.8% reduction in cycle time compared to the best of the existing CLAs featuring redundant carries with no area or power penalty. A hybrid variant involving a ripple carry adder (RCA) in the least significant stages i.e. BCLA-RCA is also considered that achieves a further 4% reduction in the forward latency and a 2.4% reduction in the cycle time compared to the proposed BCLA featuring redundant carries without area or power penalties.


title: Majority and Minority Voted Redundancy for Safety-Critical Applications

abstract: A new majority and minority voted redundancy (MMR) scheme is proposed that can provide the same degree of fault tolerance as N-modular redundancy (NMR) but with fewer function units and a less sophisticated voting logic. Example NMR and MMR circuits were implemented using a 32/28nm CMOS process and compared. The results show that MMR circuits dissipate less power, occupy less area, and encounter less critical path delay than the corresponding NMR circuits while providing the same degree of fault tolerance. Hence the MMR is a promising alternative to the NMR to efficiently implement high levels of redundancy in safety-critical applications.


title: A Structured Approach to the Analysis of Remote Sensing Images

abstract: The number of studies for the analysis of remote sensing images has been growing exponentially in the last decades. Many studies, however, only report results---in the form of certain performance metrics---by a few selected algorithms on a training and testing sample. While this often yields valuable insights, it tells little about some important aspects. For example, one might be interested in understanding the nature of a study by the interaction of algorithm, features, and the sample as these collectively contribute to the outcome; among these three, which would be a more productive direction in improving a study; how to assess the sample quality or the value of a set of features etc. With a focus on land-use classification, we advocate the use of a structured analysis. The output of a study is viewed as the result of the interplay among three input dimensions: feature, sample, and algorithm. Similarly, another dimension, the error, can be decomposed into error along each input dimension. Such a structural decomposition of the inputs or error could help better understand the nature of the problem and potentially suggest directions for improvement. We use the analysis of a remote sensing image at a study site in Guangzhou, China, to demonstrate how such a structured analysis could be carried out and what insights it generates. The structured analysis could be applied to a new study, or as a diagnosis to an existing one. We expect this will inform practice in the analysis of remote sensing images, and help advance the state-of-the-art of land-use classification.


title: Signal Shaping for Generalized Spatial Modulation and Generalized Quadrature Spatial Modulation

abstract: This paper investigates generic signal shaping methods for multiple-data-stream generalized spatial modulation (GenSM) and generalized quadrature spatial modulation (GenQSM) based on the maximizing the minimum Euclidean distance (MMED) criterion. Three cases with different channel state information at the transmitter (CSIT) are considered, including no CSIT, statistical CSIT and perfect CSIT. A unified optimization problem is formulated to find the optimal transmit vector set under size, power and sparsity constraints. We propose an optimization-based signal shaping (OBSS) approach by solving the formulated problem directly and a codebook-based signal shaping (CBSS) approach by finding sub-optimal solutions in discrete space. In the OBSS approach, we reformulate the original problem to optimize the signal constellations used for each transmit antenna combination (TAC). Both the size and entry of all signal constellations are optimized. Specifically, we suggest the use of a recursive design for size optimization. The entry optimization is formulated as a non-convex large-scale quadratically constrained quadratic programming (QCQP) problem and can be solved by existing optimization techniques with rather high complexity. To reduce the complexity, we propose the CBSS approach using a codebook generated by quadrature amplitude modulation (QAM) symbols and a low-complexity selection algorithm to choose the optimal transmit vector set. Simulation results show that the OBSS approach exhibits the optimal performance in comparison with existing benchmarks. However, the OBSS approach is impractical for large-size signal shaping and adaptive signal shaping with instantaneous CSIT due to the demand of high computational complexity. As a low-complexity approach, CBSS shows comparable performance and can be easily implemented in large-size systems.


title: Localizations for quiver Hecke algebras

abstract: We provide the localization procedure for monoidal categories by a real commuting family of braiders. For an element $w$ of the Weyl group, $\mathscr{C}_w$ is a subcategory of modules over quiver Hecke algebra which categorifies the quantum unipotent coordinate algebra $A_q[\mathfrak{n}(w)]$. We construct the localization $\widetilde{\mathscr{C}_w}$ of $\mathscr{C}_w$ by adding the inverses of simple modules which correspond to the frozen variables in the quantum cluster algebra $A_q[\mathfrak{n}(w)]$. The localization $\widetilde{\mathscr{C}_w}$ is left rigid and we expect that it is rigid.


title: Radiative heat transfer between metallic nanoparticle clusters in both near field and far field

abstract: Micro-nanoparticle systems have wide applications in thermal science and technology. In dense particulate system, the particle separation distance may be less than the characteristic thermal wavelength and near field effect will be significant and become a key factor to influence thermal radiation transfer in the system. In this study, radiative heat transfer (RHT) between two metallic nanoparticles clusters are explored using many-body radiative heat transfer theory implemented with the coupled electric and magnetic dipole (CEMD) approach, which effectively takes into account the contribution of magnetic polarization of metallic nanoparticles on heat exchange. As the focus, the effects of magnetic polarization and many-body interaction (MBI) on RHT were analyzed. The effects of fractal dimension and relative orientation of the clusters were also analyzed. Results show that the contribution of magnetically polarized eddy-current Joule dissipation dominates the RHT between Ag nanoparticle clusters. If only electric polarization (EP approach) is considered, the heat conductance will be underestimated as compared with the CEMD approach in both near field and far field regime. The effect of MBI on the RHT between Ag nanoparticle clusters is unobvious at room temperature, which is quite different from the SiC nanoparticle clusters. For the latter, MBI tends to suppress RHT significantly. The relative orientation has remarkable effect on radiative heat flux for clusters with lacy structure when the separation distance is in the near field. While for the separation distance in far field, both the relative orientation and the fractal dimension has a weak influence on radiative heat flux. This work will help the understanding of thermal transport in dense particulate system.


title: Fixup Initialization: Residual Learning Without Normalization

abstract: Normalization layers are a staple in state-of-the-art deep neural network architectures. They are widely believed to stabilize training, enable higher learning rate, accelerate convergence and improve generalization, though the reason for their effectiveness is still an active research topic. In this work, we challenge the commonly-held beliefs by showing that none of the perceived benefits is unique to normalization. Specifically, we propose fixed-update initialization (Fixup), an initialization motivated by solving the exploding and vanishing gradient problem at the beginning of training via properly rescaling a standard initialization. We find training residual networks with Fixup to be as stable as training with normalization -- even for networks with 10,000 layers. Furthermore, with proper regularization, Fixup enables residual networks without normalization to achieve state-of-the-art performance in image classification and machine translation.


title: Atomic mechanisms for the Si atom dynamics in graphene: chemical transformations at the edge and in the bulk

abstract: Recent advances in scanning transmission electron microscopy (STEM) allow to observe solid-state transformations and reactions in materials induced by thermal stimulus or electron beam on the atomic level. However, despite the rate at which large volumes of data can be generated (sometimes in the gigabyte to terabyte range per single experiment), approaches for the extraction of material-specific knowledge on the kinetics and thermodynamics of these processes are still lacking. One of the critical issues lies in being able to map the evolution of various atomic structures and determine the associated transition probabilities directly from raw experimental data characterized by high levels of noise and missing structural elements. Here, we demonstrate an approach based on the combination of multiple machine learning techniques to study the dynamic behavior of e-beam irradiated Si atoms in the bulk and at the edges of single-layer graphene in STEM experiments. First, a deep learning network is used to convert experimental STEM movies into coordinates of individual Si and carbon atoms. Then, a Gaussian mixture model is further used to establish the elementary atomic configurations of the Si atoms, defining the bonding geometries and chemical species and accounting for the discrete rotational symmetry of the host lattice. Finally, the frequencies and Markov transition probabilities between these states are determined. This analysis enables insight into the thermodynamics of defect populations and chemical reaction networks from the atomically resolved STEM data. Here, we observe a clear tendency for the formation of a 1D Si crystal along zigzag direction of graphene edges and for the Si impurity coupling to topological defects in bulk graphene.


title: Prediction of Silicate Glasses' Stiffness by High-Throughput Molecular Dynamics Simulations and Machine Learning

abstract: The development by machine learning of models predicting materials' properties usually requires the use of a large number of consistent data for training. However, quality experimental datasets are not always available or self-consistent. Here, as an alternative route, we combine machine learning with high-throughput molecular dynamics simulations to predict the Young's modulus of silicate glasses. We demonstrate that this combined approach offers excellent predictions over the entire compositional domain. By comparing the performance of select machine learning algorithms, we discuss the nature of the balance between accuracy, simplicity, and interpretability in machine learning.


title: Observation of angle-dependent mode conversion and mode hopping in 2D annular antidot lattice

abstract: We report spin-wave excitations in annular antidot lattice fabricated from 15 nm-thin Ni80Fe20 flm. The nanodots of 170 nm diameters are embedded in the 350 nm (diameter) antidot lattice to form the annular antidot lattice, which is arranged in a square lattice with edge-to-edge separation of 120 nm. A strong anisotropy in the spin-wave modes are observed with the change in orientation angle ({\phi}) of the in-plane bias magnetic field by using Time-resolved Magneto-optic Kerr microscope. A fattened fourfold rotational symmetry, mode hopping and mode conversion leading to mode quenching for three prominent spin-wave modes are observed in this lattice with the variation of the bias field orientation. Micromagnetic simulations enable us to successfully reproduce the measured evolution of frequencies with the orientation of bias magnetic field, as well as to identify the spatial profles of the modes. The magnetostatic field analysis, suggest the existence of magnetostatic coupling between the dot and antidot in annular antidot sample. Further local excitations of some selective spin-wave modes using numerical simulations showed the anisotropic spin-wave propagation through the lattice.


title: Observation of band narrowing and mode conversion in two-dimensional binary magnonic crystal

abstract: We introduce a new type of binary magnonic crystal, where Ni$_{80}$Fe$_{20}$ nanodots of two different sizes are diagonally connected forming a unit and those units are arranged in a square lattice. The magnetization dynamics of the sample is measured by using time-resolved magneto-optical Kerr effect microscope with varying magnitude and in-plane orientation ($\phi$) of the bias magnetic field. Interestingly, at $\phi=0^{\circ}$, the spin-wave mode profiles show frequency selective spatial localization of spin-wave power within the array. With the variation of $\phi$ in the range $0^{\circ}<\phi\leq 45^{\circ}$, we observe band narrowing due to localized to extended spin-wave mode conversion. Upon further increase of $\phi$, the spin-wave modes slowly lose the extended nature and become fully localized again at 90$^{\circ}$. We have extensively demonstrated the role of magnetostatic stray field distribution on the rotational symmetries obtained for the spin-wave modes. From micromagnetic simulations, we find that the dipole-exchange coupling between the nano-dots leads to remarkable modifications of the spin-wave mode profiles when compared with arrays of individual small and large dots. Numerically, we also show that the physical connection between the nano-dots provides more control points over the spin-wave propagation in the lattice at different orientations of bias magnetic field. This new type of binary magnonic crystal may find potential applications in magnonic devices such as spin-wave waveguide, filter, coupler, and other on-chip microwave communication devices.


title: Value Propagation for Decentralized Networked Deep Multi-agent Reinforcement Learning

abstract: We consider the networked multi-agent reinforcement learning (MARL) problem in a fully decentralized setting, where agents learn to coordinate to achieve the joint success. This problem is widely encountered in many areas including traffic control, distributed control, and smart grids. We assume that the reward function for each agent can be different and observed only locally by the agent itself. Furthermore, each agent is located at a node of a communication network and can exchanges information only with its neighbors. Using softmax temporal consistency and a decentralized optimization method, we obtain a principled and data-efficient iterative algorithm. In the first step of each iteration, an agent computes its local policy and value gradients and then updates only policy parameters. In the second step, the agent propagates to its neighbors the messages based on its value function and then updates its own value function. Hence we name the algorithm value propagation. We prove a non-asymptotic convergence rate 1/T with the nonlinear function approximation. To the best of our knowledge, it is the first MARL algorithm with convergence guarantee in the control, off-policy and non-linear function approximation setting. We empirically demonstrate the effectiveness of our approach in experiments.


title: Out-of-Time-Order Correlators in One-Dimensional XY model

abstract: Out-of-time-order correlators (OTOC) are considered to be a promising tool to characterize chaos in quantum systems. In this paper we study OTOC in XY model. With the presence of anisotropic parameter $\gamma$ and external magnetic field $\lambda$ in the Hamiltonian, we mainly focus on their influences on OTOC in thermodynamical limit. We find that the butterfly speed $v_B$ is dependent of these two parameters, and the recent conjectured universal form which characterizes the wavefront of chaos spreading are proved to be positive with varying $v_B$ in different phases of XY model. Moreover, we also study the behaviors of OTOC with fixed location, and we find that the early-time part fully agrees with the results derived from Hausdorff-Baker-Campbell expansion. The long-time part is studied either, while in the local case $C(t)$ decay as power law $t^{-1}$, $|F(t)|$ with nonlocal operators show quite interesting and nontrivial power law decay corresponding to different choices of operators and models. At last, we observe temperature dependence for OTOC with local operators at ($\gamma=0, \lambda=1$), and divergent behavior with low temperature for nonlocal operator case at late time.


title: On Spectral Cantor-Moran measures and a variant of Bourgain's sum of sine problem

abstract: In this paper, we show that if we have a sequence of Hadamard triples $\{(N_n,B_n,L_n)\}$ with $B_n\subset \{0,1,..,N_n-1\}$ for $n=1,2,...$, except an extreme case, then the associated Cantor-Moran measure $$ \begin{aligned} \mu = \mu(N_n,B_n) =& \delta_{\frac{1}{N_1}B_1}\ast\delta_{\frac{1}{N_1N_2}B_2}\ast \delta_{\frac{1}{N_1N_2N_3}B_3}\ast...\\ =& \mu_n\ast\mu_{>n} \end{aligned} $$ with support inside $[0,1]$ always admits an exponential orthonormal basis $E(\Lambda) = \{e^{2\pi i \lambda x}:\lambda\in\Lambda\}$ for $L^2(\mu)$, where $\Lambda$ is obtained from suitably modifying $L_n$. Here, $\mu_n$ is the convolution of the first $n$ Dirac measures and $\mu_{>n}$ denotes the tail-term. We show that the completeness of $E(\Lambda)$ in general depends on the ``equi-positivity" of the sequence of the pull-backed tail of the Cantor-Moran measure $\nu_{>n}(\cdot) = \mu_{>n}((N_1...N_n)^{-1}(\cdot))$. Such equi-positivity can be analyzed by the integral periodic zero set of the weak limit of $\{\nu_{>n}\}$. This result offers a new conceptual understanding of the completeness of exponential functions and it improves significantly many partial results studied by recent research, whose focus has been specifically on $\#B_n\le 4$. Using the Bourgain's example that a sum of sine can be asymptotically small, we shows that, in the extreme case, there exists some Cantor-Moran measure such that the equi-positive condition fails and the Fourier transform of the associated $\nu_{>n}$ uniformly converges on some unbounded set.


title: Connecting implicit and explicit large eddy simulations of two-dimensional turbulence through machine learning

abstract: In this article, we utilize machine learning to dynamically determine if a point on the computational grid requires implicit numerical dissipation for large eddy simulation (LES). The decision making process is learnt through \emph{a priori} training on quantities derived from direct numerical simulation (DNS) data. In particular, we compute eddy-viscosities obtained through the coarse graining of DNS quantities and utilize their distribution to categorize areas that require dissipation. If our learning determines that closure is necessary, an upwinded scheme is utilized for computing the non-linear Jacobian. In contrast, if it is determined that closure is unnecessary, a symmetric and second-order accurate energy and enstrophy preserving Arakawa scheme is utilized instead. This results in a closure framework that precludes the specification of any model-form for the small scale contributions of turbulence but deploys an appropriate numerical dissipation from explicit closure driven hypotheses. This methodology is deployed for the Kraichnan turbulence test-case and assessed through various statistical quantities such as angle-averaged kinetic energy spectra and vorticity structure functions. Our framework thus establishes a direct link between the use of explicit LES ideologies for closure and numerical scheme-based modeling of turbulence leading to improved statistical fidelity of \emph{a posteriori} simulations.


title: Reward Shaping via Meta-Learning

abstract: Reward shaping is one of the most effective methods to tackle the crucial yet challenging problem of credit assignment in Reinforcement Learning (RL). However, designing shaping functions usually requires much expert knowledge and hand-engineering, and the difficulties are further exacerbated given multiple similar tasks to solve. In this paper, we consider reward shaping on a distribution of tasks, and propose a general meta-learning framework to automatically learn the efficient reward shaping on newly sampled tasks, assuming only shared state space but not necessarily action space. We first derive the theoretically optimal reward shaping in terms of credit assignment in model-free RL. We then propose a value-based meta-learning algorithm to extract an effective prior over the optimal reward shaping. The prior can be applied directly to new tasks, or provably adapted to the task-posterior while solving the task within few gradient updates. We demonstrate the effectiveness of our shaping through significantly improved learning efficiency and interpretable visualizations across various settings, including notably a successful transfer from DQN to DDPG.


title: Towards a sound massive cosmology

abstract: It is known that de Rham-Gabadadze-Tolley (dRGT) massive gravity does not permit a homogeneous and isotropic universe with flat or spherical spatial metrics. We demonstrate that a singular reference metric solves this problem in an economic and straightforward way. In the dRGT massive gravity with a singular reference metric, there are sound homogeneous and isotropic cosmological solutions. We investigate cosmologies with the static and dynamical singular reference metrics, respectively. The term like dark energy appears naturally and the universe accelerates itself in some late time evolution. The term simulating dark matter also naturally emerges. We make a preliminary constraint on the parameters in the dRGT massive gravity in frame of the present cosmological model by using the data of supernovae, cosmic microwave back ground radiations, and baryonic acoustic oscillations.


title: Solution of an Open Problem about Two Families of Orthogonal Polynomials

abstract: An open problem about two new families of orthogonal polynomials was posed by Alhaidari. Here we will identify one of them as Wilson polynomials. The other family seems to be new but we show that they are discrete orthogonal polynomials on a bounded countable set with one accumulation point at 0 and we give some asymptotics as the degree tends to infinity.


title: Distributed Nash Equilibrium Seeking for Games in Systems with Bounded Control Inputs

abstract: Noticing that physical limitations are ubiquitous in practical engineering systems, this paper considers Nash equilibrium seeking for games in systems where the control inputs are bounded. More specifically, first-order integrator-type systems with bounded control inputs are firstly considered and two saturated control strategies are designed to seek for the Nash equilibrium of the game. Then, we further consider the Nash equilibrium seeking problem for games in second-order integrator-type systems. As this problem has rarely been investigated, we firstly propose a centralized seeking strategy without considering the boundedness of the control inputs, followed by a distributed counterpart. By further adapting a saturation function into the Nash equilibrium seeking strategy, a new seeking strategy is then designed for the considered second-order systems with bounded controls. In the proposed distributed strategies, consensus protocols are included for information sharing and the saturation functions are utilized to construct bounded control inputs. The convergence results are established through conducting Lyapunov stability analysis. Lastly, by considering the connectivity control of mobile sensor networks, the proposed methods are numerically verified.


title: Predicting Tomorrow's Headline using Today's Twitter Deliberations

abstract: Predicting the popularity of news article is a challenging task. Existing literature mostly focused on article contents and polarity to predict popularity. However, existing research has not considered the users' preference towards a particular article. Understanding users' preference is an important aspect for predicting the popularity of news articles. Hence, we consider the social media data, from the Twitter platform, to address this research gap. In our proposed model, we have considered the users' involvement as well as the users' reaction towards an article to predict the popularity of the article. In short, we are predicting tomorrow's headline by probing today's Twitter discussion. We have considered 300 political news article from the New York Post, and our proposed approach has outperformed other baseline models.


title: Augment your batch: better training with larger batches

abstract: Large-batch SGD is important for scaling training of deep neural networks. However, without fine-tuning hyperparameter schedules, the generalization of the model may be hampered. We propose to use batch augmentation: replicating instances of samples within the same batch with different data augmentations. Batch augmentation acts as a regularizer and an accelerator, increasing both generalization and performance scaling. We analyze the effect of batch augmentation on gradient variance and show that it empirically improves convergence for a wide variety of deep neural networks and datasets. Our results show that batch augmentation reduces the number of necessary SGD updates to achieve the same accuracy as the state-of-the-art. Overall, this simple yet effective method enables faster training and better generalization by allowing more computational resources to be used concurrently.


title: The homogenization method for topology optimization of structures: old and new

abstract: These are the lecture notes of a short course on the homogenization method for topology optimization of structures, given by Gr\'egoire Allaire, during the "GSIS International Summer School 2018" at Tohoku University (Sendai, Japan). The goal of this course is to review the necessary mathematical tools of homogenization theory and apply them to topology optimization of mechanical structures. The ultimate application, targeted in this course, is the topology optimization of structures built with lattice materials. Practical and numerical exercises are given, based on the finite element free software FreeFem++.


title: On the higher Riemann-Roch without denominators

abstract: We prove two refinements of the higher Riemann-Roch without denominators: a statement for regular closed immersions between arbitrary finite dimensional noetherian schemes, with no smoothness assumptions, and a statement for the relative cohomology of a proper morphism.


title: Breaking optomechanical cooling limit by two drive fields on a membrane-in-middle system

abstract: We present a theoretical scheme for ground state cooling of a mechanical resonator in a membrane-in-middle optomechanical system (OMS) driven by two red-detuned drive fields. The details of dynamical evolution of OMS are provided, and the effect of system conditions on cooling results are systematically studied. Most importantly, the setups with two drives are found to be capable of achieving better cooling results than the theoretical cooling limit with single cavity. Even an improvement by one order of thermal phonon number is possible with proper combination of the cavity damping rate and drive intensity.


title: Heterogeneity-aware Gradient Coding for Straggler Tolerance

abstract: Gradient descent algorithms are widely used in machine learning. In order to deal with huge volume of data, we consider the implementation of gradient descent algorithms in a distributed computing setting where multiple workers compute the gradient over some partial data and the master node aggregates their results to obtain the gradient over the whole data. However, its performance can be severely affected by straggler workers. Recently, some coding-based approaches are introduced to mitigate the straggler problem, but they are efficient only when the workers are homogeneous, i.e., having the same computation capabilities. In this paper, we consider that the workers are heterogeneous which are common in modern distributed systems. We propose a novel heterogeneity-aware gradient coding scheme which can not only tolerate a predetermined number of stragglers but also fully utilize the computation capabilities of heterogeneous workers. We show that this scheme is optimal when the computation capabilities of workers are estimated accurately. A variant of this scheme is further proposed to improve the performance when the estimations of the computation capabilities are not so accurate. We conduct our schemes for gradient descent based image classification on QingCloud clusters. Evaluation results show that our schemes can reduce the whole computation time by up to $3\times$ compared with a state-of-the-art coding scheme.


title: A simple INDIUM TIN OXIDE/glass DRA

abstract: A novel Dielectric Resonator Antenna, simply made of INDIUM TIN OXIDE coated glass slides placed on a microstrip transmission line, for communication applications is presented. Changes in the bandwidth and gain of the antenna are observed by modifying the dimensions of the INDIUM TIN OXIDE coated glass slides. Changes in gain, directivity and reflection coefficient are observed. A parametric study is conducted on the size of the DRA to understand the effect on bandwidth, reflection coefficient and gain.


title: Successive minima of line bundles

abstract: We introduce and study the successive minima of line bundles on proper algebraic varieties. The first (resp. last) minima are the width (resp. Seshadri constant) of the line bundle at very general points. The volume of the line bundle is equivalent to the product of the successive minima. For line bundles on toric varieties, the successive minima are equivalent to the (reciprocal of) usual successive minima of the difference of the moment polytope.


title: On the Universality of Invariant Networks

abstract: Constraining linear layers in neural networks to respect symmetry transformations from a group $G$ is a common design principle for invariant networks that has found many applications in machine learning. In this paper, we consider a fundamental question that has received little attention to date: Can these networks approximate any (continuous) invariant function? We tackle the rather general case where $G\leq S_n$ (an arbitrary subgroup of the symmetric group) that acts on $\mathbb{R}^n$ by permuting coordinates. This setting includes several recent popular invariant networks. We present two main results: First, $G$-invariant networks are universal if high-order tensors are allowed. Second, there are groups $G$ for which higher-order tensors are unavoidable for obtaining universality. $G$-invariant networks consisting of only first-order tensors are of special interest due to their practical value. We conclude the paper by proving a necessary condition for the universality of $G$-invariant networks that incorporate only first-order tensors.


title: Highly efficient cooling of mechanical resonator with square pulse drives

abstract: Ground state cooling of mechanical resonator is a way to generate macroscopic quantum states. Here we present a study of optomechanical cooling under the drive of square pulses without smooth profile. By illustrating the dynamical processes of cooling, we show how to choose the amplitudes and durations of square pulses, as well as the intervals between them, so that a mechanical resonator can be quickly cooled down to its ground state. Compared with the cooling under a continuous-wave drive field, the ground state cooling of a mechanical resonator can be performed more efficiently and flexibly by using square pulse drives. At certain times of such cooling process, the thermal phonon number under square pulse drives can become even lower than the theoretical limit for the cooling with a continuous-wave drive field of the same amplitude.


title: Stochastic Approximation of Smooth and Strongly Convex Functions: Beyond the $O(1/T)$ Convergence Rate

abstract: Stochastic approximation (SA) is a classical approach for stochastic convex optimization. Previous studies have demonstrated that the convergence rate of SA can be improved by introducing either smoothness or strong convexity condition. In this paper, we make use of smoothness and strong convexity simultaneously to boost the convergence rate. Let $\lambda$ be the modulus of strong convexity, $\kappa$ be the condition number, $F_*$ be the minimal risk, and $\alpha>1$ be some small constant. First, we demonstrate that, in expectation, an $O(1/[\lambda T^\alpha] + \kappa F_*/T)$ risk bound is attainable when $T = \Omega(\kappa^\alpha)$. Thus, when $F_*$ is small, the convergence rate could be faster than $O(1/[\lambda T])$ and approaches $O(1/[\lambda T^\alpha])$ in the ideal case. Second, to further benefit from small risk, we show that, in expectation, an $O(1/2^{T/\kappa}+F_*)$ risk bound is achievable. Thus, the excess risk reduces exponentially until reaching $O(F_*)$, and if $F_*=0$, we obtain a global linear convergence. Finally, we emphasize that our proof is constructive and each risk bound is equipped with an efficient stochastic algorithm attaining that bound.


title: Modelling Polymerization-Induced Self Assembly (PISA)

abstract: In this work we studied polymerization-induced self-assembly by means of computer simulations. Using this model, phase diagrams of the micelle states were constructed depending on the polymer concentration and the asymmetry of the composition for various reaction conditions. We found that if the reaction is ideal controlled radical polymerization (the initiation speed is much larger than the propagation speed and there are no side reactions such as termination or chain transfer), the phase diagram is no different from that obtained for pre-synthesized monodisperse diblock-copolymers with one insoluble block. Next, we studied two cases of slow initiation. We found that the phase diagram change dramatically: upon decreasing the initiation speed, the regions of spherical and cylindrical micelles shrink, while the region of vesicles/lamellae expands. This happens because at small initiation speed there is a significant amount of chains with a very short non-soluble block (and even without such block altogether), which do not participate in the formation of micelles. Therefore, decreasing the initiation speed essentially remaps the phase diagram coordinates by making the effective concentration lower (by decreasing the number of active chains) and the effective block length ratio higher (again, because the number of active chains gets higher, while the number of monomers remains the same).


title: Concrete Autoencoders for Differentiable Feature Selection and Reconstruction

abstract: We introduce the concrete autoencoder, an end-to-end differentiable method for global feature selection, which efficiently identifies a subset of the most informative features and simultaneously learns a neural network to reconstruct the input data from the selected features. Our method is unsupervised, and is based on using a concrete selector layer as the encoder and using a standard neural network as the decoder. During the training phase, the temperature of the concrete selector layer is gradually decreased, which encourages a user-specified number of discrete features to be learned. During test time, the selected features can be used with the decoder network to reconstruct the remaining input features. We evaluate concrete autoencoders on a variety of datasets, where they significantly outperform state-of-the-art methods for feature selection and data reconstruction. In particular, on a large-scale gene expression dataset, the concrete autoencoder selects a small subset of genes whose expression levels can be use to impute the expression levels of the remaining genes. In doing so, it improves on the current widely-used expert-curated L1000 landmark genes, potentially reducing measurement costs by 20%. The concrete autoencoder can be implemented by adding just a few lines of code to a standard autoencoder.


title: Reforming the Wishart characteristic function

abstract: The literature presents the characteristic function of the Wishart distribution on m times m matrices as an inverse power of the determinant of the Fourier variable, the exponent being the positive, real shape parameter. I demonstrate that only for two times two matrices, this expression is unambiguous -- in this case the complex range of the determinant excludes the negative real line. When m greater or equals 3 the range of the determinant contains closed lines around the origin, hence a single branch of the complex logarithm does not suffice to define the determinant's power. To resolve this issue, I give the correct analytic extension of the Laplace transform, by exploiting the Fourier-Laplace transform of a Wishart process.


title: Eva-CiM: A System-Level Performance and Energy Evaluation Framework for Computing-in-Memory Architectures

abstract: Computing-in-Memory (CiM) architectures aim to reduce costly data transfers by performing arithmetic and logic operations in memory and hence relieve the pressure due to the memory wall. However, determining whether a given workload can really benefit from CiM, which memory hierarchy and what device technology should be adopted by a CiM architecture requires in-depth study that is not only time consuming but also demands significant expertise in architectures and compilers. This paper presents an energy evaluation framework, Eva-CiM, for systems based on CiM architectures. Eva-CiM encompasses a multi-level (from device to architecture) comprehensive tool chain by leveraging existing modeling and simulation tools such as GEM5, McPAT [2] and DESTINY [3]. To support high-confidence prediction, rapid design space exploration and ease of use, Eva-CiM introduces several novel modeling/analysis approaches including models for capturing memory access and dependency-aware ISA traces, and for quantifying interactions between the host CPU and CiM modules. Eva-CiM can readily produce energy estimates of the entire system for a given program, a processor architecture, and the CiM array and technology specifications. Eva-CiM is validated by comparing with DESTINY [3] and [4], and enables findings including practical contributions from CiM-supported accesses, CiM-sensitive benchmarking as well as the pros and cons of increased memory size for CiM. Eva-CiM also enables exploration over different configurations and device technologies, showing 1.3-6.0X energy improvement for SRAM and 2.0-7.9X for FeFET-RAM, respectively.


title: Large Minors in Expanders

abstract: In this paper we study expander graphs and their minors. Specifically, we attempt to answer the following question: what is the largest function $f(n,\alpha,d)$, such that every $n$-vertex $\alpha$-expander with maximum vertex degree at most $d$ contains {\bf every} graph $H$ with at most $f(n,\alpha,d)$ edges and vertices as a minor? Our main result is that there is some universal constant $c$, such that $f(n,\alpha,d)\geq \frac{n}{c\log n}\cdot \left(\frac{\alpha}{d}\right )^c$. This bound achieves a tight dependence on $n$: it is well known that there are bounded-degree $n$-vertex expanders, that do not contain any grid with $\Omega(n/\log n)$ vertices and edges as a minor. The best previous result showed that $f(n,\alpha,d) \geq \Omega(n/\log^{\kappa}n)$, where $\kappa$ depends on both $\alpha$ and $d$. Additionally, we provide a randomized algorithm, that, given an $n$-vertex $\alpha$-expander with maximum vertex degree at most $d$, and another graph $H$ containing at most $\frac{n}{c\log n}\cdot \left(\frac{\alpha}{d}\right )^c$ vertices and edges, with high probability finds a model of $H$ in $G$, in time poly$(n)\cdot (d/\alpha)^{O\left( \log(d/\alpha) \right)}$. We note that similar but stronger results were independently obtained by Krivelevich and Nenadov: they show that $f(n,\alpha,d)=\Omega \left(\frac{n\alpha^2}{d^2\log n} \right)$, and provide an efficient algorithm, that, given an $n$-vertex $\alpha$-expander of maximum vertex degree at most $d$, and a graph $H$ with $O\left( \frac{n\alpha^2}{d^2\log n} \right)$ vertices and edges, finds a model of $H$ in $G$. Finally, we observe that expanders are the `most minor-rich' family of graphs in the following sense: for every $n$-vertex and $m$-edge graph $G$, there exists a graph $H$ with $O \left( \frac{n+m}{\log n} \right)$ vertices and edges, such that $H$ is not a minor of $G$.


title: Asymptotic development of an integral operator and boundedness of the criticality of potential centers

abstract: We study the asymptotic development at infinity of an integral operator. We use this development to give sufficient conditions in order to upper bound the number of critical periodic orbits that bifurcate from the outer boundary of the period function of planar potential centers. We apply the main results to two different families: the power-like potential family $\ddot x=x^p-x^q$, $p,q\in\mathbb{R}$, $p>q$; and the family of dehomogenized Loud's centers.


title: Automated Quality Control in Image Segmentation: Application to the UK Biobank Cardiac MR Imaging Study

abstract: Background: The trend towards large-scale studies including population imaging poses new challenges in terms of quality control (QC). This is a particular issue when automatic processing tools, e.g. image segmentation methods, are employed to derive quantitative measures or biomarkers for later analyses. Manual inspection and visual QC of each segmentation isn't feasible at large scale. However, it's important to be able to automatically detect when a segmentation method fails so as to avoid inclusion of wrong measurements into subsequent analyses which could lead to incorrect conclusions. Methods: To overcome this challenge, we explore an approach for predicting segmentation quality based on Reverse Classification Accuracy, which enables us to discriminate between successful and failed segmentations on a per-cases basis. We validate this approach on a new, large-scale manually-annotated set of 4,800 cardiac magnetic resonance scans. We then apply our method to a large cohort of 7,250 cardiac MRI on which we have performed manual QC. Results: We report results used for predicting segmentation quality metrics including Dice Similarity Coefficient (DSC) and surface-distance measures. As initial validation, we present data for 400 scans demonstrating 99% accuracy for classifying low and high quality segmentations using predicted DSC scores. As further validation we show high correlation between real and predicted scores and 95% classification accuracy on 4,800 scans for which manual segmentations were available. We mimic real-world application of the method on 7,250 cardiac MRI where we show good agreement between predicted quality metrics and manual visual QC scores. Conclusions: We show that RCA has the potential for accurate and fully automatic segmentation QC on a per-case basis in the context of large-scale population imaging as in the UK Biobank Imaging Study.


title: Magic numbers in polymer phase separation -- the importance of being rigid

abstract: Cells possess non-membrane-bound bodies, many of which are now understood as phase-separated condensates. One class of such condensates is composed of two polymer species, where each consists of repeated binding sites that interact in a one-to-one fashion with the binding sites of the other polymer. Previous biologically-motivated modeling of such a two-component system surprisingly revealed that phase separation is suppressed for certain combinations of numbers of binding sites. This phenomenon, dubbed the "magic-number effect", occurs if the two polymers can form fully-bonded small oligomers by virtue of the number of binding sites in one polymer being an integer multiple of the number of binding sites of the other. Here we use lattice-model simulations and analytical calculations to show that this magic-number effect can be greatly enhanced if one of the polymer species has a rigid shape that allows for multiple distinct bonding conformations. Moreover, if one species is rigid, the effect is robust over a much greater range of relative concentrations of the two species. Our findings advance our understanding of the fundamental physics of two-component polymer-based phase-separation and suggest implications for biological and synthetic systems.


title: Subsumption of Weakly Well-Designed SPARQL Patterns is Undecidable

abstract: Weakly well-designed SPARQL patterns is a recent generalisation of well-designed patterns, which preserve good computational properties but also capture almost all patterns that appear in practice. Subsumption is one of static analysis problems for SPARQL, along with equivalence and containment. In this paper we show that subsumption is undecidable for weakly well-designed patterns, which is in stark contrast to well-designed patterns, and to equivalence and containment.


title: An introductory guide to fluid models with anisotropic temperatures Part 1 -- CGL description and collisionless fluid hierarchy

abstract: We present a detailed guide to advanced collisionless fluid models that incorporate kinetic effects into the fluid framework, and that are much closer to the collisionless kinetic description than traditional magnetohydrodynamics. Such fluid models are directly applicable to modeling turbulent evolution of a vast array of astrophysical plasmas, such as the solar corona and the solar wind, the interstellar medium, as well as accretion disks and galaxy clusters. The text can be viewed as a detailed guide to Landau fluid models and it is divided into two parts. Part 1 is dedicated to fluid models that are obtained by closing the fluid hierarchy with simple (non Landau fluid) closures. Part 2 is dedicated to Landau fluid closures. Here in Part 1, we discuss the CGL fluid model in great detail, together with fluid models that contain dispersive effects introduced by the Hall term and by the finite Larmor radius (FLR) corrections to the pressure tensor. We consider dispersive effects introduced by the non-gyrotropic heat flux vectors. We investigate the parallel and oblique firehose instability, and show that the non-gyrotropic heat flux strongly influences the maximum growth rate of these instabilities. Furthermore, we discuss fluid models that contain evolution equations for the gyrotropic heat flux fluctuations and that are closed at the 4th-moment level by prescribing a specific form for the distribution function. For the bi-Maxwellian distribution, such a closure is known as the "normal" closure. We also discuss a fluid closure for the bi-kappa distribution. Finally, by considering one-dimensional Maxwellian fluid closures at higher-order moments, we show that such fluid models are always unstable. The last possible non Landau fluid closure is therefore the "normal" closure, and beyond the 4th-order moment, Landau fluid closures are required.


title: Nearly Optimal Sparse Polynomial Multiplication

abstract: In the sparse polynomial multiplication problem, one is asked to multiply two sparse polynomials f and g in time that is proportional to the size of the input plus the size of the output. The polynomials are given via lists of their coefficients F and G, respectively. Cole and Hariharan (STOC 02) have given a nearly optimal algorithm when the coefficients are positive, and Arnold and Roche (ISSAC 15) devised an algorithm running in time proportional to the "structural sparsity" of the product, i.e. the set supp(F)+supp(G). The latter algorithm is particularly efficient when there not "too many cancellations" of coefficients in the product. In this work we give a clean, nearly optimal algorithm for the sparse polynomial multiplication problem.


title: Observation of superconductivity in bilayer graphene/hexagonal boron nitride superlattices

abstract: A class of low-dimensional superconductivity (SC), such as most of "atomic-layer" SCs, has survived only under certain circumstances, implying a role of the substrate. Moreover, in some recent SC discoveries at heterogeneous interfaces, SC was buried in bulk solids and ex situ. Genuine atomic-layer SC is difficult to access. Here we report a novel route to atomic-layer SC in graphene superlattices. Our device comprises stacked non-twisted bilayer graphene (BLG) and hexagonal boron nitride (hBN), i.e., hBN/BLG/hBN Moir\'e superlattices. Upon in situ electrostatic doping, we observe an SC dome with a critical temperature up to $T_{\rm{BKT}} = 14 \rm{K}$, corresponding to the confinement of vortices. We believe that SC via doping Dirac materials is ubiquitous in condensed matter and that this study paves a way toward the design of a new SC family.


title: End-to-End Performance Analysis of Underwater Optical Wireless Relaying and Routing Techniques Under Location Uncertainty

abstract: On the contrary of low speed and high delay acoustic systems, underwater optical wireless communication (UOWC) can deliver a high speed and low latency service at the expense of short communication ranges. Therefore, multihop communication is of utmost importance to improve degree of connectivity and overall performance of underwater optical wireless networks (UOWNs). In this regard, this paper investigates relaying and routing techniques and provides their end-to-end (E2E) performance analysis under the location uncertainty. To achieve robust and reliable links, we first consider adaptive beamwidths and derive the divergence angles under the absence and presence of a pointing-acquisitioning-and-tracking (PAT) mechanism. Thereafter, important E2E performance metrics (e.g., data rate, bit error rate, transmission power, amplifier gain, etc.) are obtained for two potential relaying techniques; decode & forward (DF) and optical amplify & forward (AF). We develop centralized routing schemes for both relaying techniques to optimize E2E rate, bit error rate, and power consumption. Alternatively, a distributed routing protocol, namely Light Path Routing (LiPaR), is proposed by leveraging the range-beamwidth tradeoff of UOWCs. LiPaR is especially shown to be favorable when there is no PAT mechanism and available network information. In order to show the benefits of multihop communications, extensive simulations are conducted to compare different routing and relaying schemes under different network parameters and underwater environments.


title: Very fast helicity injection leading to critically stable state and large eruptive activity in solar active region NOAA 12673

abstract: Using the photospheric magnetic and coronal observations of Solar Dynamics Observatory, we studied the build-up and eruption of coronal non-potential magnetic structure in emerging active region (AR) 12673. The velocity field derived from tracked vector-magnetograms indicates persistent shear and converging motions of flux regions about the polarity inversion line (PIL). A major helicity injection occurs during rapid flux emergence consistent with the very fast flux emergence phase. While this helicity flux builds-up the sigmoid by September 4, the helicity injection by the continued shear and converging motions in the later evolution contributes to sigmoid sustenance and its core field twist as a manifestation of the flux rope which erupts after exceeding critical value of twist. Moreover, the total length of sheared PIL segments correlates with the non-neutralized current and maintains a higher value in both the polarity regions as a signature of eruptive capability of the AR according to the flux rope models. The modelled magnetic field qualitatively reproduces the sigmoidal structure capturing major features like twisted core flux as flux rope, and hook-shaped parts connecting at the middle of the PIL. Study of quasi-separatrix-layers reveals that the sheared arcade, enclosing the flux rope, is stressed to a critically stable state and its coronal height becomes doubled from September 4-6. While demonstrating the fast injection of helicity per unit flux as the crucial factor for severe space-weather events, this study explains the formation of the flux rope and recurrent eruptive nature of the AR by the critically stable state of sheared arcade early on September 6.


title: Application of the reflection functor to integrable systems on quiver varieties

abstract: The theory of representations of quivers and of their preprojective algebras are reviewed. In particular, moduli spaces of representations of these algebras, quiver varieties and reflection functor are described. The proof that the bijection between moduli spaces induced by the reflection functor is an isomorphism of symplectic affine varieties is presented. Hamiltonian systems with complete flows on the quiver varieties are introduced. An application of the reflection functor to these systems is described. The results on the case of the cyclic quiver are reviewed and a role of the reflection functor in this case is discussed. In particular, "spin" integrable generalisations of the Calogero-Moser systems and their application to the generalised KP hierarchies are considered.


title: An introductory guide to fluid models with anisotropic temperatures Part 2 -- Kinetic theory, Pad\'e approximants and Landau fluid closures

abstract: In Part 2 of our guide to collisionless fluid models, we concentrate on Landau fluid closures. These closures were pioneered by Hammett and Perkins and allow for the rigorous incorporation of collisionless Landau damping into a fluid framework. It is Landau damping that sharply separates traditional fluid models and collisionless kinetic theory, and is the main reason why the usual fluid models do not converge to the kinetic description, even in the long-wavelength low-frequency limit. We start with a brief introduction to kinetic theory, where we discuss in detail the plasma dispersion function $Z(\zeta)$, and the associated plasma response function $R(\zeta)=1+\zeta Z(\zeta)=-Z'(\zeta)/2$. We then consider a 1D (electrostatic) geometry and make a significant effort to map all possible Landau fluid closures that can be constructed at the 4th-order moment level. These closures for parallel moments have general validity from the largest astrophysical scales down to the Debye length, and we verify their validity by considering examples of the (proton and electron) Landau damping of the ion-acoustic mode, and the electron Landau damping of the Langmuir mode. We proceed by considering 1D closures at higher-order moments than the 4th-order, and as was concluded in Part 1, this is not possible without Landau fluid closures. We show that it is possible to reproduce linear Landau damping in the fluid framework to any desired precision, thus showing the convergence of the fluid and collisionless kinetic descriptions. We then consider a 3D (electromagnetic) geometry in the gyrotropic (long-wavelength low-frequency) limit and map all closures that are available at the 4th-order moment level. In the Appendix A, we provide comprehensive tables with Pad\'e approximants of $R(\zeta)$ up to the 8th-pole order, with many given in an analytic form.


title: On vertex adjacencies in the polytope of pyramidal tours with step-backs

abstract: We consider the traveling salesperson problem in a directed graph. The pyramidal tours with step-backs are a special class of Hamiltonian cycles for which the traveling salesperson problem is solved by dynamic programming in polynomial time. The polytope of pyramidal tours with step-backs $PSB (n)$ is defined as the convex hull of the characteristic vectors of all possible pyramidal tours with step-backs in a complete directed graph. The skeleton of $PSB (n)$ is the graph whose vertex set is the vertex set of $PSB (n)$ and the edge set is the set of geometric edges or one-dimensional faces of $PSB (n)$. The main result of the paper is a necessary and sufficient condition for vertex adjacencies in the skeleton of the polytope $PSB (n)$ that can be verified in polynomial time.


title: Many-body dynamical localisation of coupled quantum kicked rotors

abstract: The quantum motion of $N$ coupled kicked rotors is mapped to an interacting $N$-particle Anderson-Aubry-Andr$\'e$ tight-binding problem supporting many-body localised (MBL) phases. Interactions in configuration space are known to be insufficient for destroying Anderson localisation in a system in the MBL phase. The mapping we establish here predicts that a similar effect takes place in momentum space and determines the quantum dynamics of the coupled kicked rotors. Due to the boundedness of the Floquet quasi-energy spectrum there exists limitations on the interacting lattice models that can be mapped to quantum kicked rotors; in particular, no extensive observable can be mapped in the thermodynamic limit.


title: Hereditary classes of ordered binary structures

abstract: Balogh, Bollob\'{a}s and Morris (2006) have described a threshold phenomenon in the behavior of the profile of hereditary classes of ordered graphs. In this paper, we give an other look at their result based on the notion of monomorphic decomposition of a relational structure introduced in \cite{P-T-2013}. We prove that the class $\mathfrak S$ of ordered binary structures which do not have a finite monomorphic decomposition has a finite basis (a subset $\mathfrak A$ such that every member of $\mathfrak S$ embeds some member of $\mathfrak A$). In the case of ordered reflexive directed graphs, the basis has 1242 members and the profile of their ages grows at least as the Fibonacci function. From this result, we deduce that the following dichotomy property holds for every hereditary class $\mathfrak C$ of finite ordered binary structures of a given finite type. Either there is an integer $\ell$ such that every member of $\mathfrak C$ has a monomorphic decomposition into at most $\ell$ blocks and in this case the profile of $\mathfrak C$ is bounded by a polynomial of degree $\ell-1$ (and in fact is a polynomial), or $\mathfrak C$ contains the age of a structure which does not have a finite monomorphic decomposition, in which case the profile of $\mathfrak C$ is bounded below by the Fibonacci function.


title: Resultant Based Incremental Recovery of Camera Pose from Pairwise Matches

abstract: Incremental (online) structure from motion pipelines seek to recover the camera matrix associated with an image $I_n$ given $n-1$ images, $I_1,...,I_{n-1}$, whose camera matrices have already been recovered. In this paper, we introduce a novel solution to the six-point online algorithm to recover the exterior parameters associated with $I_n$. Our algorithm uses just six corresponding pairs of 2D points, extracted each from $I_n$ and from \textit{any} of the preceding $n-1$ images, allowing the recovery of the full six degrees of freedom of the $n$'th camera, and unlike common methods, does not require tracking feature points in three or more images. Our novel solution is based on constructing a Dixon resultant, yielding a solution method that is both efficient and accurate compared to existing solutions. We further use Bernstein's theorem to prove a tight bound on the number of complex solutions. Our experiments demonstrate the utility of our approach.


title: Joint models as latent Gaussian models - not reinventing the wheel

abstract: Joint models have received increasing attention during recent years with extensions into various directions; numerous hazard functions, different association structures, linear and non-linear longitudinal trajectories amongst others. Many of these resulted in new R packages and new formulations of the joint model. However, a joint model with a linear bivariate Gaussian association structure is still a latent Gaussian model (LGM) and thus can be implemented using most existing packages for LGM's. In this paper, we will show that these joint models can be implemented from a LGM viewpoint using the R-INLA package. As a particular example, we will focus on the joint model with a non-linear longitudinal trajectory, recently developed and termed the partially linear joint model. Instead of the usual spline approach, we argue for using a Bayesian smoothing spline framework for the joint model that is stable with respect to knot selection and hence less cumbersome for practitioners.


title: 6D Object Pose Estimation Based on 2D Bounding Box

abstract: In this paper, we present a simple but powerful method to tackle the problem of estimating the 6D pose of objects from a single RGB image. Our system trains a novel convolutional neural network to regress the unit quaternion, which represents the 3D rotation, from the partial image inside the bounding box returned by 2D detection systems. Then we propose an algorithm we call Bounding Box Equation to efficiently and accurately obtain the 3D translation, using 3D rotation and 2D bounding box. Considering that the quadratic sum of the quaternion's four elements equals to one, we add a normalization layer to keep the network's output on the unit sphere and put forward a special loss function for unit quaternion regression. We evaluate our method on the LineMod dataset and experiment shows that our approach outperforms base-line and some state of the art methods.


title: A Privacy Preserving Randomized Gossip Algorithm via Controlled Noise Insertion

abstract: In this work we present a randomized gossip algorithm for solving the average consensus problem while at the same time protecting the information about the initial private values stored at the nodes. We give iteration complexity bounds for the method and perform extensive numerical experiments.


title: The distribution of oxygen at the Ni81Fe19/Ta interface

abstract: The knowledge of how oxygen atoms are distributed at a magnetic-metal / oxide, or magnetic-metal / non-magnetic-metal interface, can be an useful tool to optimize device production. Multilayered Ni81Fe19 / Ta samples consisting of 15 bilayers of 2.5 nm each, grown onto glass substrates by magnetron sputtering from Ni81Fe19 and Ta targets, have been investigated. X-ray absorption near edge structure, extended X-Ray absorption fine structure, small angle X-ray diffraction, and simulations, were used to characterize the samples. Oxygen atoms incorporated onto Ni81Fe19 films during O2 exposition are mainly bonded to Fe atoms. This partial oxidation of the Ni81Fe19 surface works as a barrier to arriving Ta atoms, preventing intermixing at the Ni81Fe19 / Ta interface. The reduction of the Ni81Fe19 surface by the formation of TaO x is observed.


title: Instabilities of Thin Black Rings: Closing the Gap

abstract: We initiate the study of dynamical instabilities of higher-dimensional black holes using the blackfold approach, focusing on asymptotically flat boosted black strings and singly-spinning black rings in $D\ge5$. We derive novel analytic expressions for the growth rate of the Gregory-Laflamme instability for boosted black strings and its onset for arbitrary boost parameter. In the case of black rings, we study their stability properties in the region of parameter space that has so far remained inaccessible to numerical approaches. In particular, we show that very thin (ultraspinning) black rings exhibit a Gregory-Laflamme instability, giving strong evidence that black rings are unstable in the entire range of parameter space. For very thin rings, we show that the growth rate of the instability increases with increasing non-axisymmetric mode $m$ while for thicker rings, there is competition between the different modes. However, up to second order in the blackfold approximation, we do not observe an elastic instability, in particular for large modes $m\gg1$, where this approximation has higher accuracy. This suggests that the Gregory-Laflamme instability is the dominant instability for very thin black rings. Additionally, we find a long-lived mode that describes a wiggly time-dependent deformation of a black ring. We comment on disagreements between our results and corresponding ones obtained from a large $D$ analysis of black ring instabilities.


title: Negative optical torque on a microsphere in optical tweezers

abstract: We show that the optical force field in optical tweezers with elliptically polarized beams has the opposite handedness for a wide range of particle sizes and for the most common configurations. Our method is based on the direct observation of the particle equilibrium position under the effect of a transverse Stokes drag force, and its rotation around the optical axis by the mechanical effect of the optical torque. We find overall agreement with theory, with no fitting, provided that astigmatism, which is characterized separately, is included in the theoretical description. Our work opens the way for characterization of the trapping parameters, such as the microsphere complex refractive index and the astigmatism of the optical system, from measurements of the microsphere rotation angle.


title: Topological games of bounded selections

abstract: We present a new variation of the classical selection principles $\mathsf{S}_\mathrm{k}(\mathcal A, \mathcal B)$ ($k\in\mathbb N$) and $\mathsf{S}_\mathrm{fin}(\mathcal A, \mathcal B)$ that formally lies between these two properties. As in the case of the classical selection principles, we also obtain a new variation of topological games and discuss how new topological properties may emerge in the specific cases of covering and tightness.


title: The Green-Tao theorem for Piatetski-Shapiro primes

abstract: Let $m\geq 3$. Suppose that $$ 1-2^{-2^{m^24^m}}<\gamma<1. $$ Then the set $$ \{p\text{ prime}:\, p=[n^{\frac1\gamma}]\text{ for some }n\in{\mathbb N}\} $$ contains infinitely many non-trivial $m$-term arithmetic progressions.


title: Mirror symmetry for K3 surfaces

abstract: For certain K3 surfaces, there are two constructions of mirror symmetry that are very different. The first, known as BHK mirror symmetry, comes from the Landau-Ginzburg model for the K3 surface; the other, known as LPK3 mirror symmetry, is based on a lattice polarization of the K3 surface in the sense of Dolgachev's definition. There is a large class of K3 surfaces for which both versions of mirror symmetry apply. In this class we consider the K3 surfaces admitting a certain purely nonsymplectic automorphism of order 4, 8, or 12, and we complete the proof that these two formulations of mirror symmetry agree for this class of K3 surfaces.


title: Authentication and Authorization for Mobile IoT Devices using Bio-features: Recent Advances and Future Trends

abstract: Bio-features are fast becoming a key tool to authenticate the IoT devices; in this sense, the purpose of this investigation is to summaries the factors that hinder biometrics models' development and deployment on a large scale, including human physiological (e.g., face, eyes, fingerprints-palm, or electrocardiogram) and behavioral features (e.g., signature, voice, gait, or keystroke). The different machine learning and data mining methods used by authentication and authorization schemes for mobile IoT devices are provided. Threat models and countermeasures used by biometrics-based authentication schemes for mobile IoT devices are also presented. More specifically, We analyze the state of the art of the existing biometric-based authentication schemes for IoT devices. Based on the current taxonomy, We conclude our paper with different types of challenges for future research efforts in biometrics-based authentication schemes for IoT devices.


title: Discrete Sturm-Liouville problems: singularity of the $n$-th eigenvalue with application to Atkinson type

abstract: In this paper, we characterize singularity of the $n$-th eigenvalue of self-adjoint discrete Sturm-Liouville problems in any dimension. For a fixed Sturm-Liouville equation, we completely characterize singularity of the $n$-th eigenvalue. For a fixed boundary condition, unlike in the continuous case, the $n$-th eigenvalue exhibits jump phenomena and we describe the singularity under a non-degenerate assumption. Compared with the continuous case in [8, 12], the singular set here is involved heavily with coefficients of the Sturm-Liouville equations. This, along with arbitrariness of the dimension, causes difficulty when dividing areas in layers of the considered space such that the $n$-th eigenvalue has the same singularity in any given area. We study the singularity by partitioning and analyzing the local coordinate systems, and provide a Hermitian matrix which can determine the areas' division. To prove the asymptotic behavior of the $n$-th eigenvalue, we generalize the method developed in [21] to any dimension. Finally, by transforming the Sturm-Liouville problem of Atkinson type in any dimension to a discrete one, we can not only determine the number of eigenvalues, but also apply our approach above to obtain the complete characterization of singularity of the $n$-th eigenvalue for the Atkinson type.


title: On Peak Age of Information in Data Preprocessing enabled IoT Networks

abstract: Internet of Things (IoT) has been emerging as one of the use cases permeating our daily lives in 5th Generation wireless networks, where status update packages are usually required to be timely delivered for many IoT based intelligent applications. Enabling the collected raw data to be preprocessed before transmitted to the destination can provider users with better context-aware services and lighten the transmission burden. However, the effect from data preprocessing on the overall information freshness is an essential yet unrevealed issue. In this work we study the joint effects of data preprocessing and transmission procedures on information freshness measured by peak age of information (PAoI). Particularity, we formulate the considered multi-source preprocessing and transmission enabled IoT system as a tandem queue where a priority M/G/1 queue is followed by a G/G/1 queue. Then, we respectively derive the closed-form and an information theoretic approximation of the expectation of waiting time for the formulated processing queue and transmission queue, and further get the analytical expressions of the average PAoI for packages from different sources. Finally, the accuracy of our analysis is verified with simulation results.


title: On the Existence of Telescopers for Rational Functions in Three Variables

abstract: Zeilberger's method of creative telescoping is crucial for the computer-generated proofs of combinatorial and special-function identities. Telescopers are linear differential or ($q$-)recurrence operators computed by algorithms for creative telescoping. For a given class of inputs, when telescopers exist and how to construct telescopers efficiently if they exist are two fundamental problems related to creative telescoping. In this paper, we solve the existence problem of telescopers for rational functions in three variables including 18 cases. We reduce the existence problem from the trivariate case to the bivariate case and some related problems. The existence criteria given in this paper enable us to determine the termination of algorithms for creative telescoping with trivariate rational inputs.


title: Evolution of superclusters in the cosmic web

abstract: Aims. We investigate how properties of the ensemble of superclusters in the cosmic web evolve with time. Methods. We perform numerical simulations of the evolution of the cosmic web using the LambdaCDM model in box sizes L0 = 1024, 512, 256 Mpc/h. We find supercluster ensembles of models for four evolutionary stages, corresponding to the present epoch z = 0, and to redshifts z = 1, z = 3, and z = 10. We calculate fitness diameters of superclusters defined from volumes of superclusters divided to filling factors of over-density regions. Geometrical and fitness diameters of largest superclusters, and the number of superclusters as functions of the threshold density are used as percolation functions to describe geometrical properties of the ensemble of superclusters in the cosmic web. We calculate distributions of geometrical and fitness diameters and luminosities of superclusters, and follow time evolution of percolation functions and supercluster distributions. We compare percolation functions and supercluster distributions of models and samples of galaxies of the Sloan Digital Sky Survey (SDSS). Results. Our analysis shows that fitness diameters of superclusters have a minimum at certain threshold density. Fitness diameters around minima almost do not change with time in co-moving coordinates. Numbers of superclusters have maxima which are approximately constant for all evolutionary epochs. Geometrical diameters of superclusters decrease during the evolution of the cosmic web; luminosities of superclusters increase during the evolution. Conclusions. Our study suggests that evolutionary changes occur inside dynamical volumes of superclusters. The stability of fitness diameters and numbers of superclusters during the evolution is an important property of the cosmic web.


title: Updated analysis of recent results on electron and positron elastic scattering on proton

abstract: We discuss recent experimental results concerning the cross section ratio of positron over electron elastic scattering on protons, and compare with the predictions of a pre-existent calculation. The deviation from unity of this ratio, $i.e.$, a charge asymmetry different from zero, is the signature of contributions beyond the Born approximation. After reviewing the published results, we compare the elastic data to a calculation which includes the diagram corresponding to two-photon exchange. It turns out that all the data on the cross section ratio, in the limit of their precision, do not show evidence of enhanced two-photon contribution beyond the expected percent level. Our results confirm that experimental evidence for a large contribution of two-photon exchange is not yet found.


title: Acceleration of the NVT-flash calculation for multicomponent mixtures using deep neural network models

abstract: Phase equilibrium calculation, also known as flash calculation, has been extensively applied in petroleum engineering, not only as a standalone application for separation process but also an integral component of compositional reservoir simulation. It is of vital importance to accelerate flash calculation without much compromise in accuracy and reliability, turning it into an active research topic in the last two decades. In this study, we establish a deep neural network model to approximate the iterative NVT-flash calculation. A dynamic model designed for NVT flash problems is iteratively solved to produce data for training the neural network. In order to test the model's capacity to handle complex fluid mixtures, three real reservoir fluids are investigated, including one Bakken oil and two Eagle Ford oils. Compared to previous studies that follow the conventional flash framework in which stability testing precedes phase splitting calculation, we incorporate stability test and phase split calculation together and accomplish both two steps by a single deep learning model. The trained model is able to identify the single vapor, single liquid and vapor-liquid state under the subcritical region of the investigated fluids. A number of examples are presented to show the accuracy and efficiency of the proposed deep neural network. It is found that the trained model makes predictions at most 244 times faster than the iterative flash calculation under the given cases. Even though training a multi-level network model does take a large amount of time that is comparable to the computational time of flash calculations, the one-time offline training process gives the deep learning model great potential to speed up compositional reservoir simulation.


title: Dual Co-Matching Network for Multi-choice Reading Comprehension

abstract: Multi-choice reading comprehension is a challenging task that requires complex reasoning procedure. Given passage and question, a correct answer need to be selected from a set of candidate answers. In this paper, we propose \textbf{D}ual \textbf{C}o-\textbf{M}atching \textbf{N}etwork (\textbf{DCMN}) which model the relationship among passage, question and answer bidirectionally. Different from existing approaches which only calculate question-aware or option-aware passage representation, we calculate passage-aware question representation and passage-aware answer representation at the same time. To demonstrate the effectiveness of our model, we evaluate our model on a large-scale multiple choice machine reading comprehension dataset (i.e. RACE). Experimental result show that our proposed model achieves new state-of-the-art results.


title: LAXPC / AstroSat Study of ~ 1 and ~ 2 mHz Quasi-periodic Oscillations in the Be/X-ray Binary 4U 0115+63 During its 2015 Outburst

abstract: The Be X-ray Binary 4U 0115+63 was observed by Large Area X-ray Proportional Counter (LAXPC) instrument on AstroSat on 2015 October 24 during the peak of a giant Type II outburst. Prominent intensity oscillations at ~ 1 and ~ 2 mHz frequency were detected during the outburst. Nuclear Spectroscopic Telescope Array (NuSTAR) observations made during the same outburst also show mHz quasi periodic oscillations (QPOs). Details of the oscillations and their characteristics deduced from LAXPC/AstroSat and NuSTAR observations are reported in this paper. Analysis of the archival Rossi X-ray Timing Explorer (RXTE)/Proportional Counter Array (PCA) data during 2001-11 also show presence of mHz QPOs during some of the outbursts and details of these QPOs are also reported. Possible models to explain the origin of the mHz oscillations are examined. Similar QPOs, albeit at higher frequencies, have been reported from other neutron star and black hole sources and both may have a common origin. Current models to explain the instability in the inner accretion disk causing the intense oscillations are discussed.


title: Cutoff on Ramanujan complexes and classical groups

abstract: The total-variation cutoff phenomenon has been conjectured to hold for simple random walk on all transitive expanders. However, very little is actually known regarding this conjecture, and cutoff on sparse graphs in general. In this paper we establish total-variation cutoff for simple random walk on Ramanujan complexes of type $\widetilde{A}_{d}$ $(d\geq1)$. As a result, we obtain explicit generators for the finite classical groups $\mathrm{PGL}_{n}(\mathbb{F}_{q})$ for which the associated Cayley graphs exhibit total-variation cutoff.