title: Interference-Aided Energy Harvesting: Cognitive Relaying with Multiple Primary Transceivers

abstract: We consider a spectrum sharing scenario where a secondary transmitter (ST) communicates with its destination via a decode-and-forward secondary relay (SR) in the presence of interference from multiple primary transmitters. The SR harvests energy from received radio-frequency signals that include primary interference and uses it to forward the information to the secondary destination. The relay adopts a time switching policy that switches between energy harvesting and information decoding over the time. Under the primary outage constraints and the peak power constraints at both ST and SR, to determine the average secondary throughput, we derive exact analytical expressions for the secondary outage probability and the ergodic capacity, which characterize the delay-limited and the delay-tolerant transmissions, respectively. We also investigate the effects of the number of primary transceivers and the peak power constraints on the optimal energy harvesting time that maximizes the secondary throughput. By utilizing the primary interference as an energy source, the secondary network achieves a better throughput performance compared to the case where the primary interference is ignored for energy harvesting purpose. Finally, we consider a case where ST also harvests energy from primary transmissions and compare its throughput performance with that of the non-energy harvesting ST case.


title: A hybrid Alternating Least Squares -- TT Cross algorithm for parametric PDEs

abstract: We consider the approximate solution of parametric PDEs using the low-rank Tensor Train (TT) decomposition. Such parametric PDEs arise for example in uncertainty quantification problems in engineering applications. We propose an algorithm that is a hybrid of the alternating least squares and the TT cross methods. It computes a TT approximation of the whole solution, which is beneficial when multiple quantities of interest are sought. This might be needed, for example, for the computation of the probability density function (PDF) via the maximum entropy method [Kavehrad and Joseph, IEEE Trans. Comm., 1986]. The new algorithm exploits and preserves the block diagonal structure of the discretized operator in stochastic collocation schemes. This disentangles computations of the spatial and parametric degrees of freedom in the TT representation. In particular, it only requires solving independent PDEs at a few parameter values, thus allowing the use of existing high performance PDE solvers. In our numerical experiments, we apply the new algorithm to the stochastic diffusion equation and compare it with preconditioned steepest descent in the TT format, as well as with (multilevel) quasi-Monte Carlo and dimension-adaptive sparse grids methods. For sufficiently smooth random fields the new approach is orders of magnitude faster.


title: Assessing the foundation of the Trojan Horse Method

abstract: We discuss the foundation of the Trojan Horse Method (THM) within the Inclusive Non-Elastic Breakup (INEB) theory. We demonstrate that the direct part of the INEB cross section, which is of two-step character, becomes, in the DWBA limit of the three-body theory with appropriate approximations and redefinitions, similar in structure to the one-step THM cross section. We also discuss the connection of the THM to the Surrogate Method (SM), which is a genuine two-step process.


title: Hamiltonian and Algebraic Theories of Gapped Boundaries in Topological Phases of Matter

abstract: We present an exactly solvable lattice Hamiltonian to realize gapped boundaries of Kitaev's quantum double models for Dijkgraaf-Witten theories. We classify the elementary excitations on the boundary, and systematically describe the bulk-to-boundary condensation procedure. We also present the parallel algebraic/categorical structure of gapped boundaries.

categories: cond-mat.str-el math-ph math.MP math.QA quant-ph


title: Widely tunable on-chip microwave circulator for superconducting quantum circuits

abstract: We report on the design and performance of an on-chip microwave circulator with a widely (GHz) tunable operation frequency. Non-reciprocity is created with a combination of frequency conversion and delay, and requires neither permanent magnets nor microwave bias tones, allowing on-chip integration with other superconducting circuits without the need for high-bandwidth control lines. Isolation in the device exceeds 20 dB over a bandwidth of tens of MHz, and its insertion loss is small, reaching as low as 0.9 dB at select operation frequencies. Furthermore, the device is linear with respect to input power for signal powers up to hundreds of fW ($\approx 10^3$ circulating photons), and the direction of circulation can be dynamically reconfigured. We demonstrate its operation at a selection of frequencies between 4 and 6 GHz.


title: Pushing the Limits of Online Auto-tuning: Machine Code Optimization in Short-Running Kernels

abstract: We propose an online auto-tuning approach for computing kernels. Differently from existing online auto-tuners, which regenerate code with long compilation chains from the source to the binary code, our approach consists on deploying auto-tuning directly at the level of machine code generation. This allows auto-tuning to pay off in very short-running applications. As a proof of concept, our approach is demonstrated in two benchmarks, which execute during hundreds of milliseconds to a few seconds only. In a CPU-bound kernel, the average speedups achieved are 1.10 to 1.58 depending on the target micro-architecture, up to 2.53 in the most favourable conditions (all run-time overheads included). In a memory-bound kernel, less favourable to our runtime auto-tuning optimizations, the average speedups are 1.04 to 1.10, up to 1.30 in the best configuration. Despite the short execution times of our benchmarks, the overhead of our runtime auto-tuning is between 0.2 and 4.2% only of the total application execution times. By simulating the CPU-bound application in 11 different CPUs, we showed that, despite the clear hardware disadvantage of In-Order (io) cores vs. Out-of-Order (ooo) equivalent cores, online auto-tuning in io CPUs obtained an average speedup of 1.03 and an energy efficiency improvement of 39~\% over the SIMD reference in ooo CPUs.


title: Factoring the Cycle Aging Cost of Batteries Participating in Electricity Markets

abstract: When participating in electricity markets, owners of battery energy storage systems must bid in such a way that their revenues will at least cover their true cost of operation. Since cycle aging of battery cells represents a substantial part of this operating cost, the cost of battery degradation must be factored in these bids. However, existing models of battery degradation either do not fit market clearing software or do not reflect the actual battery aging mechanism. In this paper we model battery cycle aging using a piecewise linear cost function, an approach that provides a close approximation of the cycle aging mechanism of electrochemical batteries and can be incorporated easily into existing market dispatch programs. By defining the marginal aging cost of each battery cycle, we can assess the actual operating profitability of batteries. A case study demonstrates the effectiveness of the proposed model in maximizing the operating profit of a battery energy storage system taking part in the ISO New England energy and reserve markets.


title: Systematic study of {\alpha} decay using different versions of proximity formalism

abstract: Finding the best model to describe the {\alpha}-decay process is an old and ongoing challenge in nuclear physics. The present work systematically studied {\alpha}-decay half-lives for the favored ground-state-to-ground-state transitions of 344 isotopes of nuclei with 52 <= Z <= 107 using 28 versions of the proximity potential model in the framework of the WKB approximation. The present study introduces the best proximity versions with the fewest deviations with respect to experimental values. The models for Prox. 77-set 4, Prox. 77-set 5, and Dutt 2011 with the root-mean-square deviations (RMSDs) of <1 were found to predict {\alpha}-decay half-lives better than the other models. Comparison with fusion studies shows that Dutt 2011 is an appropriate model both for {\alpha}-decay studies and for prediction of the barrier characteristic in heavy-ion fusion reactions. The calculation of {\alpha}-decay half-lives were repeated for even-even, even-odd, odd-even, and odd-odd nuclei. This detailed comparative study reveals that for these versions the half-lives of the even-even nuclei with RMSDs of <0.6 show less deviation than the even-odd, odd-even, and odd-odd nuclei.


title: Stability in Chaos

abstract: Intrinsic instability of trajectories characterizes chaotic dynamical systems. We report here that trajectories can exhibit a surprisingly high degree of stability, over a very long time, in a chaotic dynamical system. We provide a detailed quantitative description of this effect for a one-dimensional model of inertial particles in a turbulent flow using large-deviation theory. Specifically, the determination of the entropy function for the distribution of finite-time Lyapunov exponents reduces to the analysis of a Schr\"odinger equation, which is tackled by semi-classical methods.


title: Quantum oscillations and a non-trivial Berry phase in the noncentrosymmetric superconductor BiPd

abstract: We report the measurements of de Haas-van Alphen (dHvA) oscillations in the noncentrosymmetric superconductor BiPd. Several pieces of a complex multi-sheet Fermi surface are identified, including a small pocket (frequency 40 T) which is three dimensional and anisotropic. From the temperature dependence of the amplitude of the oscillations, the cyclotron effective mass is ($0.18$ $\pm$ 0.1) $m_e$. Further analysis showed a non-trivial $\pi$-Berry phase is associated with the 40 T pocket, which strongly supports the presence of topological states in bulk BiPd and may result in topological superconductivity due to the proximity coupling to other bands.


title: Role of atoms in atomic gravitational-wave detectors

abstract: Recently, it has been proposed that space-based atomic sensors may be used to detect gravitational waves. These proposals describe the sensors either as clocks or as atom interferometers. Here, we seek to explore the fundamental similarities and differences between the two types of proposals. We present a framework in which the fundamental mechanism for sensitivity is identical for clock and atom interferometer proposals, with the key difference being whether or not the atoms are tightly confined by an external potential. With this interpretation in mind, we propose two major enhancements to detectors using confined atoms, which allow for an enhanced sensitivity analogous to large-momentum-transfer (LMT) used in atom interferometry (though with no transfer of momentum to the atoms), and a way to extend the useful coherence time of the sensor beyond the atom's excited state lifetime.


title: Temporal Network Comparison using Graphlet-orbit Transitions

abstract: Networks are widely used to model real-world systems and uncover their topological features. Network properties such as the degree distribution and shortest path length have been computed in numerous real-world networks, and most of them have been shown to be both scale-free and small-world networks. Graphlets and network motifs are subgraph patterns that capture richer structural information than aforementioned global network properties, and these local features are often used for network comparison. However, past work on graphlets and network motifs is almost exclusively applicable only for static networks. Many systems are better represented as temporal networks which depict not only how a system was at a given stage but also how they evolved. Time-dependent information is crucial in temporal networks and, by disregarding that data, static methods can not achieve the best possible results. This paper introduces an extension of graphlets for temporal networks. Our proposed method enumerates all 4-node graphlet-orbits in each network-snapshot, building the corresponding orbit-transition matrix in the process. Our hypothesis is that networks representing similar systems have characteristic orbit transitions which better identify them than simple static patterns, and this is assessed on a set of real temporal networks split into categories. In order to perform temporal network comparison we put forward an orbit-transition-agreement metric (OTA). OTA correctly groups a set of temporal networks that both static network motifs and graphlets fail to do so adequately. Furthermore, our method produces interpretable results which we use to uncover characteristic orbit transitions, and that can be regarded as a network-fingerprint.


title: Dark matter contribution to $b\to s \mu^+ \mu^-$ anomaly in local $U(1)_{L_\mu-L_\tau}$ model

abstract: We propose a local $U(1)_{L_\mu-L_\tau}$ model to explain $b \to s \mu^+ \mu^-$ anomaly observed at the LHCb and Belle experiments. The model also has a natural dark matter candidate $N$. We introduce $SU(2)_L$-doublet colored scalar $\widetilde{q}$ to mediate $b \to s$ transition at one-loop level. The $U(1)_{L_\mu-L_\tau}$ gauge symmetry is broken spontaneously by the scalar $S$. All the new particles are charged under $U(1)_{L_\mu-L_\tau}$. We can obtain $C_9^{\mu,{\rm NP}} \sim -1$ to solve the $b \to s\mu^+\mu^-$ anomaly and can explain the correct dark matter relic density of the universe, $\Omega_{\rm DM} h^2 \approx 0.12$, simultaneously, while evading constraints from electroweak precision tests, neutrino trident experiments and other quark flavor-changing loop processes such as $b \to s \gamma$ and $B_s-\overline{B}_s$ mixing. Our model can be tested by searching for $Z'$ and new colored scalar at the LHC and $B \to K^* \nu \overline{\nu}$ process at Belle-II.


title: Asymptotic representations of augmented q-Onsager algebra and boundary K-operators related to Baxter Q-operators

abstract: We consider intertwining relations of the augmented $q$-Onsager algebra introduced by Ito and Terwilliger, and obtain generic (diagonal) boundary $K$-operators in terms of the Cartan element of $U_{q}(sl_2)$. These $K$-operators solve reflection equations. Taking appropriate limits of these $K$-operators in Verma modules, we derive $K$-operators for Baxter Q-operators and corresponding reflection equations.


title: Direct observation of exchange-driven spin interactions in one-dimensional system

abstract: We present experimental results of transverse electron focusing measurements performed on an n-type GaAs based mesoscopic device consisting of one-dimensional (1D) quantum wires as injector and detector. We show that non-adiabatic injection of 1D electrons at a conductance of e$^2$/h results in a single first focusing peak, which on gradually increasing the injector conductance up to 2e$^2$/h , produces asymmetric two sub-peaks in the first focusing peak, each sub-peak representing the population of spin-state arising from the spatially separated spins in the injector. Further increasing the conductance flips the spin-states in the 1D channel thus reversing the asymmetry in the sub-peaks. On applying a source-drain bias, the spin-gap, so obtained, can be resolved thus providing evidence of exchange interaction induced spin polarisation in the 1D systems.


title: A Class of topological nodal rings and its realization in carbon networks

abstract: Topological nodal rings can be classified into three types according to the slopes in their energy dispersion. The first two are made of type-I and II nodal points, respectively, while the third is made of both. In carbon networks, all three types can exist. Under strain, phase transitions from a topological metal to a semiconductor, take place, and at the transition points, these nodal rings shrink into type-I, II, and III semi-Dirac points. These topological features exhibit diverse electron-hole pocket patterns and Landau levels, which give rise to exotic transport properties.


title: Casimir-Lifshitz Force for Nonreciprocal Media and Applications to Photonic Topological Insulators

abstract: Based on the theory of macroscopic quantum electrodynamics, we generalize the expression of the Casimir force for nonreciprocal media. The essential ingredient of this result is the Green's tensor between two nonreciprocal semi-infinite slabs including a reflexion matrix with four coefficients that mixes optical polarizations. This Green's tensor does not obey Lorentz's reciprocity and thus violates time-reversal symmetry. The general result for the Casimir force is analyzed in the retarded and nonretarded limits, concentrating on the influences arising from reflections with or without change of polarization. In a second step we apply our general result to a photonic topological insulator whose nonreciprocity stems from an anisotropic permittivity tensor, namely InSb. We show that there is a regime for the distance between the slabs where the magnitude of the Casimir force is tunable by an external magnetic field. Furthermore the strength of this tuning depends on the orientation of the magnetic field with respect to the slab surfaces.


title: An Efficient Approach to Communication-aware Path Planning for Long-range Surveillance Missions undertaken by UAVs

abstract: While using drones for remote surveillance missions, it is mandatory to do path planning of the vehicle since these are pilot-less vehicles. Path planning, whether offline or online, entails setting up the path as a sequence of locations in the 3D Euclidean space, whose coordinates happen to be latitude, longitude and altitude. For the specific application of remote surveillance of long linear infrastructures in non-urban terrain, the continuous 3D-ESP problem practically entails two important scalar costs. The first scalar cost is the distance traveled along the planned path. Since drones are battery operated, hence it is needed that the path length between fixed start and goal locations of a mission should be minimal at all costs. The other scalar cost is the cost of transmitting the acquired video during the mission of remote surveillance, via a camera mounted in the drone's belly. Because of the length of surveillance target which is long linear infrastructure, the amount of video generated is very high and cannot be generally stored in its entirety, on board. If the connectivity is poor along certain segments of a naive path, to boost video transmission rate, the transmission power of the signal is kept high, which in turn dissipates more battery energy. Hence a path is desired that simultaneously also betters what is known as communication cost. These two costs trade-off, and hence Pareto optimization is needed for this 3D bi-objective Euclidean shortest path problem. In this report, we study the mono-objective offline path planning problem, based on the distance cost, while posing the communication cost as an upper-bounded constraint. The bi-objective path planning solution is sketched out towards the end.


title: Impact of Cosmic Ray Transport on Galactic Winds

abstract: The role of cosmic rays generated by supernovae and young stars has very recently begun to receive significant attention in studies of galaxy formation and evolution due to the realization that cosmic rays can efficiently accelerate galactic winds. Microscopic cosmic ray transport processes are fundamental for determining the efficiency of cosmic ray wind driving. Previous studies focused on modeling of cosmic ray transport either via constant diffusion coefficient or via streaming proportional to the Alfven speed. However, in predominantly cold, neutral gas, cosmic rays can propagate faster than in the ionized medium and the effective transport can be substantially larger; i.e., cosmic rays can decouple from the gas. We perform three-dimensional magnetohydrodynamical simulations of patches of galactic disks including the effects of cosmic rays. Our simulations include the decoupling of cosmic rays in the cold, neutral interstellar medium. We find that, compared to the ordinary diffusive cosmic ray transport case, accounting for the decoupling leads to significantly different wind properties such as the gas density and temperature, significantly broader spatial distribution of cosmic rays, and larger wind speed. These results have implications for X-ray, gamma-ray and radio emission, and for the magnetization and pollution of the circumgalactic medium by cosmic rays.


title: Production of vector resonances at the LHC via WZ-scattering: a unitarized EChL analysis

abstract: In the present work we study the production of vector resonances at the LHC by means of the vector boson scattering $WZ \to WZ$ and explore the sensitivities to these resonances for the expected future LHC luminosities. We are assuming that these vector resonances are generated dynamically from the self interactions of the longitudinal gauge bosons, $W_L$ and $Z_L$, and work under the framework of the electroweak chiral Lagrangian to describe in a model independent way the supposedly strong dynamics of these modes. The properties of the vector resonances, mass, width and couplings to the $W$ and $Z$ gauge bosons are derived from the inverse amplitude method approach. We implement all these features into a single model, the IAM-MC, adapted for MonteCarlo, built in a Lagrangian language in terms of the electroweak chiral Lagrangian and a chiral Lagrangian for the vector resonances, which mimics the resonant behavior of the IAM and provides unitary amplitudes. The model has been implemented in MadGraph, allowing us to perform a realistic study of the signal versus background events at the LHC. In particular, we have focused our study on the $pp\to WZjj$ type of events, discussing first on the potential of the hadronic and semileptonic channels of the final $WZ$, and next exploring in more detail the clearest signals. These are provided by the leptonic decays of the gauge bosons, leading to a final state with $\ell_1^+\ell_1^-\ell_2^+\nu jj$, $\ell=e,\mu$, having a very distinctive signature, and showing clearly the emergence of the resonances with masses in the range of 1.5-2.5 TeV, which we have explored.


title: Toric h-vectors and Chow Betti Numbers of Dual Hypersimplices

abstract: The toric h-numbers of a dual hypersimplex and the Chow Betti numbers of the normal fan of a hypersimplex are the ranks of intersection cohomology and Chow cohomology respectively of the torus orbit closure of a generic point in the Grassmannian. We give explicit formulas for these numbers. We also show that similar formulas hold for the coordinator numbers of type A^* lattices.


title: Capturing the diversity of biological tuning curves using generative adversarial networks

abstract: Tuning curves characterizing the response selectivities of biological neurons often exhibit large degrees of irregularity and diversity across neurons. Theoretical network models that feature heterogeneous cell populations or random connectivity also give rise to diverse tuning curves. However, a general framework for fitting such models to experimentally measured tuning curves is lacking. We address this problem by proposing to view mechanistic network models as generative models whose parameters can be optimized to fit the distribution of experimentally measured tuning curves. A major obstacle for fitting such models is that their likelihood function is not explicitly available or is highly intractable to compute. Recent advances in machine learning provide ways for fitting generative models without the need to evaluate the likelihood and its gradient. Generative Adversarial Networks (GAN) provide one such framework which has been successful in traditional machine learning tasks. We apply this approach in two separate experiments, showing how GANs can be used to fit commonly used mechanistic models in theoretical neuroscience to datasets of measured tuning curves. This fitting procedure avoids the computationally expensive step of inferring latent variables, e.g. the biophysical parameters of individual cells or the particular realization of the full synaptic connectivity matrix, and directly learns model parameters which characterize the statistics of connectivity or of single-cell properties. Another strength of this approach is that it fits the entire, joint distribution of experimental tuning curves, instead of matching a few summary statistics picked a priori by the user. More generally, this framework opens the door to fitting theoretically motivated dynamical network models directly to simultaneously or non-simultaneously recorded neural responses.


title: Hidden Conformal Symmetry of Smooth Braneworld Scenarios

abstract: In this paper we generalize our previous model (arXiv: 1705.09331), on a hidden conformal symmetry of smooth braneworld scenarios, to the case with two real scalar fields non-minimally coupled to gravity. The gauge condition reduces the action of the system to the action were gravity minimally couples to one of the scalar fields, plus a cosmological constant. We show that, depending on the internal symmetry of the scalar fields, the two possibilities, $SO(2)$ or $SO(1,1)$, emerge. In the $SO(2)$ case we get a ghost-like scalar field action, which can describe two models -- Standing Wave and Sine-Gordon smooth braneworlds. For the $SO(1,1)$ case we get the standard sign for the kinetic part of the scalar field. By breaking the $SO(1,1)$ symmetry (but keeping the conformal one) we are able to get two Randall-Sundrum models, with a non-minimal coupling and with a scalar field having hyperbolic potential. We conclude that this method can be seen as a solution-generating technique and a natural way to introduce non-trivial scalar fields that can provide smooth braneworld models.


title: Fast Restricted Causal Inference

abstract: Hidden variables are well known sources of disturbance when recovering belief networks from data based only on measurable variables. Hence models assuming existence of hidden variables are under development. This paper presents a new algorithm "accelerating" the known CI algorithm of Spirtes, Glymour and Scheines {Spirtes:93}. We prove that this algorithm does not produces (conditional) independencies not present in the data if statistical independence test is reliable. This result is to be considered as non-trivial since e.g. the same claim fails to be true for FCI algorithm, another "accelerator" of CI, developed in {Spirtes:93}.


title: The Reversible Residual Network: Backpropagation Without Storing Activations

abstract: Deep residual networks (ResNets) have significantly pushed forward the state-of-the-art on image classification, increasing in performance as networks grow both deeper and wider. However, memory consumption becomes a bottleneck, as one needs to store the activations in order to calculate gradients using backpropagation. We present the Reversible Residual Network (RevNet), a variant of ResNets where each layer's activations can be reconstructed exactly from the next layer's. Therefore, the activations for most layers need not be stored in memory during backpropagation. We demonstrate the effectiveness of RevNets on CIFAR-10, CIFAR-100, and ImageNet, establishing nearly identical classification accuracy to equally-sized ResNets, even though the activation storage requirements are independent of depth.


title: Discovery of new drug therapeutic indications from gene mutation information for hepatocellular carcinoma

abstract: Hepatocellular carcinoma (HCC) is the most common primary liver malignancy and is a leading cause of cancer-related death worldwide. However, cure is not possible with currently used therapies, and there is not so much approved targeted therapy for HCC despite numerous attempts and clinical trials. So, it is essential to identify additional therapeutic strategies to block the growth of HCC tumors. As a cancer disease, it is associated with aberrant genomic and transcriptional landscapes. We sought to use a systematic drug repositioning bioinformatics approach to identify novel candidate drugs to treat HCC, which considers not only aberrant genomic information, but also the changes of transcriptional landscapes. First, we screen the collection of HCC feature genes that frequently mutated in most samples of HCC based on human mutation data. Then, the gene expression data of HCC in TCGA are combined to classify the kernel genes of HCC. Finally, the therapeutic score (TS) of each drug is calculated based on the kolmogorov-smirnov statistical method. Using this strategy, we identified five drugs that associated with HCC, including three drugs that could treat HCC and two drugs that might have side-effect on HCC. In addition, we also make Connectivity Map (CMap) profiles similarity analysis and KEGG enrichment analysis on drug targets. All these findings suggest that our approach is effective for accurate discovering novel therapeutic options for HCC and easily to be extended to other tumors.


title: A dynamical characterization of acylindrically hyperbolic groups

abstract: We give a dynamical characterization of acylindrically hyperbolic groups. As an application, we prove that non-elementary convergence groups are acylindrically hyperbolic.


title: GLSR-VAE: Geodesic Latent Space Regularization for Variational AutoEncoder Architectures

abstract: VAEs (Variational AutoEncoders) have proved to be powerful in the context of density modeling and have been used in a variety of contexts for creative purposes. In many settings, the data we model possesses continuous attributes that we would like to take into account at generation time. We propose in this paper GLSR-VAE, a Geodesic Latent Space Regularization for the Variational AutoEncoder architecture and its generalizations which allows a fine control on the embedding of the data into the latent space. When augmenting the VAE loss with this regularization, changes in the learned latent space reflects changes of the attributes of the data. This deeper understanding of the VAE latent space structure offers the possibility to modulate the attributes of the generated data in a continuous way. We demonstrate its efficiency on a monophonic music generation task where we manage to generate variations of discrete sequences in an intended and playful way.


title: Game Theory for Secure Critical Interdependent Gas-Power-Water Infrastructure

abstract: A city's critical infrastructure such as gas, water, and power systems, are largely interdependent since they share energy, computing, and communication resources. This, in turn, makes it challenging to endow them with fool-proof security solutions. In this paper, a unified model for interdependent gas-power-water infrastructure is presented and the security of this model is studied using a novel game-theoretic framework. In particular, a zero-sum noncooperative game is formulated between a malicious attacker who seeks to simultaneously alter the states of the gas-power-water critical infrastructure to increase the power generation cost and a defender who allocates communication resources over its attack detection filters in local areas to monitor the infrastructure. At the mixed strategy Nash equilibrium of this game, numerical results show that the expected power generation cost deviation is 35\% lower than the one resulting from an equal allocation of resources over the local filters. The results also show that, at equilibrium, the interdependence of the power system on the natural gas and water systems can motivate the attacker to target the states of the water and natural gas systems to change the operational states of the power grid. Conversely, the defender allocates a portion of its resources to the water and natural gas states of the interdependent system to protect the grid from state deviations.


title: On the Power Spectrum of Dark Matter Substructure in Strong Gravitational Lenses

abstract: Studying the smallest self-bound dark matter structure in our Universe can yield important clues about the fundamental particle nature of dark matter. Galaxy-scale strong gravitational lensing provides a unique way to detect and characterize dark matter substructures at cosmological distances from the Milky Way. Within the cold dark matter (CDM) paradigm, the number of low-mass subhalos within lens galaxies is expected to be large, implying that their contribution to the lensing convergence field is approximately Gaussian and could thus be described by their power spectrum. We develop here a general formalism to compute from first principles the substructure convergence power spectrum for different populations of dark matter subhalos. As an example, we apply our framework to two distinct subhalo populations: a truncated Navarro-Frenk-White subhalo population motivated by standard CDM, and a truncated cored subhalo population motivated by self-interacting dark matter (SIDM). We study in detail how the subhalo abundance, mass function, internal density profile, and concentration affect the amplitude and shape of the substructure power spectrum. We determine that the power spectrum is mostly sensitive to a specific combination of the subhalo abundance and moments of the mass function, as well as to the average tidal truncation scale of the largest subhalos included in the analysis. Interestingly, we show that the asymptotic slope of the substructure power spectrum at large wave number reflects the internal density profile of the subhalos. In particular, the SIDM power spectrum exhibits a characteristic steepening at large wave number absent in the CDM power spectrum, opening the possibility of using this observable, if at all measurable, to discern between these two scenarios.


title: US Cosmic Visions: New Ideas in Dark Matter 2017: Community Report

abstract: This white paper summarizes the workshop "U.S. Cosmic Visions: New Ideas in Dark Matter" held at University of Maryland on March 23-25, 2017.


title: Local Two-Sample Testing: A New Tool for Analysing High-Dimensional Astronomical Data

abstract: Modern surveys have provided the astronomical community with a flood of high-dimensional data, but analyses of these data often occur after their projection to lower-dimensional spaces. In this work, we introduce a local two-sample hypothesis test framework that an analyst may directly apply to data in their native space. In this framework, the analyst defines two classes based on a response variable of interest (e.g. higher-mass galaxies versus lower-mass galaxies) and determines at arbitrary points in predictor space whether the local proportions of objects that belong to the two classes significantly differs from the global proportion. Our framework has a potential myriad of uses throughout astronomy; here, we demonstrate its efficacy by applying it to a sample of 2487 i-band-selected galaxies observed by the HST ACS in four of the CANDELS program fields. For each galaxy, we have seven morphological summary statistics along with an estimated stellar mass and star-formation rate. We perform two studies: one in which we determine regions of the seven-dimensional space of morphological statistics where high-mass galaxies are significantly more numerous than low-mass galaxies, and vice-versa, and another study where we use SFR in place of mass. We find that we are able to identify such regions, and show how high-mass/low-SFR regions are associated with concentrated and undisturbed galaxies while galaxies in low-mass/high-SFR regions appear more extended and/or disturbed than their high-mass/low-SFR counterparts.


title: Unveiling two types of local order in liquid water using machine learning

abstract: Machine learning methods are being explored in many areas of science, with the aim of finding solution to problems that evade traditional scientific approaches due to their complexity. In general, an order parameter capable of identifying two different phases of matter separated by a correspond- ing phase transition is constructed based on symmetry arguments. This parameter measures the degree of order as the phase transition proceeds. However, when the two distinct phases are highly disordered it is not trivial to identify broken symmetries with which to find an order parameter. This poses an excellent problem to be addressed using machine learning procedures. Room tem- perature liquid water is hypothesized to be a supercritical liquid, with fluctuations of two different molecular orders associated to two parent liquid phases, one with high density and another one with low density. The validity of this hypothesis is linked to the existence of an order parameter capable of identifying the two distinct liquid phases and their fluctuations. In this work we show how two different machine learning procedures are capable of recognizing local order in liquid water. We argue that when in order to learn relevant features from this complexity, an initial, physically motivated preparation of the available data is as important as the quality of the data set, and that machine learning can become a successful analysis tool only when coupled to high level physical information.


title: Implications of tides for life on exoplanets

abstract: As evident from the nearby examples of Proxima Centauri and TRAPPIST-1, Earth-sized planets in the habitable zone of low-mass stars are common. Here, we focus on such planetary systems and argue that their (oceanic) tides could be more prominent due to stronger tidal forces. We identify the conditions under which tides may exert a significant positive influence on biotic processes including abiogenesis, biological rhythms, nutrient upwelling and stimulating photosynthesis. We conclude our analysis with the identification of large-scale algal blooms as potential temporal biosignatures in reflectance light curves that can arise indirectly as a consequence of strong tidal forces.


title: One-loop considerations for coexisting vacua in the CP conserving 2HDM

abstract: The Two-Higgs-Doublet model (2HDM) is a simple and viable extension of the Standard Model (SM) with a scalar potential complex enough that two minima may coexist. In this work we investigate if the procedure to identify our vacuum as the global minimum by tree-level formulas carries over to the one-loop corrected potential. In the CP conserving case, we identify two distinct types of coexisting minima --- the regular ones (moderate $\tan\beta$) and the non-regular ones (small or large $\tan\beta$) --- and conclude that the tree level expectation fails only for the non-regular type of coexisting minima. For the regular type, the sign of $m^2_{12}$ already precisely indicates which minima is the global one, even at one-loop.


title: DocTag2Vec: An Embedding Based Multi-label Learning Approach for Document Tagging

abstract: Tagging news articles or blog posts with relevant tags from a collection of predefined ones is coined as document tagging in this work. Accurate tagging of articles can benefit several downstream applications such as recommendation and search. In this work, we propose a novel yet simple approach called DocTag2Vec to accomplish this task. We substantially extend Word2Vec and Doc2Vec---two popular models for learning distributed representation of words and documents. In DocTag2Vec, we simultaneously learn the representation of words, documents, and tags in a joint vector space during training, and employ the simple $k$-nearest neighbor search to predict tags for unseen documents. In contrast to previous multi-label learning methods, DocTag2Vec directly deals with raw text instead of provided feature vector, and in addition, enjoys advantages like the learning of tag representation, and the ability of handling newly created tags. To demonstrate the effectiveness of our approach, we conduct experiments on several datasets and show promising results against state-of-the-art methods.


title: Multifractal Study of Quasiparticle Localization in Disordered Superconductors

abstract: The thermal metal to thermal insulator transition due to random disorder is studied in the context of the symmetries of the Bogoliubov de Gennes Hamiltonian. We focus on a three dimensional system with gapless s-wave pairing that possesses time reversal and spin rotational symmetry. The quasiparticle excitations (bogolons) undergo a metal insulator transition as the disorder increases. We determine the critical disorder strength and correlation exponent first by the transfer matrix method (TMM). We then apply a multifractal finite sized scaling (MFSS) of the bogolon wavefunction obtained from large scale diagonalization of the Hamiltonian and obtain the critical disorder strength and exponent, in agreement with those found by TMM.


title: Nonlinear Programming Methods for Distributed Optimization

abstract: In this paper we investigate how standard nonlinear programming algorithms can be used to solve constrained optimization problems in a distributed manner. The optimization setup consists of a set of agents interacting through a communication graph that have as common goal the minimization of a function expressed as a sum of (possibly non-convex) differentiable functions. Each function in the sum corresponds to an agent and each agent has associated an equality constraint. By re-casting the distributed optimization problem into an equivalent, augmented centralized problem, we show that distributed algorithms result naturally from applying standard nonlinear programming techniques. Due to the distributed formulation, the standard assumptions and convergence results no longer hold. We emphasize what changes are necessary for convergence to still be achieved for three algorithms: two algorithms based on Lagrangian methods, and an algorithm based the method of multipliers. The changes in the convergence results are necessary mainly due to the fact that the (local) minimizers of the lifted optimization problem are not regular, as a results of the distributed formulation. Unlike the standard algorithm based on the method of multipliers, for the distributed version we cannot show that the theoretical superlinear convergence rate can be achieved.


title: Finite size analysis of measurement device independent quantum cryptography with continuous variables

abstract: We study the impact of finite-size effects on the key rate of continuous-variable (CV) measurement-device-independent (MDI) quantum key distribution (QKD). Inspired by the parameter estimation technique developed in [Rupert \textit{et al.} Phys. Rev. A \textbf{90}, 062310 (2014)]~we adapt it to study CV-MDI-QKD and, assuming realistic experimental conditions, we analyze the impact of finite-size effects on the key rate. We find that, increasing the block-size, the performance of the protocol converges towards the ideal one, and that block-sizes between $10^{6}$ and $10^{9}$ data points can already provide a key rate $\sim10^{-2}$ bit/use over metropolitan distances.


title: One Tool, Many Languages: Language-Parametric Transformation with Incremental Parametric Syntax

abstract: We present a new approach for building source-to-source transformations that can run on multiple programming languages, based on a new way of representing programs called incremental parametric syntax. We implement this approach in Haskell in our Cubix system, and construct incremental parametric syntaxes for C, Java, JavaScript, Lua, and Python. We demonstrate a whole-program refactoring tool that runs on all of them, along with three smaller transformations that each run on several. Our evaluation shows that (1) once a transformation is written, little work is required to configure it for a new language (2) transformations built this way output readable code which preserve the structure of the original, according to participants in our human study, and (3) our transformations can still handle language corner-cases, as validated on compiler test suites.


title: A "quite superfluous" characterisation of the Jacobson radical of an associative ring

abstract: It is well-known that the Jacobson radical of a unital ring is its largest superfluous right ideal. It is recorded here that the result carries over to non-unital rings provided the notion of "superfluous" is taken relative to all regular ideals.


title: An Efficient and Distribution-Free Two-Sample Test Based on Energy Statistics and Random Projections

abstract: A common disadvantage in existing distribution-free two-sample testing approaches is that the computational complexity could be high. Specifically, if the sample size is $N$, the computational complexity of those two-sample tests is at least $O(N^2)$. In this paper, we develop an efficient algorithm with complexity $O(N \log N)$ for computing energy statistics in univariate cases. For multivariate cases, we introduce a two-sample test based on energy statistics and random projections, which enjoys the $O(K N \log N)$ computational complexity, where $K$ is the number of random projections. We name our method for multivariate cases as Randomly Projected Energy Statistics (RPES). We can show RPES achieves nearly the same test power with energy statistics both theoretically and empirically. Numerical experiments also demonstrate the efficiency of the proposed method over the competitors.


title: An intercomparison of monthly surface air temperature on islands and proximate moorings across the tropical Indo-Pacific

abstract: Despite the importance and vulnerability of small island nations and their ecosystems, they frequently have insufficient observations to provide baseline or long-term perspectives on climate variability and change, and global model experiments rarely have the resolution to include them. Many studies in observational climatology, climate modeling, and paleoclimate thus depend to varying degrees on an approximation equating near-surface marine and terrestrial island climates, often with direct implications for island societies and resources. Here we investigate the validity of this approximation, and by extension the viability of offshore moorings to serve as proxies for island climate variability, by comparing monthly mean observations on a diverse set of 17 islands across the global tropics with those from ocean moorings proximate to each island. While some island-mooring pairs exhibit a mean offset in surface air temperature, these cannot be explained as a simple function of one parameter including distance from the mooring, station elevation, or island dimensions. Overall, the seasonal to interannual variability in near-surface climate at monthly time scales at meteorological stations on tropical islands is captured remarkably well by moorings positioned up to 1,000 km offshore.


title: Mesoscale simulations of confined Nafion thin films

abstract: The morphology and transport properties of thin films of the ionomer Nafion, with thicknesses on the order of the bulk cluster size, have been investigated as a model system to explain the anomalous behaviour of catalyst/electrode-polymer interfaces in membrane-electrode assemblies. We have employed dissipative particle dynamics (DPD) to investigate the interaction of water and fluorocarbon chains with carbon and quartz as confining materials for a wide range of operational water contents and film thicknesses. We found confinement-induced clustering of water perpendicular to the thin film. Hydrophobic carbon forms a water depletion zone near the film interface, whereas hydrophilic quartz results in a zone with excess water. There are, on average, oscillating water-rich and fluorocarbon-rich regions, in agreement with experimental results from neutron reflectometry. Water diffusivity shows increasing directional anisotropy of up to 30% with decreasing film thickness, depending of the confining material. The percolation analysis revealed significant differences in water clustering and connectivity with the confining material. These findings indicate the fundamentally different nature of ionomer thin films, compared to membranes, and suggest explanations for increased ionic resistances observed in the catalyst layer.


title: Multigrid-based inversion for volumetric radar imaging with asteroid interior reconstruction as a potential application

abstract: This study concentrates on advancing mathematical and computational methodology for radar tomography imaging in which the unknown volumetric velocity distribution of a wave within a bounded domain is to be reconstructed. Our goal is to enable effective simulation and inversion of a large amount of full-wave data within a realistic 2D or 3D geometry. For propagating and inverting the wave, we present a rigorous multigrid-based forward approach which utilizes the finite-difference time-domain method and a nested finite element grid structure. Based on the multigrid approach, we introduce and validate a multiresolution algorithm which allows regularization of the unknown distribution through a coarse-to-fine inversion scheme. In this approach, sparse signals can be effectively inverted, as the coarse fluctuations are reconstructed before the finer ones. Furthermore, the number of nonzero entries in the system matrix can be compressed and thus the inversion procedure can be speeded up. As a test scenario we investigate satellite-based asteroid interior reconstruction. We use both full-wave and projected wave data and estimate the accuracy of the inversion under different error sources: noise and positioning inaccuracies. The results suggest that the present full-wave inversion approach allows recovering the interior with a single satellite recording backscattering data. It seems that robust results can be achieved, when the peak-to-peak signal-to-noise ratio is above 10 dB. Furthermore, it seems that reconstructing the deep interior can be enhanced if two satellites can be utilized in the measurements.


title: The Discovery of a Gravitationally Lensed Supernova Ia at Redshift 2.22

abstract: We present the discovery and measurements of a gravitationally lensed supernova (SN) behind the galaxy cluster MOO J1014+0038. Based on multi-band Hubble Space Telescope and Very Large Telescope (VLT) photometry of the supernova, and VLT spectroscopy of the host galaxy, we find a 97.5% probability that this SN is a SN Ia, and a 2.5% chance of a CC SN. Our typing algorithm combines the shape and color of the light curve with the expected rates of each SN type in the host galaxy. With a redshift of 2.2216, this is the highest redshift SN Ia discovered with a spectroscopic host-galaxy redshift. A further distinguishing feature is that the lensing cluster, at redshift 1.23, is the most distant to date to have an amplified SN. The SN lies in the middle of the color and light-curve shape distributions found at lower redshift, disfavoring strong evolution to z = 2.22. We estimate an amplification due to gravitational lensing of 2.8+0.6-0.5 (1.10 +- 0.23 mag)---compatible with the value estimated from the weak-lensing-derived mass and the mass-concentration relation from LambdaCDM simulations---making it the most amplified SN Ia discovered behind a galaxy cluster.


title: Identification of Overlapping Communities via Constrained Egonet Tensor Decomposition

abstract: Detection of overlapping communities in real-world networks is a generally challenging task. Upon recognizing that a network is in fact the union of its egonets, a novel network representation using multi-way data structures is advocated in this contribution. The introduced sparse tensor-based representation exhibits richer structure compared to its matrix counterpart, and thus enables a more robust approach to community detection. To leverage this structure, a constrained tensor approximation framework is introduced using PARAFAC decomposition. The arising constrained trilinear optimization is handled via alternating minimization, where intermediate subproblems are solved using the alternating direction method of multipliers (ADMM) to ensure convergence. The factors obtained provide soft community memberships, which can further be exploited for crisp, and possibly-overlapping community assignments. The framework is further broadened to include time-varying graphs, where the edgeset as well as the underlying communities evolve through time. Performance of the proposed approach is assessed via tests on benchmark synthetic graphs as well as real-world networks. As corroborated by numerical tests, the proposed tensor-based representation captures multi-hop nodal connections, that is, connectivity patterns within single-hop neighbors, whose exploitation yields a more robust community identification in the presence of mixing as well as overlapping communities.


title: Symmetry energies for $A = 24$ and $48$ and the USD and KB3 shell model Hamiltonians

abstract: Calculations in the sd and pf shells reported some time ago by Satu\l a\etal\ [Phys.~Lett.~B~407, 103 (1997)] are redone for an extended analysis of the results. As in the original work, we do calculations for one mass number in each shell and consider in each case the sequence of lowest energies for isospins 0, 2, and 4, briefly the symmetry spectrum. Following further the original work we study how this spectrum changes when parts of the two-nucleon interaction are turned off. The variation of its width is explored in detail. A differential combination $\epsilon_\text{W}$ of the three energies was taken in the original work as a measure of the so-called Wigner term in semi-empirical mass formulas, and it was found to decrease drastically when the two-nucleon interaction in the channel of zero isospin is turned off. Our analysis shows that the width of the symmetry spectrum experiences an equally drastic decrease, which can be explained qualitatively in terms of schematic approximations. We therefore suggest that the decrease of $\epsilon_\text{W}$ be seen mainly as a side effect of a narrowing of the symmetry spectrum rather than an independent manifestation of the two-nucleon interaction in the channel of zero isospin.


title: Fine-grained reductions from approximate counting to decision

abstract: In this paper, we introduce a general framework for fine-grained reductions of approximate counting problems to their decision versions. (Thus we use an oracle that decides whether any witness exists to multiplicatively approximate the number of witnesses with minimal overhead.) This mirrors a foundational result of Sipser (STOC 1983) and Stockmeyer (SICOMP 1985) in the polynomial-time setting, and a similar result of M\"uller (IWPEC 2006) in the FPT setting. Using our framework, we obtain such reductions for some of the most important problems in fine-grained complexity: the Orthogonal Vectors problem, 3SUM, and the Negative-Weight Triangle problem (which is closely related to All-Pairs Shortest Path). We also provide a fine-grained reduction from approximate #SAT to SAT. Suppose the Strong Exponential Time Hypothesis (SETH) is false, so that for some $1<c<2$ and all $k$ there is an $O(c^n)$-time algorithm for k-SAT. Then we prove that for all $k$, there is an $O((c+o(1))^n)$-time algorithm for approximate #$k$-SAT. In particular, our result implies that the Exponential Time Hypothesis (ETH) is equivalent to the seemingly-weaker statement that there is no algorithm to approximate #3-SAT to within a factor of $1+\epsilon$ in time $2^{o(n)}/\epsilon^2$ (taking $\epsilon > 0$ as part of the input).


title: Cloud-based or On-device: An Empirical Study of Mobile Deep Inference

abstract: Modern mobile applications are benefiting significantly from the advancement in deep learning, e.g., implementing real-time image recognition and conversational system. Given a trained deep learning model, applications usually need to perform a series of matrix operations based on the input data, in order to infer possible output values. Because of computational complexity and size constraints, these trained models are often hosted in the cloud. To utilize these cloud-based models, mobile apps will have to send input data over the network. While cloud-based deep learning can provide reasonable response time for mobile apps, it restricts the use case scenarios, e.g. mobile apps need to have network access. With mobile specific deep learning optimizations, it is now possible to employ on-device inference. However, because mobile hardware, such as GPU and memory size, can be very limited when compared to its desktop counterpart, it is important to understand the feasibility of this new on-device deep learning inference architecture. In this paper, we empirically evaluate the inference performance of three Convolutional Neural Networks (CNNs) using a benchmark Android application we developed. Our measurement and analysis suggest that on-device inference can cost up to two orders of magnitude greater response time and energy when compared to cloud-based inference, and that loading model and computing probability are two performance bottlenecks for on-device deep inferences.


title: Temporal Simultons in Optical Parametric Oscillators

abstract: We identify and demonstrate a regime of operation in optical parametric oscillators (OPOs) in which the formation of temporal simultons produces stable half-harmonic pulses. Simultons are simultanous bright-dark solitons of a signal field at frequency $\omega$ and the pump field at $2\omega$ which form in a quadratic nonlinear medium. The formation of simultons in an OPO is evidenced by sech$^2$ spectra with broad instantaneous bandwidths which increase with pump power, and large slope efficiencies. In contrast to conventional synchronously pumped OPOs, operation in this regime is achieved by using relatively large parametric gains, and a low finesse resonator detuned to a slightly longer roundtrip time than the pump repetition period. In the experiment we achieve sub-50 fs sech$^2$ pulses with a slope efficiency of 570%, and a conversion efficiency of 55% as a result of formation of simultons in an OPO. We verify the distinct features of OPO operation in the simulton regime analytically and numerically. These results represent a fundamental shift in the understanding and design of OPOs, in which the nonlinear dynamics that arise during asynchronous operation can be used to efficiently generate few-cycle long-wavelength frequency combs.


title: Group sequential designs for negative binomial outcomes

abstract: Count data and recurrent events in clinical trials, such as the number of lesions in magnetic resonance imaging in multiple sclerosis, the number of relapses in multiple sclerosis, the number of hospitalizations in heart failure, and the number of exacerbations in asthma or in chronic obstructive pulmonary disease (COPD) are often modeled by negative binomial distributions. In this manuscript we study planning and analyzing clinical trials with group sequential designs for negative binomial outcomes. We propose a group sequential testing procedure for negative binomial outcomes based on Wald statistics using maximum likelihood estimators. The asymptotic distribution of the proposed group sequential tests statistics are derived. The finite sample size properties of the proposed group sequential test for negative binomial outcomes and the methods for planning the respective clinical trials are assessed in a simulation study. The simulation scenarios are motivated by clinical trials in chronic heart failure and relapsing multiple sclerosis, which cover a wide range of practically relevant settings. Our research assures that the asymptotic normal theory of group sequential designs can be applied to negative binomial outcomes when the hypotheses are tested using Wald statistics and maximum likelihood estimators. We also propose two methods, one based on Student's t-distribution and one based on resampling, to improve type I error rate control in small samples. The statistical methods studied in this manuscript are implemented in the R package \textit{gscounts}, which is available for download on the Comprehensive R Archive Network (CRAN).


title: Spacetime supersymmetry in low-dimensional perturbative heterotic compactifications

abstract: We study the constraints of spacetime supersymmetry for perturbative three- and two-dimensional Minkowski vacua of the critical heterotic string. Assuming a standard RNS construction of the spacetime supersymmetry generators and a compact unitary internal superconformal worldsheet theory, we describe the worldsheet structures associated to various spacetime supersymmetries. In three dimensions we show that there are no CFT surprises: each allowed spacetime supersymmetry is realized by a supergravity compactification. As a recent orbifold construction shows, in two dimensions there are more exotic possibilities, and we discuss how these fit into our analysis.


title: Sucralose Interaction with Protein Structures

abstract: Sucralose is a commonly employed artificial sweetener that appears to destabilize protein native structures. This is in direct contrast to the bio-preservative nature of its natural counterpart, sucrose, which enhances the stability of biomolecules against environmental stress. We have further explored the molecular interactions of sucralose as compared to sucrose to illuminate the origin of the differences in their bio-preservative efficacy. We show that the mode of interactions of sucralose and sucrose in bulk solution differ subtly using hydration dynamics measurement and computational simulation. Sucralose does not appear to disturb the native state of proteins for moderate concentrations (<0.2 M) at room temperature. However, as the concentration increases, or in the thermally stressed state, sucralose appears to differ in its interactions with protein leading to the reduction of native state stability. This difference in interaction appears weak. We explored the difference in the preferential exclusion model using time-resolved spectroscopic techniques and observed that both molecules appear to be effective reducers of bulk hydration dynamics. However, the chlorination of sucralose appears to slightly enhance the hydrophobicity of the molecule, which reduces the preferential exclusion of sucralose from the protein-water interface. The weak interaction of sucralose with hydrophobic pockets on the protein surface differs from the behavior of sucrose. We experimentally followed up upon the extent of this weak interaction using isothermal titration calorimetry (ITC) measurements. We propose this as a possible origin for the difference in their bio-preservative properties.


title: On the Complexity of Learning Neural Networks

abstract: The stunning empirical successes of neural networks currently lack rigorous theoretical explanation. What form would such an explanation take, in the face of existing complexity-theoretic lower bounds? A first step might be to show that data generated by neural networks with a single hidden layer, smooth activation functions and benign input distributions can be learned efficiently. We demonstrate here a comprehensive lower bound ruling out this possibility: for a wide class of activation functions (including all currently used), and inputs drawn from any logconcave distribution, there is a family of one-hidden-layer functions whose output is a sum gate, that are hard to learn in a precise sense: any statistical query algorithm (which includes all known variants of stochastic gradient descent with any loss function) needs an exponential number of queries even using tolerance inversely proportional to the input dimensionality. Moreover, this hard family of functions is realizable with a small (sublinear in dimension) number of activation units in the single hidden layer. The lower bound is also robust to small perturbations of the true weights. Systematic experiments illustrate a phase transition in the training error as predicted by the analysis.


title: Intertwining wavelets or Multiresolution analysis on graphs through random forests

abstract: We propose a new method for performing multiscale analysis of functions defined on the vertices of a finite connected weighted graph. Our approach relies on a random spanning forest to downsample the set of vertices, and on approximate solutions of Markov intertwining relation to provide a subgraph structure and a filter bank leading to a wavelet basis of the set of functions. Our construction involves two parameters q and q'. The first one controls the mean number of kept vertices in the downsampling, while the second one is a tuning parameter between space localization and frequency localization. We provide an explicit reconstruction formula, bounds on the reconstruction operator norm and on the error in the intertwining relation, and a Jackson-like inequality. These bounds lead to recommend a way to choose the parameters q and q'. We illustrate the method by numerical experiments.


title: A Real-Time Solver For Time-Optimal Control Of Omnidirectional Robots with Bounded Acceleration

abstract: We are interested in the problem of time-optimal control of omnidirectional robots with bounded acceleration (TOC-ORBA). While there exist approximate solutions for such robots, and exact solutions with unbounded acceleration, exact solvers to the TOC-ORBA problem have remained elusive until now. In this paper, we present a real-time solver for true time-optimal control of omnidirectional robots with bounded acceleration. We first derive the general parameterized form of the solution to the TOC-ORBA problem by application of Pontryagin's maximum principle. We then frame the boundary value problem of TOC-ORBA as an optimization problem over the parametrized control space. To overcome local minima and poor initial guesses to the optimization problem, we introduce a two-stage optimal control solver (TSOCS): The first stage computes an upper bound to the total time for the TOC-ORBA problem and holds the time constant while optimizing the parameters of the trajectory to approach the boundary value conditions. The second stage uses the parameters found by the first stage, and relaxes the constraint on the total time to solve for the parameters of the complete TOC-ORBA problem. We further implement TSOCS as a closed loop controller to overcome actuation errors on real robots in real-time. We empirically demonstrate the effectiveness of TSOCS in simulation and on real robots, showing that 1) it runs in real time, generating solutions in less than 0.5ms on average; 2) it generates faster trajectories compared to an approximate solver; and 3) it is able to solve TOC-ORBA problems with non-zero final velocities that were previously unsolvable in real-time.


title: Communication Lower Bounds of Bilinear Algorithms for Symmetric Tensor Contractions

abstract: We introduce a new theoretical framework for deriving lower bounds on data movement in bilinear algorithms. Bilinear algorithms are a general representation of fast algorithms for bilinear functions, which include computation of matrix multiplication, convolution, and symmetric tensor contractions. A bilinear algorithm is described by three matrices. Our communication lower bounds are based on quantifying the minimal matrix ranks of matching subsets of columns of these matrices. This infrastructure yields new communication lower bounds for symmetric tensor contraction algorithms, which provide qualitative new insights. Tensor symmetry (invariance under permutation of modes) is common to many applications of tensor computations (e.g., tensor representation of hypergraphs, analysis of high order moments in data, as well as tensors modelling interactions of electrons in computational chemistry). Tensor symmetry enables reduction in representation size as well as arithmetic cost of contractions by factors that scale with the number of equivalent permutations. However, we derive lower bounds showing that these arithmetic cost and memory reductions can necessitate increases in data movement by factors that scale with the size of the tensors.


title: Simplified Long Short-term Memory Recurrent Neural Networks: part I

abstract: We present five variants of the standard Long Short-term Memory (LSTM) recurrent neural networks by uniformly reducing blocks of adaptive parameters in the gating mechanisms. For simplicity, we refer to these models as LSTM1, LSTM2, LSTM3, LSTM4, and LSTM5, respectively. Such parameter-reduced variants enable speeding up data training computations and would be more suitable for implementations onto constrained embedded platforms. We comparatively evaluate and verify our five variant models on the classical MNIST dataset and demonstrate that these variant models are comparable to a standard implementation of the LSTM model while using less number of parameters. Moreover, we observe that in some cases the standard LSTM's accuracy performance will drop after a number of epochs when using the ReLU nonlinearity; in contrast, however, LSTM3, LSTM4 and LSTM5 will retain their performance.


title: Gaussian Estimates: A Brief History

abstract: Two-sided Gaussian estimates for the fundamental solution of a second order linear parabolic differential equation are upper and lower bounds in terms of the fundamental solution of the classical heat conduction equation. In his seminal 1958 paper Nash stated, without proof, two-sided non-Gaussian bounds for the fundamental solution of a uniformly parabolic divergence structure equation assuming only boundedness of the coefficients. In his 1967-1968 papers Aronson derived truly Gaussian estimates for the fundamental solutions of a large class of linear parabolic equations (including the divergence structure equation) under minimal non-regularity assumptions on the coefficients. Subsequently in 1986 Fabes & Stroock derived Gaussian estimates for the divergence structure equation directly from the ideas of Nash and went on to prove Nash's continuity theorem and the Harnack inequality as a consequence of their estimate. In this note I describe these results together with various extensions.


title: UAV Air-to-Ground Channel Characterization for mmWave Systems

abstract: Communication at mmWave bands carries critical importance for 5G wireless networks. In this paper, we study the characterization of mmWave air-to-ground (AG) channels for unmanned aerial vehicle (UAV) communications. In particular, we use ray tracing simulations using Remcom Wireless InSite software to study the behavior of AG mmWave bands at two different frequencies: 28~GHz and 60~GHz. Received signal strength (RSS) and root mean square delay spread (RMS-DS) of multipath components (MPCs) are analyzed for different UAV heights considering four different environments: urban, suburban, rural, and over sea. It is observed that the RSS mostly follows the two ray propagation model along the UAV flight path for higher altitudes. This two ray propagation model is affected by the presence of high rise scatterers in urban scenario. Moreover, we present details of a universal serial radio peripheral (USRP) based channel sounder that can be used for AG channel measurements for mmWave (60 GHz) UAV communications.


title: Smooth backfitting of proportional hazards with multiplicative components

abstract: Smooth backfitting has proven to have a number of theoretical and practical advantages in structured regression. Smooth backfitting projects the data down onto the structured space of interest providing a direct link between data and estimator. This paper introduces the ideas of smooth backfitting to survival analysis in a proportional hazard model, where we assume an underlying conditional hazard with multiplicative components. We develop asymptotic theory for the estimator and we use the smooth backfitter in a practical application, where we extend recent advances of in-sample forecasting methodology by allowing more information to be incorporated, while still obeying the structured requirements of in-sample forecasting.


title: Simplified Long Short-term Memory Recurrent Neural Networks: part II

abstract: This is part II of three-part work. Here, we present a second set of inter-related five variants of simplified Long Short-term Memory (LSTM) recurrent neural networks by further reducing adaptive parameters. Two of these models have been introduced in part I of this work. We evaluate and verify our model variants on the benchmark MNIST dataset and assert that these models are comparable to the base LSTM model while use progressively less number of parameters. Moreover, we observe that in case of using the ReLU activation, the test accuracy performance of the standard LSTM will drop after a number of epochs when learning parameter become larger. However all of the new model variants sustain their performance.


title: Smoothing of limit linear series on curves and metrized complexes of pseudocompact type

abstract: We investigate the connection between Osserman limit series (on curves of pseudocompact type) and Amini-Baker limit linear series (on metrized complexes with corresponding underlying curve) via a notion of pre-limit linear series on curves of the same type. Then, applying the smoothing theorems of Osserman limit linear series, we deduce that, fixing certain metrized complexes, or for certain types of Amini-Baker limit linear series, the smoothability is equivalent to a certain "weak glueing condition". Also for arbitrary metrized complexes of pseudocompact type the weak glueing condition (when it applies) is necessary for smoothability. As an application we confirm the lifting property of specific divisors on the metric graph associated to a certain regular smoothing family, and give a new proof of the main theorem in the paper of Cartwright, Jensen and Payne for vertex avoiding divisors, and generalize loc.cit. for divisors of rank one in the sense that, for the metric graph, there could be at most three edges (instead of two) between any pair of adjacent vertices.


title: Alloying strategy for two-dimensional GaN optical emitters

abstract: The recent progress in formation of two-dimensional (2D) GaN by a migration-enhanced encapsulated technique opens up new possibilities for group III-V 2D semiconductors with a band gap within the visible energy spectrum. Using first-principles calculations we explored alloying of 2D-GaN to achieve an optically active material with a tuneable band gap. The effect of isoelectronic III-V substitutional elements on the band gaps, band offsets, and spatial electron localization is studied. In addition to optoelectronic properties, the formability of alloys is evaluated using impurity formation energies. A dilute highly-mismatched solid solution 2D-GaN$_{1-x}$P$_x$ features an efficient band gap reduction in combination with a moderate energy penalty associated with incorporation of phosphorous in 2D-GaN, which is substantially lower than in the case of the bulk GaN. The group-V alloying elements also introduce significant disorder and localization at the valence band edge that facilitates direct band gap optical transitions thus implying the feasibility of using III-V alloys of 2D-GaN in light-emitting devices.


title: Simplified Long Short-term Memory Recurrent Neural Networks: part III

abstract: This is part III of three-part work. In parts I and II, we have presented eight variants for simplified Long Short Term Memory (LSTM) recurrent neural networks (RNNs). It is noted that fast computation, specially in constrained computing resources, are an important factor in processing big time-sequence data. In this part III paper, we present and evaluate two new LSTM model variants which dramatically reduce the computational load while retaining comparable performance to the base (standard) LSTM RNNs. In these new variants, we impose (Hadamard) pointwise state multiplications in the cell-memory network in addition to the gating signal networks.


title: Lacunary Eta-quotients Modulo Powers of Primes

abstract: An integral power series is called lacunary modulo $M$ if almost all of its coefficients are divisible by $M$. Motivated by the parity problem for the partition function, $p(n)$, Gordon and Ono studied the generating functions for $t$-regular partitions, and determined conditions for when these functions are lacunary modulo powers of primes. We generalize their results in a number of ways by studying infinite products called Dedekind eta-quotients and generalized Dedekind eta-quotients. We then apply our results to the generating functions for the partition functions considered by Nekrasov, Okounkov, and Han.


title: Quantum Gravity inspired nonlocal gravity model

abstract: We consider a nonlocal gravity model motivated by nonperturbative Quantum Gravity studies. This model, if correct, suggests the existence of strong IR relevant effects which can lead to an interesting late time cosmology. We implement the IR modification directly in the effective action. We show that, upon some assumptions on initial conditions, this model describes an observationally viable background cosmology being also consistent with local gravity tests.


title: Compliant Movement Primitives in a Bimanual Setting

abstract: Simultaneously achieving low trajectory errors and compliant control \emph{without} explicit models of the task was effectively addressed with Compliant Movement Primitives (CMP). For a single-robot task, this means that it is accurately following its trajectory, but also exhibits compliant behavior in case of perturbations. In this paper we extend this kind of behavior without explicit models to bimanual tasks. In the presence of an external perturbation on any of the robots, they will both move in synchrony in order to maintain their relative posture, and thus not exert force on the object they are carrying. Thus, they will act compliantly in their absolute task, but remain stiff in their relative task. To achieve compliant absolute behavior and stiff relative behavior, we combine joint-space CMPs with the well known symmetric control approach. To reduce the necessary feedback reaction of symmetric control, we further augment it with copying of a virtual force vector at the end-effector, calculated through the measured external joint torques. Real-world results on two Kuka LWR-4 robots in a bimanual setting confirm the applicability of the approach.


title: A Class of Integrable Metrics and Gauge Fields

abstract: Starting with the most general four-dimensional spacetime possessing two commuting Killing vectors and a nontrivial Killing tensor, we analytically integrate Einstein-Yang-Mills equations for a completely arbitrary gauge group. It is assumed that the gauge field inherits the symmetries of the background and is aligned with the principal null directions of the spacetime.


title: Nonperturbative quark-gluon thermodynamics at finite density

abstract: Thermodynamics of the quark-gluon plasma at finite density is studied in the framework of the Field Correlator Method, where thermodynamical effects of Polyakov loops and colormagnetic confinement are taken into account. Having found good agreement with numerical lattice data for zero density, we calculate pressure $P(T,\,\mu)$ for $0 <\mu< 400$ MeV and $150 <T< 1000$ MeV. For the first time the explicit intergral form is found in this region,demonstrating analytic structure in the complex $\mu$ plane. The resulting multiple complex branch points are found at the Roberge-Weiss values of $\operatorname{Im}\,\mu$, with $\operatorname{Re}\,\mu$ defined by the values of Polyakov lines and colormagnetic confinement.


title: Phenomenological theories of the low-temperature pseudogap: Hall number, specific heat and Seebeck coefficient

abstract: Since its experimental discovery, many phenomenological theories successfully reproduced the rapid rise from $p$ to $1+p$ found in the Hall number $n_H$ at the critical doping $p^*$ of the pseudogap in superconducting cuprates. Further comparison with experiments is now needed in order to narrow down candidates. In this paper, we consider three previously successful phenomenological theories in a unified formalism---an antiferromagnetic mean field (AF), a spiral incommensurate antiferromagnetic mean field (sAF), and the Yang-Rice-Zhang (YRZ) theory. We find a rapid rise in the specific heat and a rapid drop in the Seebeck coefficient for increasing doping across the transition in each of those models. The predicted rises and drops are locked, not to~$p^*$, but to the doping where anti-nodal electron pockets, characteristic of each model, appear at the Fermi surface shortly before~$p^*$. While such electron pockets are still to be found in experiments, we discuss how they could provide distinctive signatures for each model. We also show that the range of doping where those electron pockets would be found is strongly affected by the position of the van~Hove singularity.


title: Diagrammatic Exciton Basis Theory of the Photophysics of Pentacene Dimers

abstract: Covalently linked acene dimers are of interest as candidates for intramolecular singlet fission. We report many-electron calculations of the energies and wavefunctions of the optical singlets, the lowest triplet exciton and the triplet-triplet biexciton, as well as the final states of excited state absorptions from these states in a family of phenyl-linked pentacene dimers. While it is difficult to distinguish between the triplet and the triplet-triplet from their transient absorptions in the 500-600 nm region, by comparing theoretical transient absorption spectra against published and unpublished experimental transient absorptions in the near and mid infrared we conclude that the end product of photoexcitation in these particular bipentacenes is the bound triplet-triplet and not free triplets. We predict additional transient absorptions at even longer wavelengths, beyond 1500 nm, to the equivalent of the classic 2$^1$A$_g^-$ in linear polyenes.


title: Constructing new nonlinear evolution equations with supersymmetry

abstract: The factorisation method commonly used in linear supersymmetric quantum mechanics is extended, such that it can be applied to nonlinear quantum mechanical systems. The new method is distinguishable from the linear formalism, as the superpotential is forced to become eigenfunction-dependent. An example solution is given for the nonlinear Schroedinger equation and its supersymmetric partner equation. This method allows new nonlinear evolution equations to be constructed from the solutions of known nonlinear equations, and has the potential to be a useful tool for mathematicians and physicists working in the field of nonlinear systems, allowing the discovery of previously unknown `dualities' amongst soliton solutions and their respective equations.

categories: math-ph hep-th math.MP nlin.SI physics.optics


title: Periodically stationary multivariate autoregressive models

abstract: A class of multivariate periodic autoregressive models is proposed where coupling between time series is achieved through linear mean functions. Various response distributions with quadratic mean-variance relationships fit into the framework, including the negative binomial, gamma and Gaussian distributions. We develop an iterative algorithm to obtain unconditional means, variances and auto-/cross-covariances for models with higher order lags. Analytical solutions are given for the univariate model with lag one and multivariate models with linear mean-variance relationship. A special case of the model class is an established framework for modelling multivariate time series of counts from routine surveillance of infectious diseases. We extend this model class to allow for distributed lags and apply it to a dataset on norovirus gastroenteritis in two German states. The availability of unconditional moments and auto/cross-correlations enhances model assessment and interpretation.


title: SU(3) structures on S2 bundles over four-manifolds

abstract: We construct globally-defined $SU(3)$ structures on smooth compact toric varieties (SCTV) in the class of $\mathbb{CP}^1$ bundles over $M$, where $M$ is an arbitrary SCTV of complex dimension two. The construction can be extended to the case where the base is K\"ahler-Einstein of positive curvature, but not necessarily toric, and admits a parameter space which includes $SU(3)$ structures of LT type.


title: Impact of Bayesian prior on the characterization of binary black hole coalescences

abstract: In a regime where data are only mildly informative, prior choices can play a significant role in Bayesian statistical inference, potentially affecting the inferred physics. We show this is indeed the case for some of the parameters inferred from current gravitational-wave measurements of binary black hole coalescences. We reanalyze the first detections performed by the twin LIGO interferometers using alternative (and astrophysically motivated) prior assumptions. We find different prior distributions can introduce deviations in the resulting posteriors that impact the physical interpretation of these systems. For instance, (i) limits on the $90\%$ credible interval on the effective black hole spin $\chi_{\rm eff}$ are subject to variations of $\sim 10\%$ if a prior with black hole spins mostly aligned to the binary's angular momentum is considered instead of the standard choice of isotropic spin directions, and (ii) under priors motivated by the initial stellar mass function, we infer tighter constraints on the black hole masses, and in particular, we find no support for any of the inferred masses within the putative mass gap $M \lesssim 5 M_\odot$.


title: Predicting multicellular function through multi-layer tissue networks

abstract: Motivation: Understanding functions of proteins in specific human tissues is essential for insights into disease diagnostics and therapeutics, yet prediction of tissue-specific cellular function remains a critical challenge for biomedicine. Results: Here we present OhmNet, a hierarchy-aware unsupervised node feature learning approach for multi-layer networks. We build a multi-layer network, where each layer represents molecular interactions in a different human tissue. OhmNet then automatically learns a mapping of proteins, represented as nodes, to a neural embedding based low-dimensional space of features. OhmNet encourages sharing of similar features among proteins with similar network neighborhoods and among proteins activated in similar tissues. The algorithm generalizes prior work, which generally ignores relationships between tissues, by modeling tissue organization with a rich multiscale tissue hierarchy. We use OhmNet to study multicellular function in a multi-layer protein interaction network of 107 human tissues. In 48 tissues with known tissue-specific cellular functions, OhmNet provides more accurate predictions of cellular function than alternative approaches, and also generates more accurate hypotheses about tissue-specific protein actions. We show that taking into account the tissue hierarchy leads to improved predictive power. Remarkably, we also demonstrate that it is possible to leverage the tissue hierarchy in order to effectively transfer cellular functions to a functionally uncharacterized tissue. Overall, OhmNet moves from flat networks to multiscale models able to predict a range of phenotypes spanning cellular subsystems


title: Predictive Liability Models and Visualizations of High Dimensional Retail Employee Data

abstract: Employee theft and dishonesty is a major contributor to loss in the retail industry. Retailers have reported the need for more automated analytic tools to assess the liability of their employees. In this work, we train and optimize several machine learning models for regression prediction and analysis on this data, which will help retailers identify and manage risky employees. Since the data we use is very high dimensional, we use feature selection techniques to identify the most contributing factors to an employee's assessed risk. We also use dimension reduction and data embedding techniques to present this dataset in a easy to interpret format.


title: Exactly solvable model for drift of suspended ferromagnetic particles induced by the Magnus force

abstract: The phenomenon of drift motion of single-domain ferromagnetic particles induced by the Magnus force in a viscous fluid is studied analytically. We use a minimal set of equations to describe the translational and rotational motions of these particles subjected to a harmonic force and a non-uniformly rotating magnetic field. Assuming that the azimuthal angle of the magnetic field is a periodic triangular function, we analytically solve the rotational equation of motion in the steady state and calculate the drift velocity of particles. We study in detail the dependence of this velocity on the model parameters, discuss the applicability of the drift phenomenon for separation of particles in suspensions, and verify numerically the analytical predictions.


title: Chiral symmetry-breaking dynamics in the phase transformation of nematic droplets

abstract: Dynamic simulations of the isotropic-nematic phase transformation of liquid crystal droplets under homeotropic anchoring are found to predict chiral symmetry-breaking dynamics. These observations occur when using material parameters for pentyl-cyanobiphenyl (5CB) but not under the single elastic constant approximation of this material, frequently used in simulation. The twisting dynamic process occurs during the relaxation of the domain from an unstable radial texture to a stable uniform texture and involves simultaneous defect loop motion and twisting of the bulk nematic texture.


title: Recognizing Abnormal Heart Sounds Using Deep Learning

abstract: The work presented here applies deep learning to the task of automated cardiac auscultation, i.e. recognizing abnormalities in heart sounds. We describe an automated heart sound classification algorithm that combines the use of time-frequency heat map representations with a deep convolutional neural network (CNN). Given the cost-sensitive nature of misclassification, our CNN architecture is trained using a modified loss function that directly optimizes the trade-off between sensitivity and specificity. We evaluated our algorithm at the 2016 PhysioNet Computing in Cardiology challenge where the objective was to accurately classify normal and abnormal heart sounds from single, short, potentially noisy recordings. Our entry to the challenge achieved a final specificity of 0.95, sensitivity of 0.73 and overall score of 0.84. We achieved the greatest specificity score out of all challenge entries and, using just a single CNN, our algorithm differed in overall score by only 0.02 compared to the top place finisher, which used an ensemble approach.


title: Set-Direct Factorizations of Groups

abstract: We consider factorizations $G=XY$ where $G$ is a general group, $X$ and $Y$ are normal subsets of $G$ and any $g\in G$ has a unique representation $g=xy$ with $x\in X$ and $y\in Y$. This definition coincides with the customary and extensively studied definition of a direct product decomposition by subsets of a finite abelian group. Our main result states that a group $G$ has such a factorization if and only if $G$ is a central product of $\left\langle X\right\rangle $ and $\left\langle Y\right\rangle $ and the central subgroup $\left\langle X\right\rangle \cap \left\langle Y\right\rangle $ satisfies certain abelian factorization conditions. We analyze some special cases and give examples. In particular, simple groups have no non-trivial set-direct factorization.


title: LSQ14efd: observations of the cooling of a shock break-out event in a type Ic Supernova

abstract: We present the photometric and spectroscopic evolution of the type Ic supernova LSQ14efd, discovered by the La Silla QUEST survey and followed by PESSTO. LSQ14efd was discovered few days after explosion and the observations cover up to ~100 days. The early photometric points show the signature of the cooling of the shock break-out event experienced by the progenitor at the time of the supernova explosion, one of the first for a type Ic supernova. A comparison with type Ic supernova spectra shows that LSQ14efd is quite similar to the type Ic SN 2004aw. These two supernovae have kinetic energies that are intermediate between standard Ic explosions and those which are the most energetic explosions known (e.g. SN 1998bw). We computed an analytical model for the light-curve peak and estimated the mass of the ejecta 6.3 +/- 0.5 Msun, a synthesized nickel mass of 0.25 Msun and a kinetic energy of Ekin = 5.6 +/- 0.5 x 10^51 erg. No connection between LSQ14efd and a GRB event could be established. However we point out that the supernova shows some spectroscopic similarities with the peculiar SN-Ia 1999ac and the SN-Iax SN 2008A. A core-collapse origin is most probable considering the spectroscopic, photometric evolution and the detection of the cooling of the shock break-out.


title: Monitoring survey of pulsating giant stars in the Local Group galaxies: survey description, science goals, target selection

abstract: The population of nearby dwarf galaxies in the Local Group constitutes a complete galactic environment, perfect suited for studying the connection between stellar populations and galaxy evolution. In this study, we are conducting an optical monitoring survey of the majority of dwarf galaxies in the Local Group, with the Isaac Newton Telescope (INT), to identify long period variable stars (LPVs). These stars are at the end points of their evolution and therefore their luminosity can be directly translated into their birth masses; this enables us to reconstruct the star formation history. By the end of the monitoring survey, we will have performed observations over ten epochs, spaced approximately three months apart, and identified long period, dust-producing AGB stars; five epochs of data have been obtained already. LPVs are also the main source of dust; in combination with Spitzer Space Telescope images at mid-IR wavelengths we will quantify the mass loss, and provide a detailed map of the mass feedback into the interstellar medium. We will also use the amplitudes in different optical passbands to determine the radius variations of the stars, and relate this to their mass loss.


title: Composite images of Galois for elliptic curves over $\mathbf{Q}$ & Entanglement fields

abstract: Let $E$ be an elliptic curve defined over $\mathbf{Q}$ without complex multiplication. For each prime $\ell$, there is a representation $\rho_{E,\ell}\colon \text{Gal}(\overline{\mathbf{Q}}/\mathbf{Q}) \to \text{GL}_2(\mathbf{F}_{\ell})$ that describes the Galois action on the $\ell$-torsion points of $E$. Building on recent work of Rouse--Zureick-Brown and Zywina, we find models for composite level modular curves whose rational points classify elliptic curves over $\mathbf{Q}$ with simultaneously non-surjective, composite image of Galois. We also provably determine the rational points on almost all of these curves. Finally, we give an application of our results to the study of entanglement fields.


title: Multilayer shallow water models with locally variable number of layers and semi-implicit time discretization

abstract: We propose an extension of the discretization approaches for multilayer shallow water models, aimed at making them more flexible and efficient for realistic applications to coastal flows. A novel discretization approach is proposed, in which the number of vertical layers and their distribution are allowed to change in different regions of the computational domain. Furthermore, semi-implicit schemes are employed for the time discretization, leading to a significant efficiency improvement for subcritical regimes. We show that, in the typical regimes in which the application of multilayer shallow water models is justified, the resulting discretization does not introduce any major spurious feature and allows again to reduce substantially the computational cost in areas with complex bathymetry. As an example of the potential of the proposed technique, an application to a sediment transport problem is presented, showing a remarkable improvement with respect to standard discretization approaches.


title: Large Deviation Theory for Parameter Estimation in Simple Neuron Models

abstract: To investigate the complex dynamics of a biological neuron that is subject to small random perturbations we can use stochastic neuron models. While many techniques have already been developed to study properties of such models, especially the analysis of the (expected) first-passage time or (E)FPT remains difficult. In this thesis I apply the large deviation theory (LDT), which is already well-established in physics and finance, to the problem of determining the EFPT of the mean-reverting Ornstein-Uhlenbeck (OU) process. The OU process instantiates the Stochastic Leaky Integrate and Fire model and thus serves as an example of a biologically inspired mathematical neuron model. I derive several classical results using much simpler mathematics than the original publications from neuroscience and I provide a few conceivable interpretations and perspectives on these derivations. Using these results I explore some possible applications for parameter estimation and I provide an additional mathematical justification for using a Poisson process as a small-noise approximation of the full model. Finally I perform several simulations to verify these results and to reveal systematic biases of this estimator.


title: Flow of Carreau Fluids Through Non-Uniform Pore throats

abstract: The creeping flow of a generalized Newtonian fluid in a non--uniform pore throat is investigated analytically. The analytical solution determines the flow regimes and the transition point from Newtonian to the power law flow regime. As an application of the new model, a regular lattice--based pore network model is constructed to simulate the flows of shear--thinning and shear--thickening fluids through porous media. The pore throats have convergent--divergent geometries.


title: Time-reversal breaking topological phase without Hall electric current in a two-dimensional Dirac semimetal protected by nonsymmorphic symmetry

abstract: We investigate the topological phase derived by time-reversal breaking fields in a nonsymmorphic symmetry-protected two-dimensional Dirac semimetal. When the nonsymmorphic symmetry is preserved even in the presence of the field, the two-dimensional electronic states change into two distinct topological phases with the insulating gap. One phase is well-known as quantum Hall states with chiral edge modes accompanying the Hall current, but the other one is an unconventional topological phase with helical edge modes in the absence of time-reversal symmetry.