title: Performance of a Deep Neural Network at Detecting North Atlantic Right Whale Upcalls

abstract: Passive acoustics provides a powerful tool for monitoring the endangered North Atlantic right whale ($Eubalaena$ $glacialis$), but robust detection algorithms are needed to handle diverse and variable acoustic conditions and differences in recording techniques and equipment. Here, we investigate the potential of deep neural networks for addressing this need. ResNet, an architecture commonly used for image recognition, is trained to recognize the time-frequency representation of the characteristic North Atlantic right whale upcall. The network is trained on several thousand examples recorded at various locations in the Gulf of St.\ Lawrence in 2018 and 2019, using different equipment and deployment techniques. Used as a detection algorithm on fifty 30-minute recordings from the years 2015-2017 containing over one thousand upcalls, the network achieves recalls up to 80%, while maintaining a precision of 90%. Importantly, the performance of the network improves as more variance is introduced into the training dataset, whereas the opposite trend is observed using a conventional linear discriminant analysis approach. Our work demonstrates that deep neural networks can be trained to identify North Atlantic right whale upcalls under diverse and variable conditions with a performance that compares favorably to that of existing algorithms.


title: Semi-supervised ASR by End-to-end Self-training

abstract: While deep learning based end-to-end automatic speech recognition (ASR) systems have greatly simplified modeling pipelines, they suffer from the data sparsity issue. In this work, we propose a self-training method with an end-to-end system for semi-supervised ASR. Starting from a Connectionist Temporal Classification (CTC) system trained on the supervised data, we iteratively generate pseudo-labels on a mini-batch of unsupervised utterances with the current model, and use the pseudo-labels to augment the supervised data for immediate model update. Our method retains the simplicity of end-to-end ASR systems, and can be seen as performing alternating optimization over a well-defined learning objective. We also perform empirical investigations of our method, regarding the effect of data augmentation, decoding beamsize for pseudo-label generation, and freshness of pseudo-labels. On a commonly used semi-supervised ASR setting with the WSJ corpus, our method gives 14.4% relative WER improvement over a carefully-trained base system with data augmentation, reducing the performance gap between the base system and the oracle system by 50%.


title: Distinguishing freezing and thawing dark energy models through measurements of the fine-structure constant

abstract: Mapping the behaviour of dark energy is a pressing task for observational cosmology. Phenomenological classification divides dynamical dark energy models into freezing and thawing, depending on whether the dark energy equation of state is approaching or moving away from $w=p/\rho=-1$. Moreover, in realistic dynamical dark energy models the dynamical degree of freedom is expected to couple to the electromagnetic sector, leading to variations of the fine-structure constant $\alpha$. We discuss the feasibility of distinguishing between the freezing and thawing classes of models with current and forthcoming observational facilities and using a parametrisation of the dark energy equation of state, which can have either behaviour, introduced by Mukhanov as fiducial paradigm. We illustrate how freezing and thawing models lead to different redshift dependencies of $\alpha$, and use a combination of current astrophysical observations and local experiments to constrain this class of models, improving the constraints on the key coupling parameter by more than a factor of two, despite considering a more extended parameter space than the one used in previous studies. We also briefly discuss the improvements expected from future facilities and comment on the practical limitations of this class of parametrisations. In particular, we show that sufficiently sensitive data can distinguish between freezing and thawing models, at least if one assumes that the relevant parameter space does not include phantom dark energy models.


title: Creating Realistic Power Distribution Networks using Interdependent Road Infrastructure

abstract: It is well known that physical interdependencies exist between networked civil infrastructures such as transportation and power system networks. In order to analyze complex nonlinear correlations between such networks, datasets pertaining to such real infrastructures are required. However, such data are not readily available due to their proprietary nature. This work proposes a methodology to generate realistic synthetic power distribution networks for a given geographical region. A network generated in this manner is not the actual distribution system, but its functionality is very similar to the real distribution network. The synthetic network connects high voltage substations to individual residential consumers through primary and secondary distribution networks. Here, the distribution network is generated by solving an optimization problem which minimizes the overall length of the network subject to structural and power flow constraints. This work also incorporates identification of long high voltage feeders originating from substations and connecting remotely situated customers in rural geographic locations while maintaining voltage regulation within acceptable limits. The proposed methodology is applied to the state of Virginia and creates synthetic distribution networks which are validated by comparing them to actual power distribution networks at the same location.


title: Large-Scale Optical Reservoir Computing for Spatiotemporal Chaotic Systems Prediction

abstract: Reservoir computing is a relatively recent computational paradigm that originates from a recurrent neural network and is known for its wide range of implementations using different physical technologies. Large reservoirs are very hard to obtain in conventional computers, as both the computation complexity and memory usage grow quadratically. We propose an optical scheme performing reservoir computing over very large networks potentially being able to host several millions of fully connected photonic nodes thanks to its intrinsic properties of parallelism and scalability. Our experimental studies confirm that, in contrast to conventional computers, the computation time of our optical scheme is only linearly dependent on the number of photonic nodes of the network, which is due to electronic overheads, while the optical part of computation remains fully parallel and independent of the reservoir size. To demonstrate the scalability of our optical scheme, we perform for the first time predictions on large spatiotemporal chaotic datasets obtained from the Kuramoto-Sivashinsky equation using optical reservoirs with up to 50 000 optical nodes. Our results are extremely challenging for conventional von Neumann machines, and they significantly advance the state of the art of unconventional reservoir computing approaches, in general.


title: Study of Entropy-Diffusion Relation in a Deterministic Hamiltonian System through Microscopic Analysis

abstract: Although an intimate relation between entropy and diffusion has been advocated for many years and even seems to have been verified in theory and experiments, a quantitatively reliable study, and any derivation of an algebraic relation between the two does not seem to exist. Here we explore the nature of this entropy-diffusion relation in three deterministic systems where an accurate estimate of both can be carried out. We study three deterministic model systems, (a) the motion of a single point particle with constant energy in a two-dimensional periodic potential energy landscape, (b) the same in regular Lorentz gas where a point particle with constant energy moves between collisions with hard disc scatterers and (c) motion of a point particle among the boxes with small apertures. These models, introduced by Zwanzig, exhibit diffusive motion in the limit where ergodicity is shown to exist. We then explore the diffusion-entropy relation by an accurate calculation of both diffusion and entropy for the aforementioned model systems. We estimate the self-diffusion coefficient of the particle by employing computer simulations and entropy by quadrature using Boltzmann's formula. We observe an interesting crossover in the diffusion-entropy relation in some specific regions which is attributed to the emergence of correlated returns. The crossover could herald a breakdown of the Rosenfeld-like exponential scaling between the two, as observed at low temperatures. Later, we modify the scaling relation to account for the correlated motions and present a detailed analysis of the dynamical entropy obtained via Lyapunov exponent which is rather an important quantity in the study of deterministic systems.


title: Renormalized perturbation theory at large expansion orders

abstract: We present a general formalism that allows for the computation of large-order renormalized expansions in the spacetime representation, effectively doubling the numerically attainable perturbation order of renormalized Feynman diagrams. We show that this formulation compares advantageously to the currently standard techniques due to its high efficiency, simplicity, and broad range of applicability. Our formalism permits to easily complement perturbation theory with non-perturbative information, which we illustrate by implementing expansions renormalized by the addition of a gap or the inclusion of Dynamical Mean-Field Theory. As a result, we present numerically-exact results for the square-lattice Fermi-Hubbard model in the low temperature non-Fermi-liquid regime and show the momentum-dependent suppression of fermionic excitations in the antinodal region.


title: Touchless Typing Using Head Movement-based Gestures

abstract: In this paper, we propose a novel touchless typing interface that makes use of an on-screen QWERTY keyboard and a smartphone camera. The keyboard was divided into nine color-coded clusters. The user moved their head toward clusters, which contained the letters that they wanted to type. A front-facing smartphone camera recorded the head movements. A bidirectional GRU based model which used pre-trained embedding rich in head pose features was employed to translate the recordings into cluster sequences. The model achieved an accuracy of 96.78% and 86.81% under intra- and inter-user scenarios, respectively, over a dataset of 2234 video sequences collected from 22 users.


title: All-Digital Stokes Polarimetry with a Digital Micro-mirror Device

abstract: Stokes polarimetry is widely used to extract the polarisation structure of optical fields, typically from six measurements, although it can be extracted from only four. To measure the required intensities, most approaches are based on optical polarisation components. In this work, we present an all-digital approach that enables a rapid measure of all four intensities without any moving components. Our method employs a Polarisation Grating (PG) to simultaneously project the incoming mode into left- and right-circular polarised states, followed by a polarisation-insensitive Digital Micromirror Device (DMD), which digitally introduces a phase retardance for the acquisition of the remaining two polarisation states. We demonstrate how this technique can be applied to measuring the SoP, vectorness and intra-modal phase of optical fields, without any moving components and shows excellent agreement with theory, illustrating fast, real-time polarimetry.


title: No Routing Needed Between Capsules

abstract: Most capsule network designs rely on traditional matrix multiplication between capsule layers and computationally expensive routing mechanisms to deal with the capsule dimensional entanglement that the matrix multiplication introduces. By using Homogeneous Vector Capsules (HVCs), which use element-wise multiplication rather than matrix multiplication, the dimensions of the capsules remain unentangled. In this work, we study HVCs as applied to the highly structured MNIST dataset in order to produce a direct comparison to the capsule research direction of Geoffrey Hinton, et al. In our study, we show that a simple convolutional neural network using HVCs performs as well as the prior best performing capsule network on MNIST using 5.5x fewer parameters, 4x fewer training epochs, no reconstruction sub-network, and requiring no routing mechanism. The addition of multiple classification branches to the network establishes a new state of the art for the MNIST dataset with an accuracy of 99.87% for an ensemble of these models, as well as establishing a new state of the art for a single model (99.83% accurate).


title: On the boundary local time measure of super-Brownian motion

abstract: If $L^x$ is the total occupation local time of $d$-dimensional super-Brownian motion, $X$, for $d=2$ and $d=3$, we construct a random measure $\mathcal{L}$, called the boundary local time measure, as a rescaling of $L^x e^{-\lambda L^x} dx$ as $\lambda\to \infty$, thus confirming a conjecture of \cite{MP17} and further show that the support of $\mathcal{L}$ equals the topological boundary of the range of $X$, $\partial\mathcal{R}$. This latter result uses a second construction of a boundary local time $\widetilde{\mathcal{L}}$ given in terms of exit measures and we prove that $\widetilde{\mathcal{L}}=c\mathcal{L}$ a.s. for some constant $c>0$. We derive reasonably explicit first and second moment measures for $\mathcal{L}$ in terms of negative dimensional Bessel processes and use it with the energy method to give a more direct proof of the lower bound of the Hausdorff dimension of $\partial\mathcal{R}$ in \cite{HMP18}. The construction requires a refinement of the $L^2$ upper bounds in \cite{MP17} and \cite{HMP18} to exact $L^2$ asymptotics. The methods also refine the left tail bounds for $L^x$ in \cite{MP17} to exact asymptotics. We conjecture that the Minkowski content of $\partial\mathcal{R}$ is equal to the total mass of the boundary local time $\mathcal{L}$ up to some constant.


title: RatLesNetv2: A Fully Convolutional Network for Rodent Brain Lesion Segmentation

abstract: We present a fully convolutional neural network (ConvNet), named RatLesNetv2, for segmenting lesions in rodent magnetic resonance (MR) brain images. RatLesNetv2 architecture resembles an autoencoder and it incorporates residual blocks that facilitate its optimization. RatLesNetv2 is trained end to end on three-dimensional images and it requires no preprocessing. We evaluated RatLesNetv2 on an exceptionally large dataset composed of 916 T2-weighted rat brain MRI scans of 671 rats at nine different lesion stages that were used to study focal cerebral ischemia for drug development. In addition, we compared its performance with three other ConvNets specifically designed for medical image segmentation. RatLesNetv2 obtained similar to higher Dice coefficient values than the other ConvNets and it produced much more realistic and compact segmentations with notably fewer holes and lower Hausdorff distance. The Dice scores of RatLesNetv2 segmentations also exceeded inter-rater agreement of manual segmentations. In conclusion, RatLesNetv2 could be used for automated lesion segmentation, reducing human workload and improving reproducibility. RatLesNetv2 is publicly available at https://github.com/jmlipman/RatLesNetv2.


title: Characteristic classes and stability conditions for projective Kleinian orbisurfaces

abstract: We construct Bridgeland stability conditions on the derived category of smooth quasi-projective Deligne-Mumford surfaces whose coarse moduli spaces have ADE singularities. This unifies the construction for smooth surfaces and Bridgeland's work on Kleinian singularities. The construction hinges on an orbifold version of the Bogomolov-Gieseker inequality for slope semistable sheaves on the stack, and makes use of the To\"en-Hirzebruch-Riemann-Roch theorem.


title: Heating neutron stars with GeV dark matter

abstract: An old neutron star (NS) may capture halo dark matter (DM) and get heated up by the deposited kinetic energy, thus behaving like a thermal DM detector with sensitivity to a wide range of DM masses and a variety of DM-quark interactions. Near future infrared telescopes will measure NS temperatures down to a few thousand Kelvin and probe NS heating by DM capture. We focus on GeV-mass Dirac fermion DM (which is beyond the reach of current DM direct detection experiments) in scenarios in which the DM capture rate can saturate the geometric limit. For concreteness, we study (1) a model that invokes dark decays of the neutron to explain the neutron lifetime anomaly, and (2) a framework of DM coupled to quarks through a vector current portal. In the neutron dark decay model a NS can have a substantial DM population, so that the DM capture rate can reach the geometric limit through DM self-interactions even if the DM-neutron scattering cross section is tiny. We find NS heating to have greater sensitivity than multi-pion signatures in large underground detectors for the neutron dark decay model, and sub-GeV gamma-ray signatures for the quark vector portal model.


title: Aggregation and Data Driven Identification of Building Thermal Dynamic Model and Unmeasured Disturbance

abstract: An aggregate model is a single-zone equivalent of a multi-zone building, and is useful for many purposes, including model based control of large heating, ventilation and air conditioning (HVAC) equipment. This paper deals with the problem of simultaneously identifying an aggregate thermal dynamic model and unknown disturbances from input-output data. The unknown disturbance is a key challenge since it is not measurable but non-negligible. We first present a principled method to aggregate a multi-zone building model into a single zone model, and show the aggregation is not as trivial as it has been assumed in the prior art. We then provide a method to identify the parameters of the model and the unknown disturbance for this aggregate (single-zone) model. Finally, we test our proposed identification algorithm to data collected from a multi-zone building testbed in Oak Ridge National Laboratory. A key insight provided by the aggregation method allows us to recognize under what conditions the estimation of the disturbance signal will be necessarily poor and uncertain, even in the case of a specially designed test in which the disturbances affecting each zone are known (as the case of our experimental testbed). This insight is used to provide a heuristic that can be used to assess when the identification results are likely to have high or low accuracy.


title: Reflex motion in the Milky Way stellar halo resulting from the Large Magellanic Cloud infall

abstract: We present the results of N-body models of the Milky Way and Large Magellanic Cloud system and study the kinematic reflex motion in the stellar halo owing to the barycentre displacement of the disc. In agreement with previous studies, we find that the Milky Way disc may be moving at 40 km/s relative to the barycentre prior to the Large Magellanic Cloud infall. The resulting reflex motion is visible in tangential velocities of the stellar halo as a simple dipole. The signal is strongest for stars with long dynamical times, identifiable in position-velocity data as stars with large apocentres, whose dynamical memory is still well-represented by the unperturbed Milky Way potential. The signal varies across the sky depending on the stellar tracer and may be up to the same magnitude as the velocity of the disc centre-of-mass, making reflex motion a source of bias for Milky Way potential determinations based on the modeling of stellar streams and/or smooth halo tracers such as blue horizontal branch or RR Lyrae stars.


title: On infinite variants of De Morgan law in locale theory

abstract: A locale, being a complete Heyting algebra, satisfies De Morgan law $(a\vee b)^*=a^*\wedge b^*$ for pseudocomplements. The dual De Morgan law $(a\wedge b)^*={a^* \vee b^*}$ (here referred to as the second De Morgan law) is equivalent to, among other conditions, $(a\vee b)^{**} =a^{**}\vee b^{**}$, and characterizes the class of extremally disconnected locales. This paper presents a study of the subclasses of extremally disconnected locales determined by the infinite versions of the second De Morgan law and its equivalents.


title: Sparse Interpolation in Terms of Multivariate Chebyshev Polynomials

abstract: Sparse interpolation} refers to the exact recovery of a function as a short linear combination of basis functions from a limited number of evaluations. For multivariate functions, the case of the monomial basis is well studied, as is now the basis of exponential functions. Beyond the multivariate Chebyshev polynomial obtained as tensor products of univariate Chebyshev polynomials, the theory of root systems allows to define a variety of generalized multivariate Chebyshev polynomials that have connections to topics such as Fourier analysis and representations of Lie algebras. We present a deterministic algorithm to recover a function that is the linear combination of at most r such polynomials from the knowledge of r and an explicitly bounded number of evaluations of this function.


title: The geometric Burge correspondence and the partition function of polymer replicas

abstract: We construct a geometric lifting of the Burge correspondence as a composition of local birational maps on generic Young-diagram-shaped arrays. We establish its fundamental relation to the geometric Robinson-Schensted-Knuth correspondence and to the geometric Sch\"utzenberger involution. We also show a number of properties of the geometric Burge correspondence, specializing them to the case of symmetric input arrays. In particular, our construction shows that such a mapping is volume preserving in log-log variables. As an application, we consider a model of two polymer paths of given length constrained to have the same endpoint, known as polymer replica. We prove that the distribution of the polymer replica partition function in a log-gamma random environment is a Whittaker measure, and deduce the corresponding Whittaker integral identity. For a certain choice of the parameters, we notice a distributional identity between our model and the symmetric log-gamma polymer studied by O'Connell, Sepp\"al\"ainen, and Zygouras (2014).


title: A Combinatorial View of the Service Rates of Codes Problem, its Equivalence to Fractional Matching and its Connection with Batch Codes

abstract: We propose a novel technique for constructing a graph representation of a code through which we establish a significant connection between the service rate problem and the well-known fractional matching problem. Using this connection, we show that the service capacity of a coded storage system equals the fractional matching number in the graph representation of the code, and thus is lower bounded and upper bounded by the matching number and the vertex cover number, respectively. This is of great interest because if the graph representation of a code is bipartite, then the derived upper and lower bounds are equal, and we obtain the capacity. Leveraging this result, we characterize the service capacity of the binary simplex code whose graph representation, as we show, is bipartite. Moreover, we show that the service rate problem can be viewed as a generalization of the multiset primitive batch codes problem.


title: On $\sigma$-arithmetic graphs of finite groups

abstract: Let $G$ be a finite group and $\sigma$ a partition of the set of all? primes $\Bbb{P}$, that is, $\sigma =\{\sigma_i \mid i\in I \}$, where $\Bbb{P}=\bigcup_{i\in I} \sigma_i$ and $\sigma_i\cap \sigma_j= \emptyset $ for all $i\ne j$. If $n$ is an integer, we write $\sigma(n)=\{\sigma_i \mid \sigma_{i}\cap \pi (n)\ne \emptyset \}$ and $\sigma (G)=\sigma (|G|)$. We call a graph $\Gamma$ with the set of all vertices $V(\Gamma)=\sigma (G)$ ($G\ne 1$) a $\sigma$-arithmetic graph of $G$, and we associate with $G\ne 1$ the following three directed $\sigma$-arithmetic graphs: (1) the $\sigma$-Hawkes graph $\Gamma_{H\sigma }(G)$ of $G$ is a $\sigma$-arithmetic graph of $G$ in which $(\sigma_i, \sigma_j)\in E(\Gamma_{H\sigma }(G))$ if $\sigma_j\in \sigma (G/F_{\{\sigma_i\}}(G))$; (2) the $\sigma$-Hall graph $\Gamma_{\sigma Hal}(G)$ of $G$ in which $(\sigma_i, \sigma_j)\in E(\Gamma_{\sigma Hal}(G))$ if for some Hall $\sigma_i$-subgroup $H$ of $G$ we have $\sigma_j\in \sigma (N_{G}(H)/HC_{G}(H))$; (3) the $\sigma$-Vasil'ev-Murashko graph $\Gamma_{{\mathfrak{N}_\sigma }}(G)$ of $G$ in which $(\sigma_i, \sigma_j)\in E(\Gamma_{{\mathfrak{N}_\sigma}}(G))$ if for some ${\mathfrak{N}_{\sigma }}$-critical subgroup $H$ of $G$ we have $\sigma_i \in \sigma (H)$ and $\sigma_j\in \sigma (H/F_{\{\sigma_i\}}(H))$. In this paper, we study the structure of $G$ depending on the properties of these three graphs of $G$.


title: Learning to Catch Security Patches

abstract: Timely patching is paramount to safeguard users and maintainers against dire consequences of malicious attacks. In practice, patching is prioritized following the nature of the code change that is committed in the code repository. When such a change is labeled as being security-relevant, i.e., as fixing a vulnerability, maintainers rapidly spread the change and users are notified about the need to update to a new version of the library or of the application. Unfortunately, oftentimes, some security-relevant changes go unnoticed as they represent silent fixes of vulnerabilities. In this paper, we propose a Co-Training-based approach to catch security patches as part of an automatic monitoring service of code repositories. Leveraging different classes of features, we empirically show that such automation is feasible and can yield a precision of over 90% in identifying security patches, with an unprecedented recall of over 80%. Beyond such a benchmarking with ground truth data which demonstrates an improvement over the state-of-the-art, we confirmed that our approach can help catch security patches that were not reported as such.


title: Signals of a Quantum Universe

abstract: Structure in the Universe is widely believed to have originated from quantum fluctuations during an early epoch of accelerated expansion. Yet, the patterns we observe today do not distinguish between quantum or classical primordial fluctuations; current cosmological data is consistent with either possibility. We argue here that a detection of primordial non-Gaussianity can resolve the present situation, and provide a litmus-test for the quantum origin of cosmic structure. Unlike in quantum mechanics, vacuum fluctuations cannot arise in classical theories and therefore long-range classical correlations must result from (real) particles in the initial state. Similarly to flat-space scattering processes, we show how basic principles require these particles to manifest themselves as poles in the $n$-point functions, in the so-called folded configurations. Following this observation, and assuming fluctuations are (i) correlated over large scales, and (ii) generated by local evolution during an inflationary phase, we demonstrate that: the absence of a pole in the folded limit of non-Gaussian correlators uniquely identifies the quantum vacuum as the initial state. In the same spirit as Bell's inequalities, we discuss how this can be circumvented if locality is abandoned. We also briefly discuss the implications for simulations of a non-Gaussian universe.


title: The Orbit of Planet Nine Derived from Engineering Physics

abstract: Several papers have recently suggested the possible presence of a ninth planet (Planet X) that might explain the gravitational perturbations of a number of detached Trans-Neptunian objects. To analyze the possibility further, we have applied celestial mechanics, engineering physics and statistical analysis to develop improved estimates of the planet's primary orbital elements and mass from first engineering principles, using the orbital characteristics of both the original group of 6 objects analyzed and also a second group comprising the original 6 together with 6 additional long-period asteroids selected by the authors. We show that the driving force behind the observed clustering is gravitational torque that arranges the orbits of asteroids in a systematic, orderly manner, and we develop the associated equations of motion. As evidence we show that the expected effects are fully apparent in the orbital characteristics of the correlated bodies involved, including most strikingly regarding their orbital planes, azimuth orientations and specific relative angular momenta, which we show generates a highly unexpected form of resonance in their relative angular momenta. We further show that the coordinates of Planet 9's orbit are close to the original values proposed recently by other authors, although we prove that its period has to be dramatically smaller than that proposed in recent literature by Batygin and Brown, 2019, at about 3500 yrs, the eccentricity is near 0.65, and its mass approximately 8.4 times the Earth's mass. Given the predicted orbit we show that the planet has apparently created numerous mean motion resonances, of which seven are noted specifically. As for a possible observation, Planet X, would range between V=18.9 and 26.1, probably with a magnitude of about 24.8.


title: Integrated ridesharing services with chance-constrained dynamic pricing and demand learning

abstract: The design of integrated mobility-on-demand services requires jointly considering the interactions between traveler choice behavior and operators' operation policies to design a financially sustainable pricing scheme. However, most existing studies focus on the supply side perspective, disregarding the impact of customer choice behavior in the presence of co-existing transport networks. We propose a modeling framework for dynamic integrated mobility-on-demand service operation policy evaluation with two service options: door-to-door rideshare and rideshare with transit transfer. A new constrained dynamic pricing model is proposed to maximize operator profit, taking into account the correlated structure of different modes of transport. User willingness to pay is considered as a stochastic constraint, resulting in a more realistic ticket price setting while maximizing operator profit. Unlike most studies, which assume that travel demand is known, we propose a demand learning process to calibrate customer demand over time based on customers' historical purchase data. We evaluate the proposed methodology through simulations under different scenarios on a test network by considering the interactions of supply and demand in a multimodal market. Different scenarios in terms of customer arrival intensity, vehicle capacity, and the variance of user willingness to pay are tested. Results suggest that the proposed chance-constrained assortment price optimization model allows increasing operator profit while keeping the proposed ticket prices acceptable.


title: Distant foreground and the Planck-derived Hubble constant

abstract: It is possible to reduce the discrepancy between the local measurement of the cosmological parameter $H_0$ and the value derived from the $Planck$ measurements of the Cosmic Microwave Background (CMB) by considering contamination of the CMB by emission from some medium around distant extragalactic sources, such as extremely cold coarse-grain dust. Though being distant, such a medium would still be in the foreground with respect to the CMB, and, as any other foreground, it would alter the CMB power spectrum. This could contribute to the dispersion of CMB temperature fluctuations. By generating a few random samples of CMB with different dispersions, we have checked that the increased dispersion leads to a smaller estimated value of $H_0$, the rest of the cosmological model parameters remaining fixed. This might explain the reduced value of the $Planck$-derived parameter $H_0$ with respect to the local measurements. The signature of the distant foreground in the CMB traced by SNe was previously reported by the authors of this paper -- we found a correlation between the SN redshifts, $z_{\rm SN}$, and CMB temperature fluctuations at the SNe locations, $T_{\rm SN}$. Here we have used the slopes of the regression lines $T_{\rm SN}\,/\,z_{\rm SN}$ corresponding to different {\it Planck} wave bands in order to estimate the possible temperature of the distant extragalactic medium, which turns out to be very low, about 5\,K. The most likely ingredient of this medium is coarse-grain ($grey$) dust, which is known to be almost undetectable, except for the effect of dimming remote extragalactic sources.


title: Magnetohydrodynamic Simulations of Binary Neutron Star Mergers in General Relativity: Effects of Magnetic Field Orientation on Jet Launching

abstract: Binary neutron star (NSNS) mergers can be sources of gravitational waves coincident with electromagnetic counterpart emission. To solidify their role as multimessenger sources, we present fully 3D, general relativistic, magnetohydrodynamic simulations of spinning NSNSs initially on quasicircular orbits that merge and undergo delayed collapse to a black hole (BH). The NSNSs consist of two identical stars modeled as $\Gamma=2$ polytropes with spin $\chi_{NS}= 0.36$ aligned along the direction of the total orbital angular momentum $L$. Each star is initially threaded by a dynamical unimportant interior dipole B-field. The field is extended into the exterior where a nearly force-free magnetosphere resembles that of a pulsar. The magnetic dipole moment $\mu$ is either aligned or perpendicular to $L$ and has the same initial magnitude for each orientation. For comparison, we also impose symmetry across the orbital plane in one case where $\mu$ in both stars is aligned along $L$. We find that the lifetime of the transient hypermassive neutron star remnant, the jet launching time, and the ejecta are very sensitive to the B-field orientation. By contrast, the physical properties of the BH + disk remnant, such as the mass and spin of the BH, the accretion rate, and the electromagnetic luminosity, are roughly independent of the initial B-field orientation. In addition, we find imposing symmetry across the orbital plane does not play a significant role in the final outcome of the mergers. Our results show that an incipient jet emerges only when the seed B-field has a sufficiently large-scale poloidal component aligned to $L$. The lifetime [$\Delta t\gtrsim 140(M_{NS}/1.625M_\odot)\rm ms$] and Poynting luminosities [$L_{EM}\simeq 10^{52}$erg/s] of the jet, when it forms, are consistent with typical short gamma-ray bursts, as well as with the Blandford--Znajek mechanism for launching jets.


title: The star formation histories of z~1 post-starburst galaxies

abstract: We present the star formation histories of 39 galaxies with high quality rest-frame optical spectra at 0.5<z<1.3 selected to have strong Balmer absorption lines and/or Balmer break, and compare to a sample of spectroscopically selected quiescent galaxies at the same redshift. Photometric selection identifies a majority of objects that have clear evidence for a recent short-lived burst of star formation within the last 1.5 Gyr, i.e. "post-starburst" galaxies, however we show that good quality continuum spectra are required to obtain physical parameters such as burst mass fraction and burst age. Dust attenuation appears to be the primary cause for misidentification of post-starburst galaxies, leading to contamination in spectroscopic samples where only the [OII] emission line is available, as well as a small fraction of objects lost from photometric samples. The 31 confirmed post-starburst galaxies have formed 40-90% of their stellar mass in the last 1-1.5 Gyr. We use the derived star formation histories to find that the post-starburst galaxies are visible photometrically for 0.5-1 Gyr. This allows us to update a previous analysis to suggest that 25-50% of the growth of the red sequence at z~1 could be caused by a starburst followed by rapid quenching. We use the inferred maximum historical star formation rates of several 100-1000 Msun/yr and updated visibility times to confirm that sub-mm galaxies are likely progenitors of post-starburst galaxies. The short quenching timescales of 100-200 Myr are consistent with cosmological hydrodynamic models in which rapid quenching is caused by the mechanical expulsion of gas due to an AGN.


title: Soluble Fermionic Quantum Critical Point in Two Dimensions

abstract: We study a model for a quantum critical point in two spatial dimensions between a semimetallic phase, characterized by a stable quadratic Fermi node, and an ordered phase, in which the spectrum develops a band gap. The quantum critical behavior can be computed exactly, and we explicitly derive the scaling laws of various observables. While the order-parameter correlation function at criticality satisfies the usual power law with anomalous exponent $\eta_\phi = 2$, the correlation length and the expectation value of the order parameter exhibit essential singularities upon approaching the quantum critical point from the insulating side, akin to the Berezinskii-Kosterlitz-Thouless transition. The susceptibility, on the other hand, has a power-law divergence with non-mean-field exponent $\gamma = 2$. On the semimetallic side, the correlation length remains infinite, leading to an emergent scale invariance throughout this phase.


title: Impact of uncertainties in the halo velocity profile on direct detection of sub-GeV dark matter

abstract: We use the state-of-the-art high-resolution cosmological simulations by IllustrisTNG to derive the velocity distribution and local density of dark matter in galaxies like our Milky Way and find a substantial spread in both quantities. Next we use our findings to examine the sensitivity to the dark matter velocity profile of underground searches using electron scattering in germanium and silicon targets. We find that sub-GeV dark matter search is strongly affected by these uncertainties, unlike nuclear recoil searches for heavier dark matter, especially in multiple electron-hole modes, for which the sensitivity to the scattering cross-section is also weaker. Therefore, by improving the sensitivity to lower ionization thresholds not only projected sensitivities will be boosted but also the dependence on the astrophysical uncertainties will become significantly reduced.


title: Mass-to-Light Ratios of Spatially Resolved Stellar Populations in M31

abstract: A galaxy's stellar mass-to-light ratio ($M_\star/L$) is a useful tool for converting luminosity to stellar mass ($M_\star$). However, the practical utility of $M_\star/L$ inferred from stellar population synthesis (SPS) models is limited by mismatches between the real and assumed models for star formation history (SFH) and dust geometry, both of which vary within galaxies. Here, we measure spatial variations in $M_\star/L$ and their dependence on color, star formation history, and dust across the disk of M31, using a map of $M^\mathrm{CMD}_\star$ derived from color-magnitude diagrams of resolved stars in the Panchromatic Hubble Andromeda Treasury (PHAT) survey. First, we find comparable scatter in $M_\star/L$ for the optical and mid-IR, contrary to the common idea that $M_\star/L$ is less variable in the IR. Second, we confirm that $M_\star/L$ is correlated with color for both the optical and mid-IR and report color vs. $M_\star/L$ relations (CMLRs) in M31 for filters used in the Sloan Digital Sky Survey (SDSS) and Widefield Infrared Survey Explorer (WISE). Third, we show that the CMLR residuals correlate with recent SFH, such that quiescent regions are offset to higher $M_\star/L$ than star-forming regions at a fixed color. The mid-IR CMLR, however, is not linear due to the high scatter of $M_\star/L$ in star-forming regions. Finally, we find a flatter optical CMLR than any SPS-based CMLRs in the literature. We show this is an effect of dust geometry, which is typically neglected but should be accounted for when using optical data to map $M_\star/L$.


title: Curvature Perturbations From Stochastic Particle Production During Inflation

abstract: We calculate the curvature power spectrum sourced by spectator fields that are excited repeatedly and non-adiabatically during inflation. In the absence of detailed information of the nature of spectator field interactions, we consider an ensemble of models with intervals between the repeated interactions and interaction strengths drawn from simple probabilistic distributions. We show that the curvature power spectra of each member of the ensemble shows rich structure with many features, and there is a large variability between different realizations of the same ensemble. Such features can be probed by the cosmic microwave background (CMB) and large scale structure observations. They can also have implications for primordial black hole formation and CMB spectral distortions. The geometric random walk behavior of the spectator field allows us to calculate the ensemble-averaged power spectrum of curvature perturbations semi-analytically. For sufficiently large stochastic sourcing, the ensemble-averaged power spectrum shows a scale dependence arising from the time spent by modes outside the horizon during the period of particle production, in spite of there being no preferred scale in the underlying model. We find that the magnitude of the ensemble-averaged power spectrum overestimates the typical power spectra in the ensemble because the ensemble distribution of the power spectra is highly non-Gaussian with fat tails.


title: From the pseudogap metal to the Fermi liquid using ancilla qubits

abstract: We propose a new parton theory of the hole-doped cuprates, describing the evolution from the pseudogap metal with small Fermi surfaces to the conventional Fermi liquid with a large Fermi surface. We introduce 2 ancilla qubits per square lattice site, and employ them to obtain a variational wavefunction of a fractionalized Fermi liquid for the pseudogap metal state. We propose a multi-layer Hamiltonion for the cuprates, with the electrons residing in the 'physical' layer, and the ancilla qubits in two 'hidden' layers: the hidden layers can be decoupled from the physical layer by a canonical transformation which leaves the hidden layers in a trivial gapped state. This Hamiltonian yields an emergent gauge theory which describes not only the fractionalized Fermi liquid, but also the conventional Fermi liquid, and possible exotic intermediate phases and critical points. The fractionalized Fermi liquid has hole pockets with quasiparticle weight which is large only on "Fermi arcs", and fermionic spinon excitations which carry charges of the emergent gauge fields.


title: Exploring Primordial Black Holes from the Multiverse with Optical Telescopes

abstract: Primordial black holes (PBHs) are a viable candidate for dark matter if the PBH masses are in the currently unconstrained "sublunar" mass range. We revisit the possibility that PBHs were produced by nucleation of false vacuum bubbles during inflation. We show that this scenario can produce a population of PBHs that simultaneously accounts for all dark matter, explains the candidate event in Subaru Hyper Suprime-Cam (HSC) data, and contains both heavy black holes as observed by LIGO and very heavy seeds of supermassive black holes. We demonstrate with numerical studies that future observations of HSC, as well as other optical surveys, such as LSST, will be able to provide a definitive test for this generic PBH formation mechanism if it is the dominant source of dark matter.


title: Self-testing of a single quantum device under computational assumptions

abstract: Self-testing is a method to characterise an arbitrary quantum system based only on its classical input-output correlations, and plays an important role in device-independent quantum information processing as well as quantum complexity theory. Prior works on self-testing require the assumption that the system's state is shared among multiple parties that only perform local measurements and cannot communicate. Here, we replace the setting of multiple non-communicating parties, which is difficult to enforce in practice, by a single computationally bounded party. Specifically, we construct a protocol that allows a classical verifier to robustly certify that a single computationally bounded quantum device must have prepared a Bell pair and performed single-qubit measurements on it, up to a change of basis applied to both the device's state and measurements. This means that under computational assumptions, the verifier is able to certify the presence of entanglement, a property usually closely associated with two separated subsystems, inside a single quantum device. To achieve this, we build on techniques first introduced by Brakerski et al. (2018) and Mahadev (2018) which allow a classical verifier to constrain the actions of a quantum device assuming the device does not break post-quantum cryptography.


title: Low Mach and thin domain limit for the compressible Euler system

abstract: We consider the compressible Euler system describing the motion of an ideal fluid confined to a straight layer $\Omega_{\delta}=(0,\delta)\times\mathbb{R}^2, \ \ \delta>0$. In the framework of dissipative measure-valued solutions, we show the convergence to the strong solution of the 2D incompressible Euler system when the Mach number tends to zero and $\delta\rightarrow0$.


title: Highly-Chirped Bragg Gratings for Integrated Silica Spectrometers

abstract: A blazed chirped Bragg grating in a planar silica waveguide device was used to create an integrated diffractive element for a spectrometer. The grating diffracts light from a waveguide and creates a wavelength dependent focus in a manner similar to a bulk diffraction grating spectrometer. An external imaging system is used to analyse the light, later device iterations plan to integrate detectors to make a fully integrated spectrometer. Devices were fabricated with grating period chirp rates in excess of 100nm/mm, achieving a focal length of 5.5mm. Correction of coma aberrations resulted in a device with a footprint of 20mm x 10mm, a peak FWHM resolution of 1.8nm, a typical FWHM resolution of 2.6nm and operating with a 160nm bandwidth centered at 1550nm.


title: On-shell heavy particle effective theories

abstract: We introduce on-shell variables for Heavy Particle Effective Theories (HPETs) with the goal of extending Heavy Black Hole Effective Theory to higher spins and of facilitating its application to higher post-Minkowskian orders. These variables inherit the separation of spinless and spin-inclusive effects from the HPET fields, resulting in an explicit spin-multipole expansion of the three-point amplitude for any spin. By matching amplitudes expressed using the on-shell HPET variables to those derived from the one-particle effective action, we find that the spin-multipole expansion of a heavy spin-$s$ particle corresponds exactly to the multipole expansion (up to order $2s$) of a Kerr black hole, that is, without needing to take the infinite spin limit. Finally, we show that tree-level radiative processes with same-helicity bosons emitted from a heavy spin-$s$ particle exhibit a spin-multipole universality.


title: Collective vibrations of confined levitating droplets

abstract: We report a new type of fluid-based driven dissipative oscillator system consisting of a lattice of millimetric fluid droplets bouncing on a vertically vibrating liquid bath and bound within an annular ring. We characterize the system behavior as it is energized through a progressive increase in the bath's vibrational acceleration. Depending on the number of drops, the onset of motion of the lattice may take the form of either out-of-phase oscillations or a striking solitary wave-like instability. Theoretical modeling demonstrates that these behaviors may be attributed to different bifurcations at the onset of instability. The results presented here demonstrate the potential and utility of the walking droplet system as a platform for investigating wave-mediated, inertial, non-equilibrium particle dynamics at the macroscale.


title: Reducing bias in dual flip angle T1-mapping in human brain at 7T

abstract: Purpose: To address the systematic bias in whole-brain dual flip angle (DFA) T1-mapping at 7T by optimizing the flip angle pair and carefully selecting RF pulse shape and duration. Theory and Methods: Spoiled gradient echoes can be used to estimate whole-brain maps of T1. This can be accomplished by using only two acquisitions with different flip angles, i.e., a DFA-based approach. Although DFA-based T1-mapping is seemingly straightforward to implement, it is sensitive to bias caused by incomplete spoiling and incidental magnetization transfer (MT) effects. Further bias is introduced by the increased B0 and B1+ inhomogeneities at 7T. Experiments were performed to determine the optimal flip angle pair and appropriate RF pulse shape and duration. Obtained T1 estimates were validated using inversion recovery prepared EPI and compared to literature values. A multi-echo readout was used to increase SNR, enabling quantification of R2* and susceptibility, X. Results: Incomplete spoiling was observed above a local flip angle of approximately 20 degrees. An asymmetric gauss-filtered sinc pulse with a constant duration of 700 us showed a sufficiently flat frequency response profile to avoid incomplete excitation in areas with high B0 offsets. A pulse duration of 700 us minimized effects from incidental MT. Conclusion: When performing DFA-based T1-mapping one should (i) limit the higher flip angle to avoid incomplete spoiling, (ii) use a RF pulse shape insensitive to B0 inhomogeneities and (iii) apply a constant RF pulse duration, balanced to minimize incidental MT.


title: Subdirect products and propagating equations with an application to Moufang theorem

abstract: We introduce the concept of propagating equations and focus on the case of associativity propagating in varieties of loops. An equation $\varepsilon$ propagates in an algebra $X$ if $\varepsilon(\overrightarrow y)$ holds whenever $\varepsilon(\overrightarrow x)$ holds and the elements of $\overrightarrow y$ are contained in the subalgebra of $X$ generated by $\overrightarrow x$. If $\varepsilon$ propagates in $X$ then it propagates in all subalgebras and products of $X$ but not necessarily in all homomorphic images of $X$. If $\mathcal V$ is a variety, the propagating core $\mathcal V_{[\varepsilon]} = \{X\in\mathcal V:\varepsilon$ propagates in $X\}$ is a quasivariety but not necessarily a variety. We prove by elementary means Goursat's Lemma for loops and describe all subdirect products of $X^k$ and all finitely generated loops in $\mathbf{HSP}(X)$ for a nonabelian simple loop $X$. If $\mathcal V$ is a variety of loops in which associativity propagates, $X$ is a finite loop in which associativity propagates and every subloop of $X$ is either nonabelian simple or contained in $\mathcal V$, then associativity propagates in $\mathbf{HSP}(X)\lor\mathcal V$. We study the propagating core $\mathcal S_{[x(yz)=(xy)z]}$ of Steiner loops with respect to associativity. While this is not a variety, we exhibit many varieties contained in $\mathcal S_{[x(yz)=(xy)z]}$, each providing a solution to Rajah's problem, i.e., a variety of loops not contained in Moufang loops in which Moufang Theorem holds.


title: The Threshold Dimension of a Graph

abstract: Let $G$ be a graph, and let $u$, $v$, and $w$ be vertices of $G$. If the distance between $u$ and $w$ does not equal the distance between $v$ and $w$, then $w$ is said to resolve $u$ and $v$. The metric dimension of $G$, denoted $\beta(G)$, is the cardinality of a smallest set $W$ of vertices such that every pair of vertices of $G$ is resolved by some vertex of $W$. The threshold dimension of a graph $G$, denoted $\tau(G)$, is the minimum metric dimension among all graphs $H$ having $G$ as a spanning subgraph. In other words, the threshold dimension of $G$ is the minimum metric dimension among all graphs obtained from $G$ by adding edges. If $\beta(G) = \tau(G)$, then $G$ is said to be \emph{irreducible}; otherwise, we say that $G$ is reducible. If $H$ is a graph having $G$ as a spanning subgraph and such that $\beta(H)=\tau(G)$, then $H$ is called a threshold graph of $G$. The threshold dimension of a graph is expressed in terms of a minimum number of strong products of paths that admits a certain type of embedding of the graph. A sharp upper bound for the threshold dimension of trees is established. It is also shown that the irreducible trees are precisely those of metric dimension at most 2. Moreover, if $T$ is a tree with metric dimension 3 or 4, then $T$ has threshold dimension $2$. It is shown, in these two cases, that a threshold graph for $T$ can be obtained by adding exactly one or two edges to $T$, respectively. However, these results do not extend to trees with metric dimension $5$, i.e., there are trees of metric dimension $5$ with threshold dimension exceeding $2$.


title: Ergodic-localized junctions in a periodically-driven spin chain

abstract: We report the analogue simulation of an ergodiclocalized junction by using an array of 12 coupled superconducting qubits. To perform the simulation, we fabricated a superconducting quantum processor that is divided into two domains: a driven domain representing an ergodic system, while the second is localized under the effect of disorder. Due to the overlap between localized and delocalized states, for small disorder there is a proximity effect and localization is destroyed. To experimentally investigate this, we prepare a microwave excitation in the driven domain and explore how deep it can penetrate the disordered region by probing its dynamics. Furthermore, we performed an ensemble average over 50 realizations of disorder, which clearly shows the proximity effect. Our work opens a new avenue to build quantum simulators of driven-disordered systems with applications in condensed matter physics and material science

categories: quant-ph cond-mat.dis-nn cond-mat.stat-mech cond-mat.str-el


title: Software-Defined Location Privacy Protection for Vehicular Networks

abstract: While the adoption of connected vehicles is growing, security and privacy concerns are still the key barriers raised by society. These concerns mandate automakers and standardization groups to propose convenient solutions for privacy preservation. One of the main proposed solutions is the use of Pseudonym-Changing Strategies (PCSs). However, ETSI has recently published a technical report which highlights the absence of standardized and efficient PCSs [1]. This alarming situation mandates an innovative shift in the way that the privacy of end-users is protected during their journey. Software-Defined Networking (SDN) is emerging as a key 5G enabler to manage the network in a dynamic manner. SDN-enabled wireless networks are opening up new programmable and highly-flexible privacy-aware solutions. We exploit this paradigm to propose an innovative software-defined location privacy architecture for vehicular networks. The proposed architecture is context-aware, programmable, extensible, and able to encompass all existing and future pseudonym-changing strategies. To demonstrate the merit of our architecture, we consider a case study that involves four pseudonym-changing strategies, which we deploy over our architecture and compare with their static implementations. We also detail how the SDN controller dynamically switches between the strategies according to the context.


title: Database development and exploration of microstructure versus process relationships using variational autoencoders

abstract: The paper demonstrates the use of variational autoencoders for graphical representation of a large database containing process-microstructure relationships. Correlating microstructural features to processing is an essential first step to answer the difficult problem of process sequence design. In this paper, a large database of 346,200 orientation distribution functions resulting from a variety of process sequences is constructed, where each sequence comprises up to four stages of tension, compression and rolling along different directions in various permutations. This opensource database is constructed for collaborative development of process design algorithms. The paper demonstrates a novel application of the large database: graphical representation of texture-process relationships. A variational autoencoder is used to reduce the entire database to a two dimensional latent space where variations in processes and properties can be visualized. Using proximity analysis in this latent space, we can quickly unearth multiple process solutions to the problem of texture or property design.


title: Reflexion maps and geometry of surfaces in R^4

abstract: In this article we introduce new affinely invariant points---`special parabolic points'---on the parabolic set of a generic surface $M$ in real 4-space, associated with symmetries in the 2-parameter family of reflexions of $M$ in points of itself. The parabolic set itself is detected in this way, and each arc is given a sign, which changes at the special points, where the family has an additional degree of symmetry. Other points of $M$ which are detected by the family of reflexions include inflexion points of real and imaginary type, and the first of these is also associated with sign changes on the parabolic set. We show how to compute the special points globally for the case where $M$ is given in Monge form and give some examples illustrating the birth of special parbolic points in a 1-parameter family of surfaces. The tool we use from singularity theory is the contact classification of certain symmetric maps from the plane to the plane and we give the beginning of this classification, including versal unfoldings which we relate to the geometry of $M$.


title: Automatic MRI-Driven Model Calibration for Advanced Brain Tumor Progression Analysis

abstract: Our objective is the calibration of mathematical tumor growth models from a single multiparametric scan. The target problem is the analysis of preoperative Glioblastoma (GBM) scans. To this end, we present a fully automatic tumor-growth calibration methodology that integrates a single-species reaction-diffusion partial differential equation (PDE) model for tumor progression with multiparametric Magnetic Resonance Imaging (mpMRI) scans to robustly extract patient specific biomarkers i.e., estimates for (i) the tumor cell proliferation rate, (ii) the tumor cell migration rate, and (iii) the original, localized site(s) of tumor initiation. Our method is based on a sparse reconstruction algorithm for the tumor initial location (TIL). This problem is particularly challenging due to nonlinearity, ill-posedeness, and ill conditioning. We propose a coarse-to-fine multi-resolution continuation scheme with parameter decomposition to stabilize the inversion. We demonstrate robustness and practicality of our method by applying the proposed method to clinical data of 206 GBM patients. We analyze the extracted biomarkers and relate tumor origin with patient overall survival by mapping the former into a common atlas space. We present preliminary results that suggest improved accuracy for prediction of patient overall survival when a set of imaging features is augmented with estimated biophysical parameters. All extracted features, tumor initial positions, and biophysical growth parameters are made publicly available for further analysis. To our knowledge, this is the first fully automatic scheme that can handle multifocal tumors and can localize the TIL to a few millimeters.


title: Weakly Supervised Lesion Co-segmentation on CT Scans

abstract: Lesion segmentation in medical imaging serves as an effective tool for assessing tumor sizes and monitoring changes in growth. However, not only is manual lesion segmentation time-consuming, but it is also expensive and requires expert radiologist knowledge. Therefore many hospitals rely on a loose substitute called response evaluation criteria in solid tumors (RECIST). Although these annotations are far from precise, they are widely used throughout hospitals and are found in their picture archiving and communication systems (PACS). Therefore, these annotations have the potential to serve as a robust yet challenging means of weak supervision for training full lesion segmentation models. In this work, we propose a weakly-supervised co-segmentation model that first generates pseudo-masks from the RECIST slices and uses these as training labels for an attention-based convolutional neural network capable of segmenting common lesions from a pair of CT scans. To validate and test the model, we utilize the DeepLesion dataset, an extensive CT-scan lesion dataset that contains 32,735 PACS bookmarked images. Extensive experimental results demonstrate the efficacy of our co-segmentation approach for lesion segmentation with a mean Dice coefficient of 90.3%.


title: A hot terrestrial planet orbiting the bright M dwarf L 168-9 unveiled by TESS

abstract: We report the detection of a transiting super-Earth-sized planet (R=1.39+-0.09 Rearth) in a 1.4-day orbit around L 168-9 (TOI-134),a bright M1V dwarf (V=11, K=7.1) located at 25.15+-0.02 pc. The host star was observed in the first sector of the Transiting Exoplanet Survey Satellite (TESS) mission and, for confirmation and planet mass measurement, was followed up with ground-based photometry, seeing-limited and high-resolution imaging, and precise radial velocity (PRV) observations using the HARPS and PFS spectrographs. Combining the TESS data and PRV observations, we find the mass of L168-9 b to be 4.60+-0.56 Mearth, and thus the bulk density to be 1.74+0.44-0.33 times larger than that of the Earth. The orbital eccentricity is smaller than 0.21 (95% confidence). This planet is a Level One Candidate for the TESS Mission's scientific objective - to measure the masses of 50 small planets - and is one of the most observationally accessible terrestrial planets for future atmospheric characterization.


title: Some interpretations for algebraic invariants of edge ideals of hypergraphs via combinatorial invariants

abstract: The present work is concerned with characterizing some algebraic invariants of edge ideals of hypergraphs. To this aim, firstly, we introduce some kinds of combinatorial invariants similar to matching numbers for hypergraphs. Then we compare them to each other and to previously existing ones. These invariants are used for characterizing or bounding some algebraic invariants of edge ideals of hypergraphs such as graded Betti numbers, projective dimension and Castelnouvo-Mumford regularity.


title: Recognizing Developers' Emotions while Programming

abstract: Developers experience a wide range of emotions during programming tasks, which may have an impact on job performance. In this paper, we present an empirical study aimed at (i) investigating the link between emotion and progress, (ii) understanding the triggers for developers' emotions and the strategies to deal with negative ones, (iii) identifying the minimal set of non-invasive biometric sensors for emotion recognition during programming task. Results confirm previous findings about the relation between emotions and perceived productivity. Furthermore, we show that developers' emotions can be reliably recognized using only a wristband capturing the electrodermal activity and heart-related metrics.


title: Analyticity of the percolation density $\theta$ in all dimensions

abstract: We prove that for Bernoulli bond percolation on $\mathbb{Z}^d$, $d\geq 2$ the percolation density is an analytic function of the parameter in the supercritical interval $(p_c,1]$. This answers a question of Kesten from 1981.


title: Correlated and in-situ electrical transmission electron microscopy studies and related membrane fabrication

abstract: Understanding the interplay between the structure, composition and opto-electronic properties of semiconductor nano-objects requires combining transmission electron microscopy (TEM) based techniques with electrical and optical measurements on the very same specimen. Recent developments in TEM technologies allow not only the identification and in-situ electrical characterization of a particular object, but also the direct visualization of its modification in-situ by techniques such as Joule heating. Over the past years, we have carried out a number of studies in these fields that are reviewed in this contribution. In particular, we discuss here i) correlated studies where the same unique object is characterized electro-optically and by TEM, ii) in-situ Joule heating studies where a solid-state metal-semiconductor reaction is monitored in the TEM, and iii) in-situ biasing studies to better understand the electrical properties of contacted single nanowires. In addition, we provide detailed fabrication steps for the silicon nitride membranes crucial to these correlated and in-situ measurements.


title: Imputation for High-Dimensional Linear Regression

abstract: We study high-dimensional regression with missing entries in the covariates. A common strategy in practice is to \emph{impute} the missing entries with an appropriate substitute and then implement a standard statistical procedure acting as if the covariates were fully observed. Recent literature on this subject proposes instead to design a specific, often complicated or non-convex, algorithm tailored to the case of missing covariates. We investigate a simpler approach where we fill-in the missing entries with their conditional mean given the observed covariates. We show that this imputation scheme coupled with standard off-the-shelf procedures such as the LASSO and square-root LASSO retains the minimax estimation rate in the random-design setting where the covariates are i.i.d.\ sub-Gaussian. We further show that the square-root LASSO remains \emph{pivotal} in this setting. It is often the case that the conditional expectation cannot be computed exactly and must be approximated from data. We study two cases where the covariates either follow an autoregressive (AR) process, or are jointly Gaussian with sparse precision matrix. We propose tractable estimators for the conditional expectation and then perform linear regression via LASSO, and show similar estimation rates in both cases. We complement our theoretical results with simulations on synthetic and semi-synthetic examples, illustrating not only the sharpness of our bounds, but also the broader utility of this strategy beyond our theoretical assumptions.


title: End-to-End Vision-Based Adaptive Cruise Control (ACC) Using Deep Reinforcement Learning

abstract: This paper presented a deep reinforcement learning method named Double Deep Q-networks to design an end-to-end vision-based adaptive cruise control (ACC) system. A simulation environment of a highway scene was set up in Unity, which is a game engine that provided both physical models of vehicles and feature data for training and testing. Well-designed reward functions associated with the following distance and throttle/brake force were implemented in the reinforcement learning model for both internal combustion engine (ICE) vehicles and electric vehicles (EV) to perform adaptive cruise control. The gap statistics and total energy consumption are evaluated for different vehicle types to explore the relationship between reward functions and powertrain characteristics. Compared with the traditional radar-based ACC systems or human-in-the-loop simulation, the proposed vision-based ACC system can generate either a better gap regulated trajectory or a smoother speed trajectory depending on the preset reward function. The proposed system can be well adaptive to different speed trajectories of the preceding vehicle and operated in real-time.


title: iGLU 2.0: A new non-invasive, accurate serum glucometer for smart healthcare

abstract: To best of the authors knowledge, this article presents the first-ever non-invasive glucometer that takes into account serum glucose for high accuracy. In case of blood glucose measurement, serum glucose value has always been considered precise blood glucose value during prandial modes. Serum glucose can be measured in laboratory and more stable glucose level compare to capillary glucose. However, this invasive approach is not convenient for frequent measurement. Sometimes, Conventional invasive blood glucose measurement may be responsible for cause of trauma and chance of blood related infections. To overcome this issue, in the current paper, we propose a novel Internet-of-Medical (IoMT) enabled glucometer for non-invasive precise serum glucose measurement. In this work, a near-infrared (NIR) spectroscopic technique has been used for glucose measurement. The novel device called iGLU 2.0 is based on optical detection and precise machine learning (ML) regression models. The optimal multiple polynomial regression and deep neural network models have been presented to analyze the precise measurement. The glucose values of serum are saved on cloud through open IoT platform for endocrinologist at remote location. To validate iGLU 2.0, Mean Absolute Relative Difference (mARD) and Average Error (AvgE) are obtained 6.07% and 6.09%, respectively from predicted blood glucose values for capillary glucose. For serum glucose, mARD and AvgE are found 4.86% and 4.88%, respectively. These results represent that the proposed non-invasive glucose measurement device is more precise for serum glucose compared to capillary glucose.


title: An Overview of Fingerprint-Based Authentication: Liveness Detection and Beyond

abstract: In this paper, we provide an overview of fingerprint sensing methods used for authentication. We analyze the current fingerprint sensing technologies, from algorithmic, as well as from hardware perspectives. We then focus on methods to detect physical liveness, defined as techniques that can be used to ensure that a living human user is attempting to authenticate on a system. We analyze how effective these methods are at preventing attacks where a malicious entity tries to trick a fingerprint-based authentication system to accept a fake finger as a real one (spoofing attacks). We then identify broader attack points against biometric data, such as fingerprints. Finally, we propose novel measures to protect fingerprint data.


title: Characterization of microwave absorption in carbon nanotubes using resonance aperture transmission method

abstract: A new method to characterize microwave electromagnetic absorption of a bulk carbon nanotube material is proposed and experimentally evaluated in this paper. The method is based on the measurement of microwave transmission through a capacitive-resonator aperture in a conductive screen loaded with a CNT sample under test. This method allows to measure microwave permittivity and absorption of thin samples, several atomic layers to few micrometers thick, with linear dimensions much smaller than the wavelength of radiation in free space. This minimal sample requirement restricts the application of conventional microwave characterization methods such as free-space or waveguide permittivity characterization. It is demonstrated that the resonance E-field enhancement inside the CRA leads to strong EM interaction of the microwave E-field with the CNT sample under test thus enabling high sensitivity and dynamic range of the measurement procedure. Another advantage of the proposed technique over conventional non-resonance characterization methods is that in the resonance transmission band, the CRA operation is reflection-less which leads to a relatively simple qualitative algebraic de-embedding procedure of the material parameters based on the principle of energy conservation. The experimental microwave absorption data of the multiwall CNT samples are presented in the S frequency band (2-4GHz), demonstrating microwave absorption properties of the multiwall CNT ribbons.


title: On the dimension of groups that satisfy certain conditions on their finite subgroups

abstract: We prove that the Bredon cohomological dimension and the virtual cohomological dimension coincide for groups that admit a cocompact model for $\underline{E}G$ and satisfy properties (M) and (NM). Among the examples of groups satisfying these hypothesis are cocompact and arithmetic Fuchsian groups, one relator groups, the Hilbert modular group and $3$-manifold groups.


title: A tutorial on the range variant of asymmetric numeral systems

abstract: This paper is intended to be a brief and accessible introduction to the range variant of asymmetric numeral systems (ANS), a system for lossless compression of sequences which can be used as a drop in replacement for arithmetic coding (AC). Because of the relative simplicity of ANS, we are able to provide enough mathematical detail to rigorously prove that ANS attains a compression rate close to the Shannon limit. Pseudo-code, intuitive interpretation and diagrams are given alongside the mathematical derivations. A working Python demo which accompanies this tutorial is available at https://raw.githubusercontent.com/j-towns/ans-notes/master/rans.py.


title: Certified and fast computations with shallow covariance kernels

abstract: Many techniques for data science and uncertainty quantification demand efficient tools to handle Gaussian random fields, which are defined in terms of their mean functions and covariance operators. Recently, parameterized Gaussian random fields have gained increased attention, due to their higher degree of flexibility. However, especially if the random field is parameterized through its covariance operator, classical random field discretization techniques fail or become inefficient. In this work we introduce and analyze a new and certified algorithm for the low-rank approximation of a parameterized family of covariance operators which represents an extension of the adaptive cross approximation method for symmetric positive definite matrices. The algorithm relies on an affine linear expansion of the covariance operator with respect to the parameters, which needs to be computed in a preprocessing step using, e.g., the empirical interpolation method. We discuss and test our new approach for isotropic covariance kernels, such as Mat\'ern kernels. The numerical results demonstrate the advantages of our approach in terms of computational time and confirm that the proposed algorithm provides the basis of a fast sampling procedure for parameter dependent Gaussian random fields.


title: Ensemble Rejection Sampling

abstract: We introduce Ensemble Rejection Sampling, a scheme for exact simulation from the posterior distribution of the latent states of a class of non-linear non-Gaussian state-space models. Ensemble Rejection Sampling relies on a proposal for the high-dimensional state sequence built using ensembles of state samples. Although this algorithm can be interpreted as a rejection sampling scheme acting on an extended space, we show under regularity conditions that the expected computational cost to obtain an exact sample increases cubically with the length of the state sequence instead of exponentially for standard rejection sampling. We demonstrate this methodology by sampling exactly state sequences according to the posterior distribution of a stochastic volatility model and a non-linear autoregressive process. We also present an application to rare event simulation.


title: Learning a distance function with a Siamese network to localize anomalies in videos

abstract: This work introduces a new approach to localize anomalies in surveillance video. The main novelty is the idea of using a Siamese convolutional neural network (CNN) to learn a distance function between a pair of video patches (spatio-temporal regions of video). The learned distance function, which is not specific to the target video, is used to measure the distance between each video patch in the testing video and the video patches found in normal training video. If a testing video patch is not similar to any normal video patch then it must be anomalous. We compare our approach to previously published algorithms using 4 evaluation measures and 3 challenging target benchmark datasets. Experiments show that our approach either surpasses or performs comparably to current state-of-the-art methods.


title: Impact of ionizing radiation on superconducting qubit coherence

abstract: The practical viability of any qubit technology stands on long coherence times and high-fidelity operations, with the superconducting qubit modality being a leading example. However, superconducting qubit coherence is impacted by broken Cooper pairs, referred to as quasiparticles, with a density that is empirically observed to be orders of magnitude greater than the value predicted for thermal equilibrium by the Bardeen-Cooper-Schrieffer (BCS) theory of superconductivity. Previous work has shown that infrared photons significantly increase the quasiparticle density, yet even in the best isolated systems, it still remains higher than expected, suggesting that another generation mechanism exists. In this Letter, we provide evidence that ionizing radiation from environmental radioactive materials and cosmic rays contributes to this observed difference, leading to an elevated quasiparticle density that would ultimately limit superconducting qubits of the type measured here to coherence times in the millisecond regime. We further demonstrate that introducing radiation shielding reduces the flux of ionizing radiation and positively correlates with increased coherence time. Albeit a small effect for today's qubits, reducing or otherwise mitigating the impact of ionizing radiation will be critical for realizing fault-tolerant superconducting quantum computers.


title: Estimation of Building Rooftop Temperature from High Spatial Resolution Aerial Thermal Images

abstract: This letter presents a novel technique to calculate temperatures of building rooftops and other impervious surfaces from high spatial resolution aerial thermal images. In this study, we collected aerial radiance images of 30cm spatial resolution using a FLIR Phoenix imager in long-wave and mid-wave infrared wavelengths for the city of Cedar Falls, USA to estimate building roof temperature loss. Simultaneous ground temperature measurements were made at pre-selected ground targets and roofs using 9 Fluke 561r infrared thermometers. Atmospheric correction of aerial images was performed by Empirical Line Calibration (ELC) method. The resulting ground-leaving radiances were corrected for the emissivity of different roof types and the true kinetic temperature of the building roofs was calculated. The ELC model was observed to perform better when only impervious surface targets were used for the regression. With an R2=0.71 for ELC, the method produced a root mean squared error of 0.74{\deg}C for asphalt roofs. Further, we observed that the microclimate plays a significant role while synchronizing aerial and ground measurements.


title: Effects of Strain and Film Thickness on the Stability of the Rhombohedral Phase of HfO2

abstract: The discovery of ferroelectric polarization in HfO2-based ultrathin films has spawned a lot of interest due to their potential applications in data storage. Recently, a new R3m rhombohedral phase was proposed to be responsible for the emergence of ferroelectricity in the [111]-oriented Hf0.5Zr0.5O2 thin films, but the fundamental mechanism of ferroelectric polarization in such films remains poorly understood. In this paper, we employ density-functional-theory calculations to investigate structural and polarization properties of the R3m HfO2 phase. We find that the film thickness and in-plane compressive strain effects play a key role in stabilizing the R3m phase leading to robust ferroelectricity of [111]-oriented R3m HfO2.


title: VerSe: A Vertebrae Labelling and Segmentation Benchmark for Multi-detector CT Images

abstract: Vertebral labelling and segmentation are two fundamental tasks in an automated spine processing pipeline. Reliable and accurate processing of spine images is expected to benefit clinical decision-support systems for diagnosis, surgery planning, and population-based analysis on spine and bone health. However, designing automated algorithms for spine processing is challenging predominantly due to considerable variations in anatomy and acquisition protocols and due to a severe shortage of publicly available data. Addressing these limitations, the Large Scale Vertebrae Segmentation Challenge (VerSe) was organised in conjunction with the International Conference on Medical Image Computing and Computer Assisted Intervention (MICCAI) in 2019 and 2020, with a call for algorithms towards labelling and segmentation of vertebrae. Two datasets containing a total of 374 multi-detector CT scans from 355 patients were prepared and 4505 vertebrae have individually been annotated at voxel-level by a human-machine hybrid algorithm (https://osf.io/nqjyw/, https://osf.io/t98fz/). A total of 25 algorithms were benchmarked on these datasets. In this work, we present the the results of this evaluation and further investigate the performance-variation at vertebra-level, scan-level, and at different fields-of-view. We also evaluate the generalisability of the approaches to an implicit domain shift in data by evaluating the top performing algorithms of one challenge iteration on data from the other iteration. The principal takeaway from VerSe: the performance of an algorithm in labelling and segmenting a spine scan hinges on its ability to correctly identify vertebrae in cases of rare anatomical variations. The content and code concerning VerSe can be accessed at: https://github.com/anjany/verse.


title: Ab initio many-body GW correlations in the electronic structure of LaNiO$_2$

abstract: We present an ab initio $GW$ self-energy calculation of the electronic structure of LaNiO$_2$. With respect to density-functional theory we find that in $GW$ the La 4$f$ states undergo an important $+$2 eV upward shift from the Fermi level, while the O 2$p$ states are pulled down by $-$1.5 eV, thus reinforcing the charge-transfer character of this material. However, $GW$ many-body effects leave the $d$-like bands at the Fermi level almost unaffected, so that the Fermi-surface topology is preserved, unlike in cuprates.


title: Sheaf Representations and Duality in Logic

abstract: The fundamental duality theories relating algebra and geometry that were discovered in the mid-20th century can also be applied to logic via its algebraization under categorical logic. They thereby result in known and new completeness theorems. This idea can be taken even further via what is sometimes called ``categorification'' to establish a new connection between logic and geometry, a glimpse of which can also be had in topos theory.


title: Diffusion synthetic acceleration for heterogeneous domains, compatible with voids

abstract: A standard approach to solving the S$_N$ transport equations is to use source iteration with diffusion synthetic acceleration (DSA). Although this approach is widely used and effective on many problems, there remain some practical issues with DSA preconditioning, particularly on highly heterogeneous domains. For large-scale parallel simulation, it is critical that both (i) preconditioned source iteration converges rapidly, and (ii) the action of the DSA preconditioner can be applied using fast, scalable solvers, such as algebraic multigrid (AMG). For heterogeneous domains, these two interests can be at odds. In particular, there exist DSA diffusion discretizations that can be solved rapidly using AMG, but they do not always yield robust/fast convergence of the larger source iteration. Conversely, there exist robust DSA discretizations where source iteration converges rapidly on difficult heterogeneous problems, but fast parallel solvers like AMG tend to struggle applying the action of such operators. Moreover, very few current methods for the solution of deterministic transport are compatible with voids. This paper develops a new heterogeneous DSA preconditioner based on only preconditioning the optically thick subdomains. The resulting method proves robust on a variety of heterogeneous transport problems, including a linearized hohlraum mesh related to inertial confinement fusion. Moreover, the action of the preconditioner is easily computed using $\mathcal{O}(1)$ AMG iterations, {convergence of the transport iteration typically requires $2-5\times$ less iterations than current state-of-the-art ``full DSA,'' and the proposed method is} trivially compatible with voids. On the hohlraum problem, rapid convergence is obtained by preconditioning less than 3\% of the mesh elements with $5-10$ AMG iterations.


title: Traveling Bands, Clouds, and Vortices of Chiral Active Matter

abstract: We consider stochastic dynamics of self-propelled particles with nonlocal normalized alignment interactions subject to phase lag. The role of the lag is to indirectly generate chirality into particle motion. To understand large scale behavior, we derive a continuum description of an active Brownian particle (ABP) flow with macroscopic scaling in the form of a partial differential equation (PDE) for a one-particle probability density function (DF). Due to indirect chirality, we find a new spatially homogeneous nonstationary analytic solution for this class of equations. Our development of kinetic and hydrodynamic theories towards such a solution reveals the existence of a wide variety of spatially nonhomogeneous patterns reminiscent of the traveling bands, clouds, and vortical structures of linear active matter. Our model may thereby serve as the basis for understanding the nature of chiral active media and designing multiagent swarms with designated behavior.


title: On simulation in automata networks

abstract: An automata network is a finite graph where each node holds a state from some finite alphabet and is equipped with an update function that changes its state according to the configuration of neighboring states. More concisely, it is given by a finite map $f:Q^n\rightarrow Q^n$. In this paper we study how some (sets of) automata networks can be simulated by some other (set of) automata networks with prescribed update mode or interaction graph. Our contributions are the following. For non-Boolean alphabets and for any network size, there are intrinsically non-sequential transformations (i.e. that can not be obtained as composition of sequential updates of some network). Moreover there is no universal automaton network that can produce all non-bijective functions via compositions of asynchronous updates. On the other hand, we show that there are universal automata networks for sequential updates if one is allowed to use a larger alphabet and then use either projection onto or restriction to the original alphabet. We also characterize the set of functions that are generated by non-bijective sequential updates. Following Tchuente, we characterize the interaction graphs $D$ whose semigroup of transformations is the full semigroup of transformations on $Q^n$, and we show that they are the same if we force either sequential updates only, or all asynchronous updates.


title: A 3D Reactive Navigation Algorithm for Mobile Robots by Using Tentacle-Based Sampling

abstract: This paper introduces a reactive navigation framework for mobile robots in 3-dimensional (3D) space. The proposed approach does not rely on the global map information and achieves fast navigation by employing a tentacle based sampling and their heuristic evaluations on-the-fly. This reactive nature of the approach comes from the prior arrangement of navigation points on tentacles (parametric contours) to sample the navigation space. These tentacles are evaluated at each time-step, based on heuristic features such as closeness to the goal, previous tentacle preferences and nearby obstacles in a robot-centered 3D grid. Then, the navigable sampling point on the selected tentacle is passed to a controller for the motion execution. The proposed framework does not only extend its 2D tentacle-based counterparts into 3D, but also introduces offline and online parameters, whose tuning provides versatility and adaptability of the algorithm to work in unknown environments. To demonstrate the superior performance of the proposed algorithm over a state-of-art method, the statistical results from physics-based simulations on various maps are presented. The video of the work is available at https://youtu.be/rrF7wHCz-0M.


title: Phase diagram of the large $N$ Gross-Neveu model in a finite periodic box

abstract: We analyze the Gross-Neveu model in the limit of large number of flavors of massless fermions. We study the phase diagram in a two and three dimensional periodic box at a fixed thermal to spatial aspect ratio, $\frac{\beta}{\ell}$, with a flavor independent chemical potential. We assume the bilinear condensate, when one exists, has a specific momentum in the spatial direction(s). The main known features of the phase diagram in the $\ell\to\infty$ limit of the two dimensional model are also seen on a finite $\ell\times\beta$ torus -- a phase with a homogeneous (zero momentum) condensate; a phase with an inhomogeneous (non-zero momentum) condensate and a phase with no condensate. We observe that the inhomogeneous phase contains several sub-phases characterized by a specific spatial momentum. Unlike the two dimensional model, we do not find evidence for a phase with a inhomogeneous condensate in the three dimensional model.


title: Comparison of Syntactic and Semantic Representations of Programs in Neural Embeddings

abstract: Neural approaches to program synthesis and understanding have proliferated widely in the last few years; at the same time graph based neural networks have become a promising new tool. This work aims to be the first empirical study comparing the effectiveness of natural language models and static analysis graph based models in representing programs in deep learning systems. It compares graph convolutional networks using different graph representations in the task of program embedding. It shows that the sparsity of control flow graphs and the implicit aggregation of graph convolutional networks cause these models to perform worse than naive models. Therefore it concludes that simply augmenting purely linguistic or statistical models with formal information does not perform well due to the nuanced nature of formal properties introducing more noise than structure for graph convolutional networks.


title: No loss of derivatives for hyperbolic operators with Zygmund-continuous coefficients in time

abstract: In this note we prove a well-posedness result, without loss of derivatives, for strictly hyperbolic wave operators having coefficients which are Zygmund-continuous in the time variable and Lipschitz-continuous in the space variables. The proof is based on Tarama's idea of introducing a lower order corrector in the energy, in order to produce special algebraic cancellations when computing its time derivative, combined with paradifferential calculus with parameters, in order to handle the low regularity of the coefficients with respect to $x$.


title: Modular network for high accuracy object detection

abstract: We present a novel modular object detection convolutional neural network that significantly improves the accuracy of object detection. The network consists of two stages in a hierarchical structure. The first stage is a network that detects general classes. The second stage consists of separate networks to refine the classification and localization of each of the general classes objects. Compared to a state of the art object detection networks the classification error in the modular network is improved by approximately 3-5 times, from 12% to 2.5 %-4.5%. This network is easy to implement and has a 0.94 mAP. The network architecture can be a platform to improve the accuracy of widespread state of the art object detection networks and other kinds of deep learning networks. We show that a deep learning network initialized by transfer learning becomes more accurate as the number of classes it later trained to detect becomes smaller.


title: Propagation of $L^p_{\beta}$-norm, $ 1<p\le \infty$, for the system of Boltzmann equations for monatomic gas mixtures

abstract: With the existence and uniqueness of a vector value solution for the full non-linear homogeneous Boltzmann system of equations describing multi-component monatomic gas mixtures for binary interactions proved [8], we present in this manuscript several properties for such a solution. We start by proving the gain of integrability of the gain term of the multispecies collision operator, extending the work done previously for the single species case [3]. In addition, we study the integrability properties of the multispecies collision operator as a bilinear form, revisiting and expanding the work done for a single gas [2]. With these estimates, together with a control by below for the loss term of the collision operator as in [8], we develop the propagation for the polynomially and exponentially $\beta$-weighted $L^p_\beta$-norms for the vector value solution. Finally, we extend such $L^p_\beta$-norms propagation property to $p=\infty$.


title: Cohomology of hyperfinite Borel actions

abstract: We study cocycles of countable groups $\Gamma$ of Borel automorphisms of a standard Borel space $(X, \mathcal{B})$ taking values in a locally compact second countable group $G$. We prove that for a hyperfinite group $\Gamma$ the subgroup of coboundaries is dense in the group of cocycles. We describe all Borel cocycles of the $2$-odometer and show that any such cocycle is cohomologous to a cocycle with values in a countable dense subgroup $H$ of $G$. We also provide a Borel version of Gottschalk-Hedlund theorem.


title: Gaussian-Smooth Optimal Transport: Metric Structure and Statistical Efficiency

abstract: Optimal transport (OT), and in particular the Wasserstein distance, has seen a surge of interest and applications in machine learning. However, empirical approximation under Wasserstein distances suffers from a severe curse of dimensionality, rendering them impractical in high dimensions. As a result, entropically regularized OT has become a popular workaround. However, while it enjoys fast algorithms and better statistical properties, it looses the metric structure that Wasserstein distances enjoy. This work proposes a novel Gaussian-smoothed OT (GOT) framework, that achieves the best of both worlds: preserving the 1-Wasserstein metric structure while alleviating the empirical approximation curse of dimensionality. Furthermore, as the Gaussian-smoothing parameter shrinks to zero, GOT $\Gamma$-converges towards classic OT (with convergence of optimizers), thus serving as a natural extension. An empirical study that supports the theoretical results is provided, promoting Gaussian-smoothed OT as a powerful alternative to entropic OT.


title: Interplay of curvature and rigidity in shape-based models of confluent tissue

abstract: Rigidity transitions in simple models of confluent cells have been a powerful organizing principle in understanding the dynamics and mechanics of dense biological tissue. In this work we explore the interplay between geometry and rigidity in two-dimensional vertex models confined to the surface of a sphere. By considering shapes of cells defined by perimeters whose magnitude depends on geodesic distances and areas determined by spherical polygons, the critical shape index in such models is affected by the size of the cell relative to the radius of the sphere on which it is embedded. This implies that cells can collectively rigidify by growing the size of the sphere, i.e. by tuning the curvature of their domain. Finite-temperature studies indicate that cell motility is affected well away from the zero-temperature transition point.


title: A random-walk model for dark matter halo spins

abstract: We extend the random-walk model of Vitvitska et al. for predicting the spins of dark matter halos from their merger histories. Using updated merger rates, orbital parameter distributions, and N-body constraints we show that this model can accurately reproduce the distribution of spin parameters measured in N-body simulations when we include a weak correlation between the spins of halos and the angular momenta of infalling subhalos. We further show that this model is in approximate agreement with the correlation of the spin magnitude over time as determined from N-body simulations, while it slightly underpredicts the correlation in the direction of the spin vector measured from the same simulations. This model is useful for predicting spins from merger histories derived from non-N-body sources, thereby circumventing the need for very high resolution simulations to permit accurate measurements of spins. It may be particularly relevant to modeling systems which accumulate angular momentum from halos over time (such as galactic discs) - we show that this model makes small but significant changes in the distribution of galactic disc sizes computed using the Galacticus semi-analytic galaxy formation model.


title: Detection of Thin Boundaries between Different Types of Anomalies in Outlier Detection using Enhanced Neural Networks

abstract: Outlier detection has received special attention in various fields, mainly for those dealing with machine learning and artificial intelligence. As strong outliers, anomalies are divided into the point, contextual and collective outliers. The most important challenges in outlier detection include the thin boundary between the remote points and natural area, the tendency of new data and noise to mimic the real data, unlabelled datasets and different definitions for outliers in different applications. Considering the stated challenges, we defined new types of anomalies called Collective Normal Anomaly and Collective Point Anomaly in order to improve a much better detection of the thin boundary between different types of anomalies. Basic domain-independent methods are introduced to detect these defined anomalies in both unsupervised and supervised datasets. The Multi-Layer Perceptron Neural Network is enhanced using the Genetic Algorithm to detect newly defined anomalies with higher precision so as to ensure a test error less than that calculated for the conventional Multi-Layer Perceptron Neural Network. Experimental results on benchmark datasets indicated reduced error of anomaly detection process in comparison to baselines.