title: Fault Localization with Code Coverage Representation Learning

abstract: In this paper, we propose DeepRL4FL, a deep learning fault localization (FL) approach that locates the buggy code at the statement and method levels by treating FL as an image pattern recognition problem. DeepRL4FL does so via novel code coverage representation learning (RL) and data dependencies RL for program statements. Those two types of RL on the dynamic information in a code coverage matrix are also combined with the code representation learning on the static information of the usual suspicious source code. This combination is inspired by crime scene investigation in which investigators analyze the crime scene (failed test cases and statements) and related persons (statements with dependencies), and at the same time, examine the usual suspects who have committed a similar crime in the past (similar buggy code in the training data). For the code coverage information, DeepRL4FL first orders the test cases and marks error-exhibiting code statements, expecting that a model can recognize the patterns discriminating between faulty and non-faulty statements/methods. For dependencies among statements, the suspiciousness of a statement is seen taking into account the data dependencies to other statements in execution and data flows, in addition to the statement by itself. Finally, the vector representations for code coverage matrix, data dependencies among statements, and source code are combined and used as the input of a classifier built from a Convolution Neural Network to detect buggy statements/methods. Our empirical evaluation shows that DeepRL4FL improves the top-1 results over the state-of-the-art statement-level FL baselines from 173.1% to 491.7%. It also improves the top-1 results over the existing method-level FL baselines from 15.0% to 206.3%.


title: Quantification of Bore Path Uncertainty in Borehole Heat Exchanger Arrays

abstract: Borehole heat exchanger arrays have become a common implement for the utilization of thermal energy in the soil. Building these facilities is expensive, especially the drilling of boreholes, into which closed-pipe heat exchangers are inserted. Therefore, cost-reducing drilling methods are common practice, which can produce inaccuracies of varying degree. This brings into question how much these inaccuracies could potentially affect the performance of a planned system. In the presented case study, an uncertainty quantification for seasonally operated borehole heat exchanger arrays is performed to analyze the bore paths' deviations impact. We introduce an adaptive, anisotropic stochastic collocation method, known as the generalized Smolyak algorithm, which was previously unused in this context and apply it to a numerical model of the borehole heat exchanger array. Our results show that the borehole heat exchanger array performance is surprisingly reliable even with potentially severe implementation errors during their construction. This, coupled with the potential uses of the presented method in similar applications gives planners and investors valuable information regarding the viability of borehole heat exchanger arrays in the face of uncertainty. With this paper, we hope to provide a powerful statistical tool to the field of geothermal energy, in which uncertainty quantification methods are still rarely used at this point. The discussed case study represents a jumping-off point for further investigations on the effects of uncertainty on borehole heat exchanger arrays and borehole thermal energy storage systems.


title: Strong bulk-surface interaction dominated in-plane anisotropy of electronic structure in GaTe

abstract: Recently, intriguing physical properties have been unraveled in anisotropic layered semiconductors, in which the in-plane electronic band structure anisotropy often originates from the low crystallographic symmetry and thus a thickness-independent character emerges. Here, we apply high-resolution angle-resolved photoemission spectroscopy to directly image the in-plane anisotropic energy bands in monoclinic gallium telluride (GaTe). Our first-principles calculations reveal the in-plane anisotropic energy band structure of GaTe measured experimentally is dominated by a strong bulk-surface interaction rather than geometric factors, surface effect and quantum confinement effect. Furthermore, accompanied by the thickness of GaTe increasing from mono- to few-layers, the strong interlayer coupling of GaTe induces direct-indirect-direct band gap transitions and the in-plane anisotropy of hole effective mass is reversed. Our results shed light on the physical origins of in-plane anisotropy of electronic structure in GaTe, paving the way for the design and device applications of nanoelectronics and optoelectronics based on anisotropic layered semiconductors.


title: Singularity-aware motion planning for multi-axis additive manufacturing

abstract: Multi-axis additive manufacturing enables high flexibility of material deposition along dynamically varied directions. The Cartesian motion platforms of these machines include three parallel axes and two rotational axes. Singularity on rotational axes is a critical issue to be tackled in motion planning for ensuring high quality of manufacturing results. The highly nonlinear mapping in the singular region can convert a smooth toolpath with uniformly sampled waypoints defined in the model coordinate system into a highly discontinuous motion in the machine coordinate system, which leads to over-extrusion / under-extrusion of materials in filament-based additive manufacturing. The problem is challenging as both the maximal and the minimal speeds at the tip of a printer head must be controlled in motion. Moreover, collision may occur when sampling-based collision avoidance is employed. In this paper, we present a motion planning method to support the manufacturing realization of designed toolpaths for multi-axis additive manufacturing. Problems of singularity and collision are considered in an integrated manner to improve the motion therefore the quality of fabrication.


title: PA-ResSeg: A Phase Attention Residual Network for Liver Tumor Segmentation from Multi-phase CT Images

abstract: In this paper, we propose a phase attention residual network (PA-ResSeg) to model multi-phase features for accurate liver tumor segmentation, in which a phase attention (PA) is newly proposed to additionally exploit the images of arterial (ART) phase to facilitate the segmentation of portal venous (PV) phase. The PA block consists of an intra-phase attention (Intra-PA) module and an inter-phase attention (Inter-PA) module to capture channel-wise self-dependencies and cross-phase interdependencies, respectively. Thus it enables the network to learn more representative multi-phase features by refining the PV features according to the channel dependencies and recalibrating the ART features based on the learned interdependencies between phases. We propose a PA-based multi-scale fusion (MSF) architecture to embed the PA blocks in the network at multiple levels along the encoding path to fuse multi-scale features from multi-phase images. Moreover, a 3D boundary-enhanced loss (BE-loss) is proposed for training to make the network more sensitive to boundaries. To evaluate the performance of our proposed PA-ResSeg, we conducted experiments on a multi-phase CT dataset of focal liver lesions (MPCT-FLLs). Experimental results show the effectiveness of the proposed method by achieving a dice per case (DPC) of 0.77.87, a dice global (DG) of 0.8682, a volumetric overlap error (VOE) of 0.3328 and a relative volume difference (RVD) of 0.0443 on the MPCT-FLLs. Furthermore, to validate the effectiveness and robustness of PA-ResSeg, we conducted extra experiments on another multi-phase liver tumor dataset and obtained a DPC of 0.8290, a DG of 0.9132, a VOE of 0.2637 and a RVD of 0.0163. The proposed method shows its robustness and generalization capability in different datasets and different backbones.


title: Four point functions in CFT's with slightly broken higher spin symmetry

abstract: We compute spinning four point functions in the quasi-fermionic three dimensional conformal field theory with slightly broken higher spin symmetry at finite t'Hooft coupling. More concretely, we obtain a formula for $\langle j_s j_{\tilde{0}} j_{\tilde{0}} j_{\tilde{0}} \rangle$, where $j_s$ is a higher spin current and $j_{\tilde{0}}$ is the scalar single trace operator. Our procedure consists in writing a plausible ansatz in Mellin space and using crossing, pseudo-conservation and Regge boundedness to fix all undetermined coefficients. Our method can potentially be generalised to compute all spinning four point functions in these theories.


title: An astronomical image content-based recommendation system using combined deep learning models in a fully unsupervised mode

abstract: We have developed a method that maps large astronomical images onto a two-dimensional map and clusters them. A combination of various state-of-the-art machine learning (ML) algorithms is used to develop a fully unsupervised image quality assessment and clustering system. Our pipeline consists of a data pre-processing step where individual image objects are identified in a large astronomical image and converted to smaller pixel images. This data is then fed to a deep convolutional autoencoder jointly trained with a self-organizing map (SOM). This part can be used as a recommendation system. The resulting output is eventually mapped onto a two-dimensional grid using a second, deep, SOM. We use data taken from ground-based telescopes and, as a case study, compare the system's ability and performance with the results obtained by supervised methods presented by Teimoorinia et al. (2020). The availability of target labels in this data allowed a comprehensive performance comparison between our unsupervised and supervised methods. In addition to image-quality assessments performed in this project, our method can have various other applications. For example, it can help experts label images in a considerably shorter time with minimum human intervention. It can also be used as a content-based recommendation system capable of filtering images based on the desired content.


title: Unscented Kalman Inversion: Efficient Gaussian Approximation to the Posterior Distribution

abstract: The unscented Kalman inversion (UKI) method presented in [1] is a general derivative-free approach for the inverse problem. UKI is particularly suitable for inverse problems where the forward model is given as a black box and may not be differentiable. The regularization strategies, convergence property, and speed-up strategies [1,2] of the UKI are thoroughly studied, and the method is capable of handling noisy observation data and solving chaotic inverse problems. In this paper, we study the uncertainty quantification capability of the UKI. We propose a modified UKI, which allows to well approximate the mean and covariance of the posterior distribution for well-posed inverse problems with large observation data. Theoretical guarantees for both linear and nonlinear inverse problems are presented. Numerical results, including learning of permeability parameters in subsurface flow and of the Navier-Stokes initial condition from solution data at positive times are presented. The results obtained by the UKI require only $O(10)$ iterations, and match well with the expected results obtained by the Markov Chain Monte Carlo method.


title: Right-cancellable protomodular algebras

abstract: A new protomodular analog of the classical criterion for the existence of a group term in the algebraic theory of a variety of universal algebras is given. To this end, the notion of a right-cancellable protomodular algebra is introduced. It is proved that the algebraic theory of a variety of universal algebras contains a group term if and only if it contains protomodular terms with respect to which all algebras from the variety are right-cancellable. This, in particular, gives a partial answer to the extended version of an open problem from loop theory whether any Hausdorff topological (semi-)loop is completely regular. Moreover, the right-cancellable algebras from the simplest protomodular varieties are characterized as sets with principal group actions as well as groups with simple additional structures.


title: Field-effect chirality devices with Dirac semimetal

abstract: Charge-based field-effect transistors (FETs) greatly suffer from unavoidable carrier scattering and heat dissipation. In analogy to valley degree of freedom in semiconductors, chiral anomaly current in Weyl/Dirac semimetals is theoretically predicted to be nearly non-dissipative over long distances, but still lacks experimental ways to efficiently control its transport. Here we demonstrate field-effect chirality devices with Dirac semimetal PtSe2, in which its Fermi level is close to the Dirac point in conduction band owing to intrinsic defects. The chiral anomaly is further corroborated with nonlocal valley transport measurement, which can also be effectively modulated by external fields, showing robust nonlocal valley transport with micrometer diffusion length. Similar to charge-based FETs, the chiral conductivity in PtSe2 devices can be modulated by electrostatic gating with an ON/OFF ratio more than 103. We also demonstrate basic logic functions in the devices with electric and magnetic fields as input signals.


title: Quasistationary Distributions and Ergodic Control Problems

abstract: We introduce and study the basic properties of two ergodic stochastic control problems associated with the quasistationary distribution (QSD) of a diffusion process $X$ relative to a bounded domain. The two problems are in some sense dual, with one defined in terms of the generator associated with $X$ and the other in terms of its adjoint. Besides proving wellposedness of the associated Hamilton-Jacobi-Bellman equations, we describe how they can be used to characterize important properties of the QSD. Of particular note is that the QSD itself can be identified, up to normalization, in terms of the cost potential of the control problem associated with the adjoint.


title: Gravitational Wave Direct Detection does not Constrain the Tensor Spectral Index at CMB Scales

abstract: I discuss constraints on the power spectrum of primordial tensor perturbations from a combination of Cosmic Microwave Background (CMB) measurements and the gravitational wave direct detection experiments LIGO/Virgo and DECIGO. There are two main points: (1) Inflation predicts an approximately power-law form for the primordial tensor spectrum, but makes no prediction for its amplitude. Given that neither Planck nor LIGO/Virgo has actually detected primordial tensor modes, it is trivially true that no model-independent constraint on the slope of the tensor power spectrum is possible with current data. (2) CMB and LIGO/Virgo scales differ by more than 19 orders of magnitude, and 16 for DECIGO. I show that a power-law extrapolation from CMB to direct detection frequencies overestimates the amplitude of primordial tensor modes by as much as two orders of magnitude relative to an ensemble of realistic single-field inflation models. Moreover, the primordial tensor amplitude at direct detection scales is mostly uncorrelated with the tensor spectral index at CMB scales, and any constraint is strongly dependent on the specific form of the inflationary potential.


title: A number theoretic characterization of $E$-smooth and (FRS) morphisms: estimates on the number of $\mathbb{Z}/p^{k}\mathbb{Z}$-points

abstract: We provide uniform estimates on the number of $\mathbb{Z}/p^{k}\mathbb{Z}$-points lying on fibers of flat morphisms between smooth varieties whose fibers have rational singularities, termed (FRS) morphisms. For each individual fiber, the estimates were known by work of Avni and Aizenbud, but we render them uniform over all fibers. The proof technique for individual fibers is based on Hironaka's resolution of singularities and Denef's formula, but breaks down in the uniform case. Instead, we use recent results from the theory of motivic integration. Our estimates are moreover equivalent to the (FRS) property, just like in the absolute case by Avni and Aizenbud. In addition, we define new classes of morphisms, called $E$-smooth morphisms ($E\in\mathbb{N}$), which refine the (FRS) property, and use the methods we developed to provide uniform number-theoretic estimates as above for their fibers. Similar estimates are given for fibers of $\varepsilon$-jet flat morphisms, improving previous results by the last two authors.


title: Linearly ordered sets with only one operator have the amalgamation property

abstract: The class of linearly ordered sets with one order preserving unary operation has the Strong Amalgamation Property (SAP). The class of linearly ordered sets with one strict order preserving unary operation has AP but not SAP. The class of linearly ordered sets with two order preserving unary operations has not AP. For every set $F$, the class of linearly ordered sets with an $F$-indexed family of automorphisms has SAP. Corresponding results are proved in the case of order reversing operations. Various subclasses of the above classes are considered and some model-theoretical consequences are presented.


title: On the Initialization for Convex-Concave Min-max Problems

abstract: Convex-concave min-max problems are ubiquitous in machine learning, and people usually utilize first-order methods (e.g., gradient descent ascent) to find the optimal solution. One feature which separates convex-concave min-max problems from convex minimization problems is that the best known convergence rates for min-max problems have an explicit dependence on the size of the domain, rather than on the distance between initial point and the optimal solution. This means that the convergence speed does not have any improvement even if the algorithm starts from the optimal solution, and hence, is oblivious to the initialization. Here, we show that strict-convexity-strict-concavity is sufficient to get the convergence rate to depend on the initialization. We also show how different algorithms can asymptotically achieve initialization-dependent convergence rates on this class of functions. Furthermore, we show that the so-called "parameter-free" algorithms allow to achieve improved initialization-dependent asymptotic rates without any learning rate to tune. In addition, we utilize this particular parameter-free algorithm as a subroutine to design a new algorithm, which achieves a novel non-asymptotic fast rate for strictly-convex-strictly-concave min-max problems with a growth condition and H{\"o}lder continuous solution mapping. Experiments are conducted to verify our theoretical findings and demonstrate the effectiveness of the proposed algorithms.


title: Visual Navigation with a 2-pixel Camera---Possibilities and Limitations

abstract: Borrowing terminology from fluid mechanics, the concepts of {\em Eulerian} and {\em Lagrangian optical flow sensing} are introduced. Eulerian optical flow sensing assumes that each photoreceptor in the camera or eye can instantaneously detect feature image points and their velocities on the retina. If this assumption is satisfied, even a two pixel imaging system can provide a moving agent with information about its movement along a corridor that is sufficiently precise as to be used as a robustly reliable steering signal. Implementing Eulerian optical flow sensing poses significant challenges, however. Lagrangian optical flow, on the other hand, tracks feature image points as they move on the retina. This form of visual sensing is the basis for many standard computer vision implementations, including Lukas-Kanade and Horn-Schunck. Lagrangian optical flow has its own challenges, not least of which is that it is badly confounded by rotational components of motion. Combined steering and sensing strategies for mitigating the effects of rotational motions are considered.


title: A Novel Adaptive Deep Network for Building Footprint Segmentation

abstract: Building footprint segmentations for high resolution images are increasingly demanded for many remote sensing applications. By the emerging deep learning approaches, segmentation networks have made significant advances in the semantic segmentation of objects. However, these advances and the increased access to satellite images require the generation of accurate object boundaries in satellite images. In the current paper, we propose a novel network-based on Pix2Pix methodology to solve the problem of inaccurate boundaries obtained by converting satellite images into maps using segmentation networks in order to segment building footprints. To define the new network named G2G, our framework includes two generators where the first generator extracts localization features in order to merge them with the boundary features extracted from the second generator to segment all detailed building edges. Moreover, different strategies are implemented to enhance the quality of the proposed networks' results, implying that the proposed network outperforms state-of-the-art networks in segmentation accuracy with a large margin for all evaluation metrics. The implementation is available at https://github.com/A2Amir/A-Novel-Adaptive-Deep-Network-for-Building-Footprint-Segmentation.


title: Automated Generation of Interorganizational Disaster Response Networks through Information Extraction

abstract: When a disaster occurs, maintaining and restoring community lifelines subsequently require collective efforts from various stakeholders. Aiming at reducing the efforts associated with generating Stakeholder Collaboration Networks (SCNs), this paper proposes a systematic approach to reliable information extraction for stakeholder collaboration and automated network generation. Specifically, stakeholders and their interactions are extracted from texts through Named Entity Recognition (NER), one of the techniques in natural language processing. Once extracted, the collaboration information is transformed into structured datasets to generate the SCNs automatically. A case study of stakeholder collaboration during Hurricane Harvey was investigated and it had demonstrated the feasibility and applicability of the proposed method. Hence, the proposed approach was proved to significantly reduce practitioners' interpretation and data collection workloads. In the end, discussions and future work are provided.


title: On Optimizing the Trade-off between Privacy and Utility in Data Provenance

abstract: Organizations that collect and analyze data may wish or be mandated by regulation to justify and explain their analysis results. At the same time, the logic that they have followed to analyze the data, i.e., their queries, may be proprietary and confidential. Data provenance, a record of the transformations that data underwent, was extensively studied as means of explanations. In contrast, only a few works have studied the tension between disclosing provenance and hiding the underlying query. This tension is the focus of the present paper, where we formalize and explore for the first time the tradeoff between the utility of presenting provenance information and the breach of privacy it poses with respect to the underlying query. Intuitively, our formalization is based on the notion of provenance abstraction, where the representation of some tuples in the provenance expressions is abstracted in a way that makes multiple tuples indistinguishable. The privacy of a chosen abstraction is then measured based on how many queries match the obfuscated provenance, in the same vein as k-anonymity. The utility is measured based on the entropy of the abstraction, intuitively how much information is lost with respect to the actual tuples participating in the provenance. Our formalization yields a novel optimization problem of choosing the best abstraction in terms of this tradeoff. We show that the problem is intractable in general, but design greedy heuristics that exploit the provenance structure towards a practically efficient exploration of the search space. We experimentally prove the effectiveness of our solution using the TPC-H benchmark and the IMDB dataset.


title: Hyperpolarization of cis-15N,15N'-azobenzene by parahydrogen at ultralow magnetic fields

abstract: Development of the methods to exploit nuclear hyperpolarization and search for molecules whose nuclear spins can be efficiently hyperpolarized is an active area in nuclear magnetic resonance. Of particular interest are those molecules that have long nuclear relaxation times, making them to be suitable candidates as contrast agents in magnetic resonance imaging. In this work, we present a detailed study of SABRE SHEATH (Signal Amplification By Reversible Exchange in Shield Enabled Alignment Transfer to Heteronuclei) experiments of 15N,15N' azobenzene. In SABRE SHEATH experiments nuclear spins of the target are hyperpolarized by transfer of spin polarization from parahydrogen at ultralow fields during a reversible chemical process. The studied system is complicated, and we are concerned only about a subset of the data, presenting details for the molecules that experience fast chemical exchange at the catalytic complex and thus are involved in polarizing the free azobenzene. Azobenzene exists in two isomers trans- and cis-. We show that all nuclear spins in cis-azobenzene can be efficiently hyperpolarized by SABRE at suitable magnetic fields. Enhancement factors (relative to 9.4 T) reach several thousands of times for 15N spins and a few tens of times for the 1H spins. There are two approaches to observe either hyperpolarized magnetization of 15N/1H spins or hyperpolarized singlet order of the 15N spin pair. We compare these approaches and present the field dependencies of SABRE experiments for them. No hyperpolarization of trans 15N,15N' azobenzene was observed. The results presented here will be useful for further experiments where hyperpolarized cis-15N,15N' azobenzene is switched by light to trans 15N,15N' azobenzene for storing the produced hyperpolarization in the long-lived spin state of the 15N pair of trans-15N,15N' azobenzene.


title: Jenss-Bayley Latent Change Score Model with Individual Ratio of Growth Acceleration in the Framework of Individual Measurement Occasions

abstract: Longitudinal data analysis has been widely employed to examine between-individual differences in within-individual changes. One challenge of such analyses is that the rate-of-change is only available indirectly when change patterns are nonlinear with respect to time. Latent change score models (LCSMs), which can be employed to investigate the change in rate-of-change at the individual level, have been developed to address this challenge. We extend an existing LCSM with the Jenss-Bayley growth curve \cite[Chapter~18]{Grimm2016growth} and propose a novel expression for change scores that allows for (1) unequally-spaced study waves and (2) individual measurement occasions around each wave. We also extend the existing model to estimate the individual ratio of the growth acceleration (that largely determines the trajectory shape and is viewed as the most important parameter in the Jenss-Bayley model). We present the proposed model by a simulation study and a real-world data analysis. Our simulation study demonstrates that the proposed model can estimate the parameters unbiasedly and precisely and exhibit target confidence interval coverage. The simulation study also shows that the proposed model with the novel expression for the change scores outperforms the existing model. An empirical example using longitudinal reading scores shows that the model can estimate the individual ratio of the growth acceleration and generate individual rate-of-change in practice. We also provide the corresponding code for the proposed model.


title: Depths of Copernican craters on lunar maria and highlands

abstract: We present a study on the relationship between the ratio of the depth of a crater to its diameter and the diameter for lunar craters both on the maria and on the highlands. We consider craters younger than 1.1 billion years in age, i.e. of Copernican period. The aim of this work is to improve our understanding of such relationships based on our new estimates of the craters's depth and diameter. Previous studies considered similar relationships for much older craters (up to 3.2 billion years). We calculated the depths of craters with diameters from 10 to 100 km based on the altitude profiles derived from data obtained by the Lunar Orbiter Laser Altimeter (LOLA) onboard the Lunar Reconnaissance Orbiter (LRO). The ratio h/D of the depth h of a crater to its diameter D can diverge by up to a factor of two for craters with almost the same diameters. The linear and power approximations (regressions) of the dependence of h/D on D were made for simple and complex Copernican craters selected from the data from Mazrouei et al. (2019) and Losiak et al. (2015). For the separation of highland craters into two groups based only on their dependences of h/D on D, at D<18 km these are mostly simple craters, although some complex craters can have diameters D>16 km. Depths of mare craters with D<14 km are greater than 0.15D. Following Pike's (1981) classification, we group mare craters of D<15 km as simple craters. Mare craters with 15<D<18 km fit both approximation curves for simple and complex craters. Depths of mare craters with D>18 km are in a better agreement with the approximation curve of h/D vs. D for complex craters than for simple craters. At the same diameter, mare craters are deeper than highland craters at a diameter smaller than 30-40 km. For greater diameters, highland craters are deeper.


title: A ring accelerator? Unusual jet dynamics in the IceCube candidate PKS 1502+106

abstract: On 2019/07/30.86853 UT, IceCube detected a high-energy astrophysical neutrino candidate. The Flat Spectrum Radio Quasar PKS 1502+106 is located within the 50 percent uncertainty region of the event. Our analysis of 15 GHz Very Long Baseline Array (VLBA) and astrometric 8 GHz VLBA data, in a time span prior and after the IceCube event, reveals evidence for a radio ring structure which develops with time. Several arc-structures evolve perpendicular to the jet ridge line. We find evidence for precession of a curved jet based on kinematic modelling and a periodicity analysis. An outflowing broad line region (BLR) based on the C IV line emission (Sloan Digital Sky Survey, SDSS) is found. We attribute the atypical ring to an interaction of the precessing jet with the outflowing material. We discuss our findings in the context of a spine-sheath scenario where the ring reveals the sheath and its interaction with the surroundings (narrow line region, NLR, clouds). We find that the radio emission is correlated with the $\gamma$-ray emission, with radio lagging the $\gamma$-rays. Based on the $\gamma$-ray variability timescale, we constrain the $\gamma$-ray emission zone to the BLR (30-200 $r_{\rm g}$) and within the jet launching region. We discuss that the outflowing BLR provides the external radiation field for $\gamma$-ray production via external Compton scattering. The neutrino is most likely produced by proton-proton interaction in the blazar zone (beyond the BLR), enabled by episodic encounters of the jet with dense clouds, i.e. some molecular cloud in the NLR.


title: N-Shot Learning for Augmenting Task-Oriented Dialogue State Tracking

abstract: Augmentation of task-oriented dialogues has followed standard methods used for plain-text such as back-translation, word-level manipulation, and paraphrasing despite its richly annotated structure. In this work, we introduce an augmentation framework that utilizes belief state annotations to match turns from various dialogues and form new synthetic dialogues in a bottom-up manner. Unlike other augmentation strategies, it operates with as few as five examples. Our augmentation strategy yields significant improvements when both adapting a DST model to a new domain, and when adapting a language model to the DST task, on evaluations with TRADE and TOD-BERT models. Further analysis shows that our model performs better on seen values during training, and it is also more robust to unseen values. We conclude that exploiting belief state annotations enhances dialogue augmentation and results in improved models in n-shot training scenarios.


title: Extremal general affine surface areas

abstract: For a convex body $K$ in $\mathbb{R}^n$, we introduce and study the extremal general affine surface areas, defined by \[ {\rm IS}_{\varphi}(K):=\sup_{K^\prime\subset K}{\rm as}_{\varphi}(K),\quad {\rm os}_{\psi}(K):=\inf_{K^\prime\supset K}{\rm as}_{\psi}(K) \] where ${\rm as}_{\varphi}(K)$ and ${\rm as}_{\psi}(K)$ are the $L_\varphi$ and $L_\psi$ affine surface area of $K$, respectively. We prove that there exist extremal convex bodies that achieve the supremum and infimum, and that the functionals ${\rm IS}_{\varphi}$ and ${\rm os}_{\psi}$ are continuous. In our main results, we prove Blaschke-Santal\'o type inequalities and inverse Santal\'o type inequalities for the extremal general affine surface areas. This article may be regarded as an Orlicz extension of the recent work of Giladi, Huang, Sch\"utt and Werner (2020), who introduced and studied the extremal $L_p$ affine surface areas.


title: Nash equilibrium mapping vs Hamiltonian dynamics vs Darwinian evolution for some social dilemma games in the thermodynamic limit

abstract: How cooperation evolves and manifests itself in the thermodynamic or infinite player limit of social dilemma games is a matter of intense speculation. Various analytical methods have been proposed to analyze the thermodynamic limit of social dilemmas. In this work, we compare two analytical methods, i.e., Darwinian evolution and Nash equilibrium mapping, with a numerical agent-based approach. For completeness, we also give results for another analytical method, Hamiltonian dynamics. In contrast to Hamiltonian dynamics, which involves the maximization of payoffs of all individuals, in Darwinian evolution, the payoff of a single player is maximized with respect to its interaction with the nearest neighbor. While the Hamiltonian dynamics method utterly fails as compared to Nash equilibrium mapping, the Darwinian evolution method gives a false positive for game magnetization -- the net difference between the fraction of cooperators and defectors -- when payoffs obey the condition a + d = b + c, wherein a,d represents the diagonal elements and b,c the off-diagonal elements in a symmetric social dilemma game payoff matrix. When either a + d =/= b + c or when one looks at the average payoff per player, the Darwinian evolution method fails, much like the Hamiltonian dynamics approach. On the other hand, the Nash equilibrium mapping and numerical agent-based method agree well for both game magnetization and average payoff per player for the social dilemmas in question, i.e., the Hawk-Dove game and the Public goods game. This paper thus brings to light the inconsistency of the Darwinian evolution method vis-a-vis both Nash equilibrium mapping and a numerical agent-based approach.

categories: cond-mat.stat-mech econ.TH physics.soc-ph q-bio.PE


title: Hyperuniformity and phase enrichment in vortex and rotor assemblies

abstract: Ensembles of particles rotating in a two-dimensional fluid can exhibit chaotic dynamics yet develop signatures of hidden order. Such "rotors" are found in the natural world spanning vastly disparate length scales - from the rotor proteins in cellular membranes to models of atmospheric dynamics. Here we show that an initially random distribution of either ideal vortices in an inviscid fluid, or driven rotors in a viscous membrane, spontaneously self assembles. Despite arising from drastically different physics, these systems share a Hamiltonian structure that sets geometrical conservation laws resulting in distinct structural states. We find that the rotationally invariant interactions isotropically suppress long wavelength fluctuations - a hallmark of a disordered hyperuniform material. With increasing area fraction, the system orders into a hexagonal lattice. In mixtures of two co-rotating populations, the stronger population will gain order from the other and both will become phase enriched. Finally, we show that classical 2D point vortex systems arise as exact limits of the experimentally accessible microscopic membrane rotors, yielding a new system through which to study topological defects.


title: Unrealizable Cores for Reactive Systems Specifications

abstract: One of the main challenges of reactive synthesis, an automated procedure to obtain a correct-by-construction reactive system, is to deal with unrealizable specifications. One means to deal with unrealizability, in the context of GR(1), an expressive assume-guarantee fragment of LTL that enables efficient synthesis, is the computation of an unrealizable core, which can be viewed as a fault-localization approach. Existing solutions, however, are computationally costly, are limited to computing a single core, and do not correctly support specifications with constructs beyond pure GR(1) elements. In this work we address these limitations. First, we present QuickCore, a novel algorithm that accelerates unrealizable core computations by relying on the monotonicity of unrealizability, on an incremental computation, and on additional properties of GR(1) specifications. Second, we present Punch, a novel algorithm to efficiently compute all unrealizable cores of a specification. Finally, we present means to correctly handle specifications that include higher-level constructs beyond pure GR(1) elements. We implemented our ideas on top of Spectra, an open-source language and synthesis environment. Our evaluation over benchmarks from the literature shows that QuickCore is in most cases faster than previous algorithms, and that its relative advantage grows with scale. Moreover, we found that most specifications include more than one core, and that Punch finds all the cores significantly faster than a competing naive algorithm.


title: Observing quantum coherence from photons scattered in free-space

abstract: Quantum channels in free-space, an essential prerequisite for fundamental tests of quantum mechanics and quantum technologies in open space, have so far been based on direct line-of-sight because the predominant approaches for photon-encoding, including polarization and spatial modes, are not compatible with randomly scattered photons. Here we demonstrate a novel approach to transfer and recover quantum coherence from scattered, non-line-of-sight photons analyzed in a multimode and imaging interferometer for time-bins, combined with photon detection based on a 8x8 single-photon-detector-array. The observed time-bin visibility for scattered photons remained at a high $95\%$ over a wide scattering angle range of -45 degree to +45 degree, while the individual pixels in the detector array resolve or track an image in its field of view of ca. 0.5 degrees. Using our method we demonstrate the viability of two novel applications. Firstly, using scattered photons as an indirect channel for quantum communication thereby enabling non-line-of-sight quantum communication with background suppression, and secondly, using the combined arrival time and quantum coherence to enhance the contrast of low-light imaging and laser ranging under high background light. We believe our method will instigate new lines for research and development on applying photon coherence from scattered signals to quantum sensing, imaging, and communication in free-space environments.


title: Primal-Dual Stochastic Mirror Descent for MDPs

abstract: We consider the problem of learning the optimal policy for infinite-horizon Markov decision processes (MDPs). For this purpose, some variant of Stochastic Mirror Descent is proposed for convex programming problems with Lipschitz-continuous functionals. An important detail is the ability to use inexact values of functional constraints and compute the value of dual variables. We analyze this algorithm in a general case and obtain an estimate of the convergence rate that does not accumulate errors during the operation of the method. Using this algorithm, we get the first parallel algorithm for mixing average-reward MDPs with a generative model without reduction to discounted MDP. One of the main features of the presented method is low communication costs in a distributed centralized setting, even with very large networks.


title: Observation of Possible Nonlinear Anomalous Hall Effect in Organic Two-Dimensional Dirac Fermion System

abstract: We report the observation of nonlinear anomalous Hall effect (NLAHE) in the multilayered organic conductor $\alpha$-(BEDT-TTF)$_2$I$_3$ in the charge order (CO) insulating phase just under the critical pressure for transition into two-dimensional (2D) massless Dirac fermion (DF) phase. We successfully extracted the finite nonlinear Hall voltage proportional to square current at zero magnetic field. The observed NLAHE features, current direction dependence and correlation with CO, are consistent with the previous estimation assuming 2D massive DF with a pair of tilted Dirac cones. This is the first observation of topological transport in organic conductors, and also the first example of NLAHE in the electronic phase with spontaneous symmetry breaking.


title: Spline parameterization of neural network controls for deep learning

abstract: Based on the continuous interpretation of deep learning cast as an optimal control problem, this paper investigates the benefits of employing B-spline basis functions to parameterize neural network controls across the layers. Rather than equipping each layer of a discretized ODE-network with a set of trainable weights, we choose a fixed number of B-spline basis functions whose coefficients are the trainable parameters of the neural network. Decoupling the trainable parameters from the layers of the neural network enables us to investigate and adapt the accuracy of the network propagation separated from the optimization learning problem. We numerically show that the spline-based neural network increases robustness of the learning problem towards hyperparameters due to increased stability and accuracy of the network propagation. Further, training on B-spline coefficients rather than layer weights directly enables a reduction in the number of trainable parameters.


title: Automatic evaluation of human oocyte developmental potential from microscopy images

abstract: Infertility is becoming an issue for an increasing number of couples. The most common solution, in vitro fertilization, requires embryologists to carefully examine light microscopy images of human oocytes to determine their developmental potential. We propose an automatic system to improve the speed, repeatability, and accuracy of this process. We first localize individual oocytes and identify their principal components using CNN (U-Net) segmentation. Next, we calculate several descriptors based on geometry and texture. The final step is an SVM classifier. Both the segmentation and classification training is based on expert annotations. The presented approach leads to a classification accuracy of 70%.


title: The Poisson equation involving surface measures

abstract: We prove the (optimal) $W^{1,\infty}$-regularity of weak solutions to the equation $-\Delta u = Q \; \mathcal{H}^{n-1} \llcorner \Gamma$ in a domain $\Omega \subset \mathbb{R}^n$ with Dirichlet boundary conditions, where $\Gamma \subset \subset \Omega$ is a compact (Lipschitz) manifold and $Q \in L^\infty(\Gamma)$. We also discuss optimality and necessity of the assumptions on $Q$ and $\Gamma$. Our findings can be applied to study the regularity of solutions for several free boundary problems, in particular the biharmonic Alt-Caffarelli Problem.


title: Instrumental variables, spatial confounding and interference

abstract: Unobserved spatial confounding variables are prevalent in environmental and ecological applications where the system under study is complex and the data are often observational. Instrumental variables (IVs) are a common way to address unobserved confounding; however, the efficacy of using IVs on spatial confounding is largely unknown. This paper explores the effectiveness of IVs in this situation -- with particular attention paid to the spatial scale of the instrument. We show that, in case of spatially-dependent treatments, IVs are most effective when they vary at a finer spatial resolution than the treatment. We investigate IV performance in extensive simulations and apply the model in the example of long term trends in the air pollution and cardiovascular mortality in the United States over 1990-2010. Finally, the IV approach is also extended to the spatial interference setting, in which treatments can affect nearby responses.


title: Lepton Flavour Universality tests in electroweak penguin decays at LHCb

abstract: The coupling of the electroweak gauge bosons of the Standard Model (SM) to leptons is flavour universal. Extensions of the SM do not necessarily have this property. Rare decays of heavy flavour are suppressed in the SM and new particles may give sizeable contributions to these processes, therefore, their precise study allows for sensitive tests of lepton flavour universality. Of particular interest are rare b to sll decays that are well accessible at the LHCb experiment. Recent results from LHCb on lepton flavour universality in rare b to sll decays are discussed.


title: Checking the admissibility of odd-vertex pairings is hard

abstract: Nash-Williams proved that every graph has a well-balanced orientation. A key ingredient in his proof is admissible odd-vertex pairings. We show that for two slightly different definitions of admissible odd-vertex pairings, deciding whether a given odd-vertex pairing is admissible is co-NP-complete. This resolves a question of Frank. We also show that deciding whether a given graph has an orientation that satisfies arbitrary local arc-connectivity requirements is NP-complete.


title: Two infinite families of elliptic curves with rank greater than one

abstract: We prove, using elementary methods, that each member of the infinite families of elliptic curves given by $E_m \colon y^2=x^3 - x + m^6$ and $E_m' \colon y^2=x^3 + x - m^6$ have rank at least $2$ and 3, respectively, under mild restrictions on $m$. We also prove stronger results for $E_m$ and $E_m'$ using more technical machinery.


title: Evaluation of local stress state due to grain-boundary sliding during creep within a crystal plasticity finite element multi-scale framework

abstract: Previous studies demonstrate that grain-boundary sliding could accelerate creep rate and give rise to large internal stresses that can lead to damage development, e.g. formation of wedge cracks. The present study provides more insight into the effects of grain-boundary sliding (GBS) on the deformation behaviour of realistic polycrystalline aggregates during creep, through the development of a computational framework which combines: i) the use of interface elements for sliding at grain boundaries, and ii) special triple point (in 2D) or triple line (in 3D) elements to prevent artificial dilation at these locations in the microstructure with iii) a physically-based crystal plasticity constitutive model for time-dependent inelastic deformation of the individual grains. Experimental data at various scales is used to calibrate the framework and compare with model predictions. We pay particular consideration to effects of grain boundary sliding during creep of Type 316 stainless steel, which is used extensively in structural components of the UK fleet of Advanced Gas Cooled Nuclear Reactors (AGRs). It is found that the anisotropic deformation of the grains and the mismatch in crystallographic orientation between neighbouring grains play a significant role in determining the extent of sliding on a given boundary. Their effect on the development of stress within the grains, particularly at triple grain junctions, and the increase in axial stress along transverse boundaries are quantified. The article demonstrates that the magnitude of the stress along the facets is highly-dependent on the crystallographic orientations of the neighbouring grains and the relative amount of sliding. Boundaries, transverse to the applied load tend to carry higher normal stresses of the order of 100-180 MPa, in cases where the neighbouring grains consist of plastically-harder crystallographic orientations.


title: Efficient Frequency Doubling with Active Stabilization on Chip

abstract: Thin-film lithium niobate (TFLN) is superior for integrated nanophotonics due to its outstanding properties in nearly all aspects: strong second-order nonlinearity, fast and efficient electro-optic effects, wide transparency window, and little two photon absorption and free carrier scattering. Together, they permit highly integrated nanophotonic circuits capable of complex photonic processing by incorporating disparate elements on the same chip. Yet, there has to be a demonstration that synergizes those superior properties for system advantage. Here we demonstrate such a chip that capitalizes on TFLNs favorable ferroelectricity, high second-order nonlinearity, and strong electro-optic effects. It consists of a monolithic circuit integrating a Z-cut, quasi-phase matched microring with high quality factor and a phase modulator used in active feedback control. By Pound-Drever-Hall locking, it realizes stable frequency doubling at about 50% conversion with only milliwatt pump, marking the highest by far among all nanophotonic platforms with milliwatt pumping. Our demonstration addresses a long-outstanding challenge facing cavity-based optical processing, including frequency conversion, frequency comb generation, and all-optical switching, whose stable performance is hindered by photorefractive or thermal effects. Our results further establish TFLN as an excellent material capable of optical multitasking, as desirable to build multi-functional chip devices.


title: Alternative proof of upper bound of spanning trees in a graph

abstract: We give a proof for sharp estimate for the number of spanning trees using linear algebra and generalize this bound to multigraphs. In addition, we show that this bound is tight for complete graphs. In addition, we give estimates for number of spanning trees in specific type of graphs such as product of graphs or Cartesian product of two graphs.


title: Lanczos potential of Weyl field: interpretations and applications

abstract: An attempt is made to uncover the physical meaning and significance of the obscure Lanczos tensor field which is regarded as a potential of the Weyl field. Despite being a fundamental building block of any metric theory of gravity, the Lanczos tensor has not been paid proper attention as it deserves. By providing an elucidation on this tensor field through its derivation in some particularly chosen spacetimes, we try to find its adequate interpretation. Though the Lanczos field is traditionally introduced as a gravitational analogue of the electromagnetic 4-potential field, the performed study unearths its another feature - a relativistic analogue of the Newtonian gravitational force field. A new domain of applicability of the Lanczos tensor is introduced which corroborates this new feature of the tensor.


title: On the nature of exotic matter in Morris-Thorne wormholes

abstract: It is well known that phantom energy, which is characterized by the equation of state $p=\omega\rho$, $\omega <-1$, can support Morris-Thorne wormholes since it leads to a violation of the null energy condition. The purpose of this note is to show that the converse is also true in the following sense: for a typical shape function, the equation of state of exotic matter in the vicinity of the throat is given by $p_r=\omega\rho$, $\omega <-1$, where $p_r$ is the radial pressure. Some implications thereof are also considered.


title: Virasoro constraints in quantum singularity theories

abstract: We introduce Virasoro operators for any Landau-Ginzburg pair (W, G) where W is a non-degenerate quasi-homogeneous polynomial and G is a certain group of diagonal symmetries. We propose a conjecture that the total ancestor potential of the FJRW theory of the pair (W,G) is annihilated by these Virasoro operators. We prove the conjecture in various cases, including: (1) invertible polynomials with the maximal group, (2) some two-variable polynomials with the minimal group, (3) certain Calabi-Yau polynomials with groups. We also discuss the connections among Virasoro constraints, mirror symmetry of Landau-Ginzburg models, and Landau-Ginzburg/Calabi-Yau correspondence.


title: Lagging Heat Models in Thermodynamics and Bioheat Transfer: a Critical Review

abstract: The accuracy of the classical heat conduction model, known as Fourier's law, is highly questioned, dealing with the micro and nanosystems and biological tissues. In other words, the results obtained from the classical equations deviate from the available experimental data. It means that the continuum heat diffusion equation is insufficient and inappropriate for modeling heat transport in these cases. There are several techniques for modeling non-Fourier heat conduction. In the present paper, we place our focus on the dual-phase-lag (DPL) approach. The DPL model, as a popular modification of Fourier's law, has already been utilized in numerous situations, such as simulating ultrafast laser heating and heat conduction in carbon nanotubes. There has been a sharp increase in research on non-Fourier heat conduction in recent years. Several studies have been performed in the fields of thermoelasticity, thermodynamics, transistor modeling, and bioheat transport. This review presents the most recent non-Fourier bioheat conduction works and the related thermodynamics background. The various mathematical tools, modeling different thermal therapies, and relevant criticisms and disputes are discussed. Finally, the novel and other possible studies are also presented to provide a better overview, and the roadmap to the future research and challenges ahead is drawn up.


title: Time-Varying Coefficient Model Estimation Through Radial Basis Functions

abstract: In this paper we estimate the dynamic parameters of a time-varying coefficient model through radial kernel functions in the context of a longitudinal study. Our proposal is based on a linear combination of weighted kernel functions involving a bandwidth, centered around a given set of time points. In addition, we study different alternatives of estimation and inference including a Frequentist approach using weighted least squares along with bootstrap methods, and a Bayesian approach through both Markov chain Monte Carlo and variational methods. We compare the estimation strategies mention above with each other, and our radial kernel functions proposal with an expansion based on regression spline, by means of an extensive simulation study considering multiples scenarios in terms of sample size, number of repeated measurements, and subject-specific correlation. Our experiments show that the capabilities of our proposal based on radial kernel functions are indeed comparable with or even better than those obtained from regression splines. We illustrate our methodology by analyzing data from two AIDS clinical studies.


title: Gravitational lensing in 4-D Einstein-Gauss-Bonnet gravity in the presence of plasma

abstract: In this paper we have assumed a weak-field regime to explore the gravitational lensed photons in a 4 dimensional Einstein-Gauss-Bonnet gravity, which is very much in the limelight these days. The investigation is conducted in three distinct paradigms: uniform plasma, singular isothermal sphere and a non-singular isothermal sphere. The lensing angle associated with the distribution factor of the medium is individually derived for each case and further utilized to study the magnification of the image source brightness, selectively for the uniform plasma and singular isothermal sphere. The attained results are brought forth in contrast with the standard Schwarzschild geometry


title: Ion-Beam Modification of Metastable Gallium Oxide Polymorphs

abstract: Gallium oxide with a corundum structure ({\alpha}-Ga2O3) has recently attracted great attention in view of electronic and photonic applications due to its unique properties including a wide band gap exceeding that of the most stable beta phase (\b{eta}-Ga2O3). However, the lower thermal stability of the {\alpha}-phase at ambient conditions in comparison with the \b{eta}-phase requires careful investigation of its resistance to other external influences such as ion irradiation, ion doping, etc. In this work, the structural changes under the action of Al+ ion irradiation have been investigated for a polymorphic gallium oxide layers grown by hydride vapor phase epitaxy on c-plane sapphire and consisting predominantly of {\alpha}-phase with inclusions of {\alpha}(\k{appa})-phase. It is established by the X-ray diffraction technique that inclusions of {\alpha}(\k{appa})-phase in the irradiated layer undergo the expansion along the normal to the substrate surface, while there is no a noticeable deformation for the {\alpha}-phase. This speaks in favor of the different radiation tolerance of various Ga2O3 polymorphs, especially the higher radiation tolerance of the {\alpha}-phase. This fact should be taken into account when utilizing ion implantation to modify gallium oxide properties in terms of development of efficient doping strategies.


title: Machine Learning Techniques to Construct Patched Analog Ensembles for Data Assimilation

abstract: Using generative models from the machine learning literature to create artificial ensemble members for use within data assimilation schemes has been introduced in [Grooms QJRMS, 2020] as constructed analog ensemble optimal interpolation (cAnEnOI). Specifically, we study general and variational autoencoders for the machine learning component of this method, and combine the ideas of constructed analogs and ensemble optimal interpolation in the data assimilation piece. To extend the scalability of cAnEnOI for use in data assimilation on complex dynamical models, we propose using patching schemes to divide the global spatial domain into digestible chunks. Using patches makes training the generative models possible and has the added benefit of being able to exploit parallelism during the generative step. Testing this new algorithm on a 1D toy model, we find that larger patch sizes make it harder to train an accurate generative model (i.e. a model whose reconstruction error is small), while conversely the data assimilation performance improves at larger patch sizes. There is thus a sweet spot where the patch size is large enough to enable good data assimilation performance, but not so large that it becomes difficult to train an accurate generative model. In our tests the new patched cAnEnOI method outperforms the original (unpatched) cAnEnOI, as well as the ensemble square root filter results from [Grooms QJRMS, 2020].


title: Discontinuous yielding of pristine micro-crystals

abstract: We study the mechanical response of a dislocation-free 2D crystal under homogenous shear using a new mesoscopic approach to crystal plasticity, a Landau-type theory, accounting for the global invariance of the energy in the space of strain tensors while operating with an infinite number of equivalent energy wells. The advantage of this approach is that it eliminates arbitrariness in dealing with topological transitions involved, for instance, in nucleation and annihilation of dislocations. We use discontinuous yielding of pristine micro-crystals as a benchmark problem for the new theory and show that the nature of the catastrophic instability, which in this setting inevitably follows the standard affine response, depends not only on lattice symmetry but also on the orientation of the crystal in the loading device. The ensuing dislocation avalanche involves cooperative dislocation nucleation, resulting in the formation of complex microstructures controlled by a nontrivial self-induced coupling between different plastic mechanisms.


title: Creep-fatigue interactions in a polycrystalline structural material under typical high-temperature power plant operating conditions

abstract: A micromechanical model at the microscale within a crystal plasticity self-consistent model, is used to analyse loading histories in Type 316H stainless steel, common to structural components in high-temperature power plants. The study compares the SCM predictions on changes in mechanical behaviour and creep properties to analyses via the UK R5 structural assessment code. Plant-relevant cyclic-creep histories in Type 316H stainless steel at 550C are examined with focus on the estimation of accumulated creep strain. This aims to quantify how better understanding of creep deformation under cyclic loading can inform R5 code life assessment methods, which are strain based. The effect of cyclic plasticity on primary and secondary creep is quantified. The levels of creep strain accumulated during different dwells, following loading from different macroscopic stress states, is also evaluated and link is made to experimental observations in explaining the response. The impact of displacement- and load-controlled dwells on the evolution of the cyclic hysteresis loop is examined which revealed how different assumptions for creep strain accumulation could employ information from SCM results to reduce uncertainty. It also provides guidance in the development of interpolation frameworks for structural parameters, based on the SCM results which could be used directly in structural assessment and design approaches to estimate the accumulation of inelastic strain with a lower degree of conservatism.


title: One-Point Gradient-Free Methods for Smooth and Non-Smooth Saddle-Point Problems

abstract: In this paper, we analyze gradient-free methods with one-point feedback for stochastic saddle point problems $\min_{x}\max_{y} \varphi(x, y)$. For non-smooth and smooth cases, we present analysis in a general geometric setup with arbitrary Bregman divergence. For problems with higher-order smoothness, the analysis is carried out only in the Euclidean case. The estimates we have obtained repeat the best currently known estimates of gradient-free methods with one-point feedback for problems of imagining a convex or strongly convex function. The paper uses three main approaches to recovering the gradient through finite differences: standard with a random direction, as well as its modifications with kernels and residual feedback. We also provide experiments to compare these approaches for the matrix game.


title: Compressible fluid inside a linear oscillator

abstract: We use Feireisl-Lions theory to deduce the existence of weak solutions to a system describing the dynamics of a linear oscillator containing a Newtonian compressible fluid. The appropriate Navier-Stokes equation is considered on a domain whose movement has one degree of freedom. The equation is paired with the Newton law and we assume a no-slip boundary condition.


title: Pricing Exchange Rate Options and Quanto Caps in the Cross-Currency Random Field LIBOR Market Model

abstract: We develop an arbitrage-free random field LIBOR market model to price cross-currency derivatives. The uncertainty of the forward LIBOR rates of our cross-currency model is driven by a two time parameter random field instead of a finite dimensional Brownian motion. To demonstrate the applications of this model, we develop an approximate closed-form pricing formula for Quanto caps and cross-currency swaps. Further, we derive an exact pricing formula for an exchange rate option in the random field setting.


title: Exploiting ultrasound tongue imaging for the automatic detection of speech articulation errors

abstract: Speech sound disorders are a common communication impairment in childhood. Because speech disorders can negatively affect the lives and the development of children, clinical intervention is often recommended. To help with diagnosis and treatment, clinicians use instrumented methods such as spectrograms or ultrasound tongue imaging to analyse speech articulations. Analysis with these methods can be laborious for clinicians, therefore there is growing interest in its automation. In this paper, we investigate the contribution of ultrasound tongue imaging for the automatic detection of speech articulation errors. Our systems are trained on typically developing child speech and augmented with a database of adult speech using audio and ultrasound. Evaluation on typically developing speech indicates that pre-training on adult speech and jointly using ultrasound and audio gives the best results with an accuracy of 86.9%. To evaluate on disordered speech, we collect pronunciation scores from experienced speech and language therapists, focusing on cases of velar fronting and gliding of /r/. The scores show good inter-annotator agreement for velar fronting, but not for gliding errors. For automatic velar fronting error detection, the best results are obtained when jointly using ultrasound and audio. The best system correctly detects 86.6% of the errors identified by experienced clinicians. Out of all the segments identified as errors by the best system, 73.2% match errors identified by clinicians. Results on automatic gliding detection are harder to interpret due to poor inter-annotator agreement, but appear promising. Overall findings suggest that automatic detection of speech articulation errors has potential to be integrated into ultrasound intervention software for automatically quantifying progress during speech therapy.


title: Simulation of runaway electron generation in fusion grade tokamak and suppression by impurity injection

abstract: During disruptions in fusion-grade tokamaks like ITER, large electric fields are induced following the thermal quench (TQ) period which can generate a substantial amount of Runaway Electrons (REs) that can carry up to 10 MA current with energies as high as several tens of MeV [1-3] in current quench phase (CQ). These runaway electrons can cause significant damage to the plasma-facing-components due to their localized energy deposition. To mitigate these effects, impurity injections of high-Z atoms have been proposed [1-3]. In this paper, we use a self-consistent 0D tokamak disruption model as implemented in PREDICT code [6] which has been upgraded to take into account the effect of impurity injections on RE dynamics as suggested in [4-5]. Dominant RE generation mechanisms such as the secondary avalanche mechanism as well as primary RE-generation mechanisms namely Dreicer, hot-tail, tritium decay and Compton scattering (from {\gamma}-rays emitted from activated walls) have been taken into account. In these simulations, the effect of impurities is taken into account considering collisions of REs with free and bound electrons as well as scattering from full and partially-shielded nuclear charge. These corrections were also implemented in the relativistic test particle model to simulate RE-dynamics in momentum space. We show that the presence of impurities has a non-uniform effect on the Runaway Electron Distribution function. We also show that the combined effect of pitch-angle scattering induced by the collisions with impurity ions and synchrotron emission loss results in the faster dissipation of RE-energy distribution function [7]. The variation of different RE generation mechanisms during different phases of the disruption, mainly before and after impurity injections is reported.


title: A Resolvent Criterion Approach to Strong Decay of a Multilayered Lam\'e-Heat System

abstract: We consider a multilayer hyperbolic-parabolic PDE system which constitutes a coupling of 3D thermal - 2D elastic - 3D elastic dynamics, in which the boundary interface coupling between 3D fluid and 3D structure is realized via a 2D elastic equation. Our main result here is one of strong decay for the given multilayered - heat system. That is, the solution to this composite PDE system is stabilized asymptotically to the zero state. Our proof of strong stability takes place in the "frequency domain" and ultimately appeals to the pointwise resolvent condition introduced by Tomilov [45]. This very useful result, however, requires that the semigroup associated with our multilayered FSI system be completely non-unitary (c.n.u). Accordingly, we firstly establish that the semigroup $\{e^{\mathcal {A}t}\}_{t\geq 0}$ is indeed c.n.u., in part by invoking relatively recent results of global uniqueness for overdetermined Lam\'e systems on nonsmooth domains. Although the entire proof also requires higher regularity results for some trace terms, this \textit{"resolvent criterion approach"} allows us to establish a "classially soft" proof of strong decay. In particular, it avoids the sort of technical PDE multipliers invoked in [9].


title: Bounded Exhaustive Search of Alloy Specification Repairs

abstract: The rising popularity of declarative languages and the hard to debug nature thereof have motivated the need for applicable, automated repair techniques for such languages. However, despite significant advances in the program repair of imperative languages, there is a dearth of repair techniques for declarative languages. This paper presents BeAFix, an automated repair technique for faulty models written in Alloy, a declarative language based on first-order relational logic. BeAFix is backed with a novel strategy for bounded exhaustive, yet scalable, exploration of the spaces of fix candidates and a formally rigorous, sound pruning of such spaces. Moreover, different from the state-of-the-art in Alloy automated repair, that relies on the availability of unit tests, BeAFix does not require tests and can work with assertions that are naturally used in formal declarative languages. Our experience with using BeAFix to repair thousands of realworld faulty models, collected by other researchers, corroborates its ability to effectively generate correct repairs and outperform the state-of-the-art.


title: Terahertz-Band Joint Ultra-Massive MIMO Radar-Communications: Model-Based and Model-Free Hybrid Beamforming

abstract: Wireless communications and sensing at terahertz (THz) band are increasingly investigated as promising short-range technologies because of the availability of high operational bandwidth at THz. In order to address the extremely high attenuation at THz, ultra-massive multiple-input multiple-output (MIMO) antenna systems have been proposed for THz communications to compensate propagation losses. However, the cost and power associated with fully digital beamformers of these huge antenna arrays are prohibitive. In this paper, we develop wideband hybrid beamformers based on both model-based and model-free techniques for a new group-of-subarrays (GoSA) ultra-massive MIMO structure in low-THz band. Further, driven by the recent developments to save the spectrum, we propose beamformers for a joint ultra-massive MIMO radar-communications system, wherein the base station serves multi-antenna user equipment (RX), and tracks radar targets by generating multiple beams toward both RX and the targets. We formulate the GoSA beamformer design as an optimization problem to provide a trade-off between the unconstrained communications beamformers and the desired radar beamformers. To mitigate the beam split effect at THz band arising from frequency-independent analog beamformers, we propose a phase correction technique to align the beams of multiple subcarriers toward a single physical direction. To further decrease the ultra-massive MIMO computational complexity and enhance robustness, we also implement deep learning solutions to the proposed model-based hybrid beamformers. Numerical experiments demonstrate that both techniques outperform the conventional approaches in terms of spectral efficiency and radar beampatterns, as well as exhibiting less hardware cost and computation time.


title: Optimal control of point-to-point navigation in turbulent time-dependent flows using Reinforcement Learning

abstract: We present theoretical and numerical results concerning the problem to find the path that minimizes the time to navigate between two given points in a complex fluid under realistic navigation constraints. We contrast deterministic Optimal Navigation (ON) control with stochastic policies obtained by Reinforcement Learning (RL) algorithms. We show that Actor-Critic RL algorithms are able to find quasi-optimal solutions in the presence of either time-independent or chaotically evolving flow configurations. For our application, ON solutions develop unstable behavior within the typical duration of the navigation process, and are therefore not useful in practice. We first explore navigation of turbulent flow using a constant propulsion speed. Based on a discretized phase-space, the propulsion direction is adjusted with the aim to minimize the time spent to reach the target. Further, we explore a case where additional control is obtained by allowing the engine to power off. Exploiting advection of the underlying flow, allows the target to be reached with less energy consumption. In this case, we optimize a linear combination between the total navigation time and the total time the engine is switched off. Our approach can be generalized to other setups, for example, navigation under imperfect environmental forecast or with different models for the moving vessel.


title: Thermodynamics of Riemannian Kerr-AdS black holes in Poincar\'e gauge theory

abstract: A Hamiltonian approach to black hole entropy is used to study Riemannian Kerr-AdS solutions in the general, parity-violating Poincar\'e gauge theory. Entropy and the asymptotic charges are entirely determined by the parity-even sector of the theory, whereas the parity-odd contributions vanish. Entropy is found to be proportional to the horizon area, and the first law of black hole thermodynamics is confirmed.


title: CP-MDP: A CANDECOMP-PARAFAC Decomposition Approach to Solve a Markov Decision Process Multidimensional Problem

abstract: Markov Decision Process (MDP) is the underlying model for optimal planning for decision-theoretic agents in stochastic environments. Although much research focuses on solving MDP problems both in tabular form or using factored representations, none focused on tensor decomposition methods. Solving MDPs using tensor algebra offers the prospect of leveraging advances in tensor-based computations to further increase solver efficiency. In this paper, we develop an MDP solver for a multidimensional problem using a tensor decomposition method to compress the transition models and optimize the value iteration and policy iteration algorithms. We empirically evaluate our approach against tabular methods and show our approach can compute much larger problems using substantially less memory, opening up new possibilities for tensor-based approaches in stochastic planning


title: Hubble Space Telescope observations of the old pulsar PSR J0108-1431

abstract: We present results of optical-UV observations of the 200 Myr old rotation-powered radio pulsar J0108$-$1431 with the Hubble Space Telescope. We found a putative candidate for the far-UV (FUV) pulsar counterpart, with the flux density $f_\nu = 9.0\pm 3.2$ nJy at $\lambda = 1528$ \AA. The pulsar was not detected, however, at longer wavelengths, with $3\sigma$ upper limits of 52, 37, and 87 nJy at $\lambda =$ 4326, 3355, and 2366 \AA, respectively. Assuming that the pulsar counterpart was indeed detected in FUV, and the previously reported marginal $U$ and $B$ detections with the Very Large Telescope were real, the optical-UV spectrum of the pulsar can be described by a power-law model with a nearly flat $f_\nu$ spectrum. Similar to younger pulsars detected in the optical, the slope of the nonthermal spectrum steepens in the X-ray range. The pulsar's luminosity in the 1500--6000 \AA wavelength range, $L \sim 1.2\times 10^{27} (d/210\,{\rm pc})^2$ erg s$^{-1}$, corresponds to a high efficiency of conversion of pulsar rotation energy loss rate $\dot {E}$ to the optical-UV radiation, $\eta = L/\dot{E} \sim (1$--$6)\times 10^{-4}$, depending on somewhat uncertain values of distance and spectral slope. The brightness temperature of the bulk neutron star surface does not exceed 59,000 K ($3\sigma$ upper bound), as seen by a distant observer. If we assume that the FUV flux is dominated by a thermal component, then the surface temperature can be in the range of 27,000--55,000 K, requiring a heating mechanism to operate in old neutron stars.


title: Silent versus modal multi-speaker speech recognition from ultrasound and video

abstract: We investigate multi-speaker speech recognition from ultrasound images of the tongue and video images of the lips. We train our systems on imaging data from modal speech, and evaluate on matched test sets of two speaking modes: silent and modal speech. We observe that silent speech recognition from imaging data underperforms compared to modal speech recognition, likely due to a speaking-mode mismatch between training and testing. We improve silent speech recognition performance using techniques that address the domain mismatch, such as fMLLR and unsupervised model adaptation. We also analyse the properties of silent and modal speech in terms of utterance duration and the size of the articulatory space. To estimate the articulatory space, we compute the convex hull of tongue splines, extracted from ultrasound tongue images. Overall, we observe that the duration of silent speech is longer than that of modal speech, and that silent speech covers a smaller articulatory space than modal speech. Although these two properties are statistically significant across speaking modes, they do not directly correlate with word error rates from speech recognition.


title: BiconNet: An Edge-preserved Connectivity-based Approach for Salient Object Detection

abstract: Salient object detection (SOD) is viewed as a pixel-wise saliency modeling task by traditional deep learning-based methods. A limitation of current SOD models is insufficient utilization of inter-pixel information, which usually results in imperfect segmentation near edge regions and low spatial coherence. As we demonstrate, using a saliency mask as the only label is suboptimal. To address this limitation, we propose a connectivity-based approach called bilateral connectivity network (BiconNet), which uses connectivity masks together with saliency masks as labels for effective modeling of inter-pixel relationships and object saliency. Moreover, we propose a bilateral voting module to enhance the output connectivity map, and a novel edge feature enhancement method that efficiently utilizes edge-specific features. Through comprehensive experiments on five benchmark datasets, we demonstrate that our proposed method can be plugged into any existing state-of-the-art saliency-based SOD framework to improve its performance with negligible parameter increase.


title: Expectation of the Site Frequency Spectrum

abstract: The site frequency spectrum describes variation among a set of n DNA sequences. Its i'th entry (i=1,2,...,n-1) is the number of nucleotide sites at which the mutant allele is present in i copies. Under selective neutrality, random mating, and constant population size, the expected value of the spectrum is well known but somewhat puzzling. Each additional sequence added to a sample adds an entry to the end of the expected spectrum but does not affect existing entries. This note reviews the reasons for this behavior.


title: Transformers with Competitive Ensembles of Independent Mechanisms

abstract: An important development in deep learning from the earliest MLPs has been a move towards architectures with structural inductive biases which enable the model to keep distinct sources of information and routes of processing well-separated. This structure is linked to the notion of independent mechanisms from the causality literature, in which a mechanism is able to retain the same processing as irrelevant aspects of the world are changed. For example, convnets enable separation over positions, while attention-based architectures (especially Transformers) learn which combination of positions to process dynamically. In this work we explore a way in which the Transformer architecture is deficient: it represents each position with a large monolithic hidden representation and a single set of parameters which are applied over the entire hidden representation. This potentially throws unrelated sources of information together, and limits the Transformer's ability to capture independent mechanisms. To address this, we propose Transformers with Independent Mechanisms (TIM), a new Transformer layer which divides the hidden representation and parameters into multiple mechanisms, which only exchange information through attention. Additionally, we propose a competition mechanism which encourages these mechanisms to specialize over time steps, and thus be more independent. We study TIM on a large-scale BERT model, on the Image Transformer, and on speech enhancement and find evidence for semantically meaningful specialization as well as improved performance.


title: Extreme Superposition: High-Order Fundamental Rogue Waves in the Far-Field Regime

abstract: We study fundamental rogue-wave solutions of the focusing nonlinear Schr\"odinger equation in the limit that the order of the rogue wave is large and the independent variables $(x,t)$ are proportional to the order (the far-field limit). We first formulate a Riemann-Hilbert representation of these solutions that allows the order to vary continuously rather than by integer increments. The intermediate solutions in this continuous family include also soliton solutions for zero boundary conditions spectrally encoded by a single complex-conjugate pair of poles of arbitrary order, as well as other solutions having nonzero boundary conditions matching those of the rogue waves albeit with far slower decay as $x\to\pm\infty$. The large-order far-field asymptotic behavior of the solution depends on which of three disjoint regions $\mathcal{C}$, $\mathcal{S}$, and $\mathcal{E}$ contains the rescaled variables. On the regions $\mathcal{C}$ and $\mathcal{S}$ we show that the asymptotic behavior is the same for all continuous orders, while in the region $\mathcal{E}$ the discrete sequence of rogue-wave orders produces distinctive asymptotic behavior that is different from other cases.


title: Quasi-stationary distribution for the Langevin process in cylindrical domains, part II: overdamped limit

abstract: Consider the Langevin process, described by a vector (positions and momenta) in $\mathbb{R}^{d}\times\mathbb{R}^d$. Let $\mathcal O$ be a $\mathcal{C}^2$ open bounded and connected set of $\mathbb{R}^d$. Recent works showed the existence of a unique quasi-stationary distribution (QSD) of the Langevin process on the domain $D:=\mathcal{O}\times\mathbb{R}^d$. In this article, we study the overdamped limit of this QSD, i.e. when the friction coefficient goes to infinity. In particular, we show that the marginal law in position of the overdamped limit is the QSD of the overdamped Langevin process on the domain $\mathcal{O}$.


title: Covalent 2D Cr$_2$Te$_3$ ferromagnet

abstract: To broaden the scope of van der Waals 2D magnets, we report the synthesis and magnetism of covalent 2D magnetic Cr$_2$Te$_3$ with a thickness down to one-unit-cell. The 2D Cr$_2$Te$_3$ crystals exhibit robust ferromagnetism with a Curie temperature of 180 K, a large perpendicular anisotropy of 7*105 J m-3, and a high coercivity of ~ 4.6 kG at 20 K. First-principles calculations further show a transition from canted to collinear ferromagnetism, a transition from perpendicular to in-plane anisotropy, and emergent half-metallic behavior in atomically-thin Cr$_2$Te$_3$, suggesting its potential application for injecting carriers with high spin polarization into spintronic devices.


title: Nonlocal doubly nonlinear diffusion problems with nonlinear boundary conditions

abstract: We study the existence and uniqueness of mild and strong solutions of nonlocal nonlinear diffusion problems of $p$-Laplacian type with nonlinear boundary conditions posed in metric random walk spaces. These spaces include, among others, weighted discrete graphs and $\mathbb{R}^N$ with a random walk induced by a nonsingular kernel. We also study the case of nonlinear dynamical boundary conditions. The generality of the nonlinearities considered allow us to cover the nonlocal counterparts of a large scope of local diffusion problems: Stefan problems, Hele-Shaw problems, diffusion in porous media problems, obstacle problems, and more. Nonlinear semigroup theory is the basis for this study.


title: $L(2,1)$-Labeling of the iterated Mycielski of graphs and some related to matching problems

abstract: In this paper, we study the $L(2, 1)$-Labeling of the Mycielski and the iterated Mycielski of graphs in general. For a graph $G$ and all $t\geq 1$, we give sharp bounds for $\lambda(M^t(G))$ the $L(2, 1)$-labeling number of the $t$-th iterated Mycielski in terms of the number of iterations $t$, the order $n$, the maximum degree $\bigtriangleup$, and $\lambda(G)$ the $L(2, 1)$-labeling number of $G$. For $t=1$, we present necessary and sufficient conditions between the $4$-star matching number of the complement graph and $\lambda(M(G))$ the $L(2, 1)$-labeling number of the Mycielski of a graph, with some applications to special graphs. For all $t\geq 2$, we prove that for any graph $G$ of order $n$, we have $2^{t-1}(n+2)-2\leq \lambda(M^t(G))\leq 2^{t}(n+1)-2$. Thereafter, we characterize the graphs achieving the upper bound $2^t(n+1)-2$, then by using the Marriage Theorem and Tutte's characterization of graphs with a perfect $2$-matching, we characterize all graphs without isolated vertices achieving the lower bound $2^{t-1}(n+2)-2$. We determine the $L(2, 1)$-labeling number for the Mycielski and the iterated Mycielski of some graph classes.


title: Constrained Differentially Private Federated Learning for Low-bandwidth Devices

abstract: Federated learning becomes a prominent approach when different entities want to learn collaboratively a common model without sharing their training data. However, Federated learning has two main drawbacks. First, it is quite bandwidth inefficient as it involves a lot of message exchanges between the aggregating server and the participating entities. This bandwidth and corresponding processing costs could be prohibitive if the participating entities are, for example, mobile devices. Furthermore, although federated learning improves privacy by not sharing data, recent attacks have shown that it still leaks information about the training data. This paper presents a novel privacy-preserving federated learning scheme. The proposed scheme provides theoretical privacy guarantees, as it is based on Differential Privacy. Furthermore, it optimizes the model accuracy by constraining the model learning phase on few selected weights. Finally, as shown experimentally, it reduces the upstream and downstream bandwidth by up to 99.9% compared to standard federated learning, making it practical for mobile systems.


title: Bar pattern speeds in CALIFA galaxies III. Solving the puzzle of ultrafast bars

abstract: More than 10% of the barred galaxies with a direct measurement of the bar pattern speed host an ultrafast bar. These bars extend beyond the corotation radius and challenge our understanding of the orbital structure of barred galaxies. Most of them are found in spiral galaxies, rather than in lenticular ones. We analysed the properties of the ultrafast bars detected in the CALIFA Survey to investigate whether they are an artefact resulting from an overestimation of the bar radius and/or an underestimation of the corotation radius or a new class of bars, whose orbital structure has not yet been understood. We revised the available measurements of the bar radius based on ellipse fitting and Fourier analysis and of the bar pattern speed from the Tremaine-Weinberg method. In addition, we measured the bar radius from the analysis of the maps tracing the transverse-to-radial force ratio, which we obtained from the deprojected i-band images of the galaxies retrieved from the SDSS Survey. We found that nearly all the sample galaxies are spirals with an inner ring or pseudo-ring circling the bar and/or strong spiral arms, which hamper the measurement of the bar radius from the ellipse fitting and Fourier analysis. According to these methods, the bar ends overlap the ring or the spiral arms making the adopted bar radius unreliable. On the contrary, the bar radius from the ratio maps are shorter than the corotation radius. This is in agreement with the theoretical predictions and findings of numerical simulations about the extension and stability of the stellar orbits supporting the bars. We conclude that ultrafast bars are no longer observed when the correct measurement of the bar radius is adopted. Deriving the bar radius in galaxies with rings and strong spiral arms is not straightforward and a solid measurement method based on both photometric and kinematic data is still missing.


title: Extra gauge bosons and lepton flavor universality violation in $\Upsilon$ and $B$ meson decays

abstract: Lepton flavor universality can be tested through the ratio of semileptonic $B$ meson decays and leptonic $\Upsilon$ meson decays, with $\Upsilon \equiv \Upsilon(nS)$ ($n=1,2,3$). For the charged-current transitions $b \to c\tau\bar{\nu}_\tau$, discrepancies between the experiment and the Standard Model (SM) have been observed in recent years by different flavor facilities such as BABAR, Belle, and LHCb. While for the neutral-current transitions $b \bar{b} \to \tau\bar{\tau}$, the BABAR experiment reported recently a new measurement of leptonic decay ratio $R_{\Upsilon(3S)} = {\rm BR}(\Upsilon(3S) \to \tau^+\tau^-)/{\rm BR}(\Upsilon(3S) \to \mu^+\mu^-)$, showing an agreement with the SM at the $1.8 \sigma$ level. In light of this new BABAR result and regarding the connection between new physics (NP) interpretations to the charged-current $b \to c \tau \bar{\nu}_{\tau}$ anomalies and neutral-current $b \bar{b} \to \tau \bar{\tau}$ processes, in this study, we revisit the NP consequences of this measurement within a simplified model with extra massive gauge bosons that coupled predominantly to left-handed leptons of the third-generation. We show that the BABAR measurement of $R_{\Upsilon(3S)}$ cannot easily be accommodated (within its experimental $1\sigma$ range) together with the other $b \to c\tau\bar{\nu}_\tau$ data, hinting toward a new anomalous observable.


title: End-to-end Uncertainty-based Mitigation of Adversarial Attacks to Automated Lane Centering

abstract: In the development of advanced driver-assistance systems (ADAS) and autonomous vehicles, machine learning techniques that are based on deep neural networks (DNNs) have been widely used for vehicle perception. These techniques offer significant improvement on average perception accuracy over traditional methods, however, have been shown to be susceptible to adversarial attacks, where small perturbations in the input may cause significant errors in the perception results and lead to system failure. Most prior works addressing such adversarial attacks focus only on the sensing and perception modules. In this work, we propose an end-to-end approach that addresses the impact of adversarial attacks throughout perception, planning, and control modules. In particular, we choose a target ADAS application, the automated lane centering system in OpenPilot, quantify the perception uncertainty under adversarial attacks, and design a robust planning and control module accordingly based on the uncertainty analysis. We evaluate our proposed approach using both the public dataset and production-grade autonomous driving simulator. The experiment results demonstrate that our approach can effectively mitigate the impact of adversarial attacks and can achieve 55% to 90% improvement over the original OpenPilot.


title: Luttinger sum rules and spin fractionalization in the SU(N) Kondo Lattice

abstract: We show how Oshikawa's theorem for the Fermi surface volume of the Kondo lattice can be extended to the SU$(N)$ symmetric case. By extending the theorem, we are able to show that the mechanism of Fermi surface expansion seen in the large $N$ mean-field theory is directly linked to the expansion of the Fermi surface in a spin-$1/2$ Kondo lattice. This linkage enables us to interpret the expansion of the Fermi surface in a Kondo lattice as a fractionalization of the local moments into heavy electrons. Our method allows extension to a pure U(1) spin liquid, where we find the volume of the spinon Fermi surface by applying a spin-twist, analogous to Oshikawa's flux insertion. Lastly, we discuss the possibility of interpreting the FL$^*$ phase characterised by a small Fermi surface in the absence of symmetry breaking, as a non-topological coexistence of such a U(1) spin liquid and an electronic Fermi liquid.


title: Better Together? How Externalities of Size Complicate Notions of Solidarity and Actuarial Fairness

abstract: Consider a cost-sharing game with players of different contribution to the total cost: an example might be an insurance company calculating premiums for a population of mixed-risk individuals. Two natural and competing notions of fairness might be to a) charge each individual the same price or b) charge each individual according to the cost that they bring to the pool. In the insurance literature, these general approaches are referred to as "solidarity" and "actuarial fairness" and are commonly viewed as opposites. However, in insurance (and many other natural settings), the cost-sharing game also exhibits "externalities of size": all else being equal, larger groups have lower average cost. In the insurance case, we analyze a model with externalities of size due to a reduction in the variability of losses. We explore how this complicates traditional understandings of fairness, drawing on literature in cooperative game theory. First, we explore solidarity: we show that it is possible for both groups (high and low risk) to strictly benefit by joining an insurance pool where costs are evenly split, as opposed to being in separate risk pools. We build on this by producing a pricing scheme that maximally subsidizes the high risk group, while maintaining an incentive for lower risk people to stay in the insurance pool. Next, we demonstrate that with this new model, the price charged to each individual has to depend on the risk of other participants, making naive actuarial fairness inefficient. Furthermore, we prove that stable pricing schemes must be ones where players have the anti-social incentive of desiring riskier partners, contradicting motivations for using actuarial fairness. Finally, we describe how these results relate to debates about fairness in machine learning and potential avenues for future research.


title: The variational approach of an elliptic problem and its solution by finite elements

abstract: In this paper, we deal with an elliptic problem with the Dirichlet boundary condition. We operate in Sobolev spaces and the main analytic tool we use is the Lax-Milgram lemma. First, we present the variational approach of the problem which allows us to apply different functional analysis techniques. Then we study thoroughly the well-posedness of the problem. We conclude our work with a solution of the problem using numerical analysis techniques and the free software freefem++.


title: High-Dimensional Bayesian Optimization with Sparse Axis-Aligned Subspaces

abstract: Bayesian optimization (BO) is a powerful paradigm for efficient optimization of black-box objective functions. High-dimensional BO presents a particular challenge, in part because the curse of dimensionality makes it difficult to define -- as well as do inference over -- a suitable class of surrogate models. We argue that Gaussian process surrogate models defined on sparse axis-aligned subspaces offer an attractive compromise between flexibility and parsimony. We demonstrate that our approach, which relies on Hamiltonian Monte Carlo for inference, can rapidly identify sparse subspaces relevant to modeling the unknown objective function, enabling sample-efficient high-dimensional BO. In an extensive suite of experiments comparing to existing methods for high-dimensional BO we demonstrate that our algorithm, Sparse Axis-Aligned Subspace BO (SAASBO), achieves excellent performance on several synthetic and real-world problems without the need to set problem-specific hyperparameters.


title: Proton irradiation effects on metal-YBCO interfaces

abstract: 10 MeV proton-irradiation effects on a YBCO-based test structure were analyzed by measuring its current-voltage (IV) characteristics for different cumulated fluences. For fluences of up to $\sim$80$\cdot$10$^9$~p/cm$^2$ no changes in the electrical behavior of the device were observed, while for a fluence of $\sim$~300$\cdot$10$^9~$ p/cm$^2$ it becomes less conducting. A detailed analysis of the room temperature IV characteristics based on the $\gamma$ power exponent parameter [$\gamma=dLn(I)/dLn(V)$] allowed us to reveal the main conduction mechanisms as well as to establish the equivalent circuit model of the device. The changes produced in the electrical behavior, in accordance with Monte Carlo TRIM simulations, suggest that the main effect induced by protons is the displacement of oxygen atoms within the YBCO lattice, particularly from oxygen-rich to oxygen-poor areas, where they become trapped.