title: Turing-like patterns in an asymmetric dynamic Ising model

abstract: To investigate novel aspects of pattern formation in spin systems, we use a mapping between reactive concentrations in a reaction-diffusion system and spin orientations in a dynamic multiple-spin Ising model. While pattern formation in Ising models always rely on infinite-range interactions, this mapping allows us to design an finite-range interactions Ising model that can produce patterns observed in reaction diffusion systems including Turing patterns with a tunable typical length scale. This model has asymmetric interactions and several spin types coexisting at a site. While we use the example of genetic regulation during embryo-genesis to build our model, it can be used to study the behavior of other complex systems of interacting agents.


title: Safe and adaptive importance sampling: a mixture approach

abstract: This paper investigates adaptive importance sampling algorithms for which the policy, the sequence of distributions used to generate the particles, is a mixture distribution between a flexible kernel density estimate (based on the previous particles), and a "safe" heavy-tailed density. When the share of samples generated according to the safe density goes to zero but not too quickly, two results are established: (i) uniform convergence rates are derived for the policy toward the target density; (ii) a central limit theorem is obtained for the resulting integral estimates. The fact that the asymptotic variance is the same as the variance of an "oracle" procedure with variance-optimal policy, illustrates the benefits of the approach. In addition, a subsampling step (among the particles) can be conducted before constructing the kernel estimate in order to decrease the computational effort without altering the performance of the method. The practical behavior of the algorithms is illustrated in a simulation study.


title: An ordered set of power-efficiency trade-offs

abstract: In this paper, we derive a number of inequalities which express power-efficiency trade-offs that hold generally for thermodynamic machines operating in non-equilibrium stationary states. One of these inequalities concerns the output power, which is bounded by a quadratic function of the thermodynamic efficiency multiplied by a factor. Different factors can be obtained according to the level of knowledge one has about the underlying dynamics of the machine, they can depend for instance on the covariance of the input flux, the dynamical activity, or the non-equilibrium conductance.


title: A Bayesian Nonparametric Approach for Evaluating the Causal Effect of Treatment in Randomized Trials with Semi-Competing Risks

abstract: We develop a Bayesian nonparametric (BNP) approach to evaluate the causal effect of treatment in a randomized trial where a nonterminal event may be censored by a terminal event, but not vice versa (i.e., semi-competing risks). Based on the idea of principal stratification, we define a novel estimand for the causal effect of treatment on the nonterminal event. We introduce identification assumptions, indexed by a sensitivity parameter, and show how to draw inference using our BNP approach. We conduct simulation studies and illustrate our methodology using data from a brain cancer trial.


title: Topological Data Analysis in Information Space

abstract: Various kinds of data are routinely represented as discrete probability distributions. Examples include text documents summarized by histograms of word occurrences and images represented as histograms of oriented gradients. Viewing a discrete probability distribution as a point in the standard simplex of the appropriate dimension, we can understand collections of such objects in geometric and topological terms. Importantly, instead of using the standard Euclidean distance, we look into dissimilarity measures with information-theoretic justification, and we develop the theory needed for applying topological data analysis in this setting. In doing so, we emphasize constructions that enable usage of existing computational topology software in this context.


title: Early Detection of Mental Stress Using Advanced Neuroimaging and Artificial Intelligence

abstract: While different neuroimaging modalities have been proposed to detect mental stress, each modality experiences certain limitations. This study proposed novel approaches to detect stress based on fusion of EEG and fNIRS signals in the feature-level using joint independent component analysis (jICA) and canonical correlation analysis method (CCA) and predected the level of stress using machine-learning approach. The jICA and CCA were then developed to combine the features to detect mental stress. The jICA fusion scheme discovers relationships between modalities by utilizing ICA to identify sources from each modality that modulate in the same way across subjects. The CCA fuse information from two sets of features to discover the associations across modalities and to ultimately estimate the sources responsible for these associations. The study further explored the functional connectivity (FC) and evaluated the performance of the fusion methods based on their classification performance and compared it with the result obtained by each individual modality. The jICA fusion technique significantly improved the classification accuracy, sensitivity and specificity on average by +3.46% compared to the EEG and +11.13% compared to the fNIRS. Similarly, CCA method improved the classification accuracy, sensitivity and specificity on average by +8.56% compared to the EEG and +13.03% compared to the fNIRS, respectively. The overall performance of the proposed fusion methods significantly improved the detection rate of mental stress, p<0.05. The FC significantly reduced under stress and suggested EEG and fNIRS as a potential biomarker of stress.


title: Spin dynamics of antiferromagnetically coupled ferromagnetic bilayers -- the case of Cr$_2$WO$_6$ and Cr$_2$MoO$_6$

abstract: Recent inelastic neutron diffraction measurements on Cr$_2$(Te, W, Mo)O$_6$ have revealed that these systems consist of bilayers of spin-3/2 Cr$^{3+}$ ions with strong antiferromagnetic inter-bilayer coupling and tuneable intra-bilayer coupling from ferro (for W and Mo) to antiferro (for Te). These measurements have determined the ground state spin structure and the values of sublattice magnetization, which shows significant reduction of sublattice magnetization from the atomic spin value of 3.0$\mu_B$ for Cr$^{3+}$ atoms. In an earlier paper we theoretically investigated the low temperature spin dynamics of Cr$_2$TeO$_6$ bilayer system where both the intra and inter-bilayer couplings are antiferromagnetic. In this paper we investigate Cr$_2$WO$_6$ and Cr$_2$MoO$_6$ systems where intra-bilayer exchange couplings are ferromagnetic but the inter-bilayer exchange couplings are antiferromagnetic. We obtain the magnon dispersion, sublattice magnetization, two-magnon density of states, longitudinal spin-spin correlation function, and its powder average and compare the results for these systems with results for Cr$_2$TeO$_6$.


title: Real order (an)-isotropic total variation in image processing - Part II: Learning of optimal structures

abstract: This article is the second work in our series of papers dedicated to image processing models based on the fractional order total variation $TV^r$. In our first work of this series, we studied key analytic properties of these semi-norms. Here we focus on the more applied aspects of such models: first, in order to obtain a better reconstructed image, we propose several extensions of the fractional order total variation. Such generalizations, collectively denoted by RVL, will be modular, i.e. the parameters therein are mutually independent, and can be fine tuned to the particular task. Then, we will study the bilevel training schemes based on RVL, and show that such schemes are well defined, i.e. they admit minimizers. Finally, we will provide some numerical examples, showing that training schemes based on RVL are effectively better than those based on classical regularizer $TV$ and $TGV^2$.


title: A Novel Monocular Disparity Estimation Network with Domain Transformation and Ambiguity Learning

abstract: Convolutional neural networks (CNN) have shown state-of-the-art results for low-level computer vision problems such as stereo and monocular disparity estimations, but still, have much room to further improve their performance in terms of accuracy, numbers of parameters, etc. Recent works have uncovered the advantages of using an unsupervised scheme to train CNN's to estimate monocular disparity, where only the relatively-easy-to-obtain stereo images are needed for training. We propose a novel encoder-decoder architecture that outperforms previous unsupervised monocular depth estimation networks by (i) taking into account ambiguities, (ii) efficient fusion between encoder and decoder features with rectangular convolutions and (iii) domain transformations between encoder and decoder. Our architecture outperforms the Monodepth baseline in all metrics, even with a considerable reduction of parameters. Furthermore, our architecture is capable of estimating a full disparity map in a single forward pass, whereas the baseline needs two passes. We perform extensive experiments to verify the effectiveness of our method on the KITTI dataset.


title: The Steep-Bounce Zeta Map in Parabolic Cataland

abstract: As a classical object, the Tamari lattice has many generalizations, including $\nu$-Tamari lattices and parabolic Tamari lattices. In this article, we unify these generalizations in a bijective fashion. We first prove that parabolic Tamari lattices are isomorphic to $\nu$-Tamari lattices for bounce paths $\nu$. We then introduce a new combinatorial object called `left-aligned colorable tree', and show that it provides a bijective bridge between various parabolic Catalan objects and certain nested pairs of Dyck paths. As a consequence, we prove the Steep-Bounce Conjecture using a generalization of the famous zeta map in $q,t$-Catalan combinatorics. A generalization of the zeta map on parking functions, which arises in the theory of diagonal harmonics, is also obtained as a labeled version of our bijection.


title: Implications of chiral symmetry on $S$-wave pionic resonances and the scalar charmed mesons

abstract: The chiral symmetry of QCD requires energy-dependent pionic strong interactions at low energies. This constraint, however, is not fulfilled by the usual Breit--Wigner parameterization of pionic resonances, leading to masses larger than the real ones. We derive relations between nonleptonic three-body decays of the $B$-meson into a $D$-meson and a pair of light pseudoscalar mesons based on SU(3) chiral symmetry. Employing effective field theory methods, we demonstrate that taking into account the final-state interactions, the experimental data of the decays $B^-\to D^+\pi^-\pi^-$, $B_s^0\to \bar{D}^0K^-\pi^+$, $B^0\to\bar{D}^0\pi^-\pi^+$, $B^-\to D^+\pi^-K^-$ and $B^0\to\bar{D}^0\pi^-K^+$ can all be described by the nonperturbative $\pi/\eta/K$-$D/D_s$ scattering amplitudes previously obtained from a combination of chiral effective field theory and lattice QCD calculations. The results provide a strong support of the scenario that the broad scalar charmed meson $D^\ast_0(2400)$ should be replaced by two states, the lower one of which has a mass of around 2.1 GeV, much smaller than that extracted from experimental data using a Breit--Wigner parameterization.


title: Estimating the density of resident coastal fish using underwater cameras: accounting for individual detectability

abstract: Technological advances in underwater video recording are opening novel opportunities for monitoring wild fish. However, extracting data from videos is often challenging. Nevertheless, it has been recently demonstrated that accurate and precise estimates of density for animals (whose normal activities are restricted to a bounded area or home range) can be obtained from counts averaged across a relatively low number of video frames. The method, however, requires that individual detectability (PID, the probability of detecting a given animal provided that it is actually within the area surveyed by a camera) has to be known. Here we propose a Bayesian implementation for estimating PID after combining counts from cameras with counts from any reference method. The proposed framework was demonstrated using Serranus scriba as a case-study, a widely distributed and resident coastal fish. Density and PID were calculated after combining fish counts from unbaited remote underwater video (RUV) and underwater visual censuses (UVC) as reference method. The relevance of the proposed framework is that after estimating PID, fish density can be estimated accurately and precisely at the UVC scale (or at the scale of the preferred reference method) using RUV only. This key statement has been extensively demonstrated using computer simulations yielded by real empirical data. Finally, we provide a simulation tool-kit for comparing the expected precision attainable for different sampling effort and for species with different levels of PID. Overall, the proposed method may contribute to substantially enlarge the spatio-temporal scope of density monitoring programs for many resident fish.


title: Suppressing cosmic variance with paired-and-fixed cosmological simulations: average properties and covariances of dark matter clustering statistics

abstract: Making cosmological inferences from the observed galaxy clustering requires accurate predictions for the mean clustering statistics and their covariances. Those are affected by cosmic variance -- the statistical noise due to the finite number of harmonics. The cosmic variance can be suppressed by fixing the amplitudes of the harmonics instead of drawing them from a Gaussian distribution predicted by the inflation models. Initial realizations also can be generated in pairs with 180 degrees flipped phases to further reduce the variance. Here, we compare the consequences of using paired-and-fixed vs Gaussian initial conditions on the average dark matter clustering and covariance matrices predicted from N-body simulations. As in previous studies, we find no measurable differences between paired-and-fixed and Gaussian simulations for the average density distribution function, power spectrum and bispectrum. Yet, the covariances from paired-and-fixed simulations are suppressed in a complicated scale- and redshift-dependent way. The situation is particularly problematic on the scales of Baryon Acoustic Oscillations where the covariance matrix of the power spectrum is lower by only 20% compared to the Gaussian realizations, implying that there is not much of a reduction of the cosmic variance. The non-trivial suppression, combined with the fact that paired-and-fixed covariances are noisier than from Gaussian simulations, suggests that there is no path towards obtaining accurate covariance matrices from paired-and-fixed simulations. Because the covariances are crucial for the observational estimates of galaxy clustering statistics and cosmological parameters, paired-and-fixed simulations, though useful for some applications, cannot be used for the production of mock galaxy catalogs.


title: Topology-based Representative Datasets to Reduce Neural Network Training Resources

abstract: One of the main drawbacks of the practical use of neural networks is the long time required in the training process. Such a training process consists of an iterative change of parameters trying to minimize a loss function. These changes are driven by a dataset, which can be seen as a set of labelled points in an n-dimensional space. In this paper, we explore the concept of are representative dataset which is a dataset smaller than the original one, satisfying a nearness condition independent of isometric transformations. Representativeness is measured using persistence diagrams (a computational topology tool) due to its computational efficiency. We prove that the accuracy of the learning process of a neural network on a representative dataset is "similar" to the accuracy on the original dataset when the neural network architecture is a perceptron and the loss function is the mean squared error. These theoretical results accompanied by experimentation open a door to reducing the size of the dataset to gain time in the training process of any neural network.


title: A control problem related to the parabolic dominative $p$-Laplace equation

abstract: We show that value functions of a certain time-dependent control problem in $\Omega\times (0,T)$, with a continuous payoff $F$ on the parabolic boundary, converge uniformly to the viscosity solution of the parabolic dominative $p$-Laplace equation $$2(n+p)u_t=\Delta u+(p-2)\lambda_n(D^2 u),$$ with the boundary data $F$. Here $2\leq p< \infty$, and $\lambda_n(D^2 u)$ is the largest eigenvalue of the Hessian $D^2 u$.


title: The solar fast dynamo action extended to A, Ap and Am type stars

abstract: Initially, we present our original fast dynamo mechanism being potentially at work concerning the Sun and the solar type stars. Then, based on this new prototype model, a dynamo action model explaining the generation of the magnetic fields for A, Ap and Am type stars will be derived. We argue that their magnetic behavior is essentially due to locally exotic plasma properties leading to large scale charge separation in their interior. Their powerful magnetic fields are generated by the revolution speed of these excess positive or negative charges tightly associated with the formation of one or two so called Torus structures. The lifespan of the Torus charge defines the stellar activity cycle. In the solar case, two Torus structures are formed preserving their charge for 11 years. The Torus charging results from two processes: First, the stellar differential rotation with its shear layers and second, the local plasma property that the plasma rotational speed (U) is equal or greater than the local light speed (v) in the medium. If U is greater than v, then two nearby pocket charges (initially formed by the stellar differential rotation) will be mutually attracted. That is, the attractive magnetic force will be greater than the repulsive electric one. In addition, the same poloidal field, in the subphotosphere shear layer, produces an excess charge generating the starspot activity. The solar birthplace of Torus, within the tachocline, is positioned at 45 deg. latitudes. Each Torus drifts equatorward increasing steadily its own rotational speed and, consequently, its own strength of poloidal magnetic field. A single Torus, over the equatorial plane, is probably formed in the interior of an Ap type star producing its strong and global dipole like magnetic field.


title: Limits of flexural wave absorption by open lossy resonators: reflection and transmission problems

abstract: The limits of flexural wave absorption by open lossy resonators are analytically and numerically reported in this work for both the reflection and transmission problems. An experimental validation for the reflection problem is presented. The reflection and transmission of flexural waves in 1D resonant thin beams are analyzed by means of the transfer matrix method. The hypotheses, on which the analytical model relies, are validated by experimental results. The open lossy resonator, consisting of a finite length beam thinner than the main beam, presents both energy leakage due to the aperture of the resonators to the main beam and inherent losses due to the viscoelastic damping. Wave absorption is found to be limited by the balance between the energy leakage and the inherent losses of the open lossy resonator. The perfect compensation of these two elements is known as the critical coupling condition and can be easily tuned by the geometry of the resonator. On the one hand, the scattering in the reflection problem is represented by the reflection coefficient. A single symmetry of the resonance is used to obtain the critical coupling condition. Therefore the perfect absorption can be obtained in this case. On the other hand, the transmission problem is represented by two eigenvalues of the scattering matrix, representing the symmetric and anti-symmetric parts of the full scattering problem. In the geometry analyzed in this work, only one kind of symmetry can be critically coupled, and therefore, the maximal absorption in the transmission problem is limited to 0.5. The results shown in this work pave the way to the design of resonators for efficient flexural wave absorption.


title: Ontology of Card Sleights

abstract: We present a machine-readable movement writing for sleight-of-hand moves with cards -- a "Labanotation of card magic." This scheme of movement writing contains 440 categories of motion, and appears to taxonomize all card sleights that have appeared in over 1500 publications. The movement writing is axiomatized in $\mathcal{SROIQ}$(D) Description Logic, and collected formally as an Ontology of Card Sleights, a computational ontology that extends the Basic Formal Ontology and the Information Artifact Ontology. The Ontology of Card Sleights is implemented in OWL DL, a Description Logic fragment of the Web Ontology Language. While ontologies have historically been used to classify at a less granular level, the algorithmic nature of card tricks allows us to transcribe a performer's actions step by step. We conclude by discussing design criteria we have used to ensure the ontology can be accessed and modified with a simple click-and-drag interface. This may allow database searches and performance transcriptions by users with card magic knowledge, but no ontology background.


title: Robot mirroring: A framework for self-tracking feedback through empathy with an artificial agent representing the self

abstract: Current technologies have enabled us to track and quantify our physical state and behavior. Self-tracking aims to achieve increased awareness to decrease undesired behaviors and lead to a healthier lifestyle. However, inappropriately communicated self-tracking results might cause the opposite effect. In this work, we propose a subtle self-tracking feedback by mirroring the self's state into an artificial agent. By eliciting empathy towards the artificial agent and fostering helping behaviors, users would help themselves as well. Finally, we reflected on the implications of this design framework, and the methodology to design and implement it. A series of interviews to expert designers pointed out to the importance of having multidisciplinary teams working in parallel. Moreover, an agile methodology with a sprint zero for the initial design, and shifted user research, design, and implementation sprints were proposed. Similar systems with data flow and hardware dependencies would also benefit from the proposed agile design process.


title: Interplay between Loewner and Dirichlet energies via conformal welding and flow-lines

abstract: The Loewner energy of a Jordan curve is the Dirichlet energy of its Loewner driving term. It is finite if and only if the curve is a Weil-Petersson quasicircle. In this paper, we describe cutting and welding operations on finite Dirichlet energy functions defined in the plane, allowing expression of the Loewner energy in terms of Dirichlet energy dissipation. We show that the Loewner energy of a unit vector field flow-line is equal to the Dirichlet energy of the harmonically extended winding. We also give an identity involving a complex-valued function of finite Dirichlet energy that expresses the welding and flow-line identities simultaneously. As applications, we prove that arclength isometric welding of two domains is sub-additive in the energy, and that the energy of equipotentials in a simply connected domain is monotone. Our main identities can be viewed as action functional analogs of both the welding and flow-line couplings of Schramm-Loewner evolution curves with the Gaussian free field.


title: On the convergence of a finite volume method for the Navier-Stokes-Fourier system

abstract: We study convergence of a finite volume scheme for the Navier-Stokes-Fourier system describing the motion of compressible viscous and heat conducting fluids. The numerical flux uses upwinding with an additional numerical diffusion of order $\mathcal{O} (h^{ \varepsilon+1})$, $0<\varepsilon<1$. The approximate solutions are piecewise constant functions with respect to the underlying mesh. We show that any uniformly bounded sequence of numerical solutions converges unconditionally to the solution of the Navier-Stokes-Fourier system. In particular, the existence of the solution to the Navier-Stokes-Fourier system is not a priori assumed.


title: Accurate 3D Face Reconstruction with Weakly-Supervised Learning: From Single Image to Image Set

abstract: Recently, deep learning based 3D face reconstruction methods have shown promising results in both quality and efficiency.However, training deep neural networks typically requires a large volume of data, whereas face images with ground-truth 3D face shapes are scarce. In this paper, we propose a novel deep 3D face reconstruction approach that 1) leverages a robust, hybrid loss function for weakly-supervised learning which takes into account both low-level and perception-level information for supervision, and 2) performs multi-image face reconstruction by exploiting complementary information from different images for shape aggregation. Our method is fast, accurate, and robust to occlusion and large pose. We provide comprehensive experiments on three datasets, systematically comparing our method with fifteen recent methods and demonstrating its state-of-the-art performance.


title: On a model of axisymmetric flows in a free boundary domain

abstract: We consider a model of axisymmetric flows for a free boundary vortex embedded in a statically stable fluid at rest. We identify the boundary of the vortex by solving a variational problem. Then, we reduce the analysis of the dynamics of the vortex to the study of a class of continuity equations for which we construct a solution.


title: Computational modeling of active deformable membranes embedded in 3D flows

abstract: Active gel theory has recently been very successful in describing biologically active materials such as actin filaments or moving bacteria in temporally fixed and simple geometries such as cubes or spheres. Here we develop a computational algorithm to compute the dynamic evolution of an arbitrarily shaped, deformable thin membrane of active material embedded in a 3D flowing liquid. For this, our algorithm combines active gel theory with the classical theory of thin elastic shells. To compute the actual forces resulting from active stresses, we apply a parabolic fitting procedure to the triangulated membrane surface. Active forces are then dynamically coupled via an Immersed-Boundary method to the surrounding fluid whose dynamics can be solved by any standard, e.g. Lattice-Boltzmann, flow solver. We validate our algorithm using the Green's functions of [Berthoumieux et al. New J. Phys. 16, 065005 (2014)] for an active cylindrical membrane subjected (i) to a locally increased active stress and (ii) to a homogeneous active stress. For the latter scenario, we predict in addition a so far unobserved non-axisymmetric instability. We highlight the versatility of our method by analyzing the flow field inside an actively deforming cell embedded in external shear flow. Further applications may be cytoplasmic streaming or active membranes in blood flows.


title: Periodic Dirac operator on the half-line

abstract: We consider the Dirac operator with a periodic potential on the half-line with the Dirichlet boundary condition at zero. Its spectrum consists of an absolutely continuous part plus at most one eigenvalue in each open gap. The Dirac resolvent admits a meromorphic continuation onto a two-sheeted Riemann surface with a unique simple pole on each open gap: on the first sheet (an eigenvalue) or on the second sheet (a resonance). These poles are called states and there are no other poles. If the potential is shifted by real parameter t, then the continuous spectrum does not change but the states can change theirs positions. We prove that each state is smooth and in general, non-monotonic function of t. We prove that a state is a strictly monotone function of t for a specific potential. Using these results we obtain formulas to recover potentials of special forms.


title: Prioritized Process Test: An Alternative to Current Process Testing Strategies

abstract: Testing processes and workflows in information and Internet of Things systems is a major part of the typical software testing effort. Consistent and efficient path-based test cases are desired to support these tests. Because certain parts of software system workflows have a higher business priority than others, this fact has to be involved in the generation of test cases. In this paper, we propose a Prioritized Process Test (PPT), which is a model-based test case generation algorithm that represents an alternative to currently established algorithms that use directed graphs and test requirements to model the system under test. The PPT accepts a directed multigraph as a model to express priorities, and edge weights are used instead of test requirements. To determine the test-coverage level of test cases, a test-depth-level concept is used. We compared the presented PPT with five alternatives (i.e., the Process Cycle Test, a naive reduction of test set created by the Process Cycle Test, Brute Force algorithm, Set-covering Based Solution and Matching-based Prefix Graph Solution) for edge coverage and edge-pair coverage. To assess the optimality of the path-based test cases produced by these strategies, we used fourteen metrics based on the properties of these test cases and 59 models that were created for three real-world systems. For all edge coverage, the PPT produced more optimal test cases than the alternatives in terms of the majority of the metrics. For edge-pair coverage, the PPT strategy yielded similar results to those of the alternatives. Thus, the PPT strategy is an applicable alternative, as it reflects both the required test coverage level and the business priority in parallel.


title: Spaces of small cellularity have nowhere constant continuous images of small weight

abstract: We call a continuous map $f : X \to Y$ nowhere constant if it is not constant on any non-empty open subset of its domain $X$. Clearly, this is equivalent with the assumption that every fiber $f^{-1}(y)$ of $f$ is nowhere dense in $X$. We call the continuous map $f : X \to Y$ pseudo-open if for each nowhere dense $Z \subset Y$ its inverse image $f^{-1}(Z)$ is nowhere dense in $X$. Clearly, if $Y$ is crowded, i.e. has no isolated points, then $f$ is nowhere constant. The aim of this paper is to study the following, admittedly imprecise, question: How "small" nowhere constant, resp. pseudo-open continuous images can "large" spaces have? Our main results yield the following two precise answers to this question, explaining also our title. Both of them involve the cardinal function $\widehat{c}(X)$, the "hat version" of cellularity, which is defined as the smallest cardinal $\kappa$ such that there is no $\kappa$-sized disjoint family of open sets in $X$. Thus, for instance, $\widehat{c}(X) = \omega_1$ means that $X$ is CCC. THEOREM A. Any crowded Tychonov space $X$ has a crowded Tychonov nowhere constant continuous image $Y$ of weight $w(Y) \le \widehat{c}(X)$. Moreover, in this statement $\le$ may be replaced with $<$ iff there are no $\widehat{c}(X)$-Suslin lines (or trees). THEOREM B. Any crowded Tychonov space $X$ has a crowded Tychonov pseudo-open continuous image $Y$ of weight $w(Y) \le 2^{<\widehat{c}(X)}$. If Martin's axiom holds then there is a CCC crowded Tychonov space $X$ such that for any crowded Hausdorff pseudo-open continuous image $Y$ of $X$ we have $w(Y) \ge \mathfrak{c}\,( = 2^{< \omega_1})$.


title: Competition between Phase Separation and Spin Density Wave or Charge Density Wave Order: Role of Long-Range Interactions

abstract: Recent studies of pairing and charge order in materials such as FeSe, SrTiO$_3$, and 2H-NbSe$_2$ have suggested that momentum dependence of the electron-phonon coupling plays an important role in their properties. Initial attempts to study Hamiltonians which either do not include or else truncate the range of Coulomb repulsion have noted that the resulting spatial non-locality of the electron-phonon interaction leads to a dominant tendency to phase separation. Here we present Quantum Monte Carlo results for such models in which we incorporate both on-site and intersite electron-electron interactions. We show that these can stabilize phases in which the density is homogeneous and determine the associated phase boundaries. As a consequence, the physics of momentum dependent electron-phonon coupling can be determined outside of the trivial phase separated regime.


title: Homogenization of time-fractional diffusion equations with periodic coefficients

abstract: We consider the initial boundary value problem for the time-fractional diffusion equation with a homogeneous Dirichlet boundary condition and an inhomogeneous initial data $a(x)\in L^{2}(D)$ in a bounded domain $D\subset \mathbb{R}^d$ with a sufficiently smooth boundary. We analyze the homogenized solution under the assumption that the diffusion coefficient $\kappa^{\epsilon}(x)$ is smooth and periodic with the period $\epsilon>0$ being sufficiently small. We derive that its first order approximation has a convergence rate of $\mathcal{O}(\epsilon^{1/2})$ when the dimension $d\leq 2$ and $\mathcal{O}(\epsilon^{1/6})$ when $d=3$. Several numerical tests are presented to show the performance of the first order approximation.


title: The VHE Gamma-Ray View of the FSRQ PKS~1510-089

abstract: The flat spectrum radio quasar PKS 1510-089 is a monitored target in many wavelength bands due to its high variability. It was detected as a very-high-energy (VHE) $\gamma$-ray emitter with H.E.S.S. in 2009, and has since been a regular target of VHE observations by the imaging Cherenkov observatories H.E.S.S. and MAGIC. In this paper, we summarize the current state of results focusing on the monitoring effort with H.E.S.S. and the discovery of a particularly strong VHE flare in 2016 with H.E.S.S. and MAGIC. While the source has now been established as a weak, but regular emitter at VHE, no correlation with other energy bands has been established. This is underlined by the 2016 VHE flare, where the detected optical and high-energy $\gamma$-ray counterparts evolve differently than the VHE flux.


title: Segmentation-Based Deep-Learning Approach for Surface-Defect Detection

abstract: Automated surface-anomaly detection using machine learning has become an interesting and promising area of research, with a very high and direct impact on the application domain of visual inspection. Deep-learning methods have become the most suitable approaches for this task. They allow the inspection system to learn to detect the surface anomaly by simply showing it a number of exemplar images. This paper presents a segmentation-based deep-learning architecture that is designed for the detection and segmentation of surface anomalies and is demonstrated on a specific domain of surface-crack detection. The design of the architecture enables the model to be trained using a small number of samples, which is an important requirement for practical applications. The proposed model is compared with the related deep-learning methods, including the state-of-the-art commercial software, showing that the proposed approach outperforms the related methods on the specific domain of surface-crack detection. The large number of experiments also shed light on the required precision of the annotation, the number of required training samples and on the required computational cost. Experiments are performed on a newly created dataset based on a real-world quality control case and demonstrates that the proposed approach is able to learn on a small number of defected surfaces, using only approximately 25-30 defective training samples, instead of hundreds or thousands, which is usually the case in deep-learning applications. This makes the deep-learning method practical for use in industry where the number of available defective samples is limited. The dataset is also made publicly available to encourage the development and evaluation of new methods for surface-defect detection.


title: Microscopically based energy density functionals for nuclei using the density matrix expansion

abstract: While ab initio many-body techniques have been able to successfully describe the properties of light and intermediate mass nuclei based on chiral effective field theory interactions, neutron-rich nuclei still remain out of reach for these methods. Conversely, energy density functional approaches can be used to calculate properties of heavy nuclei but rely mostly on phenomenological interactions. A usable form of the nuclear energy density functional that is rooted in the modern theory of nuclear forces was presented recently. The first component of this new set of functionals corresponds to the direct part (Hartree term) of the expectation value of local chiral potentials on a Slater determinant. The exchange term, which is a functional of the non-local density, is transformed into a local functional by applying the density matrix expansion. In order to reduce the computational cost due to the direct implementation of non-separable, local interactions in the Hartree term, we use an approximation to represent the regularized Yukawa functions in terms of a sum of (separable) Gaussian functions. These proceedings analyze the accuracy of such an approximation in terms of the number of Gaussian functions and look for an optimal value that gives an acceptable level of accuracy while maintaining the computational memory requirements in a many-body calculation as low as possible.


title: On symmetrizing the ultraspherical spectral method for self-adjoint problems

abstract: A mechanism is described to symmetrize the ultraspherical spectral method for self-adjoint problems. The resulting discretizations are symmetric and banded. An algorithm is presented for an adaptive spectral decomposition of self-adjoint operators. Several applications are explored to demonstrate the properties of the symmetrizer and the adaptive spectral decomposition.


title: Alexandrov geometry: foundations

abstract: Alexandrov spaces are defined via axioms similar to those given by Euclid. The Alexandrov axioms replace certain equalities with inequalities. Depending on the signs of the inequalities, we obtain Alexandrov spaces with curvature bounded above and curvature bounded below. The definitions of the two classes of spaces are similar, but their properties and known applications are quite different. Our approach is novel in its attention to the interrelatedness of the two fields, and its emphasis on the way each illuminates the other. The goal of this book is to give a comprehensive exposition of the structure theory of Alexandrov spaces with curvature bounded above and below. It includes all the basic material as well as selected topics inspired by considering the two contexts simultaneously. We only consider the intrinsic theory, leaving applications aside. This book includes material up to the definition of dimension. Another volume still in preparation will cover further topics.


title: On the asymptotics of the rescaled Appell polynomials

abstract: We introduce a new representation for the rescaled Appell polynomials and use it to obtain asymptotic expansions to arbitrary order. This representation consists of a finite sum and an integral over a universal contour (i.e. independent of the particular polynomials considered within the Appell family). We illustrate our method by studying the zero attractors for rescaled Appell polynomials. We also discuss the asymptotics to arbitrary order of the rescaled Bernoulli polynomials.


title: High performance photonic microwave filters based on a 50GHz optical soliton crystal Kerr micro-comb

abstract: We demonstrate a photonic radio frequency (RF) transversal filter based on an integrated optical micro-comb source featuring a record low free spectral range of 49 GHz yielding 80 micro-comb lines across the C-band. This record-high number of taps, or wavelengths for the transversal filter results in significantly increased performance including a QRF factor more than four times higher than previous results. Further, by employing both positive and negative taps, an improved out-of-band rejection of up to 48.9 dB is demonstrated using Gaussian apodization, together with a tunable centre frequency covering the RF spectra range, with a widely tunable 3-dB bandwidth and versatile dynamically adjustable filter shapes. Our experimental results match well with theory, showing that our transversal filter is a competitive solution to implement advanced adaptive RF filters with broad operational bandwidths, high frequency selectivity, high reconfigurability, and potentially reduced cost and footprint. This approach is promising for applications in modern radar and communications systems.


title: Learning Gentle Object Manipulation with Curiosity-Driven Deep Reinforcement Learning

abstract: Robots must know how to be gentle when they need to interact with fragile objects, or when the robot itself is prone to wear and tear. We propose an approach that enables deep reinforcement learning to train policies that are gentle, both during exploration and task execution. In a reward-based learning environment, a natural approach involves augmenting the (task) reward with a penalty for non-gentleness, which can be defined as excessive impact force. However, augmenting with only this penalty impairs learning: policies get stuck in a local optimum which avoids all contact with the environment. Prior research has shown that combining auxiliary tasks or intrinsic rewards can be beneficial for stabilizing and accelerating learning in sparse-reward domains, and indeed we find that introducing a surprise-based intrinsic reward does avoid the no-contact failure case. However, we show that a simple dynamics-based surprise is not as effective as penalty-based surprise. Penalty-based surprise, based on predicting forceful contacts, has a further benefit: it encourages exploration which is contact-rich yet gentle. We demonstrate the effectiveness of the approach using a complex, tendon-powered robot hand with tactile sensors. Videos are available at http://sites.google.com/view/gentlemanipulation.


title: Optimizing thermodynamic trajectories using evolutionary and gradient-based reinforcement learning

abstract: Using a model heat engine, we show that neural network-based reinforcement learning can identify thermodynamic trajectories of maximal efficiency. We consider both gradient and gradient-free reinforcement learning. We use an evolutionary learning algorithm to evolve a population of neural networks, subject to a directive to maximize the efficiency of a trajectory composed of a set of elementary thermodynamic processes; the resulting networks learn to carry out the maximally-efficient Carnot, Stirling, or Otto cycles. When given an additional irreversible process, this evolutionary scheme learns a previously unknown thermodynamic cycle. Gradient-based reinforcement learning is able to learn the Stirling cycle, whereas an evolutionary approach achieves the optimal Carnot cycle. Our results show how the reinforcement learning strategies developed for game playing can be applied to solve physical problems conditioned upon path-extensive order parameters.

categories: cs.NE cond-mat.stat-mech cs.LG physics.comp-ph


title: Beam-helicity asymmetries for single-hadron production in semi-inclusive deep-inelastic scattering from unpolarized hydrogen and deuterium targets

abstract: A measurement of beam-helicity asymmetries for single-hadron production in deep-inelastic scattering is presented. Data from the scattering of 27.6 GeV electrons and positrons off gaseous hydrogen and deuterium targets were collected by the HERMES experiment. The asymmetries are presented separately as a function of the Bjorken scaling variable, the hadron transverse momentum, and the fractional energy for charged pions and kaons as well as for protons and anti-protons. These asymmetries are also presented as a function of the three aforementioned kinematic variables simultaneously.


title: Quantum information with general quantum variables: a formalism encompassing qubits, qudits, and quantum continuous variables

abstract: This note presents a simple and unified formulation of the most fundamental structures used in quantum information with qubits, arbitrary dimension qudits, and quantum continuous variables. This \emph{general quantum variables} construction provides a succinct language for formulating many results in quantum computation and information so that they are applicable in all dimensions. The structures included within this formalism include: a generalization to arbitrary dimension of the three Pauli operators, and the associated mutually unbiased bases; the Pauli and Clifford groups; many important quantum gates; standard sets of generators for the Clifford group; and simple universal gate sets. This formalism provides a convenient, intuitive and extensible language for easily generalizing results that were originally derived for a single type of system (often qubits or quantum continuous variables), and that rely on only those structures listed above, to apply in all dimensions.


title: How to observe the vacuum decay in low-energy heavy-ion collisions

abstract: In slow collisions of two bare nuclei with the total charge larger than the critical value $Z_{\rm cr} \approx 173$, the initially neutral vacuum can spontaneously decay into the charged vacuum and two positrons. Detection of the spontaneous emission of positrons would be the direct evidence of this fundamental phenomenon. However, the spontaneously produced particles are indistinguishable from the dynamical background in the positron spectra. We show that the vacuum decay can nevertheless be observed via impact-sensitive measurements of pair-production probabilities. Possibility of such observation is demonstrated using numerical calculations of pair production in low-energy collisions of heavy nuclei.


title: Attractive interactions between like-oriented surface steps from an ab initio perspective: Role of the elastic and electrostatic contributions

abstract: At low temperatures the interactions between like-oriented steps on a surface are believed to be dominated by elastic repulsions. This belief is based on the results of the classical continuum field theories of elasticity. Electrostatic interactions are also believed to have a role but no theory has yet been developed for their description. In our paper we critically revisit the assumptions of these theories, and found that the ab initio approaches, based on the Born-Oppenheimer distinction between the valence electrons and the ions degrees of freedom, allow to identify distinctly both elastic and electrostatic contributions to the step energy and the step interaction energy. By using ab-initio calculations and specifying to the case of the technologically important (001)surface of GaAs, we show that the elastic interactions are not always repulsive and attractive step interactions between like-oriented steps are also possible. Furthermore, the repulsiveness and attractiveness of the step interactions depend on the step atomic structure.


title: Learning Convolutional Transforms for Lossy Point Cloud Geometry Compression

abstract: Efficient point cloud compression is fundamental to enable the deployment of virtual and mixed reality applications, since the number of points to code can range in the order of millions. In this paper, we present a novel data-driven geometry compression method for static point clouds based on learned convolutional transforms and uniform quantization. We perform joint optimization of both rate and distortion using a trade-off parameter. In addition, we cast the decoding process as a binary classification of the point cloud occupancy map. Our method outperforms the MPEG reference solution in terms of rate-distortion on the Microsoft Voxelized Upper Bodies dataset with 51.5% BDBR savings on average. Moreover, while octree-based methods face exponential diminution of the number of points at low bitrates, our method still produces high resolution outputs even at low bitrates. Code and supplementary material are available at https://github.com/mauriceqch/pcc_geo_cnn .


title: Convolutional Sparse Coding for Compressed Sensing CT Reconstruction

abstract: Over the past few years, dictionary learning (DL)-based methods have been successfully used in various image reconstruction problems. However, traditional DL-based computed tomography (CT) reconstruction methods are patch-based and ignore the consistency of pixels in overlapped patches. In addition, the features learned by these methods always contain shifted versions of the same features. In recent years, convolutional sparse coding (CSC) has been developed to address these problems. In this paper, inspired by several successful applications of CSC in the field of signal processing, we explore the potential of CSC in sparse-view CT reconstruction. By directly working on the whole image, without the necessity of dividing the image into overlapped patches in DL-based methods, the proposed methods can maintain more details and avoid artifacts caused by patch aggregation. With predetermined filters, an alternating scheme is developed to optimize the objective function. Extensive experiments with simulated and real CT data were performed to validate the effectiveness of the proposed methods. Qualitative and quantitative results demonstrate that the proposed methods achieve better performance than several existing state-of-the-art methods.


title: OCGAN: One-class Novelty Detection Using GANs with Constrained Latent Representations

abstract: We present a novel model called OCGAN for the classical problem of one-class novelty detection, where, given a set of examples from a particular class, the goal is to determine if a query example is from the same class. Our solution is based on learning latent representations of in-class examples using a denoising auto-encoder network. The key contribution of our work is our proposal to explicitly constrain the latent space to exclusively represent the given class. In order to accomplish this goal, firstly, we force the latent space to have bounded support by introducing a tanh activation in the encoder's output layer. Secondly, using a discriminator in the latent space that is trained adversarially, we ensure that encoded representations of in-class examples resemble uniform random samples drawn from the same bounded space. Thirdly, using a second adversarial discriminator in the input space, we ensure all randomly drawn latent samples generate examples that look real. Finally, we introduce a gradient-descent based sampling technique that explores points in the latent space that generate potential out-of-class examples, which are fed back to the network to further train it to generate in-class examples from those points. The effectiveness of the proposed method is measured across four publicly available datasets using two one-class novelty detection protocols where we achieve state-of-the-art results.


title: Open quantum systems in Heisenberg picture

abstract: In the framework of the Heisenberg picture, an alternative derivation of the reduced density matrix of a driven dissipative quantum harmonic oscillator as the prototype of an open quantum system is investigated. The reduced density matrix for different initial states of the combined system is obtained from a general formula, and different limiting cases are studied. Exact expressions for the corresponding characteristic function in quantum thermodynamics and Wigner quasi distribution function are found. A possible generalization based on the Magnus expansion of the evolution operator is presented.


title: Traversing the noise of dynamic mini-batch sub-sampled loss functions: A visual guide

abstract: Mini-batch sub-sampling in neural network training is unavoidable, due to growing data demands, memory-limited computational resources such as graphical processing units (GPUs), and the dynamics of on-line learning. In this study we specifically distinguish between static mini-batch sub-sampled loss functions, where mini-batches are intermittently fixed during training, resulting in smooth but biased loss functions; and the dynamic sub-sampling equivalent, where new mini-batches are sampled at every loss evaluation, trading bias for variance in sampling induced discontinuities. These render automated optimization strategies such as minimization line searches ineffective, since critical points may not exist and function minimizers find spurious, discontinuity induced minima. This paper suggests recasting the optimization problem to find stochastic non-negative associated gradient projection points (SNN-GPPs). We demonstrate that the SNN-GPP optimality criterion is less susceptible to sub-sampling induced discontinuities than critical points or minimizers. We conduct a visual investigation, comparing local minimum and SNN-GPP optimality criteria in the loss functions of a simple neural network training problem for a variety of popular activation functions. Since SNN-GPPs better approximate the location of true optima, particularly when using smooth activation functions with high curvature characteristics, we postulate that line searches locating SNN-GPPs can contribute significantly to automating neural network training


title: Ultrafast band-gap renormalization and build-up of optical gain in monolayer MoTe$_2$

abstract: The dynamics of band-gap renormalization and gain build-up in monolayer MoTe$_2$ is investigated by evaluating the non-equilibrium Dirac-Bloch equations with the incoherent carrier-carrier and carrier-phonon scattering treated via quantum-Boltzmann type scattering equations. For the case where an approximately $300$ fs-long high intensity optical pulse generates charge-carrier densities in the gain regime, the strong Coulomb coupling leads to a relaxation of excited carriers on a few fs time scale. The pump-pulse generation of excited carriers induces a large band-gap renormalization during the time scale of the pulse. Efficient phonon coupling leads to a subsequent carrier thermalization within a few ps, which defines the time scale for the optical gain build-up energetically close to the low-density exciton resonance.


title: A local version of Einstein's formula for the effective viscosity of suspensions

abstract: We prove a local variant of Einstein's formula for the effective viscosity of dilute suspensions, that is $\mu^\prime=\mu (1+\frac 5 2\phi+o(\phi))$, where $\phi$ is the volume fraction of the suspended particles. Up to now rigorous justifications have only been obtained for dissipation functionals of the flow field. We prove that the formula holds on the level of the Stokes equation (with variable viscosity). We consider a regime where the number $N$ of particles suspended in the fluid goes to infinity while their size $R$ and the volume fraction $\phi=NR^3$ approach zero. We establish $L^\infty$ and $L^p$ estimates for the difference of the microscopic solution to the solution of the homogenized equation. Here we assume that the particles are contained in a bounded region and are well separated in the sense that the minimal distance is comparable to the average one. The main tools for the proof are a dipole approximation of the flow field of the suspension together with the so-called method of reflections and a coarse graining of the volume density.


title: Multiscale population dynamics in reproductive biology: singular perturbation reduction in deterministic and stochastic models

abstract: In this study, we describe different modeling approaches for ovarian follicle population dynamics, based on either ordinary (ODE), partial (PDE) or stochastic (SDE) differential equations, and accounting for interactions between follicles. We put a special focus on representing the population-level feedback exerted by growing ovarian follicles onto the activation of quiescent follicles. We take advantage of the timescale difference existing between the growth and activation processes to apply model reduction techniques in the framework of singular perturbations. We first study the linear versions of the models to derive theoretical results on the convergence to the limit models. In the nonlinear cases, we provide detailed numerical evidence of convergence to the limit behavior. We reproduce the main semi-quantitative features characterizing the ovarian follicle pool, namely a bimodal distribution of the whole population, and a slope break in the decay of the quiescent pool with aging.


title: Thermodynamics of Non-Commutative Scalar-Tensor-Vector Gravity Black holes

abstract: In this paper, we analyze the thermodynamic stability of Schwarzschild Modified Gravity (MOG) black holes in a non-commutative framework. We show that, unlike a commutative MOG black hole, in the coherent state picture of non-commutativity MOG black holes are thermodynamically stable. At the final stage of evaporation, a stable remnant with zero temperature and finite entropy is left in this noncommutative framework. Also, we consider Parikh-Wilczek tunneling mechanism of massive particles from non-commutative MOG black holes and demonstrate that information leaks out of non-commutative MOG black holes in the form of some non-thermal correlations.


title: Theory for Glassy Behavior of Supercooled Liquid Mixtures

abstract: We present a model for glassy dynamics in supercooled liquid mixtures. Given the relaxation behavior of individual supercooled liquids, the model predicts the relaxation times of their mixtures as temperature is decreased. The model is based on dynamical facilitation theory for glassy dynamics, which provides a physical basis for relaxation and vitrification of a supercooled liquid. This is in contrast to empirical linear interpolations such as the Gordon-Taylor equation typically used to predict glass transition temperatures of liquid mixtures. To understand the behavior of supercooled liquid mixtures we consider a multi-component variant of the kinetically constrained East model in which components have a different energy scale and can also diffuse when locally mobile regions, i.e., excitations, are present. Using a variational approach we determine an effective single component model with a single effective energy scale that best approximates a mixture. When scaled by this single effective energy, we show that experimental relaxation times of many liquid mixtures all collapse onto the 'parabolic law' predicted by dynamical facilitation theory. The model can be used to predict transport properties and glass transition temperatures of mixtures of glassy materials, with implications in atmospheric chemistry, biology, and pharmaceuticals.


title: Single Image Deraining: A Comprehensive Benchmark Analysis

abstract: We present a comprehensive study and evaluation of existing single image deraining algorithms, using a new large-scale benchmark consisting of both synthetic and real-world rainy images.This dataset highlights diverse data sources and image contents, and is divided into three subsets (rain streak, rain drop, rain and mist), each serving different training or evaluation purposes. We further provide a rich variety of criteria for dehazing algorithm evaluation, ranging from full-reference metrics, to no-reference metrics, to subjective evaluation and the novel task-driven evaluation. Experiments on the dataset shed light on the comparisons and limitations of state-of-the-art deraining algorithms, and suggest promising future directions.


title: Representations of real numbers induced by probability distributions on $\mathbb{N}$

abstract: We observe that a probability distribution supported by $\mathbb{N}$, induces a representation of real numbers in [0, 1) with digits in $\mathbb{N}$. We first study the Hausdorff dimension of sets with prescribed digits with respect to these representations. Than we determine the prevalent frequency of digits and the Hausdorff dimension of sets with prescribed frequencies of digits. As examples we consider the geometric distribution, the Poisson distribution and the zeta distribution.


title: Surprises in High-Dimensional Ridgeless Least Squares Interpolation

abstract: Interpolators -- estimators that achieve zero training error -- have attracted growing attention in machine learning, mainly because state-of-the art neural networks appear to be models of this type. In this paper, we study minimum $\ell_2$ norm ("ridgeless") interpolation in high-dimensional least squares regression. We consider two different models for the feature distribution: a linear model, where the feature vectors $x_i \in {\mathbb R}^p$ are obtained by applying a linear transform to a vector of i.i.d. entries, $x_i = \Sigma^{1/2} z_i$ (with $z_i \in {\mathbb R}^p$); and a nonlinear model, where the feature vectors are obtained by passing the input through a random one-layer neural network, $x_i = \varphi(W z_i)$ (with $z_i \in {\mathbb R}^d$, $W \in {\mathbb R}^{p \times d}$ a matrix of i.i.d. entries, and $\varphi$ an activation function acting componentwise on $W z_i$). We recover -- in a precise quantitative way -- several phenomena that have been observed in large-scale neural networks and kernel machines, including the "double descent" behavior of the prediction risk, and the potential benefits of overparametrization.


title: Sequential Optimization of Speed, Thermal Load, and Power Split in Connected HEVs

abstract: The emergence of connected and automated vehicles (CAVs) provides an unprecedented opportunity to capitalize on these technologies well beyond their original designed intents. While abundant evidence has been accumulated showing substantial fuel economy improvement benefits achieved through advanced powertrain control, the implications of the CAV operation on power and thermal management have not been fully investigated. In this paper, in order to explore the opportunities for the coordination between the onboard thermal management and the power split control, we present a sequential optimization solution for eco-driving speed trajectory planning, air conditioning (A/C) thermal load planning (eco-cooling), and powertrain control in hybrid electric CAVs to evaluate the individual as well as the collective energy savings through proactive usage of traffic data for vehicle speed prediction. Simulation results over a real-world driving cycle show that compared to a baseline non-CAV, 11.9%, 14.2%, and 18.8% energy savings can be accumulated sequentially through speed, thermal load, and power split optimizations, respectively.


title: Cosmological aspects of a hyperbolic solution in $f(R,T)$ gravity

abstract: This article deals with a cosmological scenario in $ f(R,T) $ gravity for a flat FLRW model of the universe. We consider the $ f(R,T) $ function as $ f(R)+f(T) $ which starts with a quadratic correction of the geometric term $ f(R) $ having structure $ f(R)=R+\alpha R^2 $, and a linear matter term $ f(T)=2\lambda T $. To achieve the solution of the gravitational field equations in the $ f(R,T) $ formalism, we take the form of a geometrical parameter, i.e. scale factor $ a(t)= sinh^{\frac{1}{n}}(\beta t) $ \cite{cha}, where $ \beta $ and $ n $ are model parameters. An eternal acceleration can be predicted by the model for $ 0<n<1 $, while the cosmic transition from the early decelerated phase to the present accelerated epoch can be anticipated for $ n\geq 1 $. The obtained model facilitate the formation of structure in the Universe according to the Jeans instability condition as our model transits from radiation dominated era to matter dominated era. We study the varying role of the equation of state parameter $ \omega $. We analyze our model by studying the behavior of the scalar field and discuss the energy conditions on our achieved solution. We examine the validity of our model via Jerk parameter, Om diagnostic, Velocity of sound and Statefinder diagnostic tools. We investigate the constraints on the model parameter $ n $ and $ H_0 $ (Hubble constant) using some observational datasets: $SNeIa$ dataset, $ H(z)$ (Hubble parameter) dataset, $ BAO $ (Baryon Acoustic Oscillation data) and their combinations as joint observational datasets $ H(z)$ + $ SNeIa $ and $ H(z)$ + $ SNeIa $ + $ BAO $. It is testified that the present study is well consistent with these observations. We also perform some cosmological tests and a detailed discussion of the model.


title: Bonnor-Melvin universe with a cosmological constant

abstract: We generalize the well-known Bonnor-Melvin solution of the Einstein-Maxwell equations to the case of a non-vanishing cosmological constant. The spacetime is again cylindrically symmetric and static but, unlike the original solution, it truly represents a homogeneous magnetic field.


title: Incomplete Nutation Diffusion Imaging: an ultrafast, single-scan approach for diffusion mapping

abstract: Purpose: Diffusion Magnetic Resonance Imaging (dMRI) is confounded by its long acquisition duration, thereby thwarting the detection of rapid microstructural changes, especially when diffusivity variations are accompanied by rapid changes in T2. The purpose of the present study is to accelerate dMRI to a single scan acquisition, and to enable a more accurate estimation of diffusivity as function of time. Methods: A general methodology termed Incomplete Initial Nutation Diffusion Imaging (INDI) capturing two diffusion contrasts in a single shot, is presented. INDI creates a longitudinal magnetization reservoir that facilitates the successive acquisition of two images separated by only a few milliseconds. INDI's theory is presented, followed by proof-of-concept ex- and in-vivo experiments at 16.4 T and 9.4 T. Results: Mean diffusivities (MDs) extracted from INDI were comparable with Diffusion Tensor Imaging (DTI) and the two-shot IDE in the water phantom. As expected in the brain tissues, DTI provided lower MD than UF-IDE and IDE, but IDE and UF-IDE were in excellent agreement. Simulations are presented for identifying the regimes where INDI is most beneficial. Conclusions: INDI accelerates dMRI acquisition to single-shot mode, which can be of great importance for mapping dynamic microstructural properties in-vivo without T2 bias.


title: BOLD-fMRI in the mouse auditory pathway

abstract: The auditory pathway is widely distributed throughout the brain, and is perhaps one of the most interesting networks in the context of neuroplasticity. Accurate mapping of neural activity in the entire pathway, preferably noninvasively, and with high resolution, could be instrumental for understanding such longitudinal processes. Functional magnetic resonance imaging (fMRI) has clear advantages for such characterizations, as it is noninvasive, provides relatively high spatial resolution and lends itself for repetitive studies, albeit relying on an indirect neurovascular coupling to deliver its information. Indeed, fMRI has been previously used to characterize the auditory pathway in humans and in rats. In the mouse, however, the auditory pathway has insofar only been mapped using manganese-enhanced MRI. Here, we describe a novel setup specifically designed for high-resolution mapping of the mouse auditory pathway using high-field fMRI. Robust and consistent Blood-Oxygenation-Level-Dependent (BOLD) responses were documented along nearly the entire auditory pathway, from the cochlear nucleus (CN), through the superior olivary complex (SOC), nuclei of the lateral lemniscus (LL), inferior colliculus (IC) and the medial geniculate body (MGB). By contrast, clear BOLD responses were not observed in auditory cortex (AC) in this study. Diverse BOLD latencies were mapped ROI- and pixel-wise using coherence analysis, evidencing different averaged BOLD time courses in different auditory centers. Some degree of tonotopy was identified in the IC, SOC, and MGB in the pooled dataset though it could not be assessed per subject due to a lack of statistical power. Given the importance of the mouse model in plasticity studies, animal models, and optogenetics, and fMRI's potential to map dynamic responses to specific cues, this first fMRI study of the mouse auditory pathway paves the way for future (...)


title: Joint Data compression and Computation offloading in Hierarchical Fog-Cloud Systems

abstract: Data compression has the potential to significantly improve the computation offloading performance in hierarchical fog-cloud systems. However, it remains unknown how to optimally determine the compression ratio jointly with the computation offloading decisions and the resource allocation. This joint optimization problem is studied in the current paper where we aim to minimize the maximum weighted energy and service delay cost (WEDC) of all users. First, we consider a scenario where data compression is performed only at the mobile users. We prove that the optimal offloading decisions have a threshold structure. Moreover, a novel three-step approach employing convexification techniques is developed to optimize the compression ratios and the resource allocation. Then, we address the more general design where data compression is performed at both the mobile users and the fog server. We propose three efficient algorithms to overcome the strong coupling between the offloading decisions and resource allocation. We show that the proposed optimal algorithm for data compression at only the mobile users can reduce the WEDC by a few hundred percent compared to computation offloading strategies that do not leverage data compression or use sub-optimal optimization approaches. Besides, the proposed algorithms for additional data compression at the fog server can further reduce the WEDC.


title: High-fidelity quantum tomography with imperfect measurements

abstract: In the current work we address the problem of quantum process tomography (QPT) in the case of imperfect preparation and measurement of the states which are used for QPT. The fuzzy measurements approach which helps us to efficiently take these imperfections into account is considered. However, to implement such a procedure one should have a detailed information about the errors. An approach for obtaining the partial information about them is proposed. It is based on the tomography of the ideal identity gate. This gate could be implemented by performing the measurement right after the initial state preparation. By using the result of the identity gate tomography we were able to significantly improve further QPT procedures. The proposed approach has been tested experimentally on the IBM superconducting quantum processor. As a result, we have obtained an increase in fidelity from 89% to 98% for Hadamard transformation and from 77% to 95% for CNOT gate.


title: Rapid Convergence of the Unadjusted Langevin Algorithm: Isoperimetry Suffices

abstract: We study the Unadjusted Langevin Algorithm (ULA) for sampling from a probability distribution $\nu = e^{-f}$ on $\mathbb{R}^n$. We prove a convergence guarantee in Kullback-Leibler (KL) divergence assuming $\nu$ satisfies a log-Sobolev inequality and the Hessian of $f$ is bounded. Notably, we do not assume convexity or bounds on higher derivatives. We also prove convergence guarantees in R\'enyi divergence of order $q > 1$ assuming the limit of ULA satisfies either the log-Sobolev or Poincar\'e inequality. We also prove a bound on the bias of the limiting distribution of ULA assuming third-order smoothness of $f$, without requiring isoperimetry.


title: The resolution of the universal Abel map via tropical geometry and applications

abstract: Let $g$ and $n$ be nonnegative integers and $\mathcal A=(a_0,\dots,a_n)$ a sequence of $n+1$ integers summing up to $d$. Let $\overline{\mathcal M}_{g,n+1}$ be the moduli space of $(n+1)$-pointed stable curves of genus $g$ and $\overline{\mathcal J}_{\mu,g}\rightarrow \overline{\mathcal M}_{g,1}$ be the Esteves' universal Jacobian, where $\mu$ is a universal genus-$g$ polarization of degree $d$. We give an explicit resolution of the universal Abel map $\alpha_{\mathcal A,\mu}\colon \overline{\mathcal M}_{g,n+1}\dashrightarrow \overline{\mathcal J}_{\mu,g}$, taking a pointed curve $(X,p_0,\dots,p_n)$ to $\mathcal{O}_X(\sum_{0\le i\le n} a_ip_i)$. The blowup of $\overline{\mathcal M}_{g,n+1}$ giving rise to the resolution is inspired by the resolution of the tropical analogue of the map $\alpha_{\mathcal A,\mu}$ (in the category of generalized cone complexes). As an application, we describe the double ramification cycle in terms of the universal sheaf inducing the resolution of the map $\alpha_{\mathcal A,\mu}$.


title: Accurate and Infinite Prime Prediction from Novel Quasi-Prime Analytical Methodology

abstract: It is known that prime numbers occupy specific geometrical patterns or moduli when numbers from one to infinity are distributed around polygons having sides that are integer multiple of number 6. In this paper, we will show that not only prime numbers occupy these moduli, but non-prime numbers sharing these same moduli have unique prime-ness properties. When utilizing digital root methodologies, these non-prime numbers provide a novel method to accurately identify prime numbers and prime factors without trial division or probabilistic-based methods. We will also show that the icositetragon (24-sided regular polygon) is a unique polygon pertaining to prime numbers and their ultimate incidence and distribution.


title: Identifying Maximal Non-Redundant Integer Cone Generators

abstract: A non-redundant integer cone generator (NICG) of dimension $d$ is a set $S$ of vectors from $\{0,1\}^d$ whose vector sum cannot be generated as a positive integer linear combination of a proper subset of $S$. The largest possible cardinality of NICG of a dimension $d$, denoted by $N(d)$, provides an upper bound on the sparsity of systems of integer equations with a large number of integer variables. A better estimate of $N(d)$ means that we can consider smaller sub-systems of integer equations when solving systems with many integer variables. Furthermore, given that we can reduce constraints on set algebra expressions to constraints on cardinalities of Venn regions, tighter upper bound on $N(d)$ yields a more efficient decision procedure for a logic of sets with cardinality constraints (BAPA), which has been applied in software verification. Previous attempts to compute $N(d)$ using SAT solvers have not succeeded even for $d=3$. The only known values were computed manually: $N(d)=d$ for $d < 4$ and $N(4) > 4$. We provide the first exact values for $d > 3$, namely, $N(4)=5$, $N(5)=7$, and $N(6)=9$, which is a significant improvement of the known asymptotic bound (which would give only e.g. $N(6) \le 29$, making a decision procedure impractical for $d=6$). We also give lower bounds for $N(7)$, $N(8)$, $N(9)$, and $N(10)$, which are: $11$, $13$, $14$, and $16$, respectively. We describe increasingly sophisticated specialized search algorithms that we used to explore the space of non-redundant generators and obtain these results.


title: Temperature-dependent Fowler-Nordheim electron barrier height in SiO2/4H-SiC MOS capacitors

abstract: This paper reports on the physical and temperature-dependent electrical characterizations of the oxide/semiconductor interface in MOS capacitors with a SiO2 layer deposited on 4H-SiC using dichlorosilane and nitrogen-based vapor precursors. The capacitors, subjected to a standard post deposition annealing process in N2O, exhibited an interface state density Dit = 9.0 x 1011cm-2eV-1 below the conduction band edge. At room temperature, a barrier height (conduction band offset) of 2.8 eV was observed, along with the presence of negative charges in the insulator. The SiO2 insulating properties were evaluated by studying the experimental temperature-dependence of the gate current. In particular, the temperature-dependent electrical measurements showed a negative temperature coefficient of the Fowler-Nordheim electron barrier height (dFB/dT = - 0.98 meV/{\deg}C), which was very close to the expected value for an ideal SiO2/4H-SiC system and much lower compared to the values reported for thermally grown SiO2. This smaller dependence of FB on the temperature and the increase of the current level with temperature in the transcharacteristics measured in the relative fabricated MOSFETs represents a clear advantage of our deposited SiO2 for the operation of MOSFET devices at high temperatures.


title: Box-constrained monotone $L_\infty$-approximations to Lipschitz regularizations, with applications to robust testing

abstract: Tests of fit to exact models in statistical analysis often lead to rejections even when the model is a useful approximate description of the random generator of the data. Among possible relaxations of a fixed model, the one defined by contamination neighbourhoods, namely, $\mathcal{V}_\alpha(P_0)=\{(1-\alpha)P_0+\alpha Q: Q \in \mathcal{P}\}$, where $\mathcal{P}$ is the set of all probabilities in the sample space, has received much attention, from its central role in Robust Statistics. For probabilities on the real line, consistent tests of fit to $\mathcal{V}_\alpha(P_0)$ can be based on $d_K(P_0,R_\alpha(P))$, the minimal Kolmogorov distance between $P_0$ and the set of trimmings of $P$, $R_\alpha(P)=\big\{\tilde P\in\mathcal{P}:\tilde P\ll P,\,{\textstyle \frac{d\tilde P}{dP}\leq\frac{1}{1-\alpha}}\, P\text{-a.s.}\big\}$. We show that this functional admits equivalent formulations in terms of, either best approximation in uniform norm by $L$-Lipschitz functions satisfying a box constraint, or as the best monotone approximation in uniform norm to the $L$-Lipschitz regularization, which is seen to be expressable in terms of the average of the Pasch-Hausdorff envelopes. This representation for the solution of the variational problem allows to obtain results showing stability of the functional $d_K(P_0,R_\alpha(P))$, as well as directional differentiability, providing the basis for a Central Limit Theorem for that functional.


title: Cooling rate effect on thermoremanent magnetization in archaeological baked clays: an experimental study on modern bricks

abstract: The influence of cooling rate on the intensity of thermoremanent magnetization (TRM) and the necessity to correct archaeo/palaeointensities for this effect have long been recognized. However the reliability of the correction is still questioned. We studied 35 bricks baked in two modern kilns (SK and BK) in known experimental conditions and with measurements of the direction and intensity of the geomagnetic field at the site. The smallest kiln (SK, 0.2 m 3) cooled in around 12 hours and the biggest kiln (BK, 8 m 3) in around 40 hours. Thermomagnetic, hysteresis and backfield curves indicated that the main magnetic carriers were Ti-poor titanomagnetites and Tipoortitanohematites. The fraction of the TRM carried by Ti-poor titanohematites is the maindifference between the two sets of bricks. This fraction is around 5-10% in bricks from BK kilnand up to 40% in those from SK kiln. Intensities of the Earth's magnetic field were determinedusing the original Thellier-Thellier protocol with correction of TRM anisotropy. The averageintensities overestimate the expected field intensity by 5% (SK) and 6% (BK). This resultemphasizes the necessity of the cooling rate correction. In order to have a detailed evaluation ofthe cooling rate effect, we used several slow cooling rates: 0.8, 0.4, 0.2 and 0.1{\textdegree}C/min. Thecorrection factors obtained with the 0.8{\textdegree}C/min cooling ranged between -2% and 21% and wereproportional to the TRM fraction carried by Ti-poor titanohematite. The higher proportion ofthese grains in bricks from SK kiln led to an overestimate of the correction factor and anunderestimate of the intensity by 7%. However, the expected intensity is recovered whentemperature steps higher than 580{\textdegree}C (i.e. in the range of Ti-poor titanohematites unblockingtemperatures) were excluded from the calculation of archaeointensity and cooling rate correction.In the case of the BK kiln bricks, for which Ti-poor titanohematites does not contributesignificantly to the TRM, all tested cooling rates give average intensities close to the expectedvalue. Incorrectly estimating the duration of the archaeological cooling has therefore a lowimpact on the accuracy of the archaeointensity data on these kinds of material.


title: Reachability in Vector Addition Systems is Primitive-Recursive in Fixed Dimension

abstract: The reachability problem in vector addition systems is a central question, not only for the static verification of these systems, but also for many inter-reducible decision problems occurring in various fields. The currently best known upper bound on this problem is not primitive-recursive, even when considering systems of fixed dimension. We provide significant refinements to the classical decomposition algorithm of Mayr, Kosaraju, and Lambert and to its termination proof, which yield an ACKERMANN upper bound in the general case, and primitive-recursive upper bounds in fixed dimension. While this does not match the currently best known TOWER lower bound for reachability, it is optimal for related problems.


title: Hydrodynamic limit for a diffusive system with boundary conditions

abstract: We study the hydrodynamic limit for the isothermal dynamics of an anharmonic chain under hyperbolic space-time scaling and with nonvanishing viscosity. The temperature is kept constant by a contact with a heat bath, realised via a stochastic momentum-preserving noise added to the dynamics. The noise is designed so it contributes to the macroscopic limit. Dirichlet boundary conditions are also considered: one end of the chain is kept fixed, while a time-varying tension is applied to the other end. Moreover, Neumann boundary conditions are added in such a way that the system produces the correct thermodynamic entropy in the macroscopic limit. We show that the volume stretch and momentum converge (in an appropriate sense) to a smooth solution of a system of parabolic conservation laws (isothermal Navier-Stokes equations in Lagrangian coordinates) with boundary conditions. Finally, changing the external tension allows us to define thermodynamic isothermal transformations between equilibrium states. We use this to deduce the first and the second law of Thermodynamics for our model.


title: Deep Learning for the Degraded Broadcast Channel

abstract: Machine learning has shown promising results for communications system problems. We present results on the use of deep auto-encoders in order to learn a transceiver for the multiuser degraded broadcast channel, and see that the auto encoder is able to learn to communicate on this channel using superposition coding. Additionally, the deep neural net is able to determine a bit labeling and optimize the per user power allocation that depends on the per user SNR.


title: The surface complex of Seifert fibered spaces

abstract: We define the surface complex for $3$-manifolds and embark on a case study in the arena of Seifert fibered spaces. The base orbifold of a Seifert fibered space captures some of the topology of the Seifert fibered space, so, not surprisingly, the surface complex of a Seifert fibered space always contains a subcomplex isomorphic to the curve complex of the base orbifold.


title: Electromagnetic Window into the Dawn of Black Holes

abstract: Massive 10^6-10^10 Msun black holes (BHs) are ubiquitous in local galactic nuclei. They were common by the time the Universe is several Gyr old, and many of them were in place within the first 1~Gyr after the Big Bang. Their quick assembly has been attributed to mechanisms such as the rapid collapse of gas into the nuclei of early protogalaxies, accretion and mergers of stellar-mass BHs accompanying structure formation at early times, and the runaway collapse of early, ultra-dense stellar clusters. The origin of the early massive BHs remains an intriguing and long-standing unsolved puzzle in astrophysics. Here we discuss strategies for discerning between BH seeding models using electromagnetic observations. We argue that the most direct answers will be obtained through detection of BHs with masses M<10^5 Msun at redshifts z>10, where we expect them to first form. Reaching out to these redshifts and down to these masses is crucial, because BHs are expected to lose the memory of their initial assembly by the time they grow well above 10^5 Msun and are incorporated into higher-mass galaxies. The best way to detect 10^4-10^5 Msun BHs at high redshifts is by a sensitive X-ray survey. Critical constraining power is augmented by establishing the properties and the environments of their host galaxies in deep optical/IR imaging surveys. Required OIR data can be obtained with the JWST and WFIRST missions. The required X-ray flux limits (down to 10^{-19} erg/s/cm^2) are accessible only with a next-generation X-ray observatory which has both high (sub-1") angular resolution and high throughput. A combination of deep X-ray and OIR surveys will be capable of probing several generic markers of the BH seed scenarios, and resolving the long-stanging puzzle of their origin. These electromagnetic observations are also highly synergistic with the information from LISA on high-z BH mergers.


title: The volumetric rate of normal type Ia supernovae in the local universe discovered by the Palomar Transient Factory

abstract: We present the volumetric rate of normal type Ia supernovae (SNe Ia) discovered by the Palomar Transient Factory (PTF). Using strict data-quality cuts, and considering only periods when the PTF maintained a regular cadence, PTF discovered 90 SNe Ia at $z\le0.09$ in a well-controlled sample over three years of operation (2010-2012). We use this to calculate the volumetric rate of SN Ia events by comparing this sample to simulations of hundreds of millions of SN Ia light curves produced in statistically representative realisations of the PTF survey. This quantifies the recovery efficiency of each PTF SN Ia event, and thus the relative weighting of each event. From this, the volumetric SN Ia rate was found to be $r_v=2.43\pm0.29\,\text{(stat)}_{-0.19}^{+0.33}\text{(sys)}\times10^{-5}\,\text{SNe yr}^{-1}\,\text{Mpc}^{-3}\, h_{70}^{3}$. This represents the most precise local measurement of the SN Ia rate. We fit a simple SN Ia delay-time distribution model, $\propto\mathrm{t}^{-\beta}$, to our PTF rate measurement combined with a literature sample of rate measurements from surveys at higher-redshifts. We find $\beta{\sim}1$, consistent with a progenitor channel governed by the gravitational in-spiral of binary white dwarfs.


title: Probing the superconducting ground state of ZrIrSi: A $\mu$SR study

abstract: The superconducting ground state of newly reported ZrIrSi is probed by means of $\mu$SR technique along with resistivity measurement. The occurrence of superconductivity at $T_\mathrm{C}$ = 1.7 K is confirmed by resistivity measurement. ZF-$\mu$SR study revealed that below $T_\mathrm{C}$, there is no spontaneous magnetic field in the superconducting state, indicates TRS is preserved in case of ZrIrSi. From TF-$\mu$SR measurement, we have estimated the superfluid density as a function of temperature, which is described by an isotropic $s-$wave model with a superconducting gap $2\Delta(0)/k_\mathrm{B}T_\mathrm{C}$ = 5.1, indicates the presence of strong spin-orbit coupling. {\it Ab-initio} electronic structure calculation indicates that there are four bands passing through the Fermi level, forming four Fermi surface pockets. We find that the low-energy bands are dominated by the $4d$-orbitals of transition metal Zr, with substantially lesser weight from the $5d$-orbitals of the Ir-atoms.


title: Notes on the analytic solution of box model equations for gravity-driven particle currents with constant volume

abstract: We summarize the physical equations and analytic solutions of three versions of the box model equations, suitable for the integral formulation of axisymmetric gravity-driven particle currents with constant volume. The first model is based on a simple constant resisting stress, while the second and third models assume flow dilution by particle deposition. The third model is characterized by assuming an interstitial fluid lighter than the ambient fluid. All the calculations are performed on a flat topography. Ambient fluid entrainment and cooling effects are not considered. All particles are assumed to deposit at the same velocity.


title: Data Augmentation for Leaf Segmentation and Counting Tasks in Rosette Plants

abstract: Deep learning techniques involving image processing and data analysis are constantly evolving. Many domains adapt these techniques for object segmentation, instantiation and classification. Recently, agricultural industries adopted those techniques in order to bring automation to farmers around the globe. One analysis procedure required for automatic visual inspection in this domain is leaf count and segmentation. Collecting labeled data from field crops and greenhouses is a complicated task due to the large variety of crops, growth seasons, climate changes, phenotype diversity, and more, especially when specific learning tasks require a large amount of labeled data for training. Data augmentation for training deep neural networks is well established, examples include data synthesis, using generative semi-synthetic models, and applying various kinds of transformations. In this paper we propose a method that preserves the geometric structure of the data objects, thus keeping the physical appearance of the data-set as close as possible to imaged plants in real agricultural scenes. The proposed method provides state of the art results when applied to the standard benchmark in the field, namely, the ongoing Leaf Segmentation Challenge hosted by Computer Vision Problems in Plant Phenotyping.


title: Bifurcation analysis and chaos control of periodically driven discrete fractional order memristive Duffing Oscillator

abstract: Discrete fractional order chaotic systems extends the memory capability to capture the discrete nature of physical systems. In this research, the memristive discrete fractional order chaotic system is introduced. The dynamics of the system was studied using bifurcation diagrams and phase space construction. The system was found chaotic with fractional order $0.465<n<0.562$. The dynamics of the system under different values makes it useful as a switch. Controllers were developed for the tracking control of the two systems to different trajectories. The effectiveness of the designed controllers were confirmed using simulations


title: Numerical Simulations of Gravitational Waves from Early-Universe Turbulence

abstract: We perform direct numerical simulations of magnetohydrodynamic turbulence in the early universe and numerically compute the resulting stochastic background of gravitational waves and relic magnetic fields. These simulations do not make the simplifying assumptions of earlier analytic work. If the turbulence is assumed to have an energy-carrying scale that is about a hundredth of the Hubble radius at the time of generation, as expected in a first-order phase transition, the peak of gravitational wave power will be in the mHz frequency range for a signal produced at the electroweak scale. The efficiency of gravitational wave (GW) production varies significantly with how the turbulence is driven. Detectability of turbulence at the electroweak scale by the planned Laser Interferometer Space Antenna (LISA) requires anywhere from 0.1% to 10% of the thermal plasma energy density to be in plasma motions or magnetic fields, depending on the model of the driving process. Our results predict a new universal form below the spectral peak frequency that is shallower than previously thought. This implies larger values of the GW energy spectra in the low-frequency range. This extends the range where turbulence is detectable with LISA to lower frequencies, corresponding to higher energy scales than the assumed energy-carrying scale.


title: A procedure for making high dynamic-range radio images: Deep imaging of the kiloparsec-scale radio structures of a distant blazar, NRAO 530, with JVLA data

abstract: Using JVLA data of Sgr A* and NRAO 530 (as a calibrator) at 5.5, 9 and 33 GHz, we developed a procedure for the reduction of wideband data. We have demonstrated that, correcting for residual interferometer errors such as residual delays, astronomers can now achieve high-fidelity radio images with a dynamic range exceeding 1,000,000:1. We outline the detailed procedure, noting that it can have broad application to the analysis of broadband continuum observations. We apply this procedure to observations of a distant blazar, NRAO 530, revealing its structures in unprecedented detail. Our 5.5-GHz image shows that the structure of NRAO 530 is prominently characterized by a moderately curved western jet terminating at a hot spot. Close to the radio core, an abrupt bending of the jet is revealed in the high-resolution (<100 mas) images at 33 GHz, showing an evolution of the position angle of the jet from the north at the VLBI scale (50 mas, or ~400 pc projected), increasing toward the west at the larger VLA scales (1 arcsec, or ~10 kpc). The continuation of the jet axis drift forms the curved western jet extending out to 20 arcsec. In contrast, a faint and broad counter-jet is present on the eastern side with a curvature antisymmetric to the western jet. The eastern jet terminates at a bright hotspot, forming an edge-brightened diffuse lobe. The observed contrast in brightness between the western and eastern jets suggests that the jets on the VLA scale are mildly relativistic. The radiation from the western jet is boosted while the radiation from the receding eastern jet is plausibly suppressed owing to the relativistic Doppler effect.


title: Reliability Maximization in Uncertain Graphs

abstract: Network reliability measures the probability that a target node is reachable from a source node in an uncertain graph, i.e., a graph where every edge is associated with a probability of existence. In this paper, we investigate the novel and fundamental problem of adding a small number of edges in the uncertain network for maximizing the reliability between a given pair of nodes. We study the NP-hardness and the approximation hardness of our problem, and design effective, scalable solutions. Furthermore, we consider extended versions of our problem (e.g., multiple source and target nodes can be provided as input) to support and demonstrate a wider family of queries and applications, including sensor network reliability maximization and social influence maximization. Experimental results validate the effectiveness and efficiency of the proposed algorithms.


title: A Polynomial-time Solution for Robust Registration with Extreme Outlier Rates

abstract: We propose a robust approach for the registration of two sets of 3D points in the presence of a large amount of outliers. Our first contribution is to reformulate the registration problem using a Truncated Least Squares (TLS) cost that makes the estimation insensitive to a large fraction of spurious point-to-point correspondences. The second contribution is a general framework to decouple rotation, translation, and scale estimation, which allows solving in cascade for the three transformations. Since each subproblem (scale, rotation, and translation estimation) is still non-convex and combinatorial in nature, out third contribution is to show that (i) TLS scale and (component-wise) translation estimation can be solved exactly and in polynomial time via an adaptive voting scheme, (ii) TLS rotation estimation can be relaxed to a semidefinite program and the relaxation is tight in practice, even in the presence of an extreme amount of outliers. We validate the proposed algorithm, named TEASER (Truncated least squares Estimation And SEmidefinite Relaxation), in standard registration benchmarks showing that the algorithm outperforms RANSAC and robust local optimization techniques, and favorably compares with Branch-and-Bound methods, while being a polynomial-time algorithm. TEASER can tolerate up to 99% outliers and returns highly-accurate solutions.


title: DC-SPP-YOLO: Dense Connection and Spatial Pyramid Pooling Based YOLO for Object Detection

abstract: Although the YOLOv2 method is extremely fast on object detection, its detection accuracy is restricted due to the low performance of its backbone network and the underutilization of multi-scale region features. Therefore, a dense connection (DC) and spatial pyramid pooling (SPP) based YOLO (DC-SPP-YOLO) method for ameliorating the object detection accuracy of YOLOv2 is proposed in this paper. Specifically, the dense connection of convolution layers is employed in the backbone network of YOLOv2 to strengthen the feature extraction and alleviate the vanishing-gradient problem. Moreover, an improved spatial pyramid pooling is introduced to pool and concatenate the multi-scale region features, so that the network can learn the object features more comprehensively. The DC-SPP-YOLO model is established and trained based on a new loss function composed of MSE (mean square error) loss and cross-entropy loss. The experimental results indicated that the mAP (mean Average Precision) of DC-SPP-YOLO is higher than that of YOLOv2 on the PASCAL VOC datasets and the UA-DETRAC datasets. The effectiveness of DC-SPP-YOLO method proposed is demonstrated.