title: Noisy Concurrent Training for Efficient Learning under Label Noise

abstract: Deep neural networks (DNNs) fail to learn effectively under label noise and have been shown to memorize random labels which affect their generalization performance. We consider learning in isolation, using one-hot encoded labels as the sole source of supervision, and a lack of regularization to discourage memorization as the major shortcomings of the standard training procedure. Thus, we propose Noisy Concurrent Training (NCT) which leverages collaborative learning to use the consensus between two models as an additional source of supervision. Furthermore, inspired by trial-to-trial variability in the brain, we propose a counter-intuitive regularization technique, target variability, which entails randomly changing the labels of a percentage of training samples in each batch as a deterrent to memorization and over-generalization in DNNs. Target variability is applied independently to each model to keep them diverged and avoid the confirmation bias. As DNNs tend to prioritize learning simple patterns first before memorizing the noisy labels, we employ a dynamic learning scheme whereby as the training progresses, the two models increasingly rely more on their consensus. NCT also progressively increases the target variability to avoid memorization in later stages. We demonstrate the effectiveness of our approach on both synthetic and real-world noisy benchmark datasets.


title: LAAT: Locally Aligned Ant Technique for discovering multiple faint low dimensional structures of varying density

abstract: Dimensionality reduction and clustering are often used as preliminary steps for many complex machine learning tasks. The presence of noise and outliers can deteriorate the performance of such preprocessing and therefore impair the subsequent analysis tremendously. In manifold learning, several studies indicate solutions for removing background noise or noise close to the structure when the density is substantially higher than that exhibited by the noise. However, in many applications, including astronomical datasets, the density varies alongside manifolds that are buried in a noisy background. We propose a novel method to extract manifolds in the presence of noise based on the idea of Ant colony optimization. In contrast to the existing random walk solutions, our technique captures points that are locally aligned with major directions of the manifold. Moreover, we empirically show that the biologically inspired formulation of ant pheromone reinforces this behavior enabling it to recover multiple manifolds embedded in extremely noisy data clouds. The algorithm performance in comparison to state-of-the-art approaches for noise reduction in manifold detection and clustering is demonstrated, on several synthetic and real datasets, including an N-body simulation of a cosmological volume.


title: Berrut Approximated Coded Computing: Straggler Resistance Beyond Polynomial Computing

abstract: One of the major challenges in using distributed learning to train complicated models with large data sets is to deal with stragglers effect. As a solution, coded computation has been recently proposed to efficiently add redundancy to the computation tasks. In this technique, coding is used across data sets, and computation is done over coded data, such that the results of an arbitrary subset of worker nodes with a certain size are enough to recover the final results. The major challenges with those approaches are (1) they are limited to polynomial function computations, (2) the size of the subset of servers that we need to wait for grows with the multiplication of the size of the data set and the model complexity (the degree of the polynomial), which can be prohibitively large, (3) they are not numerically stable for computation over real numbers. In this paper, we propose Berrut Approximated Coded Computing (BACC), as an alternative approach, which is not limited to polynomial function computation. In addition, the master node can approximately calculate the final results, using the outcomes of any arbitrary subset of available worker nodes. The approximation approach is proven to be numerically stable with low computational complexity. In addition, the accuracy of the approximation is established theoretically and verified by simulation results in different settings such as distributed learning problems. In particular, BACC is used to train a deep neural network on a cluster of servers, which outperforms repetitive computation (repetition coding) in terms of the rate of convergence.


title: Review: Deep Learning in Electron Microscopy

abstract: Deep learning is transforming most areas of science and technology, including electron microscopy. This review paper offers a practical perspective aimed at developers with limited familiarity. For context, we review popular applications of deep learning in electron microscopy. Afterwards, we discuss hardware and software needed to get started with deep learning and interface with electron microscopes. We then review neural network components, popular architectures, and their optimization. Finally, we discuss future directions of deep learning in electron microscopy.


title: Novel Parameter Estimation and Radar Detection Approaches for Multiple Point-like Targets: Designs and Comparisons

abstract: In this work, we develop and compare two innovative strategies for parameter estimation and radar detection of multiple point-like targets. The first strategy, which appears here for the first time, jointly exploits the maximum likelihood approach and Bayesian learning to estimate targets' parameters including their positions in terms of range bins. The second strategy relies on the intuition that for high signal-to-interference plus-noise ratio values, the energy of data containing target components projected onto the nominal steering direction should be higher than the energy of data affected by interference only. The adaptivity with respect to the interference covariance matrix is also considered exploiting a training data set collected in the proximity of the window under test. Finally, another important innovation aspect concerns the adaptive estimation of the unknown number of targets by means of the model order selection rules.


title: More Embeddings, Better Sequence Labelers?

abstract: Recent work proposes a family of contextual embeddings that significantly improves the accuracy of sequence labelers over non-contextual embeddings. However, there is no definite conclusion on whether we can build better sequence labelers by combining different kinds of embeddings in various settings. In this paper, we conduct extensive experiments on 3 tasks over 18 datasets and 8 languages to study the accuracy of sequence labeling with various embedding concatenations and make three observations: (1) concatenating more embedding variants leads to better accuracy in rich-resource and cross-domain settings and some conditions of low-resource settings; (2) concatenating additional contextual sub-word embeddings with contextual character embeddings hurts the accuracy in extremely low-resource settings; (3) based on the conclusion of (1), concatenating additional similar contextual embeddings cannot lead to further improvements. We hope these conclusions can help people build stronger sequence labelers in various settings.


title: Central exclusive production of $W$ boson pairs in $pp$ collisions at the LHC in hadronic and semi-leptonic final states

abstract: We present a phenomenology study on central exclusive production of $W^+W^-$ boson pairs in proton-proton collisions at the Large Hadron Collider at 14 TeV using the forward proton detectors, such as the ATLAS Forward Proton or the CMS-TOTEM Precision Proton Spectrometer detectors. Final states where at least one of the $W$ bosons decay hadronically in a large-radius jet are considered. The latter extends previous efforts that consider solely leptonic final states. A measurement of exclusive $W^+W^-$ also allows us to further constrain anomalous quartic gauge boson interactions between photons and $W$ bosons. Expected limits on anomalous quartic gauge couplings $a_{0,C}^W$ associated to dimension-six effective operators are derived for the hadronic, semi-leptonic, and leptonic final states. It is found that the couplings can be probed down to one-dimensional values of $a_{0}^W = 3.7\times 10^{-7}$ GeV$^{-2}$ and $a_{C}^W = 9.2 \times 10^{-7}$ GeV$^{-2}$ at $95\%$ CL at an integrated luminosity of 300 fb$^{-1}$ by combining all final states, compared to values of about $a_{0}^W = 4\times 10^{-6}$ GeV$^{-2}$ and $a_{C}^W = 1\times 10^{-5}$ GeV$^{-2}$ at 95\% CL expected for the leptonic channel alone.


title: Simulation studies on regional predictive control

abstract: We apply recently proposed regional predictive control approaches to several examples to evaluate their efficiency. We present numerical results as well as hardware-in-the-loop (HIL) results.


title: Ozone-based sequential infiltration synthesis of Al2O3 nanostructures in symmetric block copolymer

abstract: Sequential infiltration synthesis (SIS) provides an original strategy to grow inorganic materials by infiltrating gaseous precursors in polymeric films. Combined with micro-phase separated nanostructures resulting from block copoly-mer (BCP) self-assembly, SIS selectively binds the precursors to only one domain mimicking the morphology of the original BCP template. This methodology represents a smart solution for the fabrication of inorganic nanostructures starting from self-assembled BCP thin films, in view of advanced lithographic application and of functional nanostructure synthesis. The SIS process using Trimethylaluminum (TMA) and H2O precursors in self-assembled PS-b-PMMA BCP thin films established as a model system, where the PMMA phase is selectively infiltrated. However the temperature range allowed by polymeric material restricts the available precursors to highly reactive reagents, such as TMA. In order to extend the SIS methodology and access a wide library of materials, a crucial step is the implementation of processes using reactive reagents that are fully compatible with the initial polymeric template. This work reports a comprehensive morphological (SEM, SE, AFM) and physico-chemical (XPS) investigation of alumina nanostructures synthesized by means of a SIS process using O3 as oxygen precursor in self-assembled PS-b-PMMA thin films with lamellar morphology. The comparison with the H2O-based SIS pro-cess validates the possibility to use O3 as oxygen precursor expanding the possible range of precursors for the fabrication of inorganic nanostructures.


title: Temporal array with superconducting nanowire single-photon detectors for photon-number-resolution

abstract: We present an experimental realization of a 16 element, temporal-array, photon-number-resolving (PNR) detector, which is a multiplexed single-photon detector that splits an input signal over multiple time-bins, and the time-bins are detected using two superconducting nanowire single-photon detectors (SNSPD). A theoretical investigation of the PNR capabilities of the detector is performed and it is concluded that compared to a single-photon detector, our array detector can resolve one order of magnitude higher mean photon numbers, given the same number of input pulses to measure. This claim is experimentally verified and we show that the detector can accurately predict photon numbers between $10^{-3}$ to $10^{2}$. Our present detector is incapable of single-shot photon-number measurements with high precision since its effective quantum efficiency is $49\,\%$. Using SNSPDs with a higher quantum efficiency the PNR performance will improve, but the photon-number resolution will still be limited by the array size.


title: Distributed Prediction-Correction ADMM for Time-Varying Convex Optimization

abstract: This paper introduces a dual-regularized ADMM approach to distributed, time-varying optimization. The proposed algorithm is designed in a prediction-correction framework, in which the computing nodes predict the future local costs based on past observations, and exploit this information to solve the time-varying problem more effectively. In order to guarantee linear convergence of the algorithm, a regularization is applied to the dual, yielding a dual-regularized ADMM. We analyze the convergence properties of the time-varying algorithm, as well as the regularization error of the dual-regularized ADMM. Numerical results show that in time-varying settings, despite the regularization error, the performance of the dual-regularized ADMM can outperform inexact gradient-based methods, as well as exact dual decomposition techniques, in terms of asymptotical error and consensus constraint violation.


title: Getting the model right; an information criterion for spectroscopy

abstract: Robust model-fitting to spectroscopic transitions is a requirement across many fields of science. The corrected Akaike and Bayesian information criteria (AICc and BIC) are most frequently used to select the optimal number of fitting parameters. In general, AICc modelling is thought to overfit (too many model parameters) and BIC underfits. For spectroscopic modelling, both AICc and BIC lack in two important respects: (a) no penalty distinction is made according to line strength such that parameters of weak lines close to the detection threshold are treated with equal importance as strong lines and (b) no account is taken of the way in which spectral lines impact on narrow data regions. In this paper we introduce a new information criterion that addresses these shortcomings, the "Spectral Information Criterion" (SpIC). Spectral simulations are used to compare performances. The main findings are (i) SpIC clearly outperforms AICc for high signal to noise data, (ii) SpIC and AICc work equally well for lower signal to noise data, although SpIC achieves this with fewer parameters, and (iii) BIC does not perform well (for this application) and should be avoided. The new method should be of broader applicability (beyond spectroscopy), wherever different model parameters influence separated small ranges within a larger dataset and/or have widely varying sensitivities.


title: Histopathology for Mohs Micrographic Surgery with Photoacoustic Remote Sensing Microscopy

abstract: Mohs micrographic surgery (MMS) is a precise oncological technique where layers of tissue are resected and examined with intraoperative histopathology to minimize the removal of normal tissue while completely excising the cancer. To achieve intraoperative pathology, the tissue is frozen, sectioned and stained over a 20- to 60-minute period, then analyzed by the MMS surgeon. Surgery is continued one layer at a time until no cancerous cells remain, meaning MMS can take several hours to complete. Ideally, it would be desirable to circumvent or augment frozen sectioning methods and directly visualize subcellular morphology on the unprocessed excised tissues. Employing photoacoustic remote sensing (PARS) microscopy, we present a non-contact label-free reflection-mode method of performing such visualizations in frozen sections of human skin. PARS leverages endogenous optical absorption contrast within cell nuclei to provide visualizations reminiscent of histochemical staining techniques. Presented here, is the first true one to one comparison between PARS microscopy and standard histopathological imaging in human tissues. We demonstrate the ability of PARS microscopy to provide large grossing scans (>1 cm2, sufficient to visualize entire MMS sections) and regional scans with subcellular lateral resolution (~300 nm).


title: A planetary system with two transiting mini-Neptunes near the radius valley transition around the bright M dwarf TOI-776

abstract: We report the discovery and characterization of two transiting planets around the bright M1 V star LP 961-53 (TOI-776, J = 8.5 mag, M = 0.54+-0.03 Msun) detected during Sector 10 observations of the Transiting Exoplanet Survey Satellite (TESS). Combining the TESS photometry with HARPS radial velocities, as well as ground-based follow-up transit observations from MEarth and LCOGT telescopes, we measured for the inner planet, TOI-776 b, a period of 8.25 d, a radius of 1.85+-0.13 Re, and a mass of 4.0+-0.9 Me; and for the outer planet, TOI-776 c, a period of 15.66 d, a radius of 2.02+-0.14 Re, and a mass of 5.3+-1.8 Me. The Doppler data shows one additional signal, with a period of 34 d, associated with the rotational period of the star. The analysis of fifteen years of ground-based photometric monitoring data and the inspection of different spectral line indicators confirm this assumption. The bulk densities of TOI-776 b and c allow for a wide range of possible interior and atmospheric compositions. However, both planets have retained a significant atmosphere, with slightly different envelope mass fractions. Thanks to their location near the radius gap for M dwarfs, we can start to explore the mechanism(s) responsible for the radius valley emergence around low-mass stars as compared to solar-like stars. While a larger sample of well-characterized planets in this parameter space is still needed to draw firm conclusions, we tentatively estimate that the stellar mass below which thermally-driven mass loss is no longer the main formation pathway for sculpting the radius valley is between 0.63 and 0.54 Msun. Due to the brightness of the star, the TOI-776 system is also an excellent target for the James Webb Space Telescope, providing a remarkable laboratory to break the degeneracy in planetary interior models and to test formation and evolution theories of small planets around low-mass stars.


title: Error protected qubits in a silicon photonic chip

abstract: General purpose quantum computers can, in principle, entangle a number of noisy physical qubits to realise composite qubits protected against errors. Architectures for measurement-based quantum computing intrinsically support error-protected qubits and are the most viable approach for constructing an all-photonic quantum computer. Here we propose and demonstrate an integrated silicon photonic architecture that both entangles multiple photons, and encodes multiple physical qubits on individual photons, to produce error-protected qubits. We realise reconfigurable graph states to compare several schemes with and without error-correction encodings and implement a range of quantum information processing tasks. We observe a success rate increase from 62.5% to 95.8% when running a phase estimation algorithm without and with error protection, respectively. Finally, we realise hypergraph states, which are a generalised class of resource states that offer protection against correlated errors. Our results show how quantum error-correction encodings can be implemented with resource-efficient photonic architectures to improve the performance of quantum algorithms.


title: Complex-Valued vs. Real-Valued Neural Networks for Classification Perspectives: An Example on Non-Circular Data

abstract: The contributions of this paper are twofold. First, we show the potential interest of Complex-Valued Neural Network (CVNN) on classification tasks for complex-valued datasets. To highlight this assertion, we investigate an example of complex-valued data in which the real and imaginary parts are statistically dependent through the property of non-circularity. In this context, the performance of fully connected feed-forward CVNNs is compared against a real-valued equivalent model. The results show that CVNN performs better for a wide variety of architectures and data structures. CVNN accuracy presents a statistically higher mean and median and lower variance than Real-Valued Neural Network (RVNN). Furthermore, if no regularization technique is used, CVNN exhibits lower overfitting. The second contribution is the release of a Python library (Barrachina 2019) using Tensorflow as back-end that enables the implementation and training of CVNNs in the hopes of motivating further research on this area.


title: Powers of binomial edge ideals with quadratic Gr\"obner bases

abstract: We study powers of binomial edge ideals associated with closed and block graphs.


title: Metastability in graded and step like variation of field and anisotropy of the Blume-Capel ferromagnet

abstract: The metastable behaviour of two dimensional anisotropic Blume-Capel ferromagnet, under the influence of graded and stepped like variations of magnetic field, has been investigated by extensive Monte Carlo simulation using Metropolis single spin flip algorithm. Starting from the initial perfectly ordered state, reversal of magnetisation has been studied in presence of the field. Metastable lifetime or the reversal time of magnetisation for the system having uniform anisotropy and in the presence of both graded and stepped field has been studied. Also the same has been explored for a {\it graded and stepped anisotropic system} in presence of a uniform field. Finite size effect is analyzed for the variation of reversal time with gradient of field ($G_h$) and a scaling relation $\langle \tau \rangle \sim L^{-\beta} f(G_h L^\alpha)$ is proposed. Spatial variation of density profile for the projection states (i.e., +1, 0 and -1) has also been studied at the moment of reversal of magnetisation. The spatial variation of the number of spin flips per site is studied. Motion of the interface (or domain wall) with the gradient of field and gradient of anisotropy are investigated and are found to follow the hyperbolic tangent like behaviours in both cases. Since, both the anisotropy and the applied field have significant impact on the metastable lifetime, an interesting competitive scenario is observed for a {\it graded anisotropic system in graded field} and a {\it stepped anisotropic system in stepped field}. A line of marginal competition has been found out for both the cases separating the regions of field dominated reversal and anisotropy dominated reversal. The decay of the metastable volume fraction was found to follow the Avrami's law. The mean reversal time was observed to decrease monotonically across the phase boundary of the ferro-para transition.


title: Dispersive optomechanics of supercavity modes in high-index disks

abstract: In this work, we study the dispersive coupling between optical quasi-bound-states in the continuum at telecom wavelengths and GHz-mechanical modes in high-index wavelength-sized disks. We show that such cavities can display values of the optomechanical coupling rate on par with optomechanical crystal cavities ($g_{0}/2\pi\backsimeq$ 800 kHz). Interestingly, optomechanical coupling of optical resonances with mechanical modes at frequencies well above 10 GHz seems attainable. We also show that mechanical leakage in the substrate can be extremely reduced by placing the disk over a thin silica pedestal. Our results suggest a new route for ultra compact optomechanical cavities which can potentially be arranged in massive arrays forming optomechanical metasurfaces.


title: Trimethylaluminum Diffusion in PMMA Thin Films during Sequential Infiltration Synthesis: In Situ Dynamic Spectroscopic Ellipsometric Investigation

abstract: Sequential infiltration synthesis (SIS) provides a successful route to grow inorganic materials into polymeric films by penetrating of gaseous precursors into the polymer, both in order to enhance the functional properties of the polymer creating an organic-inorganic hybrid material, and to fabricate inorganic nanostructures when infiltrating in patterned polymer films or in selfassembled block copolymers. A SIS process consists in a controlled sequence of metal organic precursor and co-reactant vapor exposure cycles of the polymer films in an atomic layer deposition (ALD) reactor. In this work, we present a study of the SIS process of alumina using trimethylaluminum (TMA) and H2O in various polymer films using in situ dynamic spectroscopic ellipsometry (SE). In situ dynamic SE enables time-resolved monitoring of the swelling of the polymer, which is relevant to the diffusion and retain of the metal precursor into the polymer itself. Different swelling behaviour of poly(methylmethacrylate) (PMMA) and polystyrene (PS) was observed when exposed to TMA vapor. PMMA films swell more significantly than PS films do, resulting in very different infiltrated Al2O3 thickness after polymer removal in O2 plasma. PMMA films reach different swollen states upon TMA exposure and reaction with H2O, depending on the TMA dose and on the purge duration after TMA exposure, which correspond to different amounts of metal precursor retained inside the polymer and converted to alumina. Diffusion coefficients of TMA in PMMA were extracted investigating the swelling of pristine PMMA films during TMA infiltration and shown to be dependent on polymer molecular weight. In situ dynamic SE monitoring allows to control the SIS process tuning it from an ALD-like process for long purge to a chemical vapour deposition - like process selectively confined inside the polymer films


title: Deciphering the mechanism of near-threshold $J/\psi$ photoproduction

abstract: The photoproduction of the $J/\psi$ off the proton is believed to deepen our understanding of various physics issues. On the one hand, it is proposed to provide access to the origin of the proton mass, based on the QCD multipole expansion. On the other hand, it can be employed in a study of pentaquark states. The process is usually assumed to proceed through vector-meson dominance, that is the photon couples to a $J/\psi$ which rescatters with the proton to give the $J/\psi p$ final state. In this Letter, we provide a compelling hint for and propose measurements necessary to confirm a novel production mechanism via the $\Lambda_c \bar D^{(*)}$ intermediate states. In particular, there must be cusp structures at the $\Lambda_c \bar D^{(*)}$ thresholds in the energy dependence of the $J/\psi$ photoproduction cross section. The same mechanism also implies the $J/\psi$-nucleon scattering lengths of order 1 mfm. Given this, one expects only a minor contribution of charm quarks to the nucleon mass.


title: Knowledge-Assisted Deep Reinforcement Learning in 5G Scheduler Design: From Theoretical Framework to Implementation

abstract: In this paper, we develop a knowledge-assisted deep reinforcement learning (DRL) algorithm to design wireless schedulers in the fifth-generation (5G) cellular networks with time-sensitive traffic. Since the scheduling policy is a deterministic mapping from channel and queue states to scheduling actions, it can be optimized by using deep deterministic policy gradient (DDPG). We show that a straightforward implementation of DDPG converges slowly, has a poor quality-of-service (QoS) performance, and cannot be implemented in real-world 5G systems, which are non-stationary in general. To address these issues, we propose a theoretical DRL framework, where theoretical models from wireless communications are used to formulate a Markov decision process in DRL. To reduce the convergence time and improve the QoS of each user, we design a knowledge-assisted DDPG (K-DDPG) that exploits expert knowledge of the scheduler design problem, such as the knowledge of the QoS, the target scheduling policy, and the importance of each training sample, determined by the approximation error of the value function and the number of packet losses. Furthermore, we develop an architecture for online training and inference, where K-DDPG initializes the scheduler off-line and then fine-tunes the scheduler online to handle the mismatch between off-line simulations and non-stationary real-world systems. Simulation results show that our approach reduces the convergence time of DDPG significantly and achieves better QoS than existing schedulers (reducing 30% ~ 50% packet losses). Experimental results show that with off-line initialization, our approach achieves better initial QoS than random initialization and the online fine-tuning converges in few minutes.


title: Pionless Effective Field Theory Evaluation of Nuclear Polarizability in Muonic Deuterium

abstract: We calculate the longitudinal structure function of the deuteron up through next-to-next-to-leading order in the framework of pionless effective field theory. We use these results to compute the two-photon polarizability contribution to Lamb shift in muonic deuterium, which can be utilized to extract the nuclear charge radius of the deuteron. We present analytical expressions order-by-order for the relevant transition matrix elements and the longitudinal structure function, and we give numerical results for the corresponding contributions to the Lamb shift. We also discuss the impact of relativistic and other higher-order effects. We find agreement with previous calculations and explain the accuracy of our calculation.


title: S2SD: Simultaneous Similarity-based Self-Distillation for Deep Metric Learning

abstract: Deep Metric Learning (DML) provides a crucial tool for visual similarity and zero-shot applications by learning generalizing embedding spaces, although recent work in DML has shown strong performance saturation across training objectives. However, generalization capacity is known to scale with the embedding space dimensionality. Unfortunately, high dimensional embeddings also create higher retrieval cost for downstream applications. To remedy this, we propose \emph{Simultaneous Similarity-based Self-distillation (S2SD). S2SD extends DML with knowledge distillation from auxiliary, high-dimensional embedding and feature spaces to leverage complementary context during training while retaining test-time cost and with negligible changes to the training time. Experiments and ablations across different objectives and standard benchmarks show S2SD offers notable improvements of up to 7% in Recall@1, while also setting a new state-of-the-art. Code available at https://github.com/MLforHealth/S2SD.


title: Richardson-Gaudin geminal wavefunctions in a Slater determinant basis

abstract: Geminal wavefunctions have been employed to model strongly-correlated electrons. These wavefunctions represent products of weakly-correlated pairs of electrons and reasonable approximations are computable with polynomial cost. In particular, Richardson-Gaudin states have recently been employed as a variational ansatz. This contribution serves to explain the Richardson-Gaudin wavefunctions in the conventional language of quantum chemistry.


title: nazgul: A statistical approach to gamma-ray burst localization. Triangulation via non-stationary time-series models

abstract: Context. Gamma-ray bursts can be located via arrival time signal triangulation using gamma-ray detectors in orbit throughout the solar system. The classical approach based on cross-correlations of binned light curves ignores the Poisson nature of the time-series data, and is unable to model the full complexity of the problem. Aims. To present a statistically proper and robust GRB timing/triangulation algorithm as a modern update to the original procedures used for the Interplanetary Network (IPN). Methods. A hierarchical Bayesian forward model for the unknown temporal signal evolution is learned via random Fourier features (RFF) and fitted to each detector's time-series data with time-differences that correspond to GRB's position on the sky via the appropriate Poisson likelihood. Results. Our novel method can robustly estimate the position of a GRB as verified via simulations. The uncertainties generated by the method are robust and in many cases more precise compared to the classical method. Thus, we have a method that can become a valuable tool for gravitational wave follow-up. All software and analysis scripts are made publicly available here (https://github.com/grburgess/nazgul) for the purpose of replication.


title: Intrinsic dissipative Floquet superconductors beyond mean-field theory

abstract: We study the intrinsic superconductivity in a dissipative Floquet electronic system in the presence of attractive interactions. Based on the functional Keldysh theory beyond the mean-field treatment, we find that the system shows a time-periodic bosonic condensation and reaches an intrinsic dissipative Floquet superconducting (SC) phase. Due to the interplay between dissipations and periodic modulations, the Floquet SC gap becomes "soft" and contains the diffusive fermionic modes with finite lifetimes. However, bosonic modes of the bosonic condensation are still propagating even in the presence of dissipations.


title: Regional predictive control with suboptimally extended regions of validity

abstract: Model predictive control (MPC) is based on perpetually solving optimization problems. The solution of the optimization is usually interpreted as the optimal input for the current state. However, the solution of the optimization does not just provide an optimal input, but an entire optimal affine feedback law and a polytope on which this law is optimal. We recently proposed to use this feedback law as long as the system remains in its polytope. This can be interpreted as an event-based approach, where leaving the current polytope is the event that triggers the next optimization. This approach is especially appropriate for a networked control setting since the feedback laws and their polytopes can be evaluated with a low computational effort on lean local nodes. In this article the region of validity for a feedback law is extended. More precisely, the optimal polytopes are extended to nonlinearly bounded regions of validity resulting from the intersection of a feasibility and stability region. As a result, fewer quadratic programs need to be solved compared to the optimal approach. The new validity regions are still suitable for the evaluation on a lean local node. Moreover, the regions can be adjusted for a desired closed-loop performance.


title: Sparsification Lower Bounds for List $H$-Coloring

abstract: We investigate the List $H$-Coloring problem, the generalization of graph coloring that asks whether an input graph $G$ admits a homomorphism to the undirected graph $H$ (possibly with loops), such that each vertex $v \in V(G)$ is mapped to a vertex on its list $L(v) \subseteq V(H)$. An important result by Feder, Hell, and Huang [JGT 2003] states that List $H$-Coloring is polynomial-time solvable if $H$ is a so-called bi-arc graph, and NP-complete otherwise. We investigate the NP-complete cases of the problem from the perspective of polynomial-time sparsification: can an $n$-vertex instance be efficiently reduced to an equivalent instance of bitsize $O(n^{2-\varepsilon})$ for some $\varepsilon > 0$? We prove that if $H$ is not a bi-arc graph, then List $H$-Coloring does not admit such a sparsification algorithm unless $NP \subseteq coNP/poly$. Our proofs combine techniques from kernelization lower bounds with a study of the structure of graphs $H$ which are not bi-arc graphs.


title: Feature Engineering for Data-driven Traffic State Forecast in Urban Road Networks

abstract: Most traffic state forecast algorithms when applied to urban road networks consider only the links in close proximity to the target location. However, for longer-term forecasts also the traffic state of more distant links or regions of the network are expected to provide valuable information for a data-driven algorithm. This paper studies these expectations of using a network clustering algorithm and one year of Floating Car (FCD) collected by a large fleet of vehicles. First, a clustering algorithm is applied to the data in order to extract congestion-prone regions in the Munich city network. The level of congestion inside these clusters is analyzed with the help of statistical tools. Clear spatio-temporal congestion patterns and correlations between the clustered regions are identified. These correlations are integrated into a K- Nearest Neighbors (KNN) travel time prediction algorithm. In a comparison with other approaches, this method achieves the best results. The statistical results and the performance of the KNN predictor indicate that the consideration of the network-wide traffic is a valuable feature for predictors and a promising way to develop more accurate algorithms in the future.


title: Effect of the Density of Reactive Sites in P(S-r-MMA) Film During Al2O3 Growth by Sequential Infiltration Synthesis

abstract: Sequential infiltration synthesis (SIS) consists in a controlled sequence of metal organic precursors and co-reactant vapor exposure cycles of polymer films. Two aspects characterize a SIS process: precursor molecule diffusion within the polymer matrix and precursor molecule entrapment into polymer films via chemical reaction. In this paper, we investigated SIS process for the alumina synthesis using trimethylaluminum (TMA) and H2O in thin films of poly(styrene-random-methyl methacrylate) (P(S-r-MMA)) with variable MMA content. The amount of alumina grown in the P(S-r-MMA) films linearly depends on MMA content. A relatively low concentration of MMA in the copolymer matrix is enough to guarantee the volumetric growth of alumina in the polymer film. In pure polystyrene, metal oxide seeds grow in the sub-surface region of the film. In-situ dynamic spectroscopic ellipsometry (SE) analyses provide quantitative information about TMA diffusivity in pristine P(S-r-MMA) matrices as a function of MMA fraction, allowing further insight into the process kinetics as a function of the density of reactive sites in the polymer film. This work improves the understanding of infiltration synthesis mechanism and provides a practical approach to potentially expand the library of polymers that can be effectively infiltrated by introducing reactive sites in the polymer chain.


title: Stability analysis of a signaling circuit with dual species of GTPase switches

abstract: GTPases are molecular switches that regulate a wide range of cellular processes, such as organelle biogenesis, position, shape, and signal transduction. These enzymes operate by toggling between an active ("ON") guanosine triphosphate (GTP)-bound state and an inactive ("OFF") guanosine diphosphate (GDP)-bound state; such a toggle is regulated by GEFs (guanine nucleotide exchange factors) and GAPs (GTPase activating proteins). Here we dissect a network motif between monomeric (m) and trimeric (t) GTPases assembled exclusively in eukaryotic cells of multicellular organisms. To this end, we develop a system of ordinary differential equations in which these two classes of GTPases are interlinked conditional to their ON/OFF states within a motif through feedforward and feedback loops. We provide formulas for the steady states of the system and perform local stability analysis to investigate the role of the different connections between the GTPase switches. A feedforward from the active mGTPase to the GEF of the tGTPase was sufficient to provide two locally stable states: one where both active/inactive forms of the mGTPase can be interpreted as having low concentrations and the other where both m- and tGTPase have high concentrations. When a feedback loop from the GEF of the tGTPase to the GAP of the mGTPase was added to the feedforward system, two other locally stable states emerged, both having the tGTPase inactivated and being interpreted as having low active tGTPase concentrations. Finally, the addition of a second feedback loop, from the active tGTPase to the GAP of the mGTPase, gives rise to a family of steady states parametrized by the inactive tGTPase concentrations. Our findings reveal that the coupling of these two different GTPase motifs can dramatically change their steady state behaviors and shed light on how such coupling may impact information processing in eukaryotic cells.


title: Mobility edge of Stark many-body localization

abstract: We investigate many-body localization of interacting spinless fermions in a one-dimensional disordered and tilted lattice. The fermions undergo energy-dependent transitions from ergodic to Stark many-body localization driven by the tilted potential, which are manifested by the appearance of mobility edges between delocalized states and Stark many-body localized states even when the disorder is weak. We can concretely diagnose these transitions rather than crossovers by finite-size scaling of energy-level statistics. Moreover, in the Stark many-body localization, the entanglement entropy obeys the area law scaling, in analogy to that in the conventional many-body localization.


title: Weyl's theorem for commuting tuple of paranormal and $\ast$-paranormal operators

abstract: In this article, we show that a commuting pair $T=(T_1,T_2)$ of $\ast$-paranormal operators $T_1$ and $T_2$ with quasitriangular property satisfy the Weyl's theorem-I, that is $$\sigma_T(T)\setminus\sigma_{T_W}(T)=\pi_{00}(T)$$ and a commuting pair of paranormal operators satisfy Weyl's theorem-II, that is $$\sigma_T(T)\setminus\omega(T)=\pi_{00}(T),$$ where $\sigma_T(T),\, \sigma_{T_W}(T),\,\omega(T)$ and $\pi_{00}(T)$ are the Taylor spectrum, the Taylor Weyl spectrum, the joint Weyl spectrum and the set consisting of isolated eigenvalues of $T$ with finite multiplicity, respectively. Moreover, we prove that Weyl's theorem-II holds for $f(T)$, where $T$ is a commuting pair of paranormal operators and $f$ is an analytic function in a neighbourhood of $\sigma_T(T)$.


title: Transition scenario in hypersonic axisymmetrical compression ramp flow

abstract: A high-fidelity simulation of the shock/transitional boundary layer interaction caused by a 15-degrees axisymmetrical compression ramp is performed at a freestream Mach number of 5 and a transitional Reynolds number. The inlet of the computational domain is perturbed with a white noise in order to excite convective instabilities. Coherent structures are extracted using Spectral Proper Orthogonal Decomposition (SPOD), which gives a mathematically optimal decomposition of spatio-temporally correlated structures within the flow. The mean flow is used to perform a resolvent analysis in order to study non-normal linear amplification mechanisms. The comparison between the resolvent analysis and the SPOD results provides insight on both the linear and non-linear mechanisms at play in the flow. To carry out the analysis, the flow is separated into three main regions of interest: the attached boundary layer, the mixing layer and the reattachment region. The observed transition process is dependent on the linear amplification of oblique modes in the boundary layer over a broad range of frequencies. These modes interact nonlinearly to create elongated streamwise structures which are then amplified by a linear mechanism in the rest of the domain until they break down in the reattachment region. The early nonlinear interaction is found to be essential for the transition process.


title: Improved Quantum Boosting

abstract: Boosting is a general method to convert a weak learner (which generates hypotheses that are just slightly better than random) into a strong learner (which generates hypotheses that are much better than random). Recently, Arunachalam and Maity gave the first quantum improvement for boosting, by combining Freund and Schapire's AdaBoost algorithm with a quantum algorithm for approximate counting. Their booster is faster than classical boosting as a function of the VC-dimension of the weak learner's hypothesis class, but worse as a function of the quality of the weak learner. In this paper we give a substantially faster and simpler quantum boosting algorithm, based on Servedio's SmoothBoost algorithm.


title: Formulog: Datalog for SMT-Based Static Analysis (Extended Version)

abstract: Satisfiability modulo theories (SMT) solving has become a critical part of many static analyses, including symbolic execution, refinement type checking, and model checking. We propose Formulog, a domain-specific language that makes it possible to write a range of SMT-based static analyses in a way that is both close to their formal specifications and amenable to high-level optimizations and efficient evaluation. Formulog extends the logic programming language Datalog with a first-order functional language and mechanisms for representing and reasoning about SMT formulas; a novel type system supports the construction of expressive formulas, while ensuring that neither normal evaluation nor SMT solving goes wrong. Our case studies demonstrate that a range of SMT-based analyses can naturally and concisely be encoded in Formulog, and that -- thanks to this encoding -- high-level Datalog-style optimizations can be automatically and advantageously applied to these analyses.


title: Dynamics of delayed neural field models in two-dimensional spatial domains

abstract: Delayed neural field models can be viewed as a dynamical system in an appropriate functional analytic setting. On two dimensional rectangular space domains, and for a special class of connectivity and delay functions, we describe the spectral properties of the linearized equation. We transform the characteristic integral equation for the delay differential equation (DDE) into a linear partial differential equation (PDE) with boundary conditions. We demonstrate that finding eigenvalues and eigenvectors of the DDE is equivalent with obtaining nontrivial solutions of this boundary value problem (BVP). When the connectivity kernel consists of a single exponential, we construct a basis of the solutions of this BVP that forms a complete set in $L^2$. This gives a complete characterization of the spectrum and is used to construct a solution to the resolvent problem. As an application we give an example of a Hopf bifurcation and compute the first Lyapunov coefficient.


title: Functional data analysis: An application to COVID-19 data in the United States

abstract: The COVID-19 pandemic so far has caused huge negative impacts on different areas all over the world, and the United States (US) is one of the most affected countries. In this paper, we use methods from the functional data analysis to look into the COVID-19 data in the US. We explore the modes of variation of the data through a functional principal component analysis (FPCA), and study the canonical correlation between confirmed and death cases. In addition, we run a cluster analysis at the state level so as to investigate the relation between geographical locations and the clustering structure. Lastly, we consider a functional time series model fitted to the cumulative confirmed cases in the US, and make forecasts based on the dynamic FPCA. Both point and interval forecasts are provided, and the methods for assessing the accuracy of the forecasts are also included.


title: The Bi-Laplacian with Wentzell boundary conditions on Lipschitz domains

abstract: We investigate the Bi-Laplacian with Wentzell boundary conditions in a bounded domain $\Omega\subseteq\mathbb{R}^d$ with Lipschitz boundary $\Gamma$. More precisely, using form methods, we show that the associated operator on the ground space $L^2(\Omega)\times L^2(\Gamma)$ has compact resolvent and generates a holomorphic and strongly continuous real semigroup of self-adjoint operators. Furthermore, we give a full characterization of the domain in terms of Sobolev spaces, also proving H\"older regularity of solutions, allowing classical interpretation of the boundary condition. Finally, we investigate spectrum and asymptotic behavior of the semigroup, as well as eventual positivity.


title: Exact saturation in pseudo-elementary classes for simple and stable theories

abstract: We study PC-exact saturation for stable and simple theories. Among other results, we show that PC-exact saturation characterizes the stability cardinals of size at least continuum of a countable stable theory and, additionally, that simple unstable theories have PC-exact saturation at singular cardinals, satisfying mild set-theoretic hypotheses, which had previously been open even for the random graph. We characterize supersimplicity of countable theories in terms of having PC-exact saturation at singular cardinals of countable cofinality. We also consider the local analogue of PC-exact saturation, showing that local PC-exact saturation for singular cardinals of countable cofinality characterizes supershort theories.


title: GraphCodeBERT: Pre-training Code Representations with Data Flow

abstract: Pre-trained models for programming language have achieved dramatic empirical improvements on a variety of code-related tasks such as code search, code completion, code summarization, etc. However, existing pre-trained models regard a code snippet as a sequence of tokens, while ignoring the inherent structure of code, which provides crucial code semantics and would enhance the code understanding process. We present GraphCodeBERT, a pre-trained model for programming language that considers the inherent structure of code. Instead of taking syntactic-level structure of code like abstract syntax tree (AST), we use data flow in the pre-training stage, which is a semantic-level structure of code that encodes the relation of "where-the-value-comes-from" between variables. Such a semantic-level structure is neat and does not bring an unnecessarily deep hierarchy of AST, the property of which makes the model more efficient. We develop GraphCodeBERT based on Transformer. In addition to using the task of masked language modeling, we introduce two structure-aware pre-training tasks. One is to predict code structure edges, and the other is to align representations between source code and code structure. We implement the model in an efficient way with a graph-guided masked attention function to incorporate the code structure. We evaluate our model on four tasks, including code search, clone detection, code translation, and code refinement. Results show that code structure and newly introduced pre-training tasks can improve GraphCodeBERT and achieves state-of-the-art performance on the four downstream tasks. We further show that the model prefers structure-level attentions over token-level attentions in the task of code search.


title: On the EIT problem for nonorientable surfaces

abstract: Let $(\Omega,g)$ be a smooth compact two-dimensional Riemannian manifold with boundary, $\Lambda_g: f\mapsto \partial_\nu u|_{\partial\Omega}$ its DN map, where $u$ obeys $\Delta_g u=0$ in $\Omega$ and $u|_{\partial \Omega}=f$. The Electric Impedance Tomography problem is to determine $\Omega$ from $\Lambda_g$. A criterion is proposed that enables one to detect (via $\Lambda_g$) whether $\Omega$ is orientable or not. The algebraic version of the BC-method is applied to solve the EIT problem for the Moebius band. The main instrument is the algebra of holomorphic functions on the double covering ${\mathbb M}$ of $M$, which is determined by $\Lambda_g$ up to an isometric isomorphism. Its Gelfand spectrum (the set of characters) plays the role of the material for constructing a relevant copy $(M',g')$ of $(M,g)$. This copy is conformally equivalent to the original, provides $\partial M'=\partial M,\,\,\Lambda_{g'}=\Lambda_g$, and thus solves the problem.


title: A new front-tracking Lagrangian model for the modeling of dynamic and post-dynamic recrystallization

abstract: A new method for the simulation of evolving multi-domains problems has been introduced in previous works (RealIMotion), Florez et al. (2020) and further developed in parallel in the context of isotropic Grain Growth (GG) with no consideration for the effects of the Stored Energy (SE) due to dislocations. The methodology consists in a new front-tracking approach where one of the originality is that not only interfaces between grains are discretized but their bulks are also meshed and topological changes of the domains are driven by selective local remeshing operations performed on the Finite Element (FE) mesh. In this article, further developments and studies of the model will be presented, mainly on the development of a model taking into account grain boundary migration by (GBM) SE. Further developments for the nucleation of new grains will be presented, allowing to model Dynamic Recrystallization (DRX) and Post-Dynamic Recrystallization (PDRX) phenomena. The accuracy and the performance of the numerical algorithms have been proven to be very promising in Florez et al. (2020). Here the results for multiple test cases will be given in order to validate the accuracy of the model taking into account GG and SE. The computational performance will be evaluated for the DRX and PDRX mechanisms and compared to a classical Finite Element (FE) framework using a Level-Set (LS) formulation.


title: Face Mask Detection using Transfer Learning of InceptionV3

abstract: The world is facing a huge health crisis due to the rapid transmission of coronavirus (COVID-19). Several guidelines were issued by the World Health Organization (WHO) for protection against the spread of coronavirus. According to WHO, the most effective preventive measure against COVID-19 is wearing a mask in public places and crowded areas. It is very difficult to monitor people manually in these areas. In this paper, a transfer learning model is proposed to automate the process of identifying the people who are not wearing mask. The proposed model is built by fine-tuning the pre-trained state-of-the-art deep learning model, InceptionV3. The proposed model is trained and tested on the Simulated Masked Face Dataset (SMFD). Image augmentation technique is adopted to address the limited availability of data for better training and testing of the model. The model outperformed the other recently proposed approaches by achieving an accuracy of 99.9% during training and 100% during testing.


title: Self-attenuation of extreme events in Navier-Stokes turbulence

abstract: Turbulent fluid flows are ubiquitous in nature and technology, and are mathematically described by the incompressible Navier-Stokes equations (INSE). A hallmark of turbulence is spontaneous generation of intense whirls, resulting from amplification of the fluid rotation-rate (vorticity) by its deformation-rate (strain). This interaction, encoded in the non-linearity of INSE, is non-local, i.e., depends on the entire state of the flow, constituting a serious hindrance in turbulence theory and in establishing regularity of INSE. Here, we unveil a novel aspect of this interaction, by separating strain into local and non-local contributions utilizing the Biot-Savart integral of vorticity in a sphere of radius R. Analyzing highly-resolved numerical turbulent solutions to INSE, we find that when vorticity becomes very large, the local strain over small R surprisingly counteracts further amplification. This uncovered self-attenuation mechanism is further shown to be connected to local Beltramization of the flow, and could provide a direction in establishing the regularity of INSE.

categories: physics.flu-dyn cs.NA math.NA physics.comp-ph


title: Microtubule Tracking in Electron Microscopy Volumes

abstract: We present a method for microtubule tracking in electron microscopy volumes. Our method first identifies a sparse set of voxels that likely belong to microtubules. Similar to prior work, we then enumerate potential edges between these voxels, which we represent in a candidate graph. Tracks of microtubules are found by selecting nodes and edges in the candidate graph by solving a constrained optimization problem incorporating biological priors on microtubule structure. For this, we present a novel integer linear programming formulation, which results in speed-ups of three orders of magnitude and an increase of 53% in accuracy compared to prior art (evaluated on three 1.2 x 4 x 4$\mu$m volumes of Drosophila neural tissue). We also propose a scheme to solve the optimization problem in a block-wise fashion, which allows distributed tracking and is necessary to process very large electron microscopy volumes. Finally, we release a benchmark dataset for microtubule tracking, here used for training, testing and validation, consisting of eight 30 x 1000 x 1000 voxel blocks (1.2 x 4 x 4$\mu$m) of densely annotated microtubules in the CREMI data set (https://github.com/nilsec/micron).


title: A Principle of Least Action for the Training of Neural Networks

abstract: Neural networks have been achieving high generalization performance on many tasks despite being highly over-parameterized. Since classical statistical learning theory struggles to explain this behavior, much effort has recently been focused on uncovering the mechanisms behind it, in the hope of developing a more adequate theoretical framework and having a better control over the trained models. In this work, we adopt an alternate perspective, viewing the neural network as a dynamical system displacing input particles over time. We conduct a series of experiments and, by analyzing the network's behavior through its displacements, we show the presence of a low kinetic energy displacement bias in the transport map of the network, and link this bias with generalization performance. From this observation, we reformulate the learning problem as follows: finding neural networks which solve the task while transporting the data as efficiently as possible. This offers a novel formulation of the learning problem which allows us to provide regularity results for the solution network, based on Optimal Transport theory. From a practical viewpoint, this allows us to propose a new learning algorithm, which automatically adapts to the complexity of the given task, and leads to networks with a high generalization ability even in low data regimes.


title: Modeling human visual search: A combined Bayesian searcher and saliency map approach for eye movement guidance in natural scenes

abstract: Finding objects is essential for almost any daily-life visual task. Saliency models have been useful to predict fixation locations in natural images, but are static, i.e., they provide no information about the time-sequence of fixations. Nowadays, one of the biggest challenges in the field is to go beyond saliency maps to predict a sequence of fixations related to a visual task, such as searching for a given target. Bayesian observer models have been proposed for this task, as they represent visual search as an active sampling process. Nevertheless, they were mostly evaluated on artificial images, and how they adapt to natural images remains largely unexplored. Here, we propose a unified Bayesian model for visual search guided by saliency maps as prior information. We validated our model with a visual search experiment in natural scenes recording eye movements. We show that, although state-of-the-art saliency models perform well in predicting the first two fixations in a visual search task, their performance degrades to chance afterward. This suggests that saliency maps alone are good to model bottom-up first impressions, but are not enough to explain the scanpaths when top-down task information is critical. Thus, we propose to use them as priors of Bayesian searchers. This approach leads to a behavior very similar to humans for the whole scanpath, both in the percentage of target found as a function of the fixation rank and the scanpath similarity, reproducing the entire sequence of eye movements.


title: A Glimpse of the First Eight Months of the COVID-19 Literature on Microsoft Academic Graph: Themes, Citation Contexts, and Uncertainties

abstract: As scientists worldwide search for answers to the overwhelmingly unknown behind the deadly pandemic, the literature concerning COVID-19 has been growing exponentially. Keeping abreast of the body of literature at such a rapidly advancing pace poses significant challenges not only to active researchers but also to the society as a whole. Although numerous data resources have been made openly available, the analytic and synthetic process that is essential in effectively navigating through the vast amount of information with heightened levels of uncertainty remains a significant bottleneck. We introduce a generic method that facilitates the data collection and sense-making process when dealing with a rapidly growing landscape of a research domain such as COVID-19 at multiple levels of granularity. The method integrates the analysis of structural and temporal patterns in scholarly publications with the delineation of thematic concentrations and the types of uncertainties that may offer additional insights into the complexity of the unknown. We demonstrate the application of the method in a study of the COVID-19 literature.


title: (Mis-)Matching Type-B Anomalies on the Higgs Branch

abstract: Building on arXiv:1911.05827, we uncover new properties of type-B conformal anomalies for Coulomb-branch operators in continuous families of 4D $\mathcal{N}=2$ SCFTs. We study a large class of such anomalies on the Higgs branch, where conformal symmetry is spontaneously broken, and compare them with their counterpart in the CFT phase. In Lagrangian theories, the non-perturbative matching of the anomalies can be determined with a weak coupling Feynman diagram computation involving massive multi-loop banana integrals. We extract the part corresponding to the anomalies of interest. Our calculations support the general conjecture that the Coulomb-branch type-B conformal anomalies always match on the Higgs branch when the IR Coulomb-branch chiral ring is empty. In the opposite case, there are anomalies that do not match. An intriguing implication of the mismatch is the existence of a second covariantly constant metric on the conformal manifold (other than the Zamolodchikov metric), which imposes previously unknown restrictions on its holonomy group.


title: Exploring the di-$J/\psi$ resonances based on {\it{ab initio}} perturbative QCD

abstract: We extend the nonrelativistic QCD framework to explore the nature of the newly discovered di-$J/\psi$ resonances. Assuming them as either molecule-like states or tetraquarks, we calculated their hadroproduction cross sections at the LHC. We find that the observed resonances are most likely spin-2 particles, and there should exist their spin-0 counterparts near these resonances. The ratio of production cross sections of the observed resonances to the latent spin-0 ones are also presented, which might help to distinguish molecule-like states from tetraquarks.


title: A factorization theorem for harmonic maps

abstract: Let $f$ be a harmonic map from a Riemann surface to a Riemannian $n$-manifold. We prove that if there is a holomorphic diffeomorphism $h$ between open subsets of the surface such that $f\circ h = f$, then $f$ factors through a holomorphic map onto another Riemann surface. If such $h$ is anti-holomorphic, we obtain an analogous statement. For minimal maps, this result is well known and is a consequence of the theory of branched immersions of surfaces due to Gulliver-Osserman-Royden. Our proof relies on various geometric properties of the Hopf differential.


title: Event-Based Backpropagation can compute Exact Gradients for Spiking Neural Networks

abstract: Spiking neural networks combine analog computation with event-based communication using discrete spikes. While the impressive advances of deep learning are enabled by training non-spiking artificial neural networks using the backpropagation algorithm, applying this algorithm to spiking networks was previously hindered by the existence of discrete spike events and discontinuities. For the first time, this work derives the backpropagation algorithm for a continuous-time spiking neural network and a general loss function by applying the adjoint method together with the proper partial derivative jumps, allowing for backpropagation through discrete spike events without approximations. This algorithm, EventProp, backpropagates errors at spike times in order to compute the exact gradient in an event-based, temporally and spatially sparse fashion. We use gradients computed via EventProp to train networks on the Yin-Yang and MNIST datasets using either a spike time or voltage based loss function and report competitive performance. Our work supports the rigorous study of gradient-based learning algorithms in spiking neural networks and provides insights toward their implementation in novel brain-inspired hardware.


title: Lifshitz phase transitions in one-dimensional Gamma model

abstract: In this paper, we study quantum phase transitions and magnetic properties of a one-dimensional spin-1/2 Gamma model, which describes the off-diagonal exchange interactions between edge-shared octahedra with strong spin-orbit couplings along the sawtooth chain. The competing exchange interactions between the nearest neighbors and the second neighbors stabilize semimetallic ground state in terms of spinless fermions, and give rise to a rich phase diagram, which consists of three gapless phases. We find distinct phases are characterized by the number of Weyl nodes in the momentum space, and such changes in the topology of the Fermi surface without symmetry breaking produce a variety of Lifshitz transitions, in which the Weyl nodes situating at $k=\pi$ interchange from type I to type II. A coexistence of type-I and type-II Weyl nodes is found in phase II. The information measures including concurrence, entanglement entropy and relative entropy can effectively signal the second-order transitions. The results indicate that the Gamma model can act as an exactly solvable model to describe Lifshitz phase transitions in correlated electron systems.


title: Evaluating Interactive Summarization: an Expansion-Based Framework

abstract: Allowing users to interact with multi-document summarizers is a promising direction towards improving and customizing summary results. Different ideas for interactive summarization have been proposed in previous work but these solutions are highly divergent and incomparable. In this paper, we develop an end-to-end evaluation framework for expansion-based interactive summarization, which considers the accumulating information along an interactive session. Our framework includes a procedure of collecting real user sessions and evaluation measures relying on standards, but adapted to reflect interaction. All of our solutions are intended to be released publicly as a benchmark, allowing comparison of future developments in interactive summarization. We demonstrate the use of our framework by evaluating and comparing baseline implementations that we developed for this purpose, which will serve as part of our benchmark. Our extensive experimentation and analysis of these systems motivate our design choices and support the viability of our framework.


title: Demultiplexing infrasound phonons with tunable magnetic lattices

abstract: Controlling infrasound signals is crucial to many processes ranging from predicting atmospheric events and seismic activities to sensing nuclear detonations. These waves can be manipulated through phononic crystals and acoustic metamaterials. However, at such ultra-low frequencies, the size (usually on the order of meters) and the mass (usually on the order of many kilograms) of these materials can hinder its potential applications in the infrasonic domain. Here, we utilize tunable lattices of repelling magnets to guide and sort infrasound waves into different channels based on their frequencies. We construct our lattices by confining meta-atoms (free-floating macroscopic disks with embedded magnets) within a magnetic boundary. By changing the confining boundary, we control the meta-atoms' spacing and therefore the intensity of their coupling potentials and wave propagation characteristics. As a demonstration of principle, we present the first experimental realization of an infrasound phonon demultiplexer (i.e., guiding ultra-low frequency waves into different channels based on their frequencies). The realized platform can be utilized to manipulate ultra-low frequency waves, within a relatively small volume, while utilizing negligible mass. In addition, the self-assembly nature of the meta-atoms can be key in creating re-programmable materials with exceptional nonlinear properties.


title: Reconfigurable Artificial Microswimmers with Internal Feedback

abstract: Micron-size self-propelling particles are often proposed as synthetic models for biological microswimmers, yet they lack internally regulated adaptation, which is central to the autonomy of their biological counterparts. Conversely, adaptation and autonomy can be encoded in larger-scale soft-robotic devices, but transferring these capabilities to the colloidal scale remains elusive. Here, we create a new class of responsive microswimmers, powered by induced-charge electrophoresis, which can adapt their motility to external stimuli via an internal feedback. Using sequential capillary assembly, we fabricate deterministic colloidal clusters comprising soft thermoresponsive microparticles, which, upon spontaneous reconfiguration, induce motility changes, such as adaptation of the clusters' propulsion velocity and reversal of its direction. We rationalize the response in terms of a coupling between self-propulsion and variations of particle shape and dielectric properties. Harnessing those allows for strategies to achieve local dynamical control with simple illumination patterns, revealing exciting opportunities for the development of new tactic active materials.


title: The EDGE-CALIFA survey: exploring the role of the molecular gas on the galaxy star formation quenching

abstract: Understanding how galaxies cease to form stars represents an outstanding challenge for galaxy evolution theories. This process of "star formation quenching" has been related to various causes, including Active Galactic Nuclei (AGN) activity, the influence of large-scale dynamics, and the environment in which galaxies live. In this paper, we present the first results from a follow-up of CALIFA survey galaxies with observations of molecular gas obtained with the APEX telescope. Together with EDGE survey CARMA observations, we collect $^{12}$CO observations that cover approximately one effective radius in 472 CALIFA galaxies. We observe that the deficit of galaxy star formation with respect to the star formation main sequence (SFMS) increases with the absence of molecular gas and with a reduced efficiency of conversion of molecular gas into stars, in line with results of other integrated studies. However, by dividing the sample into galaxies dominated by star formation and galaxies quenched in their centres (as indicated by the average value of the H$\alpha$ equivalent width), we find that this deficit increases sharply once a certain level of gas consumption is reached, indicating that different mechanisms drive separation from the SFMS in star-forming and quenched galaxies. Our results indicate that differences in the amount of molecular gas at a fixed stellar mass are the primary driver for the dispersion in the SFMS, and the most likely explanation for the start of star-formation quenching. However, once a galaxy is quenched, changes in star formation efficiency drive how much a retired galaxy separates in star formation rate from star-forming ones of similar masses. In other words, once a paucity of molecular gas has significantly reduced star formation, changes in the star formation efficiency are what drives a galaxy deeper into the red cloud, retiring it.


title: Sharp rigidity estimates for incompatible fields as consequence of the Bourgain Brezis div-curl result

abstract: In this note we show that a sharp rigidity estimate and a sharp Korn's inequality for matrix-valued fields whose incompatibility is a bounded measure can be obtained as a consequence of a Hodge decomposition with critical integrability due to Bourgain and Brezis.


title: Computational models in Electroencephalography

abstract: Computational models lie at the intersection of basic neuroscience and healthcare applications because they allow researchers to test hypotheses \textit{in silico} and predict the outcome of experiments and interactions that are very hard to test in reality. Yet, what is meant by "computational model" is understood in many different ways by researchers in different fields of neuroscience and psychology, hindering communication and collaboration. In this review, we point out the state of the art of computational modeling in Electroencephalography (EEG) and outline how these models can be used to integrate findings from electrophysiology, network-level models, and behavior. On the one hand, computational models serve to investigate the mechanisms that generate brain activity, for example measured with EEG, such as the transient emergence of oscillations at different frequency bands and/or with different spatial topographies. On the other hand, computational models serve to design experiments and test hypotheses \emph{in silico}. The final purpose of computational models of EEG is to obtain a comprehensive understanding of the mechanisms that underlie the EEG signal. This is crucial for an accurate interpretation of EEG measurements that may ultimately serve in the development of novel clinical applications.


title: Photonic and plasmonic transition radiation from graphene

abstract: In the framework of full Maxwell equations, we systematically study the electromagnetic radiation, namely the transition radiation, when a swift electron perpendicularly crosses a monolayer graphene. Based on the plane wave expansion and the Sommerfeld integration, we demonstrate the spatial distribution of this free-electron radiation process in the frequency domain, which clearly shows the broadband excitation of both the photons and graphene plasmons. Moreover, the radiation spectra for the excited photons and graphene plasmons are analytically derived. We find that the excitation of photons and graphene plasmons favours different velocities of the swift electron. To be specific, a higher electron velocity can give rise to the excitation of photons with better directivity and higher intensity, while a lower electron velocity can enable the efficient excitation of graphene plasmons in a broader frequency range. Our work indicates that the interaction between swift charged particles and ultrathin 2D materials or metasurfaces is promising for the design of terahertz on-chip radiation sources.


title: Towards Stable Imbalanced Data Classification via Virtual Big Data Projection

abstract: Virtual Big Data (VBD) proved to be effective to alleviate mode collapse and vanishing generator gradient as two major problems of Generative Adversarial Neural Networks (GANs) very recently. In this paper, we investigate the capability of VBD to address two other major challenges in Machine Learning including deep autoencoder training and imbalanced data classification. First, we prove that, VBD can significantly decrease the validation loss of autoencoders via providing them a huge diversified training data which is the key to reach better generalization to minimize the over-fitting problem. Second, we use the VBD to propose the first projection-based method called cross-concatenation to balance the skewed class distributions without over-sampling. We prove that, cross-concatenation can solve uncertainty problem of data driven methods for imbalanced classification.


title: Transfer Graph Neural Networks for Pandemic Forecasting

abstract: The recent outbreak of COVID-19 has affected millions of individuals around the world and has posed a significant challenge to global healthcare. From the early days of the pandemic, it became clear that it is highly contagious and that human mobility contributes significantly to its spread. In this paper, we study the impact of population movement on the spread of COVID-19, and we capitalize on recent advances in the field of representation learning on graphs to capture the underlying dynamics. Specifically, we create a graph where nodes correspond to a country's regions and the edge weights denote human mobility from one region to another. Then, we employ graph neural networks to predict the number of future cases, encoding the underlying diffusion patterns that govern the spread into our learning model. Furthermore, to account for the limited amount of training data, we capitalize on the pandemic's asynchronous outbreaks across countries and use a model-agnostic meta-learning based method to transfer knowledge from one country's model to another's. We compare the proposed approach against simple baselines and more traditional forecasting techniques in 3 European countries. Experimental results demonstrate the superiority of our method, highlighting the usefulness of GNNs in epidemiological prediction. Transfer learning provides the best model, highlighting its potential to improve the accuracy of the predictions in case of secondary waves, if data from past/parallel outbreaks is utilized.


title: Conformal welding of quantum disks

abstract: Two-pointed quantum disks with a weight parameter $W > 0$ are a family of finite-area random surfaces that arise naturally in Liouville quantum gravity. In this paper we show that conformally welding two quantum disks according to their boundary lengths gives another quantum disk decorated with an independent chordal $\mathrm{SLE}_\kappa(\rho_-;\rho_+)$ curve. This is the finite-volume counterpart of the classical result of Sheffield (2010) and Duplantier-Miller-Sheffield (2014) on the welding of infinite-area two-pointed quantum surfaces called quantum wedges, which is fundamental to the mating-of-trees theory. Our results can be used to give unified proofs of the mating-of-trees theorems for the quantum disk and the quantum sphere, in addition to a mating-of-trees description of the weight $W = \frac{\gamma^2}{2}$ quantum disk. Moreover, it serves as a key ingredient in our companion work [AHS21], which proves an exact formula for $\mathrm{SLE}_\kappa(\rho_-;\rho_+)$ using conformal welding of random surfaces and a conformal welding result giving the so-called SLE loop.


title: Detecci\'on de comunidades en redes: Algoritmos y aplicaciones

abstract: This master's thesis work has the objective of performing an analysis of the methods for detecting communities in networks. As an initial part, I study of the main features of graph theory and communities, as well as common measures in this problem. Subsequently, I was performed a review of the main methods of detecting communities, developing a classification, taking into account its characteristics and computational complexity for the detection of strengths and weaknesses in the methods, as well as later works. Then, study the problem of classification of a clustering method, this in order to evaluate the quality of the communities detected by analyzing different measures. Finally conclusions are elaborated and possible lines of work that can be derived.


title: The variance of relative surprisal as single-shot quantifier

abstract: The variance of (relative) surprisal, also known as varentropy, so far mostly plays a role in information theory as quantifying the leading order corrections to asymptotic i.i.d.~limits. Here, we comprehensively study the use of it to derive single-shot results in (quantum) information theory. We show that it gives genuine sufficient and necessary conditions for approximate state-transitions between pairs of quantum states in the single-shot setting, without the need for further optimization. We also clarify its relation to smoothed min- and max-entropies, and construct a monotone for resource theories using only the standard (relative) entropy and variance of (relative) surprisal. This immediately gives rise to enhanced lower bounds for entropy production in random processes. We establish certain properties of the variance of relative surprisal which will be useful for further investigations, such as uniform continuity and upper bounds on the violation of sub-additivity. Motivated by our results, we further derive a simple and physically appealing axiomatic single-shot characterization of (relative) entropy which we believe to be of independent interest. We illustrate our results with several applications, ranging from interconvertibility of ergodic states, over Landauer erasure to a bound on the necessary dimension of the catalyst for catalytic state transitions and Boltzmann's H-theorem.


title: Impact and dynamics of hate and counter speech online

abstract: Citizen-generated counter speech is a promising way to fight hate speech and promote peaceful, non-polarized discourse. However, there is a lack of large-scale longitudinal studies of its effectiveness for reducing hate speech. To this end, we perform an exploratory analysis of the effectiveness of counter speech using several different macro- and micro-level measures to analyze 180,000 political conversations that took place on German Twitter over four years. We report on the dynamic interactions of hate and counter speech over time and provide insights into whether, as in `classic' bullying situations, organized efforts are more effective than independent individuals in steering online discourse. Taken together, our results build a multifaceted picture of the dynamics of hate and counter speech online. While we make no causal claims due to the complexity of discourse dynamics, our findings suggest that organized hate speech is associated with changes in public discourse and that counter speech -- especially when organized -- may help curb hateful rhetoric in online discourse.


title: Non-equilibrium Properties of an Active Nanoparticle in a Harmonic Potential

abstract: Active particles break out of thermodynamic equilibrium thanks to their directed motion, which leads to complex and interesting behaviors in the presence of confining potentials. When dealing with active nanoparticles, however, the overwhelming presence of rotational diffusion hinders directed motion, leading to an increase of their effective temperature, but otherwise masking the effects of self-propulsion. Here, we demonstrate an experimental system where an active nanoparticle immersed in a critical solution and held in an optical harmonic potential features far-from-equilibrium behavior beyond an increase of its effective temperature. When increasing the laser power, we observe a cross-over from a Boltzmann distribution to a non-equilibrium state, where the particle performs fast orbital rotations about the beam axis. These findings are rationalized by solving the Fokker-Planck equation for the particle's position and orientation in terms of a moment expansion. The proposed self-propulsion mechanism results from the particle's non-sphericity and the lower critical point of the solute.


title: Authoritarianism vs. democracy: Simulating responses to disease outbreaks

abstract: Disease outbreaks force the governments to rapid decisions to deal with. However, the rapid stream of decision-making could be costly in terms of the democratic representativeness. The aim of the paper is to investigate the trade-off between pluralism of preferences and the time required to approach a decision. To this aim we develop and test a modified version of the Hegselmann and Krause (2002) model to capture these two characteristics of the decisional process in different institutional contexts. Using a twofold geometrical institutional setting, we simulate the impact of disease outbreaks to check whether countries exhibits idiosyncratic effects, depending on their institutional frameworks. Main findings show that the degree of pluralism in political decisions is not necessarily associated with worse performances in managing emergencies, provided that the political debate is mature enough.


title: A Multimodal Memes Classification: A Survey and Open Research Issues

abstract: Memes are graphics and text overlapped so that together they present concepts that become dubious if one of them is absent. It is spread mostly on social media platforms, in the form of jokes, sarcasm, motivating, etc. After the success of BERT in Natural Language Processing (NLP), researchers inclined to Visual-Linguistic (VL) multimodal problems like memes classification, image captioning, Visual Question Answering (VQA), and many more. Unfortunately, many memes get uploaded each day on social media platforms that need automatic censoring to curb misinformation and hate. Recently, this issue has attracted the attention of researchers and practitioners. State-of-the-art methods that performed significantly on other VL dataset, tends to fail on memes classification. In this context, this work aims to conduct a comprehensive study on memes classification, generally on the VL multimodal problems and cutting edge solutions. We propose a generalized framework for VL problems. We cover the early and next-generation works on VL problems. Finally, we identify and articulate several open research issues and challenges. This is the first study that presents the generalized view of the advanced classification techniques concerning memes classification to the best of our knowledge. We believe this study presents a clear road-map for the Machine Learning (ML) research community to implement and enhance memes classification techniques.


title: Feasibility Analysis and Preliminary Design of ChipSat Entry for In-situ Investigation of the Atmosphere of Venus

abstract: Recent miniaturization of electronics in very small, low-cost and low-power configurations suitable for use in spacecraft have inspired innovative small-scale satellite concepts, such as ChipSats, centimeter-scale satellites with a mass of a few grams. These extremely small spacecraft have the potential to usher in a new age of space science accessibility. Due to their low ballistic coefficient, ChipSats can potentially be used in a swarm constellation for extended surveys of planetary atmospheres, providing large amounts of data with high reliability and redundancy. We present a preliminary feasibility analysis of a ChipSat planetary atmospheric entry mission with the purpose of searching for traces of microscopic lifeforms in the atmosphere of Venus. Indeed, the lower cloud layer of the Venusian atmosphere could be a good target for searching for microbial lifeforms, due to the favourable atmospheric conditions and the presence of micron-sized sulfuric acid aerosols. A numerical model simulating the planetary entry of a spacecraft of specified geometry, applicable to any atmosphere for which sufficient atmospheric data are available, is implemented and verified. The results are used to create a high-level design of a ChipSat mission cruising in the Venusian atmosphere at altitudes favorable for the existence of life. The paper discusses the ChipSat mission concept and considerations about the spacecraft preliminary design at system level, including the selection of a potential payload.


title: High-energy gamma-ray study of the dynamically young SNR G150.3+4.5

abstract: The supernova remnant (SNR) G150.3+4.5 was recently discovered in the radio band; it exhibits a shell-like morphology with an angular size of $\sim 3^{\circ}$, suggesting either an old or a nearby SNR. Extended $\gamma$-ray emission spatially coincident with the SNR was reported in the Fermi Galactic Extended Source Catalog, with a power-law spectral index of $\Gamma$ = 1.91 $\pm$ 0.09. Studying particle acceleration in SNRs through their $\gamma$-ray emission is of primary concern to assess the nature of accelerated particles and the maximum energy they can reach. Using more than ten years of Fermi-LAT data, we investigate the morphological and spectral properties of the SNR G150.3+4.5 from 300 MeV to 3 TeV. We use the latest releases of the Fermi-LAT catalog, the instrument response functions and the Galactic and isotropic diffuse emissions. We use ROSAT all-sky survey data to assess any thermal and nonthermal X-ray emission, and we derive minimum and maximum distance to G150.3+4.5. We describe the $\gamma$-ray emission of G150.3+4.5 by an extended component which is found to be spatially coincident with the radio SNR. The spectrum is hard and the detection of photons up to hundreds of GeV points towards an emission from a dynamically young SNR. The lack of X-ray emission gives a tight constraint on the ambient density $n_0 \leq 3.6 \times 10^{-3}$ cm$^{-3}$. Since G150.3+4.5 is not reported as a historical SNR, we impose a lower limit on its age of $t$ = 1 kyr. We estimate its distance to be between 0.7 and 4.5 kpc. We find that G150.3+4.5 is spectrally similar to other dynamically young and shell-type SNRs, such as RX J1713.7$-$3946 or Vela Junior. The broadband nonthermal emission is explained with a leptonic scenario, implying a downstream magnetic field of $B = 5$ $\mu$G and acceleration of particles up to few TeV energies.


title: Restored strange metal phase through suppression of charge density waves in underdoped YBa$_2$Cu$_3$O$_{7-\delta}$

abstract: The normal state of optimally doped cuprates is dominated by the "strange metal" phase that shows a linear temperature ($T$) dependence of the resistivity persisting down to the lowest $T$. For underdoped cuprates this behavior is lost below the pseudogap temperature $T^*$, where Charge Density Waves (CDW) together with other intertwined local orders characterize the ground state. Here we show that the $T$-linear resistivity of highly strained, ultrathin and underdoped YBa$_2$Cu$_3$O$_{7-\delta}$ films is restored when the CDW amplitude, detected by Resonant Inelastic X-ray scattering, is suppressed. This observation points towards an intimate connection between the onset of CDW and the departure from $T$-linear resistivity in underdoped cuprates. Our results illustrate the potential of using strain control to manipulate the ground state of quantum materials.


title: A sharp upper bound for the $2$-torsion of class groups of multiquadratic fields

abstract: Let $K$ be a multiquadratic extension of $\mathbb{Q}$ and let $\text{Cl}^{+}(K)$ be its narrow class group. Recently, the authors \cite{KP} gave a bound for $|\text{Cl}^{+}(K)[2]|$ only in terms of the degree of $K$ and the number of ramifying primes. In the present work we show that this bound is sharp in a wide number of cases. Furthermore, we extend this to ray class groups.


title: Solar Flare Energy Partitioning and Transport -- the Impulsive Phase (a Heliophysics 2050 White Paper)

abstract: Solar flares are a fundamental component of solar eruptive events (SEEs; along with solar energetic particles, SEPs, and coronal mass ejections, CMEs). Flares are the first component of the SEE to impact our atmosphere, which can set the stage for the arrival of the associated SEPs and CME. Magnetic reconnection drives SEEs by restructuring the solar coronal magnetic field, liberating a tremendous amount of energy which is partitioned into various physical manifestations: particle acceleration, mass and magnetic-field eruption, atmospheric heating, and the subsequent emission of radiation as solar flares. To explain and ultimately predict these geoeffective events, the heliophysics community requires a comprehensive understanding of the processes that transform and distribute stored magnetic energy into other forms, including the broadband radiative enhancement that characterises flares. This white paper, submitted to the Heliophysics 2050 Workshop, discusses the flare impulsive phase part of SEEs, setting out the questions that need addressing via a combination of theoretical, modelling, and observational research. In short, by 2050 we must determine the mechanisms of particle acceleration and propagation, and must push beyond the paradigm of energy transport via nonthermal electron beams, to also account for accelerated protons & ions and downward directed Alfven waves.


title: Password similarity using probabilistic data structures

abstract: Passwords should be easy to remember, yet expiration policies mandate their frequent change. Caught in the crossfire between these conflicting requirements, users often adopt creative methods to perform slight variations over time. While easily fooling the most basic checks for similarity, these schemes lead to a substantial decrease in actual security, because leaked passwords, albeit expired, can be effectively exploited as seeds for crackers. This work describes an approach based on Bloom filters to detect password similarity, which can be used to discourage password reuse habits. The proposed scheme intrinsically obfuscates the stored passwords to protect them in case of database leaks, and can be tuned to be resistant to common cryptanalytic techniques, making it suitable for usage on exposed systems.


title: Deuterated Polycyclic Aromatic Hydrocarbons in the Interstellar Medium: The C--D Band Strengths of Mono-Deuterated Species

abstract: Deuterium (D) is one of the light elements created in the big bang. As the Galaxy evolves, the D/H abundance in the interstellar medium (ISM) decreases from its primordial value due to "astration". However, the observed gas-phase D/H abundances of some interstellar sightlines are substantially lower than the expected reduction by astration alone. The missing D could have been depleted onto polycyclic aromatic hydrocarbon (PAH) molecules which are ubiquitous and abundant in interstellar regions. To quantitatively explore the hypothesis of PAHs as a possible reservoir of interstellar D, we compute quantum-chemically the infrared vibrational spectra of mono-deuterated PAHs (and their cations) of various sizes. We find that, as expected, when H in PAHs is replaced by D, the C-H stretching and bending modes at 3.3, 8.6 and 11.3 $\mu$m respectively shift to longer wavelengths at $\sim$4.4, 11.4 and 15.4 $\mu$m by a factor of $\sim$$\sqrt{13/7}$, the difference in reduced mass between the C-H and C-D oscillators. We derive from the computed spectra the mean intrinsic strengths of the 3.3 $\mu$m C-H stretch and 4.4 $\mu$m C--D stretch to be $\langle A_{3.3} \rangle \sim 13.4$ km/mol and $\langle A_{4.4} \rangle \sim 7.4$ km/mol for neutral deuterated PAHs which would dominate the interstellar 3.3 and 4.4 $\mu$m emission. By comparing the computationally-derived mean ratio of $\langle A_{4.4}/A_{3.3} \rangle \sim 0.56$ for neutral PAHs with the mean ratio of the observed intensities of $\langle (I_{4.4}/I_{3.3})_{\rm obs} \rangle \sim 0.019$, we estimate the degree of deuteration (i.e., the fraction of peripheral atoms attached to C in the form of D) to be ~2.4\%, corresponding to a D-enrichment of a factor of ~1200 with respect to the interstellar D/H abundance.


title: Evolutionary Selective Imitation: Interpretable Agents by Imitation Learning Without a Demonstrator

abstract: We propose a new method for training an agent via an evolutionary strategy (ES), in which we iteratively improve a set of samples to imitate: Starting with a random set, in every iteration we replace a subset of the samples with samples from the best trajectories discovered so far. The evaluation procedure for this set is to train, via supervised learning, a randomly initialised neural network (NN) to imitate the set and then execute the acquired policy against the environment. Our method is thus an ES based on a fitness function that expresses the effectiveness of imitating an evolving data subset. This is in contrast to other ES techniques that iterate over the weights of the policy directly. By observing the samples that the agent selects for learning, it is possible to interpret and evaluate the evolving strategy of the agent more explicitly than in NN learning. In our experiments, we trained an agent to solve the OpenAI Gym environment Bipedalwalker-v3 by imitating an evolutionarily selected set of only 25 samples with a NN with only a few thousand parameters. We further test our method on the Procgen game Plunder and show here as well that the proposed method is an interpretable, small, robust and effective alternative to other ES or policy gradient methods.


title: Toward a Comprehensive Model of Snow Crystal Growth: 8. Characterizing Structure-Dependent Attachment Kinetics near -14 C

abstract: In this paper I examine snow crystal growth near -14 C in comparison with a comprehensive model that includes Structure-Dependent Attachment Kinetics (SDAK). Analyzing a series of ice-growth observations in air, I show that the data strongly support the model, which stipulates that basal growth is described by classical terrace nucleation on faceted surfaces in this temperature region. In contrast, prism growth exhibits a pronounced "SDAK dip" that substantially reduces the nucleation barrier on narrow prism facets (relative to that found on broad prism facets). I use these measurements to further characterize and refine the SDAK model, which effectively explains the robust formation of platelike snow crystals in air near 14 C.


title: Bayesian Matrix Completion for Hypothesis Testing

abstract: We aim to infer bioactivity of each chemical by assay endpoint combination, addressing sparsity of toxicology data. We propose a Bayesian hierarchical framework which borrows information across different chemicals and assay endpoints, facilitates out-of-sample prediction of activity for chemicals not yet assayed, quantifies uncertainty of predicted activity, and adjusts for multiplicity in hypothesis testing. Furthermore, this paper makes a novel attempt in toxicology to simultaneously model heteroscedastic errors and a nonparametric mean function, leading to a broader definition of activity whose need has been suggested by toxicologists. Real application identifies chemicals most likely active for neurodevelopmental disorders and obesity.


title: On the Absorption Properties of Metallic Needles

abstract: Needle-like metallic particles have been suggested to explain a wide variety of astrophysical phenomena, ranging from the mid-infrared interstellar extinction to the thermalization of starlight to generate the cosmic microwave background. These suggestions rely on the amplitude and the wavelength dependence of the absorption cross sections of metallic needles. On the absence of an exact solution to the absorption properties of metallic needles, their absorption cross sections are often derived from the antenna approximation. However, it is shown here that the antenna approximation is not an appropriate representation since it violates the Kramers-Kronig relation. Stimulated by the recent discovery of iron whiskers in asteroid Itokawa and graphite whiskers in carbonaceous chondrites, we call for rigorous calculations of the absorption cross sections of metallic needle-like particles, presumably with the discrete dipole approximation. We also call for experimental studies of the formation and growth mechanisms of metallic needle-like particles as well as experimental measurements of the absorption cross sections of metallic needles of various aspect ratios over a wide wavelength range to bound theoretical calculations.