title: Optimal bounds and extremal trajectories for time averages in nonlinear dynamical systems

abstract: For any quantity of interest in a system governed by ordinary differential equations, it is natural to seek the largest (or smallest) long-time average among solution trajectories, as well as the extremal trajectories themselves. Upper bounds on time averages can be proved a priori using auxiliary functions, the optimal choice of which is a convex optimization problem. We prove that the problems of finding maximal trajectories and minimal auxiliary functions are strongly dual. Thus, auxiliary functions provide arbitrarily sharp upper bounds on time averages. Moreover, any nearly minimal auxiliary function provides phase space volumes in which all nearly maximal trajectories are guaranteed to lie. For polynomial equations, auxiliary functions can be constructed by semidefinite programming, which we illustrate using the Lorenz system.


title: Infinite dimensional semiclassical analysis and applications to a model in NMR

abstract: We are interested in this paper with the connection between the dynamics of a model related to Nuclear Magnetic Resonance (NMR) in Quantum Field Theory (QFT) with its classical counterpart known as the Maxwell-Bloch equations. The model in QFT is a model of Quantum Electrodynamics (QED) considering fixed spins interacting with the quantized electromagnetic field in an external constant magnetic field. This model is close to the common spin-boson model. The classical model goes back to F. Bloch [15] in 1946. Our goal is not only to study the derivation of the Maxwell-Bloch equations but to also establish a semiclassical asymptotic expansion of arbitrary high orders with control of the error terms of this standard nonlinear classical motion equations. This provides therefore quantum corrections of any order in powers of the semiclassical parameter of the Bloch equations. Besides, the asymptotic expansion for the photon number is also analyzed and a law describing the photon number time evolution is written down involving the radiation field polarization. Since the quantum photon state Hilbert space (radiation field) is infinite dimensional we are thus concerned in this article with the issue of semiclassical calculus in an infinite dimensional setting. In this regard, we are studying standard notions as Wick and anti-Wick quantizations, heat operator, Beals characterization theorem and compositions of symbols in the infinite dimensional context which can have their own interest.


title: Exoplanet Biosignatures: Observational Prospects

abstract: Exoplanet hunting efforts have revealed the prevalence of exotic worlds with diverse properties, including Earth-sized bodies, which has fueled our endeavor to search for life beyond the Solar System. Accumulating experiences in astrophysical, chemical, and climatological characterization of uninhabitable planets are paving the way to characterization of potentially habitable planets. In this paper, we review our possibilities and limitations in characterizing temperate terrestrial planets with future observational capabilities through 2030s and beyond, as a basis of a broad range of discussions on how to advance "astrobiology" with exoplanets. We discuss the observability of not only the proposed biosignature candidates themselves, but also of more general planetary properties that provide circumstantial evidence, since the evaluation of any biosignature candidate relies on their context. Characterization of temperate Earth-size planets in the coming years will focus on those around nearby late-type stars. JWST and later 30 meter-class ground-based telescopes will empower their chemical investigations. Spectroscopic studies of potentially habitable planets around solar-type stars will likely require a designated spacecraft mission for direct imaging, leveraging technologies that are already being developed and tested as part of the WFIRST mission. Successful initial characterization of a few nearby targets will be an important touchstone toward a more detailed scrutiny and a larger survey that are envisioned beyond 2030. The broad outlook this paper presents may help develop new observational techniques to detect relevant features as well as frameworks to diagnose planets based on the observables.


title: Machine learning for classification and quantification of monoclonal antibody preparations for cancer therapy

abstract: Monoclonal antibodies constitute one of the most important strategies to treat patients suffering from cancers such as hematological malignancies and solid tumors. In order to guarantee the quality of those preparations prepared at hospital, quality control has to be developed. The aim of this study was to explore a noninvasive, nondestructive, and rapid analytical method to ensure the quality of the final preparation without causing any delay in the process. We analyzed four mAbs (Inlfiximab, Bevacizumab, Ramucirumab and Rituximab) diluted at therapeutic concentration in chloride sodium 0.9% using Raman spectroscopy. To reduce the prediction errors obtained with traditional chemometric data analysis, we explored a data-driven approach using statistical machine learning methods where preprocessing and predictive models are jointly optimized. We prepared a data analytics workflow and submitted the problem to a collaborative data challenge platform called Rapid Analytics and Model Prototyping (RAMP). This allowed to use solutions from about 300 data scientists during five days of collaborative work. The prediction of the four mAbs samples was considerably improved with a misclassification rate and the mean error rate of 0.8% and 4%, respectively.


title: New Cataclysmic Variables and other Exotic Binaries in the Globular Cluster 47 Tucanae

abstract: We present 22 new (+3 confirmed) cataclysmic variables (CVs) in the non core-collapsed globular cluster 47 Tucanae (47 Tuc). The total number of CVs in the cluster is now 43, the largest sample in any globular cluster so far. For the identifications we used near-ultraviolet (NUV) and optical images from the Hubble Space Telescope, in combination with X-ray results from the Chandra X-ray Observatory. This allowed us to build the deepest NUV CV luminosity function of the cluster to date. We found that the CVs in 47 Tuc are more concentrated towards the cluster center than the main sequence turnoff stars. We compared our results to the CV populations of the core-collapsed globular clusters NGC 6397 and NGC 6752. We found that 47 Tuc has fewer bright CVs per unit mass than those two other clusters. That suggests that dynamical interactions in core-collapsed clusters play a major role creating new CVs. In 47 Tuc, the CV population is probably dominated by primordial and old dynamically formed systems. We estimated that the CVs in 47 Tuc have total masses of approx. 1.4 M_sun. We also found that the X-ray luminosity function of the CVs in the three clusters is bimodal. Additionally, we discuss a possible double degenerate system and an intriguing/unclassified object. Finally, we present four systems that could be millisecond pulsar companions given their X-ray and NUV/optical colors. For one of them we present very strong evidence for being an ablated companion. The other three could be CO- or He-WDs.


title: Feature-rich bifurcations in a simple electronic circuit

abstract: A simple electronic circuit with a voltage controlled current source is investigated. The circuit exhibits rich dynamics upon varying the circuit elements such as L,C and R, and the control factor of the current source. Among several other interesting features, the circuit demonstrates two local bifurcations, namely, node to spiral and Hopf bifurcation, and a global homoclinic bifurcation. Phase-portraits corresponding to these bifurcations are presented and the implications of these bifurcations on system stability are discussed. In particular, the circuit parameters corresponding to the onset of Hopf bifurcation may be exploited to design an oscillator with stable frequency and amplitude. The circuit may be easily implemented with nonlinear resistive elements such as diodes or transistors in saturation and a gyrator block as the voltage controlled current source.


title: The time measurement problem in quantum cosmology

abstract: In the canonical approach to quantization of gravity, one often uses relational clock variables and an interpretation in terms of conditional probabilities to overcome the problem of time. In this essay we show that these suffer from serious conceptual issues.


title: Survey of Water and Ammonia in Nearby galaxies (SWAN): Resolved Ammonia Thermometry, Water and Methanol Masers in the Nuclear Starburst of NGC 253

abstract: We present Karl G Jansky Very Large Array molecular line observations of the nearby starburst galaxy NGC 253, from SWAN: "Survey of Water and Ammonia in Nearby galaxies". SWAN is a molecular line survey at centimeter wavelengths designed to reveal the physical conditions of star forming gas over a range of star forming galaxies. NGC 253 has been observed in four 1GHz bands from 21 to 36 GHz at 6" ($\sim100$pc) spatial and 3.5 km s$^{-1}$ spectral resolution. In total we detect 19 transitions from seven molecular and atomic species. We have targeted the metastable inversion transitions of ammonia (NH$_{3}$) from (1,1) to (5,5) and the (9,9) line, the 22.2 GHz water (H$_2$O) ($6_{16}-5_{23}$) maser, and the 36.1 GHz methanol (CH$_3$OH) ($4_{-1}-3_{0}$) maser. Utilizing NH$_{3}$ as a thermometer, we present evidence for uniform heating over the central kpc of NGC 253. The molecular gas is best described by a two kinetic temperature model with a warm 130K and a cooler 57K component. A comparison of these observations with previous ALMA results suggests that the molecular gas is not heated in photon dominated regions or shocks. It is possible that the gas is heated by turbulence or cosmic rays. In the galaxy center we find evidence for NH$_{3}$(3,3) masers. Furthermore we present velocities and luminosities of three water maser features related to the nuclear starburst. We partially resolve CH$_3$OH masers seen at the edges of the bright molecular emission, which coincides with expanding molecular superbubbles. This suggests that the masers are pumped by weak shocks in the bubble surfaces.


title: Efficient Learning of Harmonic Priors for Pitch Detection in Polyphonic Music

abstract: Automatic music transcription (AMT) aims to infer a latent symbolic representation of a piece of music (piano-roll), given a corresponding observed audio recording. Transcribing polyphonic music (when multiple notes are played simultaneously) is a challenging problem, due to highly structured overlapping between harmonics. We study whether the introduction of physically inspired Gaussian process (GP) priors into audio content analysis models improves the extraction of patterns required for AMT. Audio signals are described as a linear combination of sources. Each source is decomposed into the product of an amplitude-envelope, and a quasi-periodic component process. We introduce the Mat\'ern spectral mixture (MSM) kernel for describing frequency content of singles notes. We consider two different regression approaches. In the sigmoid model every pitch-activation is independently non-linear transformed. In the softmax model several activation GPs are jointly non-linearly transformed. This introduce cross-correlation between activations. We use variational Bayes for approximate inference. We empirically evaluate how these models work in practice transcribing polyphonic music. We demonstrate that rather than encourage dependency between activations, what is relevant for improving pitch detection is to learnt priors that fit the frequency content of the sound events to detect.


title: The Bag Semantics of Ontology-Based Data Access

abstract: Ontology-based data access (OBDA) is a popular approach for integrating and querying multiple data sources by means of a shared ontology. The ontology is linked to the sources using mappings, which assign views over the data to ontology predicates. Motivated by the need for OBDA systems supporting database-style aggregate queries, we propose a bag semantics for OBDA, where duplicate tuples in the views defined by the mappings are retained, as is the case in standard databases. We show that bag semantics makes conjunctive query answering in OBDA coNP-hard in data complexity. To regain tractability, we consider a rather general class of queries and show its rewritability to a generalisation of the relational calculus to bags.


title: Elementary example of energy and momentum of an extended physical system in special relativity

abstract: An instructive paradox concerning classical description of energy and momentum of extended physical systems in special relativity theory is explained using an elementary example of two point-like massive bodies rotating on a circle in their center-of-mass frame of reference, connected by an arbitrarily light and infinitesimally thin string. Namely, from the point of view of the inertial observers who move with respect to the rotating system, the sums of energies and momenta of the two bodies oscillate, instead of being constant in time. This result is understood in terms of the mechanism that binds the bodies: the string contributes to the system total energy and momentum no matter how light it is. Its contribution eliminates the unphysical oscillations from the system total four-momentum. Generality of the relativistic approach, applied here to the rotor example, suggests that in every extended physical system its binding mechanism contributes to its total energy and momentum.


title: Gradient Estimators for Implicit Models

abstract: Implicit models, which allow for the generation of samples but not for point-wise evaluation of probabilities, are omnipresent in real-world problems tackled by machine learning and a hot topic of current research. Some examples include data simulators that are widely used in engineering and scientific research, generative adversarial networks (GANs) for image synthesis, and hot-off-the-press approximate inference techniques relying on implicit distributions. The majority of existing approaches to learning implicit models rely on approximating the intractable distribution or optimisation objective for gradient-based optimisation, which is liable to produce inaccurate updates and thus poor models. This paper alleviates the need for such approximations by proposing the Stein gradient estimator, which directly estimates the score function of the implicitly defined distribution. The efficacy of the proposed estimator is empirically demonstrated by examples that include meta-learning for approximate inference, and entropy regularised GANs that provide improved sample diversity.


title: Snapshot Difference Imaging using Time-of-Flight Sensors

abstract: Computational photography encompasses a diversity of imaging techniques, but one of the core operations performed by many of them is to compute image differences. An intuitive approach to computing such differences is to capture several images sequentially and then process them jointly. Usually, this approach leads to artifacts when recording dynamic scenes. In this paper, we introduce a snapshot difference imaging approach that is directly implemented in the sensor hardware of emerging time-of-flight cameras. With a variety of examples, we demonstrate that the proposed snapshot difference imaging technique is useful for direct-global illumination separation, for direct imaging of spatial and temporal image gradients, for direct depth edge imaging, and more.


title: Deep adversarial neural decoding

abstract: Here, we present a novel approach to solve the problem of reconstructing perceived stimuli from brain responses by combining probabilistic inference with deep learning. Our approach first inverts the linear transformation from latent features to brain responses with maximum a posteriori estimation and then inverts the nonlinear transformation from perceived stimuli to latent features with adversarial training of convolutional neural networks. We test our approach with a functional magnetic resonance imaging experiment and show that it can generate state-of-the-art reconstructions of perceived faces from brain activations.


title: Cup products in the etale cohomology of number fields

abstract: This paper concerns cup product pairings in \'etale cohomology related to work of M. Kim and of W. McCallum and R. Sharifi. We will show that by considering Ext groups rather than cohomology groups, one arrives at a pairing which combines invariants defined by Kim with a pairing defined by McCallum and Sharifi. We also prove a formula for Kim's invariant in terms of Artin maps in the case of cyclic unramified Kummer extensions. One consequence is that for all $n > 1$, there are infinitely many number fields $F$ over which there are both trivial and non-trivial Kim invariants associated to cyclic groups of order $n$.


title: The Kernel Mixture Network: A Nonparametric Method for Conditional Density Estimation of Continuous Random Variables

abstract: This paper introduces the kernel mixture network, a new method for nonparametric estimation of conditional probability densities using neural networks. We model arbitrarily complex conditional densities as linear combinations of a family of kernel functions centered at a subset of training points. The weights are determined by the outer layer of a deep neural network, trained by minimizing the negative log likelihood. This generalizes the popular quantized softmax approach, which can be seen as a kernel mixture network with square and non-overlapping kernels. We test the performance of our method on two important applications, namely Bayesian filtering and generative modeling. In the Bayesian filtering example, we show that the method can be used to filter complex nonlinear and non-Gaussian signals defined on manifolds. The resulting kernel mixture network filter outperforms both the quantized softmax filter and the extended Kalman filter in terms of model likelihood. Finally, our experiments on generative models show that, given the same architecture, the kernel mixture network leads to higher test set likelihood, less overfitting and more diversified and realistic generated samples than the quantized softmax approach.


title: Fast Singular Value Shrinkage with Chebyshev Polynomial Approximation Based on Signal Sparsity

abstract: We propose an approximation method for thresholding of singular values using Chebyshev polynomial approximation (CPA). Many signal processing problems require iterative application of singular value decomposition (SVD) for minimizing the rank of a given data matrix with other cost functions and/or constraints, which is called matrix rank minimization. In matrix rank minimization, singular values of a matrix are shrunk by hard-thresholding, soft-thresholding, or weighted soft-thresholding. However, the computational cost of SVD is generally too expensive to handle high dimensional signals such as images; hence, in this case, matrix rank minimization requires enormous computation time. In this paper, we leverage CPA to (approximately) manipulate singular values without computing singular values and vectors. The thresholding of singular values is expressed by a multiplication of certain matrices, which is derived from a characteristic of CPA. The multiplication is also efficiently computed using the sparsity of signals. As a result, the computational cost is significantly reduced. Experimental results suggest the effectiveness of our method through several image processing applications based on matrix rank minimization with nuclear norm relaxation in terms of computation time and approximation precision.


title: Well-conditioned frames for high order finite element methods

abstract: The purpose of this paper is to discuss representations of high order $C^0$ finite element spaces on simplicial meshes in any dimension. When computing with high order piecewise polynomials the conditioning of the basis is likely to be important. The main result of this paper is a construction of representations by frames such that the associated $L^2$ condition number is bounded independently of the polynomial degree. To our knowledge, such a representation has not been presented earlier. The main tools we will use for the construction is the bubble transform, introduced previously in [Falk and Winther, Found Comput Math (2016) 16: 297], and properties of Jacobi polynomials on simplexes in higher dimensions. We also include a brief discussion of preconditioned iterative methods for the finite element systems in the setting of representations by frames.


title: A Comparison of Reinforcement Learning Techniques for Fuzzy Cloud Auto-Scaling

abstract: A goal of cloud service management is to design self-adaptable auto-scaler to react to workload fluctuations and changing the resources assigned. The key problem is how and when to add/remove resources in order to meet agreed service-level agreements. Reducing application cost and guaranteeing service-level agreements (SLAs) are two critical factors of dynamic controller design. In this paper, we compare two dynamic learning strategies based on a fuzzy logic system, which learns and modifies fuzzy scaling rules at runtime. A self-adaptive fuzzy logic controller is combined with two reinforcement learning (RL) approaches: (i) Fuzzy SARSA learning (FSL) and (ii) Fuzzy Q-learning (FQL). As an off-policy approach, Q-learning learns independent of the policy currently followed, whereas SARSA as an on-policy always incorporates the actual agent's behavior and leads to faster learning. Both approaches are implemented and compared in their advantages and disadvantages, here in the OpenStack cloud platform. We demonstrate that both auto-scaling approaches can handle various load traffic situations, sudden and periodic, and delivering resources on demand while reducing operating costs and preventing SLA violations. The experimental results demonstrate that FSL and FQL have acceptable performance in terms of adjusted number of virtual machine targeted to optimize SLA compliance and response time.


title: Multi-Task Learning Using Uncertainty to Weigh Losses for Scene Geometry and Semantics

abstract: Numerous deep learning applications benefit from multi-task learning with multiple regression and classification objectives. In this paper we make the observation that the performance of such systems is strongly dependent on the relative weighting between each task's loss. Tuning these weights by hand is a difficult and expensive process, making multi-task learning prohibitive in practice. We propose a principled approach to multi-task deep learning which weighs multiple loss functions by considering the homoscedastic uncertainty of each task. This allows us to simultaneously learn various quantities with different units or scales in both classification and regression settings. We demonstrate our model learning per-pixel depth regression, semantic and instance segmentation from a monocular input image. Perhaps surprisingly, we show our model can learn multi-task weightings and outperform separate models trained individually on each task.


title: Homogeneous Yang-Baxter deformations as generalized diffeomorphisms

abstract: Yang-Baxter (YB) deformations of string sigma model provide deformed target spaces. We propose that homogeneous YB deformations always lead to a certain class of $\beta$-twisted backgrounds and represent the bosonic part of the supergravity fields in terms of the classical r-matrix associated with the YB deformation. We then show that various $\beta$-twisted backgrounds can be realized by considering generalized diffeomorphisms in the undeformed background. Our result extends the notable relation between the YB deformations and (non-commuting) TsT transformations. We also discuss more general deformations beyond the YB deformations.


title: Deep Learning as a Tool to Predict Flow Patterns in Two-Phase Flow

abstract: In order to better model complex real-world data such as multiphase flow, one approach is to develop pattern recognition techniques and robust features that capture the relevant information. In this paper, we use deep learning methods, and in particular employ the multilayer perceptron, to build an algorithm that can predict flow pattern in twophase flow from fluid properties and pipe conditions. The preliminary results show excellent performance when compared with classical methods of flow pattern prediction.


title: Evaluation of Direct Haptic 4D Volume Rendering of Partially Segmented Data for Liver Puncture Simulation

abstract: This work presents an evaluation study using a force feedback evaluation framework for a novel direct needle force volume rendering concept in the context of liver puncture simulation. PTC/PTCD puncture interventions targeting the bile ducts have been selected to illustrate this concept. The haptic algorithms of the simulator system are based on (1) partially segmented patient image data and (2) a non-linear spring model effective at organ borders. The primary aim is to quantitatively evaluate force errors caused by our patient modeling approach, in comparison to haptic force output obtained from using gold-standard, completely manually-segmented data. The evaluation of the force algorithms compared to a force output from fully manually segmented gold-standard patient models, yields a low mean of 0.12 N root mean squared force error and up to 1.6 N for systematic maximum absolute errors. Force errors were evaluated on 31,222 preplanned test paths from 10 patients. Only twelve percent of the emitted forces along these paths were affected by errors. This is the first study evaluating haptic algorithms with deformable virtual patients in silico. We prove haptic rendering plausibility on a very high number of test paths. Important errors are below just noticeable differences for the hand-arm system.


title: On convex closed planar curves as equidistant sets

abstract: The equidistant set of two nonempty subsets $K$ and $L$ in the Euclidean plane is a set all of whose points have the same distance from $K$ and $L$. Since the classical conics can be also given in this way, equidistant sets can be considered as a kind of their generalizations: $K$ and $L$ are called the focal sets. In their paper \cite{PS} the authors posed the problem of the characterization of closed subsets in the Euclidean plane that can be realized as the equidistant set of two connected disjoint closed sets. We prove that any convex closed planar curve can be given as an equidistant set, i.e. the set of equidistant curves contains the entire class of convex closed planar curves. In this sense the equidistancy is a generalization of the convexity.


title: VAE with a VampPrior

abstract: Many different methods to train deep generative models have been introduced in the past. In this paper, we propose to extend the variational auto-encoder (VAE) framework with a new type of prior which we call "Variational Mixture of Posteriors" prior, or VampPrior for short. The VampPrior consists of a mixture distribution (e.g., a mixture of Gaussians) with components given by variational posteriors conditioned on learnable pseudo-inputs. We further extend this prior to a two layer hierarchical model and show that this architecture with a coupled prior and posterior, learns significantly better models. The model also avoids the usual local optima issues related to useless latent dimensions that plague VAEs. We provide empirical studies on six datasets, namely, static and binary MNIST, OMNIGLOT, Caltech 101 Silhouettes, Frey Faces and Histopathology patches, and show that applying the hierarchical VampPrior delivers state-of-the-art results on all datasets in the unsupervised permutation invariant setting and the best results or comparable to SOTA methods for the approach with convolutional networks.


title: BAMHealthCloud: A Biometric Authentication and Data Management System for Healthcare Data in Cloud

abstract: Advancements in healthcare industry with new technology and population growth has given rise to security threat to our most personal data. The healthcare data management system consists of records in different formats such as text, numeric, pictures and videos leading to data which is big and unstructured. Also, hospitals have several branches at different locations throughout a country and overseas. In view of these requirements a cloud based healthcare management system can be an effective solution for efficient health care data management. One of the major concerns of a cloud based healthcare system is the security aspect. It includes theft to identity, tax fraudulence, insurance frauds, medical frauds and defamation of high profile patients. Hence, a secure data access and retrieval is needed in order to provide security of critical medical records in health care management system. Biometric authentication mechanism is suitable in this scenario since it overcomes the limitations of token theft and forgetting passwords in conventional token id-password mechanism used for providing security. It also has high accuracy rate for secure data access and retrieval. In this paper we propose BAMHealthCloud which is a cloud based system for management of healthcare data, it ensures security of data through biometric authentication. It has been developed after performing a detailed case study on healthcare sector in a developing country. Training of the signature samples for authentication purpose has been performed in parallel on hadoop MapReduce framework using Resilient Backpropagation neural network. From rigorous experiments it can be concluded that it achieves a speedup of 9x, Equal error rate (EER) of 0.12, sensitivity of 0.98 and specificity of 0.95 as compared to other approaches existing in literature.


title: Noncommutative Blackwell-Ross martingale inequality

abstract: We establish a noncommutative Blackwell--Ross inequality for supermartingales under a suitable condition which generalize Khan's works to the noncommutative setting. We then employ it to deduce an Azuma-type inequality.


title: Projective quantum modules and projective ideals of C*-algebras

abstract: We introduce in non-coordinate presentation the notions of a quantum algebra and of a quantum module over such an algebra. Then we give the definition of a projective quantum module and of a free quantum module, the latter as a particular case of the notion of a free object in a rigged category. (Here we say "quantum" instead of frequently used protean adjective "operator"). After this we discuss the general connection between projectivity and freeness. Then we show that for a Banach quantum algebra A and a Banach quantum space E the Banach quantum A-module $A\widehat\otimes_{op}E$ is free, where " $\widehat\otimes_{op}$ " denotes the operator-projective tensor product of Banach quantum spaces. This is used in the proof of the following theorem: all closed left ideals in a separable C*-algebra, endowed with the standard quantization, are projective left quantum modules over this algebra.


title: Characterizations of operator Birkhoff-James orthogonality

abstract: In this paper, we obtain some characterizations of the (strong) Birkhoff--James orthogonality for elements of Hilbert $C^*$-modules and certain elements of $\mathbb{B}(\mathscr{H})$. Moreover, we obtain a kind of Pythagorean relation for bounded linear operators. In addition, for $T\in \mathbb{B}(\mathscr{H})$ we prove that if the norm attaining set $\mathbb{M}_T$ is a unit sphere of some finite dimensional subspace $\mathscr{H}_0$ of $\mathscr{H}$ and $\|T\|_{{{\mathscr{H}}_0}^\perp} < \|T\|$, then for every $S\in\mathbb{B}(\mathscr{H})$, $T$ is the strong Birkhoff--James orthogonal to $S$ if and only if there exists a unit vector $\xi\in {\mathscr{H}}_0$ such that $\|T\|\xi = |T|\xi$ and $S^*T\xi = 0$. Finally, we introduce a new type of approximate orthogonality and investigate this notion in the setting of inner product $C^*$-modules.


title: Are we seeing accretion flows in a 250kpc-sized Ly-alpha halo at z=3?

abstract: Using MUSE on the ESO-VLT, we obtained a 4 hour exposure of the z=3.12 radio galaxy MRC0316-257. We detect features down to ~10^-19 erg/s/cm^2/arcsec^2 with the highest surface brightness regions reaching more than a factor of 100 higher. We find Ly-alpha emission out to ~250 kpc in projection from the active galactic nucleus (AGN). The emission shows arc-like morphologies arising at 150-250 kpc from the nucleus in projection with the connected filamentary structures reaching down into the circum-nuclear region. The most distant arc is offset by 700 km/s relative to circum-nuclear HeII 1640 emission, which we assume to be at the systemic velocity. As we probe emission closer to the nucleus, the filamentary emission narrows in projection on the sky, the relative velocity decreases to ~250 km/s, and line full-width at half maximum range from 300-700 km/s. From UV line ratios, the emission on scales of 10s of kpc from the nucleus along a wide angle in the direction of the radio jets is clearly excited by the radio jets and ionizing radiation of the AGN. Assuming ionization equilibrium, the more extended emission outside of the axis of the jet direction would require 100% or more illumination to explain the observed surface brightness. High speed (>300 km/s) shocks into rare gas would provide sufficiently high surface brightness. We discuss the possibility that the arcs of Ly-alpha emission represent accretion shocks and the filamentary emission represent gas flows into the halo, and compare our results with gas accretion simulations.


title: Robust doubly charged nodal lines and nodal surfaces in centrosymmetric systems

abstract: Weyl points in three spatial dimensions are characterized by a $\mathbb{Z}$-valued charge -- the Chern number -- which makes them stable against a wide range of perturbations. A set of Weyl points can mutually annihilate only if their net charge vanishes, a property we refer to as robustness. While nodal loops are usually not robust in this sense, it has recently been shown using homotopy arguments that in the centrosymmetric extension of the $\textrm{AI}$ symmetry class they nevertheless develop a $\mathbb{Z}_2$ charge analogous to the Chern number. Nodal loops carrying a non-trivial value of this $\mathbb{Z}_2$ charge are robust, i.e. they can be gapped out only by a pairwise annihilation and not on their own. As this is an additional charge independent of the Berry $\pi$-phase flowing along the band degeneracy, such nodal loops are, in fact, doubly charged. In this manuscript, we generalize the homotopy discussion to the centrosymmetric extensions of all Atland-Zirnbauer classes. We develop a taylored mathematical framework dubbed the AZ+$\mathcal{I}$ classification and show that in three spatial dimensions such robust and multiply charged nodes appear in four of such centrosymmetric extensions, namely AZ+$\mathcal{I}$ classes $\textrm{CI}$ and $\textrm{AI}$ lead to doubly charged nodal lines, while $\textrm{D}$ and $\textrm{BDI}$ support doubly charged nodal surfaces. We remark that no further crystalline symmetries apart from the spatial inversion are necessary for their stability. We provide a description of the corresponding topological charges, and develop simple tight-binding models of various semimetallic and superconducting phases that exhibit these nodes. We also indicate how the concept of robust and multiply charged nodes generalizes to other spatial dimensions.


title: The footprint of cometary dust analogs: I. Laboratory experiments of low-velocity impacts and comparison with Rosetta data

abstract: Cometary dust provides a unique window on dust growth mechanisms during the onset of planet formation. Measurements by the Rosetta spacecraft show that the dust in the coma of comet 67P/Churyumov-Gerasimenko has a granular structure at size scales from sub-um up to several hundreds of um, indicating hierarchical growth took place across these size scales. However, these dust particles may have been modified during their collection by the spacecraft instruments. Here we present the results of laboratory experiments that simulate the impact of dust on the collection surfaces of COSIMA and MIDAS, instruments onboard the Rosetta spacecraft. We map the size and structure of the footprints left by the dust particles as a function of their initial size (up to several hundred um) and velocity (up to 6 m/s). We find that in most collisions, only part of the dust particle is left on the target; velocity is the main driver of the appearance of these deposits. A boundary between sticking/bouncing and fragmentation as an outcome of the particle-target collision is found at v ~ 2 m/s. For velocities below this value, particles either stick and leave a single deposit on the target plate, or bounce, leaving a shallow footprint of monomers. At velocities > 2 m/s and sizes > 80 um, particles fragment upon collision, transferring up to 50 per cent of their mass in a rubble-pile-like deposit on the target plate. The amount of mass transferred increases with the impact velocity. The morphologies of the deposits are qualitatively similar to those found by the COSIMA instrument.


title: Quantifying tidal stream disruption in a simulated Milky Way

abstract: Simulations of tidal streams show that close encounters with dark matter subhalos induce density gaps and distortions in on-sky path along the streams. Accordingly, observing disrupted streams in the Galactic halo would substantiate the hypothesis that dark matter substructure exists there, while in contrast, observing collimated streams with smoothly varying density profiles would place strong upper limits on the number density and mass spectrum of subhalos. Here, we examine several measures of stream "disruption" and their power to distinguish between halo potentials with and without substructure and with different global shapes. We create and evolve a population of 1280 streams on a range of orbits in the Via Lactea II simulation of a Milky Way-like halo, replete with a full mass range of {\Lambda}CDM subhalos, and compare it to two control stream populations evolved in smooth spherical and smooth triaxial potentials, respectively. We find that the number of gaps observed in a stellar stream is a poor indicator of the halo potential, but that (i) the thinness of the stream on-sky, (ii) the symmetry of the leading and trailing tails, and (iii) the deviation of the tails from a low-order polynomial path on-sky ("path regularity") distinguish between the three potentials more effectively. We find that globular cluster streams on low-eccentricity orbits far from the galactic center (apocentric radius ~ 30-80 kpc) are most powerful in distinguishing between the three potentials. If they exist, such streams will shortly be discoverable and mapped in high dimensions with near-future photometric and spectroscopic surveys.


title: LongHCPulse: Long Pulse Heat Capacity on a Quantum Design PPMS

abstract: This paper presents LongHCPulse: software which enables heat capacity to be collected on a Quantum Design PPMS using a long-pulse method. This method, wherein heat capacity is computed from the time derivative of sample temperature over long (30 min) measurement times, is necessary for probing first order transitions and shortens the measurement time by a factor of five. LongHCPulse also includes plotting utilities based on the Matplotlib library. I illustrate the use of LongHCPulse with the example of data taken on ${\rm Yb_{2}Ti_{2}O_{7}}$, and compare the results to the standard semi-adiabatic method.


title: Power-law liquid in cuprate superconductors from fermionic unparticles

abstract: Recent photoemission spectroscopy measurements [arXiv:1509.01611] on cuprate superconductors have inferred that over a wide range of doping, the imaginary part of the electron self-energy scales as $\Sigma^{\prime\prime}\sim(\omega^2+\pi^2T^2)^a$ with $a=1$ in the overdoped Fermi-liquid state and $a<0.5$ in the optimal to underdoped regime. We show that this non-Fermi-liquid scaling behavior can naturally be explained by the presence of a scale-invariant state of matter known as unparticles. We evaluate analytically the electron self-energy due to interactions with fermionic unparticles. We find that, in agreement with experiments, the imaginary part of the self-energy scales with respect to temperature and energy as $\Sigma^{\prime\prime}\sim T^{2+2\alpha}$ and $\omega^{2+2\alpha}$, where $\alpha$ is the anomalous dimension of the unparticle propagator. In addition, the calculated occupancy and susceptibility of fermionic unparticles, unlike those of normal fermions, have significant spectral weights even at high energies. This unconventional behavior is attributed to the branch cut in the unparticle propagator which broadens the unparticle spectral function over a wide energy range and non-trivially alters the scattering phase space by enhancing (suppressing) the intrinsic susceptibility at low energies for negative (positive) $\alpha$. Our work presents new evidence suggesting that unparticles might be important low-energy degrees of freedom in strongly coupled systems such as the cuprate superconductors.


title: Streaming Sparse Gaussian Process Approximations

abstract: Sparse pseudo-point approximations for Gaussian process (GP) models provide a suite of methods that support deployment of GPs in the large data regime and enable analytic intractabilities to be sidestepped. However, the field lacks a principled method to handle streaming data in which both the posterior distribution over function values and the hyperparameter estimates are updated in an online fashion. The small number of existing approaches either use suboptimal hand-crafted heuristics for hyperparameter learning, or suffer from catastrophic forgetting or slow updating when new data arrive. This paper develops a new principled framework for deploying Gaussian process probabilistic models in the streaming setting, providing methods for learning hyperparameters and optimising pseudo-input locations. The proposed framework is assessed using synthetic and real-world datasets.


title: Deep Convolutional Neural Networks as strong gravitational lens detectors

abstract: Future large-scale surveys with high resolution imaging will provide us with a few $10^5$ new strong galaxy-scale lenses. These strong lensing systems however will be contained in large data amounts which are beyond the capacity of human experts to visually classify in a unbiased way. We present a new strong gravitational lens finder based on convolutional neural networks (CNNs). The method was applied to the Strong Lensing challenge organised by the Bologna Lens Factory. It achieved first and third place respectively on the space-based data-set and the ground-based data-set. The goal was to find a fully automated lens finder for ground-based and space-based surveys which minimizes human inspect. We compare the results of our CNN architecture and three new variations ("invariant" "views" and "residual") on the simulated data of the challenge. Each method has been trained separately 5 times on 17 000 simulated images, cross-validated using 3 000 images and then applied to a 100 000 image test set. We used two different metrics for evaluation, the area under the receiver operating characteristic curve (AUC) score and the recall with no false positive ($\mathrm{Recall}_{\mathrm{0FP}}$). For ground based data our best method achieved an AUC score of $0.977$ and a $\mathrm{Recall}_{\mathrm{0FP}}$ of $0.50$. For space-based data our best method achieved an AUC score of $0.940$ and a $\mathrm{Recall}_{\mathrm{0FP}}$ of $0.32$. On space-based data adding dihedral invariance to the CNN architecture diminished the overall score but achieved a higher no contamination recall. We found that using committees of 5 CNNs produce the best recall at zero contamination and consistenly score better AUC than a single CNN. We found that for every variation of our CNN lensfinder, we achieve AUC scores close to $1$ within $6\%$.


title: Giant Planets Can Act As Stabilizing Agents on Debris Disks

abstract: We have explored the evolution of a cold debris disk under the gravitational influence of dwarf planet sized objects (DPs), both in the presence and absence of an interior giant planet. Through detailed long-term numerical simulations, we demonstrate that, when the giant planet is not present, DPs can stir the eccentricities and inclinations of disk particles, in linear proportion to the total mass of the DPs; on the other hand, when the giant planet is included in the simulations, the stirring is approximately proportional to the mass squared. This creates two regimes: below a disk mass threshold (defined by the total mass of DPs), the giant planet acts as a stabilizing agent of the orbits of cometary nucleii, diminishing the effect of the scatterers; above the threshold, the giant contributes to the dispersion of the particles.


title: Jordan frame no-hair for spherical scalar-tensor black holes

abstract: A no-hair theorem for spherical black holes in scalar-tensor gravity is presented. Contrary to the existing theorems, which are proved in the Einstein conformal frame, this proof is performed entirely in the Jordan frame. The theorem is limited to spherical symmetry (instead of axisymmetry) but holds for non-constant Brans-Dicke couplings.


title: Reducing the Top Quark Mass Uncertainty with Jet Grooming

abstract: The measurement of the top quark mass has large systematic uncertainties coming from the Monte Carlo simulations that are used to match theory and experiment. We explore how much that uncertainty can be reduced by using jet grooming procedures. We estimate the inherent ambiguity in what is meant by Monte Carlo mass to be around 530 MeV without any corrections. This uncertainty can be reduced by 60% to 200 MeV by calibrating to the W mass and a further 33% to 140 MeV by applying soft-drop jet grooming (or by 20% more to 170 MeV with trimming). At e+e- colliders, the associated uncertainty is around 110 MeV, reducing to 50 MeV after calibrating to the W mass. By analyzing the tuning parameters, we conclude that the importance of jet grooming after calibrating to the W mass is to reduce sensitivity to the underlying event.


title: Softmax Q-Distribution Estimation for Structured Prediction: A Theoretical Interpretation for RAML

abstract: Reward augmented maximum likelihood (RAML), a simple and effective learning framework to directly optimize towards the reward function in structured prediction tasks, has led to a number of impressive empirical successes. RAML incorporates task-specific reward by performing maximum-likelihood updates on candidate outputs sampled according to an exponentiated payoff distribution, which gives higher probabilities to candidates that are close to the reference output. While RAML is notable for its simplicity, efficiency, and its impressive empirical successes, the theoretical properties of RAML, especially the behavior of the exponentiated payoff distribution, has not been examined thoroughly. In this work, we introduce softmax Q-distribution estimation, a novel theoretical interpretation of RAML, which reveals the relation between RAML and Bayesian decision theory. The softmax Q-distribution can be regarded as a smooth approximation of the Bayes decision boundary, and the Bayes decision rule is achieved by decoding with this Q-distribution. We further show that RAML is equivalent to approximately estimating the softmax Q-distribution, with the temperature $\tau$ controlling approximation error. We perform two experiments, one on synthetic data of multi-class classification and one on real data of image captioning, to demonstrate the relationship between RAML and the proposed softmax Q-distribution estimation method, verifying our theoretical analysis. Additional experiments on three structured prediction tasks with rewards defined on sequential (named entity recognition), tree-based (dependency parsing) and irregular (machine translation) structures show notable improvements over maximum likelihood baselines.


title: Deep De-Aliasing for Fast Compressive Sensing MRI

abstract: Fast Magnetic Resonance Imaging (MRI) is highly in demand for many clinical applications in order to reduce the scanning cost and improve the patient experience. This can also potentially increase the image quality by reducing the motion artefacts and contrast washout. However, once an image field of view and the desired resolution are chosen, the minimum scanning time is normally determined by the requirement of acquiring sufficient raw data to meet the Nyquist-Shannon sampling criteria. Compressive Sensing (CS) theory has been perfectly matched to the MRI scanning sequence design with much less required raw data for the image reconstruction. Inspired by recent advances in deep learning for solving various inverse problems, we propose a conditional Generative Adversarial Networks-based deep learning framework for de-aliasing and reconstructing MRI images from highly undersampled data with great promise to accelerate the data acquisition process. By coupling an innovative content loss with the adversarial loss our de-aliasing results are more realistic. Furthermore, we propose a refinement learning procedure for training the generator network, which can stabilise the training with fast convergence and less parameter tuning. We demonstrate that the proposed framework outperforms state-of-the-art CS-MRI methods, in terms of reconstruction error and perceptual image quality. In addition, our method can reconstruct each image in 0.22ms--0.37ms, which is promising for real-time applications.


title: Exoplanet Predictions Based on Harmonic Orbit Resonances

abstract: The current exoplanet database includes 5454 confirmed planets and candidate planets observed with the KEPLER mission. We find 932 planet pairs from which we extract distance and orbital period ratios. While earlier studies used the Titius-Bode law or a generalized version with logarithmic spacing, which both lack a physical model, we employ here the theory of harmonic orbit resonances, which contains quantized ratios instead, to explain the observed planet distance ratios and to predict undetected exoplanets. We find that the most prevailing harmonic ratios are (2:1), (3:2), and (5:3), in 73\% of the cases, while alternative harmonic ratios of (5:4), (4:3), (5:2), (3:1) occur in 27\% of the other cases. Our orbital predictions includes 171 exoplanets, 2 Jupiter moons, one Saturn moon, 3 Uranus moons, and 4 Neptune moons. The accuracy of the predicted planet distances amounts to a few percent, which fits the data significantly better than the Titius-Bode law or a logarithmic spacing. This information may be useful for targeted exoplanet searches with Kepler data and to estimate the number of live-carrying planets in habitable zones.


title: Asymmetry and non-dispersivity in the Aharonov-Bohm effect

abstract: Decades ago, Aharonov and Bohm showed that electrons are affected by electromagnetic potentials in the absence of forces due to fields. Zeilinger's theorem describes this absence of classical force in quantum terms as the "dispersionless" nature of the Aharonov-Bohm effect. Shelankov predicted the presence of a quantum "force" for the same Aharonov-Bohm physical system as elucidated by Berry. Here, we report an experiment designed to test Shelankov's prediction and we provide a theoretical analysis that is intended to elucidate the relation between Shelankov's prediction and Zeilinger's theorem. The experiment consists of the Aharonov--Bohm physical system; free electrons pass a magnetized nanorod and far--field electron diffraction is observed. The diffraction pattern is asymmetric confirming one of Shelankov's predictions and giving indirect experimental evidence for the presence of a quantum "force". Our theoretical analysis shows that Zeilinger's theorem and Shelankov's result are both special cases of one theorem.


title: A Fast Frequent Directions Algorithm for Low Rank Approximation

abstract: Recently a deterministic method, frequent directions (FD) is proposed to solve the high dimensional low rank approximation problem. It works well in practice, but experiences high computational cost. In this paper, we establish a fast frequent directions algorithm for the low rank approximation problem, which implants a randomized algorithm, sparse subspace embedding (SpEmb) in FD. This new algorithm makes use of FD's natural block structure and sends more information through SpEmb to each block in FD. We prove that our new algorithm produces a good low rank approximation with a sketch of size linear on the rank approximated. Its effectiveness and efficiency are demonstrated by the experimental results on both synthetic and real world datasets, as well as applications in network analysis.


title: Coboundary categories and local rules

abstract: First we develop the theory of local rules for coboundary categories. Then we describe the local rules in two main cases. First for the quantum groups in general and in the seminormal representations of the Hecke algebras. Then for crystals in general and specifically for crystals of minuscule representations. Finally we show how growth diagrams can be extended to construct the action of the cactus group on highest weight words.


title: Simultaneous Multiple Surface Segmentation Using Deep Learning

abstract: The task of automatically segmenting 3-D surfaces representing boundaries of objects is important for quantitative analysis of volumetric images, and plays a vital role in biomedical image analysis. Recently, graph-based methods with a global optimization property have been developed and optimized for various medical imaging applications. Despite their widespread use, these require human experts to design transformations, image features, surface smoothness priors, and re-design for a different tissue, organ or imaging modality. Here, we propose a Deep Learning based approach for segmentation of the surfaces in volumetric medical images, by learning the essential features and transformations from training data, without any human expert intervention. We employ a regional approach to learn the local surface profiles. The proposed approach was evaluated on simultaneous intraretinal layer segmentation of optical coherence tomography (OCT) images of normal retinas and retinas affected by age related macular degeneration (AMD). The proposed approach was validated on 40 retina OCT volumes including 20 normal and 20 AMD subjects. The experiments showed statistically significant improvement in accuracy for our approach compared to state-of-the-art graph based optimal surface segmentation with convex priors (G-OSC). A single Convolution Neural Network (CNN) was used to learn the surfaces for both normal and diseased images. The mean unsigned surface positioning errors obtained by G-OSC method 2.31 voxels (95% CI 2.02-2.60 voxels) was improved to $1.27$ voxels (95% CI 1.14-1.40 voxels) using our new approach. On average, our approach takes 94.34 s, requiring 95.35 MB memory, which is much faster than the 2837.46 s and 6.87 GB memory required by the G-OSC method on the same computer system.


title: A New 3D Segmentation Methodology for Lumbar Vertebral Bodies for the Measurement of BMD and Geometry

abstract: In this paper a new technique is presented that extracts the geometry of lumbar vertebral bodies from spiral CT scans. Our new multi-step segmentation approach yields highly accurate and precise measurement of the bone mineral density (BMD) in different volumes of interest which are defined relative to a local anatomical coordinate systems. The approach also enables the analysis of the geometry of the relevant vertebrae. Intra- and inter operator precision for segmentation, BMD measurement and position of the coordinate system are below 1.5% in patient data, accuracy errors are below 1.5% for BMD and below 4% for volume in phantom data. The long-term goal of the approach is to improve fracture prediction in osteoporosis.


title: Sparse Coding on Stereo Video for Object Detection

abstract: Deep Convolutional Neural Networks (DCNN) require millions of labeled training examples for image classification and object detection tasks, which restrict these models to domains where such datasets are available. In this paper, we explore the use of unsupervised sparse coding applied to stereo-video data to help alleviate the need for large amounts of labeled data. We show that replacing a typical supervised convolutional layer with an unsupervised sparse-coding layer within a DCNN allows for better performance on a car detection task when only a limited number of labeled training examples is available. Furthermore, the network that incorporates sparse coding allows for more consistent performance over varying initializations and ordering of training examples when compared to a fully supervised DCNN. Finally, we compare activations between the unsupervised sparse-coding layer and the supervised convolutional layer, and show that the sparse representation exhibits an encoding that is depth selective, whereas encodings from the convolutional layer do not exhibit such selectivity. These result indicates promise for using unsupervised sparse-coding approaches in real-world computer vision tasks in domains with limited labeled training data.


title: Duality for spherical representations in exceptional theta correspondences

abstract: We study the exceptional theta correspondence for real groups obtained by restricting the minimal representation of the split exceptional group of the type E_n, to a split dual pair where one member is the exceptional group of the type G_2. We prove that the correspondence gives a bijection between spherical representations if n=6,7, and a slightly weaker statement if n=8.


title: A New 3D Method to Segment the Lumbar Vertebral Bodies and to Determine Bone Mineral Density and Geometry

abstract: In this paper we present a new 3D segmentation approach for the vertebrae of the lower thoracic and the lumbar spine in spiral computed tomography datasets. We implemented a multi-step procedure. Its main components are deformable models, volume growing, and morphological operations. The performance analysis that included an evaluation of accuracy using the European Spine Phantom, and of intra-operator precision using clinical CT datasets from 10 patients highlight the potential for clinical use. The intra-operator precision of the segmentation procedure was better than 1% for Bone Mineral Density (BMD) and better than 1.8% for volume. The long-term goal of this work is to enable better fracture prediction and improved patient monitoring in the field of osteoporosis. A true 3D segmentation also enables an accurate measurement of geometrical parameters that can augment the classical measurement of BMD.


title: Quaternion algebra for Stokes-Mueller formalism

abstract: It is shown that the Stokes-Mueller formalism can be reformulated in terms of quaternions, and the quaternion approach is more suitable for the formalism of Mueller-Jones states that we have recently described. In terms of quaternions it can be shown that the vector and matrix states and the Jones matrix associated to nondepolarizing optical systems are different representations isomorphic to the same quaternion state, and this quaternion state turns out to be the rotator of the Stokes quaternion. It is also shown that the coherent linear combination of nondepolarizing optical media states and depolarization phenomena can be reformulated in terms of quaternion states.


title: Ballistic near-field heat transport in dense many-body systems

abstract: Radiative heat-transport mediated by near-field interactions is known to be superdiffusive in dilute, many-body systems. In this Letter we use a generalized Landauer theory of radiative heat transfer in many-body planar systems to demonstrate a nonmonotonic transition from superdiffusive to ballistic transport in dense systems. We show that such a transition is associated to a change of the polarization of dominant modes, leading to dramatically different thermal relaxation dynamics spanning over three orders of magnitude. This result could have important consequences on thermal management at nanoscale of many-body systems.


title: Local Information with Feedback Perturbation Suffices for Dictionary Learning in Neural Circuits

abstract: While the sparse coding principle can successfully model information processing in sensory neural systems, it remains unclear how learning can be accomplished under neural architectural constraints. Feasible learning rules must rely solely on synaptically local information in order to be implemented on spatially distributed neurons. We describe a neural network with spiking neurons that can address the aforementioned fundamental challenge and solve the L1-minimizing dictionary learning problem, representing the first model able to do so. Our major innovation is to introduce feedback synapses to create a pathway to turn the seemingly non-local information into local ones. The resulting network encodes the error signal needed for learning as the change of network steady states caused by feedback, and operates akin to the classical stochastic gradient descent method.


title: State Complexity of Reversals of Deterministic Finite Automata with Output

abstract: We investigate the worst-case state complexity of reversals of deterministic finite automata with output (DFAOs). In these automata, each state is assigned some output value, rather than simply being labelled final or non-final. This directly generalizes the well-studied problem of determining the worst-case state complexity of reversals of ordinary deterministic finite automata. If a DFAO has $n$ states and $k$ possible output values, there is a known upper bound of $k^n$ for the state complexity of reversal. We show this bound can be reached with a ternary input alphabet. We conjecture it cannot be reached with a binary input alphabet except when $k = 2$, and give a lower bound for the case $3 \le k < n$. We prove that the state complexity of reversal depends solely on the transition monoid of the DFAO and the mapping that assigns output values to states.


title: The c=1 String Theory S-Matrix Revisited

abstract: We revisit the perturbative S-matrix of c=1 string theory from the worldsheet perspective. We clarify the origin of the leg pole factors, the non-analyticity of the string amplitudes, and the validity as well as limitations of earlier computations based on resonance momenta. We compute the tree level 4-point amplitude and the genus one 2-point reflection amplitude by numerically integrating Virasoro conformal blocks with DOZZ structure constants on the sphere and on the torus, with sufficiently generic complex Liouville momenta, and find agreement with known answers from the c=1 matrix model.


title: Data-driven Optimal Cost Selection for Distributionally Robust Optimization

abstract: Recently, (Blanchet, Kang, and Murhy 2016, and Blanchet, and Kang 2017) showed that several machine learning algorithms, such as square-root Lasso, Support Vector Machines, and regularized logistic regression, among many others, can be represented exactly as distributionally robust optimization (DRO) problems. The distributional uncertainty is defined as a neighborhood centered at the empirical distribution. We propose a methodology which learns such neighborhood in a natural data-driven way. We show rigorously that our framework encompasses adaptive regularization as a particular case. Moreover, we demonstrate empirically that our proposed methodology is able to improve upon a wide range of popular machine learning estimators.


title: Exciton spectrum in two-dimensional transition metal dichalcogenides: The role of Diracness

abstract: The physics of excitons, electron-hole pairs that are bound together by their mutual Coulomb attraction, can to great extent be understood in the framework of the quantum-mechanical hydrogen model. This model has recently been challenged by spectroscopic measurements on two-dimensional transition-metal dichalchogenides that unveil strong deviations from a hydrogenic spectrum. Here, we show that this deviation is due to the particular relativistic character of electrons in this class of materials. Indeed, their electrons are no longer described in terms of a Schroedinger but a massive Dirac equation that intimately links electrons to holes. Dirac excitons therefore inherit a relativistic quantum spin-1/2 that contributes to the angular momentum and thus the exciton spectrum. Most saliently, the level spacing is strongly reduced as compared to the hydrogen model, in agreement with spectroscopic measurements and ab-initio calculations.


title: Demonstration of a quantum key distribution network in urban fibre-optic communication lines

abstract: We report the results of the implementation of a quantum key distribution (QKD) network using standard fibre communication lines in Moscow. The developed QKD network is based on the paradigm of trusted repeaters and allows a common secret key to be generated between users via an intermediate trusted node. The main feature of the network is the integration of the setups using two types of encoding, i.e. polarisation encoding and phase encoding. One of the possible applications of the developed QKD network is the continuous key renewal in existing symmetric encryption devices with a key refresh time of up to 14 s.


title: Compressing Over-the-Counter Markets

abstract: Over-the-counter markets are at the center of the postcrisis global reform of the financial system. We show how the size and structure of such markets can undergo rapid and extensive changes when participants engage in portfolio compression, a post-trade netting technology. Tightly-knit and concentrated trading structures, as featured by many large over-the-counter markets, are especially susceptible to reductions of notional and reconfigurations of network structure resulting from compression activities. Using transaction-level data on credit-default-swaps markets, we estimate reduction levels consistent with the historical development observed in these markets since the Global Financial Crisis. Finally, we study the effect of a mandate to centrally clear over-the-counter markets. When participants engage in both central clearing and portfolio compression, we find large netting failures if clearinghouses proliferate. Allowing for compression across clearinghouses by-and-large offsets this adverse effect.


title: Exotic looped trajectories in double-slit experiments with matter waves

abstract: We study the observation of exotic looped trajectories in double-slit experiments with matter waves. We consider the relative intensity at $x=0$ as a function of the time-of-flight from the double-slit to the screen inside the interferometer. This allows us to define a fringe visibility associated to the contribution to the interference pattern given by exotic lopped trajectories. We demonstrate that the Sorkin parameter is given in terms of this visibility and of the axial phases which include the Gouy phase. We verify how this parameter can be obtained by measuring the relative intensity at the screen. We show that the effect of exotic looped trajectories can be significantly increased by simply adjusting the parameters of the double-slit apparatus. Applying our results to the case of neutron interferometry, we obtain a maximum Sorkin parameter of the order of $|\kappa_{max}|\approx0.2$, which is the value of the fringe visibility.


title: Clustering under Local Stability: Bridging the Gap between Worst-Case and Beyond Worst-Case Analysis

abstract: Recently, there has been substantial interest in clustering research that takes a beyond worst-case approach to the analysis of algorithms. The typical idea is to design a clustering algorithm that outputs a near-optimal solution, provided the data satisfy a natural stability notion. For example, Bilu and Linial (2010) and Awasthi et al. (2012) presented algorithms that output near-optimal solutions, assuming the optimal solution is preserved under small perturbations to the input distances. A drawback to this approach is that the algorithms are often explicitly built according to the stability assumption and give no guarantees in the worst case; indeed, several recent algorithms output arbitrarily bad solutions even when just a small section of the data does not satisfy the given stability notion. In this work, we address this concern in two ways. First, we provide algorithms that inherit the worst-case guarantees of clustering approximation algorithms, while simultaneously guaranteeing near-optimal solutions when the data is stable. Our algorithms are natural modifications to existing state-of-the-art approximation algorithms. Second, we initiate the study of local stability, which is a property of a single optimal cluster rather than an entire optimal solution. We show our algorithms output all optimal clusters which satisfy stability locally. Specifically, we achieve strong positive results in our local framework under recent stability notions including metric perturbation resilience (Angelidakis et al. 2017) and robust perturbation resilience (Balcan and Liang 2012) for the $k$-median, $k$-means, and symmetric/asymmetric $k$-center objectives.


title: Improved Very-short-term Wind Forecasting using Atmospheric Classification

abstract: We present a regime-switching vector-autoregressive method for very-short-term wind speed forecasting at multiple locations with regimes based on large-scale meteorological phenomena. Statistical methods short-term wind forecasting out-perform numerical weather prediction for forecast horizons up to a few hours, and the spatio-temporal interdependency between geographically dispersed locations may be exploited to improve forecast skill. Here we show that conditioning spatio-temporal interdependency on `atmospheric modes' can further improve forecast performance. The modes are defined from the atmospheric classification of wind and pressure fields at the surface level, and the geopotential height field at the 500hPa level. The data fields are extracted from the MERRA-2 reanalysis dataset with an hourly temporal resolution over the UK, atmospheric patterns are classified using self-organising maps and then clustered to optimise forecast performance. In a case study based on 6 years of measurements from 23 weather stations in the UK, a set of three atmospheric modes are found to be optimal for forecast performance. The skill in the one- to six-hour-ahead forecasts is improved at all sites compared to persistence and competitive benchmarks. Across the 23 test sites, one-hour-ahead root mean squared error is reduced by between 0.3% and 4.1% compared to the best performing benchmark, and by and average of 1.6% over all sites; the six-hour-ahead accuracy is improved by an average of 3.1%.


title: Multiphoton Effects Enhanced Due to Ultrafast Photon-Number Fluctuations

abstract: Multi-photon processes are the essence of nonlinear optics. Optical harmonics generation and multi-photon absorption, ionization, polymerization or spectroscopy are widely used in practical applications. Generally, the rate of an n-photon effect scales as the n-th order autocorrelation function of the incident light, which is high for light with strong photon-number fluctuations. Therefore `noisy' light sources are much more efficient for multi-photon effects than coherent sources with the same mean power, pulse duration and repetition rate. Here we generate optical harmonics of order 2-4 from bright squeezed vacuum (BSV), a state of light consisting of only quantum noise with no coherent component. We observe up to two orders of magnitude enhancement in the generation of optical harmonics due to ultrafast photon-number fluctuations. This feature is especially important for the nonlinear optics of fragile structures where the use of a `noisy' pump can considerably increase the effect without overcoming the damage threshold.


title: Theoretical and computational aspects of entanglement

abstract: We show that the two notions of entanglement: the maximum of the geometric measure of entanglement and the maximum of the nuclear norm is attained for the same states. We affirm the conjecture of Higuchi-Sudberry on the maximum entangled state of four qubits. We introduce the notion of d-density tensor for mixed d-partite states. We show that d-density tensor is separable if and only if its nuclear norm is $1$. We suggest an alternating method for computing the nuclear norm of tensors. We apply the above results to symmetric tensors. We give many numerical examples.


title: Statistical physics of human cooperation

abstract: Extensive cooperation among unrelated individuals is unique to humans, who often sacrifice personal benefits for the common good and work together to achieve what they are unable to execute alone. The evolutionary success of our species is indeed due, to a large degree, to our unparalleled other-regarding abilities. Yet, a comprehensive understanding of human cooperation remains a formidable challenge. Recent research in social science indicates that it is important to focus on the collective behavior that emerges as the result of the interactions among individuals, groups, and even societies. Non-equilibrium statistical physics, in particular Monte Carlo methods and the theory of collective behavior of interacting particles near phase transition points, has proven to be very valuable for understanding counterintuitive evolutionary outcomes. By studying models of human cooperation as classical spin models, a physicist can draw on familiar settings from statistical physics. However, unlike pairwise interactions among particles that typically govern solid-state physics systems, interactions among humans often involve group interactions, and they also involve a larger number of possible states even for the most simplified description of reality. The complexity of solutions therefore often surpasses that observed in physical systems. Here we review experimental and theoretical research that advances our understanding of human cooperation, focusing on spatial pattern formation, on the spatiotemporal dynamics of observed solutions, and on self-organization that may either promote or hinder socially favorable states.

categories: physics.soc-ph cond-mat.stat-mech cs.SI nlin.AO q-bio.PE


title: A Lightweight Approach for On-the-Fly Reflectance Estimation

abstract: Estimating surface reflectance (BRDF) is one key component for complete 3D scene capture, with wide applications in virtual reality, augmented reality, and human computer interaction. Prior work is either limited to controlled environments (\eg gonioreflectometers, light stages, or multi-camera domes), or requires the joint optimization of shape, illumination, and reflectance, which is often computationally too expensive (\eg hours of running time) for real-time applications. Moreover, most prior work requires HDR images as input which further complicates the capture process. In this paper, we propose a lightweight approach for surface reflectance estimation directly from $8$-bit RGB images in real-time, which can be easily plugged into any 3D scanning-and-fusion system with a commodity RGBD sensor. Our method is learning-based, with an inference time of less than 90ms per scene and a model size of less than 340K bytes. We propose two novel network architectures, HemiCNN and Grouplet, to deal with the unstructured input data from multiple viewpoints under unknown illumination. We further design a loss function to resolve the color-constancy and scale ambiguity. In addition, we have created a large synthetic dataset, SynBRDF, which comprises a total of $500$K RGBD images rendered with a physically-based ray tracer under a variety of natural illumination, covering $5000$ materials and $5000$ shapes. SynBRDF is the first large-scale benchmark dataset for reflectance estimation. Experiments on both synthetic data and real data show that the proposed method effectively recovers surface reflectance, and outperforms prior work for reflectance estimation in uncontrolled environments.


title: Towards Communication-Aware Robust Topologies

abstract: We currently witness the emergence of interesting new network topologies optimized towards the traffic matrices they serve, such as demand-aware datacenter interconnects (e.g., ProjecToR) and demand-aware overlay networks (e.g., SplayNets). This paper introduces a formal framework and approach to reason about and design such topologies. We leverage a connection between the communication frequency of two nodes and the path length between them in the network, which depends on the entropy of the communication matrix. Our main contribution is a novel robust, yet sparse, family of network topologies which guarantee an expected path length that is proportional to the entropy of the communication patterns.


title: Relaxed Wasserstein with Applications to GANs

abstract: Wasserstein Generative Adversarial Networks (WGANs) provide a versatile class of models, which have attracted great attention in various applications. However, this framework has two main drawbacks: (i) Wasserstein-1 (or Earth-Mover) distance is restrictive such that WGANs cannot always fit data geometry well; (ii) It is difficult to achieve fast training of WGANs. In this paper, we propose a new class of \textit{Relaxed Wasserstein} (RW) distances by generalizing Wasserstein-1 distance with Bregman cost functions. We show that RW distances achieve nice statistical properties while not sacrificing the computational tractability. Combined with the GANs framework, we develop Relaxed WGANs (RWGANs) which are not only statistically flexible but can be approximated efficiently using heuristic approaches. Experiments on real images demonstrate that the RWGAN with Kullback-Leibler (KL) cost function outperforms other competing approaches, e.g., WGANs, even with gradient penalty.


title: Detection of the Ultraviolet Spectrum of the Hot Subdwarf Companion of 60 Cygni (B1 Ve) from a Survey of IUE Spectra of Be Stars

abstract: We used archival International Ultraviolet Explorer (IUE) high-dispersion, short wavelength spectra data to search for evidence of the spectra of hot subdwarf companions of six rapidly rotating Be stars in binary systems. We searched for the signature of a hot companion through an analysis of the cross-correlation functions of observed and model spectra that were separated into primary and secondary components using a Doppler tomography algorithm and adopted spectroscopic orbital solutions. A positive detection of the flux from a hot companion was made for the reconstructed secondary cross-correlation function of just one target, 60 Cygni (B1 Ve). We estimate that the companion of the Be star in 60 Cygni has $T_{\rm eff} = 42 \pm 4$ kK, mass ratio $M_2/M_1 = 0.15 \pm 0.02$, and monochromatic flux ratio ${f_2}/{f_1} = 0.034 \pm 0.002$ in the spectral region near 1525 \AA. If the companions of the other target Be stars are also hot, then they must be faint and contribute less than $\approx 1\%$ of the UV flux ($<0.6 \%$ in the case of $\gamma$ Cas). We also discuss in an appendix a shell episode of Pleione (28 Tau) recorded in the IUE spectra.


title: Competing Bose-Glass physics with disorder-induced Bose-Einstein condensation in the doped $S=1$ antiferromagnet Ni(Cl$_{1-x}$Br$_x$)$_2$-4SC(NH$_2$)$_2$ at high magnetic fields

abstract: We study the interplay between disorder and interactions for emergent bosonic degrees of freedom induced by an external magnetic field in the Br-doped spin-gapped antiferromagnetic material Ni(Cl$_{1-x}$Br$_x$)$_2$-4SC(NH$_2$)$_2$ (DTNX). Building on nuclear magnetic resonance experiments at high magnetic field [A. Orlova et al., Phys. Rev. Lett. 118, 067203 (2017)], we describe the localization of isolated impurity states, providing a realistic theoretical modeling for DTNX. Going beyond single impurity localization we use quantum Monte Carlo simulations to explore many-body effects from which pairwise effective interactions lead to a (impurity-induced) BEC revival [M. Dupont, S. Capponi, and N. Laflorencie, Phys. Rev. Lett. 118, 067204 (2017)]. We further address the question of the existence of a many-body localized Bose-glass (BG) phase in DTNX, which is found to compete with a series of a new kind of BEC regimes made out of the multi-impurity states. The global magnetic field-temperature phase diagram of DTNX reveals a very rich structure for low impurity concentration, with consecutive disorder-induced BEC mini-domes separated by intervening many-body localized BG regimes. Upon increasing the impurity level, multiple mini-BEC phases start to overlap, while intermediate BG regions vanish.


title: Connecting Different TMD Factorization Formalisms in QCD

abstract: In the original Collins-Soper-Sterman (CSS) presentation of the results of transverse-momentum-dependent (TMD) factorization for the Drell-Yan process, results for perturbative coefficients can be obtained from calculations for collinear factorization. Here we show how to use these results, plus known results for the quark form factor, to obtain coefficients for TMD factorization in more recent formulations, e.g., that due to Collins, and apply them to known results at order $\alpha_s^2$ and $\alpha_s^3$. We also show that the "non-perturbative" functions as obtained from fits to data are equal in the two schemes. We compile the higher-order perturbative inputs needed for the updated CSS scheme by appealing to results obtained in a variety of different formalisms. In addition, we derive the connection between both versions of the CSS formalism and several formalisms based in soft-collinear effective theory (SCET). Our work uses some important new results for factorization for the quark form factor, which we derive.


title: Doubly Robust Data-Driven Distributionally Robust Optimization

abstract: Data-driven Distributionally Robust Optimization (DD-DRO) via optimal transport has been shown to encompass a wide range of popular machine learning algorithms. The distributional uncertainty size is often shown to correspond to the regularization parameter. The type of regularization (e.g. the norm used to regularize) corresponds to the shape of the distributional uncertainty. We propose a data-driven robust optimization methodology to inform the transportation cost underlying the definition of the distributional uncertainty. We show empirically that this additional layer of robustification, which produces a method we called doubly robust data-driven distributionally robust optimization (DD-R-DRO), allows to enhance the generalization properties of regularized estimators while reducing testing error relative to state-of-the-art classifiers in a wide range of data sets.


title: Hall-Littlewood RSK field

abstract: We introduce a randomized Hall-Littlewood RSK algorithm and study its combinatorial and probabilistic properties. On the probabilistic side, a new model --- the Hall-Littlewood RSK field --- is introduced. Its various degenerations contain known objects (the stochastic six vertex model, the asymmetric simple exclusion process) as well as a variety of new ones. We provide formulas for a rich class of observables of these models, extending existing results about Macdonald processes. On the combinatorial side, we establish analogs of properties of the classical RSK algorithm: invertibility, symmetry, and a "bijectivization" of the skew-Cauchy identity.


title: Shadow wave function with a symmetric kernel

abstract: A shadow wave function with an explicit symmetric kernel is introduced. As a consequence the atoms exchange in the system is enhanced. Basic properties of this class of trial functions are kept and quantities it can describe are easily estimated. The effectiveness of this approach is analized by computing properties of interest in a system formed from 4He atoms.


title: Nestrov's Acceleration For Second Order Method

abstract: Optimization plays a key role in machine learning. Recently, stochastic second-order methods have attracted much attention due to their low computational cost in each iteration. However, these algorithms might perform poorly especially if it is hard to approximate the Hessian well and efficiently. As far as we know, there is no effective way to handle this problem. In this paper, we resort to Nestrov's acceleration technique to improve the convergence performance of a class of second-order methods called approximate Newton. We give a theoretical analysis that Nestrov's acceleration technique can improve the convergence performance for approximate Newton just like for first-order methods. We accordingly propose an accelerated regularized sub-sampled Newton. Our accelerated algorithm performs much better than the original regularized sub-sampled Newton in experiments, which validates our theory empirically. Besides, the accelerated regularized sub-sampled Newton has good performance comparable to or even better than state-of-art algorithms.


title: A New Photomechanical Molecular Switch Based on a Linear {\pi}-Conjugated System

abstract: We report the electronic transport properties of a new photo-addressable molecular switch. The switching process relies on a new concept based on linear {\pi}-conjugated dynamic systems, in which the geometry and hence the electronic properties of an oligothiophene chain can be reversibly modified by the photochemical trans-cis isomerization of an azobenzene unit fixed in a lateral loop. Electron transport measurements through self-assembled monolayers on gold, contacted with eGaIn top contact, show switching with a conductance ratio up to 1E3. Ab initio calculations have been used to identify the most energetically stable conformations of the molecular switch, the corresponding calculated conductances qualitatively explain the trend observed in the photo-switching experiments.


title: On Irregular Singularity Wave Functions and Superconformal Indices

abstract: We generalize, in a manifestly Weyl-invariant way, our previous expressions for irregular singularity wave functions in two-dimensional SU(2) q-deformed Yang-Mills theory to SU(N). As an application, we give closed-form expressions for the Schur indices of all (A_{N-1}, A_{N(n-1)-1}) Argyres-Douglas (AD) superconformal field theories (SCFTs), thus completing the computation of these quantities for the (A_N, A_M) SCFTs. With minimal effort, our wave functions also give new Schur indices of various infinite sets of "Type IV" AD theories. We explore the discrete symmetries of these indices and also show how highly intricate renormalization group (RG) flows from isolated theories and conformal manifolds in the ultraviolet to isolated theories and (products of) conformal manifolds in the infrared are encoded in these indices. We compare our flows with dimensionally reduced flows via a simple "monopole vev RG" formalism. Finally, since our expressions are given in terms of concise Lie algebra data, we speculate on extensions of our results that might be useful for probing the existence of hypothetical SCFTs based on other Lie algebras. We conclude with a discussion of some open problems.


title: Impurity-induced magnetic droplet in unconventional superconductors near magnetic instability: Application to Nd doped $\mathbf{CeCoIn_5}$

abstract: Beside the spin density wave (SDW) order inside the superconducting phase in $\mathrm{CeCoIn_5}$ at high magnetic fields, recent neutron scattering measurements have found Bragg peaks in $5\%$ Nd doped $\mathrm{CeCoIn_5}$ at low fields. The intensity of Bragg peaks in low fields is suppressed by increasing field. Based on the phenomenological and microscopic modeling, we show that for the Pauli limited $d$-wave superconductors in the vicinity of SDW instability relevant for $\mathrm{CeCoIn_5}$, magnetic impurities locally induce droplets of SDW order. Because of the strong anisotropy in the momentum space in the spin fluctuations guaranteed by the $d$-wave pairing symmetry, sharp peaks in spin structure factor at $\mathbf{Q}$s are produced by the impurities, even when the droplets of SDW do not order. At zero field, the Nd impurity spins are along the $c$ axis due to the coupling to the conduction electrons with an easy $c$ axis, besides their own crystal field effect. The in-plane magnetic field cants the impurity moments toward the $ab$ plane, which suppress the droplets of SDW order. At high fields, the long-range SDW inside the superconducting phase is stabilized as a consequence of magnon condensation. Our results are consistent with the recent neutron scattering and thermal conductivity measurements.


title: Espresso: Efficient Forward Propagation for BCNNs

abstract: There are many applications scenarios for which the computational performance and memory footprint of the prediction phase of Deep Neural Networks (DNNs) needs to be optimized. Binary Neural Networks (BDNNs) have been shown to be an effective way of achieving this objective. In this paper, we show how Convolutional Neural Networks (CNNs) can be implemented using binary representations. Espresso is a compact, yet powerful library written in C/CUDA that features all the functionalities required for the forward propagation of CNNs, in a binary file less than 400KB, without any external dependencies. Although it is mainly designed to take advantage of massive GPU parallelism, Espresso also provides an equivalent CPU implementation for CNNs. Espresso provides special convolutional and dense layers for BCNNs, leveraging bit-packing and bit-wise computations for efficient execution. These techniques provide a speed-up of matrix-multiplication routines, and at the same time, reduce memory usage when storing parameters and activations. We experimentally show that Espresso is significantly faster than existing implementations of optimized binary neural networks ($\approx$ 2 orders of magnitude). Espresso is released under the Apache 2.0 license and is available at http://github.com/fpeder/espresso.


abstract: This paper considers the distributed optimization problem over a network, where the objective is to optimize a global function formed by a sum of local functions, using only local computation and communication. We develop an Accelerated Distributed Nesterov Gradient Descent (Acc-DNGD) method. When the objective function is convex and $L$-smooth, we show that it achieves a $O(\frac{1}{t^{1.4-\epsilon}})$ convergence rate for all $\epsilon\in(0,1.4)$. We also show the convergence rate can be improved to $O(\frac{1}{t^2})$ if the objective function is a composition of a linear map and a strongly-convex and smooth function. When the objective function is $\mu$-strongly convex and $L$-smooth, we show that it achieves a linear convergence rate of $O([ 1 - C (\frac{\mu}{L})^{5/7} ]^t)$, where $\frac{L}{\mu}$ is the condition number of the objective, and $C>0$ is some constant that does not depend on $\frac{L}{\mu}$.


title: Model-Based Planning with Discrete and Continuous Actions

abstract: Action planning using learned and differentiable forward models of the world is a general approach which has a number of desirable properties, including improved sample complexity over model-free RL methods, reuse of learned models across different tasks, and the ability to perform efficient gradient-based optimization in continuous action spaces. However, this approach does not apply straightforwardly when the action space is discrete. In this work, we show that it is in fact possible to effectively perform planning via backprop in discrete action spaces, using a simple paramaterization of the actions vectors on the simplex combined with input noise when training the forward model. Our experiments show that this approach can match or outperform model-free RL and discrete planning methods on gridworld navigation tasks in terms of performance and/or planning time while using limited environment interactions, and can additionally be used to perform model-based control in a challenging new task where the action space combines discrete and continuous actions. We furthermore propose a policy distillation approach which yields a fast policy network which can be used at inference time, removing the need for an iterative planning procedure.


title: Accelerated Parallel Non-conjugate Sampling for Bayesian Non-parametric Models

abstract: Inference of latent feature models in the Bayesian nonparametric setting is generally difficult, especially in high dimensional settings, because it usually requires proposing features from some prior distribution. In special cases, where the integration is tractable, we can sample new feature assignments according to a predictive likelihood. We present a novel method to accelerate the mixing of latent variable model inference by proposing feature locations based on the data, as opposed to the prior. First, we introduce an accelerated feature proposal mechanism that we show is a valid MCMC algorithm for posterior inference. Next, we propose an approximate inference strategy to perform accelerated inference in parallel. A two-stage algorithm that combines the two approaches provides a computationally attractive method that can quickly reach local convergence to the posterior distribution of our model, while allowing us to exploit parallelization.


title: Invariants of maximal tori and unipotent constituents of some quasi-projective characters for finite classical groups

abstract: We study the decomposition of certain reducible characters of classical groups as the sum of irreducible ones. Let ${\mathbf G}$ be an algebraic group of classical type with defining characteristic $p>0$, $\mu$ a dominant weight and $W$ the Weyl group of ${\mathbf G}$. Let $G=G(q)$ be a finite classical group, where $q$ is a $p$-power. For a weight $\mu$ of ${\mathbf G}$ the sum $s_\mu$ of distinct weights $w(\mu)$ with $w\in W$ viewed as a function on the semisimple elements of $G$ is known to be a generalized Brauer character of $G$ called an orbit character of $G$. We compute, for certain orbit characters and every maximal torus $T$ of $G$, the multiplicity of the trivial character $1_T$ of $T$ in $s_\mu$. The main case is where $\mu=(q-1)\omega$ and $\omega$ is a fundamental weight of ${\mathbf G}$. Let $St$ denote the Steinberg character of $G$. Then we determine the unipotent characters occurring as constituents of $s_\mu\cdot St$ defined to be 0 at the $p$-singular elements of $G$. Let $\beta_\mu$ denote the Brauer character of a representation of $SL_{n}(q)$ arising from an irreducible representation of ${\mathbf G}$ with highest weight $\mu$. Then we determine the unipotent constituents of the characters $\beta_\mu\cdot St$ for $\mu=(q-1)\omega$, and also for some other $\mu$ (called strongly $q$-restricted). In addition, for strongly restricted weights $\mu$, we compute the \mult of $1_T$ in the restriction $\beta_\mu|_T$ for every maximal torus $T$ of $G$.


title: The Bi-Lipschitz Equisingularity of Essentially Isolated Determinantal Singularities

abstract: The bi-Lipschitz geometry is one of the main subjects in the modern approach of Singularity Theory. However, it rises from works of important mathematicians of the last century, especially Zariski. In this work we investigate the Bi-Lipschitz equisingularity of families of Essentially Isolated Determinantal Singularities inspired by the approach of Mostowski and Gaffney.


title: Mittag-Leffler functions and the truncated $\mathcal{V}$-fractional derivative

abstract: We introduce a new derivative, the so-called truncated $\mathcal{V}$-fractional derivative for $\alpha$-differentiable functions through the six parameters truncated Mittag-Leffler function, which generalizes different fractional derivatives, recently introduced: conformable fractional derivatives, alternative fractional derivative, truncated alternative fractional derivative, $M$-fractional derivative and truncated $M$-fractional derivative. This new truncated $\mathcal{V}$-fractional derivative satisfies properties of the entire order calculus, among them: linearity, product rule, quotient rule, function composition, and chain rule. Also, as in the case of the Caputo derivative, the derivative of a constant is zero. Since the six parameters Mittag-Leffler function is a generalization of Mittag-Leffler functions of one, two, three, four, and five parameters, we can extend some of the classic results of the entire order calculus, namely: Rolle's theorem, the mean value theorem and its extension. In addition, we present the theorem involving the law of exponents for derivatives and we calculated the truncated $\mathcal{V}$-fractional derivative of the two parameters Mittag-Leffler function. Finally, we present the $\mathcal{V}$-fractional integral from which, as a natural consequence, new results appear as applications. Specifically, we generalize inverse property, the fundamental theorem of calculus, a theorem associated with classical integration by parts, and the mean value theorem for integrals. Also, we calculate the $\mathcal{V}$-fractional integral of the two parameters Mittag-Leffler function. Further, through the truncated $\mathcal{V}$-fractional derivative and the $\mathcal{V}$-fractional integral, we obtain a relation with the fractional derivative and integral in the Riemann-Liouville sense, in the case $0<\alpha<1$.


title: A Novel Simplified Model for Blood Coagulation: A piecewise dynamical model for thrombin with robust predictive capabilities

abstract: Realistic description of patient-specific mechanical properties of clotting dynamics presents a major challenge. Available patient-specific data falls short of robustly characterizing myriads of complex dynamic interactions that happen during clotting. We propose a simplified switching model for a key part of the coagulation cascade that describes dynamics of just four variables. The model correctly predicts prolonged activity of thrombin, an important enzyme in the clotting process, in certain plasma factor compositions. The activity sustains beyond the time which is conventionally considered to be the end of clotting. This observation along with the simplified model is hypothesized as a necessary step towards effectively studying patient-specific properties of clotting dynamics in realistic geometries.