title: Projection-Based Iterative Mode Scheduling for Switched Systems

abstract: This paper describes a method for scheduling the events of a switched system to achieve an optimal performance. The approach has guarantees on convergence and computational complexity that parallel derivative-based iterative optimization but in the infinite dimensional, integer constrained setting of mode scheduling. In comparison to methods relying on mixed integer programming, the presented approach does not require a priori discretizations of time or state. Furthermore, in comparison to embedding and relaxation methods, every iteration of the algorithm returns a dynamically feasible solution. A large class of problems call for optimal mode scheduling. This paper considers a vehicle tracking problem and a high dimensional multimachine power network synchronization problem. For the power network example, both single horizon and receding horizon approaches prevent instability of the network, and the receding horizon approach does so at near real-time speeds on a single processor.


title: Active correction of aperture discontinuities - optimized stroke minimization I: a new adaptive interaction matrix algorithm

abstract: Future searches for biomarkers on habitable exoplanets will rely on telescope instruments that achieve extremely high contrast at small planet-to-star angular separations. Coronagraphy is a promising starlight suppression technique, providing excellent contrast and throughput for off-axis sources on clear apertures. However, the complexity of space- and ground-based telescope apertures goes on increasing over time, owing to the combination of primary mirror segmentation, secondary mirror, and support structures. These discontinuities in the telescope aperture limit the coronagraph performance. In this paper, we present ACAD-OSM, a novel active method to correct for the diffractive effects of aperture discontinuities in the final image plane of a coronagraph. Active methods use one or several deformable mirrors that are controlled with an interaction matrix to correct for the aberrations in the pupil. However, they are often limited by the amount of aberrations introduced by aperture discontinuities. This algorithm relies on the recalibration of the interaction matrix during the correction process to overcome this limitation. We first describe the ACAD-OSM technique and compare it to the previous active methods for the correction of aperture discontinuities. We then show its performance in terms of contrast and off-axis throughput for static aperture discontinuities (segmentation, struts) and for some aberrations evolving over the life of the instrument (residual phase aberrations, artifacts in the aperture, misalignments in the coronagraph design). This technique can now obtain the earth-like planet detection threshold of 10^(-10) contrast on any given aperture over at least a 10% spectral bandwidth, with several coronagraph designs.


title: Quantum qubit switch: entropy and entanglement

abstract: The quantum entanglement is considered as one of the most important notions of quantum computing. The entanglement is a feature of quantum systems and it is used as a basis for many quantum algorithms and protocols. In this paper we analyze the level of entanglement for the quantum switch, during its work. The level of entanglement during the correct operating may be compared with the situation when a noise is present in the analyzed system. The noise changes the level of quantum entanglement and we utilize this observation to evaluate if the switch works properly. We present the formulae describing the value of entanglement with use of the Concurrence measure and also the value of entropy. Both mentioned values may be utilized to detect the noise during the operating of quantum switch.


title: Notes on Cofinality Spectrum Problems

abstract: These notes are based on Appalachian Set Theory lectures given by M. Malliaris on November 5, 2016 with D. Casey as the official scribe. The aim of the lectures was to present the setup and some key arguments of "Cofinality spectrum problems in model theory, set theory and general topology" by Malliaris and Shelah. This provides a sketch of the proof that $\mathfrak{p} = \mathfrak{t}$ and that $SOP_2$ theories are maximal in Keisler's Order.


title: Apparatus to control and visualize the impact of a high-energy laser pulse on a liquid target

abstract: We present an experimental apparatus to control and visualize the response of a liquid target to a laser-induced vaporization. We use a millimeter-sized drop as target and present two liquid-dye solutions that allow a variation of the absorption coefficient of the laser light in the drop by seven orders of magnitude. The excitation source is a Q-switched Nd:YAG laser at its frequency-doubled wavelength emitting nanosecond pulses with energy densities above the local vaporization threshold. The absorption of the laser energy leads to a large-scale liquid motion at timescales that are separated by several orders of magnitude, which we spatiotemporally resolve by a combination of ultra-high-speed and stroboscopic high-resolution imaging in two orthogonal views. Surprisingly, the large-scale liquid motion at upon laser impact is completely controlled by the spatial energy distribution obtained by a precise beam-shaping technique. The apparatus demonstrates the potential for accurate and quantitative studies of laser-matter interactions.


title: Large Scale Brownian Dynamics of Confined Suspensions of Rigid Particles

abstract: We introduce methods for large scale Brownian Dynamics (BD) simulation of many rigid particles of arbitrary shape suspended in a fluctuating fluid. Our method adds Brownian motion to the rigid multiblob method at a cost comparable to the cost of deterministic simulations. We demonstrate that we can efficiently generate deterministic and random displacements for many particles using preconditioned Krylov iterative methods, if kernel methods to efficiently compute the action of the Rotne-Prager-Yamakawa (RPY) mobility matrix and it "square" root are available for the given boundary conditions. We address a major challenge in large-scale BD simulations, capturing the stochastic drift term that arises because of the configuration-dependent mobility. Unlike the widely-used Fixman midpoint scheme, our methods utilize random finite differences and do not require the solution of resistance problems or the computation of the action of the inverse square root of the RPY mobility matrix. We construct two temporal schemes which are viable for large scale simulations, an Euler-Maruyama traction scheme and a Trapezoidal Slip scheme, which minimize the number of mobility solves per time step while capturing the required stochastic drift terms. We validate and compare these schemes numerically by modeling suspensions of boomerang shaped particles sedimented near a bottom wall. Using the trapezoidal scheme, we investigate the steady-state active motion in a dense suspensions of confined microrollers, whose height above the wall is set by a combination of thermal noise and active flows. We find the existence of two populations of active particles, slower ones closer to the bottom and faster ones above them, and demonstrate that our method provides quantitative accuracy even with relatively coarse resolutions of the particle geometry.


title: Factorization tests and algorithms arising from counting modular forms and automorphic representations

abstract: A theorem of Gekeler compares the number of non-isomorphic automorphic representations associated with the space of cusp forms of weight $k$ on $\Gamma_0(N)$ to a simpler function of $k$ and $N$, showing that the two are equal whenever $N$ is squarefree. We prove the converse of this theorem (with one small exception), thus providing a characterization of squarefree integers. We also establish a similar characterization of prime numbers in terms of the number of Hecke newforms of weight $k$ on $\Gamma_0(N)$. It follows that a hypothetical fast algorithm for computing the number of such automorphic representations for even a single weight $k$ would yield a fast test for whether $N$ is squarefree. We also show how to obtain bounds on the possible square divisors of a number $N$ that has been found to not be squarefree via this test, and we show how to probabilistically obtain the complete factorization of the squarefull part of $N$ from the number of such automorphic representations for two different weights. If in addition we have the number of such Hecke newforms for even a single weight $k$, then we show how to probabilistically factor $N$ entirely. All of these computations could be performed quickly in practice, given the number(s) of automorphic representations and modular forms as input.


title: A secular increase in continental crust nitrogen during the Precambrian

abstract: Recent work indicates the presence of substantial geologic nitrogen reservoirs in the mantle and continental crust. Importantly, this geologic nitrogen has exchanged between the atmosphere and the solid Earth over time. Changes in atmospheric nitrogen (i.e. atmospheric mass) have direct effects on climate and biological productivity. It is difficult to constrain, however, the evolution of the major nitrogen reservoirs through time. Here we show a secular increase in continental crust nitrogen through Earth history recorded in glacial tills (2.9 Ga to modern), which act as a proxy for average upper continental crust composition. Archean and earliest Palaeoproterozoic tills contain 66 $\pm$ 100 ppm nitrogen, whereas Neoproterozoic and Phanerozoic tills contain 290 $\pm$ 165 ppm nitrogen, whilst the isotopic composition has remained constant at ~4$\permil$. Nitrogen has accumulated in the continental crust through time, likely sequestered from the atmosphere via biological fixation. Our findings support dynamic, non-steady state behaviour of nitrogen through time, and are consistent with net transfer of atmospheric N to geologic reservoirs over time.


title: Efficient polarisation-preserving frequency conversion from a trapped-ion-compatible wavelength to the telecom C band

abstract: We demonstrate polarisation-preserving frequency conversion of single-photon-level light at 854 nm, resonant with a trapped-ion transition and qubit, to the 1550-nm telecom C band. A total photon in / fiber-coupled photon out efficiency of $\sim$ 30 % is achieved, for a free-running photon noise rate of $\sim$ 60 Hz. This performance would enable telecom conversion of 854-nm polarisation qubits, produced in existing trapped-ion systems, with a signal-to-noise ratio greater than 1. In combination with near-future trapped-ion systems, our converter would enable the observation of entanglement between an ion and a photon that has travelled more than 100 km in optical fiber: three orders of magnitude further than the state-of-the-art.


title: Automated Dyadic Data Recorder (ADDR) Framework and Analysis of Facial Cues in Deceptive Communication

abstract: We developed an online framework that can automatically pair two crowd-sourced participants, prompt them to follow a research protocol, and record their audio and video on a remote server. The framework comprises two web applications: an Automatic Quality Gatekeeper for ensuring only high quality crowd-sourced participants are recruited for the study, and a Session Controller which directs participants to play a research protocol, such as an interrogation game. This framework was used to run a research study for analyzing facial expressions during honest and deceptive communication using a novel interrogation protocol. The protocol gathers two sets of nonverbal facial cues in participants: features expressed during questions relating to the interrogation topic and features expressed during control questions. The framework and protocol were used to gather 151 dyadic conversations (1.3 million video frames). Interrogators who were lied to expressed the smile-related lip corner puller cue more often than interrogators who were being told the truth, suggesting that facial cues from interrogators may be useful in evaluating the honesty of witnesses in some contexts. Overall, these results demonstrate that this framework is capable of gathering high quality data which can identify statistically significant results in a communication study.


title: Two coexisting families of compact stars: observational implications for millisecond pulsars

abstract: It is usually thought that a single equation of state (EoS) model "correctly" represents cores of all compact stars. Here we emphasize that two families of compact stars, viz., neutron stars and strange stars, can coexist in nature, and that neutron stars can get converted to strange stars through the nucleation process of quark matter in the stellar center. From our fully general relativistic numerical computations of the structures of fast-spinning compact stars, known as millisecond pulsars, we find that such a stellar conversion causes a simultaneous spin-up and decrease in gravitational mass of these stars. This is a new type of millisecond pulsar evolution through a new mechanism, which gives rise to relatively lower mass compact stars with higher spin rates. This could have implication for the observed mass and spin distributions of millisecond pulsars. Such a stellar conversion can also rescue some massive, spin-supported millisecond pulsars from collapsing into black holes. Besides, we extend the concept of critical mass $M_{\rm cr}$ for the neutron star sequence (Berezhiani et al. 2003; Bombaci et al. 2004) to the case of fast-spinning neutron stars, and point out that neutron star EoS models cannot be ruled out by the stellar mass measurement alone. Finally, we emphasize the additional complexity for constraining EoS models, for example, by stellar radius measurements using X-ray observations, if two families of compact stars coexist.


title: A sharp bound for winning within a proportion of the maximum of a sequence

abstract: This note considers a variation of the full-information secretary problem where the random variables to be observed are independent and identically distributed. Consider $X_1,\dots,X_n$ to be an independent sequence of random variables, let $M_n:=\max\{X_1,\dots,X_n\}$, and the objective is to select the maximum of the sequence. What is the maximum probability of "stopping at the maximum"? That is, what is the stopping time $\tau$ adapted to $X_1,...,X_n$ that maximizes $P(X_{\tau}=M_n)$? This problem was examined by Gilbert and Mosteller \cite{GilMost} when in addition the common distribution is continuous. The optimal win probability in this case is denoted by $v_{n,max}^*$. What if it is desired to "stop within a proportion of the maximum"? That is, for $0<\alpha<1$, what is the stopping rule $\tau$ that maximizes $P(X_{\tau} \geq \alpha M_n)$? In this note both problems are treated as games, it is proven that for any continuous random variable $X$, if $\tau^*$ is the optimal stopping rule then $P(X_{\tau^*} \geq \alpha M_n)\geq v_{n,max}^*$, and this lower bound is sharp. Some examples and another interesting result are presented.


title: Assimilation of nearly turbulent Rayleigh-B\'enard flow through vorticity or local circulation measurements: a computational study

abstract: We introduce a continuous (downscaling) data assimilation algorithm for the 2D B\'enard convection problem using vorticity or local circulation measurements only. In this algorithm, a nudging term is added to the vorticity equation to constrain the model. Our numerical results indicate that the approximate solution of the algorithm is converging to the unknown reference solution (vorticity and temperature) corresponding to the measurements of the 2D B\'enard convection problem when only spatial coarse-grain measurements of vorticity are assimilated. Moreover, this convergence is realized using data which is much more coarse than the resolution needed to satisfy rigorous analytical estimates.


title: How Does Knowledge of the AUC Constrain the Set of Possible Ground-truth Labelings?

abstract: Recent work on privacy-preserving machine learning has considered how data-mining competitions such as Kaggle could potentially be "hacked", either intentionally or inadvertently, by using information from an oracle that reports a classifier's accuracy on the test set. For binary classification tasks in particular, one of the most common accuracy metrics is the Area Under the ROC Curve (AUC), and in this paper we explore the mathematical structure of how the AUC is computed from an n-vector of real-valued "guesses" with respect to the ground-truth labels. We show how knowledge of a classifier's AUC on the test set can constrain the set of possible ground-truth labelings, and we derive an algorithm both to compute the exact number of such labelings and to enumerate efficiently over them. Finally, we provide empirical evidence that, surprisingly, the number of compatible labelings can actually decrease as n grows, until a test set-dependent threshold is reached.


title: Plasmon-Driven Acceleration in a Photo-Excited Nanotube

abstract: A plasmon-assisted channeling acceleration can be realized with a large channel, possibly at the nanometer scale. Carbon nanotubes (CNTs) are the most typical example of nano-channels that can confine a large number of channeled particles in a photon-plasmon coupling condition. This paper presents a theoretical and numerical study on the concept of high-field charge acceleration driven by photo-excited Luttinger-liquid plasmons (LLP) in a nanotube. An analytic description of the plasmon-assisted laser acceleration is detailed with practical acceleration parameters, in particular with specifications of a typical tabletop femtosecond laser system. The maximally achievable acceleration gradients and energy gains within dephasing lengths and CNT lengths are discussed with respect to laser-incident angles and CNT-filling ratios.


title: Relatively Hyperbolic Groups have Semistabile Fundamental Group at Infinity

abstract: Suppose $G$ is a 1-ended finitely generated group that is hyperbolic relative to P a finite collection of 1-ended finitely generated subgroups. Our main theorem states that if the boundary $\partial (G, P)$ has no cut point, then $G$ has semistable fundamental group at $\infty$. Under mild conditions on $G$ and the members of P the 1-ended hypotheses and the no cut point condition can be eliminated to obtain the same semistability conclusion. We give an example that shows our main result is somewhat optimal. Finally, we improve a "double dagge" result of F. Dahmani and D. Groves.


title: A classifier for gravitational-wave inspiral signals in non-ideal single-detector data

abstract: We describe a multivariate classifier for candidate events in a templated search for gravitational-wave (GW) inspiral signals from neutron-star--black-hole (NS-BH) binaries, in data from ground-based detectors where sensitivity is limited by non-Gaussian noise transients. The standard signal-to-noise ratio (SNR) and chi-squared test for inspiral searches use only properties of a single matched filter at the time of an event; instead, we propose a classifier using features derived from a bank of inspiral templates around the time of each event, and also from a search using approximate sine-Gaussian templates. The classifier thus extracts additional information from strain data to discriminate inspiral signals from noise transients. We evaluate a Random Forest classifier on a set of single-detector events obtained from realistic simulated advanced LIGO data, using simulated NS-BH signals added to the data. The new classifier detects a factor of 1.5 -- 2 more signals at low false positive rates as compared to the standard 're-weighted SNR' statistic, and does not require the chi-squared test to be computed. Conversely, if only the SNR and chi-squared values of single-detector events are available, Random Forest classification performs nearly identically to the re-weighted SNR.


title: Ageing in Mortal Superdiffusive L\'evy Walkers

abstract: A growing body of literature examines the effects of superdiffusive subballistic movement pre-measurement (ageing or time lag) on observations arising from single-particle tracking. A neglected aspect is the finite lifetime of these L\'{e}vy walkers, be they proteins, cells or larger structures. We examine the effects of ageing on the motility of mortal walkers, and discuss the means by which permanent stopping of walkers may be categorised as arising from `natural' death or experimental artefacts such as low photostability or radiation damage. This is done by comparison of the walkers' mean squared displacement (MSD) with the front velocity of propagation of a group of walkers, which is found to be invariant under time lags. For any running time distribution of a mortal random walker, the MSD is tempered by the stopping rate $\theta$. This provides a physical interpretation for truncated heavy-tailed diffusion processes and serves as a tool by which to better classify the underlying running time distributions of random walkers. Tempering of aged MSDs raises the issue of misinterpreting superdiffusive motion which appears Brownian or subdiffusive over certain time scales.


title: Characterizing bad semidefinite programs: normal forms and short proofs

abstract: Semidefinite programs (SDPs) -- some of the most useful and versatile optimization problems of the last few decades -- are often pathological: the optimal values of the primal and dual problems may differ and may not be attained. Such SDPs are both theoretically interesting and often impossible to solve; yet, the pathological SDPs in the literature look strikingly similar. Based on our recent work \cite{Pataki:17} we characterize pathological semidefinite systems by certain {\em excluded matrices}, which are easy to spot in all published examples. Our main tool is a normal (canonical) form of semidefinite systems, which makes their pathological behavior easy to verify. The normal form is constructed in a surprisingly simple fashion, using mostly elementary row operations inherited from Gaussian elimination. The proofs are elementary and can be followed by a reader at the advanced undergraduate level. As a byproduct, we show how to transform any linear map acting on symmetric matrices into a normal form, which allows us to quickly check whether the image of the semidefinite cone under the map is closed. We can thus introduce readers to a fundamental issue in convex analysis: the linear image of a closed convex set may not be closed, and often simple conditions are available to verify the closedness, or lack of it.


title: Duality of Floating and Illumination Bodies

abstract: We investigate a duality relation between floating and illumination bodies. The definitions of these two bodies suggest that the polar of the floating body should be similar to the illumination body of the polar. We consider this question for the class of centrally symmetric convex bodies. We provide precise estimates for $(B_p^n)$ and for centrally symmetric convex bodies with everywhere positive Gauss curvature. Our estimates show that equality of the polar of the floating body and the illumination body of the polar can only be achieved in the case of ellipsoids.


title: A Stochastic Approach to Shortcut Bridging in Programmable Matter

abstract: In a self-organizing particle system, an abstraction of programmable matter, simple computational elements called particles with limited memory and communication self-organize to solve system-wide problems of movement, coordination, and configuration. In this paper, we consider a stochastic, distributed, local, asynchronous algorithm for "shortcut bridging", in which particles self-assemble bridges over gaps that simultaneously balance minimizing the length and cost of the bridge. Army ants of the genus Eciton have been observed exhibiting a similar behavior in their foraging trails, dynamically adjusting their bridges to satisfy an efficiency trade-off using local interactions. Using techniques from Markov chain analysis, we rigorously analyze our algorithm, show it achieves a near-optimal balance between the competing factors of path length and bridge cost, and prove that it exhibits a dependence on the angle of the gap being "shortcut" similar to that of the ant bridges. We also present simulation results that qualitatively compare our algorithm with the army ant bridging behavior. Our work gives a plausible explanation of how convergence to globally optimal configurations can be achieved via local interactions by simple organisms (e.g., ants) with some limited computational power and access to random bits. The proposed algorithm also demonstrates the robustness of the stochastic approach to algorithms for programmable matter, as it is a surprisingly simple extension of our previous stochastic algorithm for compression.


title: Intelligent Disaster Response via Social Media Analysis - A Survey

abstract: The success of a disaster relief and response process is largely dependent on timely and accurate information regarding the status of the disaster, the surrounding environment, and the affected people. This information is primarily provided by first responders on-site and can be enhanced by the firsthand reports posted in real-time on social media. Many tools and methods have been developed to automate disaster relief by extracting, analyzing, and visualizing actionable information from social media. However, these methods are not well integrated in the relief and response processes and the relation between the two requires exposition for further advancement. In this survey, we review the new frontier of intelligent disaster relief and response using social media, show stages of disasters which are reflected on social media, establish a connection between proposed methods based on social media and relief efforts by first responders, and outline pressing challenges and future research directions.


title: Status Updates Through Multicast Networks

abstract: Using age of information as the freshness metric, we examine a multicast network in which real-time status updates are generated by the source and sent to a group of $n$ interested receivers. We show that in order to keep the information freshness at each receiver, the source should terminate the transmission of the current update and start sending a new update packet as soon as it receives the acknowledgements back from any $k$ out of $n$ nodes. As the source stopping threshold $k$ increases, a node is more likely to get the latest generated update, but the age of the most recent update is more likely to become outdated. We derive the age minimized stopping threshold $k$ that balances the likelihood of getting the latest update and the freshness of the latest update for shifted exponential link delay. Through numerical evaluations for different stopping strategies, we find that waiting for the acknowledgements from the earliest $k$ out of $n$ nodes leads to lower average age than waiting for a pre-selected group of $k$ nodes. We also observe that a properly chosen threshold $k$ can prevent information staleness for increasing number of nodes $n$ in the multicast network.


title: Theoretical investigations of an information geometric approach to complexity

abstract: It is known that statistical model selection as well as identification of dynamical equations from available data are both very challenging tasks. Physical systems behave according to their underlying dynamical equations which, in turn, can be identified from experimental data. Explaining data requires selecting mathematical models that best capture the data regularities. The existence of fundamental links among physical systems, dynamical equations, experimental data and statistical modeling motivate us to present in this article our theoretical modeling scheme which combines information geometry and inductive inference methods to provide a probabilistic description of complex systems in the presence of limited information. Special focus is devoted to describe the role of our entropic information geometric complexity measure. In particular, we provide several illustrative examples wherein our modeling scheme is used to infer macroscopic predictions when only partial knowledge of the microscopic nature of a given system is available. Finally, limitations, possible improvements, and future investigations are discussed.


title: Floating and Illumination Bodies for Polytopes: Duality Results

abstract: We consider the question how well a floating body can be approximated by the polar of the illumination body of the polar. We establish precise convergence results in the case of centrally symmetric polytopes. This leads to a new affine invariant which is related to the cone measure of the polytope.


title: A note about words which coincide except in one position

abstract: In this short note, we show that a result about words which coincide except in one position given as an exercise in Lothaire's Algebraic Combinatorics on Words is false. Moreover, we derive a modified statement which allows us to fix the proof of a theorem which originally used the result of this exercise.


title: Genericity of Infinite Entropy for Maps with Low Regularity

abstract: For bi-Lipschitz homeomorphisms of a compact manifold it is known that topological entropy is always finite. For compact manifolds of dimension two or greater, we show that in the closure of the space of bi-Lipschitz homeomorphisms, with respect to either the H\"older or the Sobolev topologies, topological entropy is generically infinite. We also prove versions of the $C^1$-Closing Lemma in either of these spaces. Finally, we give examples of homeomorphisms with infinite topological entropy which are H\"older and/or Sobolev of every exponent.


title: Deep Learning the Physics of Transport Phenomena

abstract: We have developed a new data-driven paradigm for the rapid inference, modeling and simulation of the physics of transport phenomena by deep learning. Using conditional generative adversarial networks (cGAN), we train models for the direct generation of solutions to steady state heat conduction and incompressible fluid flow purely on observation without knowledge of the underlying governing equations. Rather than using iterative numerical methods to approximate the solution of the constitutive equations, cGANs learn to directly generate the solutions to these phenomena, given arbitrary boundary conditions and domain, with high test accuracy (MAE$<$1\%) and state-of-the-art computational performance. The cGAN framework can be used to learn causal models directly from experimental observations where the underlying physical model is complex or unknown.


title: Addendum to: Edge-Unfolding Nearly Flat Convex Caps

abstract: This addendum to [O'R17] establishes that a nearly flat acutely triangulated convex cap in the sense of that paper can be edge-unfolded even if closed to a polyhedron by adding the convex polygonal base under the cap.


title: New limits on cosmic strings from gravitational wave observation

abstract: We combine new analysis of the stochastic gravitational wave background to be expected from cosmic strings with the latest pulsar timing array (PTA) limits to give an upper bound on the energy scale of the possible cosmic string network, $G\mu < 1.5\times 10^{-11}$ at the 95% confidence level. We also show bounds from LIGO and to be expected from LISA and BBO. Current estimates for the gravitational wave background from supermassive black hole binaries are at the level where a PTA detection is expected. But if PTAs do observe a background soon, it will be difficult in the short term to distinguish black holes from cosmic strings as the source, because the spectral indices from the two sources happen to be quite similar. If PTAs do not observe a background, then the limits on $G\mu$ will improve somewhat, but a string network with $G\mu$ substantially below $10^{-11}$ will produce gravitational waves primarily at frequencies too high for PTA observation, so significant further progress will depend on intermediate-frequency observatories such as LISA, DECIGO and BBO.


title: An Analysis of ISO 26262: Using Machine Learning Safely in Automotive Software

abstract: Machine learning (ML) plays an ever-increasing role in advanced automotive functionality for driver assistance and autonomous operation; however, its adequacy from the perspective of safety certification remains controversial. In this paper, we analyze the impacts that the use of ML as an implementation approach has on ISO 26262 safety lifecycle and ask what could be done to address them. We then provide a set of recommendations on how to adapt the standard to accommodate ML.


title: Distant clusters of galaxies in the 2XMM/SDSS footprint: follow-up observations with the LBT

abstract: Context: Galaxy clusters at high redshift are important to test cosmological models and models for the growth of structure. They are difficult to find in wide-angle optical surveys, however, leaving dedicated follow-up of X-ray selected candidates as one promising identification route. Aims: We aim to increase the number of galaxy clusters beyond the SDSS-limit, z ~ 0.75. Methods: We compiled a list of extended X-ray sources from the 2XMMp catalogue within the footprint of the Sloan Digital Sky Survey. Fields without optical counterpart were selected for further investigation. Deep optical imaging and follow-up spectroscopy were obtained with the Large Binocular Telescope, Arizona (LBT), of those candidates not known to the literature. Results: From initially 19 candidates, selected by visually screening X-ray images of 478 XMM-Newton observations and the corresponding SDSS images, 6 clusters were found in the literature. Imaging data through r,z filters were obtained for the remaining candidates, and 7 were chosen for multi-object (MOS) spectroscopy. Spectroscopic redshifts, optical magnitudes, and X-ray parameters (flux, temperature, and luminosity) are presented for the clusters with spectroscopic redshifts. The distant clusters studied here constitute one additional redshift bin for studies of the L-T relation, which does not seem to evolve from high to low redshifts. ...


title: 3D angle-of-arrival positioning using von Mises-Fisher distribution

abstract: We propose modeling an angle-of-arrival (AOA) positioning measurement as a von Mises-Fisher (VMF) distributed unit vector instead of the conventional normally distributed azimuth and elevation measurements. Describing the 2-dimensional AOA measurement with three numbers removes discontinuities and reduces nonlinearity at the poles of the azimuth-elevation coordinate system. Our computer simulations show that the proposed VMF measurement noise model based filters outperform the normal distribution based algorithms in accuracy in a scenario where close-to-pole measurements occur frequently.


title: Millimeter Wave Multi-Beam-Switching Antenna

abstract: The paper introduces a high gain end-fire printed antenna operating at 14 GHz with a bandwidth of 2 GHz. The 3x3 antenna sub-array formed of this element is used to present its beam steering and multi-beam-switching capabilities for millimeter wave Massive MIMO antenna arrays of 5G networks. The co-existence of these two techniques in a Massive MIMO sub-array, aided by the use of a simple two vector excitation matrix technique, is shown in this paper. The high gain and small form factor of the antenna is useful in applications, like 5G arrays, where space is premium and the inherent layered structure also aids to the better cooling of the antenna front end.


title: On the statistical distribution of prime numbers, a view from where the distribution of prime numbers is not erratic

abstract: Currently there is no known efficient formula for primes. Besides that, prime numbers have great importance in e.g., information technology such as public-key cryptography, and their position and possible or impossible functional generation among the natural numbers is an ancient dilemma. The properties of the functions 2ab+a+b in the domain of natural numbers are introduced, analyzed, and exhibited to illustrate how these single out all the prime numbers from the full set of odd numbers. The characterization of odd primes vs. odd non-primes can be done with 2ab+a+b among the odd natural numbers as an analogue to the other, well known type of fundamental characterization for irrational and rational numbers among the real numbers. The prime number theorem, twin primes and erratic nature of primes, are also commented upon with respect to selection, as well as with the Fermat and Euler numbers as examples.


title: Sequential Composition in the Presence of Intermediate Termination (Extended Abstract)

abstract: The standard operational semantics of the sequential composition operator gives rise to unbounded branching and forgetfulness when transparent process expressions are put in sequence. Due to transparency, the correspondence between context-free and pushdown processes fails modulo bisimilarity, and it is not clear how to specify an always terminating half counter. We propose a revised operational semantics for the sequential composition operator in the context of intermediate termination. With the revised operational semantics, we eliminate transparency, allowing us to establish a close correspondence between context-free processes and pushdown processes. Moreover, we prove the reactive Turing powerfulness of TCP with iteration and nesting with the revised operational semantics for sequential composition.


title: Encoding Multifunctional Nanostructured Metasurfaces

abstract: Metasurfaces offer unconventional control of light to shape optical wavefronts within two-dimensional nanoscale structures. A nanostructured metallic thin film can establish an exclusive relationship between its surface structure and optical spectrum. The unique structural features, continuous metal coverage and smooth surface profile, can help obtain highly coherent plasmonic resonances that strongly depend on the structural parameters of nanostructured surfaces. It allows us to encode optical responses and information into a structural format on a variety of materials. We show experimentally that the nanostructured metasurfaces can be encoded to detect and display chemical and biomolecular interactions, show colorful graphic images, and create multichannel holograms.


title: Efficient Point-Counting Algorithms for Superelliptic Curves

abstract: In this paper, we present efficient algorithms for computing the number of points and the order of the Jacobian group of a superelliptic curve over finite fields of prime order p. Our method employs the Hasse-Weil bounds in conjunction with the Hasse-Witt matrix for superelliptic curves, whose entries we express in terms of multinomial coefficients. We present a fast algorithm for counting points on specific trinomial superelliptic curves and a slower, more general method for all superelliptic curves. For the first case, we reduce the problem of simplifying the entries of the Hasse-Witt matrix modulo p to a problem of solving quadratic Diophantine equations. For the second case, we extend Bostan et al.'s method for hyperelliptic curves to general superelliptic curves. We believe the methods we describe are asymptotically the most efficient known point-counting algorithms for certain families of trinomial superelliptic curves.


title: Fractional cable equation for general geometry, a model of axons with swellings and anomalous diffusion

abstract: Different experimental studies have reported anomalous diffusion in brain tissues and notably this anomalous diffusion is expressed through fractional derivatives. Axons are important to understand neurodegenerative diseases such as multiple sclerosis, Alzheimer's disease and Parkinson's disease. Indeed, abnormal accumulation of proteins and organelles in axons is a hallmark feature of these diseases. The diffusion in the axons can become to anomalous as a result from this abnormality. In this case the voltage propagation in axons is affected. Another hallmark feature of different neurodegenerative diseases is given by discrete swellings along the axon. In order to model the voltage propagation in axons with anomalous diffusion and swellings, in this paper we propose a fractional cable equation for general geometry. This generalized equation depends on fractional parameters and geometric quantities such as the curvature and torsion of the cable. For a cable with a constant radius we show that the voltage decreases when the fractional effect increases. In cables with swellings we find that when the fractional effect or the swelling radius increase, the voltage decreases. A similar behavior is obtained when the number of swellings and the fractional effect increase. Moreover, we find that when the radius swelling (or the number of swellings) and the fractional effect increase at the same time, the voltage dramatically decreases.


title: Muon Pair Production in Radiative Decays of Heavy Quarkonia

abstract: Muon pair production in exclusive radiative decays $\chi_{cJ}\to J/\psi \mu\mu$, $\chi_{bJ}\to\Upsilon\mu\mu$ are considered. We present numerical values of the branching fractions and distributions over various kinematical variables.


title: Large Teams Have Developed Science and Technology; Small Teams Have Disrupted It

abstract: Teams dominate the production of high-impact science and technology. Analyzing teamwork from more than 50 million papers, patents, and software products, 1954-2014, we demonstrate across this period that larger teams developed recent, popular ideas, while small teams disrupted the system by drawing on older and less prevalent ideas. Attention to work from large teams came immediately, while advances by small teams succeeded further into the future. Differences between small and large teams magnify with impact - small teams have become known for disruptive work and large teams for developing work. Differences in topic and re- search design account for part of the relationship between team size and disruption, but most of the effect occurs within people, controlling for detailed subject and article type. Our findings suggest the importance of supporting both small and large teams for the sustainable vitality of science and technology.


title: Radiation driven instability of rapidly rotating relativistic stars: criterion and evolution equations via multipolar expansion of gravitational waves

abstract: I suggest a novel approach for deriving evolution equations for rapidly rotating relativistic stars affected by radiation-driven Chandrasekhar-Friedman-Schutz (CFS) instability. This approach is based on the multipolar expansion of gravitational wave emission and appeals to the global physical properties of the star (energy, angular momentum, and thermal state), but not to canonical energy and angular momentum, which is traditional. It leads to simple derivation of the CFS instability criterion for normal modes and the evolution equations for a star, affected by this instability. The approach also gives a precise form to simple explanation of the CFS instability: it occurs when two conditions met: (a) gravitational wave emission removes angular momentum from the rotating star (thus releasing the rotation energy) and (b) gravitational waves carry less energy, than the released amount of the rotation energy. Results are illustrated by r-mode instability in slowly rotating Newtonian stars. In resulting evolution equations emission of gravitational waves directly affects the spin frequency, being in line with the arguments by Levin & Ushomirsky (2001), but in contrast to widely accepted equations by Lindblom et al. (1998); Ho & Lai (2000), there effective spin frequency decrease is coupled with dissipation of unstable mode. This problem, initially stressed by Levin & Ushomirsky (2001), is shown to be superficial, and arises as a result of definition of the effective spin frequency by Lindblom et al. (1998); Ho & Lai (2000). Namely, it is shown, that if this definition is taken into account properly, the evolution equations coincide with obtained here in the leading order in mode amplitude. I also argue that the next-to-leading order terms in evolution equations (which differ for Owen et al. 1998 and Ho & Lai 2000) require clarification and thus it would be more self-consistent to omit them.


title: Spin-Valley Coherent Phases of the $\nu=0$ Quantum Hall State in Bilayer Graphene

abstract: Bilayer graphene (BLG) offers a rich platform for broken symmetry states stabilized by interactions. In this work we study the phase diagram of BLG in the quantum Hall regime at filling factor $\nu=0$ within the Hartree-Fock approximation. In the simplest non-interacting situation this system has eight (nearly) degenerate Landau levels near the Fermi energy, characterized by spin, valley, and orbital quantum numbers. We incorporate in our study two effects not previously considered: (i) the nonperturbative effect of trigonal warping in the single-particle Hamiltonian, and (ii) short-range SU(4) symmetry-breaking interactions that distinguish the energetics of the orbitals. We find within this model a rich set of phases, including ferromagnetic, layer-polarized, canted antiferromagnetic, Kekul\'e, a "spin-valley entangled" state, and a "broken U(1) $\times$ U(1)" phase. This last state involves independent spontaneous symmetry breaking in the layer and valley degrees of freedom, and has not been previously identified. We present phase diagrams as a function of interlayer bias $D$ and perpendicular magnetic field $B_{\perp}$ for various interaction and Zeeman couplings, and discuss which are likely to be relevant to BLG in recent measurements. Experimental properties of the various phases and transitions among them are also discussed.


title: Network Vector: Distributed Representations of Networks with Global Context

abstract: We propose a neural embedding algorithm called Network Vector, which learns distributed representations of nodes and the entire networks simultaneously. By embedding networks in a low-dimensional space, the algorithm allows us to compare networks in terms of structural similarity and to solve outstanding predictive problems. Unlike alternative approaches that focus on node level features, we learn a continuous global vector that captures each node's global context by maximizing the predictive likelihood of random walk paths in the network. Our algorithm is scalable to real world graphs with many nodes. We evaluate our algorithm on datasets from diverse domains, and compare it with state-of-the-art techniques in node classification, role discovery and concept analogy tasks. The empirical results show the effectiveness and the efficiency of our algorithm.


title: Hamiltonian approach to slip-stacking dynamics

abstract: Hamiltonian dynamics has been applied to study the slip-stacking dynamics. The canonical-perturbation method is employed to obtain the second-harmonic correction term in the slip-stacking Hamiltonian. The Hamiltonian approach provides a clear optimal method for choosing the slip-stacking parameter and improving stacking efficiency. The dynamics are applied specifically to the Fermilab Booster-Recycler complex. The dynamics can also be applied to other accelerator complexes.


title: The time-dependent Schr\"odinger equation with piecewise constant potentials

abstract: The linear Schr\"odinger equation with piecewise constant potential in one spatial dimension is a well-studied textbook problem. It is one of only a few solvable models in quantum mechanics and shares many qualitative features with physically important models. In examples such as "particle in a box" and tunneling, attention is restricted to the time-independent Schr\"odinger equation. This paper combines the Unified Transform Method and recent insights for interface problems to present fully explicit solutions for the time-dependent problem.


title: The stability index for dynamically defined Weierstrass functions

abstract: Let $\hat{T} : X \times \mathbb{R} \to X \times \mathbb{R}$ given by $\hat{T}(x,t) = (Tx, g_x(t))$ be a skew-product dynamical system where $T : X \to X$ is a mixing conformal expanding map and, for each $x \in X$, $g_x : \mathbb{R} \to \mathbb{R}$ is an affine map of the form $g_x(t)=-f(x)+\lambda(x)^{-1}t$. Under a suitable contraction hypotheses on $\lambda$ there exists a measurable function $u: X \to \mathbb{R}$ such that $\text{graph}(u) = \{(x,u(x)) \mid x\in X\}$ is $\hat{T}$-invariant and divides $X \times \mathbb{R}$ into two regions, $\mathbb{B}^+$ and $\mathbb{B}^-$, consisting of points that are repelled under iteration by $\hat{T}$ to $\pm\infty$. These two regions act as basins of attraction to $\pm \infty$ in the sense of Milnor. The two basins have a complicated local structure: a neighbourhood of a point $(x,t) \in ^+$ will typically intersect $\mathbb{B}^-$ in a set of positive measure. The stability index (as introduced by Podvigina and Ashwin \cite{podviginaashwin:11} for general Milnor attractors) is the rate of polynomial decay of the measure of this intersection. We calculate the stability index at typical points in $X \times \mathbb{R}$. We also perform a multifractal analysis of the level sets of the stability index.


title: The Muffin Problem

abstract: You have $m$ muffins and $s$ students. You want to divide the muffins into pieces and give the shares to students such that every student has $\frac{m}{s}$ muffins. Find a divide-and-distribute protocol that maximizes the minimum piece. Let $f(m,s)$ be the minimum piece in the optimal protocol. We prove that $f(m,s)$ exists, is rational, and finding it is computable (though possibly difficult). We show that $f(m,s)$ can be derived from $f(s,m)$; hence we need only consider $m\ge s$. We give a function $FC(m,s)$ such that, for $m\ge s+1$, $f(m,s)\le FC(m,s)$. It is often the case that $f(m,s)=FC(m,s)$. More formally, for all $s$, for all but a finite number of $m$, $f(m,s)=FC(m,s)$. This leads to a nice formula for $f(m,s)$, though there are exceptions to it. We give a formula $INT(m,s)$, which has 6 parts, such that for many of the exceptional $m$, $f(m,s)=INT(m,s)<FC(m,s)$. This works for most of the exceptional $m$ where ceil${2m/s}\ge 4$. There are still some exceptional $m$ with ceil${2m/s}=3$ (if its $\le 2$ then the problem is trivial). For these cases we have a way to {\it generate theorems}. For $1\le d\le 7$ we have generated formulas for $f(s+d,s)$. We do not have a theorem here but we do have a methodology which leads to, for some of the $m$, a value $BM(m,s)$ such that often $f(m,s)=BM(m,s)<INT(m,s)<FC(m,s)$. So far it seems like, for $m\ge s$, $f(m,s) = \min\{FC(m,s), INT(m,s), BM(m,s) \}$, though we have not prove this. For $1\le s\le 50$ and $s\le m\le 60$ we have obtained $f(m,s)$ for all but 20 values.


title: Magnetic Control of Transmission and Helicity of Nano-Structured Optical Beams in Magnetoplasmonic Vortex Lenses

abstract: We theoretically investigate the generation of far-field propagating optical beams with a desired orbital angular momentum by using an archetypical magnetoplasmonic tip surrounded by a gold spiral slit. The use of a magnetic material can lead to important implications once magneto-optical activity is activated through the application of an external magnetic field. The physical model and the numerical study presented here introduce the concept of magnetically tunable plasmonic vortex lens, namely a magnetoplasmonic vortex lens, which ensures a tunable selectivity in the polarization state of the generated nanostructured beam. The presented system provides a promising platform for a localized excitation of plasmonic vortices followed by their beaming in the far-field with an active modulation of both light's transmittance and helicity.


title: The Direct Effect of Toroidal Magnetic Fields on Stellar Oscillations: An Analytical Expression for the General Matrix Element

abstract: Where is the solar dynamo located and what is its modus operandi? These are still open questions in solar physics. Helio- and asteroseismology can help answer them by enabling us to study solar and stellar internal structures through global oscillations. The properties of solar and stellar acoustic modes are changing with the level of magnetic activity. However, until now, the inference on subsurface magnetic fields with seismic measures has been very limited. The aim of this paper is to develop a formalism to calculate the effect of large-scale toroidal magnetic fields on solar and stellar global oscillation eigenfunctions and eigenfrequencies. If the Lorentz force is added to the equilibrium equation of motion, stellar eigenmodes can couple. In quasi-degenerate perturbation theory, this coupling, also known as the direct effect, can be quantified by the general matrix element. We present the analytical expression of the matrix element for a superposition of subsurface zonal toroidal magnetic field configurations. The matrix element is important for forward calculations of perturbed solar and stellar eigenfunctions and frequency perturbations. The results presented here will help to ascertain solar and stellar large-scale subsurface magnetic fields, and their geometric configuration, strength, and their change over the course of activity cycles.


title: A lower bound for the principal eigenvalue of fully nonlinear elliptic operators

abstract: In this article we present a new technique to obtain a lower bound for the principal Dirichlet eigenvalue of a fully nonlinear elliptic operator. We ilustrate the construction of an appropriate radial function required to obtain the bound in several examples. In particular we use our results to prove that $\lim_{p\to \infty}\lambda_{1,p}=\lambda_{1,\infty}=\left(\frac{\pi}{2R}\right)^2$ where $\lambda_{1,p}$ and $\lambda_{1,\infty}$ are the principal eigenvalue for the homogeneous $p$-laplacian and the homogeneous infinity laplacian respectively.


title: Attack-Aware Multi-Sensor Integration Algorithm for Autonomous Vehicle Navigation Systems

abstract: In this paper, we propose a fault detection and isolation based attack-aware multi-sensor integration algorithm for the detection of cyberattacks in autonomous vehicle navigation systems. The proposed algorithm uses an extended Kalman filter to construct robust residuals in the presence of noise, and then uses a parametric statistical tool to identify cyberattacks. The parametric statistical tool is based on the residuals constructed by the measurement history rather than one measurement at a time in the properties of discrete-time signals and dynamic systems. This approach allows the proposed multi-sensor integration algorithm to provide quick detection and low false alarm rates for applications in dynamic systems. An example of INS/GNSS integration of autonomous navigation systems is presented to validate the proposed algorithm by using a software-in-the-loop simulation.


title: Reservoir of Diverse Adaptive Learners and Stacking Fast Hoeffding Drift Detection Methods for Evolving Data Streams

abstract: The last decade has seen a surge of interest in adaptive learning algorithms for data stream classification, with applications ranging from predicting ozone level peaks, learning stock market indicators, to detecting computer security violations. In addition, a number of methods have been developed to detect concept drifts in these streams. Consider a scenario where we have a number of classifiers with diverse learning styles and different drift detectors. Intuitively, the current 'best' (classifier, detector) pair is application dependent and may change as a result of the stream evolution. Our research builds on this observation. We introduce the $\mbox{Tornado}$ framework that implements a reservoir of diverse classifiers, together with a variety of drift detection algorithms. In our framework, all (classifier, detector) pairs proceed, in parallel, to construct models against the evolving data streams. At any point in time, we select the pair which currently yields the best performance. We further incorporate two novel stacking-based drift detection methods, namely the $\mbox{FHDDMS}$ and $\mbox{FHDDMS}_{add}$ approaches. The experimental evaluation confirms that the current 'best' (classifier, detector) pair is not only heavily dependent on the characteristics of the stream, but also that this selection evolves as the stream flows. Further, our $\mbox{FHDDMS}$ variants detect concept drifts accurately in a timely fashion while outperforming the state-of-the-art.


title: End-to-end Face Detection and Cast Grouping in Movies Using Erd\H{o}s-R\'{e}nyi Clustering

abstract: We present an end-to-end system for detecting and clustering faces by identity in full-length movies. Unlike works that start with a predefined set of detected faces, we consider the end-to-end problem of detection and clustering together. We make three separate contributions. First, we combine a state-of-the-art face detector with a generic tracker to extract high quality face tracklets. We then introduce a novel clustering method, motivated by the classic graph theory results of Erd\H{o}s and R\'enyi. It is based on the observations that large clusters can be fully connected by joining just a small fraction of their point pairs, while just a single connection between two different people can lead to poor clustering results. This suggests clustering using a verification system with very few false positives but perhaps moderate recall. We introduce a novel verification method, rank-1 counts verification, that has this property, and use it in a link-based clustering scheme. Finally, we define a novel end-to-end detection and clustering evaluation metric allowing us to assess the accuracy of the entire end-to-end system. We present state-of-the-art results on multiple video data sets and also on standard face databases.


title: Demonstration of broadband contrast at 1.2 $\lambda$/D and greater for the EXCEDE Starlight Suppression System

abstract: The EXoplanetary Circumstellar Environments and Disk Explorer (EXCEDE) science mission concept uses a visible-wavelength Phase-Induced Amplitude Apodization (PIAA) coronagraph to enable high-contrast imaging of circumstellar debris systems and some giant planets at angular separations reaching into the habitable zones of some of the nearest stars. We report on the experimental results obtained in the vacuum chamber at the Lockheed Martin Advanced Technology Center in 10% broadband light centered about 650 nm, with a median contrast of $1 \times 10^{-5}$ between 1.2 and 2.0 $\lambda$/D simultaneously with $3 \times 10^{-7}$ contrast between 2 and 11 $\lambda/D$ for a single-sided dark hole using a deformable mirror (DM) upstream of the PIAA coronagraph. These results are stable and repeatable as demonstrated by three measurements runs with DM settings set from scratch and maintained on the best 90% out of the 1000 collected frames. We compare the reduced experimental data with simulation results from modeling observed experimental limits; observed performance is consistent with uncorrected low-order modes not estimated by the Low Order Wavefront Sensor (LOWFS). Modeled sensitivity to bandwidth and residual tip/tilt modes is well-matched to the experiment.


title: Spontaneous avalanche dephasing in large Rydberg ensembles

abstract: Strong dipole-exchange interactions due to spontaneously produced contaminant states can trigger rapid dephasing in many-body Rydberg ensembles [E. Goldschmidt et al., PRL 116, 113001 (2016)]. Such broadening has serious implications for many proposals to coherently use Rydberg interactions, particularly Rydberg dressing proposals. The dephasing arises as a runaway process where the production of the first contaminant atoms facilitates the creation of more contaminant atoms. Here we study the time dependence of this process with stroboscopic approaches. Using a pump-probe technique, we create an excess "pump" Rydberg population and probe its effect with a different "probe" Rydberg transition. We observe a reduced resonant pumping rate and an enhancement of the excitation on both sides of the transition as atoms are added to the pump state. We also observe a timescale for population growth significantly shorter than predicted by homogeneous mean-field models, as expected from a clustered growth mechanism where high-order correlations dominate the dynamics. These results support earlier works and confirm that the time scale for the onset of dephasing is reduced by a factor which scales as the inverse of the atom number. In addition, we discuss several approaches to minimize these effects of spontaneous broadening, including stroboscopic techniques and operating at cryogenic temperatures. It is challenging to avoid the unwanted broadening effects, but under some conditions they can be mitigated.


title: Distance to the Nearest Stable Metzler Matrix

abstract: This paper considers the non-convex problem of finding the nearest Metzler matrix to a given possibly unstable matrix. Linear systems whose state vector evolves according to a Metzler matrix have many desirable properties in analysis and control with regard to scalability. This motivates the question, how close (in the Frobenius norm of coefficients) to the nearest Metzler matrix are we? Dropping the Metzler constraint, this problem has recently been studied using the theory of dissipative Hamiltonian (DH) systems, which provide a helpful characterization of the feasible set of stable matrices. This work uses the DH theory to provide a block coordinate descent algorithm consisting of a quadratic program with favourable structural properties and a semidefinite program for which recent diagonal dominance results can be used to improve tractability.


title: Uniform level set estimates for ground state eigenfunctions

abstract: We study the behaviour of the first eigenfunction of the Dirichlet Laplacian on a planar convex domain near its maximum. We show that the eccentricity and orientation of the superlevel sets of the eigenfunction stabilise as they approach the maximum, uniformly with respect to the eccentricity of the domain itself. This is achieved by obtaining quantitatively sharp second derivative estimates, which are consistent with the shape of the superlevel sets. In particular, we prove that the eigenfunction is concave (rather than merely log concave) in an entire superlevel set near its maximum. By estimating the mixed second and third derivatives partial derivatives of the eigenfunction, we also determine the rate at which a degree 2 Taylor polynomial approximates the eigenfunction itself.


title: Local Neighborhood Intensity Pattern: A new texture feature descriptor for image retrieval

abstract: In this paper, a new texture descriptor based on the local neighborhood intensity difference is proposed for content based image retrieval (CBIR). For computation of texture features like Local Binary Pattern (LBP), the center pixel in a 3*3 window of an image is compared with all the remaining neighbors, one pixel at a time to generate a binary bit pattern. It ignores the effect of the adjacent neighbors of a particular pixel for its binary encoding and also for texture description. The proposed method is based on the concept that neighbors of a particular pixel hold a significant amount of texture information that can be considered for efficient texture representation for CBIR. Taking this into account, we develop a new texture descriptor, named as Local Neighborhood Intensity Pattern (LNIP) which considers the relative intensity difference between a particular pixel and the center pixel by considering its adjacent neighbors and generate a sign and a magnitude pattern. Since sign and magnitude patterns hold complementary information to each other, these two patterns are concatenated into a single feature descriptor to generate a more concrete and useful feature descriptor. The proposed descriptor has been tested for image retrieval on four databases, including three texture image databases - Brodatz texture image database, MIT VisTex database and Salzburg texture database and one face database AT&T face database. The precision and recall values observed on these databases are compared with some state-of-art local patterns. The proposed method showed a significant improvement over many other existing methods.


title: The diffuse molecular component in the nuclear bulge of the Milky Way

abstract: Context: The bulk of the Molecular gas in the Central Molecular Zone (CMZ) of the Galactic center region shows warm kinetic temperatures, ranging from $>20$ K in the coldest and densest regions (n$\sim 10^{4-5}$ cm$^{-3}$) up to more than 100 K for densities of about n$\sim 10^3$ cm$^{-3}$. Recently, a more diffuse, hotter ($n \sim 100$ cm$^{-3}$, $T\sim 250$ K) gas component was discovered through absorption observations of H$_3^+$. This component may be widespread in the Galactic center, and low density gas detectable in absorption may be present even outside the CMZ along sightlines crossing the extended bulge of the Galaxy. Aims: We aim to observe and characterize diffuse and low density gas using observations of 3-mm molecular transitions seen in absorption. Methods: Using the Atacama Large (sub)Millimeter Array (ALMA) we observed the absorption against the quasar J1744-312, which is located towards the Galactic bulge region at $(l,b)=(-2{^\circ}.13,-1{^\circ}.0)$, but outside the main molecular complexes. Results: ALMA observations in absorption against the J1744-312 quasar reveal a rich and complex chemistry in low density molecular and presumably diffuse clouds. We detected three velocity components, at $\sim $ 0, $-153$, and $-192$ km s$^{-1}$. The component at $\sim 0$ km s$^{-1}$ could represent gas in the Galactic disk while the velocity components at $-153$, and $-192$ km s$^{-1}$ likely originate from the Galactic bulge. We detected 12 molecules in the survey, but only 7 in the Galactic bulge gas.


title: Hadamard Full Propelinear Codes of type Q. Rank and Kernel

abstract: Hadamard full propelinear codes (HFP-codes) are introduced and their equivalence with Hadamard groups is proven (on the other hand, it is already known the equivalence of Hadamard groups with relative $(4n,2,4n,2n)$-difference sets in a group and also with cocyclic Hadamard matrices). We compute the available values for the rank and dimension of the kernel of HFP-codes of type Q and we show that the dimension of the kernel is always 1 or $2$. We also show that when the dimension of the kernel is 2 then the dimension of the kernel of the transposed code is 1 (so, both codes are not equivalent). Finally, we give a construction method such that from an HFP-code of length $4n$, dimension of the kernel $k=2$, and maximum rank $r=2n$, we obtain an HFP-code of double length $8n$, dimension of the kernel $k=2$, and maximum rank $r=4n$.


title: Invicem Lorentz Oscillator Model (ILOM)

abstract: We present detailed theoretical model of electron motion inside of an atom or molecule when subjected to electro-magnetic wave. This theoretical model, that we called Invicem Lorentz Oscillator Model (ILOM), contains Lorentz Oscillator Model (LOM) as a simplified approximation. The ILOM demonstrates artificial nature of the separation of descriptions of electromagnetic wave propagation in media into Linear and Nonlinear Optics.


title: On the classification of automorphisms of trees

abstract: We identify the complexity of the classification problem for automorphisms of a given countable regularly branching tree up to conjugacy. We consider both the rooted and unrooted cases. Additionally, we calculate the complexity of the conjugacy problem in the case of automorphisms of several non-regularly branching trees.


title: Conservative finite volume schemes for the quasi-geostrophic equation on coastal-conforming unstructured primal-dual meshes

abstract: In this paper we propose finite volume schemes for solving the inviscid and viscous quasi-geostrophic equations on coastal-conforming unstructured primal-dual meshes. Several approaches for enforcing the boundary conditions are also explored along with these schemes. The pure transport part in these schemes are shown to conserve the potential vorticity along fluid paths in an averaged sense, and conserve the potential enstrophy up to the time truncation errors. Numerical tests based on the centroidal Voronoi-Delaunay meshes are performed to confirm these properties, and to distinguish the dynamical behaviors of these schemes. Finally some potential applications of these schemes in different situations are discussed.


title: Tidal Love numbers and moment-Love relations of polytropic stars

abstract: The physical significance of tidal deformation in astronomical systems has long been known. The recently discovered universal I-Love-Q relations, which connect moment of inertia, quadrupole tidal Love number, and spin-induced quadrupole moment of compact stars, also underscore the special role of tidal deformation in gravitational wave astronomy. Motivated by the observation that such relations also prevail in Newtonian stars and crucially depend on the stiffness of a star, we consider the tidal Love numbers of Newtonian polytropic stars whose stiffness is characterised by a polytropic index $n$. We first perturbatively solve the Lane-Emden equation governing the profile of polytropic stars through the application of the scaled delta expansion method and then formulate perturbation series for the multipolar tidal Love number about the two exactly solvable cases with $n=0$ and $n=1$, respectively. Making use of these two series to form a two-point Pad\'e approximant, we find an approximate expression of the quadrupole tidal Love number, whose error is less than $2.5 \times 10^{-5}$ per cent (0.39 per cent) for $n\in[0,1]$ ($n\in[0,3]$). Similarly, we also determine the mass moments for polytropic stars accurately. Based on these findings, we are able to show that the I-Love-Q relations are in general stationary about the incompressible limit irrespective of the equation of state (EOS) of a star. Moreover, for the I-Love-Q relations, there is a secondary stationary point near $n \approx 0.4444$, thus showing the insensitivity to $n$ for $n\in[0,1]$. Our investigation clearly tracks the universality of the I-Love-Q relations from their validity for stiff stars such as neutron stars to their breakdown for soft stars.


title: Statistical Physics and Representations in Real and Artificial Neural Networks

abstract: This document presents the material of two lectures on statistical physics and neural representations, delivered by one of us (R.M.) at the Fundamental Problems in Statistical Physics XIV summer school in July 2017. In a first part, we consider the neural representations of space (maps) in the hippocampus. We introduce an extension of the Hopfield model, able to store multiple spatial maps as continuous, finite-dimensional attractors. The phase diagram and dynamical properties of the model are analyzed. We then show how spatial representations can be dynamically decoded using an effective Ising model capturing the correlation structure in the neural data, and compare applications to data obtained from hippocampal multi-electrode recordings and by (sub)sampling our attractor model. In a second part, we focus on the problem of learning data representations in machine learning, in particular with artificial neural networks. We start by introducing data representations through some illustrations. We then analyze two important algorithms, Principal Component Analysis and Restricted Boltzmann Machines, with tools from statistical physics.


title: Parameterization of SWIM Mobility Model Using Contact Traces

abstract: Opportunistic networks (OppNets) are focused to exploit direct, localised communications which occur in a peer-to-peer manner mostly based on people's movements and their contact durations. Therefore, the use of realistic mobility models is critical to evaluate the data dissemination in OppNets. One of the mobility models that is available in OMNeT++ which can be used to mimic human movement patterns is Small Worlds in Motion (SWIM). The SWIM model is based on the intuition that humans often visit nearby locations and if the visited location is far away, then it is probably due to the popularity of the location. As an alternative to mobility of a node, pairwise contact probabilities are also used to evaluate the data dissemination in OppNets. Pairwise contact probabilities can be used to predict that a node will be met by a particular node. These probabilities can be derived in many ways. One of the ways is to calculate the average probability with which a node will meet another particular node at any point of time. Another way is to calculate the probability with which a node will meet another based on the time of day. The way of calculating pairwise contact probability depends on the scenario. In this work, the pairwise contact probabilities obtained from the real traces are used to tune the parameters of the SWIM mobility model. The traces and the SWIM model are compared in terms of contact durations, inter-contact times and, number of pairwise contacts. How to decide SWIM parameters using real contact traces are being addressed.


title: Characterization of Extreme Copulas

abstract: In this paper our aim is to characterize the set of extreme points of the set of all n-dimensional copulas (n > 1). We have shown that a copula must induce a singular measure with respect to Lebesgue measure in order to be an extreme point in the set of n-dimensional copulas. We also have discovered some sufficient conditions for a copula to be an extreme copula. We have presented a construction of a small subset of n-dimensional extreme copulas such that any n-dimensional copula is a limit point of that subset with respect to weak convergence. The applications of such a theory are widespread, finding use in many facets of current mathematical research, such as distribution theory, survival analysis, reliability theory and optimization purposes. To illustrate the point further, examples of how such extremal representations can help in optimization have also been included.


title: Bispectrum Supersample Covariance

abstract: Modes with wavelengths larger than the survey window can have significant impact on the covariance within the survey window. The supersample covariance has been recognized as an important source of covariance for the power spectrum on small scales, and it can potentially be important for the bispectrum covariance as well. In this paper, using the response function formalism, we model the supersample covariance contributions to the bispectrum covariance and the cross covariance between the power spectrum and the bispectrum. The supersample covariances due to the long wavelength density and tidal perturbations are investigated, and the tidal contribution is a few orders of magnitude smaller than the density one because in configuration space the bispectrum estimator involves angular averaging and the tidal response function is anisotropic. The impact of the super-survey modes is quantified using numerical measurements with periodic box and subbox setups. For the matter bispectrum, the ratio between the supersample covariance correction and the small scale covariance, which can be computed using a periodic box, is roughly an order of magnitude smaller than that for the matter power spectrum. This is because for the bispectrum, the small scale non-Gaussian covariance is significantly larger than that for the power spectrum. For the cross covariance, the supersample covariance is as important as for the power spectrum covariance. The supersample covariance prediction with the halo model response function is in good agreement with numerical results.


title: Non-axially symmetric solutions of a mean field equation on $\mathbb{S}^2$

abstract: We prove the existence of a family of blow-up solutions of a mean field equation on sphere. The solutions blow up at four points where the minimum value of a potential energy function (involving the Green's function) is attained. The four blow-up points form a regular tetrahedron. Moreover, the solutions we build have a group of symmetry $T_d$ which is isomorphic to the symmetric group $S_4$. Other families of solutions can be similarly constructed with blow-up points at the vertices of equilateral triangles on a great circle or other inscribed platonic solids (cubes, octahedrons, icosahedrons and dodecahedrons). All of these solutions have the symmetries of the corresponding configuration, while they are non-axially symmetric.


title: Maximum independent sets near the upper bound

abstract: The size of a largest independent set of vertices in a given graph $G$ is denoted by $\alpha(G)$ and is called its independence number (or stability number). Given a graph $G$ and an integer $K,$ it is NP-complete to decide whether $\alpha(G) \geq K.$ An upper bound for the independence number $\alpha(G)$ of a given graph $G$ with $n$ vertices and $m $ edges is given by $\alpha(G) \leq p:=\lfloor\frac{1}{2} + \sqrt{\frac{1}{4} + n^2 - n - 2m}\rfloor.$ In this paper we will consider maximum independent sets near this upper bound. Our main result is the following: There exists an algorithm with time complexity $O(n^2)$ that, given as an input a graph $G$ with $n$ vertices, $m$ edges, $p:=\lfloor\frac{1}{2} + \sqrt{\frac{1}{4} + n^2 - n - 2m}\rfloor,$ and an integer $k \geq 0$ with $p \geq 2k+1,$ returns an induced subgraph $G_{p,k}$ of $G$ with $n_0 \leq p+2k+1$ vertices such that $\alpha(G) \leq p-k$ if and only if $\alpha(G_{p,k}) \leq p-k.$ Furthermore, we will show that we can decide in time $O(1.2738^{3k} + kn)$ whether $\alpha(G_{p,k}) \leq p-k.$


title: Fine-grained Recognition in the Wild: A Multi-Task Domain Adaptation Approach

abstract: While fine-grained object recognition is an important problem in computer vision, current models are unlikely to accurately classify objects in the wild. These fully supervised models need additional annotated images to classify objects in every new scenario, a task that is infeasible. However, sources such as e-commerce websites and field guides provide annotated images for many classes. In this work, we study fine-grained domain adaptation as a step towards overcoming the dataset shift between easily acquired annotated images and the real world. Adaptation has not been studied in the fine-grained setting where annotations such as attributes could be used to increase performance. Our work uses an attribute based multi-task adaptation loss to increase accuracy from a baseline of 4.1% to 19.1% in the semi-supervised adaptation case. Prior do- main adaptation works have been benchmarked on small datasets such as [46] with a total of 795 images for some domains, or simplistic datasets such as [41] consisting of digits. We perform experiments on a subset of a new challenging fine-grained dataset consisting of 1,095,021 images of 2, 657 car categories drawn from e-commerce web- sites and Google Street View.


title: Inferring Generative Model Structure with Static Analysis

abstract: Obtaining enough labeled data to robustly train complex discriminative models is a major bottleneck in the machine learning pipeline. A popular solution is combining multiple sources of weak supervision using generative models. The structure of these models affects training label quality, but is difficult to learn without any ground truth labels. We instead rely on these weak supervision sources having some structure by virtue of being encoded programmatically. We present Coral, a paradigm that infers generative model structure by statically analyzing the code for these heuristics, thus reducing the data required to learn structure significantly. We prove that Coral's sample complexity scales quasilinearly with the number of heuristics and number of relations found, improving over the standard sample complexity, which is exponential in $n$ for identifying $n^{\textrm{th}}$ degree relations. Experimentally, Coral matches or outperforms traditional structure learning approaches by up to 3.81 F1 points. Using Coral to model dependencies instead of assuming independence results in better performance than a fully supervised model by 3.07 accuracy points when heuristics are used to label radiology data without ground truth labels.


title: Closure Properties in the Class of Multiple Context Free Groups

abstract: We show that the class of groups with $k$-multiple context-free word problem is closed under graphs of groups with finite edge groups.


title: Context-Driven Exploration of Complex Chemical Reaction Networks

abstract: The construction of a reaction network containing all relevant intermediates and elementary reactions is necessary for the accurate description of chemical processes. In the case of a complex chemical reaction (involving, for instance, many reactants or highly reactive species), the size of such network may grow rapidly. Here, we present a computational protocol that constructs such reaction networks in a fully automated fashion steered in an intuitive, graph-based fashion through a single graphical user interface. Starting from a set of initial reagents new intermediates are explored through intra- and intermolecular reactions of already explored intermediates or new reactants presented to the network. This is done by assembling reactive complexes based on heuristic rules derived from conceptual electronic-structure theory and exploring the corresponding approximate reaction path. A subsequent path refinement leads to a minimum-energy path which connects the new intermediate to the existing ones to form a connected reaction network. Tree traversal algorithms are then employed to detect reaction channels and catalytic cycles. We apply our protocol to the formose reaction to study different pathways of sugar formation and to rationalize its autocatalytic nature.


title: Fine-Grained Car Detection for Visual Census Estimation

abstract: Targeted socioeconomic policies require an accurate understanding of a country's demographic makeup. To that end, the United States spends more than 1 billion dollars a year gathering census data such as race, gender, education, occupation and unemployment rates. Compared to the traditional method of collecting surveys across many years which is costly and labor intensive, data-driven, machine learning driven approaches are cheaper and faster--with the potential ability to detect trends in close to real time. In this work, we leverage the ubiquity of Google Street View images and develop a computer vision pipeline to predict income, per capita carbon emission, crime rates and other city attributes from a single source of publicly available visual data. We first detect cars in 50 million images across 200 of the largest US cities and train a model to predict demographic attributes using the detected cars. To facilitate our work, we have collected the largest and most challenging fine-grained dataset reported to date consisting of over 2600 classes of cars comprised of images from Google Street View and other web sources, classified by car experts to account for even the most subtle of visual differences. We use this data to construct the largest scale fine-grained detection system reported to date. Our prediction results correlate well with ground truth income data (r=0.82), Massachusetts department of vehicle registration, and sources investigating crime rates, income segregation, per capita carbon emission, and other market research. Finally, we learn interesting relationships between cars and neighborhoods allowing us to perform the first large scale sociological analysis of cities using computer vision techniques.


title: An explanation for dark matter and dark energy consistent with the Standard Model of particle physics and General Relativity

abstract: Analyses of internal galaxy and cluster dynamics typically employ Newton's law of gravity, which neglects the field self-interaction effects of General Relativity. This may be why dark matter seems necessary. The Universe evolution, on the other hand, is treated with the full theory, General Relativity. However, the approximations of isotropy and homogeneity, normally used to derive and solve the Universe evolution equations, effectively suppress General Relativity's field self-interaction effects and this may introduce the need for dark energy. Calculations have shown that field self-interaction increases the binding of matter inside massive systems, which may account for galaxy and cluster dynamics without invoking dark matter. In turn, energy conservation dictates that the increased binding must be balanced by an effectively decreased gravitational interaction outside the massive system. In this article, such suppression is estimated and its consequence for the Universe's evolution is discussed. Observations are reproduced without need for dark energy.


title: Scalable Annotation of Fine-Grained Categories Without Experts

abstract: We present a crowdsourcing workflow to collect image annotations for visually similar synthetic categories without requiring experts. In animals, there is a direct link between taxonomy and visual similarity: e.g. a collie (type of dog) looks more similar to other collies (e.g. smooth collie) than a greyhound (another type of dog). However, in synthetic categories such as cars, objects with similar taxonomy can have very different appearance: e.g. a 2011 Ford F-150 Supercrew-HD looks the same as a 2011 Ford F-150 Supercrew-LL but very different from a 2011 Ford F-150 Supercrew-SVT. We introduce a graph based crowdsourcing algorithm to automatically group visually indistinguishable objects together. Using our workflow, we label 712,430 images by ~1,000 Amazon Mechanical Turk workers; resulting in the largest fine-grained visual dataset reported to date with 2,657 categories of cars annotated at 1/20th the cost of hiring experts.


title: Mathematical Symbolism in a Russian Literary Masterpiece

abstract: Andrei Bely's modernist novel "Petersburg," first published in 1913, is considered a pinnacle of the Symbolist movement. Nabokov famously ranked it as one of the four greatest masterpieces of 20th-century prose. The author's father, Bugaev, was an influential mathematician and for 12 years served as the president of the Moscow Mathematical Society; he was also a source of inspiration for one of the main characters in the son's novel. While the philosophical views and political leanings of the mathematicians surrounding Bely, and their impact on "Petersburg," have been a topic of recent academic interest, there has not yet been a direct investigation of the surprisingly frequent and sophisticated mathematical passages in the book itself. We attempt here to rectify this gap in the scholarly literature, and in doing so find a rich tapestry of mathematical ideas and allusions.


title: Exact coherent states with hairpin-like vortex structure in channel flow

abstract: Hairpin vortices are widely studied as an important structural aspect of wall turbulence. The present work describes, for the first time, nonlinear traveling wave solutions to the Navier--Stokes equations in the channel flow geometry -- exact coherent states (ECS) -- that display hairpin-like vortex structure. This solution family comes into existence at a saddle-node bifurcation at Reynolds number Re=666. At the bifurcation, the solution has a highly symmetric quasistreamwise vortex structure similar to that reported for previously studied ECS. With increasing distance from the bifurcation, however, both the upper and lower branch solutions develop a vortical structure characteristic of hairpins: a spanwise-oriented "head" near the channel centerplane where the mean shear vanishes connected to counter-rotating quasistreamwise "legs" that extend toward the channel wall. At Re=1800, the upper branch solution has mean and Reynolds shear-stress profiles that closely resemble those of turbulent mean profiles in the same domain.


title: On the height of solutions to norm form equations

abstract: Let $k$ be a number field. We consider norm form equations associated to a full $O_k$-module contained in a finite extension field $l$. It is known that the set of solutions is naturally a union of disjoint equivalence classes of solutions. We prove that each nonempty equivalence class of solutions contains a representative with Weil height bounded by an expression that depends on parameters defining the norm form equation.


title: A General Regularized Continuous Formulation for the Maximum Clique Problem

abstract: In this paper, we develop a general regularization-based continuous optimization framework for the maximum clique problem. In particular, we consider a broad class of regularization terms that can be included in the classic Motzkin-Strauss formulation and we develop conditions that guarantee the equivalence between the continuous regularized problem and the original one in both a global and a local sense. We further analyze, from a computational point of view, two different regularizers that satisfy the general conditions.


title: Type II supernova spectral diversity I: Observations, sample characterization and spectral line evolution

abstract: We present 888 visual-wavelength spectra of 122 nearby type II supernovae (SNe II) obtained between 1986 and 2009, and ranging between 3 and 363 days post explosion. In this first paper, we outline our observations and data reduction techniques, together with a characterization based on the spectral diversity of SNe~II. A statistical analysis of the spectral matching technique is discussed as an alternative to non-detection constraints for estimating SN explosion epochs. The time evolution of spectral lines is presented and analysed in terms of how this differs for SNe of different photometric, spectral, and environmental properties: velocities, pseudo-equivalent widths, decline rates, magnitudes, time durations, and environment metallicity. Our sample displays a large range in ejecta expansion velocities, from $\sim9600$ to $\sim1500$ km s$^{-1}$ at 50 days post explosion with a median H$_{\alpha}$ value of 7300 km s$^{-1}$. This is most likely explained through differing explosion energies. Significant diversity is also observed in the absolute strength of spectral lines, characterised through their pseudo-equivalent widths. This implies significant diversity in both temperature evolution (linked to progenitor radius) and progenitor metallicity between different SNe~II. Around 60\% of our sample show an extra absorption component on the blue side of the H$_{\alpha}$ P-Cygni profile ("Cachito" feature) between 7 and 120 days since explosion. Studying the nature of Cachito, we conclude that these features at early times (before $\sim35$ days) are associated with \ion{Si}{2} $\lambda6355$, while past the middle of the plateau phase they are related to high velocity (HV) features of hydrogen lines.


title: Stochastic basis adaptation and spatial domain decomposition for PDEs with random coefficients

abstract: We present a novel uncertainty quantification approach for high-dimensional stochastic partial differential equations that reduces the computational cost of polynomial chaos methods by decomposing the computational domain into non-overlapping subdomains and adapting the stochastic basis in each subdomain so the local solution has a lower dimensional random space representation. The local solutions are coupled using the Neumann-Neumann algorithm, where we first estimate the interface solution then evaluate the interior solution in each subdomain using the interface solution as a boundary condition. The interior solutions in each subdomain are computed independently of each other, which reduces the operation count from $O(N^\alpha)$ to $O(M^\alpha),$ where $N$ is the total number of degrees of freedom, $M$ is the number of degrees of freedom in each subdomain, and the exponent $\alpha>1$ depends on the uncertainty quantification method used. In addition, the localized nature of solutions makes the proposed approach highly parallelizable. We illustrate the accuracy and efficiency of the approach for linear and nonlinear differential equations with random coefficients.


title: BlockSci: Design and applications of a blockchain analysis platform

abstract: Analysis of blockchain data is useful for both scientific research and commercial applications. We present BlockSci, an open-source software platform for blockchain analysis. BlockSci is versatile in its support for different blockchains and analysis tasks. It incorporates an in-memory, analytical (rather than transactional) database, making it several hundred times faster than existing tools. We describe BlockSci's design and present four analyses that illustrate its capabilities. This is a working paper that accompanies the first public release of BlockSci, available at https://github.com/citp/BlockSci. We seek input from the community to further develop the software and explore other potential applications.


title: Exploiting Problem Structure in Optimization under Uncertainty via Online Convex Optimization

abstract: In this paper, we consider two paradigms that are developed to account for uncertainty in optimization models: robust optimization (RO) and joint estimation-optimization (JEO). We examine recent developments on efficient and scalable iterative first-order methods for these problems, and show that these iterative methods can be viewed through the lens of online convex optimization (OCO). The standard OCO framework has seen much success for its ability to handle decision-making in dynamic, uncertain, and even adversarial environments. Nevertheless, our applications of interest present further flexibility in OCO via three simple modifications to standard OCO assumptions: we introduce two new concepts of weighted regret and online saddle point problems and study the possibility of making lookahead (anticipatory) decisions. Our analyses demonstrate that these flexibilities introduced into the OCO framework have significant consequences whenever they are applicable. For example, in the strongly convex case, minimizing unweighted regret has a proven optimal bound of $O(\log(T)/T)$, whereas we show that a bound of $O(1/T)$ is possible when we consider weighted regret. Similarly, for the smooth case, considering $1$-lookahead decisions results in a $O(1/T)$ bound, compared to $O(1/\sqrt{T})$ in the standard OCO setting. Consequently, these OCO tools are instrumental in exploiting structural properties of functions and resulting in improved convergence rates for RO and JEO. In certain cases, our results for RO and JEO match the best known or optimal rates in the corresponding problem classes without data uncertainty.


title: Pair-breaking of multi-gap superconductivity under parallel magnetic fields in electric-field-induced surface metallic state

abstract: Roles of paramagnetic and diamagnetic pair-breaking effects in superconductivity in electric-field-induced surface metallic state are studied by Bogoliubov-de Gennes equation, when magnetic fields are applied parallel to the surface. The multi-gap states of sub-bands are related to the depth dependence and the magnetic field dependence of superconductivity. In the Fermi-energy density of states and the spin density, sub-band contributions successively appear from higher-level sub-bands with increasing magnetic fields. The characteristic magnetic field dependence may be a key feature to identify the multi-gap structure of the surface superconductivity.


title: Metric thickenings of Euclidean submanifolds

abstract: Given a sample $Y$ from an unknown manifold $X$ embedded in Euclidean space, it is possible to recover the homology groups of $X$ by building a Vietoris--Rips or \v{C}ech simplicial complex on top of the vertex set $Y$. However, these simplicial complexes need not inherit the metric structure of the manifold, in particular when $Y$ is infinite. Indeed, a simplicial complex is not even metrizable if it is not locally finite. We instead consider metric thickenings, called the \emph{Vietoris--Rips} and \emph{\v{C}ech thickenings}, which are equipped with the 1-Wasserstein metric in place of the simplicial complex topology. We show that for Euclidean subsets $X$ with positive reach, the thickenings satisfy metric analogues of Hausmann's theorem and the nerve lemma (the metric Vietoris--Rips and \v{C}ech thickenings of $X$ are homotopy equivalent to $X$ for scale parameters less than the reach). To our knowledge this is the first version of Hausmann's theorem for Euclidean submanifolds (as opposed to Riemannian manifolds), and our result also extends to non-manifold shapes (as not all sets of positive reach are manifolds). In contrast to Hausmann's original proof, our homotopy equivalence is a deformation retraction, is realized by canonical maps in both directions, and furthermore can be proven to be a homotopy equivalence via simple linear homotopies from the map compositions to the corresponding identity maps.


title: An accurate and efficient method for calculating fluid exchange between fractures and matrix with a non-conforming mesh

abstract: I propose the 1D Subgrid Method for calculating fluid exchange between fractures and matrix with a non-conforming mesh. The method is demonstrated on linear and radial flow problems. The method has negligible computational cost and is accurate at all levels of mesh refinement.