title: Exploiting NOMA for UAV Communications in Large-Scale Cellular Networks

abstract: This paper advocates a pair of strategies in non-orthogonal multiple access (NOMA) in unmanned aerial vehicles (UAVs) communications, where multiple UAVs play as new aerial communications platforms for serving terrestrial NOMA users. A new multiple UAVs framework with invoking stochastic geometry technique is proposed, in which a pair of practical strategies are considered: 1) the UAV-centric strategy for offloading actions and 2) the user-centric strategy for providing emergency communications. In order to provide practical insights for the proposed NOMA assisted UAV framework, an imperfect successive interference cancelation (ipSIC) scenario is taken into account. For both UAV-centric strategy and user-centric strategy, we derive new exact expressions for the coverage probability. We also derive new analytical results for orthogonal multiple access (OMA) for providing a benchmark scheme. The derived analytical results in both user-centric strategy and UAV-centric strategy explicitly indicate that the ipSIC coefficient is a dominant component in terms of coverage probability. Numerical results are provided to confirm that i) for both user-centric strategy and UAV-centric strategy, NOMA assisted UAV cellular networks is capable of outperforming OMA by setting power allocation factors and targeted rate properly; and ii) the coverage probability of NOMA assisted UAV cellular framework is affected to a large extent by ipSIC coefficient, target rates and power allocations factors of paired NOMA users.


title: Ideal zeta functions associated to a family of class-2-nilpotent Lie rings

abstract: We produce explicit formulae for various ideal zeta functions associated to the members of an infinite family of class-$2$-nilpotent Lie rings, introduced in [1], in terms of Igusa functions. As corollaries we obtain information about analytic properties of global ideal zeta functions, local functional equations, topological, reduced, and graded ideal zeta functions, as well as representation zeta functions for the unipotent group schemes associated to the Lie rings in question.


title: Colored five-vertex models and Demazure atoms

abstract: Type A Demazure atoms are pieces of Schur functions, or sets of tableaux whose weights sum to such functions. Inspired by colored vertex models of Borodin and Wheeler, we will construct solvable lattice models whose partition functions are Demazure atoms; the proof of this makes use of a Yang-Baxter equation for a colored five-vertex model. As a biproduct, we construct Demazure atoms on Kashiwara's $\mathcal{B}_\infty$ crystal and give new algorithms for computing Lascoux-Sch\"utzenberger keys.


title: Tipping Phenomena and Points of No Return in Ecosystems: Beyond Classical Bifurcations

abstract: We discuss tipping phenomena (critical transitions) in nonautonomous systems using an example of a bistable ecosystem model with environmental changes represented by time-varying parameters [Scheffer et al., Ecosystems, 11 (2008), pp. 275--279]. We give simple testable criteria for the occurrence of nonautonomous tipping from the herbivore-dominating equilibrium to the plant-only equilibrium using global properties of the autonomous frozen system with fixed-in-time parameters. To begin with, we use classical autonomous bifurcation analysis to identify a codimension-three degenerate Bogdanov-Takens bifurcation: the source of a dangerous subcritical Hopf bifurcation and the organizing center for bifurcation-induced tipping (B-tipping). Then, we introduce the concept of basin instability for equilibria to identify parameter paths along which genuine nonautonomous rate-induced tipping (R-tipping) occurs without crossing any classical autonomous bifurcations. We explain nonautonomous R-tipping in terms of maximal canard trajectories and produce nonautonomous tipping diagrams in the plane of the magnitude and rate of a parameter shift to uncover intriguing R-tipping tongues and wiggling tipping-tracking bifurcation curves. Discussion of nontrivial dynamics arising from the interaction between B-tipping and R-tipping identifies "points of no return" where tipping cannot be prevented by the parameter trend reversal and "points of return tipping" where tipping is inadvertently induced by the parameter trend reversal. Our results give new insight into the sensitivity of ecosystems to the magnitudes and rates of environmental change. Finally, a comparison between "tilted" saddle-node and subcritical Hopf normal forms reveals some universal tipping properties due to basin instability, a generic dangerous bifurcation, or the combination of both.


title: From Heegaard splittings to trisections; porting 3-dimensional ideas to dimension 4

abstract: These notes summarize and expand on a mini-course given at CIRM in February 2018 as part of Winter Braids VIII. We somewhat obsessively develop the slogan `Trisections are to 4-manifolds as Heegaard splittings are to 3-manifolds', focusing on and clarifying the distinction between three ways of thinking of things: the basic definitions as decompositions of manifolds, the Morse theoretic perspective and descriptions in terms of diagrams. We also lay out these themes in two important relative settings: 4-manifolds with boundary and 4-manifolds with embedded 2-dimensional submanifolds.


title: Nanoscale Substrate Roughness Hinders Domain Formation in Supported Lipid Bilayers

abstract: Supported Lipid Bilayers (SLBs) are model membranes formed at solid substrate surfaces. This architecture renders the membrane experimentally accessible to surface sensitive techniques used to study their properties, including Atomic Force Microscopy (AFM), optical fluorescence microscopy, Quartz Crystal Microbalance (QCM) and X-Ray/Neutron Reflectometry, and allows integration with technology for potential biotechnological applications such as drug screening devices. The experimental technique often dictates substrate choice or treatment, and it is anecdotally recognised that certain substrates are suitable for the particular experiment, but the exact influence of the substrate has not been comprehensively investigated. Here, we study the behavior of a simple model bilayer, phase separating on a variety of commonly used substrates, including glass, mica, silicon and quartz, with drastically different results. The distinct micron scale domains observed on mica, identical to those seen in free-floating Giant Unilamellar Vesicles (GUVs), are reduced to nanometer scale domains on glass and quartz. The mechanism for the arrest of domain formation is investigated, and the most likely candidate is nanoscale surface roughness, acting as a drag on the hydrodynamic motion of small domains during phase separation. Evidence was found that the physico-chemical properties of the surface have a mediating effect, most likely due to changes in the lubricating interstitial water layer between surface and bilayer.


title: Deep Convolutional Neural Network for Automated Detection of Mind Wandering using EEG Signals

abstract: Mind wandering (MW) is a ubiquitous phenomenon which reflects a shift in attention from task-related to task-unrelated thoughts. There is a need for intelligent interfaces that can reorient attention when MW is detected due to its detrimental effects on performance and productivity. In this paper, we propose a deep learning model for MW detection using Electroencephalogram (EEG) signals. Specifically, we develop a channel-wise deep convolutional neural network (CNN) model to classify the features of focusing state and MW extracted from EEG signals. This is the first study that employs CNN to automatically detect MW using only EEG data. The experimental results on the collected dataset demonstrate promising performance with 91.78% accuracy, 92.84% sensitivity, and 90.73% specificity.


title: Empathic Robot for Group Learning: A Field Study

abstract: This work explores a group learning scenario with an autonomous empathic robot. We address two research questions: (1) Can an autonomous robot designed with empathic competencies foster collaborative learning in a group context? (2) Can an empathic robot sustain positive educational outcomes in long-term collaborative learning interactions with groups of students? To answer these questions, we developed an autonomous robot with empathic competencies that is able to interact with a group of students in a learning activity about sustainable development. Two studies were conducted. The first study compares learning outcomes in children across 3 conditions: learning with an empathic robot; learning with a robot without empathic capabilities; and learning without a robot. The results show that the autonomous robot with empathy fosters meaningful discussions about sustainability, which is a learning outcome in sustainability education. The second study features groups of students who interact with the robot in a school classroom for two months. The long-term educational interaction did not seem to provide significant learning gains, although there was a change in game-actions to achieve more sustainability during game-play. This result reflects the need to perform more long-term research in the field of educational robots for group learning.


title: C^1 Deformations of almost-Grassmannian structures with strongly essential symmetry

abstract: We construct a family of $(2,n)$-almost Grassmannian structures of regularity $C^1$, each admitting a one-parameter group of strongly essential automorphisms, and each not flat on any neighborhood of the higher-order fixed point. This shows that Theorem 1.3 of [9] does not hold assuming only $C^1$ regularity of the structure (see also [2, Prop 3.5]).


title: How should you discount your backtest PnL?

abstract: In-sample overfitting is a drawback of any backtest-based investment strategy. It is thus of paramount importance to have an understanding of why and how the in-sample overfitting occurs. In this article we propose a simple framework that allows one to model and quantify in-sample PnL overfitting. This allows us to compute the factor appropriate for discounting PnLs of in-sample investment strategies.


title: Entropy formula of N-body system

abstract: We prove a proposition that the entropy of the system composed of finite $N$ molecules of ideal gas is the $q$-entropy or Havrda-Charv\'at-Tsallis entropy, which is also known as Tsallis entropy, with the entropic index $q=\frac{D(N-1)-4}{D(N-1)-2}$ in $D$-dimensional space. The indispensable infinity assumption used by Boltzmann and others in their derivation of entropy formulae is not involved in our derivation, therefore our derived formula is exact. The analogy of the $N$-body system brings us to obtain the entropic index of a combined system $q_C$ formed from subsystems having different entropic indexes $q_A$ and $q_B$ as $\frac{1}{1-q_C}=\frac{1}{1-q_A}+\frac{1}{1-q_B}+\frac{D+2}{2}$. It is possible to use the number $N$ for the physical measure of deviation from Boltzmann entropy.


title: The suppressed radiative recombination rate in a quantum photocell with three electron donors

abstract: The radiative recombination of electron-hole pairs represents a great challenge to the photon-to-charge efficiency in the photocell. In this paper, we investigate how to suppress radiative recombination rate (RRR) in a proposed quantum photocell with three dipole-dipole coupled and uncoupled electron donors. The results showed that the RRR could be suppressed in this photocell with three uncoupled electron donors but be enhanced with three dipole-dipole coupled electron donors by the ambient circumstance temperatures, and the increasing energy gap in the donors, the decreasing gap between the donors and acceptor inhabited the RRR with three dipole-dipole both coupled and uncoupled electron donors. When the photocell was manipulated by the electrostatic dipole-dipole coupling strength J at room temperature, the RRR was suppressed to a smaller minimum by the gap between the donors and acceptor than those by different gaps in the donors. These suppressed strategies for RRR point out some significant ways to increase the photon-to-charge efficiency and deserve the further experimental verification.


title: A classical field theory formulation for the numerical solution of time harmonic electromagnetic fields

abstract: Finite element representations of Maxwell's equations pose unusual challenges inherent to the variational representation of the `curl-curl' equation for the fields. We present a variational formulation based on classical field theory. Borrowing from QED, we modify the Lagrangian by adding an implicit gauge-fixing term. Our formulation, in the language of differential geometry, shows that conventional edge elements should be replaced by the simpler nodal elements for time-harmonic problems. We demonstrate how this formulation, adhering to the deeper underlying symmetries of the four-dimensional covariant field description, provides a highly general, robust numerical framework.


title: Helium enrichment in intermediate-age Magellanic Clouds clusters: towards an ubiquity of multiple stellar populations?

abstract: Intermediate-age star clusters in the Magellanic Clouds harbour signatures of the multiple stellar populations long thought to be restricted to old globular clusters. We compare synthetic horizontal branch models with Hubble Space Telescope photometry of clusters in the Magellanic Clouds, with age between ~2 and ~10 Gyr, namely NGC 121, Lindsay 1, NGC 339, NGC 416, Lindsay 38, Lindsay 113, Hodge 6 and NGC 1978. We find a clear signature of initial helium abundance spreads (delta(Y)) in four out of these eight clusters (NGC 121, Lindsay 1, NGC 339, NGC 416) and we quantify the value of delta(Y). For two clusters (Lindsay 38, Lindsay 113) we can only determine an upper limit for delta(Y), whilst for the two youngest clusters in our sample (Hodge 6 and NGC 1978) no conclusion about the existence of an initial He spread can be reached. Our delta(Y) estimates are consistent with the correlation between maximum He abundance spread and mass of the host cluster found in Galactic globular clusters. This result strengthens the emerging view that the formation of multiple stellar populations is a standard process in massive star clusters, not limited to a high redshift environment.


title: Communication Enhancement Through Quantum Coherent Control of $N$ Channels in an Indefinite Causal-order Scenario

abstract: In quantum Shannon theory, transmission of information is enhanced by quantum features. Up to very recently, the trajectories of transmission remained fully classical. Recently, a new paradigm was proposed by playing quantum tricks on two completely depolarizing quantum channels i.e. using coherent control in space or time of the two quantum channels. We extend here this control to the transmission of information through a network of an arbitrary number $N$ of channels with arbitrary individual capacity i.e. information preservation characteristics in the case of indefinite causal order. We propose a formalism to assess information transmission in the most general case of $N$ channels in an indefinite causal order scenario yielding the output of such transmission. Then we explicitly derive the quantum switch output and the associated Holevo limit of the information transmission for $N=2$, $N=3$ as a function of all involved parameters. We find in the case $N=3$ that the transmission of information for three channels is twice of transmission of the two channel case when a full superposition of all possible causal orders is used.


title: A Bootstrap Test for the Existence of Moments for GARCH Processes

abstract: This paper studies the joint inference on conditional volatility parameters and the innovation moments by means of bootstrap to test for the existence of moments for GARCH(p,q) processes. We propose a residual bootstrap to mimic the joint distribution of the quasi-maximum likelihood estimators and the empirical moments of the residuals and also prove its validity. A bootstrap-based test for the existence of moments is proposed, which provides asymptotically correctly-sized tests without losing its consistency property. It is simple to implement and extends to other GARCH-type settings. A simulation study demonstrates the test's size and power properties in finite samples and an empirical application illustrates the testing approach.


title: A Note on the Modified Albertson Index

abstract: The modified Albertson index, denoted by $A\!^*\!$, of a graph $G$ is defined as $A\!^*\!(G)=\sum_{uv\in E(G)} |(d_{u})^{2}- (d_{v})^{2}|$, where $d_u$, $d_v$ denote the degrees of the vertices $u$, $v$, respectively, of $G$ and $E(G)$ is the edge set of $G$. In this note, a sharp lower bound of $A\!^*$ in terms of the maximum degree for the case of trees is derived. The $n$-vertex trees having maximal and minimal $A\!^*$ values are also characterized here. Moreover, it is shown that $A\!^*\!(G)$ is non-negative even integer for every graph $G$ and that there exist infinitely many connected graphs whose $A\!^*$ value is $2t$ for every integer $t\in\{0,3,4,5\}\cup\{8,9,10,\cdots\}$.


title: Exact Markov Chain-based Runtime Analysis of a Discrete Particle Swarm Optimization Algorithm on Sorting and OneMax

abstract: Meta-heuristics are powerful tools for solving optimization problems whose structural properties are unknown or cannot be exploited algorithmically. We propose such a meta-heuristic for a large class of optimization problems over discrete domains based on the particle swarm optimization (PSO) paradigm. We provide a comprehensive formal analysis of the performance of this algorithm on certain "easy" reference problems in a black-box setting, namely the sorting problem and the problem OneMax. In our analysis we use a Markov model of the proposed algorithm to obtain upper and lower bounds on its expected optimization time. Our bounds are essentially tight with respect to the Markov model. We show that for a suitable choice of algorithm parameters the expected optimization time is comparable to that of known algorithms and, furthermore, for other parameter regimes, the algorithm behaves less greedy and more explorative, which can be desirable in practice in order to escape local optima. Our analysis provides a precise insight on the tradeoff between optimization time and exploration. To obtain our results we introduce the notion of indistinguishability of states of a Markov chain and provide bounds on the solution of a recurrence equation with non-constant coefficients by integration.


title: Melt-Blowing of Viscoelastic Jets in Turbulent Airflows: Stochastic Modeling and Simulation

abstract: In melt-blowing processes mico- and nanofibers are produced by the extrusion of polymeric jets into a directed, turbulent high-speed airflow. Up to now the physical mechanism for the drastic jet thinning is not fully understood, since in the existing literature the numerically computed/predicted fiber thickness differs several orders of magnitude from those experimentally measured. Recent works suggest that this discrepancy might arise from the neglect of the turbulent aerodynamic fluctuations in the simulations. In this paper we confirm this suggestion numerically. Due to the complexity of the process direct numerical simulations of the multiscale-multiphase problem are not possible. Hence, we develop a numerical framework for a growing fiber in turbulent air that makes the simulation of industrial setups feasible. For this purpose we employ an asymptotic viscoelastic model for the fiber. The turbulent effects are taken into account by a stochastic aerodynamic force model where the underlying velocity fluctuations are reconstructed from a $k$-$\epsilon$ turbulence description of the airflow. Our numerical results show the significance of the turbulence on the jet thinning and give fiber diameters of realistic order of magnitude.


title: Error Rate Analysis of Amplitude-Coherent Detection over Rician Fading Channels with Receiver Diversity

abstract: Amplitude-coherent (AC) detection is an efficient detection technique that can simplify the receiver design while providing reliable symbol error rate (SER). Therefore, this work considers AC detector design and SER analysis using M-ary amplitude shift keying (MASK) modulation over Rician fading channels. More specifically, we derive the optimum, near-optimum and a suboptimum AC detectors and compare their SER to the coherent, noncoherent and the heuristic AC detectors. Moreover, the analytical SER of the heuristic detector is derived using two different approaches for single and multiple receiving antennas. One of the derived expressions is expressed in terms of a single integral that can be evaluated numerically, while the second approach gives a closed-form analytical expression for the SER, which is also used to derive a simple formula for the asymptotic SER at high signal-to-noise ratios (SNRs). The obtained analytical and simulation results show that the SER of the AC and coherent MASK detectors are comparable, particularly for high values of the Rician K-factor, and small number of receiving antennas. Moreover, the obtained results show that the SER of the optimal AC detector is equivalent to that of the coherent detector. However, the optimal AC detector complexity is prohibitively high, particularly at high SNRs. In most of the scenarios, the heuristic AC detector significantly outperforms the optimum noncoherent detector, except for the binary ASK case at low SNRs. Moreover, the obtained results show that the heuristic AC detector is immune to phase noise, and thus, it outperforms the coherent detector in scenarios where system is subject to considerable phase noise.


title: Modular Block-diagonal Curvature Approximations for Feedforward Architectures

abstract: We propose a modular extension of backpropagation for the computation of block-diagonal approximations to various curvature matrices of the training objective (in particular, the Hessian, generalized Gauss-Newton, and positive-curvature Hessian). The approach reduces the otherwise tedious manual derivation of these matrices into local modules, and is easy to integrate into existing machine learning libraries. Moreover, we develop a compact notation derived from matrix differential calculus. We outline different strategies applicable to our method. They subsume recently-proposed block-diagonal approximations as special cases, and are extended to convolutional neural networks in this work.


title: A non-iterative method for robustly computing the intersections between a line and a curve or surface

abstract: The need to compute the intersections between a line and a high-order curve or surface arises in a large number of finite element applications. Such intersection problems are easy to formulate but hard to solve robustly. We introduce a non-iterative method for computing intersections by solving a matrix singular value decomposition (SVD) and an eigenvalue problem. That is, all intersection points and their parametric coordinates are determined in one-shot using only standard linear algebra techniques available in most software libraries. As a result, the introduced technique is far more robust than the widely used Newton-Raphson iteration or its variants. The maximum size of the considered matrices depends on the polynomial degree $q$ of the shape functions and is $2q \times 3q$ for curves and $6 q^2 \times 8 q^2$ for surfaces. The method has its origin in algebraic geometry and has here been considerably simplified with a view to widely used high-order finite elements. In addition, the method is derived from a purely linear algebra perspective without resorting to algebraic geometry terminology. A complete implementation is available from http://bitbucket.org/nitro-project/.


title: Effects of dark matter in star formation

abstract: The standard model for the formation of structure assumes that there existed small fluctuations in the early universe that grew due to gravitational instability. The origins of these fluctuations are as yet unclear. In this work we propose the role of dark matter in providing the seed for star formation in the early universe. Very recent observations also support the role of dark matter in the formation of these first stars. With this we set observable constraints on luminosities, temperatures, and lifetimes of these early stars with an admixture of dark matter.


title: Conformal wave equations for the Einstein-tracefree matter system

abstract: Inspired by a similar analysis for the vacuum conformal Einstein field equations by Paetz [Ann. H. Poincar\'e 16, 2059 (2015)], in this article we show how to construct a system of quasilinear wave equations for the geometric fields associated to the conformal Einstein field equations coupled to matter models whose energy-momentum tensor has vanishing trace. In this case, the equation of conservation for the energy-momentum tensor is conformally invariant. Our analysis includes the construction of a subsidiary evolution system which allows to prove the propagation of the constraints. We discuss how the underlying structure behind these systems of equations is the integrability conditions satisfied by the conformal field equations. The main result of our analysis is that both the evolution and subsidiary equations for the geometric part of the conformal Einstein-tracefree matter field equations close without the need of any further assumption on the matter models other than the vanishing of the trace of the energy-momentum tensor. Our work is supplemented by an analysis of the evolution and subsidiary equations associated to three basic tracefree matter models: the conformally invariant scalar field, the Maxwell field and the Yang-Mills field. As an application we provide a global existence and stability result for de Sitter-like spacetimes. In particular, the result for the conformally coupled scalar field is new in the literature.


title: MIMO Capacity with Average Total and Per-Antenna Power Constraints

abstract: MIMO capacity with a joint total and per-antenna average power constraint is considered in this work. The problem arises when, besides having a limited available power at the transmitter, also the individual antennas cannot radiate power beyond the limits of their corresponding RF chains. Closed-form results are illustrated in specific cases and, in particular, in the unit-rank channel matrix case. Lower-complexity optimization problems are derived for the other cases. Numerical complexity is derived in terms of the number of equivalent real variables for the optimization problem and a general formula is provided depending on the channel matrix rank. Numerical results are included to validate and illustrate the application of the proposed optimization algorithms and also to evaluate the time complexity of its implementation.


title: Fast multigrid solvers for conforming and non-conforming multi-patch Isogeometric Analysis

abstract: Isogeometric Analysis is a high-order discretization method for boundary value problems that uses a number of degrees of freedom which is as small as for a low-order method. Standard isogeometric discretizations require a global parameterization of the computational domain. In non-trivial cases, the domain is decomposed into patches having separate parameterizations and separate discretization spaces. If the discretization spaces agree on the interfaces between the patches, the coupling can be done in a conforming way. Otherwise, non-conforming discretizations (utilizing discontinuous Galerkin approaches) are required. The author and his coworkers have previously introduced multigrid solvers for Isogeometric Analysis for the conforming case. In the present paper, these results are extended to the non-conforming case. Moreover, it is shown that the multigrid solves get even more powerful if the proposed smoother is combined with a (standard) Gauss-Seidel smoother.


title: Polynomial Invariants, Knot Homologies, and Higher Twist Numbers of Weaving Knots $W(3,n)$

abstract: We investigate several conjectures in geometric topology by assembling computer data obtained by studying weaving knots, a doubly infinite family $W(p,n)$ of examples of hyperbolic knots. In particular, we compute some important polynomial knot invariants, as well as knot homologies, for the subclass $W(3,n)$ of this family. We use these knot invariants to conclude that all knots $W(3,n)$ are fibered knots and provide estimates for some geometric invariants of these knots. Finally, we study the asymptotics of the ranks of their Khovanov homology groups. Our investigations provide evidence for our conjecture that, asymptotically as $n$ grows large, the ranks of Khovanov homology groups of $W(3,n)$ are normally distributed.


title: Ultra-recursive sequences

abstract: We study a new type of sequences whose elements are defined in terms of the position, sign and magnitude of another element of the sequence. The name ultra-recursive comes from the fact that these sequences possess terms that are generated adding either the previous (as in traditional recurrences formulas) or the next terms. They are also self-referential because the rules to generate every member alludes to some value of the sequence.


title: A Study of an Agile Methodology with Scrum Approach to the Filipino Company-Sponsored I.T. Capstone Program

abstract: Purpose - The research aims to show the relevance of company client sponsored student projects in the University of Asia and the Pacific Information Technology (UA&P IT) Capstone Program through the use ofan Agile Methodology with Scrum Approach. Method - The modified program is employed on two batches with content analysis and survey results as benchmarks. Results - Surveys at the end of the sprints for both clients and students revealed that the length of the sprint was a critical factor in the development of the information system, and that students learned from addressing additional challenges such as academic load, team pressure and communication issues. Conclusion - Over-all results showed that clients were impressed and keen to adopt the student works. Recommendations - Maintainability aspects of the research can be analyzed for future studies. Increasing the sample size with additional batches could lead to discovery of additional factors not previously seen. Research Implications - The research could help improve other Capstone Programs while improving communication with company clients.


title: Do We Preach What We Practice? Investigating the Practical Relevance of Requirements Engineering Syllabi - The IREB Case

abstract: Nowadays, there exist a plethora of different educational syllabi for Requirements Engineering (RE), all aiming at incorporating practically relevant educational units (EUs). Many of these syllabi are based, in one way or the other, on the syllabi provided by the International Requirements Engineering Board (IREB), a non-profit organisation devoted to standardised certification programs for RE. IREB syllabi are developed by RE experts and are, thus, based on the assumption that they address topics of practical relevance. However, little is known about to what extent practitioners actually perceive those contents as useful. We have started a study to investigate the relevance of the EUs included in the IREB Foundation Level certification programme. In a first phase reported in this paper, we have surveyed practitioners mainly from DACH countries (Germany, Austria and Switzerland) participating in the IREB certification. Later phases will widen the scope both by including other countries and by not requiring IREB-certified participants. The results shall foster a critical reflection on the practical relevance of EUs built upon the de-facto standard syllabus of IREB.


title: 2-universality in randomly perturbed graphs

abstract: A graph $G$ is called universal for a family of graphs $\mathcal{F}$ if it contains every element $F \in \mathcal{F}$ as a subgraph. Let $\mathcal{F}(n,2)$ be the family of all graphs with maximum degree $2$. Ferber, Kronenberg, and Luh [Optimal Threshold for a Random Graph to be 2-Universal, to appear in Transactions of the American Mathematical Society] proved that there exists a $C$ such that for $p \ge C (n^{-2/3} \log^{1/3} n )$ the random graph $G(n,p)$ a.a.s is $\mathcal{F}(n,2)$-universal, which is asymptotically optimal. For any $n$-vertex graph $G_\alpha$ with minimum degree $\delta(G_\alpha) \ge \alpha n$ Aigner and Brandt [Embedding arbitrary graphs of maximum degree two, Journal of the London Mathematical Society 48 (1993), 39-51] proved that $G_\alpha$ is $\mathcal{F}(n,2)$-universal for an optimal $\alpha \ge 2/3$. In this note, we consider the model of randomly perturbed graphs, which is the union $G_\alpha \cup G(n,p)$. We prove that $G_\alpha \cup G(n,p)$ is a.a.s. $\mathcal{F}(n,2)$-universal provided that $\alpha>0$ and $p=\omega(n^{-2/3})$. This is asymptotically optimal and improves on both results from above in the respective parameter. Furthermore, this extends a result of B\"ottcher, Montgomery, Parczyk, and Person [Embedding spanning bounded degree subgraphs in randomly perturbed graphs, arXiv:1802.04603 (2018)], who embed a given $F \in \mathcal{F}(n,2)$ at these values. We also prove variants with universality for the family $\mathcal{F}^\ell(n,2)$, all graphs from $\mathcal{F}(n,2)$ with girth at least $\ell$. For example, there exists an $\ell_0$ depending only on $\alpha$ such that for all $\ell \ge \ell_0$ already $p=\omega(1/n)$ is sufficient for $\mathcal{F}^\ell(n,2)$-universality.


title: Deep Convolutional Generative Adversarial Networks Based Flame Detection in Video

abstract: Real-time flame detection is crucial in video based surveillance systems. We propose a vision-based method to detect flames using Deep Convolutional Generative Adversarial Neural Networks (DCGANs). Many existing supervised learning approaches using convolutional neural networks do not take temporal information into account and require substantial amount of labeled data. In order to have a robust representation of sequences with and without flame, we propose a two-stage training of a DCGAN exploiting spatio-temporal flame evolution. Our training framework includes the regular training of a DCGAN with real spatio-temporal images, namely, temporal slice images, and noise vectors, and training the discriminator separately using the temporal flame images without the generator. Experimental results show that the proposed method effectively detects flame in video with negligible false positive rates in real-time.


title: Long-range magnetic ordering in rocksalt-type high-entropy oxides

abstract: We report the magnetic properties of Mg$_{0.2}$Co$_{0.2}$Ni$_{0.2}$Cu$_{0.2}$Zn$_{0.2}$O, a high-entropy oxide with rocksalt structure, and the influence of substitutions on these properties. From the magnetic susceptibility and neutron diffraction measurements, we found that this compound exhibits long-range magnetic order below 120 K despite the substantial structuraldisorder. The other rocksalt-type high-entropy oxides with various chemical substitutions were found to host either an antiferromagnetic order or spin-glass state depending on the amount of magnetic ions. The presence of magnetic order for such a disordered material potentially provide a route to explore novel magnetic properties and functions.


title: The effects of a revised $^7$Be e$^-$-capture rate on solar neutrino fluxes

abstract: The electron-capture rate on $^7$Be is the main production channel for $^7$Li in several astrophysical environments. Theoretical evaluations have to account for not only the nuclear interaction, but also the processes in the plasma where $^7$Be ions and electrons interact. In the past decades several estimates were presented, pointing out that the theoretical uncertainty in the rate is in general of few percents. In the framework of fundamental solar physics, we consider here a recent evaluation for the $^7$Be+e$^-$ rate, not used up to now in the estimate of neutrino fluxes. We analysed the effects of the new assumptions on Standard Solar Models (SSMs) and compared the results obtained by adopting the revised $^7$Be+e$^-$ rate to those obtained by the one reported in a widely used compilation of reaction rates (ADE11). We found that new SSMs yield a maximum difference in the efficiency of the $^7$Be channel of about -4\% with respect to what is obtained with the previously adopted rate. This fact affects the production of neutrinos from $^8$B, increasing the relative flux up to a maximum of 2.7\%. Negligible variations are found for the physical and chemical properties of the computed solar models. The agreement with the SNO measurements of the neutral current component of the $^8$B neutrino flux is improved.


title: Interactive molecular dynamics in virtual reality from quantum chemistry to drug binding: An open-source multi-person framework

abstract: As molecular scientists have made progress in their ability to engineer nano-scale molecular structure, we are facing new challenges in our ability to engineer molecular dynamics (MD) and flexibility. Dynamics at the molecular scale differs from the familiar mechanics of everyday objects, because it involves a complicated, highly correlated, and three-dimensional many-body dynamical choreography which is often non-intuitive even for highly trained researchers. We recently described how interactive molecular dynamics in virtual reality (iMD-VR) can help to meet this challenge, enabling researchers to manipulate real-time MD simulations of flexible structures in 3D. In this article, we outline efforts to extend immersive technologies to the molecular sciences, and we introduce 'Narupa', a flexible, open-source, multi-person iMD-VR software framework which enables groups of researchers to simultaneously cohabit real-time simulation environments to interactively visualize and manipulate the dynamics of molecular structures with atomic-level precision. We outline several application domains where iMD-VR is facilitating research, communication, and creative approaches within the molecular sciences, including training machines to learn reactive potential energy surfaces (PESs), biomolecular conformational sampling, protein-ligand binding, reaction discovery using 'on-the-fly' quantum chemistry, and transport dynamics in materials. We touch on iMD-VR's various cognitive and perceptual affordances, and how these provide research insight for molecular systems. By synergistically combining human spatial reasoning and design insight with computational automation, technologies like iMD-VR have the potential to improve our ability to understand, engineer, and communicate microscopic dynamical behavior, offering the potential to usher in a new paradigm for engineering molecules and nano-architectures.

categories: physics.chem-ph cs.HC physics.bio-ph physics.ed-ph


title: Skew-symmetric entropy stable modal discontinuous Galerkin formulations

abstract: High order entropy stable discontinuous Galerkin (DG) methods for nonlinear conservation laws satisfy an inherent discrete entropy inequality. The construction of such schemes has relied on the use of carefully chosen nodal points or volume and surface quadrature rules to produce operators which satisfy a summation-by-parts (SBP) property. In this work, we show how to construct "modal" skew-symmetric DG formulations which are entropy stable for volume and surface quadratures under which a traditional SBP property does not hold. These skew-symmetric formulations avoid the use of a "strong" matrix-based SBP property, and instead rely on the assumption that discrete operators exactly differentiate constants and satisfy a discrete form of the fundamental theorem of calculus. We conclude with numerical experiments verifying the accuracy and stability of the proposed formulations, and discuss an application of skew-symmetric formulations for entropy stable DG schemes on mixed quadrilateral-triangle meshes.


title: Hierarchical Matrix Operations on GPUs: Matrix-Vector Multiplication and Compression

abstract: Hierarchical matrices are space and time efficient representations of dense matrices that exploit the low rank structure of matrix blocks at different levels of granularity. The hierarchically low rank block partitioning produces representations that can be stored and operated on in near-linear complexity instead of the usual polynomial complexity of dense matrices. In this paper, we present high performance implementations of matrix vector multiplication and compression operations for the $\mathcal{H}^2$ variant of hierarchical matrices on GPUs. This variant exploits, in addition to the hierarchical block partitioning, hierarchical bases for the block representations and results in a scheme that requires only $O(n)$ storage and $O(n)$ complexity for the mat-vec and compression kernels. These two operations are at the core of algebraic operations for hierarchical matrices, the mat-vec being a ubiquitous operation in numerical algorithms while compression/recompression represents a key building block for other algebraic operations, which require periodic recompression during execution. The difficulties in developing efficient GPU algorithms come primarily from the irregular tree data structures that underlie the hierarchical representations, and the key to performance is to recast the computations on flattened trees in ways that allow batched linear algebra operations to be performed. This requires marshaling the irregularly laid out data in a way that allows them to be used by the batched routines. Marshaling operations only involve pointer arithmetic with no data movement and as a result have minimal overhead. Our numerical results on covariance matrices from 2D and 3D problems from spatial statistics show the high efficiency our routines achieve---over 550GB/s for the bandwidth-limited mat-vec and over 850GFLOPS/s in sustained performance for the compression on the P100 Pascal GPU.


title: Advancing the matter bispectrum estimation of large-scale structure: a comparison of dark matter codes

abstract: Cosmological information from forthcoming galaxy surveys, such as LSST and Euclid, will soon exceed that available from the CMB. Higher order correlation functions, like the bispectrum, will be indispensable for realising this potential. The interpretation of this data faces many challenges because gravitational collapse of matter is a complex non-linear process, typically modelled by computationally expensive N-body simulations. Proposed alternatives using fast dark matter codes (e.g. 2LPT or particle-mesh) are primarily evaluated on their ability to reproduce clustering statistics linked to the matter power spectrum. The accuracy of these codes can be tested in more detail by looking at higher-order statistics, and in this paper we will present an efficient and optimal methodology (MODAL-LSS) to reconstruct the full bispectrum of any 3D density field. We make quantitative comparisons between a number of fast dark matter codes and Gadget at redshift $z=0.5$. This will serve as an important diagnostic tool for dark matter/halo mock catalogues and lays the foundation for realistic high precision analysis with the galaxy bispectrum. In particular, we show that the lack of small-scale power in the bispectrum of fast codes can be ameliorated by a simple `boosting' technique for the power spectrum. We also investigate the covariance of the MODAL-LSS bispectrum estimator, demonstrating the plateauing of non-Gaussian errors in contrast to simple Gaussian extrapolations. This has important consequences for the extraction of information from the bispectrum and hence parameter estimation. Finally we make quantitative comparisons of simulation bispectra with theoretical models, discussing the initial parameters required to create mock catalogues with accurate bispectra.


title: Face Alignment using a 3D Deeply-initialized Ensemble of Regression Trees

abstract: Face alignment algorithms locate a set of landmark points in images of faces taken in unrestricted situations. State-of-the-art approaches typically fail or lose accuracy in the presence of occlusions, strong deformations, large pose variations and ambiguous configurations. In this paper we present 3DDE, a robust and efficient face alignment algorithm based on a coarse-to-fine cascade of ensembles of regression trees. It is initialized by robustly fitting a 3D face model to the probability maps produced by a convolutional neural network. With this initialization we address self-occlusions and large face rotations. Further, the regressor implicitly imposes a prior face shape on the solution, addressing occlusions and ambiguous face configurations. Its coarse-to-fine structure tackles the combinatorial explosion of parts deformation. In the experiments performed, 3DDE improves the state-of-the-art in 300W, COFW, AFLW and WFLW data sets. Finally, we perform cross-dataset experiments that reveal the existence of a significant data set bias in these benchmarks.


title: Inferring the strength of social ties: a community-driven approach

abstract: Online social networks are growing and becoming denser. The social connections of a given person may have very high variability: from close friends and relatives to acquaintances to people who hardly know. Inferring the strength of social ties is an important ingredient for modeling the interaction of users in a network and understanding their behavior. Furthermore, the problem has applications in computational social science, viral marketing, and people recommendation. In this paper we study the problem of inferring the strength of social ties in a given network. Our work is motivated by a recent approach [27], which leverages the strong triadic closure (STC) principle, a hypothesis rooted in social psychology [13]. To guide our inference process, in addition to the network structure, we also consider as input a collection of tight communities. Those are sets of vertices that we expect to be connected via strong ties. Such communities appear in different situations, e.g., when being part of a community implies a strong connection to one of the existing members. We consider two related problem formalizations that reflect the assumptions of our setting: small number of STC violations and strong-tie connectivity in the input communities. We show that both problem formulations are NP-hard. We also show that one problem formulation is hard to approximate, while for the second we develop an algorithm with approximation guarantee. We validate the proposed method on real-world datasets by comparing with baselines that optimize STC violations and community connectivity separately.


title: Flat affine symplectic Lie groups

abstract: We give a new characterization of flat affine manifolds in terms of an action of the Lie algebra of classical infinitesimal affine transformations on the bundle of linear frames. We characterize flat affine symplectic Lie groups using symplectic \'etale affine representations and as a consequence of this, we show that a flat affine symplectic Lie group with bi-invariant symplectic connection contains a nontrivial one parameter subgroup formed by central translations. We give two methods for constructing flat affine symplectic Lie groups, thus obtaining all those having bi-invariant symplectic connections. We get nontrivial examples of simply connected flat affine symplectic Lie groups in every even dimension.


title: Mean curvature versus diameter and energy quantization

abstract: We first partially extend a theorem of Topping, on the relation between mean curvature and intrinsic diameter, from immersed submanifolds of $\mathbb{R} ^{n} $ to almost everywhere immersed, closed submanifolds of a compact Riemannian manifold. We use this to prove quantization of energy for pseudo-holomorphic closed curves, of all genus, in a compact locally conformally symplectic manifold.


title: Fano 4-folds with rational fibrations

abstract: We study (smooth, complex) Fano 4-folds X having a rational contraction of fiber type, that is, a rational map X-->Y that factors as a sequence of flips followed by a contraction of fiber type. The existence of such a map is equivalent to the existence of a non-zero, non-big movable divisor on X. Our main result is that if Y is not P^1 or P^2, then the Picard number rho(X) of X is at most 18, with equality only if X is a product of surfaces. We also show that if a Fano 4-fold X has a dominant rational map X-->Z, regular and proper on an open subset of X, with dim(Z)=3, then either X is a product of surfaces, or rho(X) is at most 12. These results are part of a program to study Fano 4-folds with large Picard number via birational geometry.


title: Optimal common resource in majorization-based resource theories

abstract: We address the problem of finding the optimal common resource for an arbitrary family of target states in quantum resource theories based on majorization, that is, theories whose conversion law between resources is determined by a majorization relationship, such as it happens with entanglement, coherence or purity. We provide a conclusive answer to this problem by appealing to the completeness property of the majorization lattice. We give a proof of this property that relies heavily on the more geometric construction provided by the Lorenz curves, which allows to explicitly obtain the corresponding infimum and supremum. Our framework includes the case of possibly non-denumerable sets of target states (i.e. targets sets described by continuous parameters). In addition, we show that a notion of approximate majorization, which has recently found application in quantum thermodynamics, is in close relation with the completeness of this lattice. Finally, we provide some examples of optimal common resources within the resource theory of quantum coherence.


title: Neutron oscillations for solving neutron lifetime and dark matter puzzles

abstract: A model of $n-n'$ (neutron-mirror neutron) oscillations is proposed under the framework of the mirror matter theory with slightly broken mirror symmetry. It resolves the neutron lifetime discrepancy, i.e., the 1% difference in neutron lifetime between measurements from "beam" and "bottle" experiments. In consideration of the early universe evolution, the $n-n'$ mass difference is determined to be about $2\times 10^{-6}$ eV/c$^2$ with the $n-n'$ mixing strength of about $2\times 10^{-5}$. The picture of how the mirror-to-ordinary matter density ratio is evolved in the early universe into the observed dark-to-baryon matter density ratio of about 5.4 is presented. Reanalysis of previous data and new experiments that can be carried out under current technology are discussed and recommended to test this proposed model. Other consequences of the model on astrophysics and possible oscillations of other neutral particles are discussed as well.


title: How to "DODGE" Complex Software Analytics?

abstract: Machine learning techniques applied to software engineering tasks can be improved by hyperparameter optimization, i.e., automatic tools that find good settings for a learner's control parameters. We show that such hyperparameter optimization can be unnecessarily slow, particularly when the optimizers waste time exploring "redundant tunings"', i.e., pairs of tunings which lead to indistinguishable results. By ignoring redundant tunings, DODGE, a tuning tool, runs orders of magnitude faster, while also generating learners with more accurate predictions than seen in prior state-of-the-art approaches.


title: A Finite Volume Scheme for the Solution of a Mixed Discrete-Continuous Fragmentation Model

abstract: This paper concerns the construction and analysis of a numerical scheme for a mixed discrete-continuous fragmentation equation. A finite volume scheme is developed, based on a conservative formulation of a truncated version of the equations. The approximate solutions provided by this scheme are first shown to display conservation of mass and preservation of nonnegativity. Then, by utilising a Dunford--Pettis style argument, the sequence of approximate solutions generated is shown, under given restrictions on the model and the mesh, to converge (weakly) in an appropriate $L_1$ space to a weak solution to the problem. Additionally, by applying the methods and theory of operator semigroups, we are further able to show that weak solutions to the problem are unique and necessarily classical (differentiable) solutions. Finally, numerical simulations are performed to investigate the performance of the scheme and assess its rate of convergence.


title: A Note on the Subleading Soft Graviton

abstract: We show that the soft part of the charge generating infinitesimal superrotations can be expressed, in harmonic gauge, in terms of metric components evaluated at the boundaries of null infinity that are subleading in a large radius expansion. We then recast the spin memory observable in terms of these boundary values.


title: Seasonal evolution of temperatures in Titan's lower stratosphere

abstract: The Cassini mission offered us the opportunity to monitor the seasonal evolution of Titan's atmosphere from 2004 to 2017, i.e. half a Titan year. The lower part of the stratosphere (pressures greater than 10 mbar) is a region of particular interest as there are few available temperature measurements, and because its thermal response to the seasonal and meridional insolation variations undergone by Titan remains poorly known. In this study, we measure temperatures in Titan's lower stratosphere between 6 mbar and 25 mbar using Cassini/CIRS spectra covering the whole duration of the mission (from 2004 to 2017) and the whole latitude range. We can thus characterize the meridional distribution of temperatures in Titan's lower stratosphere, and how it evolves from northern winter (2004) to summer solstice (2017). Our measurements show that Titan's lower stratosphere undergoes significant seasonal changes, especially at the South pole, where temperature decreases by 19 K at 15 mbar in 4 years.


title: Rigorous numerics of blow-up solutions for ODEs with exponential nonlinearity

abstract: Our concerns here are blow-up solutions for ODEs with exponential nonlinearity from the viewpoint of dynamical systems and their numerical validations. As an example, the finite difference discretization of $u_t = u_{xx} + e^{u^m}$ with the homogeneous Dirichlet boundary condition is considered. Our idea is based on compactification of phase spaces and time-scale desingularization as in previous works. In the present case, treatment of exponential nonlinearity is the main issue. Fortunately, under a kind of exponential homogeneity of vector field, we can treat the problem in the same way as polynomial vector fields. In particular, we can characterize and validate blow-up solutions with their blow-up times for differential equations with such exponential nonlinearity in the similar way to previous works. A series of technical treatments of exponential nonlinearity in blow-up problems is also shown with concrete validation examples.


title: Global convergence of neuron birth-death dynamics

abstract: Neural networks with a large number of parameters admit a mean-field description, which has recently served as a theoretical explanation for the favorable training properties of "overparameterized" models. In this regime, gradient descent obeys a deterministic partial differential equation (PDE) that converges to a globally optimal solution for networks with a single hidden layer under appropriate assumptions. In this work, we propose a non-local mass transport dynamics that leads to a modified PDE with the same minimizer. We implement this non-local dynamics as a stochastic neuronal birth-death process and we prove that it accelerates the rate of convergence in the mean-field limit. We subsequently realize this PDE with two classes of numerical schemes that converge to the mean-field equation, each of which can easily be implemented for neural networks with finite numbers of parameters. We illustrate our algorithms with two models to provide intuition for the mechanism through which convergence is accelerated.


title: Hausdorff dimension of limit sets for projective Anosov representations

abstract: We study the relation between critical exponents and Hausdorff dimensions of limit sets for projective Anosov representations. We prove that the Hausdorff dimension of the symmetric limit set in $\mathbf{P}(\mathbb{R}^{n}) \times \mathbf{P}({\mathbb{R}^{n}}^*)$ is bounded between two critical exponents associated respectively to a highest weight and a simple root.


title: Probing topological properties of 3D lattice dimer model with neural networks

abstract: Machine learning methods are being actively considered as a new tool of describing many body physics. However, so far, their capabilities has been only demonstrated in previously studied models, such as e.g. Ising model. Here, we consider a simple problem, demonstrating that neural networks can be successfully used to give new insights in statistical physics. Specifically, we consider 3D lattice dimer model, which consists of sites forming a lattice and bonds connecting the neighboring sites, in such a way that every bond can be either empty or filled with a dimer, and the total number of dimers ending at one site is fixed to be one. Dimer configurations can be viewed as equivalent if they are connected through a series of local flips, i.e. simultaneous 'rotation' of a pair of parallel neighboring dimers. It turns out that the whole set of dimer configurations on a given 3D lattice can be split into distinct topological classes, such that dimer configurations belonging to different classes are not equivalent. In this paper we identify these classes by using neural networks. More specifically, we train the neural networks to distinguish dimer configurations from two known topological classes, and after it, we test them on dimer configurations from unknown topological classes. We demonstrate that 3D lattice dimer model on a bipartite lattice can be described by an integer topological invariant ('Hopf number'), whereas lattice dimer model on a non-bipartite lattice is described by $Z_2$ invariant. Thus, we demonstrate that neural networks can be successfully used to identify new topological phases in condensed matter systems, whose existence can be later verified by other (e.g. analytical) techniques.


title: Distribution-Dependent Analysis of Gibbs-ERM Principle

abstract: Gibbs-ERM learning is a natural idealized model of learning with stochastic optimization algorithms (such as Stochastic Gradient Langevin Dynamics and ---to some extent--- Stochastic Gradient Descent), while it also arises in other contexts, including PAC-Bayesian theory, and sampling mechanisms. In this work we study the excess risk suffered by a Gibbs-ERM learner that uses non-convex, regularized empirical risk with the goal to understand the interplay between the data-generating distribution and learning in large hypothesis spaces. Our main results are distribution-dependent upper bounds on several notions of excess risk. We show that, in all cases, the distribution-dependent excess risk is essentially controlled by the effective dimension $\mathrm{tr}\left(\boldsymbol{H}^{\star} (\boldsymbol{H}^{\star} + \lambda \boldsymbol{I})^{-1}\right)$ of the problem, where $\boldsymbol{H}^{\star}$ is the Hessian matrix of the risk at a local minimum. This is a well-established notion of effective dimension appearing in several previous works, including the analyses of SGD and ridge regression, but ours is the first work that brings this dimension to the analysis of learning using Gibbs densities. The distribution-dependent view we advocate here improves upon earlier results of Raginsky et al. (2017), and can yield much tighter bounds depending on the interplay between the data-generating distribution and the loss function. The first part of our analysis focuses on the localized excess risk in the vicinity of a fixed local minimizer. This result is then extended to bounds on the global excess risk, by characterizing probabilities of local minima (and their complement) under Gibbs densities, a results which might be of independent interest.


title: Band filling and disorder effects on the normal state thermoelectric behavior in MgB2

abstract: By a combined experimental and theoretical approach, we investigate normal state thermoelectric transport in MgB2, as a probe of selective disorder and doping in the sigma and pi bands. We calculate the temperature dependent diffusive Seebeck coefficient Sdiff(T) with the Boltzmann equation resolved in relaxation time approximation, taking into account the scattering with phonons and impurities, the effect of renormalization and the effect doping in a rigid band approximation. We show that selective disorder has a sizeable effect on the Sdiff magnitude, as it tunes the relative contributions of sigma and pi bands. Disorder also affects the Sdiff temperature dependences, eventually yielding a linear Sdiff(T) behavior in the dirty limit. We also show that band filling has opposite effects on S, depending on which band dominates transport. In parallel, we carry out Seebeck effect measurements on neutron-irradiated Mg11B2, and on two series of doped samples Mg1-xAlxB2 and Mg(B1-xCx)2. From comparison of calculated Sdiff(T) and experimental S(T) curves, we demonstrate that diffusive and phonon drag terms give comparable contributions in clean samples, but the phonon drag term is progressively suppressed with increasing disorder. In C and Al doped samples we observe very different experimental behaviors in terms of sign, magnitude and temperature dependence. Indeed, notwithstanding the similar electron doping introduced by both substitutions, C or Al doping yields disorder which mainly affects either sigma or pi bands, respectively. With the help of our ab-initio approach, we are able to disentangle the several effects and prove that Seebeck coefficient is a very sensitive probe of the kind of disorder.


title: Nonparametric Finite Time LTI System Identification

abstract: We address the problem of learning the parameters of a stable linear time invariant (LTI) system or linear dynamical system (LDS) with unknown latent space dimension, or order, from a single time--series of noisy input-output data. We focus on learning the best lower order approximation allowed by finite data. Motivated by subspace algorithms in systems theory, where the doubly infinite system Hankel matrix captures both order and good lower order approximations, we construct a Hankel-like matrix from noisy finite data using ordinary least squares. This circumvents the non-convexities that arise in system identification, and allows accurate estimation of the underlying LTI system. Our results rely on careful analysis of self-normalized martingale difference terms that helps bound identification error up to logarithmic factors of the lower bound. We provide a data-dependent scheme for order selection and find an accurate realization of system parameters, corresponding to that order, by an approach that is closely related to the Ho-Kalman subspace algorithm. We demonstrate that the proposed model order selection procedure is not overly conservative, i.e., for the given data length it is not possible to estimate higher order models or find higher order approximations with reasonable accuracy.


title: Competing frogs on $\mathbb{Z}^d$

abstract: A two-type version of the frog model on $\mathbb{Z}^d$ is formulated, where active type $i$ particles move according to lazy random walks with probability $p_i$ of jumping in each time step ($i=1,2$). Each site is independently assigned a random number of particles. At time 0, the particles at the origin are activated and assigned type 1 and the particles at one other site are activated and assigned type 2, while all other particles are sleeping. When an active type $i$ particle moves to a new site, any sleeping particles there are activated and assigned type $i$, with an arbitrary tie-breaker deciding the type if the site is hit by particles of both types in the same time step. We show that the event $G_i$ that type $i$ activates infinitely many particles has positive probability for all $p_1,p_2\in(0,1]$ ($i=1,2$). Furthermore, if $p_1=p_2$, then the types can coexist in the sense that $\mathbb{P}(G_1\cap G_2)>0$. We also formulate several open problems. For instance, we conjecture that, when the initial number of particles per site has a heavy tail, the types can coexist also when $p_1\neq p_2$.


title: Reduction and simultaneous doping of graphene oxide to repel LDL in treatment of atherosclerosis disease

abstract: Atherosclerotic disease develops when there is endothelial dysfunction, making low density lipoproteins (LDL) which are cholesterol carriers traveling in the blood, very likely to penetrate the endothelium layer triggering the formation of the atheromatous plaque and leading to the obstruction of blood vessels. For some time now, devices called stents have been used to re-open the light of the plugged artery, allowing the restoration of blood flow. However, this treatment has disadvantages like developing intra-stent restenosis and progression of the atherosclerosis, clogging the blood vessel again. For these reasons it is necessary to develop new materials and strategies focused on the prevention and correction of this pathology and its subsequent complications. In the present work, a reduction and in situ nitrogen doping method in graphene oxide (N-rGO) was studied. By means of XPS a 9.57% at of nitrogen doping was demonstrated and was found that the pyrrole type nitrogen is predominant. Due to nitrogen doping and pyrrolic sites, the N-rGO surface reached high levels of negative electric charge, which was verified by Electrostatic Force Microscopy (EFM) and wettability tests between the LDL at different oxidation levels and the synthesized N-rGO; contact angles between 95.7{\deg} and 130.4{\deg} were obtained which indicates repulsion between both materials. The N-rGO was also characterized by SEM and TEM where a wrinkled morphology typical of this type of material was observed and the hexagonal atomic arrangement of GO is distinguished. Finally, cell viability tests using 3-(4,5-dimethylthiazolyl)-2,5-diphenyltetrazolium bromide (MTT) assay were carried out on human umbilical cord endothelial cells (HUVEC) to analyze the cellular response at different concentrations of N-rGO, obtaining favorable behavior up to concentrations of 100 ({\mu}g/mL).


title: Chromium Oxide Formation on Nanosecond and Femtosecond Laser Irradiated Thin Chromium Films

abstract: Thin coatings of Chromium oxide have been used for applications as absorbing material in solar cells, as protections for magnetic data recording devices and as shields in flexible solar cells. Thin coatings of pure chromium were vacuum deposited on a glass substrate using hot electrons from tungsten filament. These coatings were then treated with a nanosecond and femtosecond laser in ambient conditions. The microstructure, morphology and the color of the coatings treated with laser sources were modified and there was a formation of an oxide layer due to the heat dissipation on the chromium coating from the energetic photons. High-resolution scanning electron microscope studies showed the morphological evolution that are directly correlated with the laser fluence of both the nanosecond and femtosecond lasers. This morphological evolution was accompanied by the microstructural change as observed from the x-ray diffraction patterns, the chromaticity response of the coating was studied by UV-Vis spectrometer and the response of the coating in the visible region evolved with the laser fluences. The Rutherford backscattering depth profiling of the laser treated coatings revealed the diffusion of oxygen atoms in the coating as a result of laser treatment fluence.


title: A new equation of state for dense hydrogen-helium mixtures

abstract: We present a new equation of state (EOS) for dense hydrogen/helium mixtures which covers a range of densities from $10^{-8}$ to $10^6$ g.cm$^{-3}$, pressures from $10^{-9}$ to $10^{13}$ GPa and temperatures from $10^{2}$ to $10^{8}$ K. The calculations combine the EOS of Saumon, Chabrier & vanHorn (1995) in the low density, low temperature molecular/atomic domain, the EOS of Chabrier & Potekhin (1998) in the high-density, high-temperature fully ionized domain, the limits of which differ for H and He, and ab initio quantum molecular dynamics (QMD) calculations in the intermediate density and temperature regime, characteristic of pressure dissociation and ionization. The EOS for the H/He mixture is based on the so-called additive volume law and thus does not take into account the interactions between the two species. A major improvement of the present calculations over existing ones is that we calculate the entropy over the entire density-temperature domain, a necessary quantity for stellar or planetary evolution calculations. The EOS results are compared with existing experimental data, namely Hugoniot shock experiments for pure H and He, and with first principle numerical simulations for both the single elements and the mixture. This new EOS covers a wide range of physical and astrophysical conditions, from jovian planets to solar-type stars, and recovers the existing relativistic EOS at very high densities, in the domains of white dwarfs and neutron stars.


title: Cu$_2$O Microcrystals Grown on Silicon as Platform for Quantum-Degenerate Excitons and Rydberg States

abstract: Cuprous oxide (Cu$_2$O) is a semiconductor with large exciton binding energy and significant technological importance in applications such as photovoltaics and solar water splitting. It is also a superior material system for quantum optics that enabled the observation of two intriguing phenomena, i.e. Rydberg excitons as solid-state analogue to highly-excited atomic states and dense exciton gases showing quantum degeneracy when approaching the phase transition to Bose-Einstein condensation. Previous experiments focused on natural bulk crystals due to major difficulties in growing high-quality synthetic samples. Here, we present Cu$_2$O microcrystals with excellent optical material quality capable of hosting both quantum-degenerate excitons and excited Rydberg states. Growth of Cu$_2$O with exceedingly low point defect levels was achieved on silicon by a scalable thermal oxidation process compatible with lithographic patterning. Using the latter, we demonstrate Rydberg excitons in site-controlled Cu$_2$O microstructures, paving the way for a plethora of applications in integrated quantum photonics.


title: Vertical, electrolyte-gated organic transistors: continuous operation in the MA/cm$^2$ regime and use as low-power artificial synapses

abstract: Organic semiconductors are usually not thought to show outstanding performance in highly-integrated, sub 100 nm transistors. Consequently, single-crystalline materials such as SWCNTs, MoS2 or inorganic semiconductors are the material of choice at these nanoscopic dimensions. Here, we show that using a novel vertical field-effect transistor design with a channel length of only 40 nm and a footprint of 2 x 80 x 80 nm$^2$, high electrical performance with organic polymers can be realized when using electrolyte gating. Our organic transistors combine high on-state current densities of above 3 MA/cm$^2$, on/off current modulation ratios of up to 108 and large transconductances of up to 5000 S/m. Given the high on-state currents at yet large on/off ratios, our novel structures also show promise for use in artificial neural networks, where they could operate as memristive devices with sub 100 fJ energy usage.


title: Detection of perfusion ROI as a quality control in perfusion analysis

abstract: In perfusion analysis automated approaches for image processing is preferable due to reduce time-consuming tasks for radiologists. Assessment of perfusion results quality is important step in development of algorithms for automated processing. One of them is an assessment of perfusion maps quality based on detection of perfusion ROI.


title: Asynchronous Delay-Aware Accelerated Proximal Coordinate Descent for Nonconvex Nonsmooth Problems

abstract: Nonconvex and nonsmooth problems have recently attracted considerable attention in machine learning. However, developing efficient methods for the nonconvex and nonsmooth optimization problems with certain performance guarantee remains a challenge. Proximal coordinate descent (PCD) has been widely used for solving optimization problems, but the knowledge of PCD methods in the nonconvex setting is very limited. On the other hand, the asynchronous proximal coordinate descent (APCD) recently have received much attention in order to solve large-scale problems. However, the accelerated variants of APCD algorithms are rarely studied. In this paper, we extend APCD method to the accelerated algorithm (AAPCD) for nonsmooth and nonconvex problems that satisfies the sufficient descent property, by comparing between the function values at proximal update and a linear extrapolated point using a delay-aware momentum value. To the best of our knowledge, we are the first to provide stochastic and deterministic accelerated extension of APCD algorithms for general nonconvex and nonsmooth problems ensuring that for both bounded delays and unbounded delays every limit point is a critical point. By leveraging Kurdyka-Lojasiewicz property, we will show linear and sublinear convergence rates for the deterministic AAPCD with bounded delays. Numerical results demonstrate the practical efficiency of our algorithm in speed.


title: A new line on the wide binary test of gravity

abstract: The relative velocity distribution of wide binary (WB) stars is sensitive to the law of gravity at the low accelerations typical of galactic outskirts. I consider the feasibility of this wide binary test using the `line velocity' method. This involves considering only the velocity components along the direction within the sky plane orthogonal to the systemic proper motion of each WB. I apply this technique to the WB sample of Hernandez et. al. (2019), carefully accounting for large-angle effects at one order beyond leading. Based on Monte Carlo trials, the uncertainty in the one-dimensional velocity dispersion is $\approx 100$ m/s when using sky-projected relative velocities. Using line velocities reduces this to $\approx 30$ m/s because these are much less affected by distance uncertainties. My analysis does not support the Hernandez et. al. (2019) claim of a clear departure from Newtonian dynamics beyond a radius of $\approx 10$ kAU, partly because I use $2\sigma$ outlier rejection to clean their sample first. Nonetheless, the uncertainties are small enough that existing WB data are nearly sufficient to distinguish Newtonian dynamics from Modified Newtonian Dynamics. I estimate that $\approx 1000$ WB systems will be required for this purpose if using only line velocities. In addition to a larger sample, it will also be important to control for systematics like undetected companions and moving groups. This could be done statistically. The contamination can be minimized by considering a narrow theoretically motivated range of parameters and focusing on how different theories predict different proportions of WBs in this region.


title: Solving the quantum dimer and six vertex models one electric field line at a time

abstract: The nature and the very existence of the resonant plaquette valence bond state that separates the classical columnar phase and the Rokhsar and Kivelson point in the quantum dimer model remains unsettled. Here we take a different line of attack on this model, and on the closely related six vertex model, by exploiting the global conservation law of the number of electric field lines. This allows us to study a single fluctuating electric field line which we show maps exactly onto a one dimensional spin chain. In the case of the six vertex model, the electric field line maps onto the celebrated spin 1/2 XXZ model which can be solved exactly. In the quantum dimer model, the electric field line is mapped onto a two-leg spin 1/2 ladder, which we study using numerical exact diagonalization. Our findings are consistent with the existence of three distinct phases including a Luttinger liquid phase, the one-dimensional precursor to the two-dimensional plaquette valence bond solid. The uncanny resemblance of our quasi-one-dimensional electric field line problem to the full two-dimensional problem suggests that much of the behavior of the latter might be understood by thinking of it as a closely packed array of field lines which themselves are undergoing non-trivial phase transitions.


title: Emergent hydrodynamics in non-equilibrium quantum systems

abstract: A tremendous amount of recent attention has focused on characterizing the dynamical properties of periodically driven many-body systems. Here, we use a novel numerical tool termed `density matrix truncation' (DMT) to investigate the late-time dynamics of large-scale Floquet systems. We find that DMT accurately captures two essential pieces of Floquet physics, namely, prethermalization and late-time heating to infinite temperature. Moreover, by implementing a spatially inhomogeneous drive, we demonstrate that an interplay between Floquet heating and diffusive transport is crucial to understanding the system's dynamics. Finally, we show that DMT also provides a powerful method for quantitatively capturing the emergence of hydrodynamics in static (un-driven) Hamiltonians; in particular, by simulating the dynamics of generic, large-scale quantum spin chains (up to L = 100), we are able to directly extract the energy diffusion coefficient.

categories: cond-mat.str-el cond-mat.mes-hall cond-mat.stat-mech quant-ph


title: The dynamical evolution of molecular clouds near the Galactic Centre - II. Spatial structure and kinematics of simulated clouds

abstract: The evolution of molecular clouds in galactic centres is thought to differ from that in galactic discs due to a significant influence of the external gravitational potential. We present a set of numerical simulations of molecular clouds orbiting on the 100-pc stream of the Central Molecular Zone (the central $\sim500$ pc of the Galaxy) and characterise their morphological and kinematic evolution in response to the background potential and eccentric orbital motion. We find that the clouds are shaped by strong shear and torques, by tidal and geometric deformation, and by their passage through the orbital pericentre. Within our simulations, these mechanisms control cloud sizes, aspect ratios, position angles, filamentary structure, column densities, velocity dispersions, line-of-sight velocity gradients, spin angular momenta, and kinematic complexity. By comparing these predictions to observations of clouds on the Galactic Centre 'dust ridge', we find that our simulations naturally reproduce a broad range of key observed morphological and kinematic features, which can be explained in terms of well-understood physical mechanisms. We argue that the accretion of gas clouds onto the central regions of galaxies, where the rotation curve turns over and the tidal field is fully compressive, is accompanied by transformative dynamical changes to the clouds, leading to collapse and star formation. This can generate an evolutionary progression of cloud collapse with a common starting point, which either marks the time of accretion onto the tidally-compressive region or of the most recent pericentre passage. Together, these processes may naturally produce the synchronised starbursts observed in numerous (extra)galactic nuclei.


title: GRB 190114C: from prompt to afterglow?

abstract: GRB 190114C is the first gamma-ray burst detected at Very High Energies (VHE, i.e. >300 GeV) by the MAGIC Cherenkov telescope. The analysis of the emission detected by the Fermi satellite at lower energies, in the 10 keV -- 100 GeV energy range, up to ~ 50 seconds (i.e. before the MAGIC detection) can hold valuable information. We analyze the spectral evolution of the emission of GRB 190114C as detected by the Fermi Gamma-Ray Burst Monitor (GBM) in the 10 keV -- 40 MeV energy range up to ~60 sec. The first 4 s of the burst feature a typical prompt emission spectrum, which can be fit by a smoothly broken power-law function with typical parameters. Starting on ~4 s post-trigger, we find an additional nonthermal component, which can be fit by a power law. This component rises and decays quickly. The 10 keV -- 40 MeV flux of the power-law component peaks at ~ 6 s; it reaches a value of 1.7e-5 erg cm-2 s-1. The time of the peak coincides with the emission peak detected by the Large Area Telescope (LAT) on board Fermi. The power-law spectral slope that we find in the GBM data is remarkably similar to that of the LAT spectrum, and the GBM+LAT spectral energy distribution seems to be consistent with a single component. This suggests that the LAT emission and the power-law component that we find in the GBM data belong to the same emission component, which we interpret as due to the afterglow of the burst. The onset time allows us to estimate the initial jet bulk Lorentz factor Gamma_0 is about 500, depending on the assumed circum-burst density.


title: Atomic Absorption Line Diagnostics for the Physical Properties of Red Supergiants

abstract: Red supergiants (RSGs) are evolved massive stars that represent extremes, in both their physical sizes and their cool temperatures, of the massive star population. Effective temperature ($T_{\rm eff}$) is the most critical physical property needed to place a RSG on the Hertzsprung-Russell Diagram, due to the stars' cool temperatures and resulting large bolometric corrections. Several recent papers have examined the potential utility of atomic line equivalent widths in cool supergiant spectra for determining $T_{\rm eff}$ and other physical properties (dorda2016b, dorda2016a) and found strong correlations between Ti I and Fe I spectral features and $T_{\rm eff}$ in earlier-type cool supergiants (G and early K) but poor correlations in M-type stars, a spectral subtype that makes up a significant fraction of RSGs. We have extended this work by measuring the equivalent widths of Ti, Fe, and Ca lines in late K- and M-type RSGs in the Milky Way, Large Magellanic Cloud, and Small Magellanic Cloud, and compared these results to the predictions of the MARCS stellar atmosphere models. Our analyses show a poor correlation between $T_{\rm eff}$ and the Fe I and Ti I lines in our observations (at odds with strong correlations predicted by stellar atmosphere models), but do find statistically significant correlations between $T_{\rm eff}$ and the Ca II triplet (CaT) features of Milky Way RSGs, suggesting that this could be a potential diagnostic tool for determining $T_{\rm eff}$ in M type supergiants. We also examine correlations between these spectral features and other physical properties of RSGs (including metallicity, surface gravity, and bolometric magnitude), and consider the underlying physics driving the evolution of atomic line spectra in RSGs.


title: Deep and narrow CO absorption revealing molecular clouds in the Hydra-A brightest cluster galaxy

abstract: Active galactic nuclei play a crucial role in the accretion and ejection of gas in galaxies. Although their outflows are well studied, finding direct evidence of accretion has proved very difficult and has so far been done for very few sources. A promising way to study the significance of cold accretion is by observing the absorption of an active galactic nucleus's extremely bright radio emission by the cold gas lying along the line-of-sight. As such, we present ALMA CO(1-0) and CO(2-1) observations of the Hydra-A brightest cluster galaxy (z=0.054) which reveal the existence of cold, molecular gas clouds along the line-of-sight to the galaxy's extremely bright and compact mm-continuum source. They have apparent motions relative to the central supermassive black hole of between -43 and -4 km s$^{-1}$ and are most likely moving along stable, low ellipticity orbits. The identified clouds form part of a $\sim$$10^{9}$ $\text{M}_{\odot}$, approximately edge-on disc of cold molecular gas. With peak CO(2-1) optical depths of $\tau$=0.88 $^{+0.06}_{-0.06}$, they include the narrowest and by far the deepest absorption of this type which has been observed to date in a brightest cluster galaxy. By comparing the relative strengths of the lines for the most strongly absorbing region, we are able to estimate a gas temperature of $42^{+25}_{-11}$ K and line-of-sight column densities of $N_{CO}=2^{+3}_{-1}\times 10 ^{17} cm^{-2}$ and $N_{ H_{2} }=7^{+10}_{-4}\times 10 ^{20} cm^{-2}$.


title: Gravitational-wave sources from mergers of binary black-holes catalyzed by fly-bys interactions in the field

abstract: Several scenarios were suggested for the origins of gravitational-wave (GW) sources from mergers of stellar binary black holes (BBHs). Here we propose a novel origin through catalyzed formation of GW-sources from ultra-wide binaries in the field. Such binaries experience perturbations from random stellar fly-bys which excite their eccentricities. Once a wide-binary is driven to a sufficiently small peri-center approach, GW-emission becomes significant, and the binary inspirals and merges. We derive an analytic model and verify it with numerical calculation to compute the merger rate to be $\sim10{\rm \times f_{wide}\,{\rm Gpc}^{-3}yr^{-1}}$ ($f_{{\rm wide}}$ is the fraction of wide BH-binaries), which is comparable to the observationally inferred rate. The observational signatures from this channel include spin-orbit misalignment; preference for high mass-ratio BBH; preference for high velocity-dispersion host-galaxies; and a uniform delay-time distribution.


title: Strong Cosmic Censorship in higher-dimensional Reissner-Nordstr\"{o}m-de Sitter spacetime

abstract: It was recently shown that Strong Cosmic Censorship might be violated for near-extremally-charged black holes in 4-dimensional de Sitter space under scalar perturbations. Here, we extend the study of neutral massless scalar perturbations in higher dimensions and discuss the dimensional influence on the validity of Strong Cosmic Censorship hypothesis. By giving an elaborate description of neutral massless scalar perturbations of Reissner-Nordstr\"{o}m-de Sitter black holes in $d=4,5$ and $6$ dimensions we conclude that Strong Cosmic Censorship is violated near extremality.


title: Fast radio bursts as synchrotron maser emission from decelerating relativistic blast waves

abstract: Fast radio bursts (FRB) can arise from synchrotron maser emission at ultra-relativistic magnetized shocks, such as produced by flare ejecta from young magnetars. We combine PIC simulation results for the synchrotron maser with the dynamics of self-similar shock deceleration, as commonly applied to GRBs, to explore the implications for FRB emission. We assume the upstream environment into which the ultra-relativistic ejecta collides is a mildly relativistic baryon-loaded shell released following a previous flare, motivated by the high electron-ion injection rate Mdot ~ 1e19-1e21 g/s needed on larger scales to power the persistent radio nebula coincident with the repeating burster FRB 121102 and its high inferred rotation measure. The observed radio fluence peaks once the optical depth ahead of the shock to induced Compton scattering decreases to <~ few, a condition which places a GHz observer on the high frequency tail of the maser SED. Given intervals between ion shell ejection events ~1e5 s similar to the occurrence rate of the most powerful bursts from FRB 121102, we demonstrate the production of FRBs of frequency ~ 0.1-10 GHz, isotropic radiated energies ~1e37-1e40 erg and durations ~0.1-10 ms for flares of energy ~1e43-1e45 erg. Deceleration of the blast wave, and increasing transparency of the upstream medium, generates a temporal decay of the peak frequency, similar to the downward drift seen in the sub-bursts of FRB 121102 and FRB 180814.J0422+73. The delay >~ 1e5 s between major ion-injection events needed to clear sufficiently low densities around the engine for FRB emission could explain prolonged "dark" periods and clustered burst arrival times, and lead to stochastic variation in the dispersion measure. Thermal electrons heated at the shock generate a short-lived <~ 1 ms (1 s) synchrotron transient at gamma-ray (X-ray) energies, analogous to a GRB afterglow.


title: Non-Riemannian geometry of M-theory

abstract: We construct a background for M-theory that is moduli free. This background is then shown to be related to a topological phase of the $\mathrm{E}_{8(8)}$ exceptional field theory (ExFT). The key ingredient in the construction is the embedding of non-Riemannian geometry in ExFT. This allows one to describe non-relativistic geometries, such as Newton-Cartan or Gomis-Ooguri-type limits, using the ExFT framework originally developed to describe maximal supergravity. This generalises previous work by Morand and Park in the context of double field theory.


title: Breaking cosmic degeneracies: Disentangling neutrinos and modified gravity with kinematic information

abstract: Searches for modified gravity in the large-scale structure try to detect the enhanced amplitude of density fluctuations caused by the fifth force present in many of these theories. Neutrinos, on the other hand, suppress structure growth below their free-streaming length. Both effects take place on comparable scales, and uncertainty in the neutrino mass leads to a degeneracy with modified gravity parameters for probes that are measuring the amplitude of the matter power spectrum. We explore the possibility to break the degeneracy between modified gravity and neutrino effects in the growth of structures by considering kinematic information related to either the growth rate on large scales or the virial velocities inside of collapsed structures. In order to study the degeneracy up to fully non-linear scales, we employ a suite of $N$-body simulations including both $f(R)$ modified gravity and massive neutrinos. Our results indicate that velocity information provides an excellent tool to distinguish massive neutrinos from modified gravity. Models with different values of neutrino masses and modified gravity parameters possessing a comparable matter power spectrum at a given time have different growth rates. This leaves imprints in the velocity divergence, which is therefore better suited than the amplitude of density fluctuations to tell the models apart. In such models with a power spectrum comparable to $\Lambda$CDM today, the growth rate is strictly enhanced. We also find the velocity dispersion of virialised clusters to be well suited to constrain deviations from general relativity without being affected by the uncertainty in the sum of neutrino masses.


title: Exploring the formation by core accretion and the luminosity evolution of directly imaged planets: The case of HIP 65426 b

abstract: A low-mass companion to the two-solar mass star HIP65426 has recently been detected by SPHERE at around 100 au from its host. Explaining the presence of super-Jovian planets at large separations, as revealed by direct imaging, is currently an open question. We want to derive statistical constraints on the mass and initial entropy of HIP65426b and to explore possible formation pathways of directly imaged objects within the core-accretion paradigm, focusing on HIP65426b. Constraints on the planet's mass and post-formation entropy are derived from its age and luminosity combined with cooling models. For the first time, the results of population synthesis are also used to inform the results. Then, a formation model that includes N-body dynamics with several embryos per disc is used to study possible formation histories and the properties of possible additional companions. Finally, the outcomes of two- and three-planet scattering in the post-disc phase are analysed, taking tides into account. The mass of HIP65426b is found to be Mp = 9.9 +1.1 -1.8 MJ using the hot population and Mp = 10.9 +1.4 -2.0 MJ with the cold-nominal population. Core formation at small separations from the star followed by outward scattering and runaway accretion at a few hundred AU succeeds in reproducing the mass and separation of HIP65426b. Alternatively, systems having two or more giant planets close enough to be on an unstable orbit at disc dispersal are likely to end up with one planet on a wide HIP65426b-like orbit with a relatively high eccentricity (>~ 0.5). If this scattering scenario explains its formation, HIP65426b is predicted to have a high eccentricity and to be accompanied by one or several roughly Jovian-mass planets at smaller semi-major axes, which also could have a high eccentricity. This could be tested by further direct-imaging as well as radial-velocity observations.


title: Stabilizing Inputs to Approximated Nonlinear Functions for Inference with Homomorphic Encryption in Deep Neural Networks

abstract: Leveled Homomorphic Encryption (LHE) offers a potential solution that could allow sectors with sensitive data to utilize the cloud and securely deploy their models for remote inference with Deep Neural Networks (DNN). However, this application faces several obstacles due to the limitations of LHE. One of the main problems is the incompatibility of commonly used nonlinear functions in DNN with the operations supported by LHE, i.e. addition and multiplication. As common in LHE approaches, we train a model with a nonlinear function, and replace it with a low-degree polynomial approximation at inference time on private data. While this typically leads to approximation errors and loss in prediction accuracy, we propose a method that reduces this loss to small values or eliminates it entirely, depending on simple hyper-parameters. This is achieved by the introduction of a novel and elegantly simple Min-Max normalization scheme, which scales inputs to nonlinear functions into ranges with low approximation error. While being intuitive in its concept and trivial to implement, we empirically show that it offers a stable and effective approximation solution to nonlinear functions in DNN. In return, this can enable deeper networks with LHE, and facilitate the development of security- and privacy-aware analytics applications.


title: Dependence on parameters of CW globalizations of families of Harish-Chandra modules and the meromorphic continuation of $C^{\infty}$ Eisenstein series

abstract: The first main result is that the Casselman-Wallach Globalization of a real analytic family of Harish-Chandra modules is continuous in the parameter. Our proof of this result uses results from the thesis of Vincent van der Noort in several critical ways. In his thesis the holomorphic version of the result was proved in the case when the parameter space is a one dimensional complex manifold up to a branched covering. The second main result is a proof of the meromorphic continuation of $C^{\infty}$ Eisenstein series.using Langlands' results in the $K$ finite case as an application of the methods in the proof of the first part.


title: The Zwicky Transient Facility: Data Processing, Products, and Archive

abstract: The Zwicky Transient Facility (ZTF) is a new robotic time-domain survey currently in progress using the Palomar 48-inch Schmidt Telescope. ZTF uses a 47 square degree field with a 600 megapixel camera to scan the entire northern visible sky at rates of ~3760 square degrees/hour to median depths of g ~ 20.8 and r ~ 20.6 mag (AB, 5sigma in 30 sec). We describe the Science Data System that is housed at IPAC, Caltech. This comprises the data-processing pipelines, alert production system, data archive, and user interfaces for accessing and analyzing the products. The realtime pipeline employs a novel image-differencing algorithm, optimized for the detection of point source transient events. These events are vetted for reliability using a machine-learned classifier and combined with contextual information to generate data-rich alert packets. The packets become available for distribution typically within 13 minutes (95th percentile) of observation. Detected events are also linked to generate candidate moving-object tracks using a novel algorithm. Objects that move fast enough to streak in the individual exposures are also extracted and vetted. The reconstructed astrometric accuracy per science image with respect to Gaia is typically 45 to 85 milliarcsec. This is the RMS per axis on the sky for sources extracted with photometric S/N >= 10. The derived photometric precision (repeatability) at bright unsaturated fluxes varies between 8 and 25 millimag. Photometric calibration accuracy with respect to Pan-STARRS1 is generally better than 2%. The products support a broad range of scientific applications: fast and young supernovae, rare flux transients, variable stars, eclipsing binaries, variability from active galactic nuclei, counterparts to gravitational wave sources, a more complete census of Type Ia supernovae, and Solar System objects.


title: Dynamic hierarchies in temporal directed networks

abstract: The outcome of interactions in many real-world systems can be often explained by a hierarchy between the participants. Discovering hierarchy from a given directed network can be formulated as follows: partition vertices into levels such that, ideally, there are only forward edges, that is, edges from upper levels to lower levels. In practice, the ideal case is impossible, so instead we minimize some penalty function on the backward edges. One practical option for such a penalty is agony, where the penalty depends on the severity of the violation. In this paper we extend the definition of agony to temporal networks. In this setup we are given a directed network with time stamped edges, and we allow the rank assignment to vary over time. We propose 2 strategies for controlling the variation of individual ranks. In our first variant, we penalize the fluctuation of the rankings over time by adding a penalty directly to the optimization function. In our second variant we allow the rank change at most once. We show that the first variant can be solved exactly in polynomial time while the second variant is NP-hard, and in fact inapproximable. However, we develop an iterative method, where we first fix the change point and optimize the ranks, and then fix the ranks and optimize the change points, and reiterate until convergence. We show empirically that the algorithms are reasonably fast in practice, and that the obtained rankings are sensible.