title: HST/STIS transmission spectrum of the ultra-hot Jupiter WASP-76 b confirms the presence of sodium in its atmosphere

abstract: We present an atmospheric transmission spectrum of the ultra-hot Jupiter WASP-76 b by analyzing archival data obtained with the Space Telescope Imaging Spectrograph (STIS) on board the Hubble Space Telescope (HST). The dataset spans three transits, two with a wavelength coverage between 2900 and 5700 Armstrong, and the third one between 5250 and 10300 Armstrong. From the one-dimensional, time dependent spectra we constructed white and chromatic light curves, the latter with typical integration band widths of ~200 Armstrong. We computed the wavelength dependent planet-to-star radii ratios taking into consideration WASP-76's companion. The resulting transmission spectrum of WASP-76 b is dominated by a spectral slope of increasing opacity towards shorter wavelengths of amplitude of about three scale heights under the assumption of planetary equilibrium temperature. If the slope is caused by Rayleigh scattering, we derive a lower limit to the temperature of ~870 K. Following-up on previous detection of atomic sodium derived from high resolution spectra, we re-analyzed HST data using narrower bands centered around sodium. From an atmospheric retrieval of this transmission spectrum, we report evidence of sodium at 2.9-sigma significance. In this case, the retrieved temperature at the top of the atmosphere (10-5 bar) is 2300 +412-392 K. We also find marginal evidence for titanium hydride. However, additional high resolution ground-based data are required to confirm this discovery.


title: Counterterms, Kounterterms, and the variational problem in AdS gravity

abstract: We show that the Kounterterms for pure AdS gravity in arbitrary even dimensions coincide with the boundary counterterms obtained through holographic renormalization if and only if the boundary Weyl tensor vanishes. In particular, the Kounterterms lead to a well posed variational problem for generic asymptotically locally AdS manifolds only in four dimensions. We determine the exact form of the counterterms for conformally flat boundaries and demonstrate that, in even dimensions, the Kounterterms take exactly the same form. This agreement can be understood as a consequence of Anderson's theorem for the renormalized volume of conformally compact Einstein 4-manifolds and its higher dimensional generalizations by Albin and Chang, Qing and Yang. For odd dimensional asymptotically locally AdS manifolds with a conformally flat boundary, the Kounterterms coincide with the boundary counterterms except for the logarithmic divergence associated with the holographic conformal anomaly, and finite local terms.


title: Solvable Criterion for the Contextuality of any Prepare-and-Measure Scenario

abstract: Starting from arbitrary sets of quantum states and measurements, referred to as the prepare-and-measure scenario, an operationally noncontextual ontological model of the quantum statistics associated with the prepare-and-measure scenario is constructed. The operationally noncontextual ontological model coincides with standard Spekkens noncontextual ontological models for tomographically complete scenarios, while covering the non-tomographically complete case with a new notion of a reduced space, which we motivate following the guiding principles of noncontextuality. A mathematical criterion, called unit separability, is formulated as the relevant classicality criterion -- the name is inspired by the usual notion of quantum state separability. Using this criterion, we derive a new upper bound on the cardinality of the ontic space. Then, we recast the unit separability criterion as a (possibly infinite) set of linear constraints, from which we obtain two separate hierarchies of algorithmic tests to witness the non-classicality or certify the classicality of a scenario. Finally, we reformulate our results in the framework of generalized probabilistic theories and discuss the implications for simplex-embeddability in such theories.


title: Measuring the Tidal Response of Structure Formation: Anisotropic Separate Universe Simulations using TreePM

abstract: We present anisotropic "separate universe" simulations which modify the N-body code Gadget-4 in order to represent a large-scale tidal field through an anisotropic expansion factor. These simulations are used to measure the linear, quasi-linear and nonlinear response of the matter power spectrum to a spatially uniform trace-free tidal field up to wavenumber $k = 7 h \text{Mpc}^{-1}$. Together with the response to a large-scale overdensity measured in previous work, this completely describes the nonlinear matter bispectrum in the squeezed limit. We find that the response amplitude does not approach zero on small scales in physical coordinates, but rather a constant value at $z=0$, $R_K\approx 0.5$ for $k \geq 3 h\text{Mpc}^{-1}$ up to the scale where we consider our simulations reliable, $k \leq 7 h\text{Mpc}^{-1}$ at $z=0$. This shows that even the inner regions of haloes are affected by the large-scale tidal field. We also measure directly the alignment of halo shapes with the tidal field, finding a clear signal which increases with halo mass.


title: Towards a Resilient Machine Learning Classifier -- a Case Study of Ransomware Detection

abstract: The damage caused by crypto-ransomware, due to encryption, is difficult to revert and cause data losses. In this paper, a machine learning (ML) classifier was built to early detect ransomware (called crypto-ransomware) that uses cryptography by program behavior. If a signature-based detection was missed, a behavior-based detector can be the last line of defense to detect and contain the damages. We find that input/output activities of ransomware and the file-content entropy are unique traits to detect crypto-ransomware. A deep-learning (DL) classifier can detect ransomware with a high accuracy and a low false positive rate. We conduct an adversarial research against the models generated. We use simulated ransomware programs to launch a gray-box analysis to probe the weakness of ML classifiers and to improve model robustness. In addition to accuracy and resiliency, trustworthiness is the other key criteria for a quality detector. Making sure that the correct information was used for inference is important for a security application. The Integrated Gradient method was used to explain the deep learning model and also to reveal why false negatives evade the detection. The approaches to build and to evaluate a real-world detector were demonstrated and discussed.


title: Absorption line abundances in the SMC-like galaxy UGC 5282: evidence for ISM dilution from inflows on kiloparsec scales

abstract: We present an HST Cosmic Origins Spectrograph (COS) spectrum of the QSO SDSS J095109.12+330745.8 ($z_{em} = 0.645$) whose sightline passes through the SMC-like dwarf galaxy UGC 5282 ($M_B = -16.0$, $cz = 1577$ km s$^{-1}$ ), 1.2 kpc in projection from the central HII region of the galaxy. Damped Ly$\alpha$ (DLA) absorption is detected at the redshift of UGC 5282 with $\log$[$N$(HI) cm$^{-2}\:$] = 20.89$^{+0.12}_{-0.21}\:\:$. Analysis of the accompanying SII, PII and OI metal lines yields a neutral gas metallicity, $Z_{HI}$, of [S/H] $\simeq$ [P/H] $= -0.80\pm0.24$. The metallicity of ionized gas from the central HII region, $Z_{HII}\:$, measured from its emission lines is [O/H]$= -0.37\pm0.10$, a difference of $+0.43\pm 0.26$ from $Z_{HI}\:$. This difference $\delta$ is consistent with that seen towards HII regions in other star-forming galaxies and supports the idea that ionized gas near star forming regions shows systematically higher metallicities than exist in the rest of a galaxy's neutral interstellar medium (ISM). The positive values of $\delta$ found in UGC 5282 (and the other star forming galaxies) is likely due to infalling low metallicity gas from the intergalactic medium that mixes with the galaxy's ISM on kpc scales. This model is also consistent with broad Ly$\alpha$ emission detected at the bottom of the DLA absorption, offset by $\sim 125$ km s$^{-1}$ from the absorption velocity. Models of galaxy evolution that attempt to replicate population characteristics, such as the mass-metallicity relation, may need to start with a galaxy metallicity represented by $Z_{HI}$ rather than that measured traditionally from $Z_{HII}\:$.


title: Learning Unbiased Representations via Mutual Information Backpropagation

abstract: We are interested in learning data-driven representations that can generalize well, even when trained on inherently biased data. In particular, we face the case where some attributes (bias) of the data, if learned by the model, can severely compromise its generalization properties. We tackle this problem through the lens of information theory, leveraging recent findings for a differentiable estimation of mutual information. We propose a novel end-to-end optimization strategy, which simultaneously estimates and minimizes the mutual information between the learned representation and the data attributes. When applied on standard benchmarks, our model shows comparable or superior classification performance with respect to state-of-the-art approaches. Moreover, our method is general enough to be applicable to the problem of ``algorithmic fairness'', with competitive results.


title: Random Forest Classifier Based Prediction of Rogue waves on Deep Oceans

abstract: In this paper, we present a novel approach for the prediction of rogue waves in oceans using statistical machine learning methods. Since the ocean is composed of many wave systems, the change from a bimodal or multimodal directional distribution to unimodal one is taken as the warning criteria. Likewise, we explore various features that help in predicting rogue waves. The analysis of the results shows that the Spectral features are significant in predicting rogue waves. We find that nonlinear classifiers have better prediction accuracy than the linear ones. Finally, we propose a Random Forest Classifier based algorithm to predict rogue waves in oceanic conditions. The proposed algorithm has an Overall Accuracy of 89.57% to 91.81%, and the Balanced Accuracy varies between 79.41% to 89.03% depending on the forecast time window. Moreover, due to the model-free nature of the evaluation criteria and interdisciplinary characteristics of the approach, similar studies may be motivated in other nonlinear dispersive media, such as nonlinear optics, plasma, and solids, governed by similar equations, which will allow for the early detection of extreme waves


title: On the short term modulation of cosmic rays by high-speed streams at the Pierre Auger surface array detectors

abstract: We present an analysis of the short-term modulation (one rotation of Bartels-27 days) of the galactic cosmic rays (GCR) by the solar wind, based on the cosmic ray rates observed by the Pierre Auger Observatory (PAO) on their surface detectors in scaler mode. The incidence of GCR with energies below $\sim$ 50 TeV, at the top of the atmosphere, produces more than 90\% of the secondary particles registered at ground level, i.e., they are subject to solar modulation. The modulation is consistent with at least two components: The first is the modulation of the amplitude of the cosmic rays diurnal variation, anti-correlated with the solar-wind speed. The second one occurs during the high-speed stream (HSS), the baseline of the cosmic rays diurnal variation train falls, following the time profile of the solar-wind speed inversely. Based on the radial gradient of the cosmic ray diffusion theory and under some other premises, such as the latitude dependence on diurnal variation and the inclusion of drift processes in the propagation of GCR, a semi-empirical description of the modulation is possible to do, and it hereafter is called as Toy-model. Although the Toy-model does not include fluctuations due to propagation in the atmosphere, it provides satisfactory results when compared with the PAO scaler mode data. We present details of these observations as well as the Toy-model validation.


title: Tomaszewski's problem on randomly signed sums, revisited

abstract: Let $v_1$, $v_2$, ..., $v_n$ be real numbers whose squares add up to 1. Consider the $2^n$ signed sums of the form $S = \sum \pm v_i$. Boppana and Holzman (2017) proved that at least 13/32 of these sums satisfy $|S| \le 1$. Here we improve their bound to $0.427685$.


title: A Neural Architecture for Detecting Confusion in Eye-tracking Data

abstract: Encouraged by the success of deep learning in a variety of domains, we investigate a novel application of its methods on the effectiveness of detecting user confusion in eye-tracking data. We introduce an architecture that uses RNN and CNN sub-models in parallel to take advantage of the temporal and visuospatial aspects of our data. Experiments with a dataset of user interactions with the ValueChart visualization tool show that our model outperforms an existing model based on Random Forests resulting in a 22% improvement in combined sensitivity & specificity.


title: Spectrally Efficient Pilot Structure and Channel Estimation for Multiuser FBMC Systems

abstract: In this paper, we consider channel estimation problem in the uplink of filter bank multicarrier (FBMC) systems. We propose a pilot structure and a joint multiuser channel estimation method for FBMC. Opposed to the available solutions in the literature, our proposed technique does not rely on the flat-channel condition over each subcarrier band or any requirement for placing guard symbols between different users' pilots. Our proposed pilot structure reduces the training overhead by interleaving the users' pilots in time and frequency. Thus, we can accommodate a larger number of training signals within the same bandwidth and improve the spectral efficiency. Furthermore, this pilot structure inherently leads to a reduced peak-to-average power ratio (PAPR) compared with the solutions that use all the subcarriers for training. We analytically derive the Cramer-Rao lower bound (CRLB) and mean square error (MSE) expressions for our proposed method. We show that these expressions are the same. This confirms the optimality of our proposed method, which is numerically evaluated through simulations. Relying on its improved spectral efficiency, our proposed method can serve a large number of users and relax pilot contamination problem in FBMC-based massive MIMO systems. This is corroborated through simulations in terms of sum-rate performance for both single cell and multicell scenarios.


title: Early neutron star evolution in high-mass X-ray binaries

abstract: The application of standard accretion theory to observations of X-ray binaries provides valuable insights into neutron star properties, such as their spin period and magnetic field. However, most studies concentrate on relatively old systems, where the neutron star is in its late propeller, accretor, or nearly spin equilibrium phase. Here we use an analytic model from standard accretion theory to illustrate the evolution of high-mass X-ray binaries early in their life. We show that a young neutron star is unlikely to be an accretor because of the long duration of ejector and propeller phases. We apply the model to the recently discovered ~4000 yr old high-mass X-ray binary XMMU J051342.6-672412 and find that the system's neutron star, with a tentative spin period of 4.4 s, cannot be in the accretor phase and has a magnetic field B > (a few)x10^13 G, which is comparable to the magnetic field of many older high-mass X-ray binaries and is much higher than the spin equilibrium inferred value of (a few)x10^11 G. The observed X-ray luminosity could be the result of thermal emission from a young cooling magnetic neutron star or a small amount of accretion that can occur in the propeller phase.


title: Work as an external quantum observable and an operational quantum work fluctuation theorem

abstract: We propose a definition of externally measurable quantum work in driven systems. Work is given as a quantum observable on a control device which is forcing the system and can be determined without knowledge of the system Hamiltonian $H_\mathcal{S}$. We argue that quantum work fluctuation theorems which rely on the knowledge of $H_\mathcal{S}$ are of little practical relevance, contrary to their classical counterparts. Using our framework, we derive a fluctuation theorem which is operationally accessible and could in principle be implemented in experiments to determine bounds on free energy differences of unknown systems.


title: Dissipative Kerr Solitons in a III-V Microresonator

abstract: We demonstrate stable microresonator Kerr soliton frequency combs in a III-V platform (AlGaAs on SiO$_2$) through quenching of thermorefractive effects by cryogenic cooling to temperatures between 4~K and 20~K. This cooling reduces the resonator's thermorefractive coefficient, whose room-temperature value is an order of magnitude larger than that of other microcomb platforms like Si$_3$N$_4$, SiO$_2$, and AlN, by more than two orders of magnitude, and makes soliton states adiabatically accessible. Realizing such phase-stable soliton operation is critical for applications that fully exploit the ultra-high effective nonlinearity and high optical quality factors exhibited by this platform.


title: Mutual Information Maximization for Effective Lip Reading

abstract: Lip reading has received an increasing research interest in recent years due to the rapid development of deep learning and its widespread potential applications. One key point to obtain good performance for the lip reading task depends heavily on how effective the representation can be to capture the lip movement information and meanwhile to resist the noises resulted from the change of pose, lighting conditions, speaker's appearance and so on. Towards this target, we propose to introduce the mutual information constraints on both the local feature's level and the global sequence's level to enhance the relations of the features with the speech content. On the one hand, we constraint the features generated at each time step to enable them carry a strong relation with the speech content by imposing the local mutual information maximization constraint (LMIM), leading to improvements over the model's ability to discover fine-grained lip movements and the fine-grained differences among words with similar pronunciation, such as ``spend'' and ``spending''. On the other hand, we introduce the mutual information maximization constraint on the global sequence's level (GMIM), to make the model be able to pay more attention to discriminate key frames related with the speech content, and less to various noises appeared in the speaking process. By combining these two advantages together, the proposed method is expected to be both discriminative and robust for effective lip reading. To verify this method, we evaluate on two large-scale benchmark. We perform a detailed analysis and comparison on several aspects, including the comparison of the LMIM and GMIM with the baseline, the visualization of the learned representation and so on. The results not only prove the effectiveness of the proposed method but also report new state-of-the-art performance on both the two benchmarks.


title: Dead or Alive? Implications of the Muon Anomalous Magnetic Moment for 3-3-1 Models

abstract: We have witnessed a persistent puzzling anomaly in the muon magnetic moment that cannot be accounted for in the Standard Model even considering the large hadronic uncertainties. A new measurement is forthcoming, and it might give rise to a $5\sigma$ claim for physics beyond the Standard Model. Motivated by it, we explore the implications of this new result to five models based on the $SU(3)_C \times SU(3)_L \times U(1)_N$ gauge symmetry and put our conclusions into perspective with LHC bounds. We show that previous conclusions found in the context of such models change if there are more than one heavy particle running in the loop. Moreover, having in mind the projected precision aimed by the g-2 experiment at FERMILAB, we place lower mass bounds on the particles that contribute to muon anomalous magnetic moment assuming the anomaly is resolved otherwise. Lastly, we discuss how these models could accommodate such anomaly in agreement with existing bounds.


title: Neural Generators of Sparse Local Linear Models for Achieving both Accuracy and Interpretability

abstract: For reliability, it is important that the predictions made by machine learning methods are interpretable by human. In general, deep neural networks (DNNs) can provide accurate predictions, although it is difficult to interpret why such predictions are obtained by DNNs. On the other hand, interpretation of linear models is easy, although their predictive performance would be low since real-world data is often intrinsically non-linear. To combine both the benefits of the high predictive performance of DNNs and high interpretability of linear models into a single model, we propose neural generators of sparse local linear models (NGSLLs). The sparse local linear models have high flexibility as they can approximate non-linear functions. The NGSLL generates sparse linear weights for each sample using DNNs that take original representations of each sample (e.g., word sequence) and their simplified representations (e.g., bag-of-words) as input. By extracting features from the original representations, the weights can contain rich information to achieve high predictive performance. Additionally, the prediction is interpretable because it is obtained by the inner product between the simplified representations and the sparse weights, where only a small number of weights are selected by our gate module in the NGSLL. In experiments with real-world datasets, we demonstrate the effectiveness of the NGSLL quantitatively and qualitatively by evaluating prediction performance and visualizing generated weights on image and text classification tasks.


title: Generating sets and representability for symplectic capacities

abstract: K. Cieliebak, H. Hofer, J. Latschev, and F. Schlenk (CHLS) posed the problem of finding a minimal generating set for the (symplectic) capacities on a given symplectic category. We show that if the category contains a certain one-parameter family of objects, then every countably Borel-generating set of (normalized) capacities has cardinality (strictly) bigger than the continuum. This appears to be the first result regarding the problem of CHLS, except for two results of D. McDuff about the category of ellipsoids in dimension 4. We also prove that every finitely differentiably generating set of capacities on a given symplectic category is uncountable, provided that the category contains a one-parameter family of symplectic manifolds that is ``strictly volume-increasing'' and ``embedding-capacity-wise constant''. It follows that the Ekeland-Hofer capacities and the volume capacity do not finitely differentiably generate all generalized capacities on the category of ellipsoids. This answers a variant of a question of CHLS. In addition, we prove that if a given symplectic category contains a certain one-parameter family of objects, then almost no normalized capacity is domain- or target-representable. This provides some solutions to two central problems of CHLS.


title: Star formation law in the EoR from [CII] and CIII] lines

abstract: We present a novel method to simultaneously characterise the star formation law and the interstellar medium properties of galaxies in the Epoch of Reionization (EoR) through the combination of [CII] 158$\mu$m (and its known relation with star formation rate) and CIII]$\lambda$1909{\AA} emission line data. The method, based on a Markov Chain Monte Carlo algorithm, allows to determine the target galaxy average density, $n$, gas metallicity, $Z$, and burstiness parameter, $\kappa_s$, quantifying deviations from the Kennicutt-Schmidt relation. As an application, we consider COS-3018 (z=6.854), the only EoR Lyman Break Galaxy so far detected in both [CII] and CIII]. We show that COS-3018 is a moderate starburst ($\kappa_s \approx 3$), with $Z\approx 0.4\, Z_{\odot}$, and $n \approx 500\, {\rm cm^{-3}}$. Our method will be optimally applied to joint ALMA and JWST targets.


title: Neutrino phenomenology in a left-right $D_4$ symmetric model

abstract: We present a minimal left-right symmetric flavor model and analyze the predictions for the neutrino sector. In this scenario, the Yukawa sector is shaped by the dihedral $D_4$ symmetry which leads to correlations for the neutrino mixing parameters. We end up with four possible solutions within this model. We further analyzed the impact of the upcoming long-baseline neutrino oscillation experiment DUNE. Due to its high sensitivity, DUNE will be able to rule out two of the solutions. Finally, the prediction for the neutrinoless double beta decay for the model has also been examined.


title: Stationary black holes and light rings

abstract: The ringdown and shadow of the astrophysically significant Kerr Black Hole (BH) are both intimately connected to a special set of bound null orbits known as Light Rings (LRs). Does it hold that a generic equilibrium BH must possess such orbits? In this letter we prove the following theorem. A stationary, axi-symmetric, asymptotically flat black hole spacetime in 1+3 dimensions, with a non-extremal, topologically spherical, Killing horizon admits, at least, one standard LR outside the horizon for each rotation sense. The proof relies on a topological argument and assumes $C^2$-smoothness and circularity, but makes no use of the field equations. The argument is also adapted to recover a previous theorem establishing that a horizonless ultra-compact object must admit an even number of non-degenerate LRs, one of which is stable.


title: Modeling ultrasound propagation in the moving brain: applications to shear shock waves and traumatic brain injury

abstract: Traumatic brain injury studies on the living human brain are experimentally infeasible. We present a simulation approach that models ultrasound propagation in the human brain while it is moving due to the complex shear shock wave deformation from a traumatic impact. Finite difference simulations can model ultrasound propagation in complex media such as human tissue. Recently, we have shown that a finite difference approach can be used to represent displacements that are much smaller than the grid size, such as the motion encountered in shear wave propagation from ultrasound elastography. However, this subresolution displacement model, called impedance flow, was only implemented and validated for acoustical media composed of randomly distributed scatterers. Here we propose a generalization of the impedance flow method that describes the continuous subresolution motion of structured acoustical maps, and in particular of acoustical maps of the human brain. It is shown that the average error in the subresolution displacement method is small compared to the wavelength. The method is applied to acoustical maps of a moving human brain, imposed by the propagation of a shear shock wave. Then the Fullwave simulation tool is used to model transmit-receive imaging sequences based on an L7-4 imaging transducer. The simulated radiofrequency data is beamformed using a conventional delay-and-sum method and a normalized cross-correlation method designed for shock wave tracking is used to determine the tissue motion. It is shown that the proposed generalized impedance flow method accurately captures the shear shock wave motion. This method can lead to improvements in image sequence design that takes into account the aberration and multiple reflections from the brain and in the design of tracking algorithms that can more accurately capture the complex brain motion that occurs during a traumatic impact.


title: Towards windproofing LIGO: Reducing the effect of wind-driven floor tilt by using rotation sensors in active seismic isolation

abstract: Modern gravitational-wave observatories require robust low-frequency active seismic isolation in order to keep the interferometer at its ideal operating conditions. Seismometers are used to measure both the motion of the ground and isolated platform. These devices are susceptible to contamination from ground tilt at frequencies below 0.1 Hz, particularly arising from wind-pressure acting on building walls. Consequently, during LIGO's first observing run both observatories suffered significant downtime when wind-speeds were above 7 m/s. We describe the use of ground rotation sensors at the LIGO Hanford Observatory to correct nearby ground seismometers to produce tilt-free ground translation signals. The use of these signals for sensor correction control improved low-frequency seismic isolation and allowed the observatory to operate under wind speeds as high as $15-20$ m/s.


title: On the Generalised Langevin Equation for Simulated Annealing

abstract: In this paper, we consider the generalised (higher order) Langevin equation for the purpose of simulated annealing and optimisation of nonconvex functions. Our approach modifies the underdamped Langevin equation by replacing the Brownian noise with an appropriate Ornstein-Uhlenbeck process to account for memory in the system. Under reasonable conditions on the loss function and the annealing schedule, we establish convergence of the continuous time dynamics to a global minimum. In addition, we investigate the performance numerically and show better performance and higher exploration of the state space compared to the underdamped and overdamped Langevin dynamics with the same annealing schedule.


title: Extreme Ultraviolet and Soft X-Ray Diffraction Efficiency of a Blazed Reflection Grating Fabricated by Thermally Activated Selective Topography Equilibration

abstract: Future observatories utilizing reflection grating spectrometers for extreme ultraviolet (EUV) and soft X-ray (SXR) spectroscopy require high-fidelity gratings with both blazed groove facets and custom groove layouts that are often fanned or feature a slight curvature. While fabrication procedures centering on wet anisotropic etching in monocrystalline silicon produce highly efficient blazed gratings, the precision of a nonparallel groove layout is limited by the cubic structure of the silicon crystal. This motivates the pursuit of alternative techniques to grating manufacture, namely thermally activated selective topography equilibration (TASTE), which uses gray-scale electron-beam lithography to pattern multilevel structures in resist followed by an optimized polymer thermal reflow to smooth the 3D patterns into continuous surface relief profiles. Using TASTE, a mold for a reflection grating with a periodicity of 400 nm and grooves resembling an asymmetric sawtooth was patterned in 130 nm thick poly(methyl methacrylate) resist on a silicon substrate over a 50 mm by 7.5 mm area. This structure was coated with 15 nm of gold by electron-beam physical vapor deposition using titanium as an adhesion layer and then tested for EUV and SXR diffraction efficiency at beamline 6.3.2 of the Advanced Light Source synchrotron facility. Results demonstrate a quasi-blaze response characteristic of a 27 degree blaze angle with groove facets smooth to 1.5 nm rms. Absolute peak-order efficiency ranges from 75% to 25%, while total relative efficiency measures gap90% across the measured bandpass of 15.5 nm > lambda > 1.55 nm.


title: On bucket increasing trees, clustered increasing trees and increasing diamonds

abstract: In this work we analyze bucket increasing tree families. We introduce two simple stochastic growth processes, generating random bucket increasing trees of size $n$, complementing the earlier result of Mahmoud and Smythe for bucket recursive trees. On the combinatorial side, we define multilabelled generalizations of the tree families $d$-ary increasing trees and generalized plane-oriented recursive trees. Additionally, we introduce a clustering process for ordinary increasing trees and relate it to bucket increasing trees. We discuss in detail the bucket size two and present a bijection between such bucket increasing tree families and certain families of graphs called increasing diamonds, providing an explanation for phenomena observed by Bodini et al.


title: The GraphNet Zoo: An All-in-One Graph Based Deep Semi-Supervised Framework for Medical Image Classification

abstract: We consider the problem of classifying a medical image dataset when we have a limited amounts of labels. This is very common yet challenging setting as labelled data is expensive, time consuming to collect and may require expert knowledge. The current classification go-to of deep supervised learning is unable to cope with such a problem setup. However, using semi-supervised learning, one can produce accurate classifications using a significantly reduced amount of labelled data. Therefore, semi-supervised learning is perfectly suited for medical image classification. However, there has almost been no uptake of semi-supervised methods in the medical domain. In this work, we propose an all-in-one framework for deep semi-supervised classification focusing on graph based approaches, which up to our knowledge it is the first time that an approach with minimal labels has been shown to such an unprecedented scale with medical data. We introduce the concept of hybrid models by defining a classifier as a combination between an energy-based model and a deep net. Our energy functional is built on the Dirichlet energy based on the graph p-Laplacian. Our framework includes energies based on the $\ell_1$ and $\ell_2$ norms. We then connected this energy model to a deep net to generate a much richer feature space to construct a stronger graph. Our framework can be set to be adapted to any complex dataset. We demonstrate, through extensive numerical comparisons, that our approach readily compete with fully-supervised state-of-the-art techniques for the applications of Malaria Cells, Mammograms and Chest X-ray classification whilst using only 20% of labels.


title: How Fast Can We Insert? An Empirical Performance Evaluation of Apache Kafka

abstract: Message brokers see widespread adoption in modern IT landscapes, with Apache Kafka being one of the most employed platforms. These systems feature well-defined APIs for use and configuration and present flexible solutions for various data storage scenarios. Their ability to scale horizontally enables users to adapt to growing data volumes and changing environments. However, one of the main challenges concerning message brokers is the danger of them becoming a bottleneck within an IT architecture. To prevent this, knowledge about the amount of data a message broker using a specific configuration can handle needs to be available. In this paper, we propose a monitoring architecture for message brokers and similar Java Virtual Machine-based systems. We present a comprehensive performance analysis of the popular Apache Kafka platform using our approach. As part of the benchmark, we study selected data ingestion scenarios with respect to their maximum data ingestion rates. The results show that we can achieve an ingestion rate of about 420,000 messages/second on the used commodity hardware and with the developed data sender tool.


title: Extended Red Emission: Observational constraints for models

abstract: Extended Red Emission (ERE) is a widely observed optical emission process, present in a wide range of circumstellar and interstellar environments in the Milky Way galaxy as well as other galaxies. Definitive identifications of the ERE carriers and the ERE process are still a matter of debate. Numerous models have been proposed in recent decades, often developed without consideration of the growing body of observational constraints, which by now invalidate many of these models. This review focuses on the most well-established observational constraints which should help to delineate the way toward a generally accepted explanation of the ERE and an understanding of its place in the radiation physics of the interstellar medium.


title: A Note on Stein's Method for Heavy-Traffic Analysis

abstract: In this note, we apply Stein's method to analyze the steady-state distribution of queueing systems in the traditional heavy-traffic regime. Compared to previous methods (e.g., drift method and transform method), Stein's method allows us to establish stronger results with simple and template proofs. In particular, we consider discrete-time systems in this note. We first introduce the key ideas of Stein's method for heavy-traffic analysis through a single-server system. Then, we apply the developed template to analyze both load balancing problems and scheduling problems. All these three examples demonstrate the power and flexibility of Stein's method in heavy-traffic analysis. In particular, we can see that one appealing property of Stein's method is that it combines the advantages of both the drift method and the transform method.


title: Two-Sample High Dimensional Mean Test Based On Prepivots

abstract: Testing equality of mean vectors is a very commonly used criterion when comparing two multivariate random variables. Traditional tests such as Hotelling's T-squared become either unusable or output small power when the number of variables is greater than the combined sample size. In this paper, w}e propose a test using both prepivoting and Edgeworth expansion for testing the equality of two population mean vectors in the "large p, small n" setting. The asymptotic null distribution of the test statistic is derived and it is shown that the power of suggested test converges to one under certain alternatives when both n and p increase to infinity against sparse alternatives. Finite sample performance of the proposed test statistic is compared with other recently developed tests designed to also handle the "large p, small n" situation through simulations. The proposed test achieves competitive rates for both type I error rate and power. The usefulness of our test is illustrated by applications to two microarray gene expression data sets


title: On compact subsets of Sobolev spaces on manifolds

abstract: It is common that a Sobolev space defined on $\mathbb{R}^m$ has a non-compact embedding into an $L^p$-space, but it has subspaces for which this embedding becomes compact. There are three well known cases of such subspaces, the Rellich compactness, for a subspace of functions on a bounded domain (or an unbounded domain, sufficiently thin at infinity), the Strauss compactness, for a subspace of radially symmetric functions in $\mathbb{R}^m$, and the weighted Sobolev spaces. Known generalizations of Strauss compactness include subspaces of functions with block-radial symmetry, subspaces of functions with certain symmetries on Riemannian manifolds, as well as similar subspaces of more general Besov and Triebel-Lizorkin spaces. Presence of symmetries can be interpreted in terms of the rising critical Sobolev exponent corresponding to the smaller effective dimension of the quotient space.


title: Interactions between Hlawka Type-1 and Type-2 Quantities

abstract: The classical Hlawka inequality possesses deep connections with zonotopes and zonoids in convex geometry, and has been related to Minkowski space. We introduce Hlawka Type-1 and Type-2 quantities, and establish a Hlawka-type relation between them, which connects a vast number of strikingly different variants of the Hlawka inequalities, such as Serre's reverse Hlawka inequality in the future cone of the Minkowski space, the Hlawka inequality for subadditive function on abelian group by Ressel, and the integral analogs by Takahasi et al. Besides, we announce several enhanced results, such as the Hlawka inequality for the power of measure function. Particularly, we give a complete study of the Hlawka inequality for quadratic form which relates to a work of Serre.


title: QED at Large: A Survey of Engineering of Formally Verified Software

abstract: Development of formal proofs of correctness of programs can increase actual and perceived reliability and facilitate better understanding of program specifications and their underlying assumptions. Tools supporting such development have been available for over 40 years, but have only recently seen wide practical use. Projects based on construction of machine-checked formal proofs are now reaching an unprecedented scale, comparable to large software projects, which leads to new challenges in proof development and maintenance. Despite its increasing importance, the field of proof engineering is seldom considered in its own right; related theories, techniques, and tools span many fields and venues. This survey of the literature presents a holistic understanding of proof engineering for program correctness, covering impact in practice, foundations, proof automation, proof organization, and practical proof development.


title: Semi-relativistic (e,2e) study with twisted electron beam on Cu and Ag

abstract: In this communication, we report our calculations of Triple Differential Cross-Section (TDCS) for the relativistic (e,2e) process with a twisted electron beam on Cu and Ag atomic targets in coplanar asymmetric geometry mode. The theoretical formalism has been developed in the first Born approximation (FBA) in which we use the Dirac plane wave as well as the twisted electron wave for the incident electron beam to study the effect of various parameters of the twisted electron beam on the (e,2e) process. We use Dirac plane wave, semi relativistic Coulomb wave and Darwin wave function for the scattered, ejected and K-shell electron respectively. We compare the angular profiles of the TDCS of the twisted electron impact (e,2e) process with that of the plane wave. We segregate the TDCS for charge-charge interaction and current-current interaction with their interference term and study the effect of different parameters of the twisted electron beam on them. The study is also extended to the macroscopic Cu and Ag targets to further investigate the effect of the opening angle of the twisted electron beam on TDCS. The spin asymmetry in TDCS caused by polarized incident electron beam is also studied to elucidate the effects of the twisted electron beam on the (e,2e) process.


title: The CMS Trigger Upgrade for the HL-LHC

abstract: The CMS experiment has been designed with a two-level trigger system: the Level-1 Trigger, implemented on custom-designed electronics, and the High Level Trigger, a streamlined version of the CMS offline reconstruction software running on a computer farm. During its second phase the LHC will reach a luminosity of $7.5\times10^{34}\,\textrm{cm}^{-2}\,\textrm{s}^{-1}$ with a pileup of 200 collisions, producing integrated luminosity greater than 3000 fb$^{-1}$ over the full experimental run. To fully exploit the higher luminosity, the CMS experiment will introduce a more advanced Level-1 Trigger and increase the full readout rate from 100 kHz to 750 kHz. CMS is designing an efficient data-processing hardware trigger (Level-1) that will include tracking information and high-granularity calorimeter information. The current conceptual system design is expected to take full advantage of advances in FPGA and link technologies over the coming years, providing a high-performance, low-latency system for large throughput and sophisticated data correlation across diverse sources. The higher luminosity, event complexity and input rate present an unprecedented challenge to the High Level Trigger, that aims to achieve a similar efficiency and rejection factor as today despite the higher pileup and more pure preselection. In this presentation we will discuss the ongoing studies and prospects for the online reconstruction and selection algorithms for the high-luminosity era.


title: Exploring User Opinions of Fairness in Recommender Systems

abstract: Algorithmic fairness for artificial intelligence has become increasingly relevant as these systems become more pervasive in society. One realm of AI, recommender systems, presents unique challenges for fairness due to trade offs between optimizing accuracy for users and fairness to providers. But what is fair in the context of recommendation--particularly when there are multiple stakeholders? In an initial exploration of this problem, we ask users what their ideas of fair treatment in recommendation might be, and why. We analyze what might cause discrepancies or changes between user's opinions towards fairness to eventually help inform the design of fairer and more transparent recommendation algorithms.


title: A Persistent Homology Approach to Time Series Classification

abstract: Topological Data Analysis (TDA) is a rising field of computational topology in which the topological structure of a data set can be observed by persistent homology. By considering a sequence of sublevel sets, one obtains a filtration that tracks changes in topological information. These changes can be recorded in multi-sets known as {\it persistence diagrams}. Converting information stored in persistence diagrams into a form compatible with modern machine learning algorithms is a major vein of research in TDA. {\it Persistence curves}, a recently developed framework, provides a canonical and flexible way to encode the information presented in persistence diagrams into vectors. In this work, we propose a new set of metrics based on persistence curves. We prove the stability of the proposed metrics. Finally, we apply these metrics to the UCR Time Series Classification Archive. These empirical results show that our metrics perform better than the relevant benchmark in most cases and warrant further study.


title: Scale-specific dynamics of large-amplitude bursts in EEG capture behaviorally meaningful variability

abstract: Cascading large-amplitude bursts in neural activity, termed avalanches, are thought to provide insight into the complex spatially distributed interactions in neural systems. In human neuroimaging, for example, avalanches occurring during resting-state show scale-invariant dynamics, supporting the hypothesis that the brain operates near a critical point that enables long range spatial communication. In fact, it has been suggested that such scale-invariant dynamics, characterized by a power-law distribution in these avalanches, are universal in neural systems and emerge through a common mechanism. While the analysis of avalanches and subsequent criticality is increasingly seen as a framework for using complex systems theory to understand brain function, it is unclear how the framework would account for the omnipresent cognitive variability, whether across individuals and/or tasks. To address this, we analyzed avalanches in the EEG activity of healthy humans during rest as well as two distinct task conditions that varied in cognitive demands and produced behavioral measures unique to each individual. In both rest and task conditions we observed that avalanche dynamics demonstrate scale-invariant characteristics, but differ in their specific features, demonstrating individual variability. Using a new metric we call normalized engagement, which estimates the likelihood for a brain region to produce high-amplitude bursts, we also investigated regional features of avalanche dynamics. Normalized engagement showed not only the expected individual and task dependent variability, but also scale-specificity that correlated with individual behavior. Our findings expand our understanding of avalanche features and are supportive of the emerging theoretical idea that the dynamics of an active human brain operate close to a critical-like region and not a singular critical-state.


title: LCP: A Low-Communication Parallelization Method for Fast Neural Network Inference in Image Recognition

abstract: Deep neural networks (DNNs) have inspired new studies in myriad edge applications with robots, autonomous agents, and Internet-of-things (IoT) devices. However, performing inference of DNNs in the edge is still a severe challenge, mainly because of the contradiction between the intensive resource requirements of DNNs and the tight resource availability in several edge domains. Further, as communication is costly, taking advantage of other available edge devices by using data- or model-parallelism methods is not an effective solution. To benefit from available compute resources with low communication overhead, we propose the first DNN parallelization method for reducing the communication overhead in a distributed system. We propose a low-communication parallelization (LCP) method in which models consist of several almost-independent and narrow branches. LCP offers close-to-minimum communication overhead with better distribution and parallelization opportunities while significantly reducing memory footprint and computation compared to data- and model-parallelism methods. We deploy LCP models on three distributed systems: AWS instances, Raspberry Pis, and PYNQ boards. We also evaluate the performance of LCP models on a customized hardware (tailored for low latency) implemented on a small edge FPGA and as a 16mW 0.107mm2 ASIC @7nm chip. LCP models achieve a maximum and average speedups of 56x and 7x, compared to the originals, which could be improved by up to an average speedup of 33x by incorporating common optimizations such as pruning and quantization.


title: Optimal Stopping of Stochastic Transport Minimizing Submartingale Costs

abstract: Given a stochastic state process $(X_t)_t$ and a real-valued submartingale cost process $(S_t)_t$, we characterize optimal stopping times $\tau$ that minimize the expectation of $S_\tau$ while realizing given initial and target distributions $\mu$ and $\nu$, i.e., $X_0\sim \mu$ and $X_\tau \sim \nu$. A dual optimization problem is considered and shown to be attained under suitable conditions. The optimal solution of the dual problem then provides a contact set, which characterizes the location where optimal stopping can occur. The optimal stopping time is uniquely determined as the first hitting time of this contact set provided we assume a natural structural assumption on the pair $(X_t, S_t)_t$, which generalizes the twist condition on the cost in optimal transport theory. This paper extends the Brownian motion settings studied in [15, 16] and deals with more general costs.


title: Complex hyperbolic orbifolds and Lefschetz fibrations

abstract: A class of complex hyperbolic lattices in PU(2,1) called the Deligne-Mostow lattices has been reinterpreted by Hirzebruch and others in terms of line arrangements. They use branched covers over a suitable blow up of the complete quadrilateral arrangement of lines in projective 2-space to construct the complex hyperbolic surfaces over the orbifolds associated to the lattices. Fundamental domains for these lattices have been built by Pasquinelli. Here we show how the fundamental domains can be interpreted in terms of line arrangements as above. This parallel is then applied in two contexts. Dashyan uses Hirzebruch's construction to build infinitely many representations of 3-manifolds. Here we show that his construction can be generalised to all of the Deligne-Mostow lattices and more representations can be built. Wells shows that two of the Deligne-Mostow lattices in PU(2,1) can be seen as hybrids of lattices in PU(1,1). Here we show that he implicitly uses the line arrangement and we complete his analysis to all possible pairs of lines. In this way, we show that three more Deligne-Mostow lattices can be given as hybrids.


title: Measurement of the cross section for $\text{t}\bar{\text{t}}$ production with additional jets and b jets in pp collisions at $\sqrt{s} =$ 13 TeV

abstract: Measurements of the cross section for the production of top quark pairs in association with a pair of jets from bottom quarks ($\sigma_\mathrm{t\bar{t}b\bar{b}}$) and in association with a pair of jets from quarks of any flavor or gluons ($\sigma_\mathrm{t\bar{t}jj}$) and their ratio are presented. The data were collected in proton-proton collisions at a center-of-mass energy of 13 TeV by the CMS experiment at the LHC in 2016 and correspond to an integrated luminosity of 35.9 fb$^{-1}$. The measurements are performed in a fiducial phase space and extrapolated to the full phase space, separately in the dilepton and lepton+jets channels, where lepton corresponds to either an electron or a muon. The results for the dilepton and lepton+jets channels, respectively, are $\sigma_\mathrm{t\bar{t}jj}$ = 2.36 $\pm$ 0.02 (stat) $\pm$ 0.20 (syst) pb and 31.0 $\pm$ 0.2 (stat) $\pm$ 2.9 (syst) pb, and for the cross section ratio 0.017 $\pm$ 0.001 (stat) $\pm$ 0.001 (syst) and 0.020 $\pm$ 0.001 (stat) $\pm$ 0.001 (syst). The values of $\sigma_\mathrm{t\bar{t}b\bar{b}}$ are determined from the product of the $\sigma_\mathrm{t\bar{t}jj}$ and the cross section ratio, obtaining, respectively, 0.040 $\pm$ 0.002 (stat) $\pm$ 0.005 (syst) pb and 0.62 $\pm$ 0.03 (stat) $\pm$ 0.07 (syst) pb. These measurements are the most precise to date and are consistent, within the uncertainties, with the standard model expectations obtained using a matrix element calculation at next-to-leading order in quantum chromodynamics matched to a parton shower.


title: GeoDA: a geometric framework for black-box adversarial attacks

abstract: Adversarial examples are known as carefully perturbed images fooling image classifiers. We propose a geometric framework to generate adversarial examples in one of the most challenging black-box settings where the adversary can only generate a small number of queries, each of them returning the top-$1$ label of the classifier. Our framework is based on the observation that the decision boundary of deep networks usually has a small mean curvature in the vicinity of data samples. We propose an effective iterative algorithm to generate query-efficient black-box perturbations with small $\ell_p$ norms for $p \ge 1$, which is confirmed via experimental evaluations on state-of-the-art natural image classifiers. Moreover, for $p=2$, we theoretically show that our algorithm actually converges to the minimal $\ell_2$-perturbation when the curvature of the decision boundary is bounded. We also obtain the optimal distribution of the queries over the iterations of the algorithm. Finally, experimental results confirm that our principled black-box attack algorithm performs better than state-of-the-art algorithms as it generates smaller perturbations with a reduced number of queries.


title: Mean-field limit for a class of stochastic ergodic control problems

abstract: We study a family of McKean-Vlasov (mean-field) type ergodic optimal control problems with linear control, and quadratic dependence on control of the cost function. For this class of problems we establish existence and uniqueness of an optimal control. We propose an $N$-particles Markovian optimal control problem approximating the McKean-Vlasov one and we prove the convergence in relative entropy, total variation and Wasserstein distance of the law of the former to the law of the latter when $N$ goes to infinity. Some McKean-Vlasov optimal control problems with singular cost function and the relation of these problems with the mathematical theory of Bose-Einstein condensation is also established.


title: ${\mathbb Z}_2\times {\mathbb Z}_2$-graded mechanics: the classical theory

abstract: ${\mathbb Z}_2\times {\mathbb Z}_2$-graded mechanics admits four types of particles: ordinary bosons, two classes of fermions (fermions belonging to different classes commute among each other) and exotic bosons. In this paper we construct the basic ${\mathbb Z}_2\times {\mathbb Z}_2$-graded worldline multiplets (extending the cases of one-dimensional supersymmetry) and compute, based on a general scheme, their invariant classical actions and worldline sigma-models. The four basic multiplets contain two bosons and two fermions. They are $(2,2,0)$, with two propagating bosons and two propagating fermions, $(1,2,1)_{[00]}$ (the ordinary boson is propagating, while the exotic boson is an auxiliary field), $(1,2,1)_{[11]}$ (the converse case, the exotic boson is propagating, while the ordinary boson is an auxiliary field) and, finally, $(0,2,2)$ with two bosonic auxiliary fields. Classical actions invariant under the ${\mathbb Z}_2\times {\mathbb Z}_2$-graded superalgebra are constructed for both single multiplets and interacting multiplets. Furthermore, scale-invariant actions can possess a full ${\mathbb Z}_2\times {\mathbb Z}_2$-graded conformal invariance spanned by $10$ generators and containing an $sl(2)$ subalgebra.


title: DNN+NeuroSim V2.0: An End-to-End Benchmarking Framework for Compute-in-Memory Accelerators for On-chip Training

abstract: DNN+NeuroSim is an integrated framework to benchmark compute-in-memory (CIM) accelerators for deep neural networks, with hierarchical design options from device-level, to circuit-level and up to algorithm-level. A python wrapper is developed to interface NeuroSim with a popular machine learning platform: Pytorch, to support flexible network structures. The framework provides automatic algorithm-to-hardware mapping, and evaluates chip-level area, energy efficiency and throughput for training or inference, as well as training/inference accuracy with hardware constraints. Our prior work (DNN+NeuroSim V1.1) was developed to estimate the impact of reliability in synaptic devices, and analog-to-digital converter (ADC) quantization loss on the accuracy and hardware performance of inference engines. In this work, we further investigated the impact of the analog emerging non-volatile memory non-ideal device properties for on-chip training. By introducing the nonlinearity, asymmetry, device-to-device and cycle-to-cycle variation of weight update into the python wrapper, and peripheral circuits for error/weight gradient computation in NeuroSim core, we benchmarked CIM accelerators based on state-of-the-art SRAM and eNVM devices for VGG-8 on CIFAR-10 dataset, revealing the crucial specs of synaptic devices for on-chip training. The proposed DNN+NeuroSim V2.0 framework is available on GitHub.


title: Inducing Optimal Attribute Representations for Conditional GANs

abstract: Conditional GANs are widely used in translating an image from one category to another. Meaningful conditions to GANs provide greater flexibility and control over the nature of the target domain synthetic data. Existing conditional GANs commonly encode target domain label information as hard-coded categorical vectors in the form of 0s and 1s. The major drawbacks of such representations are inability to encode the high-order semantic information of target categories and their relative dependencies. We propose a novel end-to-end learning framework with Graph Convolutional Networks to learn the attribute representations to condition on the generator. The GAN losses, i.e. the discriminator and attribute classification losses, are fed back to the Graph resulting in the synthetic images that are more natural and clearer in attributes. Moreover, prior-arts are given priorities to condition on the generator side, not on the discriminator side of GANs. We apply the conditions to the discriminator side as well via multi-task learning. We enhanced the four state-of-the art cGANs architectures: Stargan, Stargan-JNT, AttGAN and STGAN. Our extensive qualitative and quantitative evaluations on challenging face attributes manipulation data set, CelebA, LFWA, and RaFD, show that the cGANs enhanced by our methods outperform by a large margin, compared to their counter-parts and other conditioning methods, in terms of both target attributes recognition rates and quality measures such as PSNR and SSIM.


title: Self-supervised Single-view 3D Reconstruction via Semantic Consistency

abstract: We learn a self-supervised, single-view 3D reconstruction model that predicts the 3D mesh shape, texture and camera pose of a target object with a collection of 2D images and silhouettes. The proposed method does not necessitate 3D supervision, manually annotated keypoints, multi-view images of an object or a prior 3D template. The key insight of our work is that objects can be represented as a collection of deformable parts, and each part is semantically coherent across different instances of the same category (e.g., wings on birds and wheels on cars). Therefore, by leveraging self-supervisedly learned part segmentation of a large collection of category-specific images, we can effectively enforce semantic consistency between the reconstructed meshes and the original images. This significantly reduces ambiguities during joint prediction of shape and camera pose of an object, along with texture. To the best of our knowledge, we are the first to try and solve the single-view reconstruction problem without a category-specific template mesh or semantic keypoints. Thus our model can easily generalize to various object categories without such labels, e.g., horses, penguins, etc. Through a variety of experiments on several categories of deformable and rigid objects, we demonstrate that our unsupervised method performs comparably if not better than existing category-specific reconstruction methods learned with supervision.


title: Optimizing Medical Treatment for Sepsis in Intensive Care: from Reinforcement Learning to Pre-Trial Evaluation

abstract: Our aim is to establish a framework where reinforcement learning (RL) of optimizing interventions retrospectively allows us a regulatory compliant pathway to prospective clinical testing of the learned policies in a clinical deployment. We focus on infections in intensive care units which are one of the major causes of death and difficult to treat because of the complex and opaque patient dynamics, and the clinically debated, highly-divergent set of intervention policies required by each individual patient, yet intensive care units are naturally data rich. In our work, we build on RL approaches in healthcare ("AI Clinicians"), and learn off-policy continuous dosing policy of pharmaceuticals for sepsis treatment using historical intensive care data under partially observable MDPs (POMDPs). POMPDs capture uncertainty in patient state better by taking in all historical information, yielding an efficient representation, which we investigate through ablations. We compensate for the lack of exploration in our retrospective data by evaluating each encountered state with a best-first tree search. We mitigate state distributional shift by optimizing our policy in the vicinity of the clinicians' compound policy. Crucially, we evaluate our model recommendations using not only conventional policy evaluations but a novel framework that incorporates human experts: a model-agnostic pre-clinical evaluation method to estimate the accuracy and uncertainty of clinician's decisions versus our system recommendations when confronted with the same individual patient history ("shadow mode").


title: Compressive Isogeometric Analysis

abstract: This work is motivated by the difficulty in assembling the Galerkin matrix when solving Partial Differential Equations (PDEs) with Isogeometric Analysis (IGA) using B-splines of moderate-to-high polynomial degree. To mitigate this problem, we propose a novel methodology named CossIGA (COmpreSSive IsoGeometric Analysis), which combines the IGA principle with CORSING, a recently introduced sparse recovery approach for PDEs based on compressive sensing. CossIGA assembles only a small portion of a suitable IGA Petrov-Galerkin discretization and is effective whenever the PDE solution is sufficiently sparse or compressible, i.e., when most of its coefficients are zero or negligible. The sparsity of the solution is promoted by employing a multilevel dictionary of B-splines as opposed to a basis. Thanks to sparsity and the fact that only a fraction of the full discretization matrix is assembled, the proposed technique has the potential to lead to significant computational savings. We show the effectiveness of CossIGA for the solution of the 2D and 3D Poisson equation over nontrivial geometries by means of an extensive numerical investigation.


title: Real-Time Monitoring of Area Angles with Synchrophasor Measurements

abstract: This paper develops a comprehensive framework of Area Angle Monitoring (AAM) to monitor the stress of bulk power transfer across an area of a power transmission system in real-time. Area angle is calculated from synchrophasor measurements to provide alert to system operators if the area angle exceeds pre-defined thresholds. This paper proposes general methods to identify these warning and emergency thresholds, and tests a mitigation strategy to relieve the area stress when the area angle exceeds the threshold. In order to handle the limited coverage of synchrophasor measurements, this paper proposes methods to estimate phase angles for boundary buses without synchrophasor measurements, which extends the application of AAM. AAM is verified for a power transmission area in the Western Electricity Coordinating Council system with both simulated data and synchrophasor measurements recorded from real events. A utility deployment to test the framework for monitoring area angle with live-stream and recorded synchrophasor data is described.


title: Web Performance with Android's Battery-Saver Mode

abstract: A Web browser utilizes a device's CPU to parse HTML, build a Document Object Model, a Cascading Style Sheets Object Model, and render trees, and parse, compile, and execute computationally-heavy JavaScript. A powerful CPU is required to perform these tasks as quickly as possible and provide the user with a great experience. However, increased CPU performance comes with increased power consumption and reduced battery life on mobile devices. As an option to extend battery life, Android offers a battery-saver mode that when activated, turns off the power-hungry and faster processor cores and turns on the battery-conserving and slower processor cores on the device. The transition from using faster processor cores to using slower processor cores throttles the CPU clock speed on the device, and therefore impacts the webpage load process. We utilize a large-scale data-set collected by a real user monitoring system of a major content delivery network to investigate the impact of Android's battery-saver mode on various mobile Web performance metrics. Our analysis suggests that users of select smartphones of Huawei and Sony experience a sudden or gradual degradation in Web performance when battery-saver mode is active. Battery-saver mode on newer flagship smartphones, however, does not impact the mobile Web performance. Finally, we encourage for new website design goals that treat slow (and throttled-CPU) devices kindly in favor of improving end-user experience and suggest that Web performance measurements should be aware of user device battery charge levels to correctly associate Web performance.


title: A Rate Splitting Strategy for Mitigating Intra-Cell Pilot Contamination in Massive MIMO

abstract: The spectral efficiency (SE) of Massive MIMO (MaMIMO) systems is affected by low quality channel estimates. Rate-Splitting (RS) has recently gained some interest in multiuser multiple antenna systems as an effective means to mitigate the multi-user interference due to imperfect channel state information. This paper investigates the benefits of RS in the downlink of a single-cell MaMIMO system when all the users use the same pilot sequence for channel estimation. Novel expressions for the SE achieved in the downlink by a single-layer RS strategy (that relies on a single successive interference cancellation at each user side) are derived and used to design precoding schemes and power allocation strategies for common and private messages. Numerical results are used to show that the proposed RS solution achieves higher spectral efficiency that conventional MaMIMO with maximum ratio precoding.


title: Absorbing time asymptotics in the oriented swap process

abstract: The oriented swap process is a natural directed random walk on the symmetric group that can be interpreted as a multi-species version of the Totally Asymmetric Simple Exclusion Process (TASEP) on a finite interval. An open problem from a 2008 paper of Angel, Holroyd, and Romik asks for the limiting distribution of the absorbing time of the process. We resolve this question by proving that this random variable satisfies GOE Tracy-Widom asymptotics. Our starting point is a distributional identity relating the behavior of the oriented swap process to last passage percolation, conjectured in a recent paper of Bisi, Cunden, Gibbons, and Romik. The main technical tool is a shift-invariance principle for multi-species TASEPs, obtained by exploiting recent results of Borodin, Gorin, and Wheeler for the stochastic colored six-vertex model.


title: Optoelectronic optimization of graded-bandgap thin-film AlGaAs solar cells

abstract: An optoelectronic optimization was carried out for an AlGaAs solar cell containing (i) an n-AlGaAs absorber layer with a graded bandgap and (ii) a periodically corrugated Ag backreflector combined with localized ohmic Pd-Ge-Au backcontacts. The bandgap of the absorber layer was varied either sinusoidally or linearly. An efficiency of 33.1% with the 2000-nm-thick n-AlGaAs absorber layer is predicted with linearly graded bandgap along with silver backreflector and localized ohmic backcontacts, in comparison to 27.4% efficiency obtained with homogeneous bandgap and a continuous ohmic backcontact. Sinusoidal grading of the bandgap {is predicted to enhance} the maximum efficiency to 34.5%. Thus, grading the bandgap of the absorber layer, along with a periodically corrugated Ag backreflector and localized ohmic Pd-Ge-Au backcontacts can help realize ultrathin and high-efficient AlGaAs solar cells for terrestrial applications.


title: The Cooperative Sorting Strategy for Connected and Automated Vehicle Platoons

abstract: This paper presents a "cooperative vehicle sorting" strategy that seeks to optimally sort connected and automated vehicles (CAVs) in a multi-lane platoon to reach an ideally organized platoon. In the proposed method, a CAV platoon is firstly discretized into a grid system, where a CAV moves from one cell to another in the discrete time-space domain. Then, the cooperative sorting problem is modeled as a path-finding problem in the graphic domain. The problem is solved by the deterministic Astar algorithm with a stepwise strategy, where only one vehicle can move within a movement step. The resultant shortest path is further optimized with an integer linear programming algorithm to minimize the sorting time by allowing multiple movements within a step. To improve the algorithm running time and address multiple shortest paths, a distributed stochastic Astar algorithm (DSA) is developed by introducing random disturbances to the edge costs to break uniform paths (with equal path cost). Numerical experiments are conducted to demonstrate the effectiveness of the proposed DSA method. The results report shorter sorting time and significantly improved algorithm running time due to the use of DSA. In addition, we find that the optimization performance can be further improved by increasing the number of processes in the distributed computing system.


title: Triangular resolutions and effectiveness for holomorphic subelliptic multipliers

abstract: A solution to the effectiveness problem in Kohn's algorithm for generating subelliptic multipliers is provided for domains that include those given by sums of squares of holomorphic functions (also including infinite sums). These domains are of particular interest due to their relation with complex and algebraic geometry and in particular, seem to include all previously known cases. Furthermore, combined with a recent result of M. Fassina, our effectiveness method allows establishing effective subelliptic estimates for more general classes of domains. Our main new tool, a triangular resolution, is the construction of subelliptic multipliers decomposable as $Q\circ\Gamma$, where $\Gamma$ is constructed from pre-multipliers and $Q$ is part of a triangular system. The effectiveness is proved via a sequence of newly proposed procedures, called here meta-procedures, built on top of the Kohn's procedures, where the order of subellipticity can be effectively tracked. Important sources of inspiration are algebraic-geometric techniques by Y.-T. Siu and procedures for triangular systems by D.W. Catlin and J.P. D'Angelo. The proposed procedures are purely algebraic and as such can also be of interest for geometric and computational problems involving Jacobian determinants, such as resolving singularities of holomorphic maps.


title: Realization of the kagome spin ice state in a frustrated intermetallic compound

abstract: Spin ices are exotic phases of matter characterized by frustrated spins obeying local ice rules, in analogy with the electric dipoles in water ice. In two dimensions, one can similarly define ice rules for in-plane Ising-like spins arranged on a kagome lattice. These ice rules require each triangle plaquette to have a single monopole, and can lead to various unique orders and excitations. Using experimental and theoretical approaches including magnetometry, thermodynamic measurements, neutron scattering and Monte Carlo simulations, we establish HoAgGe as a crystalline (i.e. non-artificial) system that realizes the kagome spin ice state. The system features a variety of partially and fully ordered states and a sequence of field-induced phases at low temperatures, all consistent with the kagome ice rule.


title: Toward fitting structured nonlinear systems by means of dynamic mode decomposition

abstract: The dynamic mode decomposition (DMD) is a data-driven method used for identifying the dynamics of complex nonlinear systems. It extracts important characteristics of the underlying dynamics using measured time-domain data produced either by means of experiments or by numerical simulations. In the original methodology, the measurements are assumed to be approximately related by a linear operator. Hence, a linear discrete-time system is fitted to the given data. However, often, nonlinear systems modeling physical phenomena have a particular known structure. In this contribution, we propose an identification and reduction method based on the classical DMD approach allowing to fit a structured nonlinear system to the measured data. We mainly focus on two types of nonlinearities: bilinear and quadratic-bilinear. By enforcing this additional structure, more insight into extracting the nonlinear behavior of the original process is gained. Finally, we demonstrate the proposed methodology for different examples, such as the Burgers' equation and the coupled van der Pol oscillators.


title: Robust Comparison in Population Protocols

abstract: There has recently been a surge of interest in the computational and complexity properties of the population model, which assumes $n$ anonymous, computationally-bounded nodes, interacting at random, and attempting to jointly compute global predicates. Significant work has gone towards investigating majority and consensus dynamics in this model: assuming that each node is initially in one of two states $X$ or $Y$, determine which state had higher initial count. In this paper, we consider a natural generalization of majority/consensus, which we call comparison. We are given two baseline states, $X_0$ and $Y_0$, present in any initial configuration in fixed, possibly small counts. Importantly, one of these states has higher count than the other: we will assume $|X_0| \ge C |Y_0|$ for some constant $C$. The challenge is to design a protocol which can quickly and reliably decide on which of the baseline states $X_0$ and $Y_0$ has higher initial count. We propose a simple algorithm solving comparison: the baseline algorithm uses $O(\log n)$ states per node, and converges in $O(\log n)$ (parallel) time, with high probability, to a state where whole population votes on opinions $X$ or $Y$ at rates proportional to initial $|X_0|$ vs. $|Y_0|$ concentrations. We then describe how such output can be then used to solve comparison. The algorithm is self-stabilizing, in the sense that it converges to the correct decision even if the relative counts of baseline states $X_0$ and $Y_0$ change dynamically during the execution, and leak-robust, in the sense that it can withstand spurious faulty reactions. Our analysis relies on a new martingale concentration result which relates the evolution of a population protocol to its expected (steady-state) analysis, which should be broadly applicable in the context of population protocols and opinion dynamics.


title: Recurrent convolutional neural networks for mandible segmentation from computed tomography

abstract: Recently, accurate mandible segmentation in CT scans based on deep learning methods has attracted much attention. However, there still exist two major challenges, namely, metal artifacts among mandibles and large variations in shape or size among individuals. To address these two challenges, we propose a recurrent segmentation convolutional neural network (RSegCNN) that embeds segmentation convolutional neural network (SegCNN) into the recurrent neural network (RNN) for robust and accurate segmentation of the mandible. Such a design of the system takes into account the similarity and continuity of the mandible shapes captured in adjacent image slices in CT scans. The RSegCNN infers the mandible information based on the recurrent structure with the embedded encoder-decoder segmentation (SegCNN) components. The recurrent structure guides the system to exploit relevant and important information from adjacent slices, while the SegCNN component focuses on the mandible shapes from a single CT slice. We conducted extensive experiments to evaluate the proposed RSegCNN on two head and neck CT datasets. The experimental results show that the RSegCNN is significantly better than the state-of-the-art models for accurate mandible segmentation.


title: Vector Topological Edge Solitons in Floquet Insulators

abstract: We introduce topological vector edge solitons in a Floquet insulator, consisting of two honeycomb arrays of helical waveguides with opposite directions of rotation in a focusing nonlinear optical medium. Zigzag edges of two arrays placed in contact create a zigzag-zigzag interface between two structures with different topology. A characteristic feature of such a photonic insulator is that, in the linear limit, it simultaneously supports two topologically protected chiral edge states at the interface between the two arrays. In the presence of nonlinearity, either bright or dark scalar Floquet edge soliton can bifurcate from a linear topological edge state. Such solitons are unidirectional and are localized in both directions, along the interface due to nonlinear self-action, and across the interface as being an edge state. The presence of two edge states with equal averaged group velocities enables the existence of stable topological vector edge solitons. In our case, these are nonlinearly coupled bright and dark solitons bifurcating from different branches of the topological Floquet edge states. Here we put forward a new mathematical description of scalar and vector small-amplitude Floquet envelope solitons in the above-mentioned continuous system. Importantly for the design of future photonic devices based on Floquet edge solitons, we find that the latter can be described by nonlinear Schrodinger equations for the mode envelopes obtained by averaging over one rotation period in the evolution coordinate.


title: Patch Graph Rewriting (Extended Version)

abstract: The basic principle of graph rewriting is the stepwise replacement of subgraphs inside a host graph. A challenge in such replacement steps is the treatment of the patch graph, consisting of those edges of the host graph that touch the subgraph, but are not part of it. We introduce the patch graph rewriting framework, a visual graph rewriting language with precise formal semantics. The language has rich expressive power in two ways. Firstly, rules can flexibly constrain the permitted shapes of patches touching matching subgraphs. Secondly, rules can freely transform patches. While the framework is designed to be easy to understand, it subsumes many approaches to graph rewriting.


title: Cut finite element error estimates for a class of nonlinear elliptic PDEs

abstract: Motivated by many applications in complex domains with boundaries exposed to large topological changes or deformations, fictitious domain methods regard the actual domain of interest as being embedded in a fixed Cartesian background. This is usually achieved via a geometric parameterization of its boundary via level-set functions. In this note, the a priori analysis of unfitted numerical schemes with cut elements is extended beyond the realm of linear problems. More precisely, we consider the discretization of semilinear elliptic boundary value problems of the form $- \Delta u +f_1(u)=f_2$ with polynomial nonlinearity via the cut finite element method. Boundary conditions are enforced, using a Nitsche-type approach. To ensure stability and error estimates that are independent of the position of the boundary with respect to the mesh, the formulations are augmented with additional boundary zone ghost penalty terms. These terms act on the jumps of the normal gradients at faces associated with cut elements. A-priori error estimates are derived, while numerical examples illustrate the implementation of the method and validate the theoretical findings.


title: Kummer surfaces for primality testing

abstract: We use the arithmetic of the Kummer surface associated to the Jacobian of a hyperelliptic curve to study the primality of integers of the form $4m^2 5^n-1$. We provide an algorithm capable of proving the primality or compositeness of most of the integers in these families and discuss in detail the necessary steps to implement this algorithm in a computer. Although an indetermination is possible, in which case another choice of initial parameters should be used, we prove that the probability of reaching this situation is exceedingly low and decreases exponentially with $n$.


title: Air-Transfer Production Method for Large-Area Picosecond Photodetectors

abstract: We have designed and prototyped the process steps for the batch production of large-area micro-channel-plate photomultipliers (MCP-PMT) using the "air-transfer" assembly process developed with single $LAPPD^{\text{TM}}$ modules. Results are presented addressing the challenges of designing a robust package that can transmit large numbers of electrical signals for pad or strip readout from inside the vacuum tube and hermetically sealing the large-perimeter window-body interface. We have also synthesized a photocathode in a large-area low-aspect-ratio volume, and shown that the micro-channel plates recover their functionality after cathode synthesis. The steps inform a design for a multi-module batch facility employing dual nested low-vacuum (LV) and ultra-high-vacuum (UHV) systems in a small-footprint. The facility design provides full access to multiple MCP-PMT modules prior to hermetic pinch-off for leak-checking and real-time photocathode optimization.


title: On Sufficient and Necessary Conditions in Bounded CTL: A Forgetting Approach

abstract: Computation Tree Logic (CTL) is one of the central formalisms in formal verification. As a specification language, it is used to express a property that the system at hand is expected to satisfy. From both the verification and the system design points of view, some information content of such property might become irrelevant for the system due to various reasons, e.g., it might become obsolete by time, or perhaps infeasible due to practical difficulties. Then, the problem arises on how to subtract such piece of information without altering the relevant system behaviour or violating the existing specifications over a given signature. Moreover, in such a scenario, two crucial notions are informative: the strongest necessary condition (SNC) and the weakest sufficient condition (WSC) of a given property. To address such a scenario in a principled way, we introduce a forgetting-based approach in CTL and show that it can be used to compute SNC and WSC of a property under a given model and over a given signature. We study its theoretical properties and also show that our notion of forgetting satisfies existing essential postulates of knowledge forgetting. Furthermore, we analyse the computational complexity of some basic reasoning tasks for the fragment CTL_AF in particular.


title: Distributed Stabilization of Two Interdependent Markov Jump Linear Systems with Partial Information

abstract: In this paper, we study the stabilization of two interdependent Markov jump linear systems (MJLSs) with partial information, where the interdependency arises as the transition of the mode of one system depends on the states of the other system. First, we formulate a framework for the two interdependent MJLSs to capture the interactions between various entities in the system, where the modes of the system cannot be observed directly. Instead, a signal which contains information of the modes can be obtained. Then, depending on the scope of the available system state information (global or local), we design centralized and distributed controllers, respectively, that can stochastically stabilize the overall interdependent MJLS. In addition, the sufficient stabilization conditions for the system under both types of information structure are derived. Finally, we provide a numerical example to illustrate the effectiveness of the designed controllers.


title: Technique for Generating Broadband FM Light

abstract: We suggest a technique for using off-resonance spectral comb generation to produce broadband frequency modulated, and therefore, amplitude quieted light. Results include closed-form formulae for the amplitudes and phases of all of the spectral components.


title: Human Preference-Based Learning for High-dimensional Optimization of Exoskeleton Walking Gaits

abstract: Optimizing lower-body exoskeleton walking gaits for user comfort requires understanding users' preferences over a high-dimensional gait parameter space. However, existing preference-based learning methods have only explored low-dimensional domains due to computational limitations. To learn user preferences in high dimensions, this work presents LineCoSpar, a human-in-the-loop preference-based framework that enables optimization over many parameters by iteratively exploring one-dimensional subspaces. Additionally, this work identifies gait attributes that characterize broader preferences across users. In simulations and human trials, we empirically verify that LineCoSpar is a sample-efficient approach for high-dimensional preference optimization. Our analysis of the experimental data reveals a correspondence between human preferences and objective measures of dynamicity, while also highlighting differences in the utility functions underlying individual users' gait preferences. This result has implications for exoskeleton gait synthesis, an active field with applications to clinical use and patient rehabilitation.


title: NSFnets (Navier-Stokes Flow nets): Physics-informed neural networks for the incompressible Navier-Stokes equations

abstract: We employ physics-informed neural networks (PINNs) to simulate the incompressible flows ranging from laminar to turbulent flows. We perform PINN simulations by considering two different formulations of the Navier-Stokes equations: the velocity-pressure (VP) formulation and the vorticity-velocity (VV) formulation. We refer to these specific PINNs for the Navier-Stokes flow nets as NSFnets. Analytical solutions and direct numerical simulation (DNS) databases provide proper initial and boundary conditions for the NSFnet simulations. The spatial and temporal coordinates are the inputs of the NSFnets, while the instantaneous velocity and pressure fields are the outputs for the VP-NSFnet, and the instantaneous velocity and vorticity fields are the outputs for the VV-NSFnet. These two different forms of the Navier-Stokes equations together with the initial and boundary conditions are embedded into the loss function of the PINNs. No data is provided for the pressure to the VP-NSFnet, which is a hidden state and is obtained via the incompressibility constraint without splitting the equations. We obtain good accuracy of the NSFnet simulation results upon convergence of the loss function, verifying that NSFnets can effectively simulate complex incompressible flows using either the VP or the VV formulations. We also perform a systematic study on the weights used in the loss function for the data/physics components and investigate a new way of computing the weights dynamically to accelerate training and enhance accuracy. Our results suggest that the accuracy of NSFnets, for both laminar and turbulent flows, can be improved with proper tuning of weights (manual or dynamic) in the loss function.


title: Deep Deterministic Portfolio Optimization

abstract: Can deep reinforcement learning algorithms be exploited as solvers for optimal trading strategies? The aim of this work is to test reinforcement learning algorithms on conceptually simple, but mathematically non-trivial, trading environments. The environments are chosen such that an optimal or close-to-optimal trading strategy is known. We study the deep deterministic policy gradient algorithm and show that such a reinforcement learning agent can successfully recover the essential features of the optimal trading strategies and achieve close-to-optimal rewards.


title: Explainable Deep Classification Models for Domain Generalization

abstract: Conventionally, AI models are thought to trade off explainability for lower accuracy. We develop a training strategy that not only leads to a more explainable AI system for object classification, but as a consequence, suffers no perceptible accuracy degradation. Explanations are defined as regions of visual evidence upon which a deep classification network makes a decision. This is represented in the form of a saliency map conveying how much each pixel contributed to the network's decision. Our training strategy enforces a periodic saliency-based feedback to encourage the model to focus on the image regions that directly correspond to the ground-truth object. We quantify explainability using an automated metric, and using human judgement. We propose explainability as a means for bridging the visual-semantic gap between different domains where model explanations are used as a means of disentagling domain specific information from otherwise relevant features. We demonstrate that this leads to improved generalization to new domains without hindering performance on the original domain.


title: LSCP: Enhanced Large Scale Colloquial Persian Language Understanding

abstract: Language recognition has been significantly advanced in recent years by means of modern machine learning methods such as deep learning and benchmarks with rich annotations. However, research is still limited in low-resource formal languages. This consists of a significant gap in describing the colloquial language especially for low-resourced ones such as Persian. In order to target this gap for low resource languages, we propose a "Large Scale Colloquial Persian Dataset" (LSCP). LSCP is hierarchically organized in a semantic taxonomy that focuses on multi-task informal Persian language understanding as a comprehensive problem. This encompasses the recognition of multiple semantic aspects in the human-level sentences, which naturally captures from the real-world sentences. We believe that further investigations and processing, as well as the application of novel algorithms and methods, can strengthen enriching computerized understanding and processing of low resource languages. The proposed corpus consists of 120M sentences resulted from 27M tweets annotated with parsing tree, part-of-speech tags, sentiment polarity and translation in five different languages.


title: An R Autograder for PrairieLearn

abstract: We describe how we both use and extend the PrarieLearn framework by taking advantage of its built-in support for external auto-graders. By using a custom Docker container, we can match our course requirements perfectly. Moreover, by relying on the flexibility of the interface we can customize our Docker container. A specific extension for unit testing is described which creates context-dependent difference between student answers and reference solution providing a more comprehensive response at test time.


title: Polynomial potentials and coupled quantum dots in two and three dimensions

abstract: Non-separable $D-$dimensional partial differential Schr\"{o}dinger equations are considered at $D=2$ and $D=3$, with the even-parity local potentials $V(x,y,\ldots)$ which are polynomials of degree four (cusp catastrophe resembling case) and six (butterfly resembling case). Their extremes (i.e., minima and maxima) are assumed pronounced, localized via a suitable ad hoc parametrization of the coupling constants. A non-numerical approximate construction of the low lying bound states $\psi(x,y,\ldots)$] is then found feasible in the dynamical regime simulating a coupled system of quantum dots in which the individual minima of $V(x,y,\ldots)$ are well separated, with the potential being locally approximated by the harmonic oscillator wells. The measurable characteristics (and, in particular, the topologically protected probability-density distributions) are then found bifurcating in a specific evolution scenario called a relocalization quantum catastrophe.


title: Ergodic Theorems for Dynamic Imprecise Probability Kinematics

abstract: We formulate an ergodic theory for the (almost sure) limit $\mathcal{P}^\text{co}_{\tilde{\mathcal{E}}}$ of a sequence $(\mathcal{P}^\text{co}_{\mathcal{E}_n})$ of successive dynamic imprecise probability kinematics (DIPK, introduced in Caprio and Gong, 2021) updates of a set $\mathcal{P}^\text{co}_{\mathcal{E}_0}$ representing the initial beliefs of an agent. As a consequence, we formulate a strong law of large numbers.


title: Decomposability and co-modular indices of tournaments

abstract: Given a tournament $T$, a module of $T$ is a subset $X$ of $V(T)$ such that for $x, y\in X$ and $v\in V(T)\setminus X$, $(x,v)\in A(T)$ if and only if $(y,v)\in A(T)$. The trivial modules of $T$ are $\emptyset$, $\{u\}$ $(u\in V(T))$ and $V(T)$. The tournament $T$ is indecomposable if all its modules are trivial; otherwise it is decomposable. The decomposability index of $T$, denoted by $\delta(T)$, is the smallest number of arcs of $T$ that must be reversed to make $T$ indecomposable. The first author conjectured that for $n \geq 5$, we have $\delta(n) = \left\lceil \frac{n+1}{4} \right\rceil$, where $\delta(n)$ is the maximum of $\delta(T)$ over the tournaments $T$ with $n$ vertices. In this paper we prove this conjecture by introducing the co-modular index of a tournament $T$, denoted by $\Delta(T)$, as the largest number of disjoint co-modules of $T$, where a co-module of $T$ is a subset $M$ of $V(T)$ such that $M$ or $V(T) \setminus M$ is a nontrivial module of $T$. We prove that for $n \geq 3$, we have $\Delta(n) = \left\lceil \frac{n+1}{2} \right\rceil$, where $\Delta(n)$ is the maximum of $\Delta(T)$ over the tournaments $T$ with $n$ vertices. Our main result is the following close relationship between the above two indices: for every tournament $T$ with at least $5$ vertices, we have $\delta(T) = \left\lceil \frac{\Delta(T)}{2} \right\rceil$. As a consequence, we obtain $\delta(n) = \left\lceil \frac{\Delta(n)}{2} \right\rceil = \left\lceil \frac{n+1}{4} \right\rceil$ for $n \geq 5$, and we answer some further related questions.


title: One-Shot Parameter Identification of the Thevenin's Model for Batteries: Methods and Validation

abstract: Parameter estimation is of foundational importance for various model-based battery management tasks, including charging control, state-of-charge estimation and aging assessment. However, it remains a challenging issue as the existing methods generally depend on cumbersome and time-consuming procedures to extract battery parameters from data. Departing from the literature, this paper sets the unique aim of identifying all the parameters offline in a one-shot procedure, including the resistance and capacitance parameters and the parameters in the parameterized function mapping from the state-of-charge to the open-circuit voltage. Considering the well-known Thevenin's battery model, the study begins with the parameter identifiability analysis, showing that all the parameters are locally identifiable. Then, it formulates the parameter identification problem in a prediction-error-minimization framework. As the non-convexity intrinsic to the problem may lead to physically meaningless estimates, two methods are developed to overcome this issue. The first one is to constrain the parameter search within a reasonable space by setting parameter bounds, and the other adopts regularization of the cost function using prior parameter guess. The proposed identifiability analysis and identification methods are extensively validated through simulations and experiments.


title: AutoGluon-Tabular: Robust and Accurate AutoML for Structured Data

abstract: We introduce AutoGluon-Tabular, an open-source AutoML framework that requires only a single line of Python to train highly accurate machine learning models on an unprocessed tabular dataset such as a CSV file. Unlike existing AutoML frameworks that primarily focus on model/hyperparameter selection, AutoGluon-Tabular succeeds by ensembling multiple models and stacking them in multiple layers. Experiments reveal that our multi-layer combination of many models offers better use of allocated training time than seeking out the best. A second contribution is an extensive evaluation of public and commercial AutoML platforms including TPOT, H2O, AutoWEKA, auto-sklearn, AutoGluon, and Google AutoML Tables. Tests on a suite of 50 classification and regression tasks from Kaggle and the OpenML AutoML Benchmark reveal that AutoGluon is faster, more robust, and much more accurate. We find that AutoGluon often even outperforms the best-in-hindsight combination of all of its competitors. In two popular Kaggle competitions, AutoGluon beat 99% of the participating data scientists after merely 4h of training on the raw data.