title: Accurate simulation of operating system updates in neuroimaging using Monte-Carlo arithmetic

abstract: Operating system (OS) updates introduce numerical perturbations that impact the reproducibility of computational pipelines. In neuroimaging, this has important practical implications on the validity of computational results, particularly when obtained in systems such as high-performance computing clusters where the experimenter does not control software updates. We present a framework to reproduce the variability induced by OS updates in controlled conditions. We hypothesize that OS updates impact computational pipelines mainly through numerical perturbations originating in mathematical libraries, which we simulate using Monte-Carlo arithmetic in a framework called "fuzzy libmath" (FL). We applied this methodology to pre-processing pipelines of the Human Connectome Project, a flagship open-data project in neuroimaging. We found that FL-perturbed pipelines accurately reproduce the variability induced by OS updates and that this similarity is only mildly dependent on simulation parameters. Importantly, we also found between-subject differences were preserved in both cases, though the between-run variability was of comparable magnitude for both FL and OS perturbations. We found the numerical precision in the HCP pre-processed images to be relatively low, with less than 8 significant bits among the 24 available, which motivates further investigation of the numerical stability of components in the tested pipeline. Overall, our results establish that FL accurately simulates results variability due to OS updates, and is a practical framework to quantify numerical uncertainty in neuroimaging.


title: Complex-valued Spatial Autoencoders for Multichannel Speech Enhancement

abstract: In this contribution, we present a novel online approach to multichannel speech enhancement. The proposed method estimates the enhanced signal through a filter-and-sum framework. More specifically, complex-valued masks are estimated by a deep complex-valued neural network, termed the complex-valued spatial autoencoder. The proposed network is capable of exploiting as well as manipulating both the phase and the amplitude of the microphone signals. As shown by the experimental results, the proposed approach is able to exploit both spatial and spectral characteristics of the desired source signal resulting in a physically plausible spatial selectivity and superior speech quality compared to other baseline methods.


title: COVID-Net US: A Tailored, Highly Efficient, Self-Attention Deep Convolutional Neural Network Design for Detection of COVID-19 Patient Cases from Point-of-care Ultrasound Imaging

abstract: The Coronavirus Disease 2019 (COVID-19) pandemic has impacted many aspects of life globally, and a critical factor in mitigating its effects is screening individuals for infections, thereby allowing for both proper treatment for those individuals as well as action to be taken to prevent further spread of the virus. Point-of-care ultrasound (POCUS) imaging has been proposed as a screening tool as it is a much cheaper and easier to apply imaging modality than others that are traditionally used for pulmonary examinations, namely chest x-ray and computed tomography. Given the scarcity of expert radiologists for interpreting POCUS examinations in many highly affected regions around the world, low-cost deep learning-driven clinical decision support solutions can have a large impact during the on-going pandemic. Motivated by this, we introduce COVID-Net US, a highly efficient, self-attention deep convolutional neural network design tailored for COVID-19 screening from lung POCUS images. Experimental results show that the proposed COVID-Net US can achieve an AUC of over 0.98 while achieving 353X lower architectural complexity, 62X lower computational complexity, and 14.3X faster inference times on a Raspberry Pi. Clinical validation was also conducted, where select cases were reviewed and reported on by a practicing clinician (20 years of clinical practice) specializing in intensive care (ICU) and 15 years of expertise in POCUS interpretation. To advocate affordable healthcare and artificial intelligence for resource-constrained environments, we have made COVID-Net US open source and publicly available as part of the COVID-Net open source initiative.


title: RockGPT: Reconstructing three-dimensional digital rocks from single two-dimensional slice from the perspective of video generation

abstract: Random reconstruction of three-dimensional (3D) digital rocks from two-dimensional (2D) slices is crucial for elucidating the microstructure of rocks and its effects on pore-scale flow in terms of numerical modeling, since massive samples are usually required to handle intrinsic uncertainties. Despite remarkable advances achieved by traditional process-based methods, statistical approaches and recently famous deep learning-based models, few works have focused on producing several kinds of rocks with one trained model and allowing the reconstructed samples to satisfy certain given properties, such as porosity. To fill this gap, we propose a new framework, named RockGPT, which is composed of VQ-VAE and conditional GPT, to synthesize 3D samples based on a single 2D slice from the perspective of video generation. The VQ-VAE is utilized to compress high-dimensional input video, i.e., the sequence of continuous rock slices, to discrete latent codes and reconstruct them. In order to obtain diverse reconstructions, the discrete latent codes are modeled using conditional GPT in an autoregressive manner, while incorporating conditional information from a given slice, rock type, and porosity. We conduct two experiments on five kinds of rocks, and the results demonstrate that RockGPT can produce different kinds of rocks with the same model, and the reconstructed samples can successfully meet certain specified porosities. In a broader sense, through leveraging the proposed conditioning scheme, RockGPT constitutes an effective way to build a general model to produce multiple kinds of rocks simultaneously that also satisfy user-defined properties.


title: Critical Droplets and sharp asymptotics for Kawasaki dynamics with strongly anisotropic interactions

abstract: In this paper we analyze metastability and nucleation in the context of the Kawasaki dynamics for the two-dimensional Ising lattice gas at very low temperature. Let $\Lambda\subset\mathbb{Z}^2$ be a finite box. Particles perform simple exclusion on $\Lambda$, but when they occupy neighboring sites they feel a binding energy $-U_1<0$ in the horizontal direction and $-U_2<0$ in the vertical one. Thus the Kawasaki dynamics is conservative inside the volume $\Lambda$. Along each bond touching the boundary of $\Lambda$ from the outside to the inside, particles are created with rate $\rho=e^{-\Delta\beta}$, while along each bond from the inside to the outside, particles are annihilated with rate $1$, where $\beta>0$ is the inverse temperature and $\Delta>0$ is an activity parameter. Thus, the boundary of $\Lambda$ plays the role of an infinite gas reservoir with density $\rho$. We consider the parameter regime $U_1>2U_2$ also known as the strongly anisotropic regime. We take $\Delta\in{(U_1,U_1+U_2)}$, so that the empty (respectively full) configuration is a metastable (respectively stable) configuration. We investigate how the transition from empty to full takes place with particular attention to the critical configurations that asymptotically have to be crossed with probability 1. The derivation of some geometrical properties of the saddles allows us to identify the full geometry of the minimal gates and their boundaries for the nucleation in the strongly anisotropic case. We observe very different behaviors for this case with respect to the isotropic ($U_1=U_2$) and weakly anisotropic ($U_1<2U_2$) ones. Moreover, we derive mixing time, spectral gap and sharp estimates for the asymptotic transition time for the strongly anisotropic case.


title: Improved $\Lambda p$ Elastic Scattering Cross Sections Between 0.9 and 2.0 GeV/c and Connections to the Neutron Star Equation of State

abstract: Strange matter is believed to exist in the cores of neutron stars based on simple kinematics. If this is true, then hyperon-nucleon interactions will play a significant part in the neutron star equation of state (EOS). Yet, compared to other elastic scattering processes, there is very little data on $\Lambda$-$N$ scattering. This experiment utilized the CLAS detector to study the $\Lambda p \rightarrow \Lambda p$ elastic scattering cross section in the incident $\Lambda$ momentum range 0.9-2.0 GeV/c. This is the first data on this reaction in several decades. The new cross sections have significantly better accuracy and precision than the existing world data, and the techniques developed here can also be used in future experiments.


title: Minimax Boundary Estimation and Estimation with Boundary

abstract: We derive non-asymptotic minimax bounds for the Hausdorff estimation of $d$-dimensional submanifolds $M \subset \mathbb{R}^D$ with (possibly) non-empty boundary $\partial M$. The model reunites and extends the most prevalent $\mathcal{C}^2$-type set estimation models: manifolds without boundary, and full-dimensional domains. We consider both the estimation of the manifold $M$ itself and that of its boundary $\partial M$ if non-empty. Given $n$ samples, the minimax rates are of order $O\bigl((\log n/n)^{2/d}\bigr)$ if $\partial M = \emptyset$ and $O\bigl((\log n/n)^{2/(d+1)}\bigr)$ if $\partial M \neq \emptyset$, up to logarithmic factors. In the process, we develop a Voronoi-based procedure that allows to identify enough points $O\bigl((\log n/n)^{2/(d+1)}\bigr)$-close to $\partial M$ for reconstructing it.


title: Resource-efficient dissipative entanglement of two trapped-ion qubits

abstract: We demonstrate a simplified method for dissipative generation of an entangled state of two trapped-ion qubits. Our implementation produces its target state faster and with higher fidelity than previous demonstrations of dissipative entanglement generation and eliminates the need for auxiliary ions. The entangled singlet state is generated in $\sim$7 ms with a fidelity of 0.949(4). The dominant source of infidelity is photon scattering. We discuss this error source and strategies for its mitigation.


title: Extendibility limits the performance of quantum processors

abstract: Resource theories in quantum information science are helpful for the study and quantification of the performance of information-processing tasks that involve quantum systems. These resource theories also find applications in other areas of study; e.g., the resource theories of entanglement and coherence have found use and implications in the study of quantum thermodynamics and memory effects in quantum dynamics. In this paper, we introduce the resource theory of unextendibility, which is associated to the inability of extending quantum entanglement in a given quantum state to multiple parties. The free states in this resource theory are the k-extendible states, and the free channels are k-extendible channels, which preserve the class of k-extendible states. We make use of this resource theory to derive non-asymptotic, upper bounds on the rate at which quantum communication or entanglement preservation is possible by utilizing an arbitrary quantum channel a finite number of times, along with the assistance of k-extendible channels at no cost. We then show that the bounds obtained are significantly tighter than previously known bounds for quantum communication over both the depolarizing and erasure channels.


title: Lung Ultrasound Segmentation and Adaptation between COVID-19 and Community-Acquired Pneumonia

abstract: Lung ultrasound imaging has been shown effective in detecting typical patterns for interstitial pneumonia, as a point-of-care tool for both patients with COVID-19 and other community-acquired pneumonia (CAP). In this work, we focus on the hyperechoic B-line segmentation task. Using deep neural networks, we automatically outline the regions that are indicative of pathology-sensitive artifacts and their associated sonographic patterns. With a real-world data-scarce scenario, we investigate approaches to utilize both COVID-19 and CAP lung ultrasound data to train the networks; comparing fine-tuning and unsupervised domain adaptation. Segmenting either type of lung condition at inference may support a range of clinical applications during evolving epidemic stages, but also demonstrates value in resource-constrained clinical scenarios. Adapting real clinical data acquired from COVID-19 patients to those from CAP patients significantly improved Dice scores from 0.60 to 0.87 (p < 0.001) and from 0.43 to 0.71 (p < 0.001), on independent COVID-19 and CAP test cases, respectively. It is of practical value that the improvement was demonstrated with only a small amount of data in both training and adaptation data sets, a common constraint for deploying machine learning models in clinical practice. Interestingly, we also report that the inverse adaptation, from labelled CAP data to unlabeled COVID-19 data, did not demonstrate an improvement when tested on either condition. Furthermore, we offer a possible explanation that correlates the segmentation performance to label consistency and data domain diversity in this point-of-care lung ultrasound application.


title: Explicit characterisation of the fractional power spaces of the Dirichlet Laplacian and Stokes operators

abstract: We identify explicitly the fractional power spaces for the $L^2$ Dirichlet Laplacian and Dirichlet Stokes operators using the theory of real interpolation. The results are not new, but we hope that our arguments are relatively accessible.


title: SELM: Siamese Extreme Learning Machine with Application to Face Biometrics

abstract: Extreme Learning Machine is a powerful classification method very competitive existing classification methods. It is extremely fast at training. Nevertheless, it cannot perform face verification tasks properly because face verification tasks require comparison of facial images of two individuals at the same time and decide whether the two faces identify the same person. The structure of Extreme Leaning Machine was not designed to feed two input data streams simultaneously, thus, in 2-input scenarios Extreme Learning Machine methods are normally applied using concatenated inputs. However, this setup consumes two times more computational resources and it is not optimized for recognition tasks where learning a separable distance metric is critical. For these reasons, we propose and develop a Siamese Extreme Learning Machine (SELM). SELM was designed to be fed with two data streams in parallel simultaneously. It utilizes a dual-stream Siamese condition in the extra Siamese layer to transform the data before passing it along to the hidden layer. Moreover, we propose a Gender-Ethnicity-Dependent triplet feature exclusively trained on a variety of specific demographic groups. This feature enables learning and extracting of useful facial features of each group. Experiments were conducted to evaluate and compare the performances of SELM, Extreme Learning Machine, and DCNN. The experimental results showed that the proposed feature was able to perform correct classification at 97.87% accuracy and 99.45% AUC. They also showed that using SELM in conjunction with the proposed feature provided 98.31% accuracy and 99.72% AUC. They outperformed the well-known DCNN and Extreme Leaning Machine methods by a wide margin.


title: Computing solution landscape of nonlinear space-fractional problems via fast approximation algorithm

abstract: The nonlinear space-fractional problems often allow multiple stationary solutions, which can be much more complicated than the corresponding integer-order problems. In this paper, we systematically compute the solution landscapes of nonlinear constant/variable-order space-fractional problems. A fast approximation algorithm is developed to deal with the variable-order spectral fractional Laplacian by approximating the variable-indexing Fourier modes, and then combined with saddle dynamics to construct the solution landscape of variable-order space-fractional phase field model. Numerical experiments are performed to substantiate the accuracy and efficiency of fast approximation algorithm and elucidate essential features of the stationary solutions of space-fractional phase field model. Furthermore, we demonstrate that the solution landscapes of spectral fractional Laplacian problems can be reconfigured by varying the diffusion coefficients in the corresponding integer-order problems.


title: Particle fluctuations and the failure of simple effective models for many-body localized phases

abstract: We investigate and compare the particle number fluctuations in the putative many-body localized (MBL) phase of a spinless fermion model with potential disorder and nearest-neighbor interactions with those in the non-interacting case (Anderson localization) and in effective models where only interaction terms diagonal in the Anderson basis are kept. We demonstrate that these types of simple effective models cannot account for the particle number fluctuations observed in the MBL phase of the microscopic model. This implies that assisted and pair hopping terms---generated when transforming the microscopic Hamiltonian into the Anderson basis---cannot be neglected. As a consequence, it appears questionable if the microscopic model possesses an exponential number of exactly conserved local charges. If such exactly conserved local charges do not exist, then particles are expected to ultimately delocalize for any finite disorder strength.


title: An Integrated Progressive Hedging and Benders Decomposition with Multiple Master Method to Solve the Brazilian Generation Expansion Problem

abstract: This paper exploits the decomposition structure of the large-scale hydrothermal generation expansion planning problem with an integrated modified Benders Decomposition and Progressive Hedging approach. We consider detailed and realistic data from the Brazilian power system to represent hourly chronological constraints based on typical days per month and year. Also, we represent the multistage stochastic nature of the optimal hydrothermal operational policy through co-optimized linear decision rules for individual reservoirs. Therefore, we ensure investment decisions compatible with a nonanticipative (implementable) operational policy. To solve the large-scale optimization problem, we propose an improved Benders Decomposition method with multiple instances of the master problem, each of which strengthened by primal cuts and new Benders cuts generated by each master's trial solution. Additionally, our new approach allows using Progressive Hedging penalization terms for accelerating the convergence of the method. We show that our method is 60\% faster than the benchmark. Finally, the consideration of a nonanticipative operational policy can save 7.64\% of the total cost (16.18\% of the investment costs) and significantly improve spot price profiles.


title: ELSED: Enhanced Line SEgment Drawing

abstract: Detecting local features, such as corners, segments or blobs, is the first step in the pipeline of many Computer Vision applications. Its speed is crucial for real-time applications. In this paper we present ELSED, the fastest line segment detector in the literature. The key for its efficiency is a local segment growing algorithm that connects gradient-aligned pixels in presence of small discontinuities. The proposed algorithm not only runs in devices with very low end hardware, but may also be parametrized to foster the detection of short or longer segments, depending on the task at hand. We also introduce new metrics to evaluate the accuracy and repeatability of segment detectors. In our experiments with different public benchmarks we prove that our method accounts the highest repeatability and it is the most efficient in the literature. In the experiments we quantify the accuracy traded for such gain.


title: The physics of cement cohesion

abstract: Cement is one of the most produced materials in the world. A major player in greenhouse gas emissions, it is the main binding agent in concrete, to which it provides a cohesive strength that rapidly increases during setting. Understanding how such cohesion emerges has been a major obstacle to advances in cement science and technology. Here, we combine computational statistical mechanics and theory to demonstrate how cement cohesion results from the organization of interlocked ions and water, progressively confined in nano-slits between charged surfaces of Calcium-Silicate-Hydrates. Due to the water/ions interlocking, dielectric screening is drastically reduced and ionic correlations are proven significantly stronger than previously thought, dictating the evolution of the nano-scale interactions during cement hydration. By developing a quantitative analytical prediction of cement cohesion based on Coulombic forces, we reconcile a novel fundamental understanding of cement hydration with the fully atomistic description of the solid cement paste and open new paths for science and technologies of construction materials.

categories: physics.app-ph cond-mat.mtrl-sci cond-mat.stat-mech physics.comp-ph


title: Topology Optimization Methods for 3D Structural Problems: A Comparative Study

abstract: The work provides an exhaustive comparison of some representative families of topology optimization methods for 3D structural optimization, such as the Solid Isotropic Material with Penalization (SIMP), the Level-set, the Bidirectional Evolutionary Structural Optimization (BESO), and the Variational Topology Optimization (VARTOP) methods. The main differences and similarities of these approaches are then highlighted from an algorithmic standpoint. The comparison is carried out via the study of a set of numerical benchmark cases using industrial-like fine-discretization meshes (around 1 million finite elements), and Matlab as the common computational platform, to ensure fair comparisons. Then, the results obtained for every benchmark case with the different methods are compared in terms of computational cost, topology quality, achieved minimum value of the objective function, and robustness of the computations (convergence in objective function and topology). Finally, some quantitative and qualitative results are presented, from which, an attempt of qualification of the methods, in terms of their relative performance, is done.


title: Charge-velocity-dependent one-scale linear model

abstract: We apply a recently developed formalism to study the evolution of a current-carrying string network under the simple but generic assumption of a linear equation of state. We demonstrate that the existence of a scaling solution with non-trivial current depends on the expansion rate of the universe, the initial root mean square current on the string, and the available energy loss mechanisms. We find that the fast expansion rate after radiation-matter equality will tend to rapidly dilute any pre-existing current and the network will evolve towards the standard Nambu-Goto scaling solution (provided there are no external current-generating mechanisms). During the radiation era, current growth is possible provided the initial conditions for the network generate a relatively large current and/or there is significant early string damping. The network can then achieve scaling with a stable non-trivial current, assuming large currents will be regulated by some leakage mechanism. The potential existence of current-carrying string networks in the radiation era, unlike the standard Nambu-Goto networks expected in the matter era, could have interesting phenomenological consequences.


title: Nonlinear radiation gauge for near Kerr spacetimes

abstract: In this paper, we introduce and explore the properties of a new gauge choice for the vacuum Einstein equation inspired by the ingoing and outgoing radiation gauges (IRG, ORG) for the linearized vacuum Einstein equation introduced by Chrzanowski in his work on metric reconstruction [DOI:10.1103/PhysRevD.11.2042] on the Kerr background. It has been shown by Price, Shankar and Whiting [arXiv:gr-qc/0611070] that the IRG/ORG are consistent gauges for the linearized vacuum Einstein equation on Petrov type II backgrounds. In [arXiv:1903.03859], the ORG was used in proving linearized stability for the Kerr spacetime, and the new non-linear radiation gauge introduced here is a direct generalization of that gauge condition, and is intended to be used to study the stability of Kerr black holes under the evolution generated by the vacuum Einstein equation.


title: On K-K-W Type Theorems for Conformal Perturbations of Twisted Dirac Operators

abstract: In this paper, we prove two Kastler-Kalau-Walze type theorems for conformal perturbations of twisted Dirac operators and conformal perturbations of twisted signature operators on four-dimensional manifolds with (resp. without) boundary.


title: Attainment Regions in Feature-Parameter Space for High-Level Debugging in Autonomous Robots

abstract: Understanding a controller's performance in different scenarios is crucial for robots that are going to be deployed in safety-critical tasks. If we do not have a model of the dynamics of the world, which is often the case in complex domains, we may need to approximate a performance function of the robot based on its interaction with the environment. Such a performance function gives us insights into the behaviour of the robot, allowing us to fine-tune the controller with manual interventions. In high-dimensionality systems, where the actionstate space is large, fine-tuning a controller is non-trivial. To overcome this problem, we propose a performance function whose domain is defined by external features and parameters of the controller. Attainment regions are defined over such a domain defined by feature-parameter pairs, and serve the purpose of enabling prediction of successful execution of the task. The use of the feature-parameter space -in contrast to the action-state space- allows us to adapt, explain and finetune the controller over a simpler (i.e., lower dimensional space). When the robot successfully executes the task, we use the attainment regions to gain insights into the limits of the controller, and its robustness. When the robot fails to execute the task, we use the regions to debug the controller and find adaptive and counterfactual changes to the solutions. Another advantage of this approach is that we can generalise through the use of Gaussian processes regression of the performance function in the high-dimensional space. To test our approach, we demonstrate learning an approximation to the performance function in simulation, with a mobile robot traversing different terrain conditions. Then, with a sample-efficient method, we propagate the attainment regions to a physical robot in a similar environment.


title: Full-Duplex Strategy for Video Object Segmentation

abstract: Previous video object segmentation approaches mainly focus on using simplex solutions between appearance and motion, limiting feature collaboration efficiency among and across these two cues. In this work, we study a novel and efficient full-duplex strategy network (FSNet) to address this issue, by considering a better mutual restraint scheme between motion and appearance in exploiting the cross-modal features from the fusion and decoding stage. Specifically, we introduce the relational cross-attention module (RCAM) to achieve bidirectional message propagation across embedding sub-spaces. To improve the model's robustness and update the inconsistent features from the spatial-temporal embeddings, we adopt the bidirectional purification module (BPM) after the RCAM. Extensive experiments on five popular benchmarks show that our FSNet is robust to various challenging scenarios (e.g., motion blur, occlusion) and achieves favourable performance against existing cutting-edges both in the video object segmentation and video salient object detection tasks. The project is publicly available at: https://dpfan.net/FSNet.


title: Source-Free Domain Adaptation for Image Segmentation

abstract: Domain adaptation (DA) has drawn high interest for its capacity to adapt a model trained on labeled source data to perform well on unlabeled or weakly labeled target data from a different domain. Most common DA techniques require concurrent access to the input images of both the source and target domains. However, in practice, privacy concerns often impede the availability of source images in the adaptation phase. This is a very frequent DA scenario in medical imaging, where, for instance, the source and target images could come from different clinical sites. We introduce a source-free domain adaptation for image segmentation. Our formulation is based on minimizing a label-free entropy loss defined over target-domain data, which we further guide with a domain-invariant prior on the segmentation regions. Many priors can be derived from anatomical information. Here, a class ratio prior is estimated from anatomical knowledge and integrated in the form of a Kullback Leibler (KL) divergence in our overall loss function. Furthermore, we motivate our overall loss with an interesting link to maximizing the mutual information between the target images and their label predictions. We show the effectiveness of our prior aware entropy minimization in a variety of domain-adaptation scenarios, with different modalities and applications, including spine, prostate, and cardiac segmentation. Our method yields comparable results to several state of the art adaptation techniques, despite having access to much less information, as the source images are entirely absent in our adaptation phase. Our straightforward adaptation strategy uses only one network, contrary to popular adversarial techniques, which are not applicable to a source-free DA setting. Our framework can be readily used in a breadth of segmentation problems, and our code is publicly available: https://github.com/mathilde-b/SFDA


title: Flat cotorsion modules over Noether algebras

abstract: For a module-finite algebra over a commutative noetherian ring, we give a complete description of flat cotorsion modules in terms of prime ideals of the algebra, as a generalization of Enochs' result for a commutative noetherian ring. As a consequence, we show that pointwise Matlis duality gives a bijective correspondence between the isoclasses of indecomposable injective left modules and the isoclasses of indecomposable flat cotorsion right modules. This correspondence is an explicit realization of Herzog's homeomorphism induced from elementary duality of Ziegler spectra.


title: Independence Properties of Generalized Submodular Information Measures

abstract: Recently a class of generalized information measures was defined on sets of items parametrized by submodular functions. In this paper, we propose and study various notions of independence between sets with respect to such information measures, and connections thereof. Since entropy can also be used to parametrize such measures, we derive interesting independence properties for the entropy of sets of random variables. We also study the notion of multi-set independence and its properties. Finally, we present optimization algorithms for obtaining a set that is independent of another given set, and also discuss the implications and applications of combinatorial independence.


title: Poincar\'e index formula and analogy with the Kosterlitz-Thouless transition in a non-rotated cold atom Bose-Einstein condensate

abstract: A dilute gas of Bose-Einstein condensed atoms in a non-rotated and axially symmetric harmonic trap is modelled by the time dependent Gross-Pitaevskii equation. When the angular momentum carried by the condensate does not vanish, the minimum energy state describes vortices (or antivortices) that propagate around the trap center. The number of (anti)vortices increases with the angular momentum, and they repel each other to form Abrikosov lattices. Besides vortices and antivortices there are also stagnation points where the superflow vanishes; to our knowledge the stagnation points have not been analyzed previously, in the context of the Gross-Pitaevskii equation. The Poincar\'e index formula states that the difference in the number of vortices and stagnation points can never change. When the number of stagnation points is small, they tend to aggregate into degenerate propagating structures. But when the number becomes sufficiently large, the stagnation points tend to pair up with the vortex cores, to propagate around the trap center in regular lattice arrangements. There is an analogy with the geometry of the Kosterlitz-Thouless transition, with the angular momentum of the condensate as the external control parameter instead of the temperature.


title: Substitution of the Fittest: A Novel Approach for Mitigating Disengagement in Coevolutionary Genetic Algorithms

abstract: We propose substitution of the fittest (SF), a novel technique designed to counteract the problem of disengagement in two-population competitive coevolutionary genetic algorithms. The approach presented is domain-independent and requires no calibration. In a minimal domain, we perform a controlled evaluation of the ability to maintain engagement and the capacity to discover optimal solutions. Results demonstrate that the solution discovery performance of SF is comparable with other techniques in the literature, while SF also offers benefits including a greater ability to maintain engagement and a much simpler mechanism.


title: Millikelvin scanning tunneling microscope at 20/22 T with a graphite enabled stick-slip approach and an energy resolution below 8$\mu$eV: Application to conductance quantization at 20 T in single atom point contacts of Al and Au

abstract: We describe a scanning tunneling microscope (STM) that operates at magnetic fields up to 22 Tesla and temperatures down to 80 mK. We discuss the design of the STM head, with an improved coarse approach, the vibration isolation system, and efforts to improve the energy resolution using compact filters for multiple lines. We measure the superconducting gap and Josephson effect in Aluminum and show that we can resolve features in the density of states as small as 8\,$\mu$eV. We measure the quantization of conductance in atomic size contacts and make atomic resolution and density of states images in the layered material 2H-NbSe$_2$. The latter experiments are performed by continuously operating the STM at magnetic fields of 20 T in periods of several days without interruption.


title: Low Temperature Plasma for Biology, Hygiene, and Medicine: Perspective and Roadmap

abstract: Plasma, the fourth and most pervasive state of matter in the visible universe, is a fascinating medium that is connected to the beginning of our universe itself. Man-made plasmas are at the core of many technological advances that include the fabrication of semiconductor devices, which enabled the modern computer and communication revolutions. The introduction of low temperature, atmospheric pressure plasmas to the biomedical field has ushered a new revolution in the healthcare arena that promises to introduce plasma-based therapies to combat some thorny and long-standing medical challenges. This paper presents an overview of where research is at today and discusses innovative concepts and approaches to overcome present challenges and take the field to the next level. It is written by a team of experts who took an in-depth look at the various biomedical applications, made critical analysis, and proposed ideas and concepts that should help the research community focus their efforts on clear and practical steps necessary to keep the field advancing for decades to come.


title: The stochastic X-ray variability of the accreting millisecond pulsar IGR J17062-6143

abstract: This work presents an investigation of the stochastic X-ray variability from the 164 Hz accreting millisecond pulsar IGR J17062-6143, based on regular observations collected with the Neutron Star Interior Composition Explorer between 2017 July and 2020 August. Over this period, the power density spectrum showed a stable morphology, with broad $\sim25\%$ rms band-limited noise below 16 Hz. Quasi-periodic oscillations (QPOs) were occasionally observed, with the most notably detections including a low-frequency QPO centered at 2.7 Hz and a sharp QPO centered at 115 Hz that may be a 2:3 resonance with the spin frequency. Further, the energy dependence of the band-limited noise is studied through a spectroscopic analysis of the complex covariance in two frequency intervals. It is found that the power law continuum is the primary driver for the observed variability, although the thermal (blackbody) emission also appears to be intrinsically variable in area ($5\%$ rms) and temperature ($1\%$ rms). Notably, the 1 keV emission feature seen in all X-ray spectra of IGR J17062-6143 varies with the same amplitude as the power law, but systematically lags behind that continuum emission. These results appear consistent with a scenario in which a time variable Compton-scattering corona is the primary source for the observed stochastic variability, with the variability observed in the emission feature and at the lowest photon energies being due to the disk reflection of the power law continuum.


title: Dynamical stability of giant planets: the critical adiabatic index in the presence of a solid core

abstract: The dissociation and ionization of hydrogen, during the formation of giant planets via core accretion, reduces the effective adiabatic index $\gamma$ of the gas and could trigger dynamical instability. We generalize the analysis of Chandrasekhar, who determined that the threshold for instability of a self-gravitating hydrostatic body lies at $\gamma=4/3$, to account for the presence of a planetary core, which we model as an incompressible fluid. We show that the dominant effect of the core is to stabilize the envelope to radial perturbations, in some cases completely (i.e. for all $\gamma > 1$). When instability is possible, unstable planetary configurations occupy a strip of $\gamma$ values whose upper boundary falls below $\gamma=4/3$. Fiducial evolutionary tracks of giant planets forming through core accretion appear unlikely to cross the dynamical instability strip that we define.


title: Empirically Determining Substellar Cloud Compositions in the era of JWST

abstract: Most brown dwarfs have atmospheres with temperatures cold enough to form clouds. A variety of materials likely condense, including refractory metal oxides and silicates; the precise compositions and crystal structures of predicted cloud particles depend on the modeling framework used and have not yet been empirically constrained. Spitzer has shown tentative evidence of the silicate feature in L dwarf spectra and JWST can measure these features in many L dwarfs. Here, we present new models to predict the signatures of the strongest cloud absorption features. We investigate different cloud mineral species and determine how particle size, mineralogy, and crystalline structure change spectral features. We find that silicate and refractory clouds have a strong cloud absorption feature for small particle sizes ($\leq$ 1 $\mu$m). Model spectra are compared to five brown dwarfs that show evidence of the silicate feature; models that include small particles in the upper layers of the atmosphere produce a broad cloud mineral feature, and that better match the observed spectra than the Ackerman & Marley (2001) cloud model. We simulate observations with the MIRI instrument on JWST for a range of nearby, cloudy brown dwarfs, demonstrating that these features could be readily detectable if small particles are present. Furthermore, for photometrically variable brown dwarfs, our predictions suggest that with JWST, by measuring spectroscopic variability inside and outside a mineral feature, we can establish silicate (or other) clouds as the cause of variability. Mid-infrared spectroscopy is a promising tool to empirically constrain the complex cloud condensation sequence in brown dwarf atmospheres.


title: Learning Activities in Colours and Rainbows for Programming Skill Development

abstract: We present how we have created a series of bilingual (English and Welsh) STEM activities focusing on rainbows, colours, light and optical effects. The activities were motivated by the many rainbows that appeared in windows in the UK, in support of the National Health Service at the start of the coronavirus pandemic. Rainbows are hopeful and are very fitting to be used as a positive iconic image at a time of much uncertainty. In this paper we explain how we have developed and organised the activities, focusing on colours, computer graphics and computer programming. Each lesson contains one or more activities, which enable people to take an active role in their learning. We have carefully prepared and organised several processes to guide academic colleagues to create and publish different activities in the theme. Which means that the activities appear similarly structured, can be categorised and searched in a consistent way. This structure can act as a blueprint for others to follow and apply to develop their own online course. The activities incrementally take people through learning about colour, rainbows, planning what to program, design and strategies to create colourful pictures using simple computer graphics principles based in processing.org.


title: Sedimentation effects on particle position and inertial deposition in 90{\deg} circular bends

abstract: Laminar fluid-particle flows in bend geometries are present in many industrial, pharmaceutical, and biomedical applications. Particle deposition has been studied extensively; however, little attention has been paid to the effect of particle sedimentation on particle position and deposition in different pipe geometry combinations. This study presented a comprehensive analysis of sedimentation effects on particle flow behaviour in 90{\deg} circular pipe bends of micron particles in laminar pipe flows. Pipe geometry combinations consisted of eight pipe diameters, nine bend radii, and 30 particle diameters in a range of 1 to 100 {\mu}m. The results demonstrated the locations of particles that sedimented to the bottom half of the straight pipe section, and the particle positions upstream from the pipe bend entrance, which was no longer in a fully developed profile. These new locations represent the effects of gravity, pulling the particles down. While obtaining these positions can be found through CFD analysis, we proposed an analytical solution to predict the particle trajectory from different release locations that would help to identify the initial particle distribution at locations upstream to the bend, to obviate the need for long upstream straight pipe sections in the CFD analysis.


title: RadioMic: Sound Sensing via mmWave Signals

abstract: Voice interfaces has become an integral part of our lives, with the proliferation of smart devices. Today, IoT devices mainly rely on microphones to sense sound. Microphones, however, have fundamental limitations, such as weak source separation, limited range in the presence of acoustic insulation, and being prone to multiple side-channel attacks. In this paper, we propose RadioMic, a radio-based sound sensing system to mitigate these issues and enrich sound applications. RadioMic constructs sound based on tiny vibrations on active sources (e.g., a speaker or human throat) or object surfaces (e.g., paper bag), and can work through walls, even a soundproof one. To convert the extremely weak sound vibration in the radio signals into sound signals, RadioMic introduces radio acoustics, and presents training-free approaches for robust sound detection and high-fidelity sound recovery. It then exploits a neural network to further enhance the recovered sound by expanding the recoverable frequencies and reducing the noises. RadioMic translates massive online audios to synthesized data to train the network, and thus minimizes the need of RF data. We thoroughly evaluate RadioMic under different scenarios using a commodity mmWave radar. The results show RadioMic outperforms the state-of-the-art systems significantly. We believe RadioMic provides new horizons for sound sensing and inspires attractive sensing capabilities of mmWave sensing devices


title: Well-posedness and optimal control for a Cahn-Hilliard-Oono system with control in the mass term

abstract: The paper treats the problem of optimal distributed control of a Cahn-Hilliard-Oono system in $\mathbb{R}^d$, $1\leq d\leq 3$, with the control located in the mass term and admitting general potentials that include both the case of a regular potential and the case of some singular potential. The first part of the paper is concerned with the dependence of the phase variable on the control variable. In the case of a logarithmic potential, we need to prove an ad hoc strict separation property, and for this reason we have to restrict ourselves to the case $d = 2$. In the rest of the work, we study the necessary first-order optimality conditions, which are proved under suitable compatibility conditions on the initial datum of the phase variable and the time derivative of the control, at least in case of potentials having unbounded domain. PLEASE NOTE: A revised version of this paper has been published in Discrete Contin. Dyn. Syst. Ser. S 15 (2022), 2135-2172: the authors point out that both the published version and the arXiv preprint contain a trivial mistake in formula (2.21), which leads to suitable changes in the assumption (2.10) and in some occurrences of it. Namely, (2.21) and (2.10) become $ \bar \varphi(t) = \bar \varphi_0 e^{-t} + \int_0^t e^{-(t-s)} \bar u(s) \, ds$ for every $t\in[0,T]$, and $-M - (\bar\varphi_0)^- , M + (\bar\varphi_0)^+$ belong to the interior of $D(\beta)$, respectively. The formula just after (2.21) becomes $-M - (\bar\varphi_0)^- \leq \bar \varphi(t) \leq M + (\bar\varphi_0)^+$ and a similar change must be done in (4.20). The first line of (2.24) is now implied by the new (2.21). Next, (2.39) becomes $-(M+ R) - (\bar\varphi_0)^- , (M+ R) + (\bar\varphi_0)^+$ belong to the interior of $D(\beta)$. Similar changes must be done in the interval written just two lines before (4.18) and in (4.18) itself.


title: Feature Augmented Hybrid CNN for Stress Recognition Using Wrist-based Photoplethysmography Sensor

abstract: Stress is a physiological state that hampers mental health and has serious consequences to physical health. Moreover, the COVID-19 pandemic has increased stress levels among people across the globe. Therefore, continuous monitoring and detection of stress are necessary. The recent advances in wearable devices have allowed the monitoring of several physiological signals related to stress. Among them, wrist-worn wearable devices like smartwatches are most popular due to their convenient usage. And the photoplethysmography (PPG) sensor is the most prevalent sensor in almost all consumer-grade wrist-worn smartwatches. Therefore, this paper focuses on using a wrist-based PPG sensor that collects Blood Volume Pulse (BVP) signals to detect stress which may be applicable for consumer-grade wristwatches. Moreover, state-of-the-art works have used either classical machine learning algorithms to detect stress using hand-crafted features or have used deep learning algorithms like Convolutional Neural Network (CNN) which automatically extracts features. This paper proposes a novel hybrid CNN (H-CNN) classifier that uses both the hand-crafted features and the automatically extracted features by CNN to detect stress using the BVP signal. Evaluation on the benchmark WESAD dataset shows that, for 3-class classification (Baseline vs. Stress vs. Amusement), our proposed H-CNN outperforms traditional classifiers and normal CNN by 5% and 7% accuracy, and 10% and 7% macro F1 score, respectively. Also for 2-class classification (Stress vs. Non-stress), our proposed H-CNN outperforms traditional classifiers and normal CNN by 3% and ~5% accuracy, and ~3% and ~7% macro F1 score, respectively.


title: Generalized Tensor Summation Compressive Sensing Network (GTSNET): An Easy to Learn Compressive Sensing Operation

abstract: In CS literature, the efforts can be divided into two groups: finding a measurement matrix that preserves the compressed information at the maximum level, and finding a reconstruction algorithm for the compressed information. In the traditional CS setup, the measurement matrices are selected as random matrices, and optimization-based iterative solutions are used to recover the signals. However, when we handle large signals, using random matrices become cumbersome especially when it comes to iterative optimization-based solutions. Even though recent deep learning-based solutions boost the reconstruction accuracy performance while speeding up the recovery, still jointly learning the whole measurement matrix is a difficult process. In this work, we introduce a separable multi-linear learning of the CS matrix by representing it as the summation of arbitrary number of tensors. For a special case where the CS operation is set as a single tensor multiplication, the model is reduced to the learning-based separable CS; while a dense CS matrix can be approximated and learned as the summation of multiple tensors. Both cases can be used in CS of two or multi-dimensional signals e.g., images, multi-spectral images, videos, etc. Structural CS matrices can also be easily approximated and learned in our multi-linear separable learning setup with structural tensor sum representation. Hence, our learnable generalized tensor summation CS operation encapsulates most CS setups including separable CS, non-separable CS (traditional vector-matrix multiplication), structural CS, and CS of the multi-dimensional signals. For both gray-scale and RGB images, the proposed scheme surpasses most state-of-the-art solutions, especially in lower measurement rates. Although the performance gain remains limited from tensor to the sum of tensor representation for gray-scale images, it becomes significant in the RGB case.


title: Pattern Recognition in Vital Signs Using Spectrograms

abstract: Spectrograms visualize the frequency components of a given signal which may be an audio signal or even a time-series signal. Audio signals have higher sampling rate and high variability of frequency with time. Spectrograms can capture such variations well. But, vital signs which are time-series signals have less sampling frequency and low-frequency variability due to which, spectrograms fail to express variations and patterns. In this paper, we propose a novel solution to introduce frequency variability using frequency modulation on vital signs. Then we apply spectrograms on frequency modulated signals to capture the patterns. The proposed approach has been evaluated on 4 different medical datasets across both prediction and classification tasks. Significant results are found showing the efficacy of the approach for vital sign signals. The results from the proposed approach are promising with an accuracy of 91.55% and 91.67% in prediction and classification tasks respectively.


title: Responding to Illegal Activities Along the Canadian Coastlines Using Reinforcement Learning

abstract: This article elaborates on how machine learning (ML) can leverage the solution of a contemporary problem related to the security of maritime domains. The worldwide ``Illegal, Unreported, and Unregulated'' (IUU) fishing incidents have led to serious environmental and economic consequences which involve drastic changes in our ecosystems in addition to financial losses caused by the depletion of natural resources. The Fisheries and Aquatic Department (FAD) of the United Nation's Food and Agriculture Organization (FAO) issued a report which indicated that the annual losses due to IUU fishing reached $25 Billion. This imposes negative impacts on the future-biodiversity of the marine ecosystem and domestic Gross National Product (GNP). Hence, robust interception mechanisms are increasingly needed for detecting and pursuing the unrelenting illegal fishing incidents in maritime territories. This article addresses the problem of coordinating the motion of a fleet of marine vessels (pursuers) to catch an IUU vessel while still in local waters. The problem is formulated as a pursuer-evader problem that is tackled within an ML framework. One or more pursuers, such as law enforcement vessels, intercept an evader (i.e., the illegal fishing ship) using an online reinforcement learning mechanism that is based on a value iteration process. It employs real-time navigation measurements of the evader ship as well as those of the pursuing vessels and returns back model-free interception strategies.


title: Respiratory Rate Estimation Based on WiFi Frame Capture

abstract: This paper presents a method that estimates the respiratory rate based on the frame capturing of wireless local area networks. The method uses beamforming feedback matrices (BFMs) contained in the captured frames, which is a rotation matrix of channel state information (CSI). BFMs are transmitted unencrypted and easily obtained using frame capturing, requiring no specific firmware or WiFi chipsets, unlike the methods that use CSI. Such properties of BFMs allow us to apply frame capturing to various sensing tasks, e.g., vital sensing. In the proposed method, principal component analysis is applied to BFMs to isolate the effect of the chest movement of the subject, and then, discrete Fourier transform is performed to extract respiratory rates in a frequency domain. Experimental evaluation results confirm that the frame-capture-based respiratory rate estimation can achieve estimation error lower than 3.2 breaths/minute.


title: Quantum Meets the Minimum Circuit Size Problem

abstract: In this work, we initiate the study of the Minimum Circuit Size Problem (MCSP) in the quantum setting. MCSP is a problem to compute the circuit complexity of Boolean functions. It is a fascinating problem in complexity theory -- its hardness is mysterious, and a better understanding of its hardness can have surprising implications to many fields in computer science. We first define and investigate the basic complexity-theoretic properties of minimum quantum circuit size problems for three natural objects: Boolean functions, unitaries, and quantum states. We show that these problems are not trivially in NP but in QCMA (or have QCMA protocols). Next, we explore the relations between the three quantum MCSPs and their variants. We discover that some reductions that are not known for classical MCSP exist for quantum MCSPs for unitaries and states, e.g., search-to-decision reduction and self-reduction. Finally, we systematically generalize results known for classical MCSP to the quantum setting (including quantum cryptography, quantum learning theory, quantum circuit lower bounds, and quantum fine-grained complexity) and also find new connections to tomography and quantum gravity. Due to the fundamental differences between classical and quantum circuits, most of our results require extra care and reveal properties and phenomena unique to the quantum setting. Our findings could be of interest for future studies, and we post several open problems for further exploration along this direction.


title: A General Regularized Distributed Solution for System State Estimation from Relative Measurements

abstract: This work presents a novel general regularized distributed solution for the state estimation problem in networked systems. Resting on the graph-based representation of sensor networks and adopting a multivariate least-squares approach, the designed solution exploits the set of the available inter-sensor relative measurements and leverages a general regularization framework, whose parameter selection is shown to control the estimation procedure convergence performance. As confirmed by the numerical results, this new estimation scheme allows (i) the extension of other approaches investigated in the literature and (ii) the convergence optimization in correspondence to any (undirected) graph modeling the given sensor network.


title: Incremental learning of LSTM framework for sensor fusion in attitude estimation

abstract: This paper presents a novel method for attitude estimation of an object in 3D space by incremental learning of the Long-Short Term Memory (LSTM) network. Gyroscope, accelerometer, and magnetometer are few widely used sensors in attitude estimation applications. Traditionally, multi-sensor fusion methods such as the Extended Kalman Filter and Complementary Filter are employed to fuse the measurements from these sensors. However, these methods exhibit limitations in accounting for the uncertainty, unpredictability, and dynamic nature of the motion in real-world situations. In this paper, the inertial sensors data are fed to the LSTM network which are then updated incrementally to incorporate the dynamic changes in motion occurring in the run time. The robustness and efficiency of the proposed framework is demonstrated on the dataset collected from a commercially available inertial measurement unit. The proposed framework offers a significant improvement in the results compared to the traditional method, even in the case of a highly dynamic environment. The LSTM framework-based attitude estimation approach can be deployed on a standard AI-supported processing module for real-time applications.


title: Chemistry of a light impurity in a Bose-Einstein condensate

abstract: In ultracold atomic gases, a unique interplay arises between phenomena known from condensed matter physics, few-body physics and chemistry. Similar to an electron in a solid, an impurity in an ultracold gas can get dressed by excitations from the medium, forming a quasiparticle called the polaron. We study how dressing of an impurity leads to a modification of its chemical reactivity. Using a Gaussian state variational method in the frame of the impurity, we demonstrate that three-body correlations lead to an instability of the polaron. This instability is connected to an Efimov resonance, but shifted to smaller interactions by many-body effects, showing that polaron formation stimulates Efimov physics and the associated chemistry.


title: From Efimov Physics to the Bose Polaron using Gaussian States

abstract: Since the Efimov effect was introduced in 1970, a detailed theoretical understanding of Efimov physics has been developed in the few-body context. However, it has proven to be challenging to describe the role Efimov-type correlations play in many-body systems such as quenched or collapsing Bose-Einstein condensates (BECs). To study the impact the Efimov effect can have in such scenarios, we consider a light impurity immersed in a weakly interacting BEC, forming a Bose polaron. In this case, the higher-order correlations are localized around the impurity, making it more feasible to develop a theoretical description. Specifically, we employ a Gaussian state variational Ansatz in the reference frame of the impurity, capable of both capturing the Efimov effect and the formation of the polaron cloud. We find that the Efimov effect leads to a cooperative binding of bosons to the impurity and the formation of a many-body bound state. As a result, the polaron is not the ground state, but rendered a metastable excited state which can decay into these Efimov clusters. While this decay is slow for small interaction strengths, it becomes more prominent as the attractive scattering length increases, up to the point where the polaron becomes completely unstable. This critical scattering length can be interpreted as a many-body shifted Efimov resonance, where the scattering of two excitations of the bath with the polaron can lead to bound state formation. Compared to the few-body case, the resonance is shifted to smaller attractive scattering lengths due to the participation of the polaron cloud in the cooperative binding process. This corresponds to an intriguing scenario of polaron-assisted chemistry, where many-body effects lead to enhanced signal of the chemical recombination process, which can be directly probed in state-of-the-art experiments.


title: Dynamic Control for Random Access in Deadline-Constrained Broadcasting

abstract: This paper considers random access in deadline-constrained broadcasting with frame-synchronized traffic. To enhance the maximum achievable timely delivery ratio (TDR), we define a dynamic control scheme that allows each active node to determine the transmission probability with certainty based on the current delivery urgency and the knowledge of current contention intensity. For an idealized environment where the contention intensity is completely known, we develop an analytical framework based on the theory of Markov Decision Process (MDP), which leads to an optimal scheme by applying backward induction. For a realistic environment where the contention intensity is incompletely known, we develop a framework using Partially Observable Markov Decision Process (POMDP), which can in theory be solved. We show that for both environments, there exists an optimal scheme that is optimal over all types of policies. To overcome the infeasibility in obtaining an optimal or near-optimal scheme from the POMDP framework, we investigate the behaviors of the optimal scheme for two extreme cases in the MDP framework, and leverage intuition gained from these behaviors to propose a heuristic scheme for the realistic environment with TDR close to the maximum achievable TDR in the idealized environment. In addition, we propose an approximation on the knowledge of contention intensity to further simplify this heuristic scheme. Numerical results with respect to a wide range of configurations are provided to validate our study.


title: Shift-invariant waveform learning on epileptic ECoG

abstract: Seizure detection algorithms must discriminate abnormal neuronal activity associated with a seizure from normal neural activity in a variety of conditions. Our approach is to seek spatiotemporal waveforms with distinct morphology in electrocorticographic (ECoG) recordings of epileptic patients that are indicative of a subsequent seizure (preictal) versus non-seizure segments (interictal). To find these waveforms we apply a shift-invariant k-means algorithm to segments of spatially filtered signals to learn codebooks of prototypical waveforms. The frequency of the cluster labels from the codebooks is then used to train a binary classifier that predicts the class (preictal or interictal) of a test ECoG segment. We use the Matthews correlation coefficient to evaluate the performance of the classifier and the quality of the codebooks. We found that our method finds recurrent non-sinusoidal waveforms that could be used to build interpretable features for seizure prediction and that are also physiologically meaningful.


title: Transformation-Enabled Precondition Inference

abstract: Precondition inference is a non-trivial problem with important applications in program analysis and verification. We present a novel iterative method for automatically deriving preconditions for the safety and unsafety of programs. Each iteration maintains over-approximations of the set of safe and unsafe initial states; which are used to partition the program's initial states into those known to be safe, known to be unsafe and unknown. We then construct revised programs with those unknown initial states and iterate the procedure until the approximations are disjoint or some termination criteria are met. An experimental evaluation of the method on a set of software verification benchmarks shows that it can infer precise preconditions (sometimes optimal) that are not possible using previous methods.


title: Analysis of nonconforming IFE methods and a new scheme for elliptic interface problems

abstract: In this paper, an important discovery has been found for nonconforming immersed finite element (IFE) methods using the integral values on edges as degrees of freedom for solving elliptic interface problems. We show that those IFE methods without penalties are not guaranteed to converge optimally if the tangential derivative of the exact solution and the jump of the coefficient are not zero on the interface. A nontrivial counter example is also provided to support our theoretical analysis. To recover the optimal convergence rates, we develop a new nonconforming IFE method with additional terms locally on interface edges. The new method is parameter-free which removes the limitation of the conventional partially penalized IFE method. We show the IFE basis functions are unisolvent on arbitrary triangles which is not considered in the literature. Furthermore, different from multipoint Taylor expansions, we derive the optimal approximation capabilities of both the Crouzeix-Raviart and the rotated-$Q_1$ IFE spaces via a unified approach which can handle the case of variable coefficients easily. Finally, optimal error estimates in both $H^1$- and $L^2$- norms are proved and confirmed with numerical experiments.


title: Dynamic Semantic Occupancy Mapping using 3D Scene Flow and Closed-Form Bayesian Inference

abstract: This paper reports on a dynamic semantic mapping framework that incorporates 3D scene flow measurements into a closed-form Bayesian inference model. Existence of dynamic objects in the environment can cause artifacts and traces in current mapping algorithms, leading to an inconsistent map posterior. We leverage state-of-the-art semantic segmentation and 3D flow estimation using deep learning to provide measurements for map inference. We develop a Bayesian model that propagates the scene with flow and infers a 3D continuous (i.e., can be queried at arbitrary resolution) semantic occupancy map outperforming its static counterpart. Extensive experiments using publicly available data sets show that the proposed framework improves over its predecessors and input measurements from deep neural networks consistently.


title: Localization, $\mathcal{PT}$-Symmetry Breaking and Topological Transitions in non-Hermitian Quasicrystals

abstract: According to the topological band theory of a Hermitian system, the different electronic phases are classified in terms of topological invariants, wherein the transition between the two phases characterized by a different topological invariant is the primary signature of a topological phase transition. Recently, it has been argued that the delocalization-localization transition in a quasicrystal, described by the non-Hermitian $\mathcal{PT}$-symmetric extension of the Aubry-Andr\'{e}-Harper (AAH) Hamiltonian can also be identified as a topological phase transition. Interestingly, the $\mathcal{PT}$-symmetry also breaks down at the same critical point. However, in this article, we have shown that the delocalization-localization transition and the $\mathcal{PT}$-symmetry breaking are not connected to a topological phase transition. To demonstrate this, we have studied the non-Hermitian $\mathcal{PT}$-symmetric AAH Hamiltonian in the presence of Rashba Spin-Orbit (RSO) coupling. We have obtained an analytical expression of the topological transition point and compared it with the numerically obtained critical points. We have found that, except in some special cases, the critical point and the topological transition point are not the same. In fact, the delocalization-localization transition takes place earlier than the topological transition whenever they do not coincide.


title: Single camera high repetition rate two-color formaldehyde planar laser-induced fluorescence thermometry with a wavelength-switching burst mode laser

abstract: We consider a two-color formaldehyde PLIF thermometry scheme using a wavelength-switching injection seeding Nd:YAG laser at 355 nm. The 28183.5 cm-1 and 28184.5 cm-1 peaks of formaldehyde are used to measure low temperature combustion zone. Using a burst mode amplifier and a high speed camera, high-repetition rate (20 kHz) temperature field measurement is validated on a laminar coflow diffusion flame and demonstrated on a turbulent confined jet in hot crossflow flame.


title: Mechanical Properties of Wilberforce Pendulum

abstract: This paper shows the study of interesting mechanical properties of Wilberforce pendulum. Analyzing qualitatively of the pendulum, it is able to know how the phenomenon occurs. By setting of the quantitative model, equation of the motion is derived. Considering the mass and moment of inertia of the spring, the experiment was done by changing the moment of inertia as the main parameter. The results were analyzed by defining oscillation ratio and conversion factor.


title: Helical spin order in EuFe2As2 and EuRbFe4As4 single crystals

abstract: In this work, fabrication and characterization of magnetic properties of EuFe2As2 and EuRbFe4As4 single crystals is reported. Magnetization measurements of samples with well defined thin film geometry and crystal orientation demonstrate a striking similarity in ferromagnetic properties of Eu subsystems in these two compounds. Measurements with magnetic field applied along ab crystal planes reveal meta-magnetic transition in both compounds. Numerical studies employing the Jz1-Jz2 Heisenberg model suggest that the ground state of the magnetic order in Eu subsystem for both compounds is the helical spin order with the helical angle about 2pi/5, while the meta-magnetic transition is the helix-to-fan first order phase transition.


title: Genetic optimization of quantum annealing

abstract: The study of optimal control of quantum annealing by modulating the pace of evolution and by introducing a counterdiabatic potential has gained significant attention in recent times. In this work, we present a numerical approach based on genetic algorithms to improve the performance of quantum annealing, which evades the Landau-Zener transitions to navigate to the ground state of the final Hamiltonian with high probability. We optimize the annealing schedules starting from polynomial ansatz by treating their coefficients as chromosomes of the genetic algorithm. We also explore shortcuts to adiabaticity by computing a practically feasible $k$-local optimal driving operator, showing that even for $k=1$ we achieve substantial improvement of the fidelity over the standard annealing solution. With these genetically optimized annealing schedules and/or optimal driving operators, we are able to perform quantum annealing in relatively short time-scales and with larger fidelity compared to traditional approaches.


title: Nodal quintic del Pezzo threefolds and their derived categories

abstract: We construct a Kawamata type semiorthogondal decomposition for the bounded derived category of coherent sheaves of nodal quintic del Pezzo threefolds, decomposing the bounded derived category into bounded derived categories of finite dimensional algebras. This is achieved by constructing birational maps from nodal quintic del Pezzo threefolds to quadric surface fibrations over the projective line.


title: Here and There with Arithmetic

abstract: In the theory of answer set programming, two groups of rules are called strongly equivalent if, informally speaking, they have the same meaning in any context. The relationship between strong equivalence and the propositional logic of here-and-there allows us to establish strong equivalence by deriving rules of each group from rules of the other. In the process, rules are rewritten as propositional formulas. We extend this method of proving strong equivalence to an answer set programming language that includes operations on integers. The formula representing a rule in this language is a first-order formula that may contain comparison symbols among its predicate constants, and symbols for arithmetic operations among its function constants. The paper is under consideration for acceptance in TPLP.


title: Plethysms of Chromatic and Tutte Symmetric Functions

abstract: Plethysm is a fundamental operation in symmetric function theory, derived directly from its connection with representation theory. However, it does not admit a simple combinatorial interpretation, and finding coefficients of Schur function plethysms is a major open question. In this paper, we introduce a graph-theoretic interpretation for any plethysm based on the chromatic symmetric function. We use this interpretation to give simple proofs of new and previously known plethystic identities, as well as chromatic symmetric function identities.


title: On-Shell Recursion Relations for Nonrelativistic Effective Field Theories

abstract: We derive on-shell recursion relations for nonrelativistic effective field theories (EFTs) with enhanced soft limits. The recursion relations are illustrated through analytic calculation of tree-level scattering amplitudes in theories with a complex Schr\"odinger-type field, real scalar with linear dispersion relation, and real scalar with Lifshitz-type dispersion relation. Our results show that the landscape of gapless nonrelativistic EFTs with local $S$-matrix can be constrained by soft theorems and the consistency of the low-energy $S$-matrix similarly to massless relativistic EFTs.


title: Quantum Quantile Mechanics: Solving Stochastic Differential Equations for Generating Time-Series

abstract: We propose a quantum algorithm for sampling from a solution of stochastic differential equations (SDEs). Using differentiable quantum circuits (DQCs) with a feature map encoding of latent variables, we represent the quantile function for an underlying probability distribution and extract samples as DQC expectation values. Using quantile mechanics we propagate the system in time, thereby allowing for time-series generation. We test the method by simulating the Ornstein-Uhlenbeck process and sampling at times different from the initial point, as required in financial analysis and dataset augmentation. Additionally, we analyse continuous quantum generative adversarial networks (qGANs), and show that they represent quantile functions with a modified (reordered) shape that impedes their efficient time-propagation. Our results shed light on the connection between quantum quantile mechanics (QQM) and qGANs for SDE-based distributions, and point the importance of differential constraints for model training, analogously with the recent success of physics informed neural networks.


title: Probabilistic motion planning for non-Euclidean and multi-vehicle problems

abstract: Trajectory planning tasks for non-holonomic or collaborative systems are naturally modeled by state spaces with non-Euclidean metrics. However, existing proofs of convergence for sample-based motion planners only consider the setting of Euclidean state spaces. We resolve this issue by formulating a flexible framework and set of assumptions for which the widely-used PRM*, RRT, and RRT* algorithms remain asymptotically optimal in the non-Euclidean setting. The framework is compatible with collaborative trajectory planning: given a fleet of robotic systems that individually satisfy our assumptions, we show that the corresponding collaborative system again satisfies the assumptions and therefore has guaranteed convergence for the trajectory-finding methods. Our joint state space construction builds in a coupling parameter $1\leq p\leq \infty$, which interpolates between a preference for minimizing total energy at one extreme and a preference for minimizing the travel time at the opposite extreme. We illustrate our theory with trajectory planning for simple coupled systems, fleets of Reeds-Shepp vehicles, and a highly non-Euclidean fractal space.


title: Higher-order motif analysis in hypergraphs

abstract: A deluge of new data on social, technological and biological networked systems suggests that a large number of interactions among system units are not limited to pairs, but rather involve a higher number of nodes. To properly encode such higher-order interactions, richer mathematical frameworks such as hypergraphs are needed, where hyperlinks describe connections among an arbitrary number of nodes. Here we introduce the concept of higher-order motifs, small connected subgraphs where vertices may be linked by interactions of any order. We provide lower and upper bounds on the number of higher-order motifs as a function of the motif size, and propose an efficient algorithm to extract complete higher-order motif profiles from empirical data. We identify different families of hypergraphs, characterized by distinct higher-order connectivity patterns at the local scale. We also capture evidences of structural reinforcement, a mechanism that associates higher strengths of higher-order interactions for the nodes that interact more at the pairwise level. Our work highlights the informative power of higher-order motifs, providing a first way to extract higher-order fingerprints in hypergraphs at the network microscale.


title: On the Exponential Sample Complexity of the Quantum State Sign Estimation Problem

abstract: We demonstrate that the ability to estimate the relative sign of an arbitrary $n$-qubit quantum state (with real amplitudes), given only $k$ copies of that state, would yield a $kn$-query algorithm for unstructured search. Thus the quantum sample complexity of sign estimation must be exponential: $\Omega(2^{n/2}/n)$. In particular, we show that an efficient procedure for solving the sign estimation problem would allow for a polynomial time solution to the NP-complete problem 3-SAT.


title: Almost everywhere non-uniqueness of integral curves for divergence-free Sobolev vector fields

abstract: We construct divergence-free Sobolev vector fields in C([0,1];W^{1,r}(T^d;R^d)) with r < d and d >=2 which simultaneously admit any finite number of distinct positive solutions to the continuity equation. We then show that the vector fields we produce have at least as many integral curves starting from a.e. point of T^d as the number of distinct positive solutions to the continuity equation these vector fields admit. Our work uses convex integration techniques for the continuity equation introduced in [L. Szekelyhidi and S. Modena, Annals of PDE, 2018] and [E. Brue, M. Colombo, C. De Lellis, Archive for Rational Mechanics and Analysis, 2021 ] to study non-uniqueness for positive solutions. We then infer non-uniqueness for integral curves from Ambrosio superposition principle.


title: A characterization of finite factorization positive monoids

abstract: We provide a characterization of the positive monoids (i.e., additive submonoids of the nonnegative real numbers) that satisfy the finite factorization property. As a result, we establish that positive monoids with well-ordered generating sets satisfy the finite factorization property, while positive monoids with co-well-ordered generating sets satisfy this property if and only if they satisfy the bounded factorization property.


title: Time-reversal-broken Weyl semimetal in the Hofstadter regime

abstract: We study the phase diagram for a lattice model of a time-reversal-broken three-dimensional Weyl semimetal (WSM) in an orbital magnetic field $B$ with a flux of $p/q$ per unit cell ($0\le p \le q-1$), with minimal crystalline symmetry. We find several interesting phases: (i) WSM phases with $2q$, $4q$, $6q$, and $8q$ Weyl nodes and corresponding surface Fermi arcs, (ii) a layered Chern insulating (LCI) phase, gapped in the bulk, but with gapless surface states, (iii) a phase in which some bulk bands are gapless with Weyl nodes, coexisting with others that are gapped but topologically nontrivial, adiabatically connected to an LCI phase, (iv) a new gapped trivially insulating phase (I$'$) with (non-topological) counter-propagating surface states, which could be gapped out in the absence of crystal symmetries. Importantly, we are able to obtain the phase boundaries analytically for all $p,q$. Analyzing the gaps for $p=1$ and very large $q$ enables us to smoothly take the zero-field limit, even though the phase diagrams look ostensibly very different for $q=1, B=0$, and $q\to\infty, B\to 0$.


title: The Einstein-Hilbert-Palatini formalism in Pseudo-Finsler Geometry

abstract: A systematic development of the so-called Palatini formalism is carried out for pseudo-Finsler metrics $L$ of any signature. Substituting in the classical Einstein-Hilbert-Palatini functional the scalar curvature by the Finslerian Ricci scalar constructed with an independent nonlinear connection $\mathrm{N}$, the affine and metric equations for $(\mathrm{N},L)$ are obtained. In Lorentzian signature with vanishing mean Landsberg tensor Lan$_i$, both the Finslerian Hilbert metric equation and the classical Palatini conclusions are recovered by means of a combination of techniques involving the (Riemannian) maximum principle and an original argument about divisibility and fiberwise analyticity. Some of these findings are also extended to (positive definite) Riemannian metrics by using the eigenvalues of the Laplacian. When Lan$_i\neq 0$, the Palatini conclusions fail necessarily, however, a good number of properties of the solutions remain. The framework and proofs are built up in detail.


title: Integer circulant determinants of order 15

abstract: We consider the values taken by $n\times n$ circulant determinants with integer entries when $n$ is the product of two distinct odd primes $p,q$. These correspond to the integer group determinants for $\mathbb Z_{pq}$, the cyclic group of order $pq$. We show that $p^2$ and $q^2$ are not determinants (more generally we show that the classic necessary divisibility conditions are never sufficient when $n$ contains at least two odd primes). We obtain a complete description of the integer group determinants for $\mathbb Z_{15}$ (the smallest unresolved group) and partial results for general $n=3p.$


title: A sensitivity study of VBS and diboson WW to dimension-6 EFT operators at the LHC

abstract: We present a parton-level study of electro-weak production of vector-boson pairs at the Large Hadron Collider, establishing the sensitivity to a set of dimension-six operators in the Standard Model Effective Field Theory (SMEFT). Different final states are statistically combined, and we discuss how the orthogonality and interdependence of different analyses must be considered to obtain the most stringent constraints. The main novelties of our study are the inclusion of SMEFT effects in non-resonant diagrams and in irreducible QCD backgrounds, and an exhaustive template analysis of optimal observables for each operator and process considered. We also assess for the first time the sensitivity of vector-boson-scattering searches in semileptonic final states.


title: Cycles of even-odd drop permutations and continued fractions of Genocchi numbers

abstract: Recently, Lazar and Wachs (arXiv:1910.07651) showed that the (median) Genocchi numbers play a fundamental role in the study of the homogenized Linial arrangement and obtained two new permutation models (called D-permutations and E-permutations) for (median) Genocchi numbers. They further conjecture that the distributions of cycle numbers over the two models are equal. In a follow-up, Eu et al. (arXiv:2103.09130) further proved the gamma-positivity of the descent polynomials of even-odd descent permutations, which are in bijection with E-permutations by Foata's fundamental transformation. This paper merges the above two papers by considering a general moment sequence which encompasses the number of cycles and number of drops of E-permutations. Using the combinatorial theory of continued fraction, the moment connection enables us to confirm Lazar-Wachs' conjecture and obtain a natural $(p,q)$-analogue of Eu et al's descent polynomials. Furthermore, we show that the $\gamma$-coefficients of our $(p,q)$-analogue of descent polynomials have the same factorization flavor as the $\gamma$-coeffcients of Br\"and\'en's $(p,q)$-Eulerian polynomials.


title: A novel approach to asteroid impact monitoring and hazard assessment

abstract: Orbit-determination programs find the orbit solution that best fits a set of observations by minimizing the RMS of the residuals of the fit. For near-Earth asteroids, the uncertainty of the orbit solution may be compatible with trajectories that impact Earth. This paper shows how incorporating the impact condition as an observation in the orbit-determination process results in a robust technique for finding the regions in parameter space leading to impacts. The impact pseudo-observation residuals are the b-plane coordinates at the time of close approach and the uncertainty is set to a fraction of the Earth radius. The extended orbit-determination filter converges naturally to an impacting solution if allowed by the observations. The uncertainty of the resulting orbit provides an excellent geometric representation of the virtual impactor. As a result, the impact probability can be efficiently estimated by exploring this region in parameter space using importance sampling. The proposed technique can systematically handle a large number of estimated parameters, account for nongravitational forces, deal with nonlinearities, and correct for non-Gaussian initial uncertainty distributions. The algorithm has been implemented into a new impact monitoring system at JPL called Sentry-II, which is undergoing extensive testing. The main advantages of Sentry-II over JPL's currently operating impact monitoring system Sentry are that Sentry-II can systematically process orbits perturbed by nongravitational forces and that it is generally more robust when dealing with pathological cases. The runtimes and completeness of both systems are comparable, with the impact probability of Sentry-II for 99% completeness being $3\times10^{-7}$.


title: Beam-Slicing for Jammer Mitigation in mmWave Massive MU-MIMO

abstract: Millimeter-wave (mmWave) massive multi-user multiple-input multiple-output (MU-MIMO) technology promises unprecedentedly high data rates for next-generation wireless systems. To be practically viable, mmWave massive MU-MIMO basestations (BS) must (i) rely on low-resolution data-conversion and (ii) be robust to jammer interference. This paper considers the problem of mitigating the impact of a permanently transmitting jammer during uplink transmission to a BS equipped with low-resolution analog-to-digital converters (ADCs). To this end, we propose SNIPS, short for Soft-Nulling of Interferers with Partitions in Space. SNIPS combines beam-slicing---a localized, analog spatial transform that focuses the jammer energy onto a subset of all ADCs---together with a soft-nulling data detector that exploits knowledge of which ADCs are contaminated by jammer interference. Our numerical results show that SNIPS is able to successfully serve 65% of the user equipments (UEs) for scenarios in which a conventional antenna-domain soft-nulling data detector is only able to serve 2% of the UEs.


title: Adaptive Simulated Annealing with Greedy Search for the Circle Bin Packing Problem

abstract: We introduce a new bin packing problem, termed the circle bin packing problem with circular items (CBPP-CI). The problem involves packing all the circular items into multiple identical circle bins as compact as possible with the objective of minimizing the number of used bins. We first define the tangent occupying action (TOA) and propose a constructive greedy algorithm that sequentially packs the items into places tangent to the packed items or the bin boundaries. Moreover, to avoid falling into a local minimum trap and efficiently judge whether an optimal solution has been established, we continue to present the adaptive simulated annealing with greedy search (ASA-GS) algorithm that explores and exploits the search space efficiently. Specifically, we offer two novel local perturbation strategies to jump out of the local optimum and incorporate the greedy search to achieve faster convergence. The parameters of ASA-GS are adaptive according to the number of items so that they can be size-agnostic across the problem scale. We design two sets of new benchmark instances, and the empirical results show that ASA-GS completely outperforms the constructive greedy algorithm. Moreover, the packing density of ASA-GS on the top few dense bins is much higher than that of the state-of-the-art algorithm for the single circle packing problem, inferring the high quality of the packing solutions for CBPP-CI.


title: Cosmological parameters from the likelihood analysis of the galaxy power spectrum and bispectrum in real space

abstract: We present a joint likelihood analysis of the halo power spectrum and bispectrum in real space. We take advantage of a large set of numerical simulations and of an even larger set of halo mock catalogs to provide a robust estimate of the covariance properties. We derive constraints on bias and cosmological parameters assuming a theoretical model from perturbation theory at one-loop for the power spectrum and tree-level for the bispectrum. By means of the Deviance Information Criterion, we select a reference bias model dependent on seven parameters that can describe the data up to $k_{\rm max,P}=0.3\, h \, {\rm Mpc}^{-1}$ for the power spectrum and $k_{\rm max,B}=0.09\, h \, {\rm Mpc}^{-1}$ for the bispectrum at redshift $z=1$. This model is able to accurately recover three selected cosmological parameters even for the rather extreme total simulation volume of $1000\, h^{-3} \, {\rm Gpc}^3$. With the same tools, we study how relations among bias parameters can improve the fit while reducing the parameter space. In addition, we compare common approximations to the covariance matrix against the full covariance estimated from the mocks, and quantify the (non-negligible) effect of ignoring the cross-covariance between the two statistics. Finally, we explore different selection criteria for the triangular configurations to include in the analysis, showing that excluding nearly equilateral triangles rather than simply imposing a fixed maximum $k_{\rm max,B}$ on all triangle sides can lead to a better exploitation of the information contained in the bispectrum.


title: A Generalized Space-Frequency Index Modulation Scheme for Downlink MIMO Transmissions with Improved Diversity

abstract: Multidimensional Index Modulations (IM) are a novel alternative to conventional modulations which can bring considerable benefits for future wireless networks. Within this scope, in this paper we present a new scheme, named as Precoding-aided Transmitter side Generalized Space-Frequency Index Modulation (PT-GSFIM), where part of the information bits select the active antennas and subcarriers which then carry amplitude and phase modulated symbols. The proposed scheme is designed for multiuser multiple-input multiple-output (MU-MIMO) scenarios and incorporates a precoder which removes multiuser interference (MUI) at the receivers. Furthermore, the proposed PT-GSFIM also integrates signal space diversity (SSD) techniques for tackling the typical poor performance of uncoded orthogonal frequency division multiplexing (OFDM) based schemes. By combining complex rotation matrices (CRM) and subcarrier-level interleaving, PT-GSFIM can exploit the inherent diversity in frequency selective channels and improve the performance without additional power or bandwidth. To support reliable detection of the multidimensional PT-GSFIM we also propose three different detection algorithms which can provide different tradeoffs between performance and complexity. Simulation results shows that proposed PT-GSFIM scheme, can provide significant gains over conventional MU-MIMO and GSM schemes.


title: Investigating The Impacting Factors on The Public's Attitudes Towards Autonomous Vehicles Using Sentiment Analysis from Social Media Data

abstract: The public's attitudes play a critical role in the acceptance, purchase, use, and research and development of autonomous vehicles (AVs). To date, the public's attitudes towards AVs were mostly estimated through traditional survey data with high labor costs and a low quantity of samples, which also might be one of the reasons why the influencing factors on the public's attitudes of AVs have not been studied from multiple aspects in a comprehensive way yet. To address the issue, this study aims to propose a method by using large-scale social media data to investigate key factors that affect the public's attitudes and acceptance of AVs. A total of 954,151 Twitter data related to AVs and 53 candidate independent variables from seven categories were extracted using the web scraping method. Then, sentiment analysis was used to measure the public attitudes towards AVs by calculating sentiment scores. Random forests algorithm was employed to preliminarily select candidate independent variables according to their importance, while a linear mixed model was performed to explore the impacting factors considering the unobserved heterogeneities caused by the subjectivity level of tweets. The results showed that the overall attitude of the public on AVs was slightly optimistic. Factors like "drunk", "blind spot", and "mobility" had the largest impacts on public attitudes. In addition, people were more likely to express positive feelings when talking about words such as "lidar" and "Tesla" that relate to high technologies. Conversely, factors such as "COVID-19", "pedestrian", "sleepy", and "highway" were found to have significantly negative effects on the public's attitudes. The findings of this study are beneficial for the development of AV technologies, the guidelines for AV-related policy formulation, and the public's understanding and acceptance of AVs.


title: Isoperimetric inequalities in cylinders with density

abstract: Given a compact Riemannian manifold with density $M$ without boundary and the real line $\mathbb{R}$ with constant density, we prove that isoperimetric regions of large volume in $M\times\mathbb{R}$ with the product density are slabs of the form $M\times [a,b]$. We previously prove, as a necessary step, the existence of isoperimetric regions in any manifold of density where a subgroup of the group of transformations preserving weighted perimeter and volume acts cocompactly.


title: HAWC measurements of the energy spectra of cosmic ray protons, helium and heavy nuclei in the TeV range

abstract: Current knowledge of the relative abundances and the energy spectra of the elemental mass groups of cosmic rays in the $10$ TeV - $1$ PeV interval is uncertain. This situation prevents carrying out precision tests that may lead to distinguish among the existing hypotheses on the origin and propagation of TeV cosmic rays in the galaxy. In order to learn more about the mass composition of these particles, we have employed HAWC data from hadron induced air showers in order to determine the spectra of three mass groups of cosmic rays: protons, helium and heavy nuclei with $Z > 2$. The energy spectra were estimated by using the Gold unfolding technique on the 2D distribution of the lateral shower age against the estimated primary energy of events with arrival zenith angles smaller than 45 degrees. The study was carried out based on simulations using the QGSJET-II-04 model. Results are presented for primary cosmic-ray energies from $10$ TeV to $251$ TeV. They reveal that the aforementioned cosmic ray spectra exhibit fine structures within the above primary energy range.


title: Robust fractionation in cancer radiotherapy

abstract: In cancer radiotherapy, the standard formulation of the optimal fractionation problem based on the linear-quadratic dose-response model is a non-convex quadratically constrained quadratic program (QCQP). An optimal solution for this QCQP can be derived by solving a two-variable linear program. Feasibility of this solution, however, crucially depends on the so-called alphaover-beta ratios for the organs-at-risk, whose true values are unknown. Consequently, the dosing schedule presumed optimal, in fact, may not even be feasible in practice. We address this by proposing a robust counterpart of the nominal formulation. We show that a robust solution can be derived by solving a small number of two-variable linear programs, each with a small number of constraints. We quantify the price of robustness, and compare the incidence and extent of infeasibility of the nominal and robust solutions via numerical experiments.


title: Regression Diagnostics meets Forecast Evaluation: Conditional Calibration, Reliability Diagrams, and Coefficient of Determination

abstract: Model diagnostics and forecast evaluation are two sides of the same coin. A common principle is that fitted or predicted distributions ought to be calibrated or reliable, ideally in the sense of auto-calibration, where the outcome is a random draw from the posited distribution. For binary responses, this is the universal concept of reliability. For real-valued outcomes, a general theory of calibration has been elusive, despite a recent surge of interest in distributional regression and machine learning. We develop a framework rooted in probability theory, which gives rise to hierarchies of calibration, and applies to both predictive distributions and stand-alone point forecasts. In a nutshell, a prediction - distributional or single-valued - is conditionally T-calibrated if it can be taken at face value in terms of the functional T. Whenever T is defined via an identification function - as in the cases of threshold (non) exceedance probabilities, quantiles, expectiles, and moments - auto-calibration implies T-calibration. We introduce population versions of T-reliability diagrams and revisit a score decomposition into measures of miscalibration (MCB), discrimination (DSC), and uncertainty (UNC). In empirical settings, stable and efficient estimators of T-reliability diagrams and score components arise via nonparametric isotonic regression and the pool-adjacent-violators algorithm. For in-sample model diagnostics, we propose a universal coefficient of determination, $$\text{R}^\ast = \frac{\text{DSC}-\text{MCB}}{\text{UNC}},$$ that nests and reinterprets the classical $\text{R}^2$ in least squares (mean) regression and its natural analogue $\text{R}^1$ in quantile regression, yet applies to T-regression in general, with MCB $\geq 0$, DSC $\geq 0$, and $\text{R}^\ast \in [0,1]$ under modest conditions.