title: A two-stage fourth order time-accurate discretization for Lax--Wendroff type flow solvers II. High order numerical boundary conditions

abstract: This paper serves to treat boundary conditions numerically with high order accuracy in order to match the two-stage fourth-order finite volume schemes for hyperbolic problems developed in [{\em J. Li and Z. Du, A two-stage fourth order time-accurate discretization {L}ax--{W}endroff type flow solvers, {I}. {H}yperbolic conservation laws, SIAM, J. Sci. Comput., 38 (2016), pp.~A3046--A3069}]. As such, it is significant when capturing small scale structures near physical boundaries. Different from previous contributions in literature, the current approach constructs a fourth order accurate approximation to boundary conditions by only using the Jacobian of the flux function (characteristic information) instead of its successive differentiation leading to tensors of high ranks in the inverse Lax-Wendroff method. Technically, data in several ghost cells are constructed with interpolation so that the interior scheme can be implemented over boundary cells, and theoretical boundary condition has to be modified properly at intermediate stages so as to make the two-stage scheme over boundary cells fully consistent with that over interior cells. This highlights the fact that {\em continuous boundary conditions only match continuous partial differential equations (PDEs), and they must be approximated in a consistent way (even though it could be exactly valued) when the PDEs are discretized.} Several numerical examples are provided to illustrate the performance of the current approach when dealing with general boundary conditions.


title: Reconciling Magma-Ocean Crystallization Models with the present-day Structure of the Earth's mantle

abstract: Terrestrial planets are thought to experience episode(s) of large-scale melting early in their history. Fractionation during magma-ocean freezing leads to unstable stratification within the related cumulate layers due to progressive iron enrichment upward, but the effects of incremental cumulate overturns during MO crystallization remain to be explored. Here, we use geodynamic models with a moving-boundary approach to study convection and mixing within the growing cumulate layer, and thereafter within the fully crystallized mantle. For fractional crystallization, cumulates are efficiently stirred due to subsequent incremental overturns, except for strongly iron-enriched late-stage cumulates, which persist as a stably stratified layer at the base of the mantle for billions of years. Less extreme crystallization scenarios can lead to somewhat more subtle stratification. In any case, the long-term preservation of at least a thin layer of extremely enriched cumulates with Fe#>0.4, as predicted by all our models, is inconsistent with seismic constraints. Based on scaling relationships, however, we infer that final-stage Fe-rich magma-ocean cumulates originally formed near the surface should have overturned as small diapirs, and hence undergone melting and reaction with the host rock during sinking. The resulting moderately iron-enriched metasomatized/hybrid rock assemblages should have accumulated at the base of the mantle, potentially fed an intermittent basal magma ocean, and be preserved through the present-day. Such moderately iron-enriched rock assemblages can reconcile the physical properties of the large low shear-wave velocity provinces in the presentday lower mantle. Thus, we reveal Hadean melting and rock-reaction processes by integrating magmaocean crystallization models with the seismic-tomography snapshot.


title: The Capacity on Degraded Relay Broadcast Channel

abstract: The relay broadcast channel (RBC) is considered, in which a transmitter communicates with two receivers with the assistance of a relay. Based on different degradation orders among the relay and the receivers' outputs, three types of physically degraded RBCs (PDRBCs) are introduced. Inner bounds and outer bounds are derived on the capacity region of the presented three types. The bounds are tight for two types of PDRBCs: 1) one receiver's output is a degraded form of the other receiver's output, and the relay's output is a degraded form of the weaker receiver's output; 2) one receiver's output is a degraded form of the relay's output, and the other receiver's output is a degraded form of the relay's output. For the Gaussian PDRBC, the bounds match, i.e., establish its capacity region.


title: Continuum Percolation of Polydisperse Rods in Quadrupole Fields: Theory and Simulations

abstract: We investigate percolation in mixtures of nanorods in the presence of external fields that align or disalign the particles with the field axis. Such conditions are found in the formulation and processing of nanocomposites, where the field may be electric, magnetic, or due to elongational flow. Our focus is on the effect of length polydispersity, which -- in the absence of a field -- is known to produce a percolation threshold that scales with the inverse weight average of the particle length. Using a model of non-interacting spherocylinders in conjunction with connectedness percolation theory, we show that a quadrupolar field always increases the percolation threshold and that the universal scaling with the inverse weight average no longer holds if the field couples to the particle length. Instead, the percolation threshold becomes a function of higher moments of the length distribution, where the order of the relevant moments crucially depends on the strength and type of field applied. The theoretical predictions compare well with the results of our Monte Carlo simulations, which eliminate finite size effects by exploiting the fact that the universal scaling of the wrapping probability function holds even in anisotropic systems. Theory and simulation demonstrate that the percolation threshold of a polydisperse mixture can be lower than that of the individual components, confirming recent work based on a mapping onto a Bethe lattice as well as earlier computer simulations involving dipole fields. Our work shows how the formulation of nanocomposites may be used to compensate for the adverse effects of aligning fields that are inevitable under practical manufacturing conditions.


title: Quantum uncertainty relation based on the mean deviation

abstract: Traditional forms of quantum uncertainty relations are invariably based on the standard deviation. This can be understood in the historical context of simultaneous development of quantum theory and mathematical statistics. Here, we present alternative forms of uncertainty relations, in both state dependent and state independent forms, based on the mean deviation. We illustrate the robustness of this formulation in situations where the standard deviation based uncertainty relation is inapplicable. We apply the mean deviation based uncertainty relation to detect EPR violation in a lossy scenario for a higher inefficiency threshold than that allowed by the standard deviation based approach. We demonstrate that the mean deviation based uncertainty relation can perform equally well as the standard deviation based uncertainty relation as non-linear witness for entanglement detection.


title: Tight-binding Piezoelectric Theory and Electromechanical Coupling Correlations for Transition Metal Dichalcogenide Monolayers

abstract: The lack of inversion symmetry in semiconducting transition metal dichalcogenide monolayers (TMDMs) enables a considerable molecular-level intrinsic piezoelectricity, which opens prospects for atomically-thin piezotronics and optoelectronics. Here, based on the tight-binding (TB) approach and Berry phase polarization theory, we establish an atomic-scale TB theory for demonstrating piezoelectric physics in TMDMs. Using the TB piezoelectric theory, we predict their electronic Gr\"{u}neisen parameters (EGP) which measure the electron-phonon couplings. By virtue of the constructed analytical piezoelectric model, we further reveal the correlation between the electronic contribution to piezoelectric coefficients and strain-induced pseudomagnetic gauge field (PMF). These predicted EGP and PMF for TMDMs are experimentally testable, and hence the TB piezoelectric model is an alternative theoretical framework for calculating electron-phonon interactions and PMF.


title: Non-Hermitian gauged topological laser arrays

abstract: Stable and phase-locked emission in an extended topological supermode of coupled laser arrays, based on concepts of non-Hermitian and topological photonics, is theoretically suggested. We consider a non-Hermitian Su-Schrieffer-Heeger chain of coupled microring resonators and show that application of a synthetic imaginary gauge field via auxiliary passive microrings leads to all supermodes of the chain, except one, to become edge states. The only extended supermode, that retains some topological protection, can stably oscillate suppressing all other non-topological edge supermodes. Numerical simulations based on a rate equation model of semiconductor laser arrays confirm stable anti-phase laser emission in the extended topological supermode and the role of the synthetic gauge field to enhance laser stability.


title: Expectations of the Cosmic Antideuteron Flux

abstract: The cosmic antideuteron is a promising probe for the dark matter annihilation signature. In order to determine the DM signature, the background astrophysical antideuteron flux should be carefully studied. In this work we provide a new calculation of the secondary antideuteron flux, and pay special attention to the uncertainties from hadronic interaction models by using several Monte Carlo generators. The uncertainties from propagation effects are also carefully investigated for both the astrophysical background and DM annihilation signature in several scenarios, which are constrained by the latest B/C ratio measured by AMS-02. Considering these uncertainties, we find that the secondary antideuteron flux is hard to detect in the near future detectors. However, the antideuteron signature from dark matter annihilation will be detectable even considering the constraint from the AMS-02 observation of the $\bar{p}/p$ ratio.


title: Perspectives on scaling and multiscaling in passive scalar turbulence

abstract: We revisit the well-known problem of multiscaling in substances passively advected by homogeneous and isotropic turbulent flows or passive scalar turbulence. To that end we propose a two-parameter continuum hydrodynamic model for an advected substance concentration $\theta$, parametrised jointly by $y$ and $\overline y$, that characterise the spatial scaling behaviour of the variances of the advecting stochastic velocity and the stochastic additive driving force, respectively. We analyse it within a one-loop dynamic renormalisation group method to calculate the multiscaling exponents of the equal-time structure functions of $\theta$. We show how the interplay between the advective velocity and the additive force may lead to simple scaling or multiscaling. In one limit, our results reduce to the well-known results from the Kraichnan model for passive scalar. Our framework of analysis should be of help for analytical approaches for the still intractable problem of fluid turbulence itself.


title: Unitary $n$-designs via random quenches in atomic Hubbard and Spin models: Application to the measurement of R\'enyi entropies

abstract: We present a general framework for the generation of random unitaries based on random quenches in atomic Hubbard and spin models, forming approximate unitary $n$-designs, and their application to the measurement of R\'enyi entropies. We generalize our protocol presented in [Elben2017: arXiv:1709.05060, to appear in Phys. Rev. Lett.] to a broad class of atomic and spin lattice models. We further present an in-depth numerical and analytical study of experimental imperfections, including the effect of decoherence and statistical errors, and discuss connections of our approach with many-body quantum chaos.


title: ViZDoom: DRQN with Prioritized Experience Replay, Double-Q Learning, & Snapshot Ensembling

abstract: ViZDoom is a robust, first-person shooter reinforcement learning environment, characterized by a significant degree of latent state information. In this paper, double-Q learning and prioritized experience replay methods are tested under a certain ViZDoom combat scenario using a competitive deep recurrent Q-network (DRQN) architecture. In addition, an ensembling technique known as snapshot ensembling is employed using a specific annealed learning rate to observe differences in ensembling efficacy under these two methods. Annealed learning rates are important in general to the training of deep neural network models, as they shake up the status-quo and counter a model's tending towards local optima. While both variants show performance exceeding those of built-in AI agents of the game, the known stabilizing effects of double-Q learning are illustrated, and priority experience replay is again validated in its usefulness by showing immediate results early on in agent development, with the caveat that value overestimation is accelerated in this case. In addition, some unique behaviors are observed to develop for priority experience replay (PER) and double-Q (DDQ) variants, and snapshot ensembling of both PER and DDQ proves a valuable method for improving performance of the ViZDoom Marine.


title: A Quinone Based Single-Molecule Switch as Building Block for Molecular Electronics

abstract: Azophenine has previously been identified as a controllable molecular switch when deposited on a Cu(110) surface, where it can be in two symmetry-equivalent states. Each of the two states can be set as well as read by means of a scanning tunneling microscope (STM). We propose a family of molecules based on the same quinone core, which show similar switching behavior without a supporting metal surface. Such a molecule could be an element in a molecular circuit or computer. Using the example of a simple hypothetical molecule, we show that it is possible to create molecules that show the necessary properties: two conformations with similar energy but different electric conductivity, and the possibility to switch between those by applying an external electric field.


title: Collective motion of Active Brownian Particles with polar alignment

abstract: We present a comprehensive computational study of the collective behavior emerging from the competition between self-propulsion, excluded volume interactions and velocity-alignment in a two-dimensionnal model of active particles. We consider an extension of the Active Brownian Particles model where the self-propulsion direction of the particles aligns with the one of their neighbors. We analyze the onset of collective motion (flocking) in a low-density regime (10% surface area) and show that it is mainly controlled by the strength of velocity-alignment interactions: the competition between self-propulsion and crowding effects plays a minor role in the emergence of flocking. However, above the flocking threshold, the system presents a richer pattern formation scenario than analogous models without alignment interactions (Active Brownian Particles) or excluded volume effects (Vicsek-like models). Depending on the parameter regime, the structure of the system is characterized by either a broad distribution of finite-sized polar clusters or the presence of an amorphous, highly fluctuating, large-scale traveling structure which can take a lane-like or band-like form (and usually a hybrid structure which is halfway in between both). We establish a phase diagram that summarizes collective behavior of polar Active Brownian Particles and propose a generic mechanism to describe the complexity of the large-scale structures observed in systems of repulsive self-propelled particles.


title: An Information Analysis on Modeling Interaction Effects in Logistic Regression

abstract: The Akaike information criterion (AIC) is commonly used to select a logistic regression model for optimal prediction of a binary response by a specified family of models. It however lacks a convincing method of prescribing a proper family of models using the desired predictors and their interaction effects. For an alternative approach to model selection, we propose a direct selection scheme which first identifies the indispensable regressors as main-effect predictors, then examines significant interaction effects between the selected predictors such that a logistic model is constructed. The two-step selection scheme is formulated by testing for valid information identity between the response and the predictors, from which the most parsimonious logistic model is derived from the least set of indispensable predictors and interaction effects. As a byproduct, the minimum AIC model is easily found in a neighborhood of the selected model. The scheme is employed to yield the logistic model for predicting the acquisition of professional licenses in a survey of employed youth workers.


title: Structures and magnetic properties of iron silicide from adaptive genetic algorithm and first-principles calculations

abstract: We performed a systematic search for low-energy structures of binary iron silicide over a wide range of compositions using the crystal structure prediction method based on adaptive genetic algorithm. 36 structures with formation energies within 50 meV/atom (11 of them are within 20 meV) above the convex hull formed by experimentally known stable structures are predicted. Magnetic properties of these low-energy structures are investigated. Some of these structures can be promising candidates for rare-earth-free permanent magnet.


title: An interacting galaxy pair at the origin of a light echo

abstract: In a low-density region of the Shapley supercluster we identified an interacting galaxy pair at redshift z=0.04865 in which the Seyfert-2 nucleus of the main galaxy (ShaSS 073) is exciting an extended emission line region (EELR, ~170 kpc^2) in the disk of the less massive companion (ShaSS 622). New integral-field spectroscopy and the multi-band data-set, spanning from far-ultraviolet to far-infrared and radio wavelengths, allowed us to obtain a detailed description of the ShaSS 622-073 system. The gas kinematics shows hints of interaction, although the overall velocity field shows a quite regular rotation in both galaxies, thus suggesting that we are observing their first encounter as confirmed by the estimated distance of 21 kpc between the two galaxy centers. The detected ~2-3 kpc AGN outflow and the geometry of the EELR in ShaSS 622 support the presence of a hollow bicone structure. The status and sources of the ionization across the whole system have been analysed through photoionization models and a Bayesian approach which prove a clear connection between the AGN and the EELR. The luminosity of the AGN (2.4x10^44 erg/s) is a factor 20 lower than the power needed to excite the gas in the EELR (4.6x10^45 erg/s) indicating a dramatic fading of the AGN in the last 3x10^4 yr. ShaSS 073-622 provides all the ingredients listed in the recipe of a light echo where a highly ionised region maintains memory of a preceding more energetic phase of a now faded AGN. This is the first case of a light echo observed between two galaxies.


title: The risk model with stochastic premiums, dependence and a threshold dividend strategy

abstract: The paper deals with a generalization of the risk model with stochastic premiums where dependence structures between claim sizes and inter-claim times as well as premium sizes and inter-premium times are modeled by Farlie--Gumbel--Morgenstern copulas. In addition, dividends are paid to its shareholders according to a threshold dividend strategy. We derive integral and integro-differential equations for the Gerber--Shiu function and the expected discounted dividend payments until ruin. Next, we concentrate on the detailed investigation of the model in the case of exponentially distributed claim and premium sizes. In particular, we find explicit formulas for the ruin probability in the model without either dividend payments or dependence as well as for the expected discounted dividend payments in the model without dependence. Finally, numerical illustrations are presented.


title: A Comprehensive Bayesian Treatment of the Universal Kriging model with Mat\'ern correlation kernels

abstract: The Gibbs reference posterior distribution provides an objective full-Bayesian solution to the problem of prediction of a stationary Gaussian process with Mat\'ern anisotropic kernel. A full-Bayesian approach is possible, because the posterior distribution is expressed as the invariant distribution of a uniformly ergodic Markovian kernel for which we give an explicit expression. In this paper, we show that it is appropriate for the Universal Kriging framework, that is when an unknown function is added to the stationary Gaussian process. We give sufficient conditions for the existence and propriety of the Gibbs reference posterior that apply to a wide variety of practical cases and illustrate the method with several examples. Finally, simulations of Gaussian processes suggest that the Gibbs reference posterior has good frequentist properties in terms of coverage of prediction intervals.


title: Structural and Chemical Orders in Ni64.5Zr35.5 Metallic Glass by Molecular Dynamics Simulation

abstract: The atomic structure of Ni64.5Zr35.5 metallic glass has been investigated by molecular dynamics (MD) simulations. The calculated structure factors from the MD glassy sample at room temperature agree well with the X-ray diffraction (XRD) and neutron diffraction (ND) experimental data. Using the pairwise cluster alignment and clique analysis methods, we show that there are three types dominant short-range order (SRO) motifs around Ni atoms in the glass sample of Ni64.5Zr35.5, i.e., Mixed-Icosahedron(ICO)-Cube, Twined-Cube and icosahedron-like clusters. Furthermore, chemical order and medium-range order (MRO) analysis show that the Mixed-ICO-Cube and Twined-Cube clusters exhibit the characteristics of the crystalline B2 phase. Our simulation results suggest that the weak glass-forming ability (GFA) of Ni64.5Zr35.5 can be attributed to the competition between the glass forming ICO SRO and the crystalline Mixed-ICO-Cube and Twined-Cube motifs.


title: Origin and Control of Polyacrylonitrile Alignments on Carbon Nanotube and Graphene Nanoribbon

abstract: While one of the most promising applications of carbon nanotubes (CNTs) is to enhance polymer orientation and crystallization to achieve advanced carbon fibers, the successful realization of this goal has been hindered by the insufficient atomistic understanding of polymer-CNT interfaces. We herein theoretically study polyacrylonitrile (PAN)-CNT hybrid structures as a representative example of polymer-CNT composites. Based on density-functional theory calculations, we first find that the relative orientation of polar PAN nitrile groups with respect to the CNT surface is the key factor that determines the PAN-CNT interface energetics and the lying-down PAN configurations are much more preferable than their standing-up counterparts. The CNT curvature is identified as another important factor, giving the largest binding energy in the zero-curvature graphene limit. Charge transfer analysis explains the unique tendency of linear PAN alignments on the CNT surface and the possibility of ordered PAN-PAN assembly. Next, performing large-scale molecular dynamics simulations, we show that the desirable linear PAN-CNT alignment can be achieved even for relatively large initial misorientations and further demonstrate that graphene nanoribbons are a promising carbon nano-reinforcement candidate. The microscopic understanding accumulated in this study will provide design guidelines for the development of next-generation carbon nanofibers.


title: Ion-exchange synthesis and superconductivity at 8.6 K of Na2Cr3As3 with quasi-one-dimensional crystal structure

abstract: A new Cr-based quasi-one-dimensional superconductor Na2Cr3As3 was synthesized by an ion-exchange method in sodium naphthalenide solution. The crystals are thread-like and the structure was analyzed by X-ray diffraction with a noncentrosymmetric hexagonal space group P-6m2 (No. 187), in which the (Cr3As3)2- linear chains are separated by Na+ ions, and the refined lattice parameters are a = 9.239(2) {\AA} and c = 4.209(6) {\AA}. The measurements for electrical resistivity, magnetic susceptibility, and heat capacity reveal a superconducting transition with unconventional characteristic at the Tc of 8.6 K, which exceeds that of all previously reported Cr-based superconductors.


title: Connections between a system of Forward-Backward SDEs and Backward Stochastic PDEs related to the utility maximization problem

abstract: Connections between a system of Forward-Backward SDEs and Backward Stochastic PDEs related to the utility maximiza- tion problem is established. Besides, we derive another version of FBSDE of the same problem and prove an existence of a solution


title: Graph Pattern Matching for Dynamic Team Formation

abstract: Finding a list of k teams of experts, referred to as top-k team formation, with the required skills and high collaboration compatibility has been extensively studied. However, existing methods have not considered the specific collaboration relationships among different team members, i.e., structural constraints, which are typically needed in practice. In this study, we first propose a novel graph pattern matching approach for top-k team formation, which incorporates both structural constraints and capacity bounds. Second, we formulate and study the dynamic top-k team formation problem due to the growing need of a dynamic environment. Third, we develop an unified incremental approach, together with an optimization technique, to handle continuous pattern and data updates, separately and simultaneously, which has not been explored before. Finally, using real-life and synthetic data, we conduct an extensive experimental study to show the effectiveness and efficiency of our graph pattern matching approach for (dynamic) top-k team formation.


title: Instrumental variables estimation with competing risk data

abstract: Time-to-event analyses are often plagued by both -- possibly unmeasured -- confounding and competing risks. To deal with the former, the use of instrumental variables for effect estimation is rapidly gaining ground. We show how to make use of such variables in competing risk analyses. In particular, we show how to infer the effect of an arbitrary exposure on cause-specific hazard functions under a semi-parametric model that imposes relatively weak restrictions on the observed data distribution. The proposed approach is flexible accommodating exposures and instrumental variables of arbitrary type, and enables covariate adjustment. It makes use of closed-form estimators that can be recursively calculated, and is shown to perform well in simulation studies. We also demonstrates its use in an application on the effect of mammography screening on the risk of dying from breast cancer


title: Finite asymptotic clusters of metric spaces

abstract: Let $(X, d)$ be an unbounded metric space and let $\tilde r=(r_n)_{n\in\mathbb N}$ be a sequence of positive real numbers tending to infinity. A pretangent space $\Omega_{\infty, \tilde r}^{X}$ to $(X, d)$ at infinity is a limit of the rescaling sequence $\left(X, \frac{1}{r_n}d\right).$ The set of all pretangent spaces $\Omega_{\infty, \tilde r}^{X}$ is called an asymptotic cluster of pretangent spaces. Such a cluster can be considered as a weighted graph $(G_{X, \tilde r}, \rho_{X})$ whose maximal cliques coincide with $\Omega_{\infty, \tilde r}^{X}$ and the weight $\rho_{X}$ is defined by metrics on $\Omega_{\infty, \tilde r}^{X}$. We describe the structure of metric spaces having finite asymptotic clusters of pretangent spaces and characterize the finite weighted graphs which are isomorphic to these clusters.


title: An accurate centroid algorithm for PSF reconstruction

abstract: In this work, we present a novel centroiding method based on Fourier space Phase Fitting(FPF) for Point Spread Function(PSF) reconstruction. We generate two sets of simulations to test our method. The first set is generated by GalSim with elliptical Moffat profile and strong anisotropy which shifts the center of the PSF. The second set of simulation is drawn from CFHT i band stellar imaging data. We find non-negligible anisotropy from CFHT stellar images, which leads to $\sim$0.08 scatter in unit of pixels using polynomial fitting method Vakili and Hogg (2016). And we apply FPF method to estimate the centroid in real space, this scatter reduces to $\sim$0.04 in SNR=200 CFHT like sample. In low SNR (50 and 100) CFHT like samples, the background noise dominates the shifting of the centroid, therefore the scatter estimated from different methods are similar. We compare polynomial fitting and FPF using GalSim simulation with optical anisotropy. We find that in all SNR$\sim$50, 100 and 200) samples, FPF performs better than polynomial fitting by a factor of $\sim$3. In general, we suggest that in real observations there are anisotropy which shift the centroid, and FPF method is a better way to accurately locate it.


title: Double barrier reflected BSDEs with stochastic Lipschitz coefficient

abstract: This paper proves the existence and uniqueness of a solution to doubly reflected backward stochastic differential equations where the coefficient is stochastic Lipschitz, by means of the penalization method.


title: Particle Clustering Machine: A Dynamical System Based Approach

abstract: Identification of the clusters from an unlabeled data set is one of the most important problems in Unsupervised Machine Learning. The state of the art clustering algorithms are based on either the statistical properties or the geometric properties of the data set. In this work, we propose a novel method to cluster the data points using dynamical systems theory. After constructing a gradient dynamical system using interaction potential, we prove that the asymptotic dynamics of this system will determine the cluster centers, when the dynamical system is initialized at the data points. Most of the existing heuristic-based clustering techniques suffer from a disadvantage, namely the stochastic nature of the solution. Whereas, the proposed algorithm is deterministic, and the outcome would not change over multiple runs of the proposed algorithm with the same input data. Another advantage of the proposed method is that the number of clusters, which is difficult to determine in practice, does not have to be specified in advance. Simulation results with are presented, and comparisons are made with the existing methods.


title: Pierre Auger Observatory and Telescope Array: Joint Contributions to the 35th International Cosmic Ray Conference (ICRC 2017)

abstract: Joint contributions of the Telescope Array Collaboration and the Pierre Auger Collaboration to the 35th International Cosmic Ray Conference (ICRC 2017), 12-20 July 2017, Bexco, Busan, Korea.


title: Proteomics Analysis of FLT3-ITD Mutation in Acute Myeloid Leukemia Using Deep Learning Neural Network

abstract: Deep Learning can significantly benefit cancer proteomics and genomics. In this study, we attempt to determine a set of critical proteins that are associated with the FLT3-ITD mutation in newly-diagnosed acute myeloid leukemia patients. A Deep Learning network consisting of autoencoders forming a hierarchical model from which high-level features are extracted without labeled training data. Dimensional reduction reduced the number of critical proteins from 231 to 20. Deep Learning found an excellent correlation between FLT3-ITD mutation with the levels of these 20 critical proteins (accuracy 97%, sensitivity 90%, specificity 100%). Our Deep Learning network could hone in on 20 proteins with the strongest association with FLT3-ITD. The results of this study allow a novel approach to determine critical protein pathways in the FLT3-ITD mutation, and provide proof-of-concept for an accurate approach to model big data in cancer proteomics and genomics.


title: Top Quark Pair Production: theory overview

abstract: In this talk we present the current status of the theoretical description of the production of top-quark pairs at the LHC, detailing recent progress for predictions at the stable-top level as well as for predictions that include the decay of the top quarks.


title: The Unified Astronomy Thesaurus: Semantic Metadata for Astronomy and Astrophysics

abstract: Several different controlled vocabularies have been developed and used by the astronomical community, each designed to serve a specific need and a specific group. The Unified Astronomy Thesaurus (UAT) attempts to provide a highly structured controlled vocabulary that will be relevant and useful across the entire discipline, regardless of content or platform. As two major use cases for the UAT include classifying articles and data, we examine the UAT in comparison with the Astronomical Subject Keywords used by major publications and the JWST Science Keywords used by STScI's Astronomer's Proposal Tool.


title: Shannon Capacity is Achievable for Binary Interactive First-Order Markovian Protocols

abstract: We address the problem of simulating an arbitrary binary interactive first-order Markovian protocol over a pair of binary symmetric channels with crossover probability $\varepsilon$. We are interested in the achievable rates of reliable simulation, i.e., in characterizing the smallest possible blowup in communications such that a vanishing error probability (in the protocol length) can be attained. Whereas for general interactive protocols the output of each party may depend on all previous outputs of its counterpart, in a (first-order) Markovian protocol this dependence is limited to the last observed output only. Previous works in the field discuss broader families of protocols but assess the achievable rates only at the limit where $\varepsilon\to0$. In this paper, we prove that the one-way Shannon capacity, $1-h(\varepsilon)$, can be achieved for any binary first-order Markovian protocol. This surprising result, is to the best of our knowledge, the first example in which non-trivial interactive protocol can be simulated in the Shannon capacity. We give two capacity achieving coding schemes, which both divide the protocol into independent blocks, and implement vertical block coding. The first scheme is based on a random separation into blocks with variable lengths. The second scheme is based on a deterministic separation into blocks, and efficiently predicting their last transmission. The prediction can be regarded as a binary pointer jumping game, for which we show that the final step can be calculated with $O(\log m)$ bits, where $m$ is the number of rounds in the game. We conclude the work by discussing possible extensions of the result to higher order models.


title: Calder\'on-Zygmund operators on Zygmund spaces on domains

abstract: Given a bounded Lipschitz domain $D\subset \mathbb{R}^d$ and a Calder\'on-Zygmund operator $T$, we study the relations between smoothness properties of $\partial D$ and the boundedness of $T$ on the Zydmund space $\mathcal{C}_{\omega}(D)$ defined for a general growth function $\omega$. In the proof we obtain a T(P) theorem for the Zygmund spaces, when one checks boundedness not only of the characteristic function, but a finite collection of polynomials restricted to the domain. Also, a new form of extra cancellation property of the even Calder\'on-Zygmund operators in polynomial domains is stated.


title: Relativistic quantum motion of spin-0 particles under the influence of non-inertial effects in the cosmic string space-time

abstract: We study solutions for the Klein-Gordon equation with vector and scalar potentials of the Coulomb types under the influence of non-inertial effects in the space-time of topological defects. We also investigate a quantum particle described by the Klein-Gordon oscillator in the background space-time generated by a string. An important result obtained is that the non-inertial effects restrict the physical region of the space-time where the particle can be placed. In addition, we show that these potentials can form bound states for the relativistic wave equation equation in this kind of background.


title: A Study of Bug Resolution Characteristics in Popular Programming Languages

abstract: This paper presents a large-scale study that investigates the bug resolution characteristics among popular Github projects written in different programming languages. We explore correlations but, of course, we cannot infer causation. Specifically, we analyse bug resolution data from approximately 70 million Source Line of Code, drawn from 3 million commits to 600 GitHub projects, primarily written in 10 programming languages. We find notable variations in apparent bug resolution time and patch (fix) size. While interpretation of results from such large-scale empirical studies is inherently difficult, we believe that the differences in medians are sufficiently large to warrant further investigation, replication, re-analysis and follow up research. For example, in our corpus, the median apparent bug resolution time (elapsed time from raise to resolve) for Ruby was 4X that for Go and 2.5X for Java. We also found that patches tend to touch more files for the corpus of strongly typed and for statically typed programs. However, we also found evidence for a lower elapsed resolution time for bug resolution committed to projects constructed from statically typed languages. These findings, if replicated in subsequent follow on studies, may shed further empirical light on the debate about the importance of static typing.


title: A remark on the Sibony function

abstract: We present an effective formula for the Sibony function for all Reinhardt domains.


title: Optimal density for values of generic polynomial maps

abstract: We establish that the optimal bound for the size of the smallest integral solution of the Oppenheim Diophantine approximation problem $\abs{Q(x)-\xi}< \epsilon$ for a generic ternary form $Q$ is $\abs{x}\ll \epsilon^{-1}$. We also establish an optimal rate of density for the values of polynomials maps in a number of other natural problems, including the values of linear forms restricted to suitable quadratic surfaces, and the values of the polynomial map defined by the generators of the ring of conjugation-invariant polynomials on $M_3(\C)$. These results are instances of a general approach that we develop, which considers a rational affine algebraic subvariety of Euclidean space, invariant and homogeneous under an action of a semisimple Lie group $G$. Given a polynomial map $F$ defined on the Euclidean space which is invariant under a semisimple subgroup $H$ of the acting group $G$, consider the family of its translates $F\circ g$ by elements of the group. We study the restriction of these polynomial functions to the integer points on the variety confined to a large Euclidean ball. Our main results establish an explicit rate of density for their values, for generic polynomials in the family. This problem has been extensively studied before when the polynomials in question are linear, in the context of classical Diophantine approximation, but very little was known about it for polynomial of higher degree. We formulate a heuristic pigeonhole lower bound for the density and an explicit upper bound for it, formulate a sufficient condition for the coincidence of the lower and upper bounds, and in a number of natural examples establish that they indeed match. Finally, we also establish a rate of density for values of homogeneous polynomials on homogeneous projective varieties.


title: Existence of $C^\alpha$ solutions to integro-PDEs

abstract: This paper is concerned with existence of a $C^{\alpha}$ viscosity solution of a second order non-translation invariant integro-PDE. We first obtain a weak Harnack inequality for such integro-PDE. We then use the weak Harnack inequality to prove H\"older regularity and existence of solutions of the integro-PDEs.


title: Bulk viscous embedded hybrid dark energy models

abstract: In this paper, we have constructed the cosmological model of the universe in a two fluids environment with a newly developed mathematical formalism. In order to construct the model Binachi type V (BV) space time is considered with a time varying deceleration parameter. Both the fluids, the viscous fluid and the dark energy (DE) fluid have shown their dominance respectively in early time and late time of the cosmic evolution. The scale factor that simulates the cosmic transition based on the value of the bulk viscous coefficient. Within the developed formalism, a general form of the skewness parameters is also obtained as a functional form of the scale factor. The physical parameter of the model such as equation of state (EoS) parameter is also derived and analysed. The state finder diagnostic pair are also obtained to understand the geometrical nature of the model.


title: Dissipative measure valued solutions for general conservation laws

abstract: In the last years measure-valued solutions started to be considered as a relevant notion of solutions if they satisfy the so-called measure-valued -- strong uniqueness principle. This means that they coincide with a strong solution emanating from the same initial data if this strong solution exists. This property has been examined for many systems of mathematical physics, including incompressible and compressible Euler system, compressible Navier-Stokes system et al. and there are also some results concerning general hyperbolic systems. Our goal is to provide a unified framework for general systems, that would cover the most interesting cases of systems, and most importantly, we give examples of equations, for which the aspect of measure-valued -- strong uniqueness has not been considered before, like incompressible magentohydrodynamics and shallow water magnetohydrodynamics.


title: On local stabilities of $p$-K\"ahler structures

abstract: By use of a natural extension map and a power series method, we obtain a local stability theorem for p-K\"ahler structures with the $(p,p+1)$-th mild $\partial\bar\partial$-lemma under small differentiable deformations.


title: The $C^1$ property of convex carrying simplices for competitive maps

abstract: For a class of competitive maps there is an invariant one-codimensional manifold (the carrying simplex) attracting all non-trivial orbits. In the present paper it is shown that its convexity implies that it is a $C^1$ submanifold-with-corners, neatly embedded in the non-negative orthant. The proof uses the characterization of neat embedding in terms of inequalities between Lyapunov exponents for ergodic invariant measures supported on the boundary of the carrying simplex.


title: Two-Loop master integrals for heavy-to-light form factors of two different massive fermions

abstract: We calculate the full set of the two-loop master integrals for heavy-to-light form factors of two different massive fermions for arbitrary momentum transfer in NNLO QCD or QED corrections. These integrals allow to determine the two-loop QCD or QED corrections to the amplitudes for heavy-to-light form factors of two massive fermions in a full analytical way, without any approximations. The analytical results of the master integrals are derived using the method of differential equations, along with a proper choosing of canonical basis for the master integrals. All the results of master integrals are expressed in terms of Goncharov polylogarithms.


title: Computing rotation numbers in open books

abstract: We give explicit formulas and algorithms for the computation of the rotation number of a nullhomologous Legendrian knot on a page of a contact open book. On the way, we derive new formulas for the computation of the Thurston-Bennequin invariant of such knots and the Euler class and the d3-invariant of the underlying contact structure.


title: Local probabilities of randomly stopped sums of power law lattice random variables

abstract: Let $X_1$ and $N\ge 0$ be integer valued power law random variables. For a randomly stopped sum $S_N=X_1+\cdots+X_N$ of independent and identically distributed copies of $X_1$ we establish a first order asymptotics of the local probabilities $P(S_N=t)$ as $t\to+\infty$. Using this result we show the $k^{-\delta}$, $0\le \delta\le 1$ scaling of the local clustering coefficient (of a randomly selected vertex of degree $k$) in a power law affiliation network.


title: Counter-streaming flows in a giant quiet-Sun filament observed in the extreme ultraviolet

abstract: A giant solar filament was visible on the solar surface from 2011 Nov. 8-23. The filament stretched over more than half a solar diameter. Multiwavelength data from the SDO instrument AIA (171, 193, 304, and 211\AA) were used to examine counter-streaming flows within the spine of the filament. H$\alpha$ images from the Kanzelh\"ohe Solar Observatory provided context information. We apply local correlation tracking (LCT) to a two-hour time series on 2011 Nov. 16 of the AIA images to derive horizontal flow velocities of the filament. To enhance the contrast of the AIA images, noise adaptive fuzzy equalization is employed, which allows us to identify and quantify counter-streaming flows in the filament. We detect counter-streaming flows in the filament, which are visible in the time-lapse movies in all examined AIA wavelength bands. In the time-lapse movies we see that these persistent flows lasted for at least two hours. Furthermore, by applying LCT to the images we clearly determine counter-streaming flows in time series of 171 and 193\AA\ images. In the 304\AA\ wavelength band, we only see minor indications for counter-streaming flows with LCT, while in the 211\AA\ wavelength band the counter-streaming flows are not detectable. The average horizontal flows reach mean flow speeds of 0.5 km/s. The highest horizontal flow speeds are identified in the 171\AA\ band with flow speeds of up to 2.5 km/s. The results are averaged over a time series of 90 min. Because the LCT sampling window has finite width, a spatial degradation cannot be avoided leading to lower estimates of the flow velocities as compared to feature tracking or Doppler measurements. The counter-streaming flows cover about 15-20% of the whole area of the EUV filament channel and are located in the central part of the spine. In conclusion, we confirm the omnipresence of counter-streaming flows also in giant quiet-Sun filaments.


title: Distributed Memory Techniques for Classical Simulation of Quantum Circuits

abstract: In this paper we describe, implement, and test the performance of distributed memory simulations of quantum circuits on the MSU Laconia Top500 supercomputer. Using OpenMP and MPI hybrid parallelization, we first use a distributed matrix-vector multiplication with one-dimensional partitioning and discuss the shortcomings of this method due to the exponential memory requirements in simulating quantum computers. We then describe a more efficient method that stores only the $2^n$ amplitudes of the $n$ qubit state vector $|\psi\rangle$ and optimize its single node performance. In our multi-node implementation, we use a single amplitude communication protocol that maximizes the number of qubits able to be simulated and minimizes the ratio of qubits that require communication to those that do not, and we present an algorithm for efficiently determining communication pairs among processors. We simulate up to 30 qubits on a single node and 33 qubits with the state vector partitioned across 64 nodes. Lastly, we discuss the advantages and disadvantages of our communication scheme, propose potential improvements, and describe other optimizations such as storing the state vector non-sequentially in memory to map communication requirements to idle qubits in the circuit.


title: Mechanical, optoelectronic and transport properties of single-layer Ca2N and Sr2N electrides

abstract: Electride materials offer attractive physical properties due to their loosely bound electrons. Ca2N, an electride in the two-dimensional (2D) form was successfully recently synthesized. We conducted extensive first-principles calculations to explore the mechanical, electronic, optical and transport response of single-layer and free-standing Ca2N and Sr2N electrides to external strain. We show that Ca2N and Sr2N sheets present isotropic elastic properties with positive Poisson's ratios, however, they yield around 50% higher tensile strength along the zigzag direction as compared with armchair. We also showed that the strain has negligible effect on the conductivity of the materials; the current in the system reduces by less than 32% for the structure under ultimate uniaxial strain along the armchair direction. Compressive strain always increases the electronic transport in the systems due to stronger overlap of the atomic orbitals. Our results show that the optical spectra are anisotropic for light polarization parallel and perpendicular to the plane. Interband transition contributions along in-plane polarization are not negligible, by considering this effect the optical properties of Ca2N and Sr2N sheets in the low frequency regime significantly changed. The insight provided by this study can be useful for the future application of Ca2N and Sr2N in nanodevices.


title: Computing the Inverse Mellin Transform of Holonomic Sequences using Kovacic's Algorithm

abstract: We describe how the extension of a solver for linear differential equations by Kovacic's algorithm helps to improve a method to compute the inverse Mellin transform of holonomic sequences. The method is implemented in the computer algebra package HarmonicSums.


title: Interpretation of increased energetic particle flux measurements by SEPT aboard the STEREO spacecraft and contamination

abstract: Context. Interplanetary (IP) shocks are known to be accelerators of energetic charged particles observed in-situ in the heliosphere. However, the acceleration of near-relativistic electrons by shocks in the interplanetary medium is often questioned. On 9 August 2011 a Corotating Interaction Region (CIR) passed STEREO B (STB) that resulted in a flux increase in the electron and ion channels of the Solar Electron and Proton Telescope (SEPT). Because electron measurements in the few keV to several 100 keV range rely on the so-called magnet foil technique, which is utilized by SEPT, ions can contribute to the electron channels. Aims. We aim to investigate whether the flux increase in the electron channels of SEPT during the CIR event on 9 August 2011 is caused by ion contamination only. Methods. We compute the SEPT response functions for protons and helium utilizing an updated GEANT4 model of SEPT. The CIR energetic particle ion spectra for protons and helium are assumed to follow a Band function in energy per nucleon with a constant helium to proton ratio. Results. Our analysis leads to a helium to proton ratio of 16.9% and a proton flux following a Band function with the parameters $I_0 = 1.24 \cdot 10^4$ / (cm2 s sr MeV/nuc.), $E_c = 79$ keV/nuc. and spectral indices of $\gamma_1 = -0.94$ and $\gamma_2 = -3.80$ which are in good agreement with measurements by the Suprathermal Ion Telescope (SIT) aboard STB. Conclusions. Since our results explain the SEPT measurements, we conclude that no significant amount of electrons were accelerated between $55$ keV and $425$ keV by the CIR.


title: Almost h-conformal slant submersions from almost quaternionic Hermitian manifolds

abstract: We introduce the notions of h-conformal slant submersions and almost h-conformal slant submersions from almost quaternionic Hermitian manifolds onto Riemannian manifolds as a generalization of Riemannian submersions, horizontally conformal submersions, slant submersions, h-slant submersions, almost h-slant submersion and conformal slant submersions. We investigate several properties of these including the integrability of distributions, the geometry of foliations and the conditions for such maps to be totally geodesic. Further we give some examples of such maps.


title: Highly anisotropic strain dependencies in PrIr$_2$Zn$_{20}$

abstract: We report thermal expansion and magnetostriction of the cubic non-Kramers system PrIr$_2$Zn$_{20}$ with a non-magnetic $\varGamma_{3}$ ground state doublet. In previous experiments, antiferroquadrupolar order at \hbox{$T_{\mathrm{Q}}=0.11$\,K} and a Fermi liquid state around $B_{\mathrm{c}}\approx5$\,T for \hbox{$\boldsymbol{B}\parallel[001]$}, indicative of possible ferrohastatic order, were discovered. For magnetic fields \hbox{$\boldsymbol{B}\parallel[001]$}, the low temperature longitudinal and transverse thermal expansion and magnetostriction are highly anisotropic. The resulting volume strain is very small, indicating that the Pr valence remains nearly constant as a function of magnetic field. We conclude that the Fermi liquid state around $B_{\mathrm{c}}$ forms through a very little change in c-f hybridization. This result is in sharp contrast to Ce- and Yb-based Kramers Kondo lattices which show significantly larger volume strains due to the high sensitivity of the Kondo temperature to hydrostatic pressure.


title: A and F stars as probes of outer Galactic disc kinematics

abstract: Previous studies of the rotation law in the outer Galactic disc have mainly used gas tracers or clump giants. Here, we explore A and F stars as alternatives: these provide a much denser sampling in the outer disc than gas tracers and have experienced significantly less velocity scattering than older clump giants. This first investigation confirms the suitability of A stars in this role. Our work is based on spectroscopy of $\sim$ 1300 photometrically-selected stars in the red calcium-triplet region, chosen to mitigate against the effects of interstellar extinction. The stars are located in two low Galactic latitude sightlines, at longitudes $\ell = 118^{\circ}$, sampling strong Galactic rotation shear, and $\ell = 178^{\circ}$, near the Anticentre. With the use of Markov Chain Monte Carlo parameter fitting, stellar parameters and radial velocities are measured, and distances computed. The obtained trend of radial velocity with distance is inconsistent with existing flat or slowly rising rotation laws from gas tracers (Brand & Blitz 1993; Reid et al. 2014). Instead, our results fit in with those obtained by Huang et al. (2016) from disc clump giants that favoured rising circular speeds. An alternative interpretation in terms of spiral arm perturbation is not straight forward. We assess the role that undetected binaries in the sample and distance error may have in introducing bias, and show that the former is a minor factor. The random errors in our trend of circular velocity are within $\pm 5$ km s$^{-1}$.


title: Opposing effects of stacking faults and antisite domain boundaries on the conduction band edge in kesterite quaternary semiconductors

abstract: We investigated stability and the electronic structure of extended defects including anti-site domain boundaries and stacking faults in the kesterite-structured semiconductors, Cu$_2$ZnSnS$_4$ (CZTS) and Cu$_2$ZnSnSe$_4$ (CZTSe). Our hybrid density functional theory calculations show that stacking faults in CZTS and CZTSe induce a higher conduction band edge than the bulk counterparts, and thus the stacking faults act as electron barriers. Antisite domain boundaries, however, accumulate electrons as the conduction band edge is reduced in energy, having an opposite role. An Ising model was constructed to account for the stability of stacking faults, which shows the nearest neighbour interaction is stronger in the case of the selenide.


title: The effect of metallicity on the atmospheres of exoplanets with fully coupled 3D hydrodynamics, equilibrium chemistry, and radiative transfer

abstract: In this work we have performed a series of simulations of the atmosphere of GJ~1214b assuming different metallicities using the Met Office Unified Model (UM). The UM is a general circulation model (GCM) that solves the deep, non-hydrostatic equations of motion and uses a flexible and accurate radiative transfer scheme, based on the two-stream and correlated-$k$ approximations, to calculate the heating rates. In this work we consistently couple a well-tested Gibbs energy minimisation scheme to solve for the chemical equilibrium abundances locally in each grid cell for a general set of elemental abundances, further improving the flexibility and accuracy of the model. As the metallicity of the atmosphere is increased we find significant changes in the dynamical and thermal structure, with subsequent implications for the simulated phase curve. The trends that we find are qualitatively consistent with previous works, though with quantitative differences. We investigate in detail the effect of increasing the metallicity by splitting the mechanism into constituents, involving the mean molecular weight, the heat capacity and the opacities. We find the opacity effect to be the dominant mechanism in altering the circulation and thermal structure. This result highlights the importance of accurately computing the opacities and radiative transfer in 3D GCMs.


title: Grinberg-Kazhdan theorem and Newton groupoids

abstract: We first prove the Grinberg-Kazhdan formal arc theorem without any assumptions on the characteristic. This part of the article is equivalent to arXiv:math-AG/0203263. Then we try to clarify the geometric ideas behind the proof by introducing the notion of Newton groupoid (which is related to Newton's method for finding roots). Newton groupoids are certain groupoids in the category of schemes associated to any generically etale morphism from a locally complete intersection to a smooth variety.


title: On the exact maximum induced density of almost all graphs and their inducibility

abstract: Let $H$ be a graph on $h$ vertices. The number of induced copies of $H$ in a graph $G$ is denoted by $i_H(G)$. Let $i_H(n)$ denote the maximum of $i_H(G)$ taken over all graphs $G$ with $n$ vertices. Let $f(n,h) = \Pi_{i}^h a_i$ where $\sum_{i=1}^h a_i = n$ and the $a_i$ are as equal as possible. Let $g(n,h) = f(n,h) + \sum_{i=1}^h g(a_i,h)$. It is proved that for almost all graphs $H$ on $h$ vertices it holds that $i_H(n)=g(n,h)$ for all $n \le 2^{\sqrt{h}}$. More precisely, we define an explicit graph property ${\cal P}_h$ which, when satisfied by $H$, guarantees that $i_H(n)=g(n,h)$ for all $n \le 2^{\sqrt{h}}$. It is proved, in particular, that a random graph on $h$ vertices satisfies ${\cal P}_h$ with probability $1-o_h(1)$. Furthermore, all extremal $n$-vertex graphs yielding $i_H(n)$ in the aforementioned range are determined. We also prove a stability result. For $H \in {\cal P}_h$ and a graph $G$ with $n \le 2^{\sqrt{h}}$ vertices satisfying $i_H(G) \ge f(n,h)$, it must be that $G$ is obtained from a balanced blowup of $H$ by adding some edges inside the blowup parts. The {\em inducibility} of $H$ is $i_H = \lim_{n \rightarrow \infty} i_H(n)/\binom{n}{h}$. It is known that $i_H \ge h!/(h^h-h)$ for all graphs $H$ and that a random graph $H$ satisfies almost surely that $i_H \le h^{3\log h}h!/(h^h-h)$. We improve upon this upper bound almost matching the lower bound. It is shown that a graph $H$ which satisfies ${\cal P}_h$ has $i_H =(1+O(h^{-h^{1/3}}))h!/(h^h-h)$.


title: Impact Assessment of Hypothesized Cyberattacks on Interconnected Bulk Power Systems

abstract: The first-ever Ukraine cyberattack on power grid has proven its devastation by hacking into their critical cyber assets. With administrative privileges accessing substation networks/local control centers, one intelligent way of coordinated cyberattacks is to execute a series of disruptive switching executions on multiple substations using compromised supervisory control and data acquisition (SCADA) systems. These actions can cause significant impacts to an interconnected power grid. Unlike the previous power blackouts, such high-impact initiating events can aggravate operating conditions, initiating instability that may lead to system-wide cascading failure. A systemic evaluation of "nightmare" scenarios is highly desirable for asset owners to manage and prioritize the maintenance and investment in protecting their cyberinfrastructure. This survey paper is a conceptual expansion of real-time monitoring, anomaly detection, impact analyses, and mitigation (RAIM) framework that emphasizes on the resulting impacts, both on steady-state and dynamic aspects of power system stability. Hypothetically, we associate the combinatorial analyses of steady state on substations/components outages and dynamics of the sequential switching orders as part of the permutation. The expanded framework includes (1) critical/noncritical combination verification, (2) cascade confirmation, and (3) combination re-evaluation. This paper ends with a discussion of the open issues for metrics and future design pertaining the impact quantification of cyber-related contingencies.


title: Orbital Period Increase in ES Ceti

abstract: We report a long-term study of the eclipse times in the 10-minute helium binary ES Ceti. The binary period increases rapidly, with P/P-dot = 6.2x10^6 yr. This is consistent with the assumption that gravitational radiation (GR) drives the mass transfer, and appears to be the first dynamical evidence that GR is indeed the driver of evolution in this class of very old cataclysmic variables -- the AM Canum Venaticorum stars.


title: Information Bottleneck on General Alphabets

abstract: We prove rigorously a source coding theorem that can probably be considered folklore, a generalization to arbitrary alphabets of a problem motivated by the Information Bottleneck method. For general random variables $(Y, X)$, we show essentially that for some $n \in \mathbb{N}$, a function $f$ with rate limit $\log|f| \le nR$ and $I(Y^n; f(X^n)) \ge nS$ exists if and only if there is a random variable $U$ such that the Markov chain $Y - X - U$ holds, $I(U; X) \le R$ and $I(U; Y) \ge S$. The proof relies on the well established discrete case and showcases a technique for lifting discrete coding theorems to arbitrary alphabets.


title: Spot the Difference by Object Detection

abstract: In this paper, we propose a simple yet effective solution to a change detection task that detects the difference between two images, which we call "spot the difference". Our approach uses CNN-based object detection by stacking two aligned images as input and considering the differences between the two images as objects to detect. An early-merging architecture is used as the backbone network. Our method is accurate, fast and robust while using very cheap annotation. We verify the proposed method on the task of change detection between the digital design and its photographic image of a book. Compared to verification based methods, our object detection based method outperforms other methods by a large margin and gives extra information of location. We compress the network and achieve 24 times acceleration while keeping the accuracy. Besides, as we synthesize the training data for detection using weakly labeled images, our method does not need expensive bounding box annotation.


title: Monolithic quantum-dot distributed feedback laser array on silicon

abstract: Electrically-pumped lasers directly grown on silicon are key devices interfacing silicon microelectronics and photonics. We report here, for the first time, an electrically-pumped, room-temperature, continuous-wave (CW) and single-mode distributed feedback (DFB) laser array fabricated in InAs/GaAs quantum-dot (QD) gain material epitaxially grown on silicon. CW threshold currents as low as 12 mA and single-mode side mode suppression ratios (SMSRs) as high as 50 dB have been achieved from individual devices in the array. The laser array, compatible with state-of-the-art coarse wavelength division multiplexing (CWDM) systems, has a well-aligned channel spacing of 20 0.2 nm and exhibits a record wavelength coverage range of 100 nm, the full span of the O-band. These results indicate that, for the first time, the performance of lasers epitaxially grown on silicon is elevated to a point approaching real-world CWDM applications, demonstrating the great potential of this technology.


title: Low-depth circuit ansatz for preparing correlated fermionic states on a quantum computer

abstract: Quantum simulations are bound to be one of the main applications of near-term quantum computers. Quantum chemistry and condensed matter physics are expected to benefit from these technological developments. Several quantum simulation methods are known to prepare a state on a quantum computer and measure the desired observables. The most resource economic procedure is the variational quantum eigensolver (VQE), which has traditionally employed unitary coupled cluster as the ansatz to approximate ground states of many-body fermionic Hamiltonians. A significant caveat of the method is that the initial state of the procedure is a single reference product state with no entanglement extracted from a classical Hartree-Fock calculation. In this work, we propose to improve the method by initializing the algorithm with a more general fermionic Gaussian state, an idea borrowed from the field of nuclear physics. We show how this Gaussian reference state can be prepared with a linear-depth circuit of quantum matchgates. By augmenting the set of available gates with nearest-neighbor phase coupling, we generate a low-depth circuit ansatz that can accurately prepare the ground state of correlated fermionic systems. This extends the range of applicability of the VQE to systems with strong pairing correlations such as superconductors, atomic nuclei, and topological materials.


title: Computational complexity lower bounds of certain discrete Radon transform approximations

abstract: For the computational model where only additions are allowed, the $\Omega(n^2\log n)$ lower bound on operations count with respect to image size $n\times n$ is obtained for two types of the discrete Radon transform implementations: the fast Hough transform and a generic strip pattern class which includes the classical Hough transform, implying the fast Hough transform algorithm asymptotic optimality. The proofs are based on a specific result from the boolean circuits complexity theory and are generalized for the case of boolean $\vee$ binary operation.


title: Gaps between prime numbers and tensor rank of multiplication in finite fields

abstract: We present effective upper bounds on the symmetric bilinear complexity of multiplication in extensions of a base finite field Fp2 of prime square order, obtained by combining estimates on gaps between prime numbers together with an optimal construction of auxiliary divisors for multiplication algorithms by evaluation-interpolation on curves. Most of this material dates back to a 2011 unpublished work of the author, but it still provides the best results on this topic at the present time. Then a few updates are given in order to take recent developments into account, including comparison with a similar work of Ballet and Zykin, generalization to classical bilinear complexity over Fp, and to short multiplication of polynomials, as well as a discussion of open questions on gaps between prime numbers or more generally values of certain arithmetic functions.


title: A New HDG Method for Dirichlet Boundary Control of Convection Diffusion PDEs II: Low Regularity

abstract: In the first part of this work, we analyzed a Dirichlet boundary control problem for an elliptic convection diffusion PDE and proposed a new hybridizable discontinuous Galerkin (HDG) method to approximate the solution. For the case of a 2D polygonal domain, we also proved an optimal superlinear convergence rate for the control under certain assumptions on the domain and on the target state. In this work, we revisit the convergence analysis without these assumptions; in this case, the solution can have low regularity and we use a different analysis approach. We again prove an optimal convergence rate for the control, and present numerical results to illustrate the convergence theory.


title: Full configuration interaction quantum Monte Carlo benchmark and multireference coupled cluster studies of tetramethyleneethane diradical

abstract: We have performed a FCI-quality benchmark calculation for the tetramethyleneethane molecule in cc-pVTZ basis set employing a subset of CASPT2(6,6) natural orbitals for the FCIQMC calculation. The results are in an excellent agreement with the previous large scale diffusion Monte Carlo calculations by Pozun et al. and available experimental results. Our computations verified that there is a maximum on PES of the ground singlet state ($^1\text{A}$) $45^{\circ}$ torsional angle and the corresponding vertical singlet-triplet energy gap is $0.01$ eV. We have employed this benchmark for assessment of the accuracy of MkCCSDT and DMRG-tailored CCSD (TCCSD) methods. MR MkCCSDT with CAS(2,2) model space, though giving good values for the singlet-triplet energy gap, is not able to properly describe the shape of the multireference singlet PES. Similarly, DMRG(24,25) is not able to correctly capture the shape of the singlet surface, due to the missing dynamic correlation. On the other hand, the DMRG-tailored CCSD method describes the shape of the ground singlet state with an excellent accuracy, but for the correct ordering requires computation of the zero-spin-projection component of the triplet state ($^3\text{B}_1$).


title: Polynomial-based rotation invariant features

abstract: One of basic difficulties of machine learning is handling unknown rotations of objects, for example in image recognition. A related problem is evaluation of similarity of shapes, for example of two chemical molecules, for which direct approach requires costly pairwise rotation alignment and comparison. Rotation invariants are useful tools for such purposes, allowing to extract features describing shape up to rotation, which can be used for example to search for similar rotated patterns, or fast evaluation of similarity of shapes e.g. for virtual screening, or machine learning including features directly describing shape. A standard approach are rotationally invariant cylindrical or spherical harmonics, which can be seen as based on polynomials on sphere, however, they provide very few invariants - only one per degree of polynomial. There will be discussed a general approach to construct arbitrarily large sets of rotation invariants of polynomials, for degree $D$ in $\mathbb{R}^n$ up to $O(n^D)$ independent invariants instead of $O(D)$ offered by standard approaches, possibly also a complete set: providing not only necessary, but also sufficient condition for differing only by rotation (and reflectional symmetry).


title: Slowing Down Top Trees for Better Worst-Case Bounds

abstract: We consider the top tree compression scheme introduced by Bille et al. [ICALP 2013] and construct an infinite family of trees on $n$ nodes labeled from an alphabet of size $\sigma$, for which the size of the top DAG is $\Theta(\frac{n}{\log_\sigma n}\log\log_\sigma n)$. Our construction matches a previously known upper bound and exhibits a weakness of this scheme, as the information-theoretic lower bound is $\Omega(\frac{n}{\log_\sigma n})$. This settles an open problem stated by Lohrey et al. [arXiv 2017], who designed a more involved version achieving the lower bound. We show that this can be also guaranteed by a very minor modification of the original scheme: informally, one only needs to ensure that different parts of the tree are not compressed too quickly. Arguably, our version is more uniform, and in particular, the compression procedure is oblivious to the value of $\sigma$.


title: A strictly commutative model for the cochain algebra of a space

abstract: The commutative differential graded algebra $A_{\mathrm{PL}}(X)$ of polynomial forms on a simplicial set $X$ is a crucial tool in rational homotopy theory. In this note, we construct an integral version $A^{\mathcal{I}}(X)$ of $A_{\mathrm{PL}}(X)$. Our approach uses diagrams of chain complexes indexed by the category of finite sets and injections $\mathcal{I}$ to model $E_{\infty}$ differential graded algebras by strictly commutative objects, called commutative $\mathcal{I}$-dgas. We define a functor $A^{\mathcal{I}}$ from simplicial sets to commutative $\mathcal{I}$-dgas and show that it is a commutative lift of the usual cochain algebra functor. In particular, it gives rise to a new construction of the $E_{\infty}$ dga of cochains. The functor $A^{\mathcal{I}}$ shares many properties of $A_{\mathrm{PL}}$, and can be viewed as a generalization of $A_{\mathrm{PL}}$ that works over arbitrary commutative ground rings. Working over the integers, a theorem by Mandell implies that $A^{\mathcal{I}}(X)$ determines the homotopy type of $X$ when $X$ is a nilpotent space of finite type.


title: Intrinsic Gaussian processes on complex constrained domains

abstract: We propose a class of intrinsic Gaussian processes (in-GPs) for interpolation, regression and classification on manifolds with a primary focus on complex constrained domains or irregular shaped spaces arising as subsets or submanifolds of R, R2, R3 and beyond. For example, in-GPs can accommodate spatial domains arising as complex subsets of Euclidean space. in-GPs respect the potentially complex boundary or interior conditions as well as the intrinsic geometry of the spaces. The key novelty of the proposed approach is to utilise the relationship between heat kernels and the transition density of Brownian motion on manifolds for constructing and approximating valid and computationally feasible covariance kernels. This enables in-GPs to be practically applied in great generality, while existing approaches for smoothing on constrained domains are limited to simple special cases. The broad utilities of the in-GP approach is illustrated through simulation studies and data examples.


title: Construction of Algorithms for Parallel Addition

abstract: An algebraic number $\beta \in \mathbb{C}$ with no conjugate of modulus 1 can serve as the base of a numeration system $(\beta, \mathcal{A})$ with parallel addition, i.e., the sum of two operands represented in base $\beta$ with digits from $\mathcal{A}$ is calculated in constant time, irrespective of the length of the operands. In order to allow parallel addition, sufficient level of redundancy must be given to the alphabet $\mathcal{A}$. The complexity of parallel addition algorithm depends heavily on the size $\#\mathcal{A}$ of the alphabet: the bigger alphabet is considered, the lower complexity of the parallel addition algorithm may be reached, and vice versa. Here we aim to find parallel addition algorithms on alphabets of the minimal possible size, for a given base. As the complexity of these algorithms becomes quite huge in general, we introduce a so-called Extending Window Method (EWM) - in fact an algorithm to construct parallel addition algorithms. This method can be applied on bases $\beta$ which are expanding algebraic integers, i.e., $\beta$ whose all conjugates are greater than 1 in modulus. Convergence of the EWM is not guaranteed, nevertheless, we have developed tools for revealing non-convergence, and there is a number of successful applications. Firstly, the EWM provides the same parallel addition algorithms as were previously introduced by A. Avizienis, C.Y. Chow & J.E. Robertson or B. Parhami for integer bases and alphabets. Then, by applying the EWM on selected complex bases $\beta$ with non-integer alphabets $\mathcal{A} \subset \mathbb{Z}[\beta]$, we obtain new results - parallel addition algorithms on alphabets of minimal possible size, which could not be found so far (manually). The EWM is helpful also in the case of block parallel addition.


title: The Effects of Barycentric and Asymmetric Transverse Velocities on Eclipse and Transit Times

abstract: It has long been recognized that the finite speed of light can affect the observed time of an event. For example, as a source moves radially toward or away from an observer, the path length and therefore the light travel time to the observer decreases or increases, causing the event to appear earlier or later than otherwise expected, respectively. This light travel time effect (LTTE) has been applied to transits and eclipses for a variety of purposes, including studies of eclipse timing variations (ETVs) and transit timing variations (TTVs) that reveal the presence of additional bodies in the system. Here we highlight another non-relativistic effect on eclipse or transit times arising from the finite speed of light---caused by an asymmetry in the transverse velocity of the two eclipsing objects, relative to the observer. This asymmetry can be due to a non-unity mass ratio or to the presence of external barycentric motion. Although usually constant, this barycentric and asymmetric transverse velocities (BATV) effect can vary between sequential eclipses if either the path length between the two objects or the barycentric transverse velocity varies in time. We discuss this BATV effect and estimate its magnitude for both time-dependent and time-independent cases. For the time-dependent cases, we consider binaries that experience a change in orbital inclination, eccentric systems with and without apsidal motion, and hierarchical triple systems. We also consider the time-independent case which, by affecting the primary and secondary eclipses differently, can influence the inferred system parameters, such as the orbital eccentricity.


title: Revealing the nature of excitons in liquid exfoliated monolayer WS2

abstract: Transition metal dichalcogenides hold promise for applications in novel optoelectronic devices. There is therefore a need for materials that can be obtained in large quantities and with well understood optical properties. In this report, we present a thorough photoluminescence (PL) investigations of monolayer tungsten disulphide obtained via liquid phase exfoliation. As shown by microscopy studies, the exfoliated nanosheets have dimensions of tens of nanometers and thickness of 2.5 monolayers on average. The monolayer content is about 20%. Our studies show that at low temperature the photoluminescence is dominated by excitons localized on nanosheet edges. As a consequence, the PL is strongly sensitive to environment and exhibits an enhanced splitting in magnetic field. As the temperature is increased, the excitons are thermally excited out of the defect states and the dominant transition is that of the negatively charged exciton. Furthermore, upon excitation with a circularly polarized light, the PL retains a degree of polarization reaching 50% and inherited from the valley polarized photoexcited excitons. The studies of PL dynamics reveal that the PL lifetime is on the order of 10 ps, probably limited by non-radiative processes. Our results underline the potential of liquid exfoliated TMD monolayers in large scale optoelectronic devices.


title: Negative spin-Hall angle and anisotropic spin-orbit torques in epitaxial IrMn

abstract: A spin-torque ferromagnetic resonance study is performed in epitaxial $\mathrm{Fe / Ir_{15}Mn_{85}}$ bilayers with different Fe thicknesses. We measure a negative spin-Hall angle of a few percent in the antiferromagnetic IrMn in contrast to previously reported positive values. A large spin-orbit field with Rashba symmetry opposing the Oersted field is also present. Magnitudes of measured spin-orbit torques depend on the crystallographic direction of current and are correlated with the exchange bias direction set during growth. We suggest that the uncompensated moments at the Fe / IrMn interface are responsible for the observed anisotropy. Our findings highlight the importance of crystalline and magnetic structures for the spin-Hall effect in antiferromagnets.


title: Detection of additional Be+sdO systems from IUE spectroscopy

abstract: There is growing evidence that some Be stars were spun up through mass transfer in a close binary system, leaving the former mass donor star as a hot, stripped-down object. There are five known cases of Be stars with hot subdwarf (sdO) companions that were discovered through International Ultraviolet Explorer (IUE) spectroscopy. Here we expand the search for Be+sdO candidates using archival FUV spectra from IUE. We collected IUE spectra for 264 stars and formed cross-correlation functions (CCFs) with a model spectrum for a hot subdwarf. Twelve new candidate Be+sdO systems were found, and eight of these display radial velocity variations associated with orbital motion. The new plus known Be+sdO systems have Be stars with spectral subtypes of B0 to B3, and the lack of later-type systems is surprising given the large number of cooler B-stars in our sample. We discuss explanations for the observed number and spectral type distribution of the Be+sdO systems, and we argue that there are probably many Be systems with stripped companions that are too faint for detection through our analysis.


title: Remnant radio-loud AGN in the Herschel-ATLAS field

abstract: Only a small fraction of observed Active Galactic Nuclei display large-scale radio emission associated with jets, yet these radio-loud AGN have become increasingly important in models of galaxy evolution. In determining the dynamics and energetics of the radio sources over cosmic time, a key question concerns what happens when their jets switch off. The resulting `remnant' radio-loud AGN have been surprisingly evasive in past radio surveys, and therefore statistical information on the population of radio-loud AGN in their dying phase is limited. In this paper, with the recent developments of LOFAR and the VLA, we are able to provide a systematically selected sample of remnant radio-loud AGN in the Herschel-ATLAS field. Using a simple core-detection method, we constrain the upper limit on the fraction of remnants in our radio-loud AGN sample to 9 per cent, implying that the extended lobe emission fades rapidly once the core/jets turn off. We also find that our remnant sample has a wide range of spectral indices ($-1.5\leqslant \alpha^{1400}_{150}\leqslant -0.5$), confirming that the lobes of some remnants may possess flat spectra at low frequencies just as active sources do. We suggest that, even with the unprecedented sensitivity of LOFAR, our sample may still only contain the youngest of the remnant population.


title: Importance of charge capture in inter-phase regions during readout of charge-coupled devices

abstract: The current understanding of charge transfer dynamics in Charge-Coupled Devices (CCDs) is that charge is moved so quickly from one phase to the next in a clocking sequence and with a density so low that trapping of charge in the inter-phase regions is negligible. However, new simulation capabilities developed at the Centre for Electronic Imaging, that includes direct input of electron density simulations, has made it possible to investigate this assumption further. As part of the radiation testing campaign of the Euclid CCD273 devices, data has been obtained using the trap pumping method, that can be used to identify and characterise single defects CCDs. Combining this data with simulations, we find that trapping during the transfer of charge between phases is indeed necessary in order to explain the results of the data analysis. This result could influence not only trap pumping theory and how trap pumping should be performed, but also how a radiation damaged CCD is read out in the most optimal way.


title: Theoretical links between universal and Bayesian compressed sensing algorithms

abstract: Quantized maximum a posteriori (Q-MAP) is a recently-proposed Bayesian compressed sensing algorithm that, given the source distribution, recovers $X^n$ from its linear measurements $Y^m=AX^n$, where $A\in R^{m\times n}$ denotes the known measurement matrix. On the other hand, Lagrangian minimum entropy pursuit (L-MEP) is a universal compressed sensing algorithm that aims at recovering $X^n$ from its linear measurements $Y^m=AX^n$, without having access to the source distribution. Both Q-MAP and L-MEP provably achieve the minimum required sampling rates, in noiseless cases where such fundamental limits are known. L-MEP is based on minimizing a cost function that consists of a linear combination of the conditional empirical entropy of a potential reconstruction vector and its corresponding measurement error. In this paper, using a first-order linear approximation of the conditional empirical entropy function, L-MEP is connected with Q-MAP. The established connection between L-MEP and Q-MAP leads to variants of Q-MAP which have the same asymptotic performance as Q-MAP in terms of their required sampling rates. Moreover, these variants suggest that Q-MAP is robust to small error in estimating the source distribution. This robustness is theoretically proven and the effect of a non-vanishing estimation error on the required sampling rate is characterized.


title: Real space pairwise electrostatic summation in a uniform neutralising background

abstract: Evaluating the total energy of an extended distribution of point charges, which interact through the Coulomb potential, is central to the study of condensed matter. With near ubiquity, the summation required is carried out using Ewald's method, which splits the problem into two separately convergent sums; one in real space and the other in reciprocal space. Density functional based electronic structure methods require the evaluation of the ion-ion repulsive energy, neutralised by a uniform background charge. Here a purely real-space approach is described. It is straightforward to implement, computationally efficient and offers linear scaling. When applied to the evaluation of the electrostatic energy of neutral ionic crystals, it is shown to be closely related to Wolf's method.


title: Universality and Thouless energy in the supersymmetric Sachdev-Ye-Kitaev Model

abstract: We investigate the supersymmetric Sachdev-Ye-Kitaev (SYK) model, $N$ Majorana fermions with infinite range interactions in $0+1$ dimensions. We have found that, close to the ground state $E \approx 0$, discrete symmetries alter qualitatively the spectral properties with respect to the non-supersymmetric SYK model. The average spectral density at finite $N$, which we compute analytically and numerically, grows exponentially with $N$ for $E \approx 0$. However the chiral condensate, which is normalized with respect the total number of eigenvalues, vanishes in the thermodynamic limit. Slightly above $E \approx 0$, the spectral density grows exponential with the energy. Deep in the quantum regime, corresponding to the first $O(N)$ eigenvalues, the average spectral density is universal and well described by random matrix ensembles with chiral and superconducting discrete symmetries. The dynamics for $E \approx 0$ is investigated by level fluctuations. Also in this case we find excellent agreement with the prediction of chiral and superconducting random matrix ensembles for eigenvalues separations smaller than the Thouless energy, which seems to scale linearly with $N$. Deviations beyond the Thouless energy, which describes how ergodicity is approached, are universality characterized by a quadratic growth of the number variance. In the time domain, we have found analytically that the spectral form factor $g(t)$, obtained from the connected two-level correlation function of the unfolded spectrum, decays as $1/t^2$ for times shorter but comparable to the Thouless time with $g(0)$ related to the coefficient of the quadratic growth of the number variance. Our results provide further support that quantum black holes are ergodic and therefore can be classified by random matrix theory.


title: Randomized Linear Algebra Approaches to Estimate the Von Neumann Entropy of Density Matrices

abstract: Thevon Neumann entropy, named after John von Neumann, is an extension of the classical concept of entropy to the field of quantum mechanics. From a numerical perspective, von Neumann entropy can be computed simply by computing all eigenvalues of a density matrix, an operation that could be prohibitively expensive for large-scale density matrices. We present and analyze three randomized algorithms to approximate von Neumann entropy of {real} density matrices: our algorithms leverage recent developments in the Randomized Numerical Linear Algebra (RandNLA) literature, such as randomized trace estimators, provable bounds for the power method, and the use of random projections to approximate the eigenvalues of a matrix. All three algorithms come with provable accuracy guarantees and our experimental evaluations support our theoretical findings showing considerable speedup with small loss in accuracy.


title: EXPSPACE-hardness of behavioural equivalences of succinct one-counter nets

abstract: We note that the remarkable EXPSPACE-hardness result in [G\"oller, Haase, Ouaknine, Worrell, ICALP 2010] ([GHOW10] for short) allows us to answer an open complexity question for simulation preorder of succinct one counter nets (i.e., one counter automata with no zero tests where counter increments and decrements are integers written in binary). This problem, as well as bisimulation equivalence, turn out to be EXPSPACE-complete. The technique of [GHOW10] was referred to by Hunter [RP 2015] for deriving EXPSPACE-hardness of reachability games on succinct one-counter nets. We first give a direct self-contained EXPSPACE-hardness proof for such reachability games (by adjusting a known PSPACE-hardness proof for emptiness of alternating finite automata with one-letter alphabet); then we reduce reachability games to (bi)simulation games by using a standard "defender-choice" technique.


title: Forcing large tight components in 3-graphs

abstract: Any $n$-vertex $3$-graph with minimum codegree at least $\lfloor n/3\rfloor$ must have a spanning tight component, but immediately below this threshold it is possible for no tight component to span more than $\lceil 2n/3\rceil$ vertices. Motivated by this observation, we ask which codegree forces a tight component of at least any given size. The corresponding function seems to have infinitely many discontinuities, but we provide upper and lower bounds, which asymptotically converge as the function nears the origin.


title: LIME: Live Intrinsic Material Estimation

abstract: We present the first end to end approach for real time material estimation for general object shapes with uniform material that only requires a single color image as input. In addition to Lambertian surface properties, our approach fully automatically computes the specular albedo, material shininess, and a foreground segmentation. We tackle this challenging and ill posed inverse rendering problem using recent advances in image to image translation techniques based on deep convolutional encoder decoder architectures. The underlying core representations of our approach are specular shading, diffuse shading and mirror images, which allow to learn the effective and accurate separation of diffuse and specular albedo. In addition, we propose a novel highly efficient perceptual rendering loss that mimics real world image formation and obtains intermediate results even during run time. The estimation of material parameters at real time frame rates enables exciting mixed reality applications, such as seamless illumination consistent integration of virtual objects into real world scenes, and virtual material cloning. We demonstrate our approach in a live setup, compare it to the state of the art, and demonstrate its effectiveness through quantitative and qualitative evaluation.


title: Liquid crystal elastomer coatings with programmed response of surface profile

abstract: Stimuli-responsive liquid crystal elastomers (LCEs) with a strong coupling of orientational molecular order and rubber-like elasticity, show a great potential as working elements in soft robotics, sensing, transport and propulsion systems. We demonstrate a dynamic thermal control of the surface topography of LCE coatings achieved through pre-designed patterns of in-plane molecular orientation. These patterns determine whether the LCE coating develops elevations, depressions, or in-plane deformations. The deterministic dependence of the out-of-plane dynamic surface profile on the in-plane orientational pattern is explained by activation forces. These forces are caused by two factors: (i) stretching-contraction of the polymer networks driven by temperature; (ii) spatially varying orientation of the LCE. The activation force concept brings the responsive LCEs into the domain of active matter. The demonstrated relationship can be used to design programmable coatings with functionalities that mimic biological tissues such as skin.


title: Analytic vortex solutions in generalized models of the Maxwell-Higgs type

abstract: This work deals with the presence of analytical vortex configurations in generalized models of the Maxwell-Higgs type in the three-dimensional spacetime. We implement a procedure that allows to decouple the first order equations, which we use to solve the model analytically. The approach is exemplified with three distinct models that show the robustness of the construction. In the third model, one finds analytical solutions that exhibit interesting compact vortex behavior.


title: Recent Advances in Recurrent Neural Networks

abstract: Recurrent neural networks (RNNs) are capable of learning features and long term dependencies from sequential and time-series data. The RNNs have a stack of non-linear units where at least one connection between units forms a directed cycle. A well-trained RNN can model any dynamical system; however, training RNNs is mostly plagued by issues in learning long-term dependencies. In this paper, we present a survey on RNNs and several new advances for newcomers and professionals in the field. The fundamentals and recent advances are explained and the research challenges are introduced.


title: Exact Calculation of the Mean-Square Error in the Method of Expansion of Iterated Ito Stochastic integrals Based on Generalized Multiple Fourier Series

abstract: The article is devoted to the developement of the method of expansion and mean-square approximation of iterated Ito stochastic integrals based on generalized multiple Fourier series converging in the sense of norm in the space $L_2([t, T]^k)$ ($k$ is the multiplicity of the iterated Ito stochastic integral). We obtain the exact and approximate expressions for the mean-square approximation error of iterated Ito stochastic integrals of multiplicity $k$ ($k\in\mathbb{N}$) from the stochastic Taylor-Ito expansion. As a result, we do not need to use redundant terms of expansions of iterated Ito stochastic integrals that complicate the numerical methods for Ito stochastic differential equations. Moreover, we proved the convergence with propability 1 for the method of expansion of iterated Ito stochastic integrals based on generalized multiple Fourier series for the cases of multiple Fourier-Legendre series and multiple trigonometric Fourier series. Mean-square approximation of iterated Stratonovich stochastic integrals is also considered in the article. The results of the article can be applied to the high-order strong numerical methods for Ito stochastic differential equations as well as for non-commutative semilinear stochastic partial differential equations with multiplicative trace class noise.


title: Cooperative Training of Deep Aggregation Networks for RGB-D Action Recognition

abstract: A novel deep neural network training paradigm that exploits the conjoint information in multiple heterogeneous sources is proposed. Specifically, in a RGB-D based action recognition task, it cooperatively trains a single convolutional neural network (named c-ConvNet) on both RGB visual features and depth features, and deeply aggregates the two kinds of features for action recognition. Differently from the conventional ConvNet that learns the deep separable features for homogeneous modality-based classification with only one softmax loss function, the c-ConvNet enhances the discriminative power of the deeply learned features and weakens the undesired modality discrepancy by jointly optimizing a ranking loss and a softmax loss for both homogeneous and heterogeneous modalities. The ranking loss consists of intra-modality and cross-modality triplet losses, and it reduces both the intra-modality and cross-modality feature variations. Furthermore, the correlations between RGB and depth data are embedded in the c-ConvNet, and can be retrieved by either of the modalities and contribute to the recognition in the case even only one of the modalities is available. The proposed method was extensively evaluated on two large RGB-D action recognition datasets, ChaLearn LAP IsoGD and NTU RGB+D datasets, and one small dataset, SYSU 3D HOI, and achieved state-of-the-art results.


title: High Performance Quantum Modular Multipliers

abstract: We present a novel set of reversible modular multipliers applicable to quantum computing, derived from three classical techniques: 1) traditional integer division, 2) Montgomery residue arithmetic, and 3) Barrett reduction. Each multiplier computes an exact result for all binary input values, while maintaining the asymptotic resource complexity of a single (non-modular) integer multiplier. We additionally conduct an empirical resource analysis of our designs in order to determine the total gate count and circuit depth of each fully constructed circuit, with inputs as large as 2048 bits. Our comparative analysis considers both circuit implementations which allow for arbitrary (controlled) rotation gates, as well as those restricted to a typical fault-tolerant gate set.


title: A first integrability result for Miquel dynamics

abstract: Miquel dynamics is a discrete-time dynamical system on the space of square-grid circle patterns. For biperiodic circle patterns with both periods equal to two, we show that the dynamics corresponds to translation on an elliptic curve, thus providing the first integrability result for this dynamics. The main tool is a geometric interpretation of the addition law on the normalization of binodal quartic curves.