title: Performance Analysis of Robust Stable PID Controllers Using Dominant Pole Placement for SOPTD Process Models

abstract: This paper derives new formulations for designing dominant pole placement based proportional-integral-derivative (PID) controllers to handle second order processes with time delays (SOPTD). Previously, similar attempts have been made for pole placement in delay-free systems. The presence of the time delay term manifests itself as a higher order system with variable number of interlaced poles and zeros upon Pade approximation, which makes it difficult to achieve precise pole placement control. We here report the analytical expressions to constrain the closed loop dominant and non-dominant poles at the desired locations in the complex s-plane, using a third order Pade approximation for the delay term. However, invariance of the closed loop performance with different time delay approximation has also been verified using increasing order of Pade, representing a closed to reality higher order delay dynamics. The choice of the nature of non-dominant poles e.g. all being complex, real or a combination of them modifies the characteristic equation and influences the achievable stability regions. The effect of different types of non-dominant poles and the corresponding stability regions are obtained for nine test-bench processes indicating different levels of open-loop damping and lag to delay ratio. Next, we investigate which expression yields a wider stability region in the design parameter space by using Monte Carlo simulations while uniformly sampling a chosen design parameter space. Various time and frequency domain control performance parameters are investigated next, as well as their deviations with uncertain process parameters, using thousands of Monte Carlo simulations, around the robust stable solution for each of the nine test-bench processes.


title: Vector fields on $\mathfrak{osp}_{2m-1|2n}(\mathbb C)$-flag supermanifolds

abstract: We compute the Lie superalgebras of holomorphic vector fields on isotropic flag supermanifolds of maximal type corresponding to the Lie superalgebra $\mathfrak{osp}_{2m-1|2n}(\mathbb C)$.


title: Time Constrained Continuous Subgraph Search over Streaming Graphs

abstract: The growing popularity of dynamic applications such as social networks provides a promising way to detect valuable information in real time. Efficient analysis over high-speed data from dynamic applications is of great significance. Data from these dynamic applications can be easily modeled as streaming graph. In this paper, we study the subgraph (isomorphism) search over streaming graph data that obeys timing order constraints over the occurrence of edges in the stream. We propose a data structure and algorithm to efficiently answer subgraph search and introduce optimizations to greatly reduce the space cost, and propose concurrency management to improve system throughput. Extensive experiments on real network traffic data and synthetic social streaming data confirms the efficiency and effectiveness of our solution.


title: Virtually abelian subgroups of $IA_n(Z/3)$ are abelian

abstract: When studying subgroups of $Out(F_n)$, one often replaces a given subgroup $H$ with one of its finite index subgroups $H_0$ so that virtual properties of $H$ become actual properties of $H_0$. In many cases, the finite index subgroup is $H_0 = H \cap IA_n(Z/3)$. For which properties is this a good choice? Our main theorem states that being abelian is such a property. Namely, every virtually abelian subgroup of $IA_n(Z/3)$ is abelian.


title: Joint Voxel and Coordinate Regression for Accurate 3D Facial Landmark Localization

abstract: 3D face shape is more expressive and viewpoint-consistent than its 2D counterpart. However, 3D facial landmark localization in a single image is challenging due to the ambiguous nature of landmarks under 3D perspective. Existing approaches typically adopt a suboptimal two-step strategy, performing 2D landmark localization followed by depth estimation. In this paper, we propose the Joint Voxel and Coordinate Regression (JVCR) method for 3D facial landmark localization, addressing it more effectively in an end-to-end fashion. First, a compact volumetric representation is proposed to encode the per-voxel likelihood of positions being the 3D landmarks. The dimensionality of such a representation is fixed regardless of the number of target landmarks, so that the curse of dimensionality could be avoided. Then, a stacked hourglass network is adopted to estimate the volumetric representation from coarse to fine, followed by a 3D convolution network that takes the estimated volume as input and regresses 3D coordinates of the face shape. In this way, the 3D structural constraints between landmarks could be learned by the neural network in a more efficient manner. Moreover, the proposed pipeline enables end-to-end training and improves the robustness and accuracy of 3D facial landmark localization. The effectiveness of our approach is validated on the 3DFAW and AFLW2000-3D datasets. Experimental results show that the proposed method achieves state-of-the-art performance in comparison with existing methods.


title: Motivic Intersection Complex of Certain Shimura varieties

abstract: Using a version of weight conservativity we demonstrate that for certain Shimura varieties (including all Shimura three-folds, most Shimura four-folds and the Siegel sixfold) the construction of the motivic intersection complex due to Wildeshaus compares with a motivic weight truncation in the sense of S. Morel. In particular it is defined up to a unique isomorphism, and satisfies the intrinsic characterization for an intermediate extension due to Wildeshaus.


title: A periodic solution of period two of a delay differential equation

abstract: In this paper we prove that the following delay differential equation \[ \frac{d}{dt}x(t)=rx(t)\left(1-\int_{0}^{1}x(t-s)ds\right), \] has a periodic solution of period two for $r>\frac{\pi^{2}}{2}$ (when the steady state, $x=1$, is unstable). In order to find the periodic solution, we study an integrable system of ordinary differential equations, following the idea by Kaplan and Yorke \cite{Kaplan=000026Yorke:1974}. The periodic solution is expressed in terms of the Jacobi elliptic functions.


title: Wavelet Analysis of the Besov Regularity of L\'evy White Noises

abstract: We characterize the local smoothness and the asymptotic growth rate of the L\'evy white noise. We do so by characterizing the weighted Besov spaces in which it is located. We extend known results in two ways. First, we obtain new bounds for the local smoothness via the Blumenthal-Getoor indices of the L\'evy white noise. We also deduce the critical local smoothness when the two indices coincide, which is true for symmetric-alpha-stable, compound Poisson, and symmetric-gamma white noises to name a few. Second, we express the critical asymptotic growth rate in terms of the moment properties of the L\'evy white noise. Previous analyses only provided lower bounds for both the local smoothness and the asymptotic growth rate. Showing the sharpness of these bounds requires us to determine in which Besov spaces a given L\'evy white noise is (almost surely) not. Our methods are based on the wavelet-domain characterization of Besov spaces and precise moment estimates for the wavelet coefficients of the L\'evy white noise.


title: Topological properties of convolutor spaces via the short-time Fourier transform

abstract: We discuss the structural and topological properties of a general class of weighted $L^1$ convolutor spaces. Our theory simultaneously applies to weighted $\mathcal{D}'_{L^1}$ spaces as well as to convolutor spaces of the Gelfand-Shilov spaces $\mathcal{K}\{M_p\}$. In particular, we characterize the sequences of weight functions $(M_p)_{p \in \mathbb{N}}$ for which the space of convolutors of $\mathcal{K}\{M_p\}$ is ultrabornological, thereby generalizing Grothendieck's classical result for the space $\mathcal{O}'_{C}$ of rapidly decreasing distributions. Our methods lead to the first direct proof of the completeness of the space $\mathcal{O}_{C}$ of very slowly increasing smooth functions.


title: Stability, convergence, and sensitivity analysis of the Filament Based Lamellipodium Model and the corresponding FEM

abstract: This paper focuses on the study of the Filament Based Lamellipodium Model (FBLM) and the corresponding Finite Element Method (FEM) from a numerical point of view. We study fundamental numerical properties of the FEM and justify the further use of the FBLM. We exhibit that the FEM satisfies a timestep stability condition that is consistent with the nature of the problem. We propose a particular strategy to automatically adapt the time step of the method. We show that the FEM convergences with respect to the (two-dimensional) space discretization in a series of characteristic and representative experiments. We embed and couple the FBLM with a complex extracellular environment comprised of chemical and haptic components and study their combined time evolution. Under this prism, we study the sensitivity of the FBLM on several of its controlling parameters and discuss their influence in the development of the model.


title: The point equivalence problem for ordinary differential equations of the second order

abstract: We use E. Cartan's method to solve the problem of equivalence of the second order ordinary differential equations with respect to the pseudogroup of point transformations.


abstract: Polarity imbalance of the photospheric magnetic field was studied using synoptic maps of NSO Kitt Peak (1976-2016). Imbalance of positive and negative fluxes was considered for the fields with strength B>50 G in the sunspot zone (5 deg - 40 deg) and for the fields with strength B<50 G at higher latitudes (40 deg - 90 deg). The 22-year periodicity in the imbalance of positive and negative fields was found which maintained itself during four solar cycles. While for the sunspot zone the sign of the imbalance always coincides with the northern hemisphere polarity, for the high latitudes the sign of the imbalance always coincides with the southern hemisphere polarity. Good correspondence of the flux imbalance with the quadrupole moment (g20) of the potential-field source-surface (PFSS) model was observed. The polarity imbalance of the sunspot zone correlates, on one hand, with the asymmetry of the magnetic field of the Sun-as-a-star and, on the other hand, with the sector structure of the interplanetary magnetic field. The obtained results show the close connection of the magnetic fields in active regions with the Sun's polar magnetic field. The weakest fields B<5 G represent quite a special group with the magnetic flux developing in antiphase to the fluxes of the stronger fields.


title: Virtual Breakpoints for x86/64

abstract: Efficient, reliable trapping of execution in a program at the desired location is a linchpin technique for dynamic malware analysis. The progression of debuggers and malware is akin to a game of cat and mouse - each are constantly in a state of trying to thwart one another. At the core of most efficient debuggers today is a combination of virtual machines and traditional binary modification breakpoints (int3). In this paper, we present a design for Virtual Breakpoints. a modification to the x86 MMU which brings breakpoint management into hardware alongside page tables. In this paper we demonstrate the fundamental abstraction failures of current trapping methods, and design a new mechanism from the hardware up. This design incorporates lessons learned from 50 years of virtualization and debugger design to deliver fast, reliable trapping without the pitfalls of traditional binary modification.


title: Multi-Pointer Co-Attention Networks for Recommendation

abstract: Many recent state-of-the-art recommender systems such as D-ATT, TransNet and DeepCoNN exploit reviews for representation learning. This paper proposes a new neural architecture for recommendation with reviews. Our model operates on a multi-hierarchical paradigm and is based on the intuition that not all reviews are created equal, i.e., only a select few are important. The importance, however, should be dynamically inferred depending on the current target. To this end, we propose a review-by-review pointer-based learning scheme that extracts important reviews, subsequently matching them in a word-by-word fashion. This enables not only the most informative reviews to be utilized for prediction but also a deeper word-level interaction. Our pointer-based method operates with a novel gumbel-softmax based pointer mechanism that enables the incorporation of discrete vectors within differentiable neural architectures. Our pointer mechanism is co-attentive in nature, learning pointers which are co-dependent on user-item relationships. Finally, we propose a multi-pointer learning scheme that learns to combine multiple views of interactions between user and item. Overall, we demonstrate the effectiveness of our proposed model via extensive experiments on \textbf{24} benchmark datasets from Amazon and Yelp. Empirical results show that our approach significantly outperforms existing state-of-the-art, with up to 19% and 71% relative improvement when compared to TransNet and DeepCoNN respectively. We study the behavior of our multi-pointer learning mechanism, shedding light on evidence aggregation patterns in review-based recommender systems.


title: End to End Performance Analysis of Relay Cooperative Communication Based on Parked Cars

abstract: Parking lots (PLs) are usually full with cars. If these cars are formed into a self-organizing vehicular network, they can be new kind of road side units (RSUs) in urban area to provide communication data forwarding between mobile terminals nearby and a base station. However cars in PLs can leave at any time, which is neglected in the existing studies. In this paper, we investigate relay cooperative communication based on parked cars in PLs. Taking the impact of the car's leaving behavior into consideration, we derive the expressions of outage probability in a two-hop cooperative communication and its link capacity. Finally, the numerical results show that the impact of a car's arriving time is greater than the impact of the duration the car has parked on outage probability.


title: Borosilicones and viscoelastic silicone rubbers: network liquids and network solids

abstract: Borosilicones (e.g., Silly Putty) have been known for 70 years, but their peculiar behaviors have remain unexplained. In this work, experiment and theory are used to show that they are network liquids---dynamic macromolecules that appear elastic on short timescales but exhibit flow on longer timescales. Each borosilicone is a vast covalent network of silicone polymer chains joined by trifunctional boron crosslinks. At any instant, the borosilicone is a highly-crosslinked elastic material. Because the boron crosslinks are temporary, however, the network evolves with time and the borosilicone exhibits liquid behavior. A simple borosilicone exemplifies a classic transient network model and behaves as a simple (Lodge) elastic fluid. Its measured moduli and viscosities fit those predicted by the transient network model and the Maxwell viscoelastic model: a spring in series with a dashpot, including the observed exponential relaxation processes. When a borosilicone includes permanent crosslinks, however, it no longer behaves as a simple elastic fluid. Its measured moduli and viscosities fit those predicted by the Fractional Maxwell viscoelastic model: a spring in series with a the spring-pot, including the observed slower-than-exponential relaxation processes. Beyond the gelation threshold, a borosilicone becomes a viscoelastic silicone rubber (VSR). With a permanent network that spans the material coupled to a temporary network that also spans the material, the VSR is a network liquid piggybacking on a network solid. The Fractional Zener viscoelastic model: an elastic spring in parallel to the Fractional Maxwell model, accurately predicts the measured moduli of VSRs. The temporary nature of boron crosslinks is due to exchange reactions. Because the mean lifetime of temporary crosslinks is a borosilicone's only significant timescale, it exhibits thermo-rheological simplicity.


title: SEURAT: SPH scheme extended with ultraviolet line radiative transfer

abstract: We present a novel Lyman alpha (Ly$\alpha$) radiative transfer code, SEURAT, where line scatterings are solved adaptively with the resolution of the smoothed particle hydrodynamics (SPH). The radiative transfer method implemented in SEURAT is based on a Monte Carlo algorithm in which the scattering and absorption by dust are also incorporated. We perform standard test calculations to verify the validity of the code; (i) emergent spectra from a static uniform sphere, (ii) emergent spectra from an expanding uniform sphere, and (iii) escape fraction from a dusty slab. Thereby we demonstrate that our code solves the Ly$\alpha$ radiative transfer with sufficient accuracy. We emphasise that SEURAT can treat the transfer of Ly$\alpha$ photons even in highly complex systems that have significantly inhomogeneous density fields. The high adaptivity of SEURAT is desirable to solve the propagation of Ly$\alpha$ photons in the interstellar medium of young star-forming galaxies like Ly$\alpha$ emitters (LAEs). Thus, SEURAT provides a powerful tool to model the emergent spectra of Ly$\alpha$ emission, which can be compared to the observations of LAEs.


title: Zeeman Splitting and Inverted Polarization of Biexciton Emission in Monolayer WS2

abstract: We investigate the magnetic-field-induced splitting of biexcitons in monolayer WS$_2$ using polarization-resolved photoluminescence spectroscopy in out-of-plane magnetic fields up to 30 T. The observed $g$ factor of the biexciton amounts to $-3.89$, closely matching the $g$ factor of the neutral exciton. The biexciton emission shows an inverted circular field-induced polarization upon linearly polarized excitation, i.e. it exhibits preferential emission from the high-energy peak in a magnetic field. This phenomenon is explained by taking into account the configuration of the biexciton constituents in momentum space and their respective energetic behavior in magnetic fields. Our findings reveal the critical role of dark excitons in the composition of this many-body state.


title: Study on Energy Consumption and Coverage of Hierarchical Cooperation of Small Cell Base Stations in Heterogeneous Networks

abstract: The demand for communication services in the era of intelligent terminals is unprecedented and huge. To meet such development, modern wireless communications must provide higher quality services with higher energy efficiency in terms of system capacity and quality of service (QoS), which could be achieved by the high-speed data rate, the wider coverage and the higher band utilization. In this paper, we propose a way to offload users from a macro base station(MBS) with a hierarchical distribution of small cell base stations(SBS). The connection probability is the key indicator of the implementation of the unload operation. Furthermore, we measure the service performance of the system by finding the conditional probability-coverage probability with the certain SNR threshold as the condition, that is, the probability of obtaining the minimum communication quality when the different base stations are connected to the user. Then, user-centered total energy consumption of the system is respectively obtained when the macro base station(MBS) and the small cell base stations(SBS) serve each of the users. The simulation results show that the hierarchical SBS cooperation in heterogeneous networks can provide a higher system total coverage probability for the system with a lower overall system energy consumption than MBS.


title: Volume entropy and information flow in a brain graph

abstract: Entropy is a classical measure to quantify the amount of information or complexity of a system. Various entropy-based measures such as functional and spectral entropies have been proposed in brain network analysis. However, they are less widely used than traditional graph theoretic measures such as global and local efficiencies because either they are not well-defined on a graph or difficult to interpret its biological meaning. In this paper, we propose a new entropy-based graph invariant, called volume entropy. It measures the exponential growth rate of the number of paths in a graph, which is a relevant measure if information flows through the graph forever. We model the information propagation on a graph by the generalized Markov system associated to the weighted edge-transition matrix. We estimate the volume entropy using the stationary equation of the generalized Markov system. A prominent advantage of using the stationary equation is that it assigns certain distribution of weights on the edges of the brain graph, which we call the stationary distribution. The stationary distribution shows the information capacity of edges and the direction of information flow on a brain graph. The simulation results show that the volume entropy distinguishes the underlying graph topology and geometry better than the existing graph measures. In brain imaging data application, the volume entropy of brain graphs was significantly related to healthy normal aging from 20s to 60s. In addition, the stationary distribution of information propagation gives a new insight into the information flow of functional brain graph.


title: Toward a Phenomenology of Computational Thinking in STEM Education

abstract: In this chapter, we argue for an epistemological shift from viewing coding and computational thinking as mastery over computational logic and symbolic forms, to viewing them as a more complex form of experience. Rather than viewing computing as regurgitation and production of a set of axiomatic computational abstractions, we argue that computing and computational thinking, should be viewed as discursive, perspectival, material and embodied experiences, among others. These experiences include, but are not subsumed by, the use and production of computational abstractions. We illustrate what this paradigmatic shift toward a more phenomenological account of computing can mean for teaching and learning STEM in K12 classrooms by presenting a critical review of the literature, as well as by presenting a review of several studies we have conducted in K12 educational settings grounded in this perspective. Our analysis reveals several phenomenological approaches that can be useful for framing computational thinking in K12 STEM classrooms.


title: Nonlinear refractive index of electric field aligned gold nanorods measured with a Hartmann-Shack wavefront aberrometer

abstract: The capability to dynamically control the nonlinear refractive index of plasmonic suspensions may enable innovative nonlinear sensing and signaling nanotechnologies. Here, we experimentally determine the effective nonlinear refractive index for gold nanorods suspended in an index matching oil aligned using electric fields, demonstrating an approach to modulate the nonlinear optical properties of the suspension. The nonlinear optical experiments were carried out using a Hartmann-Shack wavefront aberrometer in a collimated beam configuration with a high repetition rate femtosecond laser. The suspensions were probed at 800 nm, overlapping with the long-axis absorption peak of the nanorods. We find that the effective nonlinear refractive index of the gold nanorods suspension depends linearly on the orientational order parameter, S, which can be understood by a thermally induced nonlinear response. We also show the magnitude of the nonlinear response can be varied by ~60%.


title: Transition to turbulence in shear flows

abstract: Pipe flow and many other shear flows show a transition to turbulence at flow rates for which the laminar profile is stable against infinitesimal perturbations. In this brief review the recent progress in the understanding of this transition will be summarized, with a focus on the linear and nonlinear states that drive the transitions, the extended and localized patterns that appear, and on the spatio-temporal dynamics and their relation to directed percolation.


title: Inverse Uncertainty Quantification using the Modular Bayesian Approach based on Gaussian Process, Part 2: Application to TRACE

abstract: Inverse Uncertainty Quantification (UQ) is a process to quantify the uncertainties in random input parameters while achieving consistency between code simulations and physical observations. In this paper, we performed inverse UQ using an improved modular Bayesian approach based on Gaussian Process (GP) for TRACE physical model parameters using the BWR Full-size Fine-Mesh Bundle Tests (BFBT) benchmark steady-state void fraction data. The model discrepancy is described with a GP emulator. Numerical tests have demonstrated that such treatment of model discrepancy can avoid over-fitting. Furthermore, we constructed a fast-running and accurate GP emulator to replace TRACE full model during Markov Chain Monte Carlo (MCMC) sampling. The computational cost was demonstrated to be reduced by several orders of magnitude. A sequential approach was also developed for efficient test source allocation (TSA) for inverse UQ and validation. This sequential TSA methodology first selects experimental tests for validation that has a full coverage of the test domain to avoid extrapolation of model discrepancy term when evaluated at input setting of tests for inverse UQ. Then it selects tests that tend to reside in the unfilled zones of the test domain for inverse UQ, so that one can extract the most information for posterior probability distributions of calibration parameters using only a relatively small number of tests. This research addresses the "lack of input uncertainty information" issue for TRACE physical input parameters, which was usually ignored or described using expert opinion or user self-assessment in previous work. The resulting posterior probability distributions of TRACE parameters can be used in future uncertainty, sensitivity and validation studies of TRACE code for nuclear reactor system design and safety analysis.


title: Cell-less Communications in 5G Vehicular Networks Based on Vehicle-Installed Access Points

abstract: The development of intelligent transportation systems raises many requirements to the current vehicular networks. For instance, to ensure secure communications between vehicles, low latency, high connectivity and high data rate are required for vehicular networks. To meet such requirements, the 5G communication systems for vehicular networks should be improved accordingly. This article proposes a communication scheme for 5G vehicular networks, in which moving access points are deployed on vehicles to facilitate the access of vehicle users. Moreover, the adjacent vehicle-installed moving access points may cooperatively communicate with the vehicle users by joint transmissions and receptions. In this way, the vehicle users communicate with one or more unspecified cooperative access points in a cell-less manner instead of being associated to a single access point. To manage such cell-less networks, local moving software-defined cloudlets are adopted to perform transmission and scheduling management. Simulation results show that the proposed scheme significantly reduces the latency, while improving the connectivity of the vehicular networks, and can be considered as a research direction for the solution to 5G vehicular networks.


title: Fast hybrid density-functional computations using plane-wave basis sets

abstract: A new, very fast, implementation of the exact (Fock) exchange operator for electronic structure calculations within the plane-wave pseudopotential method is described in detail for both molecular and periodic systems, and carefully validated. Our method combines the recently proposed Adaptively Compressed Exchange approach, to reduce the number of times the exchange is evaluated in the self-consistent loop, with an orbital localization procedure that reduces the number of exchange integrals to be computed at each evaluation and potentially the compute time of each of them. The new implementation, already available in the Quantum ESPRESSO distribution, results in a speedup that is never smaller than 3-4x and that increases with the size of the system, according to various realistic benchmark calculations.


title: An Energy Stable One-Field Fictitious Domain Method for Fluid-Structure Interactions

abstract: In this article, the energy stability of a one-field fictitious domain method is proved and validated by numerical tests in two and three dimensions. The distinguishing feature of this method is that it only solves for one velocity field for the whole fluid-structure domain; the interactions remain decoupled until solving the final linear algebraic equations. To achieve this the finite element procedures are carried out separately on two different meshes for the fluid and solid respectively, and the assembly of the final linear system brings the fluid and solid parts together via an isoparametric interpolation matrix between the two meshes. The weak formulations are introduced in the continuous case and after discretization in time. Then the stability is analyzed through an energy estimate. Finally, numerical examples are presented to validate the energy stability properties.


title: The weakly dependent strong law of large numbers revisited

abstract: We give a short, self-contained, and elementary proof of the strong law of large numbers under a power law decay hypothesis for joint second moments. The result is related to the classical one by Lyons. However, we also provide a rate of convergence. Our proof does not use maximal inequalities and is instead inspired by the method of multiscale large versus small field decompositions in constructive quantum field theory.


title: Reducible subgroups of exceptional algebraic groups

abstract: Let $G$ be a simple algebraic group over an algebraically closed field. A closed subgroup $H$ of $G$ is called $G$-completely reducible ($G$-cr) if, whenever $H$ is contained in a parabolic subgroup $P$ of $G$, it is contained in a Levi factor of $P$. In this paper we complete the classification of connected $G$-cr subgroups when $G$ has exceptional type, by determining the $L_{0}$-irreducible connected reductive subgroups for each simple classical factor $L_{0}$ of a Levi subgroup of $G$. As an illustration, we determine all reducible, $G$-cr semisimple subgroups when $G$ has type $F_4$ and various properties thereof. This work complements results of Lawther, Liebeck, Seitz and Testerman, and is vital in classifying non-$G$-cr reductive subgroups, a project being undertaken by the authors elsewhere.


title: Clean Os(0001) electronic surface states: a first-principle fully relativistic investigation

abstract: We analyze the electronic structure of the Os(0001) surface by means of first-principle calculations based on Fully Relativistic (FR) Density Functional Theory (DFT) and a Projector Augmented-Wave (PAW) approach. We investigate surface states and resonances analyzing their spin-orbit induced energy splitting and their spin polarization. The results are compared with previously studied surfaces Ir(111), Pt(111), and Au(111). We do not find any surface state in the gap similar to the L-gap of the (111) fcc surfaces, but find Rashba split resonances that cross the Fermi level and, as in the recently studied Ir(111) surface, have a characteristic downward dispersion. Moreover, for some selected surface states we study the spin polarization with respect to ${\bf k}_{\parallel}$, the wave-vector parallel to the surface. In some cases, such as the Rashba split resonances, the spin polarization shows a smooth behavior with slow rotations, in others the rotation is faster, due to mixing and anti-crossing of the states.


title: A Finite Soluble Quotient Algorithm

abstract: An algorithm for computing power conjugate presentations for finite soluble quotients of predetermined structure of finitely presented groups is described. Practical aspects of an implementation are discussed.


title: Wasserstein Riemannian Geometry of Positive Definite Matrices

abstract: The Wasserstein distance on multivariate non-degenerate Gaussian densities is a Riemannian distance. After reviewing the properties of the distance and the metric geodesic, we present an explicit form of the Riemannian metrics on positive-definite matrices and compute its tensor form with respect to the trace inner product. The tensor is a matrix which is the solution to a Lyapunov equation. We compute the explicit formula for the Riemannian exponential, the normal coordinates charts and the Riemannian gradient. Finally, the Levi-Civita covariant derivative is computed in matrix form together with the differential equation for the parallel transport. While all computations are given in matrix form, nonetheless we discuss also the use of a special moving frame.


title: Duality and mapping tori in Heegaard Floer homology

abstract: We show that the graph TQFT for Heegaard Floer homology satisfies a strong version of Atiyah's duality axiom for a TQFT. As an application, we compute some Heegaard Floer mixed invariants of 4-dimensional mapping tori in terms of Lefschetz numbers on $HF^+$.


title: Deep Reinforcement Learning for Dynamic Treatment Regimes on Medical Registry Data

abstract: This paper presents the first deep reinforcement learning (DRL) framework to estimate the optimal Dynamic Treatment Regimes from observational medical data. This framework is more flexible and adaptive for high dimensional action and state spaces than existing reinforcement learning methods to model real-life complexity in heterogeneous disease progression and treatment choices, with the goal of providing doctor and patients the data-driven personalized decision recommendations. The proposed DRL framework comprises (i) a supervised learning step to predict the most possible expert actions, and (ii) a deep reinforcement learning step to estimate the long-term value function of Dynamic Treatment Regimes. Both steps depend on deep neural networks. As a key motivational example, we have implemented the proposed framework on a data set from the Center for International Bone Marrow Transplant Research (CIBMTR) registry database, focusing on the sequence of prevention and treatments for acute and chronic graft versus host disease after transplantation. In the experimental results, we have demonstrated promising accuracy in predicting human experts' decisions, as well as the high expected reward function in the DRL-based dynamic treatment regimes.


title: On the breakdown of the simple Arrhenius law in the normal liquid state

abstract: It is common practice to discuss the temperature effect on molecular dynamics of glass formers above the melting temperature in terms of the Arrhenius law. Using dielectric spectroscopy measurements of dc-conductivity and structural relaxation time on the example of the typical glass former propylene carbonate, we provide experimental evidence that this practice is not justified. Our conclusions are supported by employing thermodynamic density scaling and the occurrence of inflection points in isothermal dynamic data measured at elevated pressure. Additionally, we propose a more suitable approach to describe the dynamics both above and below the inflection point based upon a modified MYEGA model.


title: Coherence-enhanced, phase-dependent dissipation in long SNS Josephson junctions: revealing Andreev Bound States dynamics

abstract: One of the best known causes of dissipation in ac driven quantum systems stems from photon absorption. Dissipation can also be caused by the retarded response to the time-dependent excitation, and in general gives insight into the system's relaxation times and mechanisms. We address the dissipation in a mesoscopic normal wire with superconducting contacts, that sustains a supercurrent at zero frequency and that may be expected to remain dissipationless at frequency lower than the superconducting gap. We probe the high frequency linear response of a Normal/Superconductor ring to a time-dependent flux by coupling it to a highly sensitive multimode microwave resonator. Far from being the simple derivative of the current-phase relation, the ring's ac susceptibility also displays a dissipative component whose phase dependence is a signature of the dynamical processes occurring within the Andreev spectrum. We show how dissipation is driven by the competition between the two aforementioned mechanisms. Depending on the relative strength of those contributions, dissipation can be maximal at $\pi$, when the minigap closes, or can be maximal near $ \pi/2$, when the dc supercurrent is maximal. We also find that the dissipative response increases at low temperature and can even exceed the normal state conductance. The results are confronted with predictions of the Kubo linear response and time-dependent Usadel equations. This experiment shows the power of the ac susceptibility measurement of individual hybrid mesoscopic systems in probing in a controlled way the quantum dynamics of ABS. By spanning different physical regimes, our experiments provide a unique access to inelastic scattering and spectroscopy of an isolated quantum coherent system. This technique should be a tool of choice to investigate topological superconductivity and detect the topological protection of edge states.


title: Hierarchies of Length Scale Based Typology in Anisotropic Multiband Superconductor

abstract: Since Ginzburg and Landau's seminal work in 1950 superconducting states have been classified by the hierarchy of the fundamental length scales of the theory; the magnetic field penetration lengths and coherence lengths. In the simplest single-component case they form a dimensionless ratio \kappa. As pointed out by Ginzburg in 1952, in general the ratio depends on the direction of the applied magnetic field due to material anisotropies. Single component materials can therefore display type-1 superconductivity when the field is applied in one direction and type-2 when the field is applied in a different direction. In this paper we expand the above length scale analysis to anisotropic multi-component superconductors, that can have multiple coherence lengths as well as multiple magnetic field penetration lengths, leading to unconvential length scale hierarchies. We demonstrate that the anisotropies in multi-band superconductors cannot in general be rescaled, leading to new regimes with various mixed hierarchies in different directions. For example, a regime is possible, where for a field applied in a certain direction coherence lengths are smaller than the magnetic field penetration lengths in one of the perpendicular directions, where as the penetration lengths are lager in the another direction. We demonstrate a new regime, where vortices form stripes in the direction where coherence lengths exceed the magnetic field penetration length and vortex cores overlap, while the vortex stripes repel each other in the orthogonal direction where the magnetic field penetration length exceeds the coherence lengths.


title: Algebraic dependencies and PSPACE algorithms in approximative complexity

abstract: Testing whether a set $\mathbf{f}$ of polynomials has an algebraic dependence is a basic problem with several applications. The polynomials are given as algebraic circuits. Algebraic independence testing question is wide open over finite fields (Dvir, Gabizon, Wigderson, FOCS'07). The best complexity known is NP$^{\#\rm P}$ (Mittmann, Saxena, Scheiblechner, Trans.AMS'14). In this work we put the problem in AM $\cap$ coAM. In particular, dependence testing is unlikely to be NP-hard and joins the league of problems of "intermediate" complexity, eg. graph isomorphism & integer factoring. Our proof method is algebro-geometric-- estimating the size of the image/preimage of the polynomial map $\mathbf{f}$ over the finite field. A gap in this size is utilized in the AM protocols. Next, we study the open question of testing whether every annihilator of $\mathbf{f}$ has zero constant term (Kayal, CCC'09). We give a geometric characterization using Zariski closure of the image of $\mathbf{f}$; introducing a new problem called approximate polynomials satisfiability (APS). We show that APS is NP-hard and, using projective algebraic-geometry ideas, we put APS in PSPACE (prior best was EXPSPACE via Grobner basis computation). As an unexpected application of this to approximative complexity theory we get-- Over any field, hitting-set for $\overline{\rm VP}$ can be designed in PSPACE. This solves an open problem posed in (Mulmuley, FOCS'12, J.AMS 2017); greatly mitigating the GCT Chasm (exponentially in terms of space complexity).


title: Uniqueness of the blow-up at isolated singularities for the Alt-Caffarelli functional

abstract: In this paper we prove uniqueness of blow-ups and $C^{1,\log}$-regularity for the free-boundary of minimizers of the Alt-Caffarelli functional at points where one blow-up has an isolated singularity. We do this by establishing a (log-)epiperimetric inequality for the Weiss energy for traces close to that of a cone with isolated singularity, whose free-boundary is graphical and smooth over that of the cone in the sphere. With additional assumptions on the cone, we can prove a classical epiperimetric inequality which can be applied to deduce a $C^{1,\alpha}$ regularity result. We also show that these additional assumptions are satisfied by the De Silva-Jerison-type cones, which are the only known examples of minimizing cones with isolated singularity. Our approach draws a connection between epiperimetric inequalities and the \L ojasiewicz inequality, and, to our knowledge, provides the first regularity result at singular points in the one-phase Bernoulli problem.


title: Experimental demonstration of shaken lattice interferometry

abstract: We experimentally demonstrate a shaken lattice interferometer. Atoms are trapped in the ground Bloch state of a red-detuned optical lattice. Using a closed-loop optimization protocol based on the dCRAB algorithm, we phase-modulate (shake) the lattice to transform the atom momentum state. In this way, we implement an atom beamsplitter and build five interferometers of varying interrogation times $T_I$. The sensitivity of shaken lattice interferometry is shown to scale as $T_I^2$, consistent with simulation [1]. Finally, we show that we can measure the sign of an applied signal and optimize the interferometer in the presence of a bias signal. [1] C. A. Weidner, H. Yu, R. Kosloff, and D. Z. Anderson, Phys. Rev. A 95, 043624 (2017).


title: Collisional Charging of Individual Sub-Millimeter Particles: Using Ultrasonic Levitation to Initiate and Track Charge Transfer

abstract: Electrostatic charging of insulating fine particles can be responsible for numerous phenomena ranging from lightning in volcanic plumes to dust explosions. However, even basic aspects of how fine particles become charged are still unclear. Studying particle charging is challenging because it usually involves the complexities associated with many particle collisions. To address these issues we introduce a method based on acoustic levitation, which makes it possible to initiate sequences of repeated collisions of a single sub-millimeter particle with a flat plate, and to precisely measure the particle charge in-situ after each collision. We show that collisional charge transfer between insulators is dependent on the hydrophobicity of the contacting surfaces. We use glass, which we modify by attaching nonpolar molecules to the particle, the plate, or both. We find that hydrophilic surfaces develop significant positive charges after contacting hydrophobic surfaces. Moreover, we demonstrate that charging between a hydrophilic and a hydrophobic surface is suppressed in an acidic environment and enhanced in a basic one. Application of an electric field during each collision is found to modify the charge transfer, again depending on surface hydrophobicity. We discuss these results within the context of contact charging due to ion transfer and show that they lend strong support to OH- ions as the charge carriers.


title: Topological Poincar\'e type inequalities and lower bounds on the infimum of the spectrum for graphs

abstract: We study topological Poincar\'e type inequalities on general graphs. We characterize graphs satisfying such inequalities and then turn to the best constants in these inequalities. Invoking suitable metrics we can interpret these constants geometrically as diameters and inradii. Moreover, we can relate them to spectral theory of Laplacians once a probability measure on the graph is chosen. More specifically, we obtain a variational characterization of these constants as infimum over spectral gaps of all Laplacians on the graphs associated to probability measures.


title: Constraining Primordial non-Gaussianity with Bispectrum and Power Spectum from Upcoming Optical and Radio Surveys

abstract: We forecast constraints on primordial non-Gaussianity (PNG) and bias parameters from measurements of galaxy power spectrum and bispectrum in future radio continuum and optical surveys. In the galaxy bispectrum, we consider a comprehensive list of effects, including the bias expansion for non-Gaussian initial conditions up to second order, redshift space distortions, redshift uncertainties and theoretical errors. These effects are all combined in a single PNG forecast for the first time. Moreover, we improve the bispectrum modelling over previous forecasts, by accounting for trispectrum contributions. All effects have an impact on final predicted bounds, which varies with the type of survey. We find that the bispectrum can lead to improvements up to a factor $\sim 5$ over bounds based on the power spectrum alone, leading to significantly better constraints for local-type PNG, with respect to current limits from \textit{Planck}. Future radio and photometric surveys could obtain a measurement error of $\sigma(f_{\mathrm{NL}}^{\mathrm{loc}}) \approx 0.2$. In the case of equilateral PNG, galaxy bispectrum can improve upon present bounds only if significant improvements in the redshift determinations of future, large volume, photometric or radio surveys could be achieved. For orthogonal non-Gaussianity, expected constraints are generally comparable to current ones.


title: Optimizing Reweighted Belief Propagation for Distributed Likelihood Fusion Problems

abstract: Belief propagation (BP) is a powerful tool to solve distributed inference problems, though it is limited by short cycles in the corresponding factor graph. Such cycles may lead to incorrect solutions or oscillatory behavior. Only for certain types of problems are convergence properties understood. We extend this knowledge by investigating the use of reweighted BP for distributed likelihood fusion problems, which are characterized by equality constraints along possibly short cycles. Through a linear formulation of BP, we are able to analytically derive convergence conditions for certain types of graphs and optimize the convergence speed. We compare with standard belief consensus and observe significantly faster convergence.


title: Spectral Shape Preserving Approximation

abstract: We introduce an algorithm of joint approximation of a function and its first derivative by alternative orthogonal polynomials on the interval [0,1].The algorithm exhibits properties of shape preserving approximation for the function. A weak formulation of approximation is presented. An example on shape preserving extrapolation is given. The weak form is reduced for approximation on a discrete set of abscissas. Also, we introduce a new system of orthogonal functions with nice properties - structured orthogonal polynomials - and show that the system can be employed for a different kind of joint approximation of a function and its first derivative and may have property of shape preserving approximation. In addition, we show that structured orthogonal polynomials generate wavelet functions We complement these results with definition of structured semi-orthogonal polynomials and introduce wavelet basis functions.


title: Growth of mod$-2$ homology in higher rank locally symmetric spaces

abstract: Let $X$ be a higher rank symmetric space or a Bruhat-Tits building of dimension at least $2$ such that the isometry group of $X$ has property $(T)$. We prove that for every torsion free lattice $\Gamma\subset {\rm Isom} X$ any homology class in $H_1(\Gamma\backslash X,\mathbb F_2)$ has a representative cycle of total length $o_X({\rm Vol}(\Gamma\backslash X))$. As an application we show that $\dim_{\mathbb F_2} H_1(\Gamma\backslash X,\mathbb F_2)=o_X({\rm Vol}(\Gamma\backslash X)).$


title: Statistical analysis of the effect of the current, potential and proposed rules of a game in tennis

abstract: With the aid of mathematical modelling (basic tool is the random walk with absorbing barriers) we derive subsequent formulas to study the effect of different versions of possible rules. For different rules the probability of winning a game, the probability of break point occurrence, the mathematical expectation of the number of rallies (points) and, the mathematical expectation of the number of break points in a game are expressed. We check these rules against ATP statistics for the Top-200 men players. In conclusion, we suggest a slight but essential modification for the rule of a tennis game, namely , second service ( in case of a first service fault) is to be allowed only at the first three points (rallies). This would partially preserve the traditions (server has an advantage in the modern game) and at the same time it would reduce the predictability of the game, significantly increasing in this way the excitement for the spectators.


title: On a theorem of Bertolini-Darmon about rationality of Stark-Heegner points over genus fields of real quadratic fields

abstract: In this paper, we remove certain hypothesis in the theorem of Bertolini-Darmon on the rationality of Stark-Heegner points over narrow genus class fields of real quadratic fields. Along the way, we establish that certain normalized special values of $L$-functions are squares of rational numbers, a result that is of independent interest, and can be regarded as instances of the rank zero case of the Birch and Swinnerton-dyer conjecture modulo squares.


title: Low-Latency Heterogeneous Networks with Millimeter-Wave Communications

abstract: Heterogeneous network (HetNet) is a key enabler to largely boost network coverage and capacity in the forthcoming fifth-generation (5G) and beyond. To support the explosively growing mobile data volumes, wireless communications with millimeter-wave (mm-wave) radios have attracted massive attention, which is widely considered as a promising candidate in 5G HetNets. In this article, we give an overview on the end-to-end latency of HetNets with mm-wave communications. In general, it is rather challenging for formulating and optimizing the delay problem with buffers in mm-wave communications, since conventional graph-based network optimization techniques are not applicable when queues are considered. Toward this end, we develop an adaptive low-latency strategy, which uses cooperative networking to reduce the end-to-end latency. Then, we evaluate the performance of the introduced strategy. Results reveal the importance of proper cooperative networking in reducing the end-to-end latency. In addition, we have identified several challenges in future research for low-latency mm-wave HetNets.


title: Structural distortion behind the nematic superconductivity in Sr$_x$Bi$_2$Se$_3$

abstract: An archetypical layered topological insulator Bi$_2$Se$_3$ becomes superconductive upon doping with Sr, Nb or Cu. Superconducting properties of these materials in the presence of in-plane magnetic field demonstrate spontaneous symmetry breaking: 180$^\circ$-rotation symmetry of superconductivity versus 120$^\circ$-rotation symmetry of the crystal. Such behavior brilliantly confirms nematic topological superconductivity. To what extent this nematicity is due to superconducting pairing in these materials, rather than due to crystal structure distortions? This question remained unanswered, because so far no visible deviations from the 3-fold crystal symmetry were resolved in these materials. To address this question we grow high quality single crystals of Sr$_x$Bi$_2$Se$_3$, perform detailed X-ray diffraction and magnetotransport studies and reveal that the observed superconducting nematicity direction correlates with the direction of small structural distortions in these samples( $\sim 0.02$\% elongation in one crystallographic direction). Additional anisotropy comes from orientation of the crystallite axes. 2-fold symmetry of magnetoresistance observed in the most uniform crystals well above critical temperature demonstrates that these structural distortions are nevertheless strong enough. Our data in combination with strong sample-to-sample variation of the superconductive anisotropy parameter are indicative for significance of the structural factor in the apparent nematic superconductivity in Sr$_x$Bi$_2$Se$_3$.


title: Disinformation Warfare: Understanding State-Sponsored Trolls on Twitter and Their Influence on the Web

abstract: Over the past couple of years, anecdotal evidence has emerged linking coordinated campaigns by state-sponsored actors with efforts to manipulate public opinion on the Web, often around major political events, through dedicated accounts, or "trolls." Although they are often involved in spreading disinformation on social media, there is little understanding of how these trolls operate, what type of content they disseminate, and most importantly their influence on the information ecosystem. In this paper, we shed light on these questions by analyzing 27K tweets posted by 1K Twitter users identified as having ties with Russia's Internet Research Agency and thus likely state-sponsored trolls. We compare their behavior to a random set of Twitter users, finding interesting differences in terms of the content they disseminate, the evolution of their account, as well as their general behavior and use of Twitter. Then, using Hawkes Processes, we quantify the influence that trolls had on the dissemination of news on social platforms like Twitter, Reddit, and 4chan. Overall, our findings indicate that Russian trolls managed to stay active for long periods of time and to reach a substantial number of Twitter users with their tweets. When looking at their ability of spreading news content and making it viral, however, we find that their effect on social platforms was minor, with the significant exception of news published by the Russian state-sponsored news outlet RT (Russia Today).


title: Data-Driven Approximate Abstraction for Black-Box Piecewise Affine Systems

abstract: How to effectively and reliably guarantee the correct functioning of safety-critical cyber-physical systems in uncertain conditions is a challenging problem. This paper presents a data-driven algorithm to derive approximate abstractions for piecewise affine systems with unknown dynamics. It advocates a significant shift from the current paradigm of abstraction, which starts from a model with known dynamics. Given a black-box system with unknown dynamics and a linear temporal logic specification, the proposed algorithm is able to obtain an abstraction of the system with an arbitrarily small error and a bounded probability. The algorithm consists of three components, system identification, system abstraction, and active sampling. The effectiveness of the algorithm is demonstrated by a case study with a soft robot.


title: Surfactant and gravity dependent instability of two-layer channel flows: Linear theory covering all wave lengths

abstract: A linear stability analysis of a two-layer plane Couette flow of two immiscible fluid layers with different densities, viscosities and thicknesses, bounded by two infinite parallel plates moving at a constant relative velocity to each other, with an insoluble surfactant along the interface and in the presence of gravity is carried out. The normal modes approach is applied to the equations governing flow disturbances. These equations, together with boundary conditions at the plates and the interface, yield a linear eigenvalue problem. When inertia is neglected velocity amplitudes are linear combinations of hyperbolic functions, and a quadratic dispersion equation for the complex growth rate is obtained where coefficients depend on the aspect ratio, the viscosity ratio, the basic velocity shear, the Marangoni number Ma that measures the effects of surfactant, and the Bond number Bo that measures the influence of gravity. An extensive investigation is carried out that examines the stabilizing or destabilizing influences of these parameters. There are two continuous branches of the normal modes: a robust branch that exists even with no surfactant, and a surfactant branch that vanishes when Ma $\downarrow 0$. Due to the availability of the explicit forms for the growth rates, in many instances the numerical results are corroborated with analytical asymptotics. For the less unstable branch, a mid-wave interval of unstable wavenumbers (Halpern and Frenkel (2003)) sometimes co-exists with a long-wave one. We study the instability landscape, determined by the threshold curve of the long-wave instability and the critical curve of the mid-wave instability in the (Ma, Bo)-plane. The changes of the extremal points of the critical curves with the variation of the other parameters, such as the viscosity ratio, and the extrema bifurcation points are investigated.


title: ALMA spectral line and imaging survey of a low and a high mass-loss rate AGB star between 335 and 362 GHz

abstract: A spectral line and imaging survey of the low mass-loss rate AGB star R Dor (Mdot ~ 1e-7 Msun/yr) and the high mass-loss rate AGB star IK Tau (Mdot ~5e-6 Msun/yr) was made with ALMA between 335 and 362 GHz at a spatial resolution of ~150 mas, corresponding to the locus of the main dust formation region of both targets. Some 200 spectral features from 15 molecules (and their isotopologues) were observed, including rotational lines in both the ground and vibrationally excited states. Detected species include the gaseous precursors of dust grains such as SiO, AlO, AlOH, TiO, and TiO2. We present a spectral atlas for both stars and the parameters of all detected spectral features. A clear dichotomy for the sulphur chemistry is seen: while CS, SiS, SO, and SO2 are abundantly present in IK Tau, only SO and SO2 are detected in R Dor. Also other species such as NaCl, NS, AlO, and AlOH display a completely different behaviour. From some selected species, the minor isotopologues can be used to assess the isotopic ratios. The channel maps of many species prove that both large and small-scale inhomogeneities persist in the inner wind of both stars in the form of blobs, arcs, and/or a disk. The high sensitivity of ALMA allows us to spot the impact of these correlated density structures in the spectral line profiles. The spectral lines often display a half width at zero intensity much larger than expected from the terminal velocity, v_inf, previously derived for both objects (36 km/s versus v_inf ~17.7 km/s for IK Tau and 23 km/s versus v_inf ~5.5 km/s for R Dor). Both a more complex 3D morphology and a more forceful wind acceleration of the (underlying) isotropic wind can explain this trend. The formation of fractal grains in the region beyond ~400 mas can potentially account for the latter scenario. From the continuum map, we deduce a dust mass of ~3.7e-7 Msun for IK Tau and ~2e-8 Msun for R Dor.


title: Multiple Object Detection, Tracking and Long-Term Dynamics Learning in Large 3D Maps

abstract: In this work, we present a method for tracking and learning the dynamics of all objects in a large scale robot environment. A mobile robot patrols the environment and visits the different locations one by one. Movable objects are discovered by change detection, and tracked throughout the robot deployment. For tracking, we extend the Rao-Blackwellized particle filter of previous work with birth and death processes, enabling the method to handle an arbitrary number of objects. Target births and associations are sampled using Gibbs sampling. The parameters of the system are then learnt using the Expectation Maximization algorithm in an unsupervised fashion. The system therefore enables learning of the dynamics of one particular environment, and of its objects. The algorithm is evaluated on data collected autonomously by a mobile robot in an office environment during a real-world deployment. We show that the algorithm automatically identifies and tracks the moving objects within 3D maps and infers plausible dynamics models, significantly decreasing the modeling bias of our previous work. The proposed method represents an improvement over previous methods for environment dynamics learning as it allows for learning of fine grained processes.


title: Application of Kriging Models for a Drug Combination Experiment on Lung Cancer

abstract: Combinatorial drugs have been widely applied in disease treatment, especially chemotherapy for cancer, due to its improved efficacy and reduced toxicity compared with individual drugs. The study of combinatorial drugs requires efficient experimental designs and proper follow-up statistical modelling techniques. Linear and non-linear models are often used in the response surface modelling for such experiments. We propose the use of Kriging models to better depict the response surfaces of combinatorial drugs and take into account the measurement error. We further study how proper experimental designs can reduce the required number of runs. We illustrate our method via a combinatorial drug experiment on lung cancer. We demonstrate that only 27 runs are needed to predict all 512 runs in the original experiment and achieve better precision than existing analysis.


title: Angular momentum and parity of a two gluon system

abstract: In this work we obtain all the allowed values of the angular momentum and parity for two massless non-Abelian fields. In particular, we show that, unlike the two photons system, two gluons can be in a state with $J^P=1^\pm$, $(2k)^-$, and $(2k+1)^-$, $k\geq1$. This is a general result since it was obtained using only the laws of addition of the angular momentum.


title: The recovery of George Berkeley's objective science of 1710 and its implications for traditional science

abstract: In 1710, George Berkeley, philospher and scientist, discovered that the directions in which we see (sightlines) are tilted 45 degrees towards the past from the directions in which we look. Looking north from the 3-d present we see to the north-past and so see the 4-d past. The past only contains images of the material objects in the present, so the world we see contains no material objects. We do not see a tree; we see an image of the tree, or, with two mirrors, an infinite stream of images of the tree receding into the 4-d past. The science of the 4-d past described with data-defined language is objective. When the Relativistic Doppler Effect is so described, it is simple and complete, unlike Einstein's space and time version in his 1905 theory of relativity. Scientists still use space and time science though so Berkeley's discoveries have been neglected and lost. Their recovery in this paper reinststes the choice between a data-defined language description of the observed 4-d past and an idea-based language description of idea-based (3+1)-d space and time.


title: The Gaussian Double-Bubble Conjecture

abstract: We establish the Gaussian Double-Bubble Conjecture: the least Gaussian-weighted perimeter way to decompose $\mathbb{R}^n$ into three cells of prescribed (positive) Gaussian measure is to use a tripod-cluster, whose interfaces consist of three half-hyperplanes meeting along an $(n-2)$-dimensional plane at $120^{\circ}$ angles (forming a tripod or "Y" shape in the plane). Moreover, we prove that tripod-clusters are the unique isoperimetric minimizers (up to null-sets).


title: Time Coordinates and Clocks: Einstein's Struggle

abstract: In his Autobiographical Notes, Einstein mentioned that on his road to the final theory of general relativity it was a major difficulty to accustom himself to the idea that coordinates need not possess an immediate physical meaning in terms of lengths and times. This appears strange: that coordinates are conventional markers of events seems an obvious fact, already familiar from pre-relativistic physics. In this paper we explore the background of Einsteins difficulties, going from his 1905 paper on special relativity, through his 1907 and 1911 papers on the consequences of the equivalence principle, to the 1916 review paper on the general theory. As we shall argue, Einstein's problems were intimately connected to his early methodology, in which clarity achieved by concrete physical pictures played an essential role; and to the related fact that on his route to the general theory he focused on special situations that were easily accessible to physical intuition. The details of this background of Einstein's early reasoning have not always been sufficiently appreciated in modern commentaries. As we shall see, this has led to erroneous judgments about the status and validity of some of the early relativistic derivations by Einstein and others, in particular concerning the gravitational redshift.


title: Lower bound for the perimeter density at singular points of a minimizing cluster in $\mathbb R^N$

abstract: In this paper we study the blow-ups of the singular points in the boundary of a minimizing cluster lying in the interface of more than two chambers. We establish a sharp lower bound for the perimeter density at those points and we prove that this bound is rigid, namely having the lowest possible density completely characterizes the blow-up.


title: Adapting The Gibbs Sampler

abstract: The popularity of Adaptive MCMC has been fueled on the one hand by its success in applications, and on the other hand, by mathematically appealing and computationally straightforward optimisation criteria for the Metropolis algorithm acceptance rate (and, equivalently, proposal scale). Similarly principled and operational criteria for optimising the selection probabilities of the Random Scan Gibbs Sampler have not been devised to date. In the present work, we close this gap and develop a general purpose Adaptive Random Scan Gibbs Sampler that adapts the selection probabilities. The adaptation is guided by optimising the $L_2-$spectral gap for the target's Gaussian analogue, gradually, as target's global covariance is learned by the sampler. The additional computational cost of the adaptation represents a small fraction of the total simulation effort. ` We present a number of moderately- and high-dimensional examples, including truncated Gaussians, Bayesian Hierarchical Models and Hidden Markov Models, where significant computational gains are empirically observed for both, Adaptive Gibbs, and Adaptive Metropolis within Adaptive Gibbs version of the algorithm. We argue that Adaptive Random Scan Gibbs Samplers can be routinely implemented and substantial computational gains will be observed across many typical Gibbs sampling problems. We shall give conditions under which ergodicity of the adaptive algorithms can be established.


title: Adaptation of olfactory receptor abundances for efficient coding

abstract: Olfactory receptor usage is highly heterogeneous, with some receptor types being orders of magnitude more abundant than others. We propose an explanation for this striking fact: the receptor distribution is tuned to maximally represent information about the olfactory environment in a regime of efficient coding that is sensitive to the global context of correlated sensor responses. This model predicts that in mammals, where olfactory sensory neurons are replaced regularly, receptor abundances should continuously adapt to odor statistics. Experimentally, increased exposure to odorants leads variously, but reproducibly, to increased, decreased, or unchanged abundances of different activated receptors. We demonstrate that this diversity of effects is required for efficient coding when sensors are broadly correlated, and provide an algorithm for predicting which olfactory receptors should increase or decrease in abundance following specific environmental changes. Finally, we give simple dynamical rules for neural birth and death processes that might underlie this adaptation.


title: Model-theoretic Elekes-Szab\'o in the strongly minimal case

abstract: We prove a generalizations of the Elekes-Szab\'o theorem for relations definable in strongly minimal structures that are interpretable in distal structures.


title: Constraints on metal oxide and metal hydroxide abundances in the winds of AGB stars - Potential detection of FeO in R Dor

abstract: Using ALMA, we observed the stellar wind of two oxygen-rich Asymptotic Giant Branch (AGB) stars, IK Tau and R Dor, between 335 and 362 GHz. One aim was to detect metal oxides and metal hydroxides (AlO, AlOH, FeO, MgO, MgOH), some of which are thought to be direct precursors of dust nucleation and growth. We report on the potential first detection of FeO (v=0, Omega=4, J=11-10) in RDor (mass-loss rate, Mdot, ~1e-7 Msun/yr). The presence of FeO in IK Tau (Mdot~5e-6 Msun/yr) cannot be confirmed due to a blend with 29SiS, a molecule that is absent in R Dor. The detection of AlO in R Dor and of AlOH in IK Tau was reported earlier by Decin et al. (2017). All other metal oxides and hydroxides, as well as MgS, remain undetected. We derive a column density N(FeO) of 1.1+/-0.9e15 cm^{-2} in R Dor, or a fractional abundance [FeO/H]~1.5e-8 accounting for non-LTE effects. The derived fractional abundance [FeO/H] is a factor ~20 larger than conventional gas-phase chemical kinetic predictions. This discrepancy may be partly accounted for by the role of vibrationally excited OH in oxidizing Fe, or may be evidence for other currently unrecognised chemical pathways producing FeO. Assuming a constant fractional abundance w.r.t. H_2, the upper limits for the other metals are [MgO/H_2] <5.5e-10 (R Dor) and <7e-11 (IK Tau), [MgOH/H_2] <9e-9 (R Dor) and <1e-9 (IK Tau), [CaO/H_2] <2.5e-9 (R Dor) and <1e-10 (IK Tau), [CaOH/H_2] <6.5e-9 (R Dor) and <9e-10 (IK Tau), and [MgS/H_$] <4.5e-10 (R Dor) and <6e-11 (IK Tau). The retrieved upper limit abundances for these latter molecules are in accord with the chemical model predictions.


title: HONE: Higher-Order Network Embeddings

abstract: This paper describes a general framework for learning Higher-Order Network Embeddings (HONE) from graph data based on network motifs. The HONE framework is highly expressive and flexible with many interchangeable components. The experimental results demonstrate the effectiveness of learning higher-order network representations. In all cases, HONE outperforms recent embedding methods that are unable to capture higher-order structures with a mean relative gain in AUC of $19\%$ (and up to $75\%$ gain) across a wide variety of networks and embedding methods.


title: Spiral Galaxy Rotation Curves Without Dark Matter or MOND -- Two Conjectures

abstract: Usual explanations of spiral galaxy rotation curves assume circular orbits of stars. The consequences of giving up this assumption are investigated here. In particular, hyperbolic stellar trajectories are found to be interesting. The two suggested models for the production of such trajectories will also explain the observed flat rotation curves without the postulation of dark matter or MOND. It is suggested that spiral galaxies may have started as compact objects with significant angular momenta and then disintegrated. The first model conjectures the existence of a spinning hot disk around a spherical galactic core. The disk is held together by local gravity and electromagnetic scattering forces. However, it disintegrates at the edge producing fragments that form stars. Once separated from the disk, the stars experience only the centrally directed gravitational force due to the massive core and remaining disk. A numerical simulation shows that a high enough angular velocity of the disk produces hyperbolic stellar trajectories that agree with the observed rotation curves. The second model conjectures a significant initial thermonuclear event that produces a dust plume along with large stars. This dust plume is made of ordinary matter. However, it acts like the postulated dark matter in producing initial circular trajectories. Unlike dark matter, the plume can be shown to escape the galaxy rapidly causing the star trajectories to evolve to hyperbolic shapes. This process can be seen to produce the observed rotation curves due to the initial circular orbits. Also, as the plume dissipates rapidly it does not obfuscate the stars from view. Both models have weaknesses as do the currently known models using dark matter or MOND.


title: Von Neumann's Impossibility Proof: Mathematics in the Service of Rhetorics

abstract: According to what has become a standard history of quantum mechanics, in 1932 von Neumann persuaded the physics community that hidden variables are impossible as a matter of principle, after which leading proponents of the Copenhagen interpretation put the situation to good use by arguing that the completeness of quantum mechanics was undeniable. This state of affairs lasted, so the story continues, until Bell in 1966 exposed von Neumann's proof as obviously wrong. The realization that von Neumann's proof was fallacious then rehabilitated hidden variables and made serious foundational research possible again. It is often added in recent accounts that von Neumann's error had been spotted almost immediately by Grete Hermann, but that her discovery was of no effect due to the dominant Copenhagen Zeitgeist. We shall attempt to tell a story that is more historically accurate and less ideologically charged. Most importantly, von Neumann never claimed to have shown the impossibility of hidden variables tout court, but argued that hidden-variable theories must possess a structure that deviates fundamentally from that of quantum mechanics. Both Hermann and Bell appear to have missed this point, moreover, both raised unjustified technical objections to the proof. Von Neumann's argument was basically that hidden-variables schemes must violate the "quantum principle" that physical quantities are to be represented by operators in a Hilbert space. As a consequence, hidden-variables schemes, though possible in principle, necessarily exhibit a certain kind of contextuality.


title: Optimal Beam Sweeping and Communication in Mobile Millimeter-Wave Networks

abstract: Millimeter-wave (mm-wave) communications incur a high beam alignment cost in mobile scenarios such as vehicular networks. Therefore, an efficient beam alignment mechanism is required to mitigate the resulting overhead. In this paper, a one-dimensional mobility model is proposed where a mobile user (MU), such as a vehicle, moves along a straight road with time-varying and random speed, and communicates with base stations (BSs) located on the roadside over the mm-wave band. To compensate for location uncertainty, the BS widens its transmission beam and, when a critical beamwidth is achieved, it performs beam-sweeping to refine the MU position estimate, followed by data communication over a narrow beam. The average rate and average transmission power are computed in closed form and the optimal beamwidth for communication, number of sweeping beams, and transmission power allocation are derived so as to maximize the average rate under an average power constraint. Structural properties of the optimal design are proved, and a bisection algorithm to determine the optimal sweeping -- communication parameters is designed. It is shown numerically that an adaptation of the IEEE 802.11ad standard to the proposed model exhibits up to 90% degradation in spectral efficiency compared to the proposed scheme.


title: Quantum Mechanics and Perspectivalism

abstract: Experimental evidene of the last decades has made the status of "collapses of the wave function" even more shaky than it already was on conceptual grounds: interference effects turn out to be detectable even when collapses are typically expected to occur. Non-collapse interpretations should consequently be taken seriously. In this paper we argue that such interpretations suggest a perspectivalism according to which quantum objects are not characterized by monadic properties, but by relations to other systems. Accordingly, physical systems may possess different properties with respect to different "reference systems". We discuss some of the relevant arguments, and argue that perspectivalism both evades recent arguments that single-world interpretations are inconsistent and eliminates the need for a privileged rest frame in the relativistic case.


title: A Survey on Behavioral Biometric Authentication on Smartphones

abstract: Recent research has shown the possibility of using smartphones' sensors and accessories to extract some behavioral attributes such as touch dynamics, keystroke dynamics and gait recognition. These attributes are known as behavioral biometrics and could be used to verify or identify users implicitly and continuously on smartphones. The authentication systems that have been built based on these behavioral biometric traits are known as active or continuous authentication systems. This paper provides a review of the active authentication systems. We present the components and the operating process of the active authentication systems in general, followed by an overview of the state-of-the-art behavioral biometric traits that used to develop an active authentication systems and their evaluation on smartphones. We discuss the issues, strengths and limitations that associated with each behavioral biometric trait. Also, we introduce a comparative summary between them. Finally, challenges and open research problems are presented in this research field.


title: Air Markov Chain Monte Carlo

abstract: We introduce a class of Adapted Increasingly Rarely Markov Chain Monte Carlo (AirMCMC) algorithms where the underlying Markov kernel is allowed to be changed based on the whole available chain output but only at specific time points separated by an increasing number of iterations. The main motivation is the ease of analysis of such algorithms. Under the assumption of either simultaneous or (weaker) local simultaneous geometric drift condition, or simultaneous polynomial drift we prove the $L_2-$convergence, Weak and Strong Laws of Large Numbers (WLLN, SLLN), Central Limit Theorem (CLT), and discuss how our approach extends the existing results. We argue that many of the known Adaptive MCMC algorithms may be transformed into the corresponding Air versions, and provide an empirical evidence that performance of the Air version stays virtually the same.


title: Metastable decoherence-free subspace and pointer states in mesoscopic quantum systems

abstract: Two initially correlated coherent states, each interacting with its own independent dissipative environment exhibit a sudden transition from classical to quantum decoherence. This change in the dynamics is a turning point in the decoherence, in the sense that depending on the average number of photons of each cavity, decoherence can even be suppressed. Indeed, the quantum state is time-independent for a time span in the mesoscopic regime, revealing a decoherence-free subspace. Furthermore, the absence of decoherence is manifested in the apparition of a metastable pointer state basis.


title: On the combinatorics of faces of trees and anodyne extensions of dendroidal sets

abstract: We discuss the combinatorics of faces of trees in the context of dendroidal sets and develop a systematic treatment of dendroidal anodyne extensions. As the main example and our motivation, we prove the pushout-product property for the stable model structure on dendroidal sets.


title: Tunable photonic heat transport in a quantum heat valve

abstract: Quantum thermodynamics is emerging both as a topic of fundamental research and as means to understand and potentially improve the performance of quantum devices. A prominent platform for achieving the necessary manipulation of quantum states is superconducting circuit quantum electrodynamics (QED). In this platform, thermalization of a quantum system can be achieved by interfacing the circuit QED subsystem with a thermal reservoir of appropriate Hilbert dimensionality. Here we study heat transport through an assembly consisting of a superconducting qubit capacitively coupled between two nominally identical coplanar waveguide resonators, each equipped with a heat reservoir in the form of a normal-metal mesoscopic resistor termination. We report the observation of tunable photonic heat transport through the resonator-qubit-resonator assembly, showing that the reservoir-to-reservoir heat flux depends on the interplay between the qubit-resonator and the resonator-reservoir couplings, yielding qualitatively dissimilar results in different coupling regimes. Our quantum heat valve is relevant for the realisation of quantum heat engines and refrigerators, that can be obtained, for example, by exploiting the time-domain dynamics and coherence of driven superconducting qubits. This effort would ultimately bridge the gap between the fields of quantum information and thermodynamics of mesoscopic systems.


title: Mixing and CPV in charm hadrons at LHCb

abstract: LHCb continues to expand its world-leading sample of charmed hadrons collected during LHC's Run 1 (2010-2012) and Run 2 (2015-present). This sample is yielding some of the most stringent tests of the Standard Model understanding of charm physics. This includes precise measurements of the neutral D-meson mixing parameters and some of the most sensitive searches for direct and indirect CP violation in charm interactions.


title: Asymmetric dark matter, baryon asymmetry and lepton number violation

abstract: We study the effect of lepton number violation (LNV) on baryon asymmetry, generated in the early Universe, in the presence of a dark sector with a global symmetry $U(1)_X$, featuring asymmetric dark matter (ADM). We show that in general LNV, observable at the LHC or in neutrinoless double beta decay experiments, cannot wash out a baryon asymmetry generated at higher scales, unlike in scenarios without such dark sector. An observation of LNV at the TeV scale may thus support ADM scenarios. Considering several models with different types of dark matter (DM), we find that the DM mass is of the order of a few GeV or below in our scenario.


title: A representative agent model based on risk-neutral prices

abstract: In this paper, we determine a representative agent model based on risk-neutral information. The main idea is that the pricing kernel is transition independent, which is supported by the well-known capital asset pricing theory. Determining the representative agent model is closely related to the eigenpair problem of a second-order differential operator. The purpose of this paper is to find all such eigenpairs which are financially or economically meaningful. We provide a necessary and sufficient condition for the existence of such pairs, and prove that that all the possible eignepairs can be expressed as a one-parameter family. Finally, we find a representative agent model derived from the eigenpairs.


title: Gelfand-Tsetlin Theory for Rational Galois Algebras

abstract: In the present paper we study Gelfand-Tsetlin modules defined in terms of BGG differential operators. The structure of these modules is described with the aid of the Postnikov-Stanley polynomials introduced in [PS09]. These polynomials are used to identify the action of the Gelfand-Tsetlin subalgebra on the BGG operators. We also provide explicit bases of the corresponding Gelfand-Tsetlin modules and prove a simplicity criterion for these modules. The results hold for modules defined over standard Galois orders of type $A$ - a large class of rings that include the universal enveloping algebra of $\mathfrak{gl} (n)$ and the finite $W$-algebras of type $A$.


title: A Cyber Science Based Ontology for Artificial General Intelligence Containment

abstract: The development of artificial general intelligence is considered by many to be inevitable. What such intelligence does after becoming aware is not so certain. To that end, research suggests that the likelihood of artificial general intelligence becoming hostile to humans is significant enough to warrant inquiry into methods to limit such potential. Thus, containment of artificial general intelligence is a timely and meaningful research topic. While there is limited research exploring possible containment strategies, such work is bounded by the underlying field the strategies draw upon. Accordingly, we set out to construct an ontology to describe necessary elements in any future containment technology. Using existing academic literature, we developed a single domain ontology containing five levels, 32 codes, and 32 associated descriptors. Further, we constructed ontology diagrams to demonstrate intended relationships. We then identified humans, AGI, and the cyber world as novel agent objects necessary for future containment activities. Collectively, the work addresses three critical gaps: (a) identifying and arranging fundamental constructs; (b) situating AGI containment within cyber science; and (c) developing scientific rigor within the field.


title: Towards monitoring conformational changes of the GPCR neurotensin receptor 1 by single-molecule FRET

abstract: Neurotensin receptor 1 (NTSR1) is a G protein-coupled receptor that is important for signaling in the brain and the gut. Its agonist ligand neurotensin (NTS), a 13-amino-acid peptide, binds with nanomolar affinity from the extracellular side to NTSR1 and induces conformational changes that trigger intracellular signaling processes. Our goal is to monitor the conformational dynamics of single fluorescently labeled NTSR1. For this, we fused the fluorescent protein mNeonGreen to the C terminus of NTSR1, purified the receptor fusion protein from E. coli membranes, and reconstituted NTSR1 into liposomes with E. coli polar lipids. Using single-molecule anisotropy measurements, NTSR1 was found to be monomeric in liposomes, with a small fraction being dimeric and oligomeric, showing homoFRET. Similar results were obtained for NTSR1 in detergent solution. Furthermore, we demonstrated agonist binding to NTSR1 by time-resolved single-molecule F\"orster resonance energy transfer (smFRET), using neurotensin labeled with the fluorophore ATTO594.


title: Less is more: sampling chemical space with active learning

abstract: The development of accurate and transferable machine learning (ML) potentials for predicting molecular energetics is a challenging task. The process of data generation to train such ML potentials is a task neither well understood nor researched in detail. In this work, we present a fully automated approach for the generation of datasets with the intent of training universal ML potentials. It is based on the concept of active learning (AL) via Query by Committee (QBC), which uses the disagreement between an ensemble of ML potentials to infer the reliability of the ensemble's prediction. QBC allows the presented AL algorithm to automatically sample regions of chemical space where the ML potential fails to accurately predict the potential energy. AL improves the overall fitness of ANAKIN-ME (ANI) deep learning potentials in rigorous test cases by mitigating human biases in deciding what new training data to use. AL also reduces the training set size to a fraction of the data required when using naive random sampling techniques. To provide validation of our AL approach we develop the COMP6 benchmark (publicly available on GitHub), which contains a diverse set of organic molecules. Through the AL process, it is shown that the AL-based potentials perform as well as the ANI-1 potential on COMP6 with only 10% of the data, and vastly outperforms ANI-1 with 25% the amount of data. Finally, we show that our proposed AL technique develops a universal ANI potential (ANI-1x) that provides accurate energy and force predictions on the entire COMP6 benchmark. This universal ML potential achieves a level of accuracy on par with the best ML potentials for single molecule or materials, while remaining applicable to the general class of organic molecules comprised of the elements CHNO.

categories: physics.comp-ph cs.LG physics.chem-ph stat.ML


title: Square Sierpi\'nski carpets and Latt\`es maps

abstract: We prove that every quasisymmetric homeomorphism of a standard square Sierpi\'nski carpet $S_p$, $p\ge 3$ odd, is an isometry. This strengthens and completes earlier work by the authors. We also show that a similar conclusion holds for quasisymmetries of the double of $S_p$ across the outer peripheral circle. Finally, as an application of the techniques developed in this paper, we prove that no standard square carpet $S_p$ is quasisymmetrically equivalent to the Julia set of a postcritically-finite rational map.


title: Document Image Classification with Intra-Domain Transfer Learning and Stacked Generalization of Deep Convolutional Neural Networks

abstract: In this work, a region-based Deep Convolutional Neural Network framework is proposed for document structure learning. The contribution of this work involves efficient training of region based classifiers and effective ensembling for document image classification. A primary level of `inter-domain' transfer learning is used by exporting weights from a pre-trained VGG16 architecture on the ImageNet dataset to train a document classifier on whole document images. Exploiting the nature of region based influence modelling, a secondary level of `intra-domain' transfer learning is used for rapid training of deep learning models for image segments. Finally, stacked generalization based ensembling is utilized for combining the predictions of the base deep neural network models. The proposed method achieves state-of-the-art accuracy of 92.2% on the popular RVL-CDIP document image dataset, exceeding benchmarks set by existing algorithms.


title: Benchmarking Clinical Decision Support Search

abstract: Finding relevant literature underpins the practice of evidence-based medicine. From 2014 to 2016, TREC conducted a clinical decision support track, wherein participants were tasked with finding articles relevant to clinical questions posed by physicians. In total, 87 teams have participated over the past three years, generating 395 runs. During this period, each team has trialled a variety of methods. While there was significant overlap in the methods employed by different teams, the results were varied. Due to the diversity of the platforms used, the results arising from the different techniques are not directly comparable, reducing the ability to build on previous work. By using a stable platform, we have been able to compare different document and query processing techniques, allowing us to experiment with different search parameters. We have used our system to reproduce leading teams runs, and compare the results obtained. By benchmarking our indexing and search techniques, we can statistically test a variety of hypotheses, paving the way for further research.


title: Momentum exchange between an electromagnetic wave and a dispersive medium

abstract: We present an elementary discussion of the momentum transferred by an electromagnetic wave propagating in a dispersive medium. Our analysis is based on Minkowski's electromagnetic momentum density which have been recently seen to be consistent with a fully covariant expression of the energy-momentum tensor of the electromagnetic field in a dispersive medium and with all the experimental evidence. We show that the medium may be either pulled or pushed as the electromagnetic pulse enters in it depending on the value of the frequency.


title: Strong error analysis for stochastic gradient descent optimization algorithms

abstract: Stochastic gradient descent (SGD) optimization algorithms are key ingredients in a series of machine learning applications. In this article we perform a rigorous strong error analysis for SGD optimization algorithms. In particular, we prove for every arbitrarily small $\varepsilon \in (0,\infty)$ and every arbitrarily large $p\in (0,\infty)$ that the considered SGD optimization algorithm converges in the strong $L^p$-sense with order $\frac{1}{2}-\varepsilon$ to the global minimum of the objective function of the considered stochastic approximation problem under standard convexity-type assumptions on the objective function and relaxed assumptions on the moments of the stochastic errors appearing in the employed SGD optimization algorithm. The key ideas in our convergence proof are, first, to employ techniques from the theory of Lyapunov-type functions for dynamical systems to develop a general convergence machinery for SGD optimization algorithms based on such functions, then, to apply this general machinery to concrete Lyapunov-type functions with polynomial structures, and, thereafter, to perform an induction argument along the powers appearing in the Lyapunov-type functions in order to achieve for every arbitrarily large $ p \in (0,\infty) $ strong $ L^p $-convergence rates. This article also contains an extensive review of results on SGD optimization algorithms in the scientific literature.


abstract: We propose $S_4$ symmetric four-generation models for charged leptons. Although an $S_4$ symmetric four-generation model has been already proposed, there are some additional symmetries in the model. We construct four-generation models for charged leptons with only requirement of exact $S_4$ symmetry. It turned out that at least one of the models is consistent with observations of charged lepton masses and predicts the mass of the charged lepton of the fourth generation to be 556 GeV.


title: Sparse and Low-rank Tensor Estimation via Cubic Sketchings

abstract: In this paper, we propose a general framework for sparse and low-rank tensor estimation from cubic sketchings. A two-stage non-convex implementation is developed based on sparse tensor decomposition and thresholded gradient descent, which ensures exact recovery in the noiseless case and stable recovery in the noisy case with high probability. The non-asymptotic analysis sheds light on an interplay between optimization error and statistical error. The proposed procedure is shown to be rate-optimal under certain conditions. As a technical by-product, novel high-order concentration inequalities are derived for studying high-moment sub-Gaussian tensors. An interesting tensor formulation illustrates the potential application to high-order interaction pursuit in high-dimensional linear regression.


title: Configuration spaces of $\mathbb{C} \setminus k$

abstract: In this note, we collect mostly known formulas and methods to compute the standard and virtual Poincar\'e polynomials of the configuration spaces of the plane $\mathbb{C} \setminus k$ with $k$ deleted points and compare the answers.


title: A Comparison of Genetic Algorithms and Simulated Annealing in Maximizing the Thermal Conductivity of Discrete Massive Chains

abstract: Functions of chemical composition are complex and discrete in nature making it impossible to optimize them with gradient methods. Genetic algorithms, which do not use derivative information, are used to maximize the thermal conductivity of one-dimensional classical harmonic oscillators made from a fixed library of randomly generated molecular units. The ability for the genetic algorithm to build structures with components having no physical increment is important in optimizing molecules with a library of unrelated polymer units. The performance of genetic algorithms in this problem is compared with simulated annealing. Hyper-parameters for these routines are selected from a grid search in order to optimize their expected solution strength. The solutions found via the genetic algorithm consistently outperform those of simulated annealing at the cost of longer computer time. Together, these algorithms are able to find thermal conductor candidates that mirror solutions in continuous media.


title: Chemically stabilized epitaxial wurtzite-BN thin film

abstract: We report on the chemically stabilized epitaxial w-BN thin film grown on c-plane sapphire by pulsed laser deposition under slow kinetic condition. Traces of no other allotropes such as cubic (c) or hexagonal (h) BN phases are present. Sapphire substrate plays a significant role in stabilizing the metastable w-BN from h-BN target under unusual PLD growth condition involving low temperature and pressure and is explained based on density functional theory calculation. The hardness and the elastic modulus of the w-BN film are 37 & 339 GPa, respectively measured by indentation along <0001> direction. The results are extremely promising in advancing the microelectronic and mechanical tooling industry.