title: Scattering of mesons in quantum simulators

abstract: Simulating real-time evolution in theories of fundamental interactions represents one of the central challenges in contemporary theoretical physics. Cold-atom platforms stand as promising candidates to realize quantum simulations of non-perturbative phenomena in gauge theories, such as vacuum decay and hadron collisions, in prohibitive conditions for direct experiments. In this work, we demonstrate that present-day quantum simulators can imitate linear particle accelerators, giving access to S-matrix measurements of elastic and inelastic meson collisions in low-dimensional Abelian gauge theories. Considering for definiteness a $(1+1)$-dimensional $\mathbb{Z}_2$-lattice gauge theory realizable with Rydberg-atom arrays, we present protocols to observe and measure selected meson-meson scattering processes. We provide a benchmark theoretical study of scattering amplitudes in the regime of large fermion mass, including an exact solution valid for arbitrary coupling strength. This allows us to discuss the occurrence of inelastic scattering channels, featuring the production of new mesons with different internal structures. We present numerical simulations of realistic wavepacket collisions, which reproduce the predicted cross section peaks. This work highlights the potential of quantum simulations to give unprecedented access to real-time scattering dynamics.

categories: cond-mat.quant-gas cond-mat.stat-mech cond-mat.str-el hep-lat quant-ph


abstract: The emergence of passive galaxies in the early Universe results from the interplay among the processes responsible for their rapid assembly and for the abrupt shut-down of their SF. Investigating the individual properties and demographics of early passive galaxies will improve our understanding of these mechanisms. In this work we present a follow-up analysis of the z>3 passive galaxy candidates selected by Merlin et al. (2019) in the CANDELS fields. We begin by first confirming the accuracy of their passive classification by exploiting their sub-mm emission to demonstrate the lack of ongoing SF. Using archival ALMA observations we are able to confirm at least 61% of the observed candidates as passive. While the remainder lack sufficiently deep data for confirmation, we are able to validate the entire sample in a statistical sense. We then estimate the Stellar Mass Function (SMF) of all 101 passive candidates in three redshift bins from z=5 to z=3. We adopt a stepwise approach that has the advantage of taking into account photometric errors, observational incompleteness, and the Eddington bias without any a-posteriori correction. We observe a pronounced evolution in the SMF around z~4, indicating that we are witnessing the emergence of the passive population at this epoch. Massive (M>10^11Msun) passive galaxies, only accounting for a small (<10%) fraction of galaxies at z>4, become dominant at later epochs. Thanks to a combination of photometric quality, sample selection and methodology, we overall find a higher density of passive galaxies than previous works. The comparison with theoretical predictions, despite a qualitative agreement, denotes a still incomplete understanding of the physical processes responsible for the formation of these galaxies. Finally, we extrapolate our results to predict the number of early passive galaxies expected in surveys carried out with future facilities.


title: Systematic improvement of molecular excited state calculations by inclusion of nuclear quantum motion: a mode-resolved picture and the effect of molecular size

abstract: The energies of molecular excited states arise as solutions to the electronic Schr\"{o}dinger equation and are often compared to experiment. At the same time, nuclear quantum motion is known to be important and to induce a red-shift of excited state energies. However, it is thus far unclear whether incorporating nuclear quantum motion in molecular excited state calculations leads to a systematic improvement of their predictive accuracy, making further investigation necessary. Here we present such an investigation by employing two first-principles methods for capturing the effect of quantum fluctuations on excited state energies, which we apply to the Thiel set of organic molecules. We show that accounting for zero-point motion leads to much improved agreement with experiment, compared to `static' calculations which only account for electronic effects, and the magnitude of the red-shift can become as large as 1.36 eV. Moreover, we show that the effect of nuclear quantum motion on excited state energies largely depends on the molecular size, with smaller molecules exhibiting larger red-shifts. Our methodology also makes it possible to analyze the contribution of individual vibrational normal modes to the red-shift of excited state energies, and in several molecules we identify a limited number of modes dominating this effect. Overall, our study provides a foundation for systematically quantifying the shift of excited state energies due to nuclear quantum motion, and for understanding this effect at a microscopic level.


title: Mass Gap in Infinite Derivative Non-local Higgs: Dyson-Schwinger Approach

abstract: We investigate the non-perturbative degrees of freedom in the class of non-local Higgs theories that have been proposed as an ultraviolet completion 4-D Quantum Field Theory (QFT) generalizing the kinetic energy operators to an infinite series of higher derivatives inspired by string field theory. At the perturbative level, the degrees of freedom of non-local Higgs are the same of the local theory. We prove that, at the non-perturbative level, the physical spectrum of the Higgs mass is actually corrected from the "infinite number of derivatives" present in the action. The technique we use permits to derive the set of Dyson-Schwinger equations in differential form. This proves essentially useful when exact solutions to the local equations are known. We show that all the formalism of the local theory involving the Dyson-Schwinger approach extends quite naturally to the non-local case. Using these methods, the spectrum of non-local theories become accessible and predictable in the non-perturbative regimes. We calculate the N-point correlation functions and predict the mass-gap in the spectrum arising purely from the self-interaction and the non-local scale M. The mass gap generated gets damped in the UV and it reaches conformal limit. We discuss some implications of our result in particle physics and cosmology.


title: Is TiO emission present in the ultra-hot Jupiter WASP-33b? A reassessment using the improved ExoMol Toto line list

abstract: [abridged] Efficient absorption of stellar UV and visible radiation by TiO and VO is predicted to drive temperature inversions in the upper atmospheres of hot Jupiters. However, few inversions or detections of TiO or VO have been reported, and results are often contradictory. Using the improved ExoMol Toto line list, we searched for TiO emission in the dayside spectrum of WASP-33b using the same data in which the molecule was previously detected with an older line list at 4.8$\sigma$. We intended to confirm the molecular detection and quantify the signal improvement offered by the ExoMol Toto line list. Data from the High Dispersion Spectrograph on the Subaru Telescope was extracted and reduced in an identical manner to the previous study. Stellar and telluric contamination were then removed. High-resolution TiO emission models of WASP-33b were created using the radiative transfer code petitRADTRANS and cross-correlated with the data. We measure a 4.3$\sigma$ TiO emission signature using the ExoMol Toto models, corresponding to a WASP-33b orbital velocity semi-amplitude of $K_\mathrm{p}=252.9^{+5.0}_{-5.3}\ \mathrm{km\ s^{-1}}$ and a system velocity of $v_\mathrm{sys}=-23.0^{+4.7}_{-4.6}\ \mathrm{km\ s^{-1}}$. Injection-recovery tests using models based on the new and earlier line lists indicate that if the new models provide a perfect match to the planet spectrum, the significance of the TiO detection should have increased by a factor of $\sim$2. Although the TiO signal we find is statistically significant, comparison with previous works makes our result too ambiguous to claim a clear-cut detection. Unexpectedly, the new ExoMol Toto models provide a weaker signal than that found previously, which is offset in $K_\mathrm{p}$-$v_\mathrm{sys}$ space. This sheds some doubt on both detections, especially in light of a recently published TiO non-detection using a different dataset.


title: Luminous Red Nova AT 2019zhd, a new merger in M 31

abstract: We present the follow-up campaign of the luminous red nova (LRN) AT~2019zhd, the third event of this class observed in M 31. The object was followed by several sky surveys for about five months before the outburst, during which it showed a slow luminosity rise. In this phase, the absolute magnitude ranged from M_r=-2.8+-0.2 mag to M_r=-5.6+-0.1 mag. Then, over a four-five day period, AT 2019zhd experienced a major brightening, reaching at peak M_r=-9.61+-0.08 mag, and an optical luminosity of 1.4x10^39 erg/s. After a fast decline, the light curve settled onto a short-duration plateau in the red bands. Although less pronounced, this feature is reminiscent of the second red maximum observed in other LRNe. This phase was followed by a rapid linear decline in all bands. At maximum, the spectra show a blue continuum with prominent Balmer emission lines. The post-maximum spectra show a much redder continuum, resembling that of an intermediate-type star. In this phase, Halpha becomes very weak, Hbeta is no longer detectable and a forest of narrow absorption metal lines now dominate the spectrum. The latest spectra, obtained during the post-plateau decline, show a very red continuum (T_eff ~ 3000 K) with broad molecular bands of TiO, similar to those of M-type stars. The long-lasting, slow photometric rise observed before the peak resembles that of LRN V1309 Sco, which was interpreted as the signature of the common-envelope ejection. The subsequent outburst is likely due to the gas outflow following a stellar merging event. The inspection of archival HST images taken 22 years before the LRN discovery reveals a faint red source (M_F555W=0.21+-0.14 mag, with F555W-F814W = 2.96+-0.12 mag) at the position of AT 2019zhd, which is the most likely quiescent precursor. The source is consistent with expectations for a binary system including a predominant M5-type star.


title: CompModels: A suite of computer model test functions for Bayesian optimization

abstract: The CompModels package for R provides a suite of computer model test functions that can be used for computer model prediction/emulation, uncertainty quantification, and calibration, but in particular, the sequential optimization of computer models. The package is a mix of real-world physics problems, known mathematical functions, and black-box functions that have been converted into computer models with the goal of Bayesian (i.e., sequential) optimization in mind. Likewise, the package contains computer models that represent either the constrained or unconstrained optimization case, each with varying levels of difficulty. In this paper, we illustrate the use of the package with both real-world examples and black-box functions by solving constrained optimization problems via Bayesian optimization. Ultimately, the package is shown to provide users with a source of computer model test functions that are reproducible, shareable, and that can be used for benchmarking of novel optimization methods.


title: The luminous red nova variety: AT 2020hat and AT 2020kog

abstract: We present the results of our monitoring campaigns of the luminous red novae (LRNe) AT 2020hat in NGC 5068 and AT 2020kog in NGC 6106. The two objects were imaged (and detected) before their discovery by routine survey operations. They show a general trend of slow luminosity rise, lasting at least a few months. The subsequent major LRN outbursts were extensively followed in photometry and spectroscopy. The light curves present an initial short-duration peak, followed by a redder plateau phase. AT 2020kog is a moderately luminous event peaking at ~7 x 10^40 erg/s, while AT 2020hat is almost one order of magnitude fainter than AT 2020kog, although it is still more luminous than V838 Mon. In analogy with other LRNe, the spectra of AT 2020kog change significantly with time. They resemble those of type IIn supernovae at early phases, then they become similar to those of K-type stars during the plateau, and to M-type stars at very late phases. In contrast, AT 2020hat already shows a redder continuum at early epochs, and its spectrum shows the late appearance of molecular bands. A moderate-resolution spectrum of AT 2020hat taken at +37 d after maximum shows a forest of narrow P Cygni lines of metals with velocities of 180 km/s, along with an Halpha emission with a full-width at half-maximum velocity of 250 km/s. For AT 2020hat, a robust constraint on its quiescent progenitor is provided by archival images of the Hubble Space Telescope. The progenitor is clearly detected as a mid-K type star, with an absolute magnitude of MF606W = -3.33+-0.09 mag and a colour of F606W-F814W = 1.14+-0.05 mag, which are inconsistent with the expectations from a massive star that could later produce a core-collapse supernova. Although quite peculiar, the two objects nicely match the progenitor versus light curve absolute magnitude correlations discussed in the literature.


title: A Moment for Random Measurements

abstract: Quantum entanglement is one of the core features of quantum theory. While it is typically revealed by measurements along carefully chosen directions, here we review different methods based on so-called random or randomized measurements. Although this approach might seem inefficient at first, sampling correlations in various random directions is a powerful tool to study properties which are invariant under local-unitary transformations. Based on random measurements, entanglement can be detected and characterized without a shared reference frame between the observers or even if local reference frames cannot be defined. This overview article discusses different methods using random measurements to detect genuine multipartite entanglement and to distinguish SLOCC classes. Furthermore, it reviews how measurement directions can efficiently be obtained based on spherical designs.


title: Variational design principles for nonequilibrium colloidal assembly

abstract: Using large deviation theory and principles of stochastic optimal control, we show that rare molecular dynamics trajectories conditioned on assembling a specific target structure encode a set of interactions and external forces that lead to enhanced stability of that structure. Such a relationship can be formulated into a variational principle, for which we have developed an associated optimization algorithm and have used it to determine optimal forces for targeted self-assembly within nonequilibrium steady-states. We illustrate this perspective on inverse design in a model of colloidal cluster assembly within linear shear flow. We find that colloidal clusters can be assembled with high yield using specific short-range interactions of tunable complexity. Shear decreases the yields of rigid clusters, while small values of shear increase the yields of nonrigid clusters. The enhancement or suppression of the yield due to shear is rationalized with a generalized linear response theory. By studying 21 unique clusters made of 6, 7 or 8 particles, we uncover basic design principles for targeted assembly out of equilibrium.


title: Whistler Waves in the foot of Quasi-Perpendicular Super-Critical Shocks

abstract: Whistler waves are thought to play an essential role in the dynamics of collisionless shocks. We use the magnetospheric multiscale (MMS) spacecraft to study whistler waves around the lower hybrid frequency, upstream of 11 quasi-perpendicular super-critical shocks. We apply the 4-spacecraft timing method to unambiguously determine the wave vector $\mathbf{k}$ of whistler waves. We find that the waves are oblique to the background magnetic field with a wave-normal angle between $20^{\circ}$ and $42^{\circ}$, a wavelength around 100 km which is close to the ion inertial length. We also find that $\mathbf{k}$ is predominantly in the same plane as the magnetic field and the normal to the shock. By combining this precise knowledge of $\mathbf{k}$ with high-resolution measurements of the 3D ion velocity distribution we show that a reflected ion beam is in resonance with the waves, opening up the possibility for wave-particle interaction between the reflected ions and the observed whistlers. The linear stability analysis of a system mimicking the observed distribution, suggests that such a system can produce the observed waves.


title: Neutrino mass constraints beyond linear order: cosmology dependence and systematic biases

abstract: We demonstrate the impact on forecasted neutrino mass constraints of extending galaxy clustering and CMB lensing predictions from linear to next-to-leading-order power spectra. The redshift-space 1-loop power spectrum model we adopt requires an additional four free bias parameters, a velocity bias parameter and two new stochastic parameters. These additional nuisance parameters appreciably weaken the constraints on $M_\nu$. CMB lensing plays a significant role in helping to alleviate these degeneracies and tighten the final constraints. The constraint on the optical depth to reionisation $\tau$ has a strong effect on the constraint on $M_\nu$, but only when CMB lensing is included in the analysis to keep the degeneracies with the nuisance parameters under control. We also extract constraints when 1) using the BAO signature only as a distance probe, and 2) isolating the scale-dependence of the power spectrum, which, as shown in previous work, provides a cosmology-independent probe of $M_\nu$. All constraints except the latter remain strongly sensitive to the assumption of a flat $\Lambda$CDM universe. We perform an analysis of the magnitude of the shift introduced in the inferred $M_\nu$ value when neglecting nonlinear corrections, and show that, for a Euclid-like survey, this shift becomes roughly equal to the 1$\sigma$ constraint itself even with a conservative cut-off scale of $k_{max} = 0.1~h~{\rm Mpc}^{-1}$. We also perform a calculation of the appropriate expected bias in neutrino mass caused by not including the next, 2-loop order and expect a shift of only about 20% of the 1$\sigma$ error for $k_{max}=0.2~h~{\rm Mpc}^{-1}$ in a Euclid-like survey.


title: How network properties and epidemic parameters influence stochastic SIR dynamics on scale-free random networks

abstract: With the premise that social interactions are described by power-law distributions, we study a SIR stochastic dynamic on a static scale-free random network generated via configuration model. We verify our model with respect to deterministic considerations and provide a theoretical result on the probability of the extinction of the disease. Based on this calibration, we explore the variability in disease spread by stochastic simulations. In particular, we demonstrate how important epidemic indices change as a function of the contagiousness of the disease and the connectivity of the network. Our results quantify the role of starting node degree in determining these indices, commonly used to describe epidemic spread.


title: The Impact of Data on the Stability of Learning-Based Control- Extended Version

abstract: Despite the existence of formal guarantees for learning-based control approaches, the relationship between data and control performance is still poorly understood. In this paper, we propose a Lyapunov-based measure for quantifying the impact of data on the certifiable control performance. By modeling unknown system dynamics through Gaussian processes, we can determine the interrelation between model uncertainty and satisfaction of stability conditions. This allows us to directly asses the impact of data on the provable stationary control performance, and thereby the value of the data for the closed-loop system performance. Our approach is applicable to a wide variety of unknown nonlinear systems that are to be controlled by a generic learning-based control law, and the results obtained in numerical simulations indicate the efficacy of the proposed measure.


title: Synchronization within synchronization: transients and intermittency in ecological networks

abstract: Transients are fundamental to ecological systems with significant implications to management, conservation, and biological control. We uncover a type of transient synchronization behavior in spatial ecological networks whose local dynamics are of the chaotic, predator-prey type. In the parameter regime where there is phase synchronization among all the patches, complete synchronization (i.e., synchronization in both phase and amplitude) can arise in certain pairs of patches as determined by the network symmetry - henceforth the phenomenon of "synchronization within synchronization." Distinct patterns of complete synchronization coexist but, due to intrinsic instability or noise, each pattern is a transient and there is random, intermittent switching among the patterns in the course of time evolution. The probability distribution of the transient time is found to follow an algebraic scaling law with a divergent average transient lifetime. Based on symmetry considerations, we develop a stability analysis to understand these phenomena. The general principle of symmetry can also be exploited to explain previously discovered, counterintuitive synchronization behaviors in ecological networks.


title: On tameness of zonoids

abstract: We prove that in a globally subanalytic family of convex bodies the set of zonoids is log-analytic, and in particular it is definable in the o-minimal structure generated by globally subanalytic sets and the graph of the exponential function.


title: Measurements of Nuclear Modification Factors of $B^{0}_{s}$ and $B^+$ Mesons in PbPb Collisions with the CMS Experiment

abstract: Beauty quarks are considered as one of the best probes of the strongly interacting medium created in relativistic heavy-ion collisions because they are predominantly produced via initial hard scatterings. Measurements of B meson production provide information about the diffusion of beauty quarks and the flavor dependence of in-medium energy loss. In these studies, clarifying the hadronization mechanism is crucial for understanding the transport properties of beauty quarks. Measurements of $B^{0}_{s}$ production can shed light on the mechanisms of beauty recombination in the medium and provide information about strangeness enhancement in the quark-gluon plasma. In this talk, we will present a new measurement of the ratio of $B^{0}_{s}$ to $B^+$ mesons in PbPb collisions at 5.02 TeV with the CMS detector, using data recorded in 2018.


title: ATSal: An Attention Based Architecture for Saliency Prediction in 360 Videos

abstract: The spherical domain representation of 360 video/image presents many challenges related to the storage, processing, transmission and rendering of omnidirectional videos (ODV). Models of human visual attention can be used so that only a single viewport is rendered at a time, which is important when developing systems that allow users to explore ODV with head mounted displays (HMD). Accordingly, researchers have proposed various saliency models for 360 video/images. This paper proposes ATSal, a novel attention based (head-eye) saliency model for 360\degree videos. The attention mechanism explicitly encodes global static visual attention allowing expert models to focus on learning the saliency on local patches throughout consecutive frames. We compare the proposed approach to other state-of-the-art saliency models on two datasets: Salient360! and VR-EyeTracking. Experimental results on over 80 ODV videos (75K+ frames) show that the proposed method outperforms the existing state-of-the-art.


title: Entanglement Spectrum Crossings Reveal non-Hermitian Dynamical Topology

abstract: The development of non-Hermitian topological band theory has led to observations of novel topological phenomena in effectively classical, driven and dissipative systems. However, for open quantum many-body systems, the absence of a ground state presents a challenge to define robust signatures of non-Hermitian topology. We show that such a signature is provided by crossings in the time evolution of the entanglement spectrum. These crossings occur in quenches from the trivial to the topological phase of a driven-dissipative Kitaev chain that is described by a Markovian quantum master equation in Lindblad form. At the topological transition, which can be crossed either by changing parameters of the Hamiltonian of the system or by increasing the strength of dissipation, the time scale at which the first entanglement spectrum crossing occurs diverges with a dynamical critical exponent of $\epsilon = 1/2$. We corroborate these numerical findings with an exact analytical solution of the quench dynamics for a spectrally flat postquench Liouvillian. This exact solution suggests an interpretation of the topological quench dynamics as a fermion parity pump. Our work thus reveals signatures of non-Hermitian topology which are unique to quantum many-body systems and cannot be emulated in classical simulators of non-Hermitian wave physics.


title: LSTM-based Traffic Load Balancing and Resource Allocation for an Edge System

abstract: The massive deployment of small cell Base Stations (SBSs) empowered with computing capabilities presents one of the most ingenious solutions adopted for 5G cellular networks towards meeting the foreseen data explosion and the ultra-low latency demanded by mobile applications. This empowerment of SBSs with Multi-access Edge Computing (MEC) has emerged as a tentative solution to overcome the latency demands and bandwidth consumption required by mobile applications at the network edge. The MEC paradigm offers a limited amount of resources to support computation, thus mandating the use of intelligence mechanisms for resource allocation. The use of green energy for powering the network apparatuses (e.g., Base Stations (BSs), MEC servers) has attracted attention towards minimizing the carbon footprint and network operational costs. However, due to their high intermittency and unpredictability, the adoption of learning methods is a requisite. Towards intelligent edge system management, this paper proposes a Green-based Edge Network Management (GENM) algorithm, which is a online edge system management algorithm for enabling green-based load balancing in BSs and energy savings within the MEC server. The main goal is to minimize the overall energy consumption and guarantee the Quality of Service (QoS) within the network. To achieve this, the GENM algorithm performs dynamic management of BSs, autoscaling and reconfiguration of the computing resources, and on/off switching of the fast tunable laser drivers coupled with location-aware traffic scheduling in the MEC server. The obtained simulation results validate our analysis and demonstrate the superior performance of GENM compared to a benchmark algorithm.


title: Photoionization Models for High Density Gas

abstract: Relativistically broadened and redshifted 6.4 -- 6.9 keV iron K lines are observed from many accretion powered objects, including X-ray binaries and active galactic nuclei (AGN). Existence of gas close to the central engine implies large radiation intensities and correspondingly large gas densities if the gas is to remain partially ionized. Simple estimates indicate that high gas densities are needed to allow survival of iron against ionization. These are high enough that rates for many atomic processes are affected by mechanisms related to interactions with nearby ions and electrons. Radiation intensities are high enough that stimulated processes can be important. Most models currently in use for interpreting relativistic lines use atomic rate coefficients designed for use at low densities and neglect stimulated processes. In our work so far we have presented atomic structure calculations with the goal of providing physically appropriate models at densities consistent with line-emitting gas near compact objects. In this paper we apply these rates to photoionization calculations, and produce ionization balance curves and X-ray emissivities and opacities which are appropriate for high densities and high radiation intensities. The final step in our program will be presented in a subsequent paper: Model atmosphere calculations which incorporate these rates into synthetic spectra.


title: Energetic spin-polarized proton beams from two-stage coherent acceleration in laser-driven plasma

abstract: We propose a scheme to overcome the great challenge of polarization loss in spin-polarized ion acceleration. When a petawatt laser pulse penetrates through a compound plasma target consisting of a double layer slab and prepolarized hydrogen halide gas, a strong forward moving quasistatic longitudinal electric field is constructed by the self-generated laser-driven plasma. This field with a varying drift velocity efficiently boosts the prepolarized protons via a two-stage coherent acceleration process. Its merit is not only achieving a highly energetic beam but also eliminating the undesired polarization loss of the accelerated protons. We study the proton dynamics via Hamiltonian analyses, specifically deriving the threshold of triggering the two-stage coherent acceleration. To confirm the theoretical predictions, we perform three-dimensional PIC simulations, where unprecedented proton beams with energy approximating half GeV and polarization ratio $\sim$ 94\% are obtained.


title: Nested Mixture of Experts: Cooperative and Competitive Learning of Hybrid Dynamical System

abstract: Model-based reinforcement learning (MBRL) algorithms can attain significant sample efficiency but require an appropriate network structure to represent system dynamics. Current approaches include white-box modeling using analytic parameterizations and black-box modeling using deep neural networks. However, both can suffer from a bias-variance trade-off in the learning process, and neither provides a structured method for injecting domain knowledge into the network. As an alternative, gray-box modeling leverages prior knowledge in neural network training but only for simple systems. In this paper, we devise a nested mixture of experts (NMOE) for representing and learning hybrid dynamical systems. An NMOE combines both white-box and black-box models while optimizing bias-variance trade-off. Moreover, an NMOE provides a structured method for incorporating various types of prior knowledge by training the associative experts cooperatively or competitively. The prior knowledge includes information on robots' physical contacts with the environments as well as their kinematic and dynamic properties. In this paper, we demonstrate how to incorporate prior knowledge into our NMOE in various continuous control domains, including hybrid dynamical systems. We also show the effectiveness of our method in terms of data-efficiency, generalization to unseen data, and bias-variance trade-off. Finally, we evaluate our NMOE using an MBRL setup, where the model is integrated with a model-based controller and trained online.


title: Primordial black holes and secondary gravitational waves from inflationary model with a non-canonical kinetic term

abstract: With the enhancement mechanism provided by a noncanonical kinetic term with a peak, the amplitude of primordial curvature perturbations can be enhanced by seven orders of magnitude at small scales while keeping to be consistent with observations at large scales. The peak function and inflationary potential are not restricted in this mechanism. We use the Higgs model and T-model as examples to show how abundant primordial black hole dark matter with different mass and scalar induced secondary gravitational waves with different peak frequency are generated. We also show that the enhanced power spectrum for the primordial curvature perturbations and the energy density of the scalar induced secondary gravitational waves can have either a sharp peak or a broad peak. The primordial black holes with the mass around $10^{-14}-10^{-12} M_{\odot}$ produced with the enhancement mechanism can make up almost all dark matter, and the scalar induced secondary gravitational waves accompanied with the production of primordial black holes can be tested by the pulsar timing arrays and spaced based gravitational wave observatory. Therefore, the mechanism can be tested by primordial black hole dark matter and gravitational wave observations.


title: Lightweight Data Fusion with Conjugate Mappings

abstract: We present an approach to data fusion that combines the interpretability of structured probabilistic graphical models with the flexibility of neural networks. The proposed method, lightweight data fusion (LDF), emphasizes posterior analysis over latent variables using two types of information: primary data, which are well-characterized but with limited availability, and auxiliary data, readily available but lacking a well-characterized statistical relationship to the latent quantity of interest. The lack of a forward model for the auxiliary data precludes the use of standard data fusion approaches, while the inability to acquire latent variable observations severely limits direct application of most supervised learning methods. LDF addresses these issues by utilizing neural networks as conjugate mappings of the auxiliary data: nonlinear transformations into sufficient statistics with respect to the latent variables. This facilitates efficient inference by preserving the conjugacy properties of the primary data and leads to compact representations of the latent variable posterior distributions. We demonstrate the LDF methodology on two challenging inference problems: (1) learning electrification rates in Rwanda from satellite imagery, high-level grid infrastructure, and other sources; and (2) inferring county-level homicide rates in the USA by integrating socio-economic data using a mixture model of multiple conjugate mappings.


title: Large Scale Neural Architecture Search with Polyharmonic Splines

abstract: Neural Architecture Search (NAS) is a powerful tool to automatically design deep neural networks for many tasks, including image classification. Due to the significant computational burden of the search phase, most NAS methods have focused so far on small, balanced datasets. All attempts at conducting NAS at large scale have employed small proxy sets, and then transferred the learned architectures to larger datasets by replicating or stacking the searched cells. We propose a NAS method based on polyharmonic splines that can perform search directly on large scale, imbalanced target datasets. We demonstrate the effectiveness of our method on the ImageNet22K benchmark[16], which contains 14 million images distributed in a highly imbalanced manner over 21,841 categories. By exploring the search space of the ResNet [23] and Big-Little Net ResNext [11] architectures directly on ImageNet22K, our polyharmonic splines NAS method designed a model which achieved a top-1 accuracy of 40.03% on ImageNet22K, an absolute improvement of 3.13% over the state of the art with similar global batch size [15].


title: Power Adaptation for Vector Parameter Estimation according to Fisher Information based Optimality Criteria

abstract: The optimal power adaptation problem is investigated for vector parameter estimation according to various Fisher information based optimality criteria. By considering an observation model that involves a linear transformation of the parameter vector and an additive noise component with an arbitrary probability distribution, six different optimal power allocation problems are formulated based on Fisher information based objective functions. Via optimization theoretic approaches, various closed-form solutions are derived for the proposed problems. Also, the results are extended to cases in which nuisance parameters exist in the system model or certain types of nonlinear transformations are applied on the parameter vector. Numerical examples are presented to investigate performance of the proposed power allocation strategies.


title: SReachTools Kernel Module: Data-Driven Stochastic Reachability Using Hilbert Space Embeddings of Distributions

abstract: We present algorithms for performing data-driven stochastic reachability as an addition to SReachTools, an open-source stochastic reachability toolbox. Our method leverages a class of machine learning techniques known as kernel embeddings of distributions to approximate the safety probabilities for a wide variety of stochastic reachability problems. By representing the probability distributions of the system state as elements in a reproducing kernel Hilbert space, we can learn the "best fit" distribution via a simple regularized least-squares problem, and then compute the stochastic reachability safety probabilities as simple linear operations. This technique admits finite sample bounds and has known convergence in probability. We implement these methods as part of SReachTools, and demonstrate their use on a double integrator system, on a million-dimensional repeated planar quadrotor system, and a cart-pole system with a black-box neural network controller.


title: Noether and Hilbert (metric) energy-momentum tensors are not, in general, equivalent

abstract: Multiple methods for deriving the energy-momentum tensor for a physical theory exist in the literature. The most common methods are to use Noether's first theorem with the 4-parameter Poincar\'{e} translation, or to write the action in a curved spacetime and perform variation with respect to the metric tensor, then return to a Minkowski spacetime. These are referred to as the Noether and Hilbert (metric/ curved space/ variational) energy-momentum tensors, respectively. In electrodynamics and other simple models, the Noether and Hilbert methods yield the same result. Due to this fact, it is often asserted that these methods are generally equivalent for any theory considered, and that this gives physicists a freedom in using either method to derive an energy-momentum tensor depending on the problem at hand. $\dots$ For spin-2, the ideal candidate to check this equivalence for a more complicated model, there exist many energy-momentum tensors in the literature, none of which are gauge invariant, so it is not clear which expression one hopes to obtain from the Noether and Hilbert approaches unlike in the case of e.g. electrodynamics. It has been shown, however, that the linearized Gauss-Bonnet gravity model (second order derivatives, second rank tensor potential) has an energy-momentum tensor that is unique, gauge invariant, symmetric, conserved, and trace-free when derived from Noether's first theorem (all the same properties of the physical energy-momentum tensor of electrodynamics). This makes it the ideal candidate to check if the Noether and Hilbert methods coincide for a more complicated model. It is proven here using this model as a counterexample, by direct calculation, that the Noether and Hilbert energy-momentum tensors are not, in general, equivalent.


title: Equilibrium sequences of differentially rotating stars with post-merger-like rotational profiles

abstract: We present equilibrium sequences of rotating relativistic stars, constructed with a new rotation law that was proposed by Uryu et al. (2017). We choose rotational parameters motivated by simulations of binary neutron star merger remnants, but otherwise adopt a cold, relativistic N=1 polytropic EOS, in order to perform a detailed comparison to published equilibrium sequences that used the Komatsu, Eriguchi and Hachisu (1989) rotation law. We find a small influence of the choice of rotation law on the mass of the equilibrium models and a somewhat larger influence on their radius. The versatility of the new rotation law allows us to construct models that have a similar rotational profile and axis ratio as observed for merger remnants, while at the same time being quasi-spherical. More specifically, we construct equilibrium sequence variations with different degrees of differential rotation and identify type A and type C solutions, similar to the corresponding types in the classification of Ansorg, Gondek-Rosinska and Villain (2009). While our models are highly accurate solutions of the fully general relativistic structure equations, we demonstrate that for models relevant to merger remnants the IWM-CFC approximation still maintains an acceptable accuracy.


title: Generating the Goeritz group of $S^3$

abstract: In 1980 J. Powell \cite{Po} proposed that five specific elements sufficed to generate the Goeritz group for any genus Heegaard splitting of the 3-sphere. Here we prove that a natural expansion of Powell's proposed generators, to include all eyeglass twists and all topological conjugates of Powell's generators does suffice.


title: Meta Variational Monte Carlo

abstract: An identification is found between meta-learning and the problem of determining the ground state of a randomly generated Hamiltonian drawn from a known ensemble. A model-agnostic meta-learning approach is proposed to solve the associated learning problem and a preliminary experimental study of random Max-Cut problems indicates that the resulting Meta Variational Monte Carlo accelerates training and improves convergence.


title: Adversarial Training for EM Classification Networks

abstract: We present a novel variant of Domain Adversarial Networks with impactful improvements to the loss functions, training paradigm, and hyperparameter optimization. New loss functions are defined for both forks of the DANN network, the label predictor and domain classifier, in order to facilitate more rapid gradient descent, provide more seamless integration into modern neural networking frameworks, and allow previously unavailable inferences into network behavior. Using these loss functions, it is possible to extend the concept of 'domain' to include arbitrary user defined labels applicable to subsets of the training data, the test data, or both. As such, the network can be operated in either 'On the Fly' mode where features provided by the feature extractor indicative of differences between 'domain' labels in the training data are removed or in 'Test Collection Informed' mode where features indicative of difference between 'domain' labels in the combined training and test data are removed (without needing to know or provide test activity labels to the network). This work also draws heavily from previous works on Robust Training which draws training examples from a L_inf ball around the training data in order to remove fragile features induced by random fluctuations in the data. On these networks we explore the process of hyperparameter optimization for both the domain adversarial and robust hyperparameters. Finally, this network is applied to the construction of a binary classifier used to identify the presence of EM signal emitted by a turbopump. For this example, the effect of the robust and domain adversarial training is to remove features indicative of the difference in background between instances of operation of the device - providing highly discriminative features on which to construct the classifier.


title: Bridging Physics-based and Data-driven modeling for Learning Dynamical Systems

abstract: How can we learn a dynamical system to make forecasts, when some variables are unobserved? For instance, in COVID-19, we want to forecast the number of infected and death cases but we do not know the count of susceptible and exposed people. While mechanics compartment models are widely used in epidemic modeling, data-driven models are emerging for disease forecasting. We first formalize the learning of physics-based models as AutoODE, which leverages automatic differentiation to estimate the model parameters. Through a benchmark study on COVID-19 forecasting, we notice that physics-based mechanistic models significantly outperform deep learning. Our method obtains a 57.4% reduction in mean absolute errors for 7-day ahead COVID-19 forecasting compared with the best deep learning competitor. Such performance differences highlight the generalization problem in dynamical system learning due to distribution shift. We identify two scenarios where distribution shift can occur: changes in data domain and changes in parameter domain (system dynamics). Through systematic experiments on several dynamical systems, we found that deep learning models fail to forecast well under both scenarios. While much research on distribution shift has focused on changes in the data domain, our work calls attention to rethink generalization for learning dynamical systems.


title: $Ab$ $initio$ hydrogen dynamics and the morphology of voids in amorphous silicon

abstract: This paper presents an $ab$ $initio$ study of hydrogen dynamics inside nanometer-size voids in $a$-Si within the framework of the density-functional theory for a varying hydrogen load of 10 to 30 H atoms/void at the low and high temperature of 400 K and 700 K, respectively. Using the local density approximation and its generalized-gradient counterpart, the dynamics of hydrogen atoms inside the voids are examined with an emphasis on the diffusion of H atoms/molecules, and the resulting nanostructural changes of the void surfaces. The results from simulations suggest that the microstructure of the hydrogen distribution on the void surfaces and the morphology of the voids are characterized by the presence of a significant number of monohydride Si-H bonds, along with a few dihydride Si-H$_2$ configurations. The study also reveals that a considerable number of (about 10--45 at.%) total H atoms inside voids can appear as H$_2$ molecules for a hydrogen load of 10--30 H atoms/void. The approximate shape of the voids is addressed from a knowledge of the positions of the void-surface atoms using the convex-hull approximation and the Gaussian broadening of the pseudo-atomic surfaces of Si and H atoms.


title: Gradualizing the Calculus of Inductive Constructions

abstract: We investigate gradual variations on the Calculus of Inductive Construction (CIC) for swifter prototyping with imprecise types and terms. We observe, with a no-go theorem, a crucial tradeoff between graduality and the key properties of normalization and closure of universes under dependent product that CIC enjoys. Beyond this Fire Triangle of Graduality, we explore the gradualization of CIC with three different compromises, each relaxing one edge of the Fire Triangle. We develop a parametrized presentation of Gradual CIC (GCIC) that encompasses all three variations, and develop their metatheory. We first present a bidirectional elaboration of GCIC to a dependently-typed cast calculus, CastCIC, which elucidates the interrelation between typing, conversion, and the gradual guarantees. We use a syntactic model of CastCIC to inform the design of a safe, confluent reduction, and establish, when applicable, normalization. We study the static and dynamic gradual guarantees as well as the stronger notion of graduality with embedding-projection pairs formulated by New and Ahmed, using appropriate semantic model constructions. This work informs and paves the way towards the development of malleable proof assistants and dependently-typed programming languages.


title: Finite Horizon Discrete Models for Multi-Agent Control Systems with Coupled Dynamics

abstract: The goal of this paper is to obtain online abstractions for coupled multi-agent systems in a decentralized manner. A discrete model which captures the motion capabilities of each agent is derived over a bounded time-horizon, by discretizing a corresponding overapproximation of the agent's reachable states. The individual abstractions' composition provides a correct representation of the coupled continuous system over the horizon and renders the approach appropriate for control synthesis under high-level specifications which are assigned to the agents over this time window. Sufficient conditions are also provided for the space and time discretization to guarantee the derivation of deterministic abstractions with tunable transition capabilities.


title: Anomalous decay rate of quasinormal modes in Reissner-Nordstr\"om black holes

abstract: The anomalous decay rate of the quasinormal modes occurs when the longest-lived modes are the ones with higher angular number. Such behaviour has been recently studied in different static spacetimes, for scalar and fermionic perturbations, being observed in both cases. In this work, we extend the existent studies to the charged spacetimes, namely, the Reissner-Nordstr\"om, the Reissner-Nordstr\"om-de Sitter and the Reissner-Nordstr\"om-Anti-de Sitter black holes. We show that the anomalous decay rate behaviour of the scalar field perturbations is present for every charged geometry in the photon sphere modes, with the existence of a critical scalar field mass whenever $\Lambda \geq 0$. In general, this critical value of mass increases with the raise of the black hole charge, thus rendering a minimum in the Schwarzschild limit. We also study the dominant mode/family for the massless and massive scalar field in these geometries showing a non-trivial dominance of the spectra that depends on the black hole mass and charge.


title: An exact power series expansion for the Kallen-Sabry vacuum polarization potential

abstract: In this work it is presented an exact power series formula for the Kallen-Sabry vacuum polarization potential.


title: Notes on equivariant homology with constant coefficients

abstract: In this paper, for a finite group, we discuss a method for calculating equivariant homology with constant coefficients. We apply it to completely calculate the geometric fixed points of the equivariant spectrum representing equivariant (co)homology with constant coefficients. We also treat a more complicated example of inverting the standard representation in the equivariant homology of split extraspecial groups at the prime 2.


title: On the Target Pebbling Conjecture

abstract: Graph pebbling is a network optimization model for satisfying vertex demands with vertex supplies (called pebbles), with partial loss of pebbles in transit. The pebbling number of a demand in a graph is the smallest number for which every placement of that many supply pebbles satisfies the demand. The Target Conjecture (Herscovici-Hester-Hurlbert, 2009) posits that the largest pebbling number of a demand of fixed size $t$ occurs when the demand is entirely stacked on one vertex. This truth of this conjecture could be useful for attacking many open problems in graph pebbling, including the famous conjecture of Graham (1989) involving graph products. It has been verified for complete graphs, cycles, cubes, and trees. In this paper we prove the conjecture for 2-paths and Kneser graphs over pairs.


title: Single Microhole per Pixel in CMOS Image Sensor with Enhanced Optical Sensitivity in Near-Infrared

abstract: Silicon photodiode based CMOS sensors with backside-illumination for 300 to 1000 nm wavelength range were studied. We showed that a single hole in the photodiode increases the optical efficiency of the pixel. In near-infrared wavelengths, the enhancement allows 70% absorption in a 3 microns thick Si. It is 4x better than for the flat pixel. We compared different shapes and sizes of single holes and holes arrays. We have shown that a certain size and shape in single holes pronounce better optical efficiency enhancement. The crosstalk was successfully reduced with trenches between pixels. We optimized the trenches to achieve minimal pixel separation for 1.12 microns pixel.


title: Semantic SLAM with Autonomous Object-Level Data Association

abstract: It is often desirable to capture and map semantic information of an environment during simultaneous localization and mapping (SLAM). Such semantic information can enable a robot to better distinguish places with similar low-level geometric and visual features and perform high-level tasks that use semantic information about objects to be manipulated and environments to be navigated. While semantic SLAM has gained increasing attention, there is little research on semanticlevel data association based on semantic objects, i.e., object-level data association. In this paper, we propose a novel object-level data association algorithm based on bag of words algorithm, formulated as a maximum weighted bipartite matching problem. With object-level data association solved, we develop a quadratic-programming-based semantic object initialization scheme using dual quadric and introduce additional constraints to improve the success rate of object initialization. The integrated semantic-level SLAM system can achieve high-accuracy object-level data association and real-time semantic mapping as demonstrated in the experiments. The online semantic map building and semantic-level localization capabilities facilitate semantic-level mapping and task planning in a priori unknown environment.


title: A New Window into Planet Formation and Migration: Refractory-to-Volatile Elemental Ratios in Ultra-hot Jupiters

abstract: A primary goal of exoplanet characterization is to use a planet's current composition to understand how that planet formed. For example, the C/O ratio has long been recognized as carrying important information on the chemistry of volatile species. Refractory elements, like Fe, Mg, and Si, are usually not considered in this conversation because they condense into solids like Fe(s) or MgSiO$_3$ and would be removed from the observable, gaseous atmosphere in exoplanets cooler than about 2000~K. However, planets hotter than about 2000~K, called ultra-hot Jupiters (UHJs), are warm enough to largely avoid the condensation of refractory species. In this paper, we explore the insight that the measurement of refractory abundances can provide into a planet's origins. Through refractory-to-volatile elemental abundance ratios, we can estimate a planet's atmospheric rock-to-ice fraction and constrain planet formation and migration scenarios. We first relate a planet's present-day refractory-to-volatile ratio to its rock-to-ice ratio from formation using various compositional models for the rocky and icy components of the protoplanetary disk. We discuss potential confounding factors like the sequestration of heavy metals in the core and condensation. We then show such a measurement using atmospheric retrievals of the low-resolution UV-IR transmission spectrum of WASP-121b with PETRA, from which we estimate a refractory-to-volatile ratio of 5.0$^{+6.0}_{-2.7}\times$ solar and a rock-to-ice ratio greater than 2/3. This result is consistent with significant atmospheric enrichment by rocky planetismals. Lastly, we discuss the rich future potential for measuring refractory-to-volatile ratios in ultra-hot Jupiters with the arrival of JWST and by combining observations at low- and high-resolution.


title: An unexpected high concentration for the dark substructure in the gravitational lens SDSSJ0946+1006

abstract: The presence of an invisible substructure has previously been detected in the gravitational lens galaxy SDSSJ0946+1006 through its perturbation of the lensed images. Using flexible models for the main halo and the subhalo perturbation to fit the lensed images, we demonstrate that the subhalo has an extraordinarily high central density and steep density slope. The inferred concentration for the subhalo is well above the expected scatter in concentrations for $\Lambda$CDM halos of similar mass. We robustly infer the subhalo's projected mass within 1 kpc to be $\sim 2$-$3.7\times 10^9$M$_\odot$ at $>$95% CL for all our lens models, while the average slope of the subhalo's projected density profile over the radial range 0.75-1.25 kpc is constrained to be steeper than isothermal ($\gamma_{2D} \lesssim -1$). By modeling the subhalo light directly, we infer a conservative upper bound on its luminosity $L_V < 1.2\times 10^8L_\odot$ at 95% CL, which shows that the perturber is dark matter dominated. To compare to $\Lambda$CDM expectations, we analyze subhalos within analogues of lensing galaxies in the Illustris TNG100-1 simulation over many lines of sight, and find hundreds of subhalos that achieve a projected mass within 1 kpc of $\gtrsim 2\times10^9M_\odot$. However, less than 1% of the mock observations yield a log-slope steep enough to be consistent with our lensing models, and they $all$ have stellar masses in excess of that allowed by observations by about an order of magnitude or more. We conclude that the presence of such a dark, highly concentrated subhalo is unexpected in a $\Lambda$CDM universe. Finally, we show that this tension with CDM is not significantly reduced if the perturber is assumed to be a line-of-sight structure, rather than a subhalo.


title: The Weyl-Mellin quantization map for $\kappa$-Minkowski Noncommutative Spacetime

abstract: We present a quantization of the functions of spacetime, i.e.\ a map, analog to Weyl map, which reproduces the $\kappa$-Minkowski commutation relations, and it has the desirable properties of mapping square integrable funcions into Hilbert-Schmidt operators, as well as real functions into self-adjoint operators. The map is based on Mellin transform on radial and time coordinates. The map also define a deformed $*$ product which we discuss with examples.


title: Inferring the concentration of dark matter subhalos perturbing strongly lensed images

abstract: We demonstrate that the perturbations of strongly lensed images by low-mass dark matter subhalos are significantly impacted by the concentration of the perturbing subhalo. For subhalo concentrations expected in $\Lambda$CDM, significant constraints on the concentration can be obtained at HST resolution for subhalos with masses larger than about $10^{10}M_\odot$. Constraints are also possible for lower mass subhalos, if their concentrations are higher than the expected scatter in CDM. We also find that the concentration of lower mass perturbers down to $\sim 10^8M_\odot$ can be well-constrained with a resolution of $\sim 0.01''$, which is achievable with long-baseline interferometry. Subhalo concentration also plays a critical role in the detectability of a perturbation, such that only high concentration perturbers with mass $\lesssim 10^9M_\odot$ are likely to be detected at HST resolution. If scatter in the $\Lambda$CDM mass-concentration relation is not accounted for during lens modeling, the inferred subhalo mass can be biased by up to a factor of 3(6) for subhalos of mass $10^9 M_\odot$($10^{10} M_\odot$); this bias can be eliminated if one varies both mass and concentration during lens fitting. Alternatively, one may robustly infer the projected mass within the subhalo's perturbation radius, defined by its distance to the critical curve of the lens being perturbed. With a sufficient number of detections, these strategies will make it possible to constrain the halo mass-concentration relation at low masses in addition to the mass function, offering a probe of dark matter physics as well as the small-scale primordial power spectrum.


title: Solving path dependent PDEs with LSTM networks and path signatures

abstract: Using a combination of recurrent neural networks and signature methods from the rough paths theory we design efficient algorithms for solving parametric families of path dependent partial differential equations (PPDEs) that arise in pricing and hedging of path-dependent derivatives or from use of non-Markovian model, such as rough volatility models in Jacquier and Oumgari, 2019. The solutions of PPDEs are functions of time, a continuous path (the asset price history) and model parameters. As the domain of the solution is infinite dimensional many recently developed deep learning techniques for solving PDEs do not apply. Similarly as in Vidales et al. 2018, we identify the objective function used to learn the PPDE by using martingale representation theorem. As a result we can de-bias and provide confidence intervals for then neural network-based algorithm. We validate our algorithm using classical models for pricing lookback and auto-callable options and report errors for approximating both prices and hedging strategies.


title: Counterexample to the Bohigas Conjecture for Transmission Through aOne-Dimensional Lattice

abstract: Resonances in particle transmission through a 1D finite lattice are studied in the presence of a finite number of impurities. Although this is a one-dimensional system that is classically integrable and has no chaos, studying the statistical properties of the spectrum such as the level spacing distribution and the spectral rigidity shows quantum chaos signatures. Using a dimensionless parameter that reflects the degree of state localization, we demonstrate how the transition from regularity to chaos is affected by state localization. The resonance positions are calculated using both the Wigner-Smithtime-delay and a Siegert state method, which are in good agreement. Our results give evidence for the existence of quantum chaos in one dimension which is a counter-example to the Bohigas-Giannoni-Schmit conjecture.


title: Wafer-Scale Epitaxial Modulation of Quantum Dot Density

abstract: Precise control of the properties of semiconductor quantum dots (QDs) is vital for creating novel devices for quantum photonics and advanced opto-electronics. Suitable low QD-density for single QD devices and experiments are challenging to control during epitaxy and are typically found only in limited regions of the wafer. Here, we demonstrate how conventional molecular beam epitaxy (MBE) can be used to modulate the density of optically active QDs in one- and two- dimensional patterns, while still retaining excellent quality. We find that material thickness gradients during layer-by-layer growth result in surface roughness modulations across the whole wafer. Growth on such templates strongly influences the QD nucleation probability. We obtain density modulations between 1 and 10 QDs/${\mu}m^{2}$ and periods ranging from several millimeters down to at least a few hundred microns. This novel method is universal and expected to be applicable to a wide variety of different semiconductor material systems. We apply the method to enable growth of ultra-low noise QDs across an entire 3-inch semiconductor wafer.


title: Phase Transitions in Germanium Telluride Nanoparticle Phase-Change Materials Studied by Time-Resolved X-Ray Diffraction

abstract: Germanium telluride (GeTe), a phase-change material, is known to exhibit four different structural phases: three at room temperature (one amorphous and two crystalline, $\alpha$ and $\gamma$) and one at high temperature (crystalline $\beta$). Because transitions between the amorphous and crystalline phases lead to significant changes in material properties (e.g., refractive index and resistivity), GeTe has been investigated as a phase-change material for photonics, thermoelectrics, ferroelectrics, and spintronics. Consequently, the temperature-dependent phase transitions in GeTe have been studied for bulk and thin-film GeTe, both fabricated by sputtering. Colloidal synthesis of nanoparticles offers a more flexible fabrication approach for amorphous and crystalline GeTe. These nanoparticles are known to exhibit size-dependent properties, such as an increased crystallization temperature for the amorphous-to-$\alpha$ transition in sub-10\,nm GeTe particles. The $\alpha$-to-$\beta$ phase transition is also expected to vary with size, but this effect has not yet been investigated for GeTe. Here, we report time-resolved X-ray diffraction of GeTe nanoparticles with different diameters and from different synthetic protocols. We observe a non-volatile amorphous-to-$\alpha$ transition between 210$^{\circ}$C and 240$^{\circ}$C and a volatile $\alpha$-to-$\beta$ transition between 370$^{\circ}$C and 420$^{\circ}$C. The latter transition was reversible and repeatable. While the transition temperatures are shifted relative to the values known for bulk GeTe, the nanoparticle-based samples still exhibit the same structural phases reported for sputtered GeTe. Thus, colloidal GeTe maintains the same general phase behavior as bulk GeTe while allowing for more flexible and accessible fabrication. Therefore, nanoparticle-based GeTe films show great potential for applications, such as in active photonics.


title: Distributed Robust State Estimation for Hybrid AC/DC Distribution Systems using Multi-Source Data

abstract: Hybrid AC/DC distribution systems are becoming a popular means to accommodate the increasing penetration of distributed energy resources and flexible loads. This paper proposes a distributed and robust state estimation (DRSE) method for hybrid AC/DC distribution systems using multiple sources of data. In the proposed distributed implementation framework, a unified robust linear state estimation model is derived for each AC and DC regions, where the regions are connected via AC/DC converters and only limited information exchange is needed. To enhance the estimation accuracy of the areas with low measurement coverage, a deep neural network (DNN) is used to extract hidden system statistical information and allow deriving nodal power injections that keep up with the real-time measurement update rate. This provides the way of integrating smart meter data, SCADA measurements and zero injections together for state estimation. Simulations on two hybrid AC/DC distribution systems show that the proposed DRSE has only slight accuracy loss by the linearization formulation but offers robustness of suppressing bad data automatically, as well as benefits of improving computational efficiency.


title: Landmark and IMU Data Fusion: Systematic Convergence Geometric Nonlinear Observer for SLAM and Velocity Bias

abstract: Navigation solutions suitable for cases when both autonomous robot's pose (\textit{i.e}., attitude and position) and its environment are unknown are in great demand. Simultaneous Localization and Mapping (SLAM) fulfills this need by concurrently mapping the environment and observing robot's pose with respect to the map. This work proposes a nonlinear observer for SLAM posed on the manifold of the Lie group of $\mathbb{SLAM}_{n}\left(3\right)$, characterized by systematic convergence, and designed to mimic the nonlinear motion dynamics of the true SLAM problem. The system error is constrained to start within a known large set and decay systematically to settle within a known small set. The proposed estimator is guaranteed to achieve predefined transient and steady-state performance and eliminate the unknown bias inevitably present in velocity measurements by directly using measurements of angular and translational velocity, landmarks, and information collected by an inertial measurement unit (IMU). Experimental results obtained by testing the proposed solution on a real-world dataset collected by a quadrotor demonstrate the observer's ability to estimate the six-degrees-of-freedom (6 DoF) robot pose and to position unknown landmarks in three-dimensional (3D) space. Keywords: Simultaneous Localization and Mapping, Nonlinear filter for SLAM, Nonlinear filter for SLAM on Matrix Lie group, pose, asymptotic stability, prescribed performance, adaptive estimate, feature, inertial measurement unit, inertial vision unit, IMU, SE(3), SO(3), noise.


title: A mixed parameter formulation with applications to linear viscoelasticity

abstract: In this work we propose and analyze an abstract parameter dependent model written as a mixed variational formulation based on Volterra integrals of second kind. For the analysis, we consider a suitable adaptation to the classic mixed theory in the Volterra equations setting, and prove the well posedness of the resulting mixed viscoelastic formulation. Error estimates are derived, using the available results for Volterra equations, where all the estimates are independent of the perturbation parameter. We consider an application of the developed theory in a viscoelastic Timoshenko beam, and report numerical experiments in order to assess the independence of the perturbation parameter.


title: Ergodic and non-ergodic many-body dynamics in strongly nonlinear lattices

abstract: The study of non-linear oscillator chains in classical many-body dynamics has a storied history going back to the seminal work of Fermi, Pasta, Ulam and Tsingou (FPUT). We introduce a new family of such systems which consist of chains of $N$ harmonically coupled particles with the non-linearity introduced by confining the motion of each individual particle to a box/stadium with hard walls. The stadia are arranged on a one dimensional lattice but they individually do not have to be one dimensional thus permitting the introduction of chaos already at the lattice scale. For the most part we study the case where the motion is entirely one dimensional. We find that the system exhibits a mixed phase space for any finite value of $N$. Computations of Lyapunov spectra at randomly picked phase space locations and a direct comparison between Hamiltonian evolution and phase space averages indicate that the regular regions of phase space are not significant at large system sizes. While the continuum limit of our model is itself a singular limit of the integrable sinh-Gordon theory, we do not see any evidence for the kind of non-ergodicity famously seen in the FPUT work. Finally, we examine the chain with particles confined to two dimensional stadia where the individual stadium is already chaotic, and find a much more chaotic phase space at small system sizes.


title: Convergent subseries of divergent series

abstract: Let $\mathscr{X}$ be the set of positive real sequences $x=(x_n)$ such that the series $\sum_n x_n$ is divergent. For each $x \in \mathscr{X}$, let $\mathcal{I}_x$ be the collection of all $A\subseteq \mathbf{N}$ such that the subseries $\sum_{n \in A}x_n$ is convergent. Moreover, let $\mathscr{A}$ be the set of sequences $x \in \mathscr{X}$ such that $\lim_n x_n=0$ and $\mathcal{I}_x\neq \mathcal{I}_y$ for all sequences $y=(y_n) \in \mathscr{X}$ with $\liminf_n y_{n+1}/y_n>0$. We show that $\mathscr{A}$ is comeager and that contains uncountably many sequences $x$ which generate pairwise nonisomorphic ideals $\mathcal{I}_x$. This answers, in particular, an open question recently posed by M. Filipczak and G. Horbaczewska.


title: Low coherence-induced resonance in double-layer structures having parity-time symmetry

abstract: We derive simple formulae for the transmittance $T$ and reflectance $R$ of Gaussian-Schell beams incident upon any stratified dielectric structure by using second-order classical coherence theory in the space-frequency picture. The formalism is applied to a particular structure consisting of a double-layer, with balanced gain and loss, satisfying the parity-time symmetry conditions. It is shown that sources with a low degree of spatial coherence, on the order of the wavelength, can induce large resonant peaks in the transmitted and reflected amplitudes. The resonance peaks vanish as the spatial coherence increases.


title: Assessment and Linear Programming under Fuzzy Conditions

abstract: A new fuzzy method is developed using triangular/trapezoidal fuzzy numbers for evaluating a group's mean performance, when qualitative grades instead of numerical scores are used for assessing its members' individual performance. Also, a new technique is developed for solving Linear Programming problems with fuzzy coefficients and everyday life applications are presented to illustrate our results.


title: The node cop-win reliability of unicyclic and bicyclic graphs

abstract: Various models to quantify the reliability of a network have been studied where certain components of the graph may fail at random and the probability that the remaining graph is connected is the proxy for reliability. In this work we introduce a strengthening of one of these models by considering the probability that the remaining graph is not just connected but also cop-win. A graph is cop-win if one cop can guarantee capture of a fleeing robber in the well-studied pursuit-evasion game of Cops and Robber. More precisely, for a graph $G$ with nodes that are operational independently with probability $p$ and edges that are operational if and only if both of their endpoints are operational, the node cop-win reliability of $G$, denoted $\text{NCRel}(G,p)$, is the probability that the operational nodes induce a cop-win subgraph of $G$. It is then of interest to find graphs $G$ with $n$ nodes and $m$ edges such that $\text{NCRel}(G,p)\ge\text{NCRel}(H,p)$ for all $p\in[0,1]$ and all graphs $H$ with $n$ nodes and $m$ edges. Such a graph is called uniformly most reliable. We show that uniformly most reliable graphs exist for unicyclic and bicyclic graphs, respectively. This is in contrast to the fact that there are no known sparse graphs maximizing the corresponding notion of node reliability.


title: Linearization for High-Speed Current-Steering DACs Using Neural Networks

abstract: This paper proposes a novel foreground linearization scheme for a high-speed CS-DAC. The technique leverages neural networks (NNs) to derive a LUT that maps the inverse of the DAC transfer characteristic onto the input codes. The algorithm is shown to improve conventional methods by at least 6dB in terms of intermodulation (IM) performance for frequencies up to 9GHz on a state-of-the-art 10-bit CS-DAC operating at 40.96GS/s (gigasamples-per-second) in 14nm CMOS.


title: On the Benefits of Multiple Gossip Steps in Communication-Constrained Decentralized Optimization

abstract: In decentralized optimization, it is common algorithmic practice to have nodes interleave (local) gradient descent iterations with gossip (i.e. averaging over the network) steps. Motivated by the training of large-scale machine learning models, it is also increasingly common to require that messages be {\em lossy compressed} versions of the local parameters. In this paper, we show that, in such compressed decentralized optimization settings, there are benefits to having {\em multiple} gossip steps between subsequent gradient iterations, even when the cost of doing so is appropriately accounted for e.g. by means of reducing the precision of compressed information. In particular, we show that having $O(\log\frac{1}{\epsilon})$ gradient iterations {with constant step size} - and $O(\log\frac{1}{\epsilon})$ gossip steps between every pair of these iterations - enables convergence to within $\epsilon$ of the optimal value for smooth non-convex objectives satisfying Polyak-\L{}ojasiewicz condition. This result also holds for smooth strongly convex objectives. To our knowledge, this is the first work that derives convergence results for nonconvex optimization under arbitrary communication compression.


title: Interacting hard-sphere fluids in an external field

abstract: We present a new method for studying equilibrium properties of interacting fluids in an arbitrary external field. The fluid is composed of monodisperse spherical particles with hard-core repulsion and additional interactions of arbitrary shape and limited range. Our method of analysis is exact in one dimension and provides demonstrably good approximations in higher dimensions. It can cope with homogeneous and heterogeneous environments. We derive an equation for the pair distribution function. The solution, to be evaluated numerically, in general, or analytically for special cases, enters expressions for the entropy and free energy functionals. For some one-dimensional systems, our approach yields analytic solutions, reproducing available exact results from different approaches.


title: Study of Resource Amount Configuration for Automatic Application Offloading

abstract: In recent years, utilization of heterogeneous hardware other than small core CPU such as GPU, FPGA or many core CPU is increasing. However, when using heterogeneous hardware, barriers of technical skills such as OpenMP, CUDA and OpenCL are high. Based on that, I have proposed environment-adaptive software that enables automatic conversion, configuration, and high performance operation of once written code, according to the hardware to be placed. However, although the conversion of the code according to the migration destination environment has been studied so far, there has been no research to properly set the resource amount. In this paper, as a new element of environment adaptive software, in order to operate the application with high cost performance, I study a method to optimize the resource amount of CPUs and offload devices.


title: On the Topological Complexity of Maps

abstract: We define and develop a homotopy invariant notion for the topological complexity of a map $f:X \to Y$, denoted TC($f$), that interacts with TC($X$) and TC($Y$) in the same way cat($f$) interacts with cat($X$) and cat($Y$). Furthermore, TC($f$) and cat($f$) satisfy the same inequalities as TC($X$) and cat($X$). We compare it to other invariants defined in the papers [15,16,17,18,20]. We apply TC($f$) to studying group homomorphisms $f:H\to G$.


title: What do we expect from Multiple-choice QA Systems?

abstract: The recent success of machine learning systems on various QA datasets could be interpreted as a significant improvement in models' language understanding abilities. However, using various perturbations, multiple recent works have shown that good performance on a dataset might not indicate performance that correlates well with human's expectations from models that "understand" language. In this work we consider a top performing model on several Multiple Choice Question Answering (MCQA) datasets, and evaluate it against a set of expectations one might have from such a model, using a series of zero-information perturbations of the model's inputs. Our results show that the model clearly falls short of our expectations, and motivates a modified training approach that forces the model to better attend to the inputs. We show that the new training paradigm leads to a model that performs on par with the original model while better satisfying our expectations.


title: Efficient space-time reduced order model for linear dynamical systems in Python using less than 120 lines of code

abstract: A classical reduced order model (ROM) for dynamical problems typically involves only the spatial reduction of a given problem. Recently, a novel space-time ROM for linear dynamical problems has been developed, which further reduces the problem size by introducing a temporal reduction in addition to a spatial reduction without much loss in accuracy. The authors show an order of a thousand speed-up with a relative error of less than 0.00001 for a large-scale Boltzmann transport problem. In this work, we present for the first time the derivation of the space-time Petrov-Galerkin projection for linear dynamical systems and its corresponding block structures. Utilizing these block structures, we demonstrate the ease of construction of the space-time ROM method with two model problems: 2D diffusion and 2D convection diffusion, with and without a linear source term. For each problem, we demonstrate the entire process of generating the full order model (FOM) data, constructing the space-time ROM, and predicting the reduced-order solutions, all in less than 120 lines of Python code. We compare our Petrov-Galerkin method with the traditional Galerkin method and show that the space-time ROMs can achieve O(100) speed-ups with O(0.001) to O(0.0001) relative errors for these problems. Finally, we present an error analysis for the space-time Petrov-Galerkin projection and derive an error bound, which shows an improvement compared to traditional spatial Galerkin ROM methods.


title: Adler-Bell-Jackiw anomaly in VSR electrodynamics

abstract: In this paper, we examine the problem of anomalies of the fermionic currents in the context of the very special relativity (VSR). We consider the VSR contributions to the triangle amplitude $\left\langle J_{5}^{\lambda} J^{\mu} J^{\nu} \right\rangle $, which allows the evaluation of the vector and axial Ward identities. Actually, we observe that the VSR nonlocal effects respect the vector Ward identity, and it also contributes in a very interesting and unique way for the Adler-Bell-Jackiw anomaly.


title: Very Deep VAEs Generalize Autoregressive Models and Can Outperform Them on Images

abstract: We present a hierarchical VAE that, for the first time, generates samples quickly while outperforming the PixelCNN in log-likelihood on all natural image benchmarks. We begin by observing that, in theory, VAEs can actually represent autoregressive models, as well as faster, better models if they exist, when made sufficiently deep. Despite this, autoregressive models have historically outperformed VAEs in log-likelihood. We test if insufficient depth explains why by scaling a VAE to greater stochastic depth than previously explored and evaluating it CIFAR-10, ImageNet, and FFHQ. In comparison to the PixelCNN, these very deep VAEs achieve higher likelihoods, use fewer parameters, generate samples thousands of times faster, and are more easily applied to high-resolution images. Qualitative studies suggest this is because the VAE learns efficient hierarchical visual representations. We release our source code and models at https://github.com/openai/vdvae.


title: A PDE Characterization of Anisotropic Hardy Spaces

abstract: We obtain a differential characterization for the anisotropic Hardy space $H_A^p$ by identifying it with a parabolic Hardy space associated with a general continuous group. This allows $H_A^p$ to be defined using a parabolic differential equation of Calderon and Torchinsky. We also provide a classification of dilations corresponding to equivalent anisotropic Hardy spaces with respect to linear transformations.


title: Self-Supervised learning with cross-modal transformers for emotion recognition

abstract: Emotion recognition is a challenging task due to limited availability of in-the-wild labeled datasets. Self-supervised learning has shown improvements on tasks with limited labeled datasets in domains like speech and natural language. Models such as BERT learn to incorporate context in word embeddings, which translates to improved performance in downstream tasks like question answering. In this work, we extend self-supervised training to multi-modal applications. We learn multi-modal representations using a transformer trained on the masked language modeling task with audio, visual and text features. This model is fine-tuned on the downstream task of emotion recognition. Our results on the CMU-MOSEI dataset show that this pre-training technique can improve the emotion recognition performance by up to 3% compared to the baseline.


title: Effects of Human vs. Automatic Feedback on Students' Understanding of AI Concepts and Programming Style

abstract: The use of automatic grading tools has become nearly ubiquitous in large undergraduate programming courses, and recent work has focused on improving the quality of automatically generated feedback. However, there is a relative lack of data directly comparing student outcomes when receiving computer-generated feedback and human-written feedback. This paper addresses this gap by splitting one 90-student class into two feedback groups and analyzing differences in the two cohorts' performance. The class is an intro to AI with programming HW assignments. One group of students received detailed computer-generated feedback on their programming assignments describing which parts of the algorithms' logic was missing; the other group additionally received human-written feedback describing how their programs' syntax relates to issues with their logic, and qualitative (style) recommendations for improving their code. Results on quizzes and exam questions suggest that human feedback helps students obtain a better conceptual understanding, but analyses found no difference between the groups' ability to collaborate on the final project. The course grade distribution revealed that students who received human-written feedback performed better overall; this effect was the most pronounced in the middle two quartiles of each group. These results suggest that feedback about the syntax-logic relation may be a primary mechanism by which human feedback improves student outcomes.


title: Upgraded W-Net with Attention Gates and its Application in Unsupervised 3D Liver Segmentation

abstract: Segmentation of biomedical images can assist radiologists to make a better diagnosis and take decisions faster by helping in the detection of abnormalities, such as tumors. Manual or semi-automated segmentation, however, can be a time-consuming task. Most deep learning based automated segmentation methods are supervised and rely on manually segmented ground-truth. A possible solution for the problem would be an unsupervised deep learning based approach for automated segmentation, which this research work tries to address. We use a W-Net architecture and modified it, such that it can be applied to 3D volumes. In addition, to suppress noise in the segmentation we added attention gates to the skip connections. The loss for the segmentation output was calculated using soft N-Cuts and for the reconstruction output using SSIM. Conditional Random Fields were used as a post-processing step to fine-tune the results. The proposed method has shown promising results, with a dice coefficient of 0.88 for the liver segmentation compared against manual segmentation.


title: Viability of Optical Coherence Tomography for Iris Presentation Attack Detection

abstract: In this paper, we propose the use of Optical Coherence Tomography (OCT) imaging for the problem of iris presentation attack (PA) detection. We assess its viability by comparing its performance with respect to traditional iris imaging modalities, viz., near-infrared (NIR) and visible spectrum. OCT imaging provides a cross-sectional view of an eye, whereas traditional imaging provides 2D iris textural information. PA detection is performed using three state-of-the-art deep architectures (VGG19, ResNet50 and DenseNet121) to differentiate between bonafide and PA samples for each of the three imaging modalities. Experiments are performed on a dataset of 2,169 bonafide, 177 Van Dyke eyes and 360 cosmetic contact images acquired using all three imaging modalities under intra-attack (known PAs) and cross-attack (unknown PAs) scenarios. We observe promising results demonstrating OCT as a viable solution for iris presentation attack detection.


title: Experimental demonstration of multimode microresonator sensing by machine learning

abstract: A multimode microcavity sensor based on a self-interference microring resonator is demonstrated experimentally. The proposed multimode sensing method is implemented by recording wideband transmission spectra that consist of multiple resonant modes. It is different from the previous dissipative sensing scheme, which aims at measuring the transmission depth changes of a single resonant mode in a microcavity. Here, by combining the dissipative sensing mechanism and the machine learning algorithm, the multimode sensing information extracted from a broadband spectrum can be efficiently fused to estimate the target parameter. The multimode sensing method is immune to laser frequency noises and robust against system imperfection, thus our work presents a great step towards practical applications of microcavity sensors outside the research laboratory. The voltage applied across the microheater on the chip was adjusted to bring its influence on transmittance through the thermo-optic effects. As a proof-of-principle experiment, the voltage was detected by the multimode sensing approach. The experimental results demonstrate that the limit of detection of the multimode sensing by the general regression neural network is reduced to 6.7% of that of single-mode sensing within a large measuring range.


title: Predicting Patient COVID-19 Disease Severity by means of Statistical and Machine Learning Analysis of Blood Cell Transcriptome Data

abstract: Introduction: For COVID-19 patients accurate prediction of disease severity and mortality risk would greatly improve care delivery and resource allocation. There are many patient-related factors, such as pre-existing comorbidities that affect disease severity. Since rapid automated profiling of peripheral blood samples is widely available, we investigated how such data from the peripheral blood of COVID-19 patients might be used to predict clinical outcomes. Methods: We thus investigated such clinical datasets from COVID-19 patients with known outcomes by combining statistical comparison and correlation methods with machine learning algorithms; the latter included decision tree, random forest, variants of gradient boosting machine, support vector machine, K-nearest neighbour and deep learning methods. Results: Our work revealed several clinical parameters measurable in blood samples, which discriminated between healthy people and COVID-19 positive patients and showed predictive value for later severity of COVID-19 symptoms. We thus developed a number of analytic methods that showed accuracy and precision for disease severity and mortality outcome predictions that were above 90%. Conclusions: In sum, we developed methodologies to analyse patient routine clinical data which enables more accurate prediction of COVID-19 patient outcomes. This type of approaches could, by employing standard hospital laboratory analyses of patient blood, be utilised to identify, COVID-19 patients at high risk of mortality and so enable their treatment to be optimised.


title: Lattice Quantum Electrodynamics in (3+1)-dimensions at finite density with Tensor Networks

abstract: Gauge theories are of paramount importance in our understanding of fundamental constituents of matter and their interactions. However, the complete characterization of their phase diagrams and the full understanding of non-perturbative effects are still debated, especially at finite charge density, mostly due to the sign-problem affecting Monte Carlo numerical simulations. Here, we report the Tensor Network simulation of a three dimensional lattice gauge theory in the Hamiltonian formulation including dynamical matter: Using this sign-problem-free method, we simulate the ground states of a compact Quantum Electrodynamics at zero and finite charge densities, and address fundamental questions such as the characterization of collective phases of the model, the presence of a confining phase at large gauge coupling, and the study of charge-screening effects.

categories: hep-lat cond-mat.stat-mech cond-mat.str-el hep-th quant-ph


title: Efficient stream-based Max-Min diversification with minimal failure rate

abstract: The stream-based Max-Min diversification problem concerns the task of selecting a limited number of diverse instances from a data stream. The nature of the problem demands immediate and irrevocable decisions. The set-wise diversity to be maximized is the minimum distance among any pair of the selected instances. Standard algorithmic approaches for sequential selection disregard the possibility of selection failures, which is the situation where the last instances of the stream are picked by default to prevent having an incomplete selection. This defect can be catastrophic for the Max-Min diversification objective. In this paper we present the Failure Rate Minimization (FRM) algorithm that allows the selection of a set of disparate instances while reducing significantly the probability of having failures. This is achieved by means of both analytical and empirical techniques. FRM is put in comparison with relevant algorithms from the literature through simulations on real datasets, where we demonstrate its efficiency and low time complexity.


title: Teaching Key Machine Learning Principles Using Anti-learning Datasets

abstract: Much of the teaching of machine learning focuses on iterative hill-climbing approaches and the use of local knowledge to gain information leading to local or global maxima. In this paper we advocate the teaching of alternative methods of generalising to the best possible solution, including a method called anti-learning. By using simple teaching methods, students can achieve a deeper understanding of the importance of validation on data excluded from the training process and that each problem requires its own methods to solve. We also exemplify the requirement to train a model using sufficient data by showing that different granularities of cross-validation can yield very different results.


title: A method for evaluating options for motif detection in electricity meter data

abstract: Investigation of household electricity usage patterns, and matching the patterns to behaviours, is an important area of research given the centrality of such patterns in addressing the needs of the electricity industry. Additional knowledge of household behaviours will allow more effective targeting of demand side management (DSM) techniques.


title: Resistance Scaling on $4N$-Carpets

abstract: The $4N$ carpets are a class of infinitely ramified self-similar fractals with a large group of symmetries. For a $4N$-carpet $F$, let $\{F_n\}_{n \geq 0}$ be the natural decreasing sequence of compact pre-fractal approximations with $\cap_nF_n=F$. On each $F_n$, let $\mathcal{E}(u, v) = \int_{F_N} \nabla u \cdot \nabla v \, dx$ be the classical Dirichlet form and $u_n$ be the unique harmonic function on $F_n$ satisfying a mixed boundary value problem corresponding to assigning a constant potential between two specific subsets of the boundary. Using a method introduced by Barlow and Bass (1990), we prove a resistance estimate of the following form: there is $\rho=\rho(N) > 1$ such that $\mathcal{E}(u_n, u_n)\rho^{n}$ is bounded above and below by positive constants independent of $n$. Such estimates have implications for the existence and scaling properties of Dirichlet forms on $F$.


title: Metasurface tessellation for reconfigurable electromagnetic scattering

abstract: Metasurfaces have attracted significant research interest owing to their unprecedented control over the spatial distributions of electromagnetic fields. Herein we propose the concept of metasurface tessellation to achieve reconfigurable scattering functions. Square meta-tiles, composed of identical structures, are arranged to fill a surface. The electromagnetic scattering of the tiled surface is determined by the orientation distribution of the meta-tiles. We present three typical cases of meta-tiles consisting of binary elements to realize several distinct scattering patterns. This study provides an alternative method to build reconfigurable and multi-functional metasurface devices without external stimuli and complicated fabrication.


title: An Accurate Numerical Method and Algorithm for Constructing Solutions of Chaotic Systems

abstract: In various fields of natural science, the chaotic systems of differential equations are considered more than 50 years. The correct prediction of the behaviour of solutions of dynamical model equations is important in understanding of evolution process and reduce uncertainty. However, often used numerical methods are unable to do it on large time segments. In this article, the author considers the modern numerical method and algorithm for constructing solutions of chaotic systems on the example of tumor growth model. Also a modification of Benettin's algorithm presents for calculation of Lyapunov exponents.


title: Charging-as-a-Service: On-demand battery delivery for light-duty electric vehicles for mobility service

abstract: This study presents an innovative solution for powering electric vehicles, named Charging-as-a-Service (CaaS), that concerns the potential large-scale adoption of light-duty electric vehicles (LDEV) in the Mobility-as-a-Service (MaaS) industry. Analogous to the MaaS, the core idea of the CaaS is to dispatch service vehicles (SVs) that carry modular battery units (MBUs) to provide LDEVs for mobility service with on-demand battery delivery. The CaaS system is expected to tackle major bottlenecks of a large-scale LDEV adoption in the MaaS industry due to the lack of charging infrastructure and excess waiting and charging time. A hybrid agent-based simulation model (HABM) is developed to model the dynamics of the CaaS system with SV agents, and a trip-based stationary charging probability distribution is introduced to simulate the generation of charging demand for LDEVs. Two dispatching algorithms are further developed to support the optimal operation of the CaaS. The model is validated by assuming electrifying all 13,000 yellow taxis in New York City (NYC) that follow the same daily trip patterns. Multiple scenarios are analyzed under various SV fleet sizes and dispatching strategies. The results suggest that optimal deployment of 250 SVs may serve the LDEV fleet in NYC with an average waiting time of 5 minutes, save the travel distance at over 50 miles per minute, and gain considerable profits of up to $50 per minute. This study offers significant insights into the feasibility, service efficiency, and financial sustainability for deploying city-wide CaaS systems to power the electric MaaS industry.


title: Enhancing Poaching Predictions for Under-Resourced Wildlife Conservation Parks Using Remote Sensing Imagery

abstract: Illegal wildlife poaching is driving the loss of biodiversity. To combat poaching, rangers patrol expansive protected areas for illegal poaching activity. However, rangers often cannot comprehensively search such large parks. Thus, the Protection Assistant for Wildlife Security (PAWS) was introduced as a machine learning approach to help identify the areas with highest poaching risk. As PAWS is deployed to parks around the world, we recognized that many parks have limited resources for data collection and therefore have scarce feature sets. To ensure under-resourced parks have access to meaningful poaching predictions, we introduce the use of publicly available remote sensing data to extract features for parks. By employing this data from Google Earth Engine, we also incorporate previously unavailable dynamic data to enrich predictions with seasonal trends. We automate the entire data-to-deployment pipeline and find that, with only using publicly available data, we recuperate prediction performance comparable to predictions made using features manually computed by park specialists. We conclude that the inclusion of satellite imagery creates a robust system through which parks of any resource level can benefit from poaching risks for years to come.