title: Background derivation and image flattening: getimages

abstract: Modern high-resolution images obtained with space observatories display extremely strong intensity variations across images on all spatial scales. Source extraction in such images with methods based on global thresholding may bring unacceptably large numbers of spurious sources in bright areas while failing to detect sources in low-background or low-noise areas. It would be highly beneficial to subtract background and equalize the levels of small-scale fluctuations in the images before extracting sources or filaments. This paper describes getimages, a new method of background derivation and image flattening. It is based on median filtering with sliding windows that correspond to a range of spatial scales from the observational beam size up to a maximum structure width $X_{\lambda}$. The latter is a single free parameter of getimages that can be evaluated manually from the observed image $\mathcal{I}_{\lambda}$. The median filtering algorithm provides a background image $\tilde{\mathcal{B}}_{\lambda}$ for structures of all widths below $X_{\lambda}$. The same median filtering procedure applied to an image of standard deviations $\mathcal{D}_{\lambda}$ derived from a background-subtracted image $\tilde{\mathcal{S}}_{\lambda}$ results in a flattening image $\tilde{\mathcal{F}}_{\lambda}$. Finally, a flattened detection image $\mathcal{I}_{{\lambda}\mathrm{D}}{\,=\,}\tilde{\mathcal{S}}_{\lambda}{/}\tilde{\mathcal{F}}_{\lambda}$ is computed, whose standard deviations are uniform outside sources and filaments. Detecting sources in such greatly simplified images results in much cleaner extractions that are more complete and reliable. As a bonus, getimages reduces various observational and map-making artifacts and equalizes noise levels between independent tiles of mosaicked images.


title: Heavy-quark symmetry implies stable heavy tetraquark mesons $Q_iQ_j \bar q_k \bar q_l$

abstract: For very heavy quarks $Q$, relations derived from heavy-quark symmetry predict the existence of novel narrow doubly heavy tetraquark states of the form $Q_iQ_j \bar q_k \bar q_l$ (subscripts label flavors), where $q$ designates a light quark. By evaluating finite-mass corrections, we predict that double-beauty states composed of $bb\bar u \bar d$, $bb\bar u \bar s$, and $bb\bar d \bar s$ will be stable against strong decays, whereas the double-charm states $cc \bar q_k \bar q_l$, mixed beauty+charm states $bc \bar q_k \bar q_l$, and heavier $bb \bar q_k \bar q_l$ states will dissociate into pairs of heavy-light mesons. Observation of a new double-beauty state through its weak decays would establish the existence of tetraquarks and illuminate the role of heavy color-antitriplet diquarks as hadron constituents.


title: Possible particle-hole instabilities in interacting type-II Weyl semimetals

abstract: Type II Weyl semimetal, a three dimensional gapless topological phase, has drawn enormous interest recently. These topological semimetals enjoy overtilted dispersion and Weyl nodes that separate the particle and hole pocket. Using perturbation renormalization group, we identify possible renormalization of the interaction vertices, which show a tendency toward instability. We further adopt a self-consistent mean-field approach to study possible instability of the type II Weyl semimetals under short-range electron-electron interaction. It is found that the instabilities are much easier to form in type II Weyl semimetals than the type I case. Eight different mean-field orders are identified, among which we further show that the polar charge density wave (CDW) phase exhibits the lowest energy. This CDW order is originated from the nesting of the Fermi surfaces and could be a possible ground state in interacting type II Weyl semimetals.


title: Single shot large field of view imaging with scattering media by spatial demultiplexing

abstract: Optically focusing and imaging through strongly scattering media are challenging tasks but have widespread applications from scientific research to biomedical applications and daily life. Benefiting from the memory effect (ME) for speckle intensity correlations, only one single-shot speckle pattern can be used for the high quality recovery of the objects and avoiding some complicated procedures to reduce scattering effects. In spite of all the spatial information from a large object being embedded in a single speckle image, ME gives a strict limitation to the field of view (FOV) for imaging through scattering media. Objects beyond the ME region cannot be recovered and only produce unwanted speckle patterns, causing reduction in the speckle contrast and recovery quality. Here, we extract the spatial information by utilizing these unavoidable speckle patterns, and enlarge the FOV of the optical imaging system. Regional point spreading functions (PSFs), which are fixed and only need to be recorded once for all time use, are employed to recover corresponding spatial regions of an object by deconvolution algorithm. Then an automatic weighted averaging in an iterative process is performed to obtain the object with significantly enlarged FOV. Our results present an important step toward an advanced imaging technique with strongly scattering media.

categories: physics.optics eess.IV eess.SP physics.app-ph


title: Primordial black holes and second order gravitational waves from ultra-slow-roll inflation

abstract: The next generation of space-borne gravitational wave detectors may detect gravitational waves from extreme mass-ratio inspirals with primordial black holes. To produce primordial black holes which contribute a non-negligible abundance of dark matter and are consistent with the observations, a large enhancement in the primordial curvature power spectrum is needed. For a single field slow-roll inflation, the enhancement requires a very flat potential for the inflaton, and this will increase the number of $e$-folds. To avoid the problem, an ultra-slow-roll inflation at the near inflection point is required. We elaborate the conditions to successfully produce primordial black hole dark matter from single field inflation and propose a toy model with polynomial potential to realize the big enhancement of the curvature power spectrum at small scales while maintaining the consistency with the observations at large scales. The power spectrum for the second order gravitational waves generated by the large density perturbations at small scales is consistent with the current pulsar timing array observations.


title: The SNPR neighbourhood of tree-child networks

abstract: Network rearrangement operations like SNPR (SubNet Prune and Regraft), a recent generalisation of rSPR (rooted Subtree Prune and Regraft), induce a metric on phylogenetic networks. To search the space of these networks one important property of these metrics is the sizes of the neighbourhoods, that is, the number of networks reachable by exactly one operation from a given network. In this paper, we present exact expressions for the SNPR neighbourhood of tree-child networks, which depend on both the size and the topology of a network. We furthermore give upper and lower bounds for the minimum and maximum size of such a neighbourhood.


title: Stability for implicit-explicit schemes for non-equilibrium kinetic systems in weighted spaces with symmetrization

abstract: We consider kinetic systems and prove their stability working in weighted spaces in which the systems are symmetric. We prove stability for various explicit and implicit semi-discrete and fully discrete schemes. The applications include advective and diffusive transport coupled to the accumulation of immobile components governed by non-equilibrium relationships. We also discuss extensions to nonlinear relationships and multiple species.


title: Construction of helices from Lucas and Fibonacci sequences

abstract: By means of two complex-valued functions (depending on an integer parameter P>=1) we construct helices of integer ratio R>=1 related to the so-called Binet formulae for P-Lucas and P-Fibonacci sequences. Based on these functions a new map is defined and we show that its three-dimensional representation is also a helix. After proving that the lattice points of these later helix satisfy certain diophantine Pell's equations we call it a Pell's helix. We prove that for P-Fibonacci and Pell's helices the respective ratio is an invariant, contrasting to the P-Lucas helices whose ratio depends on P. It is also shown that suitable linear combinations of certain complex-valued maps lead to new helices related to Lucas/Fibonacci/Pell numbers. Graphical examples are given in order to illustrate the underlying theory.


title: A Duality Web in Condensed Matter Systems

abstract: We study various dualities in condensed matter systems. The dualities in three dimensions could be derived from a conjecture of a duality between a Dirac fermion theory and an interacting scalar field theory at the Wilson-Fisher fixed point and zero temperature in three dimensions. We show that the dualities are not affected by the non-trivial holonomy, use a mean-field method to study, and discuss the dualities at a finite temperature. Finally, we combine a bulk theory, which is an Abelian $p$-form theory with a theta term in $2p+2$ dimensions, and a boundary theory, which is a $2p+1$ dimensional theory, to discuss constraints and difficulties of a duality web in $2p+1$ dimensions.


title: Blow-up for semilinear damped wave equations with sub-Strauss exponent in the scattering case

abstract: It is well-known that the critical exponent for semilinear damped wave equations is Fujita exponent when the damping is effective. Lai, Takamura and Wakasa in 2017 have obtained a blow-up result not only for super-Fujita exponent but also for the one closely related to Strauss exponent when the damping is scaling invariant and its constant is relatively small,which has been recently extended by Ikeda and Sobajima. Introducing a multiplier for the time-derivative of the spatial integral of unknown functions, we succeed in employing the technics on the analysis for semilinear wave equations and proving a blow-up result for semilinear damped wave equations with sub-Strauss exponent when the damping is in the scattering range.


title: Entropy decay for the Kac evolution

abstract: We consider solutions to the Kac master equation for initial conditions where $N$ particles are in a thermal equilibrium and $M\le N$ particles are out of equilibrium. We show that such solutions have exponential decay in entropy relative to the thermal state. More precisely, the decay is exponential in time with an explicit rate that is essentially independent on the particle number. This is in marked contrast to previous results which show that the entropy production for arbitrary initial conditions is inversely proportional to the particle number. The proof relies on Nelson's hypercontractive estimate and the geometric form of the Brascamp-Lieb inequalities due to Franck Barthe. Similar results hold for the Kac-Boltzmann equation with uniform scattering cross sections.


title: Virtual PET Images from CT Data Using Deep Convolutional Networks: Initial Results

abstract: In this work we present a novel system for PET estimation using CT scans. We explore the use of fully convolutional networks (FCN) and conditional generative adversarial networks (GAN) to export PET data from CT data. Our dataset includes 25 pairs of PET and CT scans where 17 were used for training and 8 for testing. The system was tested for detection of malignant tumors in the liver region. Initial results look promising showing high detection performance with a TPR of 92.3% and FPR of 0.25 per case. Future work entails expansion of the current system to the entire body using a much larger dataset. Such a system can be used for tumor detection and drug treatment evaluation in a CT-only environment instead of the expansive and radioactive PET-CT scan.


title: Lambda number of the power graph of a finite group

abstract: The power graph $\Gamma_G$ of a finite group $G$ is the graph with the vertex set $G$, where two distinct elements are adjacent if one is a power of the other. An $L(2, 1)$-labeling of a graph $\Gamma$ is an assignment of labels from nonnegative integers to all vertices of $\Gamma$ such that vertices at distance two get different labels and adjacent vertices get labels that are at least $2$ apart. The lambda number of $\Gamma$, denoted by $\lambda(\Gamma)$, is the minimum span over all $L(2, 1)$-labelings of $\Gamma$. In this paper, we obtain bounds for $\lambda(\Gamma_G)$, and give necessary and sufficient conditions when the bounds are attained. As applications, we compute the exact value of $\lambda(\Gamma_G)$ if $G$ is a dihedral group, a generalized quaternion group, a $\mathcal{P}$-group or a cyclic group of order $pq^n$, where $p$ and $q$ are distinct primes and $n$ is a positive integer.


title: Network modelling of topological domains using Hi-C data

abstract: Chromosome conformation capture experiments such as Hi-C are used to map the three-dimensional spatial organization of genomes. One specific feature of the 3D organization is known as topologically associating domains (TADs), which are densely interacting, contiguous chromatin regions playing important roles in regulating gene expression. A few algorithms have been proposed to detect TADs. In particular, the structure of Hi-C data naturally inspires application of community detection methods. However, one of the drawbacks of community detection is that most methods take exchangeability of the nodes in the network for granted; whereas the nodes in this case, i.e. the positions on the chromosomes, are not exchangeable. We propose a network model for detecting TADs using Hi-C data that takes into account this non-exchangeability. In addition, our model explicitly makes use of cell-type specific CTCF binding sites as biological covariates and can be used to identify conserved TADs across multiple cell types. The model leads to a likelihood objective that can be efficiently optimized via relaxation. We also prove that when suitably initialized, this model finds the underlying TAD structure with high probability. Using simulated data, we show the advantages of our method and the caveats of popular community detection methods, such as spectral clustering, in this application. Applying our method to real Hi-C data, we demonstrate the domains identified have desirable epigenetic features and compare them across different cell types.


title: Nonperturbative partonic quasidistributions of the pion from chiral quark models

abstract: We evaluate nonperturbatively the quark quasidistribution amplitude and the valence quark quasidistribution function of the pion in the framework of chiral quark models, namely the Nambu--Jona-Lasinio model and the Spectral Quark Model. We arrive at simple analytic expressions, where the nonperturbative dependence on the longitudinal momentum of the pion can be explicitly assessed. The model results for the quark quasidistribution amplitude of the pion compare favorably to the data obtained from the Euclidean lattice simulations. The quark distribution amplitude, arising in the limit of infinite longitudinal momentum of the pion, agrees, after suitable QCD evolution, to the recent data extracted from Euclidean lattices, as well as to the old data from transverse lattice simulations.


title: Dispersive SYK model: band structure and quantum chaos

abstract: The Sachdev-Ye-Kitaev (SYK) model is a concrete model for non-Fermi Liquid with maximally chaotic behavior in $0+1$-$d$. In order to gain some insights into real materials in higher dimensions where fermions could hop between different sites, here we consider coupling a SYK lattice by a constant hopping. We call this dispersive SYK model. Focusing on $1+1$-$d$ homogeneous hopping, by either tuning temperature or the relative strength of random interaction (hopping) and constant hopping, we find a crossover between a dispersive metal to an incoherent metal, where dynamic exponent $z$ changes from $1$ to $\infty$. We study the crossover by calculating spectral function, charge density correlator and the Lyapunov exponent. We further find the Lyapunov exponent becomes larger when the chemical potential is tuned to approach a Van Hove singularity because of the large density of states near the Fermi suface. The effect of the topological non-trivial bands is also discussed.


title: Effect of magnetic field on the burning of a neutron star

abstract: In this article, we present the effect of a strong magnetic field in the burning of a neutron star (NS). We have used relativistic magneto-hydrostatic (MHS) conservation equations for studying the PT from nuclear matter (NM) to quark matter (QM). We found that the shock-induced phase transition (PT) is likely if the density of the star core is more than three times nuclear saturation ($\rho_s$) density. The conversion process from NS to quark star (QS) is found to be an exothermic process beyond such densities. The burning process at the star center most likely starts as a deflagration process. However, there can be a small window at lower densities where the process can be a detonation one. At small enough infalling matter velocities the resultant magnetic field of the QS is lower than that of the NS. However, for a higher value of infalling matter velocities, the magnetic field of QM becomes larger. Therefore, depending on the initial density fluctuation and on whether the PT is a violent one or not the QS could be more magnetic or less magnetic. The PT also have a considerable effect on the tilt of the magnetic axis of the star. For smaller velocities and densities the magnetic angle are not affected much but for higher infalling velocities tilt of the magnetic axis changes suddenly. The magnetic field strength and the change in the tilt axis can have a significant effect on the observational aspect of the magnetars.


title: Effects of Quantum Coherence on Work Statistics

abstract: In the conventional two-point measurement scheme of quantum thermodynamics, quantum coherence is destroyed by the first measurement. But as we know the coherence really plays an important role in the quantum thermodynamics process, and how to describe the work statistics for a quantum coherent process is still an open question. In this paper, we use the full counting statistics method to investigate the effects of quantum coherence on work statistics. First, we give a general discussion and show that for a quantum coherent process, work statistics is very different from that of the two-point measurement scheme, specifically the average work is increased or decreased and the work fluctuation can be decreased by quantum coherence, which strongly depends on the relative phase, the energy level structure and the external protocol. Then, we concretely consider a quenched 1-D transverse Ising model, and show that quantum coherence has a more significant influence on work statistics in the ferromagnetism regime compared with that in the paramagnetism regime, so that due to the presence of quantum coherence the work statistics can exhibit the critical phenonmenon even at high temperature.


title: Binary Hypothesis Testing with Byzantine Sensors: Fundamental Trade-off Between Security and Efficiency

abstract: This paper studies binary hypothesis testing based on measurements from a set of sensors, a subset of which can be compromised by an attacker. The measurements from a compromised sensor can be manipulated arbitrarily by the adversary. The asymptotic exponential rate, with which the probability of error goes to zero, is adopted to indicate the detection performance of a detector. In practice, we expect the attack on sensors to be sporadic, and therefore the system may operate with all the sensors being benign for extended period of time. This motivates us to consider the trade-off between the detection performance of a detector, i.e., the probability of error, when the attacker is absent (defined as efficiency) and the worst-case detection performance when the attacker is present (defined as security). We first provide the fundamental limits of this trade-off, and then propose a detection strategy that achieves these limits. We then consider a special case, where there is no trade-off between security and efficiency. In other words, our detection strategy can achieve the maximal efficiency and the maximal security simultaneously. Two extensions of the secure hypothesis testing problem are also studied and fundamental limits and achievability results are provided: 1) a subset of sensors, namely "secure" sensors, are assumed to be equipped with better security countermeasures and hence are guaranteed to be benign, 2) detection performance with unknown number of compromised sensors. Numerical examples are given to illustrate the main results.


title: Discover and Learn New Objects from Documentaries

abstract: Despite the remarkable progress in recent years, detecting objects in a new context remains a challenging task. Detectors learned from a public dataset can only work with a fixed list of categories, while training from scratch usually requires a large amount of training data with detailed annotations. This work aims to explore a novel approach -- learning object detectors from documentary films in a weakly supervised manner. This is inspired by the observation that documentaries often provide dedicated exposition of certain object categories, where visual presentations are aligned with subtitles. We believe that object detectors can be learned from such a rich source of information. Towards this goal, we develop a joint probabilistic framework, where individual pieces of information, including video frames and subtitles, are brought together via both visual and linguistic links. On top of this formulation, we further derive a weakly supervised learning algorithm, where object model learning and training set mining are unified in an optimization procedure. Experimental results on a real world dataset demonstrate that this is an effective approach to learning new object detectors.


title: Ozone Measurements with Meteors: A Revisit

abstract: Understanding the role of ozone in the Mesosphere/Lower Thermosphere (MLT) region is essential for understanding the atmospheric processes in the upper atmosphere. Earlier studies have shown that it is possible to use overdense meteor trails to measure ozone concentration in the meteor region. Here we revisit this topic by comparing a compilation of radar observations to satellite measurements. We observe a modest agreement between the values derived from these two methods, which confirm the usefulness of the meteor trail technique for measuring ozone content at certain heights in the MLT region. Future simultaneous measurements will help quantifying the performance of this technique.


title: Heat Reduction by Thermal Wave Crystals

abstract: Non-Fourier heat conduction models assume wave-like behavior does exist in the heat conduction process. Based on this wave-like behavior, thermal conduction controlled in a one-dimensional periodical structure, named thermal wave crystal, has been demonstrated through both theoretical analysis and numerical simulation based on the Cattaneo-Vernotte (CV) heat-conduction model. The transfer matrix method and Bloch analysis have been applied to calculate the band structure of thermal wave propagating in thermal wave crystals. And the temperature responses are obtained by using the FDTD method, which is also used to verify the correctness of the band structure. The results show that band gaps do exist due to the Bragg scattering. Then, a calculation method to predict the mid-gap frequency of band gaps for the thermal wave crystal has been introduced in this Letter. And key parameters determining the band gaps have been discussed. This study shows the potential applications of this novel mechanism, such as thermal imagining, thermal diodes and thermal waveguides.


title: Pointwise estimates of solutions to nonlinear equations for nonlocal operators

abstract: We study pointwise behavior of positive solutions to nonlinear integral equations, and related inequalities, of the type \begin{equation*} u(x) - \int_\Omega G(x, y) \, g(u(y)) d \sigma (y) = h, \end{equation*} where $(\Omega, \sigma)$ is a locally compact measure space, $G(x, y)\colon \Omega\times \Omega \to [0, +\infty]$ is a kernel, $h \ge 0$ is a measurable function, and $g\colon [0, \infty)\to [0, \infty)$ is a monotone function. This problem is motivated by the semilinear fractional Laplace equation \begin{equation*} (-\Delta)^{\frac{\alpha}{2}} u - g(u) \sigma = \mu \quad \text{in} \, \, \Omega, \quad u=0 \, \, \, \text{in} \, \, \Omega^c, \end{equation*} with measure coefficients $\sigma$, $\mu$, where $g(u)=u^q$, $q \in \mathbb{R} \setminus\{0\}$, and $0<\alpha<n$, in domains $\Omega \subseteq\mathbb{R}^n$, or Riemannian manifolds, with positive Green's function $G$.


title: ScanNet: A Fast and Dense Scanning Framework for Metastatic Breast Cancer Detection from Whole-Slide Images

abstract: Lymph node metastasis is one of the most significant diagnostic indicators in breast cancer, which is traditionally observed under the microscope by pathologists. In recent years, computerized histology diagnosis has become one of the most rapidly expanding fields in medical image computing, which alleviates pathologists' workload and reduces misdiagnosis rate. However, automatic detection of lymph node metastases from whole slide images remains a challenging problem, due to the large-scale data with enormous resolutions and existence of hard mimics. In this paper, we propose a novel framework by leveraging fully convolutional networks for efficient inference to meet the speed requirement for clinical practice, while reconstructing dense predictions under different offsets for ensuring accurate detection on both micro- and macro-metastases. Incorporating with the strategies of asynchronous sample prefetching and hard negative mining, the network can be effectively trained. Extensive experiments on the benchmark dataset of 2016 Camelyon Grand Challenge corroborated the efficacy of our method. Compared with the state-of-the-art methods, our method achieved superior performance with a faster speed on the tumor localization task and surpassed human performance on the WSI classification task.


title: A sparse-grid isogeometric solver

abstract: Isogeometric Analysis (IGA) typically adopts tensor-product splines and NURBS as a basis for the approximation of the solution of PDEs. In this work, we investigate to which extent IGA solvers can benefit from the so-called sparse-grids construction in its combination technique form, which was first introduced in the early 90s in the context of the approximation of high-dimensional PDEs. The tests that we report show that, in accordance to the literature, a sparse-grid construction can indeed be useful if the solution of the PDE at hand is sufficiently smooth. Sparse grids can also be useful in the case of non-smooth solutions when some a-priori knowledge on the location of the singularities of the solution can be exploited to devise suitable non-equispaced meshes. Finally, we remark that sparse grids can be seen as a simple way to parallelize pre-existing serial IGA solvers in a straightforward fashion, which can be beneficial in many practical situations.


title: Dimensional Affect and Expression in Natural and Mediated Interaction

abstract: There is a perceived controversy as to whether the cognitive representation of affect is better modelled using a dimensional or categorical theory. This paper first suggests that these views are, in fact, compatible. The paper then discusses this theme and related issues in reference to a commonly stated application domain of research on human affect and expression: human computer interaction (HCI). The novel suggestion here is that a more realistic framing of studies of human affect in expression with reference to HCI and, particularly HCHI (Human-Computer-Human Interaction) entails some re-formulation of the approach to the basic phenomena themselves. This theme is illustrated with several examples from several recent research projects.


title: Curved momentum spaces from quantum groups with cosmological constant

abstract: We bring the concept that quantum symmetries describe theories with nontrivial momentum space properties one step further, looking at quantum symmetries of spacetime in presence of a nonvanishing cosmological constant $\Lambda$. In particular, the momentum space associated to the $\kappa$-deformation of the de Sitter algebra in (1+1) and (2+1) dimensions is explicitly constructed as a dual Poisson-Lie group manifold parametrized by $\Lambda$. Such momentum space includes both the momenta associated to spacetime translations and the `hyperbolic' momenta associated to boost transformations, and has the geometry of (half of) a de Sitter manifold. Known results for the momentum space of the $\kappa$-Poincar\'e algebra are smoothly recovered in the limit $\Lambda\to 0$, where hyperbolic momenta decouple from translational momenta. The approach here presented is general and can be applied to other quantum deformations of kinematical symmetries, including (3+1)-dimensional ones.


title: Chain Recurrence For General Spaces

abstract: The chain relation, due to Conley, and the strong chain relation, due to Easton, are well studied for continuous maps on compact metric spaces. Following Fathi and Pageault, we use barrier functions to generalize the theory to general relations on uniform spaces. In developing the theory, we indicate why the chain ideas are naturally uniform spaces concepts. We illustrate that the extension to relations is easy and is useful even for the study of the continuous map case.


title: Robust stability conditions for feedback interconnections of distributed-parameter negative imaginary systems

abstract: Sufficient and necessary conditions for the stability of positive feedback interconnections of negative imaginary systems are derived via an integral quadratic constraint (IQC) approach. The IQC framework accommodates distributed-parameter systems with irrational transfer function representations, while generalising existing results in the literature and allowing exploitation of flexibility at zero and infinite frequencies to reduce conservatism in the analysis. The main results manifest the important property that the negative imaginariness of systems gives rise to a certain form of IQCs on positive frequencies that are bounded away from zero and infinity. Two additional sets of IQCs on the DC and instantaneous gains of the systems are shown to be sufficient and necessary for closed-loop stability along a homotopy of systems.


title: Occlusion Handling using Semantic Segmentation and Visibility-Based Rendering for Mixed Reality

abstract: Real-time occlusion handling is a major problem in outdoor mixed reality system because it requires great computational cost mainly due to the complexity of the scene. Using only segmentation, it is difficult to accurately render a virtual object occluded by complex objects such as trees, bushes etc. In this paper, we propose a novel occlusion handling method for real-time, outdoor, and omni-directional mixed reality system using only the information from a monocular image sequence. We first present a semantic segmentation scheme for predicting the amount of visibility for different type of objects in the scene. We also simultaneously calculate a foreground probability map using depth estimation derived from optical flow. Finally, we combine the segmentation result and the probability map to render the computer generated object and the real scene using a visibility-based rendering method. Our results show great improvement in handling occlusions compared to existing blending based methods.


title: Elicitability and its Application in Risk Management

abstract: Elicitability is a property of $\mathbb{R}^k$-valued functionals defined on a set of distribution functions. These functionals represent statistical properties of a distribution, for instance its mean, variance, or median. They are called elicitable if there exists a scoring function such that the expected score under a distribution takes its unique minimum at the functional value of this distribution. If such a scoring function exists, it is called strictly consistent for the functional. Motivated by the recent findings of Fissler and Ziegel concerning higher order elicitability, this thesis reviews the most important results, examples, and applications which are found in the relevant literature. Moreover, we also contribute our own examples and findings in order to give the reader a well-founded overview of the topic as well as of the most used tools and techniques. We include necessary and sufficient conditions for strictly consistent scoring functions, several elicitable as well as non-elicitable functionals and the use of elicitability in forecast comparison, regression, and estimation. Special emphasis is placed on quantitative risk management and the result that Value at Risk and Expected Shortfall are jointly elicitable.


title: CNN-based Cascaded Multi-task Learning of High-level Prior and Density Estimation for Crowd Counting

abstract: Estimating crowd count in densely crowded scenes is an extremely challenging task due to non-uniform scale variations. In this paper, we propose a novel end-to-end cascaded network of CNNs to jointly learn crowd count classification and density map estimation. Classifying crowd count into various groups is tantamount to coarsely estimating the total count in the image thereby incorporating a high-level prior into the density estimation network. This enables the layers in the network to learn globally relevant discriminative features which aid in estimating highly refined density maps with lower count error. The joint training is performed in an end-to-end fashion. Extensive experiments on highly challenging publicly available datasets show that the proposed method achieves lower count error and better quality density maps as compared to the recent state-of-the-art methods.


title: Application of Superhalogens in the Design of Organic Superconductors

abstract: Bechgaard salts, (TMTSF)2X (TMTSF = tetramethyl tetraselenafulvalene and X = complex anion), form the most popular series of organic superconductors. In these salts, TMTSF molecules act as super-electron donor and X as acceptor. We computationally examine the electronic structure and properties of X in commonly used (TMTSF)2X (X = NO3, BF4, ClO4, PF6) superconductors and notice that they belong to the class of superhalogens due to their higher vertical detachment energy than halogen anions. This prompted us to choose other superhalogens such as X = BO2, BH4, B2F7, AuF6 and study their (TMTSF)2X complexes. Our findings suggest that these complexes behave more or less similar to those of established (TMTSF)2X superconductors, particularly for X = BO2 and B2F7. We, therefore, believe that the concept of superhalogen can be successfully applied in the design of novel organic superconductors.


title: Multi-second magnetic coherence in a single domain spinor Bose-Einstein condensate

abstract: We describe a compact, robust and versatile system for studying magnetic dynamics in a spinor Bose-Einstein condensate. Condensates of 87 Rb are produced by all-optical evaporation in a 1560 nm optical dipole trap, using a non-standard loading sequence that employs an auxiliary 1529 nm beam for partial compensation of the strong differential light shift induced by the dipole trap itself. We use near-resonance Faraday rotation probing to non-destructively track the condensate magnetization, and demonstrate few-Larmor-cycle tracking with no detectable degradation of the spin polarization. In the ferromagnetic F = 1 ground state, we observe magnetic T1 and T2$^*$ coherence times limited only by the several-second residence time of the atoms in the trap.


title: Controlled Modulation of Depolarization in Laser Speckle

abstract: A new technique based on superposition of two speckle patterns is proposed and demonstrated for controlled modulation of the spatial polarization distribution of the resultant speckle. It is demonstrated both theoretically and experimentally that controlled modulation of the spatial polarization distribution of laser speckle can be achieved by proper choice of the polarization states as well as the average spatial intensity of the constituent speckles. It is also shown that the proposed technique is useful to generate different speckle patterns with sinusoidal variation in their degree of polarization, which can be tuned from zero to unity. This technique can find applications in sensing, biomedical studies, and in determining the rotation of the electric field vector after passing through a scattering medium.


title: Explicit expressions for European option pricing under a generalized skew normal distribution

abstract: Under a generalized skew normal distribution we consider the problem of European option pricing. Existence of the martingale measure is proved. An explicit expression for a given European option price is presented in terms of the cumulative distribution function of the univariate skew normal and the bivariate standard normal distributions. Some special cases are investigated in a greater detail. To carry out the sensitivity of the option price to the skew parameters, numerical methods are applied. Some concluding remarks and further works are given. The results obtained are extensions of the results provided by [4].


title: Soft communities in similarity space

abstract: The $\mathbb{S}^1$ model has been a central geometric model in the development of the field of network geometry. It has been mainly studied in its homogeneous regime, in which angular coordinates are independently and uniformly scattered on the circle. We now investigate if the model can generate networks with targeted topological features and soft communities, that is, heterogeneous angular distributions. Under these circumstances, hidden degrees must depend on angular coordinates and we propose a method to estimate them. We conclude that the model can be topologically invariant with respect to the soft-community structure. Our results might have important implications, both in expanding the scope of the model beyond the independent hidden variables limit and in the embedding of real-world networks.


title: Joint Named Entity Recognition and Stance Detection in Tweets

abstract: Named entity recognition (NER) is a well-established task of information extraction which has been studied for decades. More recently, studies reporting NER experiments on social media texts have emerged. On the other hand, stance detection is a considerably new research topic usually considered within the scope of sentiment analysis. Stance detection studies are mostly applied to texts of online debates where the stance of the text owner for a particular target, either explicitly or implicitly mentioned in text, is explored. In this study, we investigate the possible contribution of named entities to the stance detection task in tweets. We report the evaluation results of NER experiments as well as that of the subsequent stance detection experiments using named entities, on a publicly-available stance-annotated data set of tweets. Our results indicate that named entities obtained with a high-performance NER system can contribute to stance detection performance on tweets.


title: Magnetic and quadrupole moments of the $Z_c(3900)$

abstract: The electromagnetic properties of the tetraquark state $Z_c(3900)$ are investigated in the diquark-antidiquark picture and its magnetic and quadrupole moments are extracted. To this end, the light-cone QCD sum rule in electromagnetic background field is used. The magnetic and quadrupole moments encode the spatial distributions of the charge and magnetization in the particle. The result obtained for the magnetic moment is quite large and can be measured in future experiments. We obtain a nonzero but small value for the quadrupole moment of $Zc(3900)$ indicating a nonspherical charge distribution.


title: Sparse Vector Recovery: Bernoulli-Gaussian Message Passing

abstract: Low-cost message passing (MP) algorithm has been recognized as a promising technique for sparse vector recovery. However, the existing MP algorithms either focus on mean square error (MSE) of the value recovery while ignoring the sparsity requirement, or support error rate (SER) of the sparse support (non-zero position) recovery while ignoring its value. A novel low-complexity Bernoulli-Gaussian MP (BGMP) is proposed to perform the value recovery as well as the support recovery. Particularly, in the proposed BGMP, support-related Bernoulli messages and value-related Gaussian messages are jointly processed and assist each other. In addition, a strict lower bound is developed for the MSE of BGMP via the genie-aided minimum mean-square-error (GA-MMSE) method. The GA-MMSE lower bound is shown to be tight in high signal-to-noise ratio. Numerical results are provided to verify the advantage of BGMP in terms of final MSE, SER and convergence speed.


title: Life efficiency does not always increase with the dissipation rate

abstract: There does not exist a general positive correlation between important life-supporting properties and the entropy production rate. The simple reason is that nondissipative and time-symmetric kinetic aspects are also relevant for establishing optimal functioning. In fact those aspects are even crucial in the nonlinear regimes around equilibrium where we find biological processing on mesoscopic scales. We make these claims specific via examples of molecular motors, of circadian cycles and of sensory adaptation, whose performance in some regimes is indeed spoiled by increasing the dissipated power. We use the relation between dissipation and the amount of time-reversal breaking to keep the discussion quantitative also in effective models where the physical entropy production is not clearly identifiable.

categories: cond-mat.stat-mech nlin.AO physics.bio-ph q-bio.SC


title: The Bose-Einstein Correlations and the strong coupling constant at low energies

abstract: It is shown that $\alpha_s(E)$, the strong coupling constant, can be determined in the non-perturbative regime from Bose-Einstein correlations (BEC). The obtained $\alpha_s(E)$ is in agreement with the prescriptions dealt with in the Analytic Perturbative Theory approach. It also extrapolates smoothly to the standard perturbative $\alpha_s(E)$ at higher energies. Our results indicate that BEC dimension can be considered as an alternative approach to the short range measure between hadrons.


title: Owl: A General-Purpose Numerical Library in OCaml

abstract: Owl is a new numerical library developed in the OCaml language. It focuses on providing a comprehensive set of high-level numerical functions so that developers can quickly build up data analytical applications. In this abstract, we will present Owl's design, core components, and its key functionality.


title: Maximum coherence in the optimal basis

abstract: The resource theoretic measure of quantum coherence is basis dependent, and the amount of coherence contained in a state is different in different bases. We obtained analytical solutions for the maximum coherence by optimizing the reference basis and highlighted the essential role of the mutually unbiased bases (MUBs) on attaining the maximum coherence. Apart from the relative entropy of coherence, we showed that the MUBs are optimal for the robustness of coherence, the coherence weight, and the modified skew information measure of coherence for any state. Moreover, the MUBs are optimal for all the faithful coherence measures if the state is pure or is of the single qubit. We also highlighted an upper bound for the $l_1$ norm of coherence and compared it with the other bounds as well as the maximum one attainable by optimizing the reference basis.


title: Critical values of homology classes of loops and positive curvature

abstract: We study compact and simply-connected Riemannian manifolds with positive sectional curvature $K\ge 1.$ For a non-trivial homology class of lowest dimension in the space of loops based at a point $p$ or in the free loop space one can define a critical length ${\sf crl}_p\left(M,g\right)$ resp. ${\sf crl}\left(M,g\right).$ Then ${\sf crl}_p\left(M,g\right)$ equals the length of a geodesic loop and ${\sf crl}\left(M,g\right)$ equals the length of a closed geodesic. This is the idea of the proof of the existence of a closed geodesic of positive length presented by Birkhoff in case of a sphere and by Lusternik and Fet in the general case. It is the main result of the paper that the numbers ${\sf crl}_p\left(M,g\right)$ resp. ${\sf crl}\left(M,g\right)$ attain its maximal value $2\pi$ only for the round metric on the $n$-sphere. Under the additional assumption $K \le 4$ this result for ${\sf crl}\left(M,g\right)$ follows from results by Sugimoto in even dimensions and Ballmann, Thorbergsson and Ziller in odd dimensions.


title: Effect of off-diagonal exciton-phonon coupling on intramolecular singlet fission

abstract: Intramolecular singlet fission (iSF) materials provide remarkable advantages in terms of tunable electronic structures, and quantum chemistry studies have indicated strong electronic coupling modulation by high frequency phonon modes. In this work, we formulate a microscopic model of iSF with simultaneous diagonal and off-diagonal coupling to high-frequency modes. A non-perturbative treatment, the Dirac-Frenkel time-dependent variational approach is adopted using the multiple Davydov trial states. It is shown that both diagonal and off-diagonal coupling can aid efficient singlet fission if excitonic coupling is weak, and fission is only facilitated by diagonal coupling if excitonic coupling is strong. In the presence of off-diagonal coupling, it is found that high frequency modes create additional fission channels for rapid iSF. Results presented here may help provide guiding principles for design of efficient singlet fission materials by directly tuning singlet-triplet interstate coupling.


title: Flat optical conductivity in ZrSiS due to two-dimensional Dirac bands

abstract: ZrSiS exhibits a frequency-independent interband conductivity $\sigma(\omega) = \rm{const}(\omega) \equiv \sigma_{\rm{flat}}$ in a broad range from 250 to 2500 cm$^{-1}$ (30 - 300 meV). This makes ZrSiS similar to (quasi)two-dimensional Dirac electron systems, such as graphite and graphene. We assign the flat optical conductivity to the transitions between quasi-two-dimensional Dirac bands near the Fermi level. In contrast to graphene, $\sigma_{\rm{flat}}$ is not supposed to be universal but related to the length of the nodal line in the reciprocal space, $k_{0}$. When $\sigma_{\rm{flat}}$ and $k_{0}$ are connected by a simple model, we find good agreement between experiment and theory. Due to the spin-orbit coupling, the discussed Dirac bands in ZrSiS possess a small gap $\Delta$, for which we determine an upper bound max($\Delta$) = 30 meV from our optical measurements. At low temperatures the momentum-relaxation rate collapses, and the characteristic length scale of momentum relaxation is of the order of microns below 50 K.


title: AGN feedback and its quenching efficiency

abstract: In the last decade, observations have accumulated on gas outflows in galaxies, and in particular massive molecular ones. The mass outflow rate is estimated between 1-5 times the star formation rate. For the highest maximal velocities, they are driven by AGN; these outflows are therefore a clear way to moderate or suppress star formation. Some of the most convincing examples at low redshift come from the radio mode, when the radio jets are inclined towards the galaxy plane, or expand in the hot intra-cluster medium, in cool core clusters. However, AGN feedback can also be positive in many occasions, and the net effect is difficult to evaluate. The quenching efficiency is discussed in view of recent observations.


title: Convergence rates of solutions for a two-species chemotaxis-Navier-Stokes sytstem with competitive kinetics

abstract: In this paper, we study the rates of convergence of supposedly given global bounded classical solutions to a two-species chemotaxis-Navier-Stokes system with Lotka-Volterra competitive kinetics. Except in one case where the rate of convergence for the fluid component is expressed in terms of the Poincare constant and the model parameters, all other rates of convergence are shown to be expressible only in terms of the model parameters and the underlying space dimension.


title: Wavelets and Lattice Field Theory

abstract: When continuous fields are expanded in a wavelet basis, a D-dimensional continuum action becomes a (D+1)-dimensional lattice action on the naively discretized Poincare-patch coordinates of an Euclidean AdS(D+1). New possible criteria for acceptable actions open up.


title: Rome teleportation experiment analysed in the Wigner representation: the role of the zeropoint fluctuations in complete one-photon polarization-momentum Bell-state analysis

abstract: The Wigner representation of parametric down conversion in the Heisenberg picture is applied to the study of the Rome teleportation experiment. We investigate the physical meaning of the zeropoint inputs at the different areas of the experimental setup. In particular, we establish a quantitative relationship between the zeropoint sets of modes that are needed for the preparation of the quantum state to be teleported, the idle channels inside the one-photon polarization-momentum Bell-state analyser, and the possibility of performing teleportation of a polarization state whith certainty.


title: Large $N$ expansion of Wilson loops in the Gross-Witten-Wadia matrix model

abstract: We study the large $N$ expansion of winding Wilson loops in the off-critical regime of the Gross-Witten-Wadia (GWW) unitary matrix model. These have been recently considered in arXiv:1705.06542 and computed by numerical methods. We present various analytical algorithms for the precise computation of both the perturbative and instanton corrections to the Wilson loops. In the gapped phase of the GWW model we present the genus five expansion of the one-cut resolvent that captures all winding loops. Then, as a complementary tool, we apply the Periwal-Shevitz orthogonal polynomial recursion to the GWW model coupled to suitable sources and show how it generates all higher genus corrections to any specific loop with given winding. The method is extended to the treatment of instanton effects including higher order $1/N$ corrections. Several explicit examples are fully worked out and a general formula for the next-to-leading correction at general winding is provided. For the simplest cases, our calculation checks exact results from the Schwinger-Dyson equations, but the presented tools have a wider range of applicability.


title: Sneutrino DM in the NMSSM with inverse seesaw mechanism

abstract: In supersymmetric theories like the Next-to-Minimal Supersymmetric Standard Model (NMSSM), the lightest neutralino with bino or singlino as its dominant component is customarily taken as dark matter (DM) candidate. Since light Higgsinos favored by naturalness can strength the couplings of the DM and thus enhance the DM-nucleon scattering rate, the tension between naturalness and DM direct detection results becomes more and more acute with the improved experimental sensitivity. In this work, we extend the NMSSM by inverse seesaw mechanism to generate neutrino mass, and show that in certain parameter space the lightest sneutrino may act as a viable DM candidate, i.e. it can annihilate by multi-channels to get correct relic density and meanwhile satisfy all experimental constraints. The most striking feature of the extension is that the DM-nucleon scattering rate can be naturally below its current experimental bounds regardless of the higgsino mass, and hence it alleviates the tension between naturalness and DM experiments. Other interesting features include that the Higgs phenomenology becomes much richer than that of the original NMSSM due to the relaxed constraints from DM physics and also due to the presence of extra neutrinos, and that the signatures of sparticles at colliders are quite different from those with neutralino as DM candidate.


title: Learning to Infer Graphics Programs from Hand-Drawn Images

abstract: We introduce a model that learns to convert simple hand drawings into graphics programs written in a subset of \LaTeX. The model combines techniques from deep learning and program synthesis. We learn a convolutional neural network that proposes plausible drawing primitives that explain an image. These drawing primitives are like a trace of the set of primitive commands issued by a graphics program. We learn a model that uses program synthesis techniques to recover a graphics program from that trace. These programs have constructs like variable bindings, iterative loops, or simple kinds of conditionals. With a graphics program in hand, we can correct errors made by the deep network, measure similarity between drawings by use of similar high-level geometric structures, and extrapolate drawings. Taken together these results are a step towards agents that induce useful, human-readable programs from perceptual input.


title: A shape theorem for the scaling limit of the IPDSAW at criticality

abstract: In this paper we give a complete characterization of the scaling limit of the critical Interacting Partially Directed Self-Avoiding Walk (IPDSAW) introduced in Zwanzig and Lauritzen (1968). As the system size $L$ diverges, we prove that the set of occupied sites, rescaled horizontally by $L^{2/3}$ and vertically by $L^{1/3}$ converges in law for the Hausdorff distance towards a non trivial random set. This limiting set is built with a Brownian motion $B$ conditioned to come back at the origin at $a_1$ the time at which its geometric area reaches $1$. The modulus of $B$ up to $a_1$ gives the height of the limiting set, while its center of mass process is an independent Brownian motion. Obtaining the shape theorem requires to derive a functional central limit theorem for the excursion of a random walk with Laplace symmetric increments conditioned on sweeping a prescribed geometric area. This result is proven in a companion paper arXiv:1709.06448.


title: Kernel Projection of Latent Structures Regression for Facial Animation Retargeting

abstract: Inspired by kernel methods that have been used extensively in achieving efficient facial animation retargeting, this paper presents a solution to retargeting facial animation in virtual character's face model based on the kernel projection of latent structure (KPLS) regression between semantically similar facial expressions. Specifically, a given number of corresponding semantically similar facial expressions are projected into the latent space. By using the Nonlinear Iterative Partial Least Square method, decomposition of the latent variables is achieved. Finally, the KPLS is achieved by solving a kernalized version of the eigenvalue problem. By evaluating our methodology with other kernel-based solutions, the efficiency of the presented methodology in transferring facial animation to face models with different morphological variations is demonstrated.


title: Constraining Spacetime Nonmetricity with Neutron Spin Rotation in Liquid He-4

abstract: General spacetime nonmetricity coupled to neutrons is studied. In this context, it is shown that certain nonmetricity components can generate a rotation of the neutron's spin. Available data on this effect obtained from slow-neutron propagation in liquid helium are used to constrain isotropic nonmetricity components at the level of $10^{-22}\,$GeV. These results represent the first limit on the nonmetricity $\zeta^{(6)}_2S_{000}$ parameter as well as the first measurement of nonmetricity inside matter.


title: Consistency of survival tree and forest models: splitting bias and correction

abstract: Random survival forest and survival trees are popular models in statistics and machine learning. However, there is a lack of general understanding regarding consistency, splitting rules and influence of the censoring mechanism. In this paper, we investigate the statistical properties of existing methods from several interesting perspectives. First, we show that traditional splitting rules with censored outcomes rely on a biased estimation of the within-node failure distribution. To exactly quantify this bias, we develop a concentration bound of the within-node estimation based on non i.i.d. samples and apply it to the entire forest. Second, we analyze the entanglement between the failure and censoring distributions caused by univariate splits, and show that without correcting the bias at an internal node, survival tree and forest models can still enjoy consistency under suitable conditions. In particular, we demonstrate this property under two cases: a finite-dimensional case where the splitting variables and cutting points are chosen randomly, and a high-dimensional case where the covariates are weakly correlated. Our results can also degenerate into an independent covariate setting, which is commonly used in the random forest literature for high-dimensional sparse models. However, it may not be avoidable that the convergence rate depends on the total number of variables in the failure and censoring distributions. Third, we propose a new splitting rule that compares bias-corrected cumulative hazard functions at each internal node. We show that the rate of consistency of this new model depends only on the number of failure variables, which improves from non-bias-corrected versions. We perform simulation studies to confirm that this can substantially benefit the prediction error.


title: Tree based weighted learning for estimating individualized treatment rules with censored data

abstract: Estimating individualized treatment rules is a central task for personalized medicine. [zhao2012estimating] and [zhang2012robust] proposed outcome weighted learning to estimate individualized treatment rules directly through maximizing the expected outcome without modeling the response directly. In this paper, we extend the outcome weighted learning to right censored survival data without requiring either an inverse probability of censoring weighting or a semiparametric modeling of the censoring and failure times as done in [zhao2015doubly]. To accomplish this, we take advantage of the tree based approach proposed in [zhu2012recursively] to nonparametrically impute the survival time in two different ways. The first approach replaces the reward of each individual by the expected survival time, while in the second approach only the censored observations are imputed by their conditional expected failure times. We establish consistency and convergence rates for both estimators. In simulation studies, our estimators demonstrate improved performance compared to existing methods. We also illustrate the proposed method on a phase III clinical trial of non-small cell lung cancer.


title: Supersymmetric AdS$_2\times \Sigma_2$ solutions from tri-sasakian truncation

abstract: A class of $AdS_2\times \Sigma_2$, with $\Sigma_2$ being a two-sphere or a hyperbolic space, solutions within four-dimensional $N=4$ gauged supergravity coupled to three-vector multiplets with dyonic gauging is identified. The gauged supergravity has non-semisimple $SO(3)\ltimes (\mathbf{T}^3,\hat{\mathbf{T}}^3)$ gauge group and can be obtained from a consistent truncation of eleven-dimensional supergravity on a tri-sasakian manifold. The maximally symmetric vacua contain $AdS_4$ geometries with $N=1,3$ supersymmetry corresponding to $N=1$ and $N=3$ superconformal field theories (SCFTs) in three dimensions. We find supersymmetric solutions of the form $AdS_2\times \Sigma_2$ preserving two supercharges. These solutions describe twisted compactifications of the dual $N=1$ and $N=3$ SCFTs and should arise as near horizon geometries of dyonic black holes in asymptotically $AdS_4$ space-time. Most solutions have hyperbolic horizons although some of them exhibit spherical horizons. These provide a new class of $AdS_2\times \Sigma_2$ geometries with known M-theory origin.


title: Relevant sampling of the short-time Fourier transform of time-frequency localized functions

abstract: We study the random sampling of the short-time Fourier transform of functions that are localized in a compact region in the time-frequency plane. We follow the approach introduced by Bass and Gr\"ochenig for band-limited functions, and show that with a high, controllable probability, a sufficiently dense set of local random samples from the region of concentration in the time-frequency plane yields a sampling inequality for the short-time Fourier transform of time-frequency localized functions on the region.


title: Metric minimizing surfaces revisited

abstract: A surface which does not admit a length nonincreasing deformation is called metric minimizing. We show that metric minimizing surfaces in CAT(0) spaces are locally CAT(0) with respect to their intrinsic metric.


title: LEARN: Learned Experts' Assessment-based Reconstruction Network for Sparse-data CT

abstract: Compressive sensing (CS) has proved effective for tomographic reconstruction from sparsely collected data or under-sampled measurements, which are practically important for few-view CT, tomosynthesis, interior tomography, and so on. To perform sparse-data CT, the iterative reconstruction commonly use regularizers in the CS framework. Currently, how to choose the parameters adaptively for regularization is a major open problem. In this paper, inspired by the idea of machine learning especially deep learning, we unfold a state-of-the-art "fields of experts" based iterative reconstruction scheme up to a number of iterations for data-driven training, construct a Learned Experts' Assessment-based Reconstruction Network ("LEARN") for sparse-data CT, and demonstrate the feasibility and merits of our LEARN network. The experimental results with our proposed LEARN network produces a competitive performance with the well-known Mayo Clinic Low-Dose Challenge Dataset relative to several state-of-the-art methods, in terms of artifact reduction, feature preservation, and computational speed. This is consistent to our insight that because all the regularization terms and parameters used in the iterative reconstruction are now learned from the training data, our LEARN network utilizes application-oriented knowledge more effectively and recovers underlying images more favorably than competing algorithms. Also, the number of layers in the LEARN network is only 12, reducing the computational complexity of typical iterative algorithms by orders of magnitude.


title: Rate of convergence for Hilbert space valued processes

abstract: Consider a stationary, linear Hilbert space valued process. We establish Berry-Essen type results with optimal convergence rates under sharp dependence conditions on the underlying coefficient sequence of the linear operators. The case of non-linear Bernoulli-shift sequences is also considered. If the sequence is $m$-dependent, the optimal rate $(n/m)^{1/2}$ is reached. If the sequence is weakly geometrically dependent, the rate $(n/\log n)^{1/2}$ is obtained.


title: Virtual Crystals and Nakajima Monomials

abstract: An explicit description of the virtualization map for the (modified) Nakajima monomial model for crystals is given. We give an explicit description of the Lusztig data for modified Nakajima monomials in type $A_n$.


title: Finding a best approximation pair of points for two polyhedra

abstract: Given two disjoint convex polyhedra, we look for a best approximation pair relative to them, i.e., a pair of points, one in each polyhedron, attaining the minimum distance between the sets. Cheney and Goldstein showed that alternating projections onto the two sets, starting from an arbitrary point, generate a sequence whose two interlaced subsequences converge to a best approximation pair. We propose a process based on projections onto the half-spaces defining the two polyhedra, which are more negotiable than projections on the polyhedra themselves. A central component in the proposed process is the Halpern--Lions--Wittmann--Bauschke algorithm for approaching the projection of a given point onto a convex set.


title: Negation of photon loss provided by negative weak value

abstract: We propose a usage of a weak value for a quantum processing between preselection and postselection. While the weak value of a projector of 1 provides a process with certainty like the probability of 1, the weak value of -1 negates the process completely. Their mutually opposite effect is approved without a conventional `weak' condition. In addition the quantum process is not limited to be unitary; in particular we consider a loss of photons and experimentally demonstrate the negation of the photon loss by using the negative weak value of -1 against the positive weak value of 1.


title: Towards Visual Explanations for Convolutional Neural Networks via Input Resampling

abstract: The predictive power of neural networks often costs model interpretability. Several techniques have been developed for explaining model outputs in terms of input features; however, it is difficult to translate such interpretations into actionable insight. Here, we propose a framework to analyze predictions in terms of the model's internal features by inspecting information flow through the network. Given a trained network and a test image, we select neurons by two metrics, both measured over a set of images created by perturbations to the input image: (1) magnitude of the correlation between the neuron activation and the network output and (2) precision of the neuron activation. We show that the former metric selects neurons that exert large influence over the network output while the latter metric selects neurons that activate on generalizable features. By comparing the sets of neurons selected by these two metrics, our framework suggests a way to investigate the internal attention mechanisms of convolutional neural networks.


title: Adaptive Performance Optimization under Power Constraint in Multi-thread Applications with Diverse Scalability

abstract: In modern data centers, energy usage represents one of the major factors affecting operational costs. Power capping is a technique that limits the power consumption of individual systems, which allows reducing the overall power demand at both cluster and data center levels. However, literature power capping approaches do not fit well the nature of important applications based on first-class multi-thread technology. For these applications performance may not grow linearly as a function of the thread-level parallelism because of the need for thread synchronization while accessing shared resources, such as shared data. In this paper we consider the problem of maximizing the application performance under a power cap by dynamically tuning the thread-level parallelism and the power state of the CPU-cores. Based on experimental observations, we design an adaptive technique that selects in linear time the optimal combination of thread-level parallelism and CPU-core power state for the specific workload profile of the multi-threaded application. We evaluate our proposal by relying on different benchmarks, configured to use different thread synchronization methods, and compare its effectiveness to different state-of-the-art techniques.


title: A Novel Approach for Image Segmentation based on Histograms computed from Hue-data

abstract: Computer Vision is growing day by day in terms of user specific applications. The first step of any such application is segmenting an image. In this paper, we propose a novel and grass-root level image segmentation algorithm for cases in which the background has uniform color distribution. This algorithm can be used for images of flowers, birds, insects and many more where such background conditions occur. By image segmentation, the visualization of a computer increases manifolds and it can even attain near-human accuracy during classification.


title: Model of energy spectrum parameters of ground level enhancement events in solar cycle 23

abstract: Mewaldt et al. 2012 fitted the observations of the ground level enhancement (GLE) events during solar cycle 23 to the double power-law equation to obtain the four energy spectra parameters, the normalization parameter $C$, low-energy power-law slope $\gamma_1$, high-energy power-law slope $\gamma_2$, and break energy $E_0$. There are 16 GLEs from which we select $13$ for study by excluding some events with complicated situation. We analyze the four parameters with conditions of the corresponding solar events. According to solar event conditions we divide the GLEs into two groups, one with strong acceleration by interplanetary (IP) shocks and another one without strong acceleration. By fitting the four parameters with solar event conditions we obtain models of the parameters for the two groups of GLEs separately. Therefore, we establish a model of energy spectrum of solar cycle 23 GLEs which may be used in prediction in the future.


title: Incommensurate magnetism near quantum criticality in CeNiAsO

abstract: Two phase transitions in the tetragonal strongly correlated electron system CeNiAsO were probed by neutron scattering and zero field muon spin rotation. For $T <T_{N1}$ = 8.7(3) K, a second order phase transition yields an incommensurate spin density wave with wave vector $\textbf{k} = (0.44(4), 0, 0)$. For $T < T_{N2}$ = 7.6(3) K, we find co-planar commensurate order with a moment of $0.37(5)~\mu_B$, reduced to $30 \%$ of the saturation moment of the $|\pm\frac{1}{2}\rangle$ Kramers doublet ground state, which we establish by inelastic neutron scattering. Muon spin rotation in $\rm CeNiAs_{1-x}P_xO$ shows the commensurate order only exists for x $\le$ 0.1 so the transition at $x_c$ = 0.4(1) is from an incommensurate longitudinal spin density wave to a paramagnetic Fermi liquid.


title: Correct Composition of Dephased Behavioural Models

abstract: Scenarios of execution are commonly used to specify partial behaviour and interactions between different objects and components in a system. To avoid overall inconsistency in specifications, various automated methods have emerged in the literature to compose (behavioural) models. In recent work, we have shown how the theorem prover Isabelle can be combined with the constraint solver Z3 to efficiently detect inconsistencies in two or more behavioural models and, in their absence, generate the composition. Here, we extend our approach further and show how to generate the correct composition (as a set of valid traces) of dephased models. This work has been inspired by a problem from a medical domain where different care pathways (for chronic conditions) may be applied to the same patient with different starting points.


title: The Born-Oppenheimer approximation in an effective field theory language

abstract: The Born--Oppenheimer approximation is the standard tool for the study of molecular systems. It is founded on the observation that the energy scale of the electron dynamics in a molecule is larger than that of the nuclei. A very similar physical picture can be used to describe QCD states containing heavy quarks as well as light-quarks or gluonic excitations. In this work, we derive the Born--Oppenheimer approximation for QED molecular systems in an effective field theory framework by sequentially integrating out degrees of freedom living at energies above the typical energy scale where the dynamics of the heavy degrees of freedom occurs. In particular, we compute the matching coefficients of the effective field theory for the case of the $H^+_2$ diatomic molecule that are relevant to compute its spectrum up to ${\cal O}(m\alpha^5)$. Ultrasoft photon loops contribute at this order, being ultimately responsible for the molecular Lamb shift. In the effective field theory the scaling of all the operators is homogeneous, which facilitates the determination of all the relevant contributions, an observation that may become useful for high-precision calculations. Using the above case as a guidance, we construct under some conditions an effective field theory for QCD states formed by a color-octet heavy quark-antiquark pair bound with a color-octet light-quark pair or excited gluonic state, highlighting the similarities and differences between the QED and QCD systems. Assuming that the multipole expansion is applicable, we construct the heavy-quark potential up to next-to-leading order in the multipole expansion in terms of nonperturbative matching coefficients to be obtained from lattice QCD.


title: A note on concordance properties of fibers in Seifert homology spheres

abstract: In this note, we collect various properties of Seifert homology spheres from the viewpoint of Dehn surgery along a Seifert fiber. We expect that many of these are known to various experts, but include them in one place which we hope to be useful in the study of concordance and homology cobordism.


title: Unconventional order-disorder phase transition in improper ferroelectric hexagonal manganites

abstract: The improper ferroelectricity in YMnO$_3$ and other related multiferroic hexagonal manganites are known to cause topologically protected ferroelectric domains that give rise to rich and diverse physical phenomena. The local structure and structural coherence across the ferroelectric transition, however, were previously not well understood. Here we reveal the evolution of the local structure with temperature in YMnO$_3$ using neutron total scattering techniques, and interpret them with the help of first-principles calculations. The results show that, at room temperature, the local and average structures are consistent with the established ferroelectric $P6_3cm$ symmetry. On heating, both local and average structural analyses show striking anomalies from $\sim 800$ K up to the Curie temperature consistent with increasing fluctuations of the order parameter angle. These fluctuations result in an unusual local symmetry lowering into a \textit{continuum of structures} on heating. This local symmetry breaking persists into the high-symmetry non-polar phase, constituting an unconventional type of order-disorder transition.


title: IFU spectroscopy of Southern Planetary Nebulae V: Low-ionization structures

abstract: In this 5th paper of the series, we examine the spectroscopy and morphology of four southern Galactic planetary nebulae Hen 2-141, NGC5307, IC 2553, and PB6 using new integral field spectroscopy data. The morphologies and ionization structures of the sample are given as a set of emission-line maps. In addition, the physical conditions, chemical compositions, and kinematical characteristics of these objects are derived. The results show that PB6 and Hen 2-141 are of very high excitation classes and IC 2553 and NGC5307 are mid to high excitation objects. The elemental abundances reveal that PB6 is of Type I, Hen 2-141 and IC 2553 are of Type IIa, and NGC5307 of Type IIb/III. The observations unveil the presence of well- defined low-ionization structures or \knots" in all objects. The diagnostic diagrams reveal that the excitation mechanism of these knots is probably by photo-ionization of dense material by the nebular central stars. The physical analysis of six of these knots show no significant differences with their surrounding nebular gas, except their lower electron densities. In spite of the enhancement of the low-ionization emission lines of these knots, their chemical abundances are nearly comparable to their surrounding nebulae, with the exception of perhaps slightly higher nitrogen abundances in the NGC5307 knots. The integrated spectrum of IC 2553 reveals that nearly all key lines that have led researchers to characterize its central star as a weak-emission line star type are in fact of nebular origin.


title: Low frequency analogue Hawking radiation: The Korteweg-de Vries model

abstract: We derive analytic expressions for the low-frequency properties of the analogue Hawking radiation in a general weak-dispersive medium. A thermal low-frequency part of the spectrum is expected even when dispersive effects become significant. We consider the two most common class of weak-dispersive media and investigate all possible anomalous scattering processes due inhomogeneous background flows. We first argue that under minimal assumptions, the scattering processes in near-critical flows are well described by a linearized Korteweg-de Vries equation. Within our theoretical model greybody factors are neglected, that is, the mode co-moving with the flow decouples from the other ones. We also exhibit a flow example with an exact expression for the effective temperature. We see that this temperature coincides with the Hawking one only when the dispersive length scale is much smaller than the flow gradient scale. We apply the same method in inhomogeneous flows without an analogue horizon. In this case, the spectrum coefficients decrease with decreasing frequencies. Our findings are in agreement with previous numerical works, generalizing their findings to arbitrary flow profiles. Our analytical expressions provide estimates to guide ongoing experimental efforts.


title: Choosing elements from finite fields

abstract: In two important papers from 1960 Graham Higman introduced the notion of PORC functions, and he proved that for any given positive integer $n$ the number of $p$-class two groups of order $p^n$ is a PORC function of $p$. A key result in his proof of this theorem is the following: "The number of ways of choosing a finite number of elements from the finite field of order $q^n$ subject to a finite number of monomial equations and inequalities between them and their conjugates over GF($q$), considered as a function of $q$, is PORC." Higman's proof of this result involves five pages of homological algebra. Here we give a short elementary proof of the result. Our proof is constructive, and gives an algorithms for computing the relevant PORC functions.


title: Ratios of Ordered Points of Point Processes with Regularly Varying Intensity Measures

abstract: We study limiting properties of ratios of ordered points of point processes whose intensity measures have regularly varying tails, giving a systematic treatment which points the way to "large-trimming" properties of extremal processes and a variety of applications. Our point process approach facilitates a connection with the negative binomial process of Gregoire (1984) and consequently to certain generalised versions of the Poisson-Dirichlet distribution.


title: Pseudo-Prikry sequences

abstract: We generalize results of Gitik, Dzamonja-Shelah, and Magidor-Sinapova on the existence of pseudo-Prikry sequences, which are sequences that approximate the behavior of the generic objects introduced by Prikry-type forcings, in outer models of set theory. Such sequences play an important role in the study of singular cardinal combinatorics by placing restrictions on the type of behavior that can consistently be obtained in outer models. In addition, we provide results about the existence of diagonal pseudo-Prikry sequences, which approximate the behavior of the generic objects introduced by diagonal Prikry-type forcings. Our proof techniques are substantially different from those of previous results and rely on an analysis of PCF-theoretic objects in the outer model.


title: The first order partial differential equations resolved with any derivatives

abstract: In this paper we discuss the first order partial differential equations resolved with any derivatives. At first, we transform the first order partial differential equation resolved with respect to a time derivative into a system of linear equations. Secondly, we convert it into a system of the first order linear partial differential equations with constant coefficients and nonlinear algebraic equations. Thirdly, we solve them by the Fourier transform and convert them into the equivalent integral equations. At last, we extend to discuss the first order partial differential equations resolved with respect to time derivatives and the general scenario resolved with any derivatives.


title: Invertibility via distance for non-centered random matrices with continuous distributions

abstract: Let $A$ be an $n\times n$ random matrix with independent rows $R_1(A),\dots,R_n(A)$, and assume that for any $i\leq n$ and any three-dimensional linear subspace $F\subset {\mathbb R}^n$ the orthogonal projection of $R_i(A)$ onto $F$ has distribution density $\rho(x):F\to{\mathbb R}_+$ satisfying $\rho(x)\leq C_1/\max(1,\|x\|_2^{2000})$ ($x\in F$) for some constant $C_1>0$. We show that for any fixed $n\times n$ real matrix $M$ we have $${\mathbb P}\{s_{\min}(A+M)\leq t n^{-1/2}\}\leq C'\, t,\quad\quad t>0,$$ where $C'>0$ is a universal constant. In particular, the above result holds if the rows of $A$ are independent centered log-concave random vectors with identity covariance matrices. Our method is free from any use of covering arguments, and is principally different from a standard approach involving a decomposition of the unit sphere and coverings, as well as an approach of Sankar-Spielman-Teng for non-centered Gaussian matrices.


title: Heat asymptotics for nonminimal Laplace type operators and application to noncommutative tori

abstract: Let $P$ be a Laplace type operator acting on a smooth hermitean vector bundle $V$ of fiber $\mathbb{C}^N$ over a compact Riemannian manifold given locally by $P= - [g^{\mu\nu} u(x)\partial_\mu\partial_\nu + v^\nu(x)\partial_\nu + w(x)]$ where $u,\,v^\nu,\,w$ are $M_N(\mathbb{C})$-valued functions with $u(x)$ positive and invertible. For any $a \in \Gamma(\text{End}(V))$, we consider the asymptotics $\text{Tr} (a e^{-tP}) \underset{t \downarrow 0^+}{\sim} \,\sum_{r=0}^\infty a_r(a, P)\,t^{(r-d)/2}$ where the coefficients $a_r(a, P)$ can be written locally as $a_r(a, P)(x) = \text{tr}[a(x) \mathcal{R}_r(x)]$. The computation of $\mathcal{R}_2$ is performed opening the opportunity to calculate the modular scalar curvature for noncommutative tori.


title: Calkin representations for $L^{p}$

abstract: We identify the weak closures of the ranges of certain Calkin representations for $L^{p}$, $1<p<\infty$. As a consequence, assuming the continuum hypothesis, we show that the commutant of $B(L^{p})$, $1<p<\infty$, in its ultrapower may or may not be trivial depending on the ultrafilter. This extends a result of Farah, Phillips and Stepr\=ans.


title: Flux Reconstruction for Goal-Oriented A Posteriori Error Estimation

abstract: We propose a new heuristic goal-oriented a posteriori error estimator that connects the dual weighted residual method with equilibrated a posteriori error estimation. Our numerical experiments demonstrate the practical reliability of the error estimator, confirming theoretical predictions, as well as optimally convergent adaptivity even over singular domains and coarse meshes. The central algorithm is a localized flux reconstruction, which has been implemented in the finite element library deal.II. For a solid preparation we assess the performance of the equilibrated a posteriori error estimator of the energy norm in numerical experiments. Moreover, we give what seems to be first rigorous discussion in the numerical literature of localized flux reconstruction over quadrilateral meshes with hanging nodes.


title: On Generated Dynamics for Open Quantum Systems: Spectral Analysis of Effective Liouville

abstract: We point out that the quantum dynamical map of an open quantum system can be generated by an effective Liouville operator. The effective Liouville shows the dynamical breaking of time reversibility. This breaking of reversibility is expressed by the effective Liouville's discrete spectrum having negative imaginary parts. This generated dynamics of open quantum systems is capable of memory effects described by a frequency dependence of the spectrum. When memory effects can be neglected or smoothed out, the effective Liouville generates the well known semi-group dynamics of open systems in the Markov approximation. The spectral analysis of the effective Liouville -with or without memory- allows to represent the quantum dynamical map and expectation values of physical quantities by metric expressions using complex eigenvalues and bi-orthonormal eigenmodes. The long time dynamics is dominated by the zero-mode, which forms the stationary state of the open system. The relaxation time scale $\tau$ is set by the inverse of the smallest negative imaginary part of the non-vanishing effective eigenvalues of the effective Liouville. In generic cases of system-environment-coupling this time scale $\tau$ is the relevant scale for both, the relaxation of diagonal elements (populations) and the decay of off-diagonal elements (coherences) of an initial density matrix in the stationary state's diagonal representation. It also sets the time scale after which the entropy tends to its stationary value. The effective Liouville and its spectral content thus give a framework to study the dynamic breaking of time reversibility, memory effects, decoherence, relaxation to a stationary state and the role of entropy in open quantum systems.


title: A Vision For Continuous Automated Game Design

abstract: ANGELINA is an automated game design system which has previously been built as a single software block which designs games from start to finish. In this paper we outline a roadmap for the development of a new version of ANGELINA, designed to iterate on games in different ways to produce a continuous creative process that will improve the quality of its work, but more importantly improve the perception of the software as being an independently creative piece of software. We provide an initial report of the system's structure here as well as results from the first working module of the system.


title: Adaptive Delivery in Caching Networks

abstract: The problem of content delivery in caching networks is investigated for scenarios where multiple users request identical files. Redundant user demands are likely when the file popularity distribution is highly non-uniform or the user demands are positively correlated. An adaptive method is proposed for the delivery of redundant demands in caching networks. Based on the redundancy pattern in the current demand vector, the proposed method decides between the transmission of uncoded messages or the coded messages of [1] for delivery. Moreover, a lower bound on the delivery rate of redundant requests is derived based on a cutset bound argument. The performance of the adaptive method is investigated through numerical examples of the delivery rate of several specific demand vectors as well as the average delivery rate of a caching network with correlated requests. The adaptive method is shown to considerably reduce the gap between the non-adaptive delivery rate and the lower bound. In some specific cases, using the adaptive method, this gap shrinks by almost 50% for the average rate.


title: Squashed Holography with Scalar Condensates

abstract: We evaluate the partition function of the free and interacting O(N) vector model on a two-parameter family of squashed three spheres in the presence of a scalar deformation. We also find everywhere regular solutions of Einstein gravity coupled to a scalar field in AdS and in dS with the same double squashed boundary geometry. Remarkably, the thermodynamic properties of the AdS solutions qualitatively agree with the behavior predicted by the free O(N) model with a real mass deformation. The dS bulk solutions specify the semiclassical `no-boundary' measure over anisotropic deformations of inflationary, asymptotic de Sitter space. Through dS/CFT the partition function of the interacting O(N) model yields a holographic toy model of the no-boundary measure. We find this yields a qualitatively similar probability distribution which is normalizable and globally peaked at the round three sphere, with a low amplitude for strong anisotropies.


title: Low frequency analogue Hawking radiation: The Bogoliubov-de Gennes model

abstract: We analytically study the low-frequency properties of the analogue Hawking effect in Bose-Einstein condensates. We show that in one-dimensional flows displaying an analogue horizon, the Hawking effect is dominant in the low-frequency regime. This happens despite non vanishing greybody factors, that is, the coupling of the Hawking mode and its partner to the mode propagating with the flow. To show this, we obtained analytical expressions for the scattering coefficients, in general flows and taking into account the full Bogoliubov dispersion relation. We discuss the obtained expressions for the greybody factors. In particular, we show that they can be significantly decreased if the flow obeys a conformal coupling condition. We argue that in the presence of a small but non-zero temperature, reducing greybody factors greatly facilitates the observation of entanglement, that is, establishing that the state of the Hawking mode and its partner is non-separable.


title: Droplet localization in the random XXZ model and its manifestations

abstract: We examine many-body localization properties for the eigenstates that lie in the droplet sector of the random-field spin-$\frac 1 2$ XXZ chain. These states satisfy a basic single cluster localization property (SCLP), derived in \cite{EKS}. This leads to many consequences, including dynamical exponential clustering, non-spreading of information under the time evolution, and a zero velocity Lieb Robinson bound. Since SCLP is only applicable to the droplet sector, our definitions and proofs do not rely on knowledge of the spectral and dynamical characteristics of the model outside this regime. Rather, to allow for a possible mobility transition, we adapt the notion of restricting the Hamiltonian to an energy window from the single particle setting to the many body context.


title: Microscopic theory of the nearest-neighbor valence bond sector of the spin-1/2 kagome antiferromagnet

abstract: The spin-1/2 Heisenberg model on the kagome lattice, which is closely realized in layered Mott insulators such as ZnCu$_3$(OH)$_6$Cl$_2$, is one of the oldest and most enigmatic spin-1/2 lattice model. While the numerical evidence has accumulated in favor of a quantum spin liquid, the debate is still open as to whether it is a $Z_2$ spin liquid with very short-range correlations (some kind of Resonating Valence Bond spin liquid), or an algebraic spin-liquid with power-law correlations. To address this issue, we have pushed the program started by Rokhsar and Kivelson in their derivation of the effective quantum dimer model description of Heisenberg models to unprecedented accuracy for the spin-1/2 kagome, by including all the most important virtual singlet contributions on top of the orthogonalization of the nearest-neighbor valence bond singlet basis. Quite remarkably, the resulting picture is a competition between a $Z_2$ spin liquid and a diamond valence bond crystal with a 12-site unit cell, as in the DMRG simulations of Yan, Huse and White. Furthermore, we found that, on cylinders of finite diameter $d$, there is a transition between the $Z_2$ spin liquid at small $d$ and the diamond valence bond crystal at large $d$, the prediction of the present microscopic description for the 2D lattice. These results show that, if the ground state of the spin-1/2 kagome antiferromagnet can be described by nearest-neighbor singlet dimers, it is a diamond valence bond crystal, and, a contrario, that, if the system is a quantum spin liquid, it has to involve long-range singlets, consistent with the algebraic spin liquid scenario.


title: The First Planetary Microlensing Event with Two Microlensed Source Stars

abstract: We present the analysis of microlensing event MOA-2010-BLG-117, and show that the light curve can only be explained by the gravitational lensing of a binary source star system by a star with a Jupiter mass ratio planet. It was necessary to modify standard microlensing modeling methods to find the correct light curve solution for this binary-source, binary-lens event. We are able to measure a strong microlensing parallax signal, which yields the masses of the host star, $M_* = 0.58\pm 0.11 M_\odot$, and planet $m_p = 0.54\pm 0.10 M_{\rm Jup}$ at a projected star-planet separation of $a_\perp = 2.42\pm 0.26\,$AU, corresponding to a semi-major axis of $a = 2.9{+1.6\atop -0.6}\,$AU. Thus, the system resembles a half-scale model of the Sun-Jupiter system with a half-Jupiter mass planet orbiting a half-solar mass star at very roughly half of Jupiter's orbital distance from the Sun. The source stars are slightly evolved, and by requiring them to lie on the same isochrone, we can constrain the source to lie in the near side of the bulge at a distance of $D_S = 6.9 \pm 0.7\,$kpc, which implies a distance to the planetary lens system of $D_L = 3.5\pm 0.4\,$kpc. The ability to model unusual planetary microlensing events, like this one, will be necessary to extract precise statistical information from the planned large exoplanet microlensing surveys, such as the WFIRST microlensing survey.


title: Handling Nested Parallelism and Extreme Load Imbalance in an Orbital Analysis Code

abstract: Nested parallelism exists in scientific codes that are searching multi-dimensional spaces. However, implementations of nested parallelism often have overhead and load balance issues. The Orbital Analysis code we present exhibits a sparse search space, significant load imbalances, and stopping when the first solution is reached. All these aspects of the algorithm exacerbate the problem of using nested parallelism effectively. In this paper, we present an inspector/executor strategy for chunking such computations into parallel wavefronts. The presented shared memory parallelization is no longer nested and exhibits significantly less load imbalance. We evaluate this approach on an Orbital analysis code, and we improve the execution time from the original implementation by an order of magnitude. As part of a Graduate Computer Science course in Parallel Programming models, we show how the approach can be implemented in parallel Perl, Python, Chapel, Pthreads, and OpenMP. Future work includes investigating how to automate and generalize the parallelization approach.


title: Scalable and Effective Deep CCA via Soft Decorrelation

abstract: Recently the widely used multi-view learning model, Canonical Correlation Analysis (CCA) has been generalised to the non-linear setting via deep neural networks. Existing deep CCA models typically first decorrelate the feature dimensions of each view before the different views are maximally correlated in a common latent space. This feature decorrelation is achieved by enforcing an exact decorrelation constraint; these models are thus computationally expensive due to the matrix inversion or SVD operations required for exact decorrelation at each training iteration. Furthermore, the decorrelation step is often separated from the gradient descent based optimisation, resulting in sub-optimal solutions. We propose a novel deep CCA model Soft CCA to overcome these problems. Specifically, exact decorrelation is replaced by soft decorrelation via a mini-batch based Stochastic Decorrelation Loss (SDL) to be optimised jointly with the other training objectives. Extensive experiments show that the proposed soft CCA is more effective and efficient than existing deep CCA models. In addition, our SDL loss can be applied to other deep models beyond multi-view learning, and obtains superior performance compared to existing decorrelation losses.


title: Topological Graph Persistence

abstract: Graphs are a basic tool for the representation of modern data. The richness of the topological information contained in a graph goes far beyond its mere interpretation as a one-dimensional simplicial complex. We show how topological constructions can be used to gain information otherwise concealed by the low-dimensional nature of graphs. We do that by extending previous work of other researchers in homological persistence, by proposing novel graph-theoretical constructions. Beyond cliques, we use independent sets, neighborhoods, enclaveless sets and a Ramsey-inspired extended persistence.


title: Anomalous eclipses of the young star RW Aur A

abstract: Results of UBVRIJHKLM photometry, VRI polarimetry and optical spectroscopy of a young star RW Aur A obtained during 2010-11 and 2014-16 dimming events are presented. During the second dimming the star decreased its brightness to \Delta V > 4.5 mag, polarization of its light in I-band was up to 30%, and color-magnitude diagram was similar to that of UX Ori type stars. We conclude that the reason of both dimmings is an eclipses of the star by dust screen, but the size of the screen is much larger than in the case of UXORs.


title: Asymptotic-Preserving Monte Carlo methods for transport equations in the diffusive limit

abstract: We develop a new Monte Carlo method that solves hyperbolic transport equations with stiff terms, characterized by a (small) scaling parameter. In particular, we focus on systems which lead to a reduced problem of parabolic type in the limit when the scaling parameter tends to zero. Classical Monte Carlo methods suffer of severe time step limitations in these situations, due to the fact that the characteristic speeds go to infinity in the diffusion limit. This makes the problem a real challenge, since the scaling parameter may differ by several orders of magnitude in the domain. To circumvent these time step limitations, we construct a new, asymptotic-preserving Monte Carlo method that is stable independently of the scaling parameter and degenerates to a standard probabilistic approach for solving the limiting equation in the diffusion limit. The method uses an implicit time discretization to formulate a modified equation in which the characteristic speeds do not grow indefinitely when the scaling factor tends to zero. The resulting modified equation can readily be discretized by a Monte Carlo scheme, in which the particles combine a finite propagation speed with a time-step dependent diffusion term. We show the performance of the method by comparing it with standard (deterministic) approaches in the literature.