title: Detection of the $5p-4f$ orbital crossing and its optical clock transition in Pr$^{9+}$

abstract: Recent theoretical works have proposed atomic clocks based on narrow optical transitions in highly charged ions. The most interesting candidates for searches of new physics are those which occur at rare orbital crossings where the shell structure of the periodic table is reordered. There are only three such crossings expected to be accessible in highly charged ions, and hitherto none have been observed as both experiment and theory have proven difficult. In this work we observe an orbital crossing in highly charged ions for the first time, in a system chosen to be tractable from both sides: Pr$^{9+}$. We present electron beam ion trap measurements of its spectra, including the inter-configuration lines that reveal the sought-after crossing. The proposed nHz-wide clock line, found to be at 452.334(1) nm, proceeds through hyperfine admixture of its upper state with an E2-decaying level. With state-of-the-art calculations we show that it has a very high sensitivity to new physics and extremely low sensitivity to external perturbations, making it a unique candidate for proposed precision studies.


title: Improved solution to data gathering with mobile mule

abstract: In this work we study the problem of collecting protected data in ad-hoc sensor network using a mobile entity called MULE. The objective is to increase information survivability in the network. Sensors from all over the network, route their sensing data through a data gathering tree, towards a particular node, called the $sink$. In case of a failed sensor, all the aggregated data from the sensor and from its children is lost. In order to retrieve the lost data, the MULE is required to travel among all the children of the failed sensor and to re-collect the data. There is a cost to travel between two points in the plane. We aim to minimize the MULE traveling cost, given that any sensor can fail. In order to reduce the traveling cost, it is necessary to find the optimal data gathering tree and the MULE location. We are considering the problem for the unit disk graphs (UDG) and Euclidean distance cost function. We propose a primal-dual algorithm that produces a $\left(20+\varepsilon\right)$-approximate solution for the problem, where $\varepsilon\rightarrow0$ as the sensor network spreads over a larger area. The algorithm requires $O\left(n^{3}\cdot\varDelta\left(G\right)\right)$ time to construct a gathering tree and to place the MULE, where $\varDelta\left(G\right)$ is the maximum degree in the graph and $n$ is the number of nodes.


title: A notion of rank for noncommutative quadratic forms on four generators

abstract: In this paper, we extend previous work, where a notion of rank, called mu-rank, was proposed for noncommutative quadratic forms on two and three generators. In particular, we provide a definition of mu-rank one and two for noncommutative quadratic forms on four generators.


title: Asymmetric completions of partial metric spaces

abstract: Ge and Lin (2015) proved the existence and the uniqueness of p-Cauchy completions of partial metric spaces under symmetric denseness. They asked if every (non-empty) partial metric space $X$ has a p-Cauchy completion $\bar{X}$ such that $X$ is dense but not symmetrically dense in $\bar{X}$. We construct asymmetric p-Cauchy completions for all non-empty partial metric spaces. This gives a positive answer to the question. We also provide a nonstandard construction of partial metric completions.


title: Ordering-Based Causal Structure Learning in the Presence of Latent Variables

abstract: We consider the task of learning a causal graph in the presence of latent confounders given i.i.d.~samples from the model. While current algorithms for causal structure discovery in the presence of latent confounders are constraint-based, we here propose a score-based approach. We prove that under assumptions weaker than faithfulness, any sparsest independence map (IMAP) of the distribution belongs to the Markov equivalence class of the true model. This motivates the \emph{Sparsest Poset} formulation - that posets can be mapped to minimal IMAPs of the true model such that the sparsest of these IMAPs is Markov equivalent to the true model. Motivated by this result, we propose a greedy algorithm over the space of posets for causal structure discovery in the presence of latent confounders and compare its performance to the current state-of-the-art algorithms FCI and FCI+ on synthetic data.


title: First-principles approach to electric polarization and dielectric constant calculations using generalized Wannier functions

abstract: We describe a method to calculate the electronic properties of an insulator under an applied electric field. It is based on the minimization of an electric enthalpy functional with respect to the orbitals, which behave as Wannier functions under crystal translations, but are not necessarily orthogonal. This paper extends the approach of Nunes and Vanderbilt (NV) [Phys. Rev. Lett. 73, 712 (1994)], who demonstrated that a Wannier function representation can be used to study insulating crystals in the presence of a finite electric field. According to a study by Fern\'{a}ndez et al. [Phys. Rev. B. 58, R7480 (1998)], first-principles implementations of the NV approach suffer from the impact of the localization constraint on the orthogonal wave functions, what affects the accuracy of the physical results. We show that because non-orthogonal generalized Wannier functions can be more localized than their orthogonal counterparts, the error due to localization constraints is reduced, thus improving the accuracy of the calculated physical quantities.


title: Generalizing the notion of rank to noncommutative quadratic forms

abstract: In 2010, Cassidy and Vancliff extended the notion of a quadratic form on n generators to the noncommutative setting. In this article, we suggest a notion of rank for such noncommutative quadratic forms, where n = 2 or 3. Since writing an arbitrary quadratic form as a sum of squares fails in this context, our methods entail rewriting an arbitrary quadratic form as a sum of products. In so doing, we find analogs for 2 x 2 minors and determinant of a 3 x 3 matrix in this noncommutative setting.


title: Demystifying Graph Databases: Analysis and Taxonomy of Data Organization, System Designs, and Graph Queries

abstract: Graph processing has become an important part of multiple areas of computer science, such as machine learning, computational sciences, medical applications, social network analysis, and many others. Numerous graphs such as web or social networks may contain up to trillions of edges. Often, these graphs are also dynamic (their structure changes over time) and have domain-specific rich data associated with vertices and edges. Graph database systems such as Neo4j enable storing, processing, and analyzing such large, evolving, and rich datasets. Due to the sheer size of such datasets, combined with the irregular nature of graph processing, these systems face unique design challenges. To facilitate the understanding of this emerging domain, we present the first survey and taxonomy of graph database systems. We focus on identifying and analyzing fundamental categories of these systems (e.g., triple stores, tuple stores, native graph database systems, or object-oriented systems), the associated graph models (e.g., RDF or Labeled Property Graph), data organization techniques (e.g., storing graph data in indexing structures or dividing data into records), and different aspects of data distribution and query execution (e.g., support for sharding and ACID). 51 graph database systems are presented and compared, including Neo4j, OrientDB, or Virtuoso. We outline graph database queries and relationships with associated domains (NoSQL stores, graph streaming, and dynamic graph algorithms). Finally, we describe research and engineering challenges to outline the future of graph databases.


title: Point modules over regular graded skew Clifford algebras

abstract: Results of Vancliff, Van Rompay and Willaert in 1998 prove that point modules over a regular graded Clifford algebra (GCA) are determined by (commutative) quadrics of rank at most two that belong to the quadric system associated to the GCA. In 2010, Cassidy and Vancliff generalized the notion of a GCA to that of a graded skew Clifford algebra (GSCA). The results in this article show that prior results may be extended, with suitable modification, to GSCAs. In particular, using the notion of {\mu}-rank introduced recently by the authors, the point modules over a regular GSCA are determined by (noncommutative) quadrics of {\mu}-rank at most two that belong to the noncommutative quadric system associated to the GSCA.


title: Generalized Uncertainty Principle in three-dimensional gravity and the BTZ black hole

abstract: We investigate the structure of the gravity-induced Generalized Uncertainty Principle in three dimensions. The subtleties of lower dimensional gravity, and its important differences with respect to four and higher dimensions, are duly taken into account, by considering different possible candidates for the gravitational radius, $R_g$, that is the minimal length/maximal resolution of the quantum mechanical localization process. We find that the event horizon of the $M \neq 0$ Ba\~{n}ados-Teitelboim-Zanelli micro black hole furnishes the most consistent $R_g$. This allows us to obtain a suitable formula for the Generalized Uncertainty Principle in three dimensions, and also to estimate the corrections induced by the latter on the Hawking temperature and Bekenstein entropy. We also point to the extremal $M=0$ case, and its natural unit of length introduced by the cosmological constant, $\ell = 1 / \sqrt{-\Lambda}$, as a possible alternative to $R_g$, and present a condensed matter analog realization of this scenario.


title: SneakySnake: A Fast and Accurate Universal Genome Pre-Alignment Filter for CPUs, GPUs, and FPGAs

abstract: Motivation: We introduce SneakySnake, a highly parallel and highly accurate pre-alignment filter that remarkably reduces the need for computationally costly sequence alignment. The key idea of SneakySnake is to reduce the approximate string matching (ASM) problem to the single net routing (SNR) problem in VLSI chip layout. In the SNR problem, we are interested in finding the optimal path that connects two terminals with the least routing cost on a special grid layout that contains obstacles. The SneakySnake algorithm quickly solves the SNR problem and uses the found optimal path to decide whether or not performing sequence alignment is necessary. Reducing the ASM problem into SNR also makes SneakySnake efficient to implement on CPUs, GPUs, and FPGAs. Results: SneakySnake significantly improves the accuracy of pre-alignment filtering by up to four orders of magnitude compared to the state-of-the-art pre-alignment filters, Shouji, GateKeeper, and SHD. For short sequences, SneakySnake accelerates Edlib (state-of-the-art implementation of Myers's bit-vector algorithm) and Parasail (state-of-the-art sequence aligner with a configurable scoring function), by up to 37.7x and 43.9x (>12x on average), respectively, with its CPU implementation, and by up to 413x and 689x (>400x on average), respectively, with FPGA and GPU acceleration. For long sequences, the CPU implementation of SneakySnake accelerates Parasail and KSW2 (sequence aligner of minimap2) by up to 979x (276.9x on average) and 91.7x (31.7x on average), respectively. As SneakySnake does not replace sequence alignment, users can still obtain all capabilities (e.g., configurable scoring functions) of the aligner of their choice, unlike existing acceleration efforts that sacrifice some aligner capabilities. Availability: https://github.com/CMU-SAFARI/SneakySnake


title: Musical Instrument Playing Technique Detection Based on FCN: Using Chinese Bowed-Stringed Instrument as an Example

abstract: Unlike melody extraction and other aspects of music transcription, research on playing technique detection is still in its early stages. Compared to existing work mostly focused on playing technique detection for individual single notes, we propose a general end-to-end method based on Sound Event Detection by FCN for musical instrument playing technique detection. In our case, we choose Erhu, a well-known Chinese bowed-stringed instrument, to experiment with our method. Because of the limitation of FCN, we present an algorithm to detect on variable length audio. The effectiveness of the proposed framework is tested on a new dataset, its categorization of techniques is similar to our training dataset. The highest accuracy of our 3 experiments on the new test set is 87.31%. Furthermore, we also evaluate the performance of the proposed framework on 10 real-world studio music (produced by midi) and 7 real-world recording samples to address the ability of generalization on our model.


title: Diverse Behavior Is What Game AI Needs: Generating Varied Human-Like Playing Styles Using Evolutionary Multi-Objective Deep Reinforcement Learning

abstract: this paper has been withdrawn


title: Ultrafast energy absorption and photoexcitation of bulk plasmon in crystalline silicon subjected to intense near-infrared ultrashort laser pulses

abstract: We investigate the non-linear response and energy absorption in bulk silicon irradiated by intense 12-fs near-infrared laser pulses. Depending on the laser intensity, we distinguish two regimes of non-linear absorption of the laser energy: for low intensities, energy deposition and photoionization involve perturbative three-photon transition through the direct bandgap of silicon. For laser intensities near and above 10$^{14}$ W/cm$^2$, corresponding to photocarrier density of order 10$^{22}$ cm$^{-3}$, we find that absorption at near-infrared wavelengths is greatly enhanced due to excitation of bulk plasmon resonance. In this regime, the energy transfer to electrons exceeds a few times the thermal melting threshold of Si. The optical reflectivity of the photoexcited solid is found in good qualitative agreement with existing experimental data. In particular, the model predicts that the main features of the reflectivity curve of photoexcited Si as a function of the laser fluence are determined by the competition between state and band filling associated with Pauli exclusion principle and Drude free-carrier response. The non-linear response of the photoexcited solid is also investigated for irradiation of silicon with a sequence of two strong and temporary non-overlapping pulses. The cumulative effect of the two pulses is non-additive in terms of deposited energy. Photoionization and energy absorption on the leading edge of the second pulse is greatly enhanced due to free carrier absorption.


title: Boosting Network Weight Separability via Feed-Backward Reconstruction

abstract: This paper proposes a new evaluation metric and boosting method for weight separability in neural network design. In contrast to general visual recognition methods designed to encourage both intra-class compactness and inter-class separability of latent features, we focus on estimating linear independence of column vectors in weight matrix and improving the separability of weight vectors. To this end, we propose an evaluation metric for weight separability based on semi-orthogonality of a matrix and Frobenius distance, and the feed-backward reconstruction loss which explicitly encourages weight separability between the column vectors in the weight matrix. The experimental results on image classification and face recognition demonstrate that the weight separability boosting via minimization of feed-backward reconstruction loss can improve the visual recognition performance, hence universally boosting the performance on various visual recognition tasks.


title: A Conformal Three-Field Formulation for Nonlinear Elasticity: From Differential Complexes to Mixed Finite Element Methods

abstract: We introduce a new class of mixed finite element methods for 2D and 3D compressible nonlinear elasticity. The independent unknowns of these conformal methods are displacement, displacement gradient, and the first Piola-Kirchhoff stress tensor. The so-called edge finite elements of the curl operator is employed to discretize the trial space of displacement gradients. Motivated by the differential complex of nonlinear elasticity, this choice guarantees that discrete displacement gradients satisfy the Hadamard jump condition for the strain compatibility. We study the stability of the proposed mixed finite element methods by deriving some inf-sup conditions. By considering 32 choices of simplicial conformal finite elements of degrees 1 and 2, we show that 10 choices are not stable as they do not satisfy the inf-sup conditions. We numerically study the stable choices and conclude that they can achieve optimal convergence rates. By solving several 2D and 3D numerical examples, we show that the proposed methods are capable of providing accurate approximations of strain and stress.


title: The Sakellaridis-Venkatesh Conjecture for $ U(2)\backslash SO_{2,3}$

abstract: This article gives the Plancherel decomposition of $L^2\left(U(2)(F)\backslash SO_{2,3}(F)\right)$, where $F$ is a local field with characteristic $0$, and the relative character identities. Finally, we obtain a factorization of the global period $P_{U(2)}$.


title: Secure Digital Administration in Medical Environment

abstract: The efficiency and service quality in a medical environment can be improved by using electronic documents (or e-docs) and digital signatures to speed up both doctors' activity and to provide in the same time easy retrieval and use of needed data without loosing convenience. Our proposed solution satisfies these needs by making use of the AIDA system and many cutting-edge techniques to build a digitalized management system. To be more specific, we present firstly the AIDA document exchange framework for the management (creation, search, storage) of e-docs expressed in XML format. The framework provides security services, like privacy, authorization, authentication and non-repudiation. We describe next the use of AIDA in a real world medical service, namely the creation of electronic medical records (e- MRs). Doctors can use mobile devices that embed a secure trustworthy environment defined also in AIDA for digitally signing the e-MRs. Technically speaking a handheld PC equipped with a smart-card reader and integrating What You See Is What You Sign (WYSIWYS) features will be used for viewing and signing the e-MRs. Furthermore the proposed system is easily integrated with the infrastructure (e.g., database system) already in use at hospital's administration site and allows easy handling and updating of data processed on the mobile devices. The use of web interface for the operations to be executed on the mobile device or for those executed on the remote part of the system makes the whole application homogeneous and easy to use.


title: Inertial focusing in two dimensional flows with sharp viscosity stratification in a microchannel

abstract: Recent experimental studies have shown that particle transfer across streamlines can be controlled passively using stratified flows of co-flowing streams at a finite Reynolds number. The stratification modifies the forces acting on particles through the curvature of the undisturbed velocity profile. In this study, we numerically analyze the particle migration in stratified flow of two liquids of different viscosities flowing parallel to each other between two infinite parallel plates. Particle migration in two different flow fields is considered: (i) Couette flow and (ii) Poiseuille flow. A numerical approach using an immersed boundary method is employed to perform two dimensional simulations and determine conditions when particle migration from one fluid to the other can occur. This has implications in separating particles from a fluid without a membrane. The effect of the viscosity ratio, flowrate ratio, Reynolds number and particle size on focusing position are analyzed to identify conditions under which the particle migrates from one fluid to the other. It is shown that the particle migrates to the fluid with a lower viscosity in case of stratified Couette flow when the holdup of the low viscous fluid is sufficiently high. In Poiseuille flow, particle migrates to the less viscous fluid beyond a critical flowrate ratio for a fixed viscosity ratio. This critical flowrate ratio increases with particle size.


title: High Energy Hadron Production, Self-Organized Criticality and Absorbing State Phase Transition

abstract: In high energy nuclear collisions, production rates of light nuclei agree with the predictions of an ideal gas at a temperature $T=155 \pm 10$ MeV. In an equilibrium hadronic medium of this temperature, light nuclei cannot survive. In this contribution, we suggest that the observed behavior is due to an evolution in global non-equilibrium, leading to self-organized criticality and to hadron formation as an absorbing state phase transition for color degrees of freedom. At the confinement point, the initial quark-gluon medium becomes quenched by the vacuum, breaking up into all allowed free hadronic and nuclear mass states, without (or with a very short-live) subsequent formation of thermal hadronic medium.


title: Supporting Multi-point Fan Design with Dimension Reduction

abstract: Motivated by the idea of turbomachinery active subspace performance maps, this paper studies dimension reduction in turbomachinery 3D CFD simulations. First, we show that these subspaces exist across different blades---under the same parametrization---largely independent of their Mach number or Reynolds number. This is demonstrated via a numerical study on three different blades. Then, in an attempt to reduce the computational cost of identifying a suitable dimension reducing subspace, we examine statistical sufficient dimension reduction methods, including sliced inverse regression, sliced average variance estimation, principal Hessian directions and contour regression. Unsatisfied by these results, we evaluate a new idea based on polynomial variable projection---a non-linear least squares problem. Our results using polynomial variable projection clearly demonstrate that one can accurately identify dimension reducing subspaces for turbomachinery functionals at a fraction of the cost associated with prior methods. We apply these subspaces to the problem of comparing design configurations across different flight points on a working line of a fan blade. We demonstrate how designs that offer a healthy compromise between performance at cruise and sea-level conditions can be easily found by visually inspecting their subspaces.


title: CAI4CAI: The Rise of Contextual Artificial Intelligence in Computer Assisted Interventions

abstract: Data-driven computational approaches have evolved to enable extraction of information from medical images with a reliability, accuracy and speed which is already transforming their interpretation and exploitation in clinical practice. While similar benefits are longed for in the field of interventional imaging, this ambition is challenged by a much higher heterogeneity. Clinical workflows within interventional suites and operating theatres are extremely complex and typically rely on poorly integrated intra-operative devices, sensors, and support infrastructures. Taking stock of some of the most exciting developments in machine learning and artificial intelligence for computer assisted interventions, we highlight the crucial need to take context and human factors into account in order to address these challenges. Contextual artificial intelligence for computer assisted intervention, or CAI4CAI, arises as an emerging opportunity feeding into the broader field of surgical data science. Central challenges being addressed in CAI4CAI include how to integrate the ensemble of prior knowledge and instantaneous sensory information from experts, sensors and actuators; how to create and communicate a faithful and actionable shared representation of the surgery among a mixed human-AI actor team; how to design interventional systems and associated cognitive shared control schemes for online uncertainty-aware collaborative decision making ultimately producing more precise and reliable interventions.


title: Molecular Dynamics Simulation of Imidazolium CnMIM-BF4 Ionic Liquids using a Coarse Grained Force-Field

abstract: Ionic Liquids feature thermophysical properties that are of interest in solvents, energy storage materials and tenable lubrication applications. Recently, a series of coarse grained (CG) models was developed to investigate 1-ethyl-3-methylimidazolium tetrafluoroborate [C2MIM][BF4] and 1-butyl-3-methylimidazolium hexafluorophosphate [C4MIM][PF6] ionic liquids. Building on these CG models we derive force-fields to investigate the [C_{2-8}MIM] [BF4] family of ionic liquids, as a starting step for systematic investigations of lubrication under nanoconfinement conditions. The simualted equation of state and diffusion coefficients are in good agreement with experimental data and with all-atom force fields. We use these model to analyze the nano nanostructuring of ILs characteristic of cations with longer aliphatic chains as well as the ILs liquidvapour interfacial structure. The CG nature of these models enables the simulation of very long time scales, which are needed to computer reliable results of dynamic and interfacial properties. For [C>4MIM] [BF4] the break in symmetry associated to the liquid-vapor interface induces nanostructuring in polar and non-polar domains in the direction perpendicular to the interface plane, hence mimicking the behavior observed in the bulk phases.


title: A note on Lagrangian submanifolds of twistor spaces and their relation to superminimal surfaces

abstract: In this paper a bijective correspondence between superminimal surfaces of an oriented Riemannian $4$-manifold and particular Lagrangian submanifolds of the twistor space over the $4$-manifold is proven. More explicitly, for every superminimal surface a submanifold of the twistor space is constructed which is Lagrangian for all the natural almost Hermitian structures on the twistor space. The twistor fibration restricted to the constructed Lagrangian gives a circle bundle over the superminimal surface. Conversely, if a submanifold of the twistor space is Lagrangian for all the natural almost Hermitian structures, then the Lagrangian projects to a superminimal surface and is is contained in the Lagrangian constructed from this surface. In particular this produces many Lagrangian submanifolds of the twistor spaces $\mathbb{C} P^3$ and $\mathbb{F}_{1,2}(\mathbb{C}^3)$ with respect to both the K\"{a}hler structure as well as the nearly K\"{a}hler structure. Moreover, it is shown that these Lagrangian submanifolds are minimal submanifolds.


title: Fast Security-Constrained Optimal Power Flow through Low-Impact and Redundancy Screening

abstract: Determining contingency aware dispatch decisions by solving a security-constrained optimal power flow (SCOPF) is challenging for real-world power systems, as the high problem dimensionality often leads to impractical computational requirements. This problem becomes more severe when the SCOPF has to be solved not only for a single instance, but for multiple periods, e.g. in the context of electricity market analyses. This paper proposes an algorithm that identifies the minimal set of constraints that exactly define the space of feasible nodal injections for a given network and contingency scenarios. By internalizing the technical limits of the nodal injections and enforcing a minimal worst-case impact of contingencies to line flows, computational effort can be further improved. The case study applies and analyzes the methods on the IEEE 118 and A&M 2000 bus systems, as well as the German and European transmission systems. In all tested cases the proposed algorithm identifies at least 95% of the network and security constraints as redundant, leading to significant SCOPF solve time reductions. Scalability and practical implementation are explicitly discussed. The code and input data of the case study is published supplementary to the paper under an open-source license.


title: Bounds on optimal transport maps onto log-concave measures

abstract: We consider strictly log-concave measures, whose bounds degenerate at infinity. We prove that the optimal transport map from the Gaussian onto such a measure is locally Lipschitz, and that the eigenvalues of its Jacobian have controlled growth at infinity.


title: Differentiable Deep Clustering with Cluster Size Constraints

abstract: Clustering is a fundamental unsupervised learning approach. Many clustering algorithms -- such as $k$-means -- rely on the euclidean distance as a similarity measure, which is often not the most relevant metric for high dimensional data such as images. Learning a lower-dimensional embedding that can better reflect the geometry of the dataset is therefore instrumental for performance. We propose a new approach for this task where the embedding is performed by a differentiable model such as a deep neural network. By rewriting the $k$-means clustering algorithm as an optimal transport task, and adding an entropic regularization, we derive a fully differentiable loss function that can be minimized with respect to both the embedding parameters and the cluster parameters via stochastic gradient descent. We show that this new formulation generalizes a recently proposed state-of-the-art method based on soft-$k$-means by adding constraints on the cluster sizes. Empirical evaluations on image classification benchmarks suggest that compared to state-of-the-art methods, our optimal transport-based approach provide better unsupervised accuracy and does not require a pre-training phase.


title: YBC, a stellar bolometric corrections database with variable extinction coefficients: an application to PARSEC isochrones

abstract: We present the \texttt{YBC} database of stellar bolometric corrections (BCs), available at \url{http://stev.oapd.inaf.it/YBC}. We homogenize widely-used theoretical stellar spectral libraries and provide BCs for many popular photometric systems, including the Gaia filters. The database can be easily extended to additional photometric systems and stellar spectral libraries. The web interface allows users to transform their catalogue of theoretical stellar parameters into magnitudes and colours of selected filter sets. The BC tables can be downloaded or also be implemented into large simulation projects using the interpolation code provided with the database. We compute extinction coefficients on a star-by-star basis, hence taking into account the effects of spectral type and non-linearity dependency on the total extinction. We illustrate the use of these BCs in \texttt{PARSEC} isochrones. We show that using spectral-type dependent extinction coefficients is quite necessary for Gaia filters whenever $\Av\gtrsim 0.5$\,mag. BC tables for rotating stars and tables of limb-darkening coefficients are also provided.


title: Testing the Spectroscopic Extraction of Suppression of Convective Blueshift

abstract: Efforts to detect low-mass exoplanets using stellar radial velocities (RVs) are currently limited by magnetic photospheric activity. Suppression of convective blueshift is the dominant magnetic contribution to RV variability in low-activity Sun-like stars. Due to convective plasma motions, the magnitude of RV contributions from the suppression of convective blueshift is related to the depth of formation of photospheric spectral lines of a given species used to compute the RV time series. Meunier et al. (2017), used this relation to demonstrate a method for spectroscopic extraction of the suppression of convective blueshift in order to isolate RV contributions, including planetary RVs, that contribute equally to the timeseries for each spectral line. Here, we extract disk-integrated solar RVs from observations over a 2.5 year time span made with the solar telescope integrated with the HARPS-N spectrograph at the Telescopio Nazionale Galileo (La Palma, Canary Islands, Spain). We apply the methods outlined by Meunier et al. (2017). We are not, however, able to isolate physically meaningful contributions of the suppression of convective blueshift from this solar dataset, potentially because our dataset is from solar minimum when the suppression of convective blueshift may not sufficiently dominate activity contributions to RVs. This result indicates that, for low-activity Sun-like stars, one must include additional RV contributions from activity sources not considered in the Meunier et al. (2017) model at different timescales as well as instrumental variation in order to reach the sub-meter per second RV sensitivity necessary to detect low-mass planets in orbit around Sun-like stars.


title: Ab initio study of lattice dynamics of dodecaborides ZrB12 and LuB12

abstract: We performed ab initio lattice-dynamics calculations of frame-cluster dodecaborides ZrB12 and LuB12. As a whole, our calculated phonon frequencies and atom-projected density of states are consistent with the results of available first-principles calculations and experimental measurements. So we conclude that the ab initio DFT approach is quite appropriate to study the sufficiently subtle physics of these compounds. Our experiment-independent calculations provide an explicit quantitative confirmation of mixing the eigenvectors of boron and metal vibrations, which was previously observed in experiments.


title: Landing Probabilities of Random Walks for Seed-Set Expansion in Hypergraphs

abstract: We describe the first known mean-field study of landing probabilities for random walks on hypergraphs. In particular, we examine clique-expansion and tensor methods and evaluate their mean-field characteristics over a class of random hypergraph models for the purpose of seed-set community expansion. We describe parameter regimes in which the two methods outperform each other and propose a hybrid expansion method that uses partial clique-expansion to reduce the projection distortion and low-complexity tensor methods applied directly on the partially expanded hypergraphs.


title: You Can Run, But You Cannot Hide: Using Elevation Profiles to Breach Location Privacy through Trajectory Prediction

abstract: The extensive use of smartphones and wearable devices has facilitated many useful applications. For example, with Global Positioning System (GPS)-equipped smart and wearable devices, many applications can gather, process, and share rich metadata, such as geolocation, trajectories, elevation, and time. For example, fitness applications, such as Strava and Runkeeper, utilize information for activity tracking, and have recently witnessed a boom in popularity. Those trackers have their own web platforms, and allow users to share activities on such platforms, or even with other social network platforms. To preserve privacy of users while allowing sharing, those platforms allow users to disclose partial information, such as the elevation profile for an activity, which supposedly will not leak the location trajectory. In this work we examine the extent to which publicly available elevation profiles can be used to predict the location trajectory of users. To tackle this problem, we devise three threat settings under which the city, borough, or even a route can be predicted. Those threat settings define the amount of information available to the adversary to launch the prediction attacks. Establishing that simple features of elevation profiles, e.g., spectral features, are insufficient, we devise both natural language processing (NLP)-inspired text-like representation and computer vision-inspired image-like representation of elevation profiles, and we convert the problem at hand into text and image classification problem. We use both traditional machine learning- and deep learning-based techniques, and achieve a prediction success rate ranging from 59.59% to 95.83%. The findings are alarming, and highlight that sharing information such as elevation profile may have significant privacy risks.


title: Periodic orbits of the retrograde coorbital problem

abstract: Asteroid (514107) Kaepaokaawela is the first example of an object in the 1/1 mean motion resonance with Jupiter with retrograde motion around the Sun. Its orbit was shown to be stable over the age of the Solar System which implies that it must have been captured from another star when the Sun was still in its birth cluster. Kaepaokaawela orbit is also located at the peak of the capture probability in the coorbital resonance. Identifying the periodic orbits that Kaepaokaawela and similar asteroids followed during their evolution is an important step towards precisely understanding their capture mechanism. Here, we find the families of periodic orbits in the two-dimensional retrograde coorbital problem and analyze their stability and bifurcations into three-dimensional periodic orbits. Our results explain the radical differences observed in 2D and 3D coorbital capture simulations. In particular, we find that analytical and numerical results obtained for planar motion are not always valid at infinitesimal deviations from the plane.


title: Learning from both experts and data

abstract: In this work we study the problem of inferring a discrete probability distribution using both expert knowledge and empirical data. This is an important issue for many applications where the scarcity of data prevents a purely empirical approach. In this context, it is common to rely first on an initial domain knowledge a priori before proceeding to an online data acquisition. We are particularly interested in the intermediate regime where we do not have enough data to do without the initial expert a priori of the experts, but enough to correct it if necessary. We present here a novel way to tackle this issue with a method providing an objective way to choose the weight to be given to experts compared to data. We show, both empirically and theoretically, that our proposed estimator is always more efficient than the best of the two models (expert or data) within a constant.


title: Maximum sparse induced subgraphs of the binomial random graph with given number of edges

abstract: We prove that a.a.s. the maximum size of an induced subtree of the binomial random graph $G(n,p)$ is concentrated in 2 consecutive points. We also prove that, given a non-negative integer-valued function $t(k)<\varepsilon k^2$, under a certain smoothness condition on this function, a.a.s. the maximum size $k$ of an induced subgraph with exactly $t(k)$ edges of $G(n,p)$ is concentrated in 2 consecutive points as well.


title: Overcoming the curse of dimensionality for some Hamilton--Jacobi partial differential equations via neural network architectures

abstract: We propose new and original mathematical connections between Hamilton-Jacobi (HJ) partial differential equations (PDEs) with initial data and neural network architectures. Specifically, we prove that some classes of neural networks correspond to representation formulas of HJ PDE solutions whose Hamiltonians and initial data are obtained from the parameters of the neural networks. These results do not rely on universal approximation properties of neural networks; rather, our results show that some classes of neural network architectures naturally encode the physics contained in some HJ PDEs. Our results naturally yield efficient neural network-based methods for evaluating solutions of some HJ PDEs in high dimension without using grids or numerical approximations. We also present some numerical results for solving some inverse problems involving HJ PDEs using our proposed architectures.


title: Detecting muscle activation using ultrasound speed of sound inversion with deep learning

abstract: Functional muscle imaging is essential for diagnostics of a multitude of musculoskeletal afflictions such as degenerative muscle diseases, muscle injuries, muscle atrophy, and neurological related issues such as spasticity. However, there is currently no solution, imaging or otherwise, capable of providing a map of active muscles over a large field of view in dynamic scenarios. In this work, we look at the feasibility of longitudinal sound speed measurements to the task of dynamic muscle imaging of contraction or activation. We perform the assessment using a deep learning network applied to pre-beamformed ultrasound channel data for sound speed inversion. Preliminary results show that dynamic muscle contraction can be detected in the calf and that this contraction can be positively assigned to the operating muscles. Potential frame rates in the hundreds to thousands of frames per second are necessary to accomplish this.


title: Optical detection of Mott and generalized Wigner crystal states in WSe2/WS2 moir\'e superlattices

abstract: Moir\'e superlattices are emerging as a new route for engineering strongly correlated electronic states in two-dimensional van der Waals heterostructures, as recently demonstrated in the correlated insulating and superconducting states in magic-angle twisted bilayer graphene and ABC trilayer graphene/boron nitride moir\'e superlattices. Transition metal dichalcogenide (TMDC) moir\'e heterostructures provide another exciting model system to explore correlated quantum phenomena, with the addition of strong light-matter interactions and large spin-orbital coupling. Here we report the optical detection of strongly correlated phases in semiconducting WSe2/WS2 moir\'e superlattices. Our sensitive optical detection technique reveals a Mott insulator state at one hole per superlattice site ({\nu} = 1), and surprising insulating phases at fractional filling factors {\nu} = 1/3 and 2/3, which we assign to generalized Wigner crystallization on an underlying lattice. Furthermore, the unique spin-valley optical selection rules of TMDC heterostructures allow us to optically create and investigate low-energy spin excited states in the Mott insulator. We reveal an especially slow spin relaxation lifetime of many microseconds in the Mott insulating state, orders-of-magnitude longer than that of charge excitations. Our studies highlight novel correlated physics that can emerge in moir\'e superlattices beyond graphene.


title: Probing the Electronic Structure of Graphene Near and Far from the Fermi Level via Planar Tunneling Spectroscopy

abstract: Scanning tunneling spectroscopy (STS) has yielded significant insight on the electronic structure of graphene and other two-dimensional (2D) materials. STS directly measures a fundamental and directly calculable quantity: the single particle density of states (SPDOS). Due to experimental setup limitations, however, STS has been unable to explore 2D materials in ultra-high magnetic fields where electron-electron interactions can drastically change the SPDOS. Recent developments in the assembly of heterostructures composed of graphene and hexagonal boron nitride have enabled a device-based alternative to potentially overcome these roadblocks. Thus far, however, these nascent efforts are incomplete in analyzing and understanding tunneling spectra and have yet to explore graphene at high magnetic fields. Here we report an experiment at magnetic fields up to 18 T that uses graphene tunneling field effect transistors (TFETs) to establish a clear benchmark for measurement and analysis of graphene planar tunneling spectroscopy. We acquire gate tunable tunneling spectra of graphene and then use these data and electrostatic arguments to develop a systematic analysis scheme. This analysis reveals that TFET devices directly probe electronic structure features near and far from the Fermi level. In particular, our study yields identification of the Dirac point and numerous Landau levels as they fill and empty with charge via application of a gate voltage. Our work demonstrates that TFET devices are a viable platform for directly probing the electronic structure of graphene and other 2D materials in high magnetic fields, where novel electronic states emerge.


title: Applying deep learning to teleseismic phase detection and picking: PcP and PKiKP cases

abstract: The availability of a tremendous amount of seismic data demands seismological researchers to analyze seismic phases efficiently. Recently, deep learning algorithms exhibit a powerful capability of detecting and picking on P- and S-wave phases. However, it is still a challenge to process teleseismic phases fast and accurately. In this study, we detect and pick the PcP and PKiKP phases from a Hinet dataset with 7386 seismograms by applying a deep-learning-based scheme. The scheme consists of three steps: first, we prepare latent phase data, which is truncated from the whole seismogram with the theoretical arrival time; second, we identify and evaluate the latent phase via a convolutional neural network; third, we pick the first break of good or fair phase with a fully convolutional neural network. The detection result shows that the scheme recognizes 92.15% and 94.13% of PcP and PKiKP phases. The picking result has an absolute mean error of 0.0742 s and 0.0636 s for the PcP and PKiKP phases, respectively. The performance of the picking algorithm is compared with the traditional approach of STA/LTA. The scheme processes all 7386 seismograms approximately in 2 hours, especially only cost about five minutes on the last two steps.


title: Topological Factors of Rank-One Subshifts

abstract: We study topological factors of rank-one subshifts and prove that those factors that are themselves subshifts are either finite or isomorphic to the original rank-one subshifts. Thus, we completely characterize the subshift factors of rank-one subshifts.


title: Combinatorial results for network-based models of metabolic origins

abstract: A key step in the origin of life is the emergence of a primitive metabolism. This requires the formation of a subset of chemical reactions that is both self-sustaining and collectively autocatalytic. A generic theory to study such processes (called 'RAF theory') has provided a precise and computationally effective way to address these questions, both on simulated data and in laboratory studies. One of the classic applications of this theory (arising from Stuart Kauffman's pioneering work in the 1980s) involves networks of polymers under cleavage and ligation reactions; in the first part of this paper, we provide the first exact description of the number of such reactions under various model assumptions. Conclusions from earlier studies relied on either approximations or asymptotic counting, and we show that the exact counts lead to similar (though not always identical) asymptotic results. In the second part of the paper, we solve some questions posed in more recent papers concerning the computational complexity of some key questions in RAF theory. In particular, although there is a fast algorithm to determine whether or not a catalytic reaction network contains a subset that is both self-sustaining and autocatalytic (and, if so, find one), determining whether or not sets exist that satisfy certain additional constraints exist turns out to be NP-complete.


title: Autonomous Control of a Quadrotor-Manipulator; Application of Extended State Disturbance Observer

abstract: In this work, the autonomous control of a quadrotor-manipulator unmanned aerial vehicle is treated using an extended dynamic model. Due to persistent aerodynamic disturbances and dynamic couplings, the control of a quadrotor-manipulator unit is cumbersome. Active disturbance rejection control (ADRC) is proposed for the cancellation of the undesirable dynamics while a PD controller is developed for set-point tracking. Parameters were tuned using Gradient Descent-Sequential Quadratic Programming in Simulink design optimization


title: Optimal Control of Wave Energy Converters Using Epsilon-Trig Regularization Method

abstract: The wave energy converter (WEC) devices provide access to a renewable energy source. Developing control strategies to harvest maximum wave energy requires solving a constrained optimal control problem. It is shown that singular control arcs may constitute part (or the entire) of extremal trajectories. Characterizing the optimal control structure, especially with the possibility of many switches between regular and singular control arcs, is challenging due to lack of \textit{a priori} information about: 1) optimal sequence as well as number of the regular and singular control arcs, and 2) the corresponding optimal switch times (from a regular to a singular arc and vice versa). This investigation demonstrates the application of a recently developed construct, the Epsilon-Trig Regularization Method (ETRM), to the problem of maximizing energy harvesting for a point-absorber WEC model in the presence of control constraints. Utility of the ETRM for the WEC problem is demonstrated by comparing its high-quality results against those in the literature for a number of test cases.


title: A theorem of Besicovitch and a generalization of the Birkhoff Ergodic Theorem

abstract: A remarkable theorem of Besicovitch is that an integrable function $f$ on $\mathbb{R}^2$ is strongly differentiable if and only if its associated strong maximal function $M_S f$ is finite a.e. We provide an analogue of Besicovitch's result in the context of ergodic theory that provides a generalization of Birkhoff's Ergodic Theorem. In particular, we show that if $f$ is a measurable function on a standard probability space and $T$ is an invertible measure-preserving transformation on that space, then the ergodic averages of $f$ with respect to $T$ converge a.e. if and only if the associated ergodic maximal function $T^*f$ is finite a.e.


title: Image recognition from raw labels collected without annotators

abstract: Image classification problems are typically addressed by first collecting examples with candidate labels, second cleaning the candidate labels manually, and third training a deep neural network on the clean examples. The manual labeling step is often the most expensive one as it requires workers to label millions of images. In this paper we propose to work without any explicitly labeled data by i) directly training the deep neural network on the noisy candidate labels, and ii) early stopping the training to avoid overfitting. With this procedure we exploit an intriguing property of standard overparameterized convolutional neural networks trained with (stochastic) gradient descent: Clean labels are fitted faster than noisy ones. We consider two classification problems, a subset of ImageNet and CIFAR-10. For both, we construct large candidate datasets without any explicit human annotations, that only contain 10%-50% correctly labeled examples per class. We show that training on the candidate examples and regularizing through early stopping gives higher test performance for both problems than when training on the original, clean data. This is possible because the candidate datasets contain a huge number of clean examples, and, as we show in this paper, the noise generated through the label collection process is not nearly as adversarial for learning as the noise generated by randomly flipping labels.


title: Amortized Rejection Sampling in Universal Probabilistic Programming

abstract: Naive approaches to amortized inference in probabilistic programs with unbounded loops can produce estimators with infinite variance. This is particularly true of importance sampling inference in programs that explicitly include rejection sampling as part of the user-programmed generative procedure. In this paper we develop a new and efficient amortized importance sampling estimator. We prove finite variance of our estimator and empirically demonstrate our method's correctness and efficiency compared to existing alternatives on generative programs containing rejection sampling loops and discuss how to implement our method in a generic probabilistic programming framework.


title: Zero-Shot Recognition via Optimal Transport

abstract: We propose an optimal transport (OT) framework for generalized zero-shot learning (GZSL), seeking to distinguish samples for both seen and unseen classes, with the assist of auxiliary attributes. The discrepancy between features and attributes is minimized by solving an optimal transport problem. {Specifically, we build a conditional generative model to generate features from seen-class attributes, and establish an optimal transport between the distribution of the generated features and that of the real features.} The generative model and the optimal transport are optimized iteratively with an attribute-based regularizer, that further enhances the discriminative power of the generated features. A classifier is learned based on the features generated for both the seen and unseen classes. In addition to generalized zero-shot learning, our framework is also applicable to standard and transductive ZSL problems. Experiments show that our optimal transport-based method outperforms state-of-the-art methods on several benchmark datasets.


title: Deep speech inpainting of time-frequency masks

abstract: Transient loud intrusions, often occurring in noisy environments, can completely overpower speech signal and lead to an inevitable loss of information. While existing algorithms for noise suppression can yield impressive results, their efficacy remains limited for very low signal-to-noise ratios or when parts of the signal are missing. To address these limitations, here we propose an end-to-end framework for speech inpainting, the context-based retrieval of missing or severely distorted parts of time-frequency representation of speech. The framework is based on a convolutional U-Net trained via deep feature losses, obtained using speechVGG, a deep speech feature extractor pre-trained on an auxiliary word classification task. Our evaluation results demonstrate that the proposed framework can recover large portions of missing or distorted time-frequency representation of speech, up to 400 ms and 3.2 kHz in bandwidth. In particular, our approach provided a substantial increase in STOI & PESQ objective metrics of the initially corrupted speech samples. Notably, using deep feature losses to train the framework led to the best results, as compared to conventional approaches.


title: A Nuclear Structure Model for Double Charge-Exchange Processes

abstract: A new model, based on the BCS approach, is specially designed to describe nuclear phenomena $(A,Z)\rightarrow (A,Z\pm 2)$ of double-charge exchange (DCE). After being proposed, and applied in the particle-hole limit, by one of the authors (F. Krmpoti\'c [1]), so far it was never been applied within the BCS mean-field framework, nor has its ability to describe DCE processes been thoroughly explored. It is a natural extension of the pn-QRPA model, developed by Halbleib and Sorensen [2] to describe the single $\beta$-decays $(A,Z)\rightarrow (A,Z\pm 1)$, to the DCE processes. As such, it exhibits several advantages over the pn-QRPA model when is used in the evaluation of the double beta decay (DBD) rates. For instance, i) the extreme sensitivity of the nuclear matrix elements (NMEs) on the model parametrization does not occur, ii) it allows to study NMEs, not only for the fundamental state in daughter nuclei, as the pn-QRPA model does, but also for all final $0^+$ and $2^+$ states, accounting at the same time their excitation energies and the corresponding DBD Q-values, iii) together with the DBD-NMEs it provides also the energy spectra of Fermi and Gamow-Teller DCE transition strengths, as well as the locations of the corresponding resonances and their sum rules, iv) the latter are relevant for both the DBD and the DCE reactions, since the involved nuclear structure is the same; this correlation does not exist within the pn-QRPA model. As an example, detailed numerical calculations are presented for the $(A,Z)\rightarrow (A,Z+ 2)$ process in $^{48}$Ca $\rightarrow ^{48}$Ti and the $(A,Z)\rightarrow (A,Z- 2)$ process in $^{96}$Ru $\rightarrow ^{96}$Mo, involving all final $0^+$ states and $2^+$ states.


title: An Online Deep Learning Approach Toward the Prediction of Power System Stresses Using Voltage Phasors

abstract: The outage of a transmission line may change the system phase angle differences to the point that the system experience stress conditions. Hence, the angle differences for post-contingency condition of a transmission lines should be predicted in real time operation. However, online line-based phase angle difference monitoring and prediction for power system stress assessment is not a universal operating practice yet. Thus, in this paper, an online power system stress assessment framework is proposed by developing a convolutional neural network (CNN) module trained through Deep Learning approach. In the proposed framework, the continuously streaming system phase angle data, driven from phasor measurement units (PMUs) or a state estimator (SE), is used to construct power system stress indices adaptive to the structure parameters of the CNN module. Using this approach, any hidden patterns between phase angles of buses and system stress conditions are revealed at low computation cost while yielding accurate stress status and the severity of the stress. The effectiveness and scalability of the proposed method has been verified on the IEEE 118-bus and more importantly, on the PJM Interconnection system. Moreover, outperformance of the proposed method is verified by comparing the results with artificial neural network (ANN) and decision tree (DT).


title: Deep Mouse: An End-to-end Auto-context Refinement Framework for Brain Ventricle and Body Segmentation in Embryonic Mice Ultrasound Volumes

abstract: High-frequency ultrasound (HFU) is well suited for imaging embryonic mice due to its noninvasive and real-time characteristics. However, manual segmentation of the brain ventricles (BVs) and body requires substantial time and expertise. This work proposes a novel deep learning based end-to-end auto-context refinement framework, consisting of two stages. The first stage produces a low resolution segmentation of the BV and body simultaneously. The resulting probability map for each object (BV or body) is then used to crop a region of interest (ROI) around the target object in both the original image and the probability map to provide context to the refinement segmentation network. Joint training of the two stages provides significant improvement in Dice Similarity Coefficient (DSC) over using only the first stage (0.818 to 0.906 for the BV, and 0.919 to 0.934 for the body). The proposed method significantly reduces the inference time (102.36 to 0.09 s/volume around 1000x faster) while slightly improves the segmentation accuracy over the previous methods using slide-window approaches.


title: Note on the Delta Method for Finite Population Inference with Applications to Causal Inference

abstract: The delta method creates more general inference results when coupled with central limit theorem results for the finite population. This opens up a range of new estimators for which we can find finite population asymptotic properties. We focus on the use of this method to derive asymptotic distributional results and variance expressions for causal estimators. We illustrate the use of the method by obtaining a finite population asymptotic distribution for a causal ratio estimator.


title: Topological Sliding Moir\'e Heterostructure

abstract: We investigate the effect of sliding motion of layers in Moir\'e heterostructures on the electronic state. We show that the sliding Moir\'e heterostructure can generate nontrivial topology characterized by the first and second Chern number in the high dimensional manifold spanned by the physical dimensions and synthetic dimensions associated with the sliding displacement. The nontrivial topology implies a topological charge pumping caused by the sliding motion. We demonstrate the nontrivial topology and charge pumping explicitly in a one dimensional bi-chain model and the small-angle twisted bilayer graphene. Contrary to the conventional belief that the interlayer sliding in incommensurate Moir\'e Heterostructures does not affect the electronic structure, our results reveal that the sliding motion can generate nontrivial topology dynamically and hence cannot be neglected in the dynamical process.


title: Exoplanet direct imaging in ground-based conditions on THD2 bench

abstract: The next generation of ground-based instruments aims to break through the knowledge we have on exoplanets by imaging circumstellar environments always closer to the stars. However, direct imaging requires an AO system and high-contrast techniques like a coronagraph to reject the diffracted light of an observed star and an additional wavefront sensor to control quasi-static aberrations, including the non common path aberrations. To observe faint objects, a focal plane wavefront sensor with a sub-nanometric wavefront control capability is required. In the past few years, we developed the THD2 bench which is a testbed for high-contrast imaging techniques, working in visible and near infrared wavelengths and currently reaching contrast levels lower than 1e-8 under space-like simulated conditions. We recently added a turbulence wheel on the optical path which simulates the residuals given by a typical extreme adaptive optics system and we tested several ways to remove quasi-statics speckles. One way to estimate the aberrations is a method called pair-wise probing where we record few images with known-shapes we apply on the adaptive optics deformable mirror. Once estimated, we seek to minimize the focal-plane electric field by an algorithm called Electric Field Conjugation. In this paper, we present the first results obtained on the THD2 bench using these two techniques together in turbulent conditions. We then compare the achieved performance with the one expected when all the quasi-static speckles are corrected.


title: Neutrinos from the gamma-ray source eHWC J1825-134: predictions for Km$^3$ detectors

abstract: The eHWC J1825-134 source is located in the southern sky and has been recently detected by the HAWC observatory. It presents an hard spectral index and its gamma-ray flux extends up to energies close to 100 TeV without significant suppression. Amongst the HAWC sources, it is the most luminous in the multi-TeV domain and therefore is one of the first that should be searched for with a neutrino telescope in the northern hemisphere. Using an updated effective area for the forthcoming KM3NeT detector, we study the possibility to detect this source within ten years of its running time. We conclude that about a 4 to 5 sigma detection has to be expected after ten years of observations, depending on the details of the considered scenario.


title: From Importance Sampling to Doubly Robust Policy Gradient

abstract: We show that on-policy policy gradient (PG) and its variance reduction variants can be derived by taking finite difference of function evaluations supplied by estimators from the importance sampling (IS) family for off-policy evaluation (OPE). Starting from the doubly robust (DR) estimator (Jiang & Li, 2016), we provide a simple derivation of a very general and flexible form of PG, which subsumes the state-of-the-art variance reduction technique (Cheng et al., 2019) as its special case and immediately hints at further variance reduction opportunities overlooked by existing literature. We analyze the variance of the new DR-PG estimator, compare it to existing methods as well as the Cramer-Rao lower bound of policy gradient, and empirically show its effectiveness.


title: A Quantitative Physical Model of the Snow Crystal Morphology Diagram

abstract: I describe a semi-empirical molecular model of the surface attachment kinetics governing ice crystal growth from water vapor as a function of temperature, supersaturation, and crystal mesostructure. An important new hypothesis in this model is surface-energy-driven molecular diffusion enabled by a leaky Ehrlich-Schwoebel barrier. The proposed surface-diffusion behavior is sensitive to facet width and surface premelting, yielding structure-dependent attachment kinetics with a complex temperature dependence. By incorporating several reasonable assumptions regarding the surface premelting behavior on basal and prism facets, this model can explain the overarching features of the snow crystal morphology diagram, which has been an enduring scientific puzzle for nearly 75 years.


title: Relational Test Tables: A Practical Specification Language for Evolution and Security

abstract: A wide range of interesting program properties are intrinsically relational, i.e., they relate two or more program traces. Two prominent relational properties are secure information flow and conditional program equivalence. By showing the absence of illegal information flow, confidentiality and integrity properties can be proved. Equivalence proofs allow using an existing (trusted) software release as specification for new revisions. Currently, the verification of relational properties is hardly accessible to practitioners due to the lack of appropriate relational specification languages. In previous work, we introduced the concept of generalised test tables: a table-based specification language for functional (non-relational) properties of reactive systems. In this paper, we present relational test tables -- a canonical extension of generalised test tables for the specification of relational properties, which refer to two or more program runs or traces. Regression test tables support asynchronous program runs via stuttering. We show the applicability of relational test tables, using them for the specification and verification of two examples from the domain of automated product systems.


title: Large sieve inequality for power moduli

abstract: In this note we give a new bound for large sieve with characters to power moduli which improves in some range of the parameters the previous bounds of Baier/Zhao and Halupczok.


title: Structured Prediction Helps 3D Human Motion Modelling

abstract: Human motion prediction is a challenging and important task in many computer vision application domains. Existing work only implicitly models the spatial structure of the human skeleton. In this paper, we propose a novel approach that decomposes the prediction into individual joints by means of a structured prediction layer that explicitly models the joint dependencies. This is implemented via a hierarchy of small-sized neural networks connected analogously to the kinematic chains in the human body as well as a joint-wise decomposition in the loss function. The proposed layer is agnostic to the underlying network and can be used with existing architectures for motion modelling. Prior work typically leverages the H3.6M dataset. We show that some state-of-the-art techniques do not perform well when trained and tested on AMASS, a recently released dataset 14 times the size of H3.6M. Our experiments indicate that the proposed layer increases the performance of motion forecasting irrespective of the base network, joint-angle representation, and prediction horizon. We furthermore show that the layer also improves motion predictions qualitatively. We make code and models publicly available at https://ait.ethz.ch/projects/2019/spl.


title: Classical algorithms, correlation decay, and complex zeros of partition functions of quantum many-body systems

abstract: In this paper, we present a quasi-polynomial time classical algorithm that estimates the partition function of quantum many-body systems at temperatures above the thermal phase transition point. It is known that in the worst case, the same problem is NP-hard below this point. Together with our work, this shows that the transition in the phase of a quantum system is also accompanied by a transition in the hardness of approximation. We also show that in a system of n particles above the phase transition point, the correlation between two observables whose distance is at least log(n) decays exponentially. We can improve the factor of log(n) to a constant when the Hamiltonian has commuting terms or is on a 1D chain. The key to our results is a characterization of the phase transition and the critical behavior of the system in terms of the complex zeros of the partition function. Our work extends a seminal work of Dobrushin and Shlosman on the equivalence between the decay of correlations and the analyticity of the free energy in classical spin models. On the algorithmic side, our result extends the scope of a recent approach due to Barvinok for solving classical counting problems to quantum many-body systems.


title: Regular Model Checking with Regular Relations

abstract: Regular model checking is an exploration technique for infinite state systems where state spaces are represented as regular languages and transition relations are expressed using rational relations over infinite (or finite) strings. We extend the regular model checking paradigm to permit the use of more powerful transition relations: the class of regular relations, of which the rational relations are a strict subset. We use the language of monadic second-order logic on infinite strings to specify such relations and adopt streaming string transducers (SSTs) as a suitable computational model. We introduce nondeterministic SSTs over infinite strings and show that they precisely capture the relations definable in MSO. We further explore theoretical properties of omega-NSSTs required to effectively carry out regular model checking. In particular, we establish that the regular type checking problem for omega-NSSTs is decidable in PSPACE. Since the post-image of a regular language under a regular relation may not be regular (or even context-free), approaches that iteratively compute the image can not be effectively carried out in this setting. Instead, we utilize the fact that regular relations are closed under composition, which, together with our decidability result, provides a foundation for regular model checking with regular relations.


title: Anthropic Selection of Physical Constants, Quantum Entanglement, and the Multiverse Falsifiability

abstract: This paper evaluates some important aspects of the multiverse concept. Firstly, the most realistic opportunity for it which is the spacetime variability of the physical constants and may deliver worlds with different physics, hopefully fulfilling the conditions of the anthropic principles. Then, more esoteric versions of the multiverse being the realisation of some abstract mathematics or even logic. Finally, it evaluates the big challenge of getting any signal from "other universes" using recent achievements of the quantum theory.


title: A survey on agile practices and challenges of a global software development team

abstract: The Agile Manifesto describes that the most efficient and effective method of conveying information to and within a development team is through face-to-face conversation. However that is not always possible when teams are working in a Global Software Development (GSD) environment. Based on this scenario, this study presents an exploratory data analysis using survey results to explore agile practices and challenges of a global software development team that uses ScaledAgile Framework (SAFe), which is designed for the need of larger organizations. The goal of this study is to understand the team's level of knowledge in some agile practices and which types of communication are usually prioritized. As in GSD environments team members are geographically spread across multiple regions and time zones, we aim to identify challenges this environment can present. As a result of this exploratory analysis, it has been identified that communication is one of main challenges in GSD environment and that phone calls are considered to be the most efficient type of communication. Additionally, we have also identified that professionals have different levels of confidence in Agile practices and concluded that knowledge transfers among the professionals could help those team members that are not confident in some agile practices to increase their overall confidence and knowledge.


title: Bruno Touschek in Germany after the War: 1945-46

abstract: Bruno Touschek was an Austrian born theoretical physicist, who proposed and built the first electron-positron collider in 1960 in the Frascati National Laboratories in Italy. In this note we reconstruct a crucial period of Bruno Touschek's life so far scarcely explored, which runs from Summer 1945 to the end of 1946. We shall describe his university studies in G\"ottingen, placing them in the context of the reconstruction of German science after 1945. The influence of Werner Heisenberg and other prominent German physicists will be highlighted. In parallel, we shall show how the decisions of the Allied powers towards restructuring science and technology in the UK after the war effort, determined Touschek's move to the University of Glasgow in 1947.


title: Poincar\'e/Koszul Duality for General Operads

abstract: We record a result concerning the Koszul dual of the arity filtration on an operad. This result is then used to give conditions under which, for a general operad, the Poincar\'e/Koszul duality arrow of Ayala-Francis is an equivalence. We discuss how the Poincar\'e/Koszul duality arrow for the little $n$-disks operad $\mathcal{E}_n$ relates to their work when combined with the self-Koszul duality of $\mathcal{E}_n$.


title: Looking Ahead: Anticipating Pedestrians Crossing with Future Frames Prediction

abstract: In this paper, we present an end-to-end future-prediction model that focuses on pedestrian safety. Specifically, our model uses previous video frames, recorded from the perspective of the vehicle, to predict if a pedestrian will cross in front of the vehicle. The long term goal of this work is to design a fully autonomous system that acts and reacts as a defensive human driver would --- predicting future events and reacting to mitigate risk. We focus on pedestrian-vehicle interactions because of the high risk of harm to the pedestrian if their actions are miss-predicted. Our end-to-end model consists of two stages: the first stage is an encoder/decoder network that learns to predict future video frames. The second stage is a deep spatio-temporal network that utilizes the predicted frames of the first stage to predict the pedestrian's future action. Our system achieves state-of-the-art accuracy on pedestrian behavior prediction and future frames prediction on the Joint Attention for Autonomous Driving (JAAD) dataset.


title: 2I/Borisov: A C$_2$ depleted interstellar comet

abstract: The discovery of the first active interstellar object 2I/Borisov provides an unprecedented opportunity to study planetary formation processes in another planetary system. In particular, spectroscopic observations of 2I allow us to constrain the composition of its nuclear ices. We obtained optical spectra of 2I with the 4.2 m William Herschel and 2.5 m Isaac Newton telescopes between 2019 September 30 and October 13, when the comet was between 2.5 au and 2.4 au from the Sun. We also imaged the comet with broadband filters on 15 nights from September 11 to October 17, as well as with a CN narrow-band filter on October 18 and 20, with the TRAPPIST-North telescope. Broadband imaging confirms that the dust coma colours (B-V=0.82$\pm$0.02, V-R=0.46$\pm$0.03, R-I=0.44$\pm$0.03, B-R=1.28$\pm$0.03) are the same as for Solar System comets. We detect CN emission in all spectra and in the TRAPPIST narrow-band images with production rates between 1.6$\times10^{24}$ and 2.1$\times10^{24}$ molec/s. No other species are detected. We determine three-sigma upper limits for C$_2$, C$_3$, and OH production rates of 6$\times10^{23}$ molec/s, 3$\times10^{23}$ molec/s and 2$\times10^{27}$ molec/s, respectively, on October 02. There is no significant increase of the CN production rate or A(0)f$\rho$ during our observing period. Finally, we place a three-sigma upper limit on the Q(C$_2$)/Q(CN) ratio of 0.3 (on October 13). From this, we conclude that 2I is highly depleted in C$_2$, and may have a composition similar to Solar System carbon-chain depleted comets.


title: Horizon-T Experiment Detailed Calibration of Cables

abstract: The ability to extract the pulse width and translate it into the actual disk width of the Extensive Air Showers (EAS) is a hard one requiring accurate knowledge of the system performance. For that, the analysis for the cable calibration for Horizon-10T detectors has been re-analyzed in a different form that allows for better signal width measurements. An innovative detector system Horizon-10T, constructed to study EAS in the energy range above 1016 eV coming from a wide range of zenith angles (0o - 85o), is located at Tien Shan high-altitude Science Station of Lebedev Physical Institute of the Russian Academy of Sciences at approximately 3340 meters above the sea level.


title: Error estimate of a bi-fidelity method for kinetic equations with random parameters and multiple scales

abstract: In this paper, we conduct uniform error estimates of the bi-fidelity method for multi-scale kinetic equations. We take the Boltzmann and the linear transport equations as important examples. The main analytic tool is the hypocoercivity analysis for kinetic equations, considering solutions in a perturbative setting close to the global equilibrium. This allows us to obtain the error estimates in both kinetic and hydrodynamic regimes.


title: An Application of Abel's Method to the Inverse Radon Transform

abstract: A method of approximating the inverse Radon transform on the plane by integrating against a smooth kernel is investigated. For piecewise smooth integrable functions, convergence theorems are proven and Gibbs phenomena are ruled out. Geometric properties of the kernel and their implications for computer implementation are discussed. Suggestions are made for applications and an example is presented.


title: Coexistence of type-I and type-II superconductivity signatures in ZrB12 probed by muon spin rotation measurements

abstract: Superconductors usually display either type-I or type-II superconductivity and the coexistence of these two types in the same material, for example at different temperatures is rare in nature. We the employed muon spin rotation (muSR) technique to unveil the superconducting phase diagram of the dodecaboride ZrB12 and obtained clear evidence of both type-I and type-II characteristics. Most importantly, we found a region showing unusual behavior where the usually mutually exclusive muSR signatures of type-I and type-II superconductivity coexist. We reproduced that behavior in theoretical modeling that required taking into account multiple bands and multiple coherence lengths, which suggests that material has one coherence length larger and another smaller than the magnetic field penetration length (the type-1.5 regime). At stronger fields, a footprint of the type-II mixed state showing square flux-line lattice was also obtained using neutron diffraction.


title: Spectral CUSUM for Online Network Structure Change Detection

abstract: Detecting abrupt changes in the community structure of a network from noisy observations is a fundamental problem in statistics and machine learning. This paper presents an online change detection algorithm called Spectral-CUSUM to detect unknown network structure changes through a generalized likelihood ratio statistic. We characterize the average run length (ARL) and the expected detection delay (EDD) of the Spectral-CUSUM procedure and prove its asymptotic optimality. Finally, we demonstrate the good performance of the Spectral-CUSUM procedure and compare it with several baseline methods using simulations and real data examples on seismic event detection using sensor network data.


title: Workshop Report: Workshop on Opportunities, Challenges, and Best Practices for Basic Plasma Science User Facilities (May 20-21, 2019, University of Maryland, College Park MD)

abstract: This is the report on the Workshop on Opportunities, Challenges, and Best Practices for Basic Plasma Science User Facilities, held at the University of Maryland, College Park, MD, on May 20-21, 2019.


title: Semantics for Global and Local Interpretation of Deep Neural Networks

abstract: Deep neural networks (DNNs) with high expressiveness have achieved state-of-the-art performance in many tasks. However, their distributed feature representations are difficult to interpret semantically. In this work, human-interpretable semantic concepts are associated with vectors in feature space. The association process is mathematically formulated as an optimization problem. The semantic vectors obtained from the optimal solution are applied to interpret deep neural networks globally and locally. The global interpretations are useful to understand the knowledge learned by DNNs. The interpretation of local behaviors can help to understand individual decisions made by DNNs better. The empirical experiments demonstrate how to use identified semantics to interpret the existing DNNs.


title: Contextual Prediction Difference Analysis for Explaining Individual Image Classifications

abstract: Much effort has been devoted to understanding the decisions of deep neural networks in recent years. A number of model-aware saliency methods were proposed to explain individual classification decisions by creating saliency maps. However, they are not applicable when the parameters and the gradients of the underlying models are unavailable. Recently, model-agnostic methods have also received attention. As one of them, \textit{Prediction Difference Analysis} (PDA), a probabilistic sound methodology, was proposed. In this work, we first show that PDA can suffer from saturated classifiers. The saturation phenomenon of classifiers exists widely in current neural network-based classifiers. To explain the decisions of saturated classifiers better, we further propose Contextual PDA, which runs hundreds of times faster than PDA. The experiments show the superiority of our method by explaining image classifications of the state-of-the-art deep convolutional neural networks.


title: Highly efficient schemes for time fractional Allen-Cahn equation using extended SAV approach

abstract: In this paper, we propose and analyze high order efficient schemes for the time fractional Allen-Cahn equation. The proposed schemes are based on the L1 discretization for the time fractional derivative and the extended scalar auxiliary variable (SAV) approach developed very recently to deal with the nonlinear terms in the equation. The main contributions of the paper consist in: 1) constructing first and higher order unconditionally stable schemes for different mesh types, and proving the unconditional stability of the constructed schemes for the uniform mesh; 2) carrying out numerical experiments to verify the efficiency of the schemes and to investigate the coarsening dynamics governed by the time fractional Allen-Cahn equation. Particularly, the influence of the fractional order on the coarsening behavior is carefully examined. Our numerical evidence shows that the proposed schemes are more robust than the existing methods, and their efficiency is less restricted to particular forms of the nonlinear potentials.


title: Equilibrium microstructures of diblock copolymers under 3D confinement

abstract: We investigate equilibrium microstructures exhibited by diblock copolymers in confined 3D geometries. We perform Self-Consistent Field Theory (SCFT) simulations using a finite-element based computational framework (Ackerman et al. 2017), that provides the flexibility to compute equilibrium structures under arbitrary geometries. We consider a sequence of 3D geometries (tetrahedron to sphere) that have the same volume but exhibit varying curvature. This allows us to study the interplay between edge and curvature effects of the 3D geometries on the equilibrium microstructures. We observe that beyond a length scale, the equilibrium structure changes from an interconnected network to a multi-layered concentric shell as the curvature of the 3D geometry is reduced. However, below this length scale the equilibrium structure remains a multi-layered concentric shell independent of curvature. We additionally explore variations in the equilibrium microstructures at a few discrete volume fractions. This study provides insight into possible frustrated phases that can arise in AB diblock systems while varying the shape of confinement geometry.

categories: cond-mat.soft cond-mat.mes-hall cond-mat.mtrl-sci physics.comp-ph


title: Decentralized Heterogeneous Multi-Player Multi-Armed Bandits with Non-Zero Rewards on Collisions

abstract: We consider a fully decentralized multi-player stochastic multi-armed bandit setting where the players cannot communicate with each other and can observe only their own actions and rewards. The environment may appear differently to different players, $\textit{i.e.}$, the reward distributions for a given arm are heterogeneous across players. In the case of a collision (when more than one player plays the same arm), we allow for the colliding players to receive non-zero rewards. The time-horizon $T$ for which the arms are played is \emph{not} known to the players. Within this setup, where the number of players is allowed to be greater than the number of arms, we present a policy that achieves near order-optimal expected regret of order $O(\log^{1 + \delta} T)$ for some $0 < \delta < 1$ over a time-horizon of duration $T$. This paper is accepted at IEEE Transactions on Information Theory.


title: A game method for improving the interpretability of convolution neural network

abstract: Real artificial intelligence always has been focused on by many machine learning researchers, especially in the area of deep learning. However deep neural network is hard to be understood and explained, and sometimes, even metaphysics. The reason is, we believe that: the network is essentially a perceptual model. Therefore, we believe that in order to complete complex intelligent activities from simple perception, it is necessary to con-struct another interpretable logical network to form accurate and reasonable responses and explanations to external things. Researchers like Bolei Zhou and Quanshi Zhang have found many explanatory rules for deep feature extraction aimed at the feature extraction stage of convolution neural network. However, although researchers like Marco Gori have also made great efforts to improve the interpretability of the fully connected layers of the network, the problem is also very difficult. This paper firstly analyzes its reason. Then a method of constructing logical network based on the fully connected layers and extracting logical relation between input and output of the layers is proposed. The game process between perceptual learning and logical abstract cognitive learning is implemented to improve the interpretable performance of deep learning process and deep learning model. The benefits of our approach are illustrated on benchmark data sets and in real-world experiments.


title: Multi-User MABs with User Dependent Rewards for Uncoordinated Spectrum Access

abstract: Multi-user multi-armed bandits have emerged as a good model for uncoordinated spectrum access problems. In this paper we consider the scenario where users cannot communicate with each other. In addition, the environment may appear differently to different users, ${i.e.}$, the mean rewards as observed by different users for the same channel may be different. With this setup, we present a policy that achieves a regret of $O (\log{T})$. This paper has been accepted at Asilomar Conference on Signals, Systems, and Computers 2019.


title: Fast Exact Matrix Completion: A Unified Optimization Framework for Matrix Completion

abstract: We formulate the problem of matrix completion with and without side information as a non-convex optimization problem. We design fastImpute based on non-convex gradient descent and show it converges to a global minimum that is guaranteed to recover closely the underlying matrix while it scales to matrices of sizes beyond $10^5 \times 10^5$. We report experiments on both synthetic and real-world datasets that show fastImpute is competitive in both the accuracy of the matrix recovered and the time needed across all cases. Furthermore, when a high number of entries are missing, fastImpute is over $75\%$ lower in MAPE and $15$ times faster than current state-of-the-art matrix completion methods in both the case with side information and without.


title: All-Action Policy Gradient Methods: A Numerical Integration Approach

abstract: While often stated as an instance of the likelihood ratio trick [Rubinstein, 1989], the original policy gradient theorem [Sutton, 1999] involves an integral over the action space. When this integral can be computed, the resulting "all-action" estimator [Sutton, 2001] provides a conditioning effect [Bratley, 1987] reducing the variance significantly compared to the REINFORCE estimator [Williams, 1992]. In this paper, we adopt a numerical integration perspective to broaden the applicability of the all-action estimator to general spaces and to any function class for the policy or critic components, beyond the Gaussian case considered by [Ciosek, 2018]. In addition, we provide a new theoretical result on the effect of using a biased critic which offers more guidance than the previous "compatible features" condition of [Sutton, 1999]. We demonstrate the benefit of our approach in continuous control tasks with nonlinear function approximation. Our results show improved performance and sample efficiency.


title: Self-supervised classification of dynamic obstacles using the temporal information provided by videos

abstract: Nowadays, autonomous driving systems can detect, segment, and classify the surrounding obstacles using a monocular camera. However, state-of-the-art methods solving these tasks generally perform a fully supervised learning process and require a large amount of training labeled data. On another note, some self-supervised learning approaches can deal with detection and segmentation of dynamic obstacles using the temporal information available in video sequences. In this work, we propose to classify the detected obstacles depending on their motion pattern. We present a novel self-supervised framework consisting of learning offline clusters from temporal patch sequences and considering these clusters as labeled sets to train a real-time image classifier. The presented model outperforms state-of-the-art unsupervised image classification methods on large-scale diverse driving video dataset BDD100K.


title: Phase noise in modular millimeter wave massive MIMO

abstract: This paper investigates the effect of oscillator phase noise on a multiuser millimeter wave (mmWave) massive MIMO uplink as we scale up the number of base station antennas, fixing the load factor, defined as the ratio of the number of simultaneous users to the number of base station antennas. We consider a modular approach in which the base station employs an array of subarrays, or "tiles." Each tile supports a fixed number of antennas, and can therefore be implemented using a separate radio frequency integrated circuit (RFIC), with synchronization across tiles accomplished by employing a phased locked loop in each tile to synthesize an on-chip oscillator at the carrier frequency by locking on to a common lower frequency reference clock. Assuming linear minimum mean squared error (LMMSE) multiuser detection, we provide an analytical framework that can be used to specify the required power spectral density (PSD) mask for phase noise for a target system performance. Our analysis for the phase noise at the output of the LMMSE receiver indicates two distinct effects: self-noise for each user which is inversely proportional to the number of tiles, and cross-talk between users which is insensitive to the number of tiles, and is proportional to the load factor. These analytical predictions, verified by simulations for a 140 GHz system targeting a per-user data rate of 10 Gbps, show that tiling is a robust approach for scaling. Numerical results for our proposed design approach yield relatively relaxed specifications for phase noise PSD masks.


title: Entanglement-based single-shot detection of a single magnon with a superconducting qubit

abstract: The recent development of hybrid systems based on superconducting circuits has opened up the possibility of engineering sensors of quanta of different degrees of freedom. Quantum magnonics, which aims to control and read out quanta of collective spin excitations in magnetically-ordered systems, furthermore provides unique opportunities for advances in both the study of magnetism and the development of quantum technologies. Using a superconducting qubit as a quantum sensor, we report the detection of a single magnon in a millimeter-sized ferromagnetic crystal with a quantum efficiency of up to~$0.71$. The detection is based on the entanglement between a magnetostatic mode and the qubit, followed by a single-shot measurement of the qubit state. This proof-of-principle experiment establishes the single-photon detector counterpart for magnonics.


title: Absolute Properties of the Eclipsing $\gamma$ Dor Star V404 Lyrae

abstract: We present the first high-resolution spectra for the eclipsing binary V404 Lyr showing $\gamma$ Dor pulsations, which we use to study its absolute properties. By fitting models to the disentangling spectrum of the primary star, we found that it has an effective temperature of $T_{\rm eff,1}=7,330 \pm 150$ K and a rotational velocity of $v_1\sin$$i=148\pm18$ km s$^{-1}$. The simultaneous analysis of our double-lined radial velocities and the pulsation-subtracted ${\it Kepler}$ data gives us accurate stellar and system parameters of V404 Lyr. The masses, radii, and luminosities are $M_1$ = 2.17$\pm$0.06 M$_\odot$, $R_1$ = 1.91$\pm$0.02 R$_\odot$, and $L_1$ = 9.4$\pm$0.8 L$_\odot$ for the primary, and $M_2$ = 1.42$\pm$0.04 M$_\odot$, $R_2$ = 1.79$\pm$0.02 R$_\odot$, and $L_2$ = 2.9$\pm$0.2 L$_\odot$ for the secondary. The tertiary component orbiting the eclipsing pair has a mass of $M_{\rm 3b}$ = 0.71$\pm$0.15 $M_\odot$ in an orbit of $P_{\rm 3b}$ = 642$\pm$3 days, $e_{\rm 3b}$ = 0.21$\pm$0.04, and $a_{\rm 3b}$ = 509$\pm$2 R$_\odot$. The third light of $l_3=4.1\pm0.2\%$ could be partly attributable to the K-type circumbinary object. By applying a multiple frequency analysis to the eclipse-subtracted light residuals, we detected 45 frequencies with signal to noise amplitude ratios larger than 4.0. Identified as independent pulsation modes, seven frequencies ($f_1-f_6$, $f_9$), their new pulsation constants, and the location in the Hertzsprung-Russell diagram indicate that the pulsating primary is a $\gamma$ Dor-type variable star.


title: Learning and Meta-Learning of Stochastic Advection-Diffusion-Reaction Systems from Sparse Measurements

abstract: Physics-informed neural networks (PINNs) were recently proposed in [1] as an alternative way to solve partial differential equations (PDEs). A neural network (NN) represents the solution while a PDE-induced NN is coupled to the solution NN, and all differential operators are treated using automatic differentiation. Here, we first employ the standard PINN and a stochastic version, sPINN, to solve forward and inverse problems governed by a nonlinear advection-diffusion-reaction (ADR) equation, assuming we have some sparse measurements of the concentration field at random or pre-selected locations. Subsequently, we attempt to optimize the hyper-parameters of sPINN by using the Bayesian optimization method (meta-learning), and compare the results with the empirically selected hyper-parameters of sPINN.


title: Log-dependent slope of scalar induced gravitational waves in the infrared regions

abstract: We analytically calculate the scalar induced gravitational waves (SIGWs) and find a log-dependent slope of SIGW in the infrared regions $(f<f_c)$, namely $n_{\mathrm{GW}}(f)=3-2/\ln(f_c/f)$, and $n_{\mathrm{GW}}(f)=2-2/\ln(f_c/f)$ near the peak if the power spectrum of scalar curvature perturbation is quite narrow, where $f_c$ is roughly the frequency at the peak of SIGW. Such a log-dependent slope can be taken as a new template for distinguishing SIGW from other sources.


title: Biologic and Prognostic Feature Scores from Whole-Slide Histology Images Using Deep Learning

abstract: Histopathology is a reflection of the molecular changes and provides prognostic phenotypes representing the disease progression. In this study, we introduced feature scores generated from hematoxylin and eosin histology images based on deep learning (DL) models developed for prostate pathology. We demonstrated that these feature scores were significantly prognostic for time to event endpoints (biochemical recurrence and cancer-specific survival) and had simultaneously molecular biologic associations to relevant genomic alterations and molecular subtypes using already trained DL models that were not previously exposed to the datasets of the current study. Further, we discussed the potential of such feature scores to improve the current tumor grading system and the challenges that are associated with tumor heterogeneity and the development of prognostic models from histology images. Our findings uncover the potential of feature scores from histology images as digital biomarkers in precision medicine and as an expanding utility for digital pathology.


title: On a class of elliptic functions associated with the even Dirichlet characters

abstract: We construct a class of companion elliptic functions associated with the even Dirichlet characters. Using the well-known properties of the classical Weierstrass elliptic function $\wp(z|\tau)$ as the blueprint, we will derive their representations in terms of $q$-series and partial fractions, explore the significance of the coefficients of their power series expansions and establish the modular properties under the actions of the arithmetic groups $\Gamma_0(N)$ and $\Gamma_1(N)$.