title: gambit -- An Open Source Name Disambiguation Tool for Version Control Systems

abstract: Name disambiguation is a complex but highly relevant challenge whenever analysing real-world user data, such as data from version control systems. We propose gambit, a rule-based disambiguation tool that only relies on name and email information. We evaluate its performance against two commonly used algorithms with similar characteristics on manually disambiguated ground-truth data from the Gnome GTK project. Our results show that gambit significantly outperforms both algorithms, achieving an F1 score of 0.985.


title: Higher rank quantum-classical correspondence

abstract: For a compact Riemannian locally symmetric space $\Gamma\backslash G/K$ of arbitrary rank we determine the location of certain Ruelle-Taylor resonances for the Weyl chamber action. We provide a Weyl-lower bound on an appropriate counting function for the Ruelle-Taylor resonances and establish a spectral gap which is uniform in $\Gamma$ if $G/K$ is irreducible of higher rank. This is achieved by proving a quantum-classical correspondence, i.e. a 1:1-correspondence between horocyclically invariant Ruelle-Taylor resonant states and joint eigenfunctions of the algebra of invariant differential operators on $G/K$.


title: Core Challenges of Social Robot Navigation: A Survey

abstract: Robot navigation in crowded public spaces is a complex task that requires addressing a variety of engineering and human factors challenges. These challenges have motivated a great amount of research resulting in important developments for the fields of robotics and human-robot interaction over the past three decades. Despite the significant progress and the massive recent interest, we observe a number of significant remaining challenges that prohibit the seamless deployment of autonomous robots in public pedestrian environments. In this survey article, we organize existing challenges into a set of categories related to broader open problems in motion planning, behavior design, and evaluation methodologies. Within these categories, we review past work, and offer directions for future research. Our work builds upon and extends earlier survey efforts by a) taking a critical perspective and diagnosing fundamental limitations of adopted practices in the field and b) offering constructive feedback and ideas that we aspire will drive research in the field over the coming decade.


title: Self-Supervision by Prediction for Object Discovery in Videos

abstract: Despite their irresistible success, deep learning algorithms still heavily rely on annotated data. On the other hand, unsupervised settings pose many challenges, especially about determining the right inductive bias in diverse scenarios. One scalable solution is to make the model generate the supervision for itself by leveraging some part of the input data, which is known as self-supervised learning. In this paper, we use the prediction task as self-supervision and build a novel object-centric model for image sequence representation. In addition to disentangling the notion of objects and the motion dynamics, our compositional structure explicitly handles occlusion and inpaints inferred objects and background for the composition of the predicted frame. With the aid of auxiliary loss functions that promote spatially and temporally consistent object representations, our self-supervised framework can be trained without the help of any manual annotation or pretrained network. Initial experiments confirm that the proposed pipeline is a promising step towards object-centric video prediction.


title: Strong coupling effects in quantum thermal transport with the reaction coordinate method

abstract: We present a semi-analytical approach for studying quantum thermal energy transport beyond the weak system-bath coupling regime. Our treatment, which results in a renormalized, effective Hamiltonian model is based on the reaction coordinate method. In our technique, applied to the nonequilibrium spin-boson model, a collective coordinate is extracted from each environment and added into the system to construct an enlarged system. After performing additional Hamiltonian's truncation and transformation, we attain an effective two-level system with renormalized parameters, which is weakly coupled to its environments, thus can be simulated using a perturbative Markovian quantum master equation approach. We compare the heat current characteristics in our method to other techniques, and demonstrate that we properly capture strong system-bath signatures such as the turnover behavior of the heat current as a function of system-bath coupling strength. We further investigate the thermal diode effect and demonstrate that strong couplings moderately improve the rectification ratio relative to the weak coupling limit. The effective Hamiltonian method that we developed here offers fundamental insight into the strong coupling behavior, and is computationally economic. Applications of the method towards studying quantum thermal machines are anticipated.


title: Noninteracting fermionic systems with localized losses: Exact results in the hydrodynamic limit

abstract: We investigate the interplay between unitary dynamics after a quantum quench and localized dissipation in a noninteracting fermionic chain. In particular, we consider the effect of gain and loss processes, for which fermions are added and removed incoherently. We focus on the hydrodynamic limit of large distances from the source of dissipation and of long times, with their ratio being fixed. In this limit, the localized dissipation gives rise to an effective delta potential (dissipative impurity), and the time-evolution of the local correlation functions admits a simple hydrodynamic description in terms of the fermionic occupations in the initial state and the reflection and transmission amplitudes of the impurity. We derive this hydrodynamic framework from the ab initio calculation of the microscopic dynamics. This allows us to analytically characterize the effect of dissipation for several theoretically relevant initial states, such as a uniform Fermi sea, homogeneous product states, or the inhomogeneous state obtained by joining two Fermi seas. In this latter setting, when both gain and loss processes are present, we observe the emergence of exotic nonequilibrium steady states with stepwise uniform density profiles. In all instances, for strong dissipation the coherent dynamics of the system is arrested, which is a manifestation of the celebrated quantum Zeno effect.

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


title: Entropy-Guided Control Improvisation

abstract: High level declarative constraints provide a powerful (and popular) way to define and construct control policies; however, most synthesis algorithms do not support specifying the degree of randomness (unpredictability) of the resulting controller. In many contexts, e.g., patrolling, testing, behavior prediction,and planning on idealized models, predictable or biased controllers are undesirable. To address these concerns, we introduce the \emph{Entropic Reactive Control Improvisation} (ERCI) framework and algorithm which supports synthesizing control policies for stochastic games that are declaratively specified by (i) a \emph{hard constraint} specifying what must occur, (ii) a \emph{soft constraint} specifying what typically occurs, and (iii) a \emph{randomization constraint} specifying the unpredictability and variety of the controller, as quantified using causal entropy. This framework, extends the state of the art by supporting arbitrary combinations of adversarial and probabilistic uncertainty in the environment. ERCI enables a flexible modeling formalism which we argue, theoretically and empirically, remains tractable.


title: u-cf2vec: Representation Learning for Personalized Algorithm Selection in Recommender Systems

abstract: Collaborative Filtering (CF) has become the standard approach to solve recommendation systems (RS) problems. Collaborative Filtering algorithms try to make predictions about interests of a user by collecting the personal interests from multiple users. There are multiple CF algorithms, each one of them with its own biases. It is the Machine Learning practitioner that has to choose the best algorithm for each task beforehand. In Recommender Systems, different algorithms have different performance for different users within the same dataset. Meta Learning (MtL) has been used to choose the best algorithm for a given problem. Meta Learning is usually applied to select algorithms for a whole dataset. Adapting it to select the to the algorithm for a single user in a RS involves several challenges. The most important is the design of the metafeatures which, in typical meta learning, characterize datasets while here, they must characterize a single user. This work presents a new meta-learning based framework named $\mu$-cf2vec to select the best algorithm for each user. We propose using Representation Learning techniques to extract the metafeatures. Representation Learning tries to extract representations that can be reused in other learning tasks. In this work we also implement the framework using different RL techniques to evaluate which one can be more useful to solve this task. In the meta level, the meta learning model will use the metafeatures to extract knowledge that will be used to predict the best algorithm for each user. We evaluated an implementation of this framework using MovieLens 20M dataset. Our implementation achieved consistent gains in the meta level, however, in the base level we only achieved marginal gains.


title: Synthesizing Computable Functions from Rational Specifications over Infinite Words

abstract: The synthesis problem asks to automatically generate, if it exists, an algorithm from a specification of correct input-output pairs. In this paper, we consider the synthesis of computable functions of infinite words, for a classical Turing computability notion over infinite inputs. We consider specifications which are rational relations of infinite words, i.e., specifications defined non-deterministic parity transducers. We prove that the synthesis problem of computable functions from rational specifications is undecidable. We provide an incomplete but sound reduction to some parity game, such that if Eve wins the game, then the rational specification is realizable by a computable function. We prove that this function is even computable by a deterministic two-way transducer. We provide a sufficient condition under which the latter game reduction is complete. This entails the decidability of the synthesis problem of computable functions, which we prove to be ExpTime-complete, for a large subclass of rational specifications, namely deterministic rational specifications. This subclass contains the class of automatic relations over infinite words, a yardstick in reactive synthesis.


title: Photometric variability as a proxy for magnetic activity and its dependence on metallicity

abstract: Understanding how the magnetic activity of low-mass stars depends on their fundamental parameters is an important goal of stellar astrophysics. Previous studies show that activity levels are largely determined by the stellar Rossby number which is defined as the rotation period divided by the convective turnover time. However, we currently have little information on the role that chemical composition plays. In this work, we investigate how metallicity affects magnetic activity using photometric variability as an activity proxy. Similarly to other proxies, we demonstrate that the amplitude of photometric variability is well parameterised by the Rossby number, although in a more complex way. We also show that variability amplitude and metallicity are generally positively correlated. This trend can be understood in terms of the effect that metallicity has on stellar structure and, hence, the convective turnover time (or, equivalently, the Rossby number). Lastly, we demonstrate that the metallicity dependence of photometric variability results in a rotation period detection bias whereby the periods of metal-rich stars are more easily recovered for stars of a given mass.


title: Formulating Intuitive Stack-of-Tasks using Visuo-Tactile Perception for Collaborative Human-Robot Fine Manipulation

abstract: Enabling robots to work in close proximity to humans necessitates a control framework that does not only incorporate multi-sensory information for autonomous and coordinated interactions but also has perceptive task planning to ensure an adaptable and flexible collaborative behaviour. In this research, an intuitive stack-of-tasks (iSoT) formulation is proposed, that defines the robot's actions by considering the human-arm postures and the task progression. The framework is augmented with visuo-tactile information to effectively perceive the collaborative environment and intuitively switch between the planned sub-tasks. The visual feedback from depth cameras monitors and estimates the objects' poses and human-arm postures, while the tactile data provides the exploration skills to detect and maintain the desired contacts to avoid object slippage. To evaluate the performance, effectiveness and usability of the proposed framework, assembly and disassembly tasks, performed by the human-human and human-robot partners, are considered and analyzed using distinct evaluation metrics i.e, approach adaptation, grasp correction, task coordination latency, cumulative posture deviation, and task repeatability.


title: SMIL: Multimodal Learning with Severely Missing Modality

abstract: A common assumption in multimodal learning is the completeness of training data, i.e., full modalities are available in all training examples. Although there exists research endeavor in developing novel methods to tackle the incompleteness of testing data, e.g., modalities are partially missing in testing examples, few of them can handle incomplete training modalities. The problem becomes even more challenging if considering the case of severely missing, e.g., 90% training examples may have incomplete modalities. For the first time in the literature, this paper formally studies multimodal learning with missing modality in terms of flexibility (missing modalities in training, testing, or both) and efficiency (most training data have incomplete modality). Technically, we propose a new method named SMIL that leverages Bayesian meta-learning in uniformly achieving both objectives. To validate our idea, we conduct a series of experiments on three popular benchmarks: MM-IMDb, CMU-MOSI, and avMNIST. The results prove the state-of-the-art performance of SMIL over existing methods and generative baselines including autoencoders and generative adversarial networks. Our code is available at https://github.com/mengmenm/SMIL.


title: Explaining dimensionality reduction results using Shapley values

abstract: Dimensionality reduction (DR) techniques have been consistently supporting high-dimensional data analysis in various applications. Besides the patterns uncovered by these techniques, the interpretation of DR results based on each feature's contribution to the low-dimensional representation supports new finds through exploratory analysis. Current literature approaches designed to interpret DR techniques do not explain the features' contributions well since they focus only on the low-dimensional representation or do not consider the relationship among features. This paper presents ClusterShapley to address these problems, using Shapley values to generate explanations of dimensionality reduction techniques and interpret these algorithms using a cluster-oriented analysis. ClusterShapley explains the formation of clusters and the meaning of their relationship, which is useful for exploratory data analysis in various domains. We propose novel visualization techniques to guide the interpretation of features' contributions on clustering formation and validate our methodology through case studies of publicly available datasets. The results demonstrate our approach's interpretability and analysis power to generate insights about pathologies and patients in different conditions using DR results.


title: Quarkyonic Effective Field Theory, Quark-Nucleon Duality and Ghosts

abstract: We present a field theoretical description of quarkyonic matter consisting of quark, nucleon and ghost fields coupling to mesonic degrees of freedom. The ghosts are present to cancel over-counting of nucleon states that are Pauli blocked by the quark Fermi sea. Such a theory becomes an effective field theory of nucleons at low baryon density, and as such will reproduce nucleonic matter phenomenology. This theory can accommodate chiral symmetry restoration and the dynamical generation of a shell of nucleons at the Fermi surface. It is valid for finite temperature and density. In such a theory, quark-nucleon duality is accomplished by inclusion of ghost fields so that the nucleons extra degrees of freedom, that are beyond those of quarks, are compensated by the ghost fields.


title: Spectral properties of three-dimensional Anderson model

abstract: The three-dimensional Anderson model represents a paradigmatic model to understand the Anderson localization transition. In this work we first review some key results obtained for this model in the past 50 years, and then study its properties from the perspective of modern numerical approaches. Our main focus is on the quantitative comparison between the level sensitivity statistics and the level statistics. While the former studies the sensitivity of Hamiltonian eigenlevels upon inserting a magnetic flux, the latter studies the properties of unperturbed eigenlevels. We define two versions of dimensionless conductance, the first corresponding to the width of the level curvature distribution relative to the mean level spacing, and the second corresponding to the ratio of the Heisenberg and the Thouless time obtained from the spectral form factor. We show that both conductances look remarkably similar around the localization transition, in particular, they predict a nearly identical critical point consistent with other measures of the transition. We then study some further properties of those quantities: for level curvatures, we discuss particular similarities and differences between the width of the level curvature distribution and the characteristic energy studied by Edwards and Thouless in their pioneering work [J. Phys. C. 5, 807 (1972)]. In the context of the spectral form factor, we show that at the critical point it enters a broad time-independent regime, in which its value is consistent with the level compressibility obtained from the level variance. Finally, we test the scaling solution of the average level spacing ratio in the crossover regime using the cost function minimization approach introduced in [Phys. Rev. B. 102, 064207 (2020)]. We find that the extracted transition point and the scaling coefficient agree with those from the literature to high numerical accuracy.


title: Resource-Aware Stochastic Self-Triggered Model Predictive Control

abstract: This paper considers the control of uncertain systems that are operated under limited resource factors, such as battery life or hardware longevity. We consider here resource-aware self-triggered control techniques that schedule system operation non-uniformly in time in order to balance performance against resource consumption. When running in an uncertain environment, unknown disturbances may deteriorate system performance by acting adversarially against the planned event triggering schedule. In this work, we propose a resource-aware stochastic predictive control scheme to tackle this challenge, where a novel zero-order hold feedback control scheme is proposed to accommodate a time-inhomogeneous predictive control update.


title: Towards Action Model Learning for Player Modeling

abstract: Player modeling attempts to create a computational model which accurately approximates a player's behavior in a game. Most player modeling techniques rely on domain knowledge and are not transferable across games. Additionally, player models do not currently yield any explanatory insight about a player's cognitive processes, such as the creation and refinement of mental models. In this paper, we present our findings with using action model learning (AML), in which an action model is learned given data in the form of a play trace, to learn a player model in a domain-agnostic manner. We demonstrate the utility of this model by introducing a technique to quantitatively estimate how well a player understands the mechanics of a game. We evaluate an existing AML algorithm (FAMA) for player modeling and develop a novel algorithm called Blackout that is inspired by player cognition. We compare Blackout with FAMA using the puzzle game Sokoban and show that Blackout generates better player models.


title: Combining Context-Free and Contextualized Representations for Arabic Sarcasm Detection and Sentiment Identification

abstract: Since their inception, transformer-based language models have led to impressive performance gains across multiple natural language processing tasks. For Arabic, the current state-of-the-art results on most datasets are achieved by the AraBERT language model. Notwithstanding these recent advancements, sarcasm and sentiment detection persist to be challenging tasks in Arabic, given the language's rich morphology, linguistic disparity and dialectal variations. This paper proffers team SPPU-AASM's submission for the WANLP ArSarcasm shared-task 2021, which centers around the sarcasm and sentiment polarity detection of Arabic tweets. The study proposes a hybrid model, combining sentence representations from AraBERT with static word vectors trained on Arabic social media corpora. The proposed system achieves a F1-sarcastic score of 0.62 and a F-PN score of 0.715 for the sarcasm and sentiment detection tasks, respectively. Simulation results show that the proposed system outperforms multiple existing approaches for both the tasks, suggesting that the amalgamation of context-free and context-dependent text representations can help capture complementary facets of word meaning in Arabic. The system ranked second and tenth in the respective sub-tasks of sarcasm detection and sentiment identification.


title: Monotonic Alpha-divergence Minimisation for Variational Inference

abstract: In this paper, we introduce a novel family of iterative algorithms which carry out $\alpha$-divergence minimisation in a Variational Inference context. They do so by ensuring a systematic decrease at each step in the $\alpha$-divergence between the variational and the posterior distributions. In its most general form, the variational distribution is a mixture model and our framework allows us to simultaneously optimise the weights and components parameters of this mixture model. Our approach permits us to build on various methods previously proposed for $\alpha$-divergence minimisation such as Gradient or Power Descent schemes and we also shed a new light on an integrated Expectation Maximization algorithm. Lastly, we provide empirical evidence that our methodology yields improved results on several multimodal target distributions and on a real data example.


title: Elastic Forces Drive Nonequilibrium Pattern Formation in a Model of Nanocrystal Ion Exchange

abstract: Chemical transformations, such as ion exchange, are commonly employed to modify nanocrystal compositions. Yet the mechanisms of these transformations, which often operate far from equilibrium and entail mixing diverse chemical species, remain poorly understood. Here, we explore an idealized model for ion exchange in which a chemical potential drives compositional defects to accumulate at a crystal's surface. These impurities subsequently diffuse inward. We find that the nature of interactions between sites in a compositionally impure crystal strongly impacts exchange trajectories. In particular, elastic deformations which accompany lattice-mismatched species promote spatially modulated patterns in the composition. These same patterns can be produced at equilibrium in core/shell nanocrystals, whose structure mimics transient motifs observed in nonequilibrium trajectories. Moreover, the core of such nanocrystals undergoes a phase transition - from modulated to unstructured - as the thickness or stiffness of the shell is decreased. Our results help explain the varied patterns observed in heterostructured nanocrystals produced by ion exchange and suggest principles for the rational design of compositionally-patterned nanomaterials.


title: Quantum interference and imaging using intense laser fields

abstract: The interference of matter waves is one of the intriguing features of quantum mechanics that has impressed researchers and laymen since it was first suggested almost a century ago. Nowadays, attosecond science tools allow us to utilize it in order to extract valuable information from electron wavepackets. Intense laser fields are routinely employed to create electron wave packets and control their motion with sub-femtosecond and sub-nanometer precision. In this perspective article, we discuss some of the peculiarities of intense light-matter interaction. We review some of the most important techniques used in attosecond imaging, namely photoelectron holography and laser-induced electron diffraction. We attempt to ask and answer a few questions that do not get asked very often. For example, if we are interested in position space information, why are measurements carried out in momentum space? How to accurately retrieve photoelecron spactra from the numerical solution of the time-dependent Schr\"odinger equation? And, what causes the different coherence properties of high-harmonic generation and above-threshold ionization?


title: Capturing Omni-Range Context for Omnidirectional Segmentation

abstract: Convolutional Networks (ConvNets) excel at semantic segmentation and have become a vital component for perception in autonomous driving. Enabling an all-encompassing view of street-scenes, omnidirectional cameras present themselves as a perfect fit in such systems. Most segmentation models for parsing urban environments operate on common, narrow Field of View (FoV) images. Transferring these models from the domain they were designed for to 360-degree perception, their performance drops dramatically, e.g., by an absolute 30.0% (mIoU) on established test-beds. To bridge the gap in terms of FoV and structural distribution between the imaging domains, we introduce Efficient Concurrent Attention Networks (ECANets), directly capturing the inherent long-range dependencies in omnidirectional imagery. In addition to the learned attention-based contextual priors that can stretch across 360-degree images, we upgrade model training by leveraging multi-source and omni-supervised learning, taking advantage of both: Densely labeled and unlabeled data originating from multiple datasets. To foster progress in panoramic image segmentation, we put forward and extensively evaluate models on Wild PAnoramic Semantic Segmentation (WildPASS), a dataset designed to capture diverse scenes from all around the globe. Our novel model, training regimen and multi-source prediction fusion elevate the performance (mIoU) to new state-of-the-art results on the public PASS (60.2%) and the fresh WildPASS (69.0%) benchmarks.


title: Electroweak phase transition triggered by fermion sector

abstract: To realize first-order electroweak phase transition, it is necessary to generate a barrier in the thermal Higgs potential, which is usually triggered by scalar degree of freedom. We instead investigate phase transition patterns in pure fermion extensions of the standard model, and find that additional fermions with mass hierarchy and mixing could develop such barrier and realize strongly first-order phase transition in such models. In the Higgs potential with polynomial parametrization, the barrier can be generated in the following two patterns: (I) positive quadratic term, negative cubic term and positive quartic term or (II) positive quadratic term, negative quartic term and positive higher dimensional term, such as dimensional 6 operator.


title: Design Principles for Data Analysis

abstract: The data science revolution has led to an increased interest in the practice of data analysis. While much has been written about statistical thinking, a complementary form of thinking that appears in the practice of data analysis is design thinking -- the problem-solving process to understand the people for whom a product is being designed. For a given problem, there can be significant or subtle differences in how a data analyst (or producer of a data analysis) constructs, creates, or designs a data analysis, including differences in the choice of methods, tooling, and workflow. These choices can affect the data analysis products themselves and the experience of the consumer of the data analysis. Therefore, the role of a producer can be thought of as designing the data analysis with a set of design principles. Here, we introduce design principles for data analysis and describe how they can be mapped to data analyses in a quantitative, objective and informative manner. We also provide empirical evidence of variation of principles within and between both producers and consumers of data analyses. Our work leads to two insights: it suggests a formal mechanism to describe data analyses based on the design principles for data analysis, and it provides a framework to teach students how to build data analyses using formal design principles.


title: Multitask 3D CBCT-to-CT Translation and Organs-at-Risk Segmentation Using Physics-Based Data Augmentation

abstract: In current clinical practice, noisy and artifact-ridden weekly cone-beam computed tomography (CBCT) images are only used for patient setup during radiotherapy. Treatment planning is done once at the beginning of the treatment using high-quality planning CT (pCT) images and manual contours for organs-at-risk (OARs) structures. If the quality of the weekly CBCT images can be improved while simultaneously segmenting OAR structures, this can provide critical information for adapting radiotherapy mid-treatment as well as for deriving biomarkers for treatment response. Using a novel physics-based data augmentation strategy, we synthesize a large dataset of perfectly/inherently registered planning CT and synthetic-CBCT pairs for locally advanced lung cancer patient cohort, which are then used in a multitask 3D deep learning framework to simultaneously segment and translate real weekly CBCT images to high-quality planning CT-like images. We compared the synthetic CT and OAR segmentations generated by the model to real planning CT and manual OAR segmentations and showed promising results. The real week 1 (baseline) CBCT images which had an average MAE of 162.77 HU compared to pCT images are translated to synthetic CT images that exhibit a drastically improved average MAE of 29.31 HU and average structural similarity of 92% with the pCT images. The average DICE scores of the 3D organs-at-risk segmentations are: lungs 0.96, heart 0.88, spinal cord 0.83 and esophagus 0.66. This approach could allow clinicians to adjust treatment plans using only the routine low-quality CBCT images, potentially improving patient outcomes. Our code, data, and pre-trained models will be made available via our physics-based data augmentation library, Physics-ArX, at https://github.com/nadeemlab/Physics-ArX.


title: Specific loss power of magnetic nanoparticles (fluid) hyperthermia in non-adiabatic conditions

abstract: We investigate the magnetic nanoparticles (fluid) hyperthermia in non-adiabatic conditions through the calorimetric method. Specifically, we propose a theoretical approach to magnetic hyperthermia from a thermodynamic point of view. To test the robustness of the approach, we perform hyperthermia experiments and analyze the thermal behavior of magnetite and magnesium ferrite magnetic nanoparticles dispersed in water submitted to an alternating magnetic field. From our findings, besides estimating the specific loss power value from a non-adiabatic process, thus enhancing the accuracy in the determination of this quantity, we provide physical meaning to parameters found in literature that still remained not fully understood, and bring to light how they can be obtained experimentally.


title: Quantifying Sufficient Randomness for Causal Inference

abstract: Spurious association arises from covariance between propensity for the treatment and individual risk for the outcome. For sensitivity analysis with stochastic counterfactuals we introduce a methodology to characterize uncertainty in causal inference from natural experiments and quasi-experiments. Our sensitivity parameters are standardized measures of variation in propensity and individual risk, and one minus their geometric mean is an intuitive measure of randomness in the data generating process. Within our latent propensity-risk model, we show how to compute from contingency table data a threshold, $T$, of sufficient randomness for causal inference. If the actual randomness of the data generating process exceeds this threshold then causal inference is warranted.


title: Non-invasive analysis of blood-brain barrier permeability based on wavelet and machine learning approaches

abstract: The blood-brain barrier plays a decisive role in protecting the brain from toxins and pathogens. The ability to analyze the BBB opening (OBBB) is crucial for the treatment of many brain diseases, but it is very difficult to noninvasively monitor OBBB. In this paper we analyse the EEG series of healthy rats in free behaviour and after music-induced OBBB. The research is performed using two completely different methods based on wavelet analysis and machine learning approach. Both methods enable us to recognize OBBB and are in a good agreement with each other. The comparative analysis was carried out using F-measures and ROC-curves.


title: Exploiting Asynchronous Priority Scheduling in Parallel Eikonal Solvers

abstract: Numerical solutions to the Eikonal equation are computed using variants of the fast marching method, the fast sweeping method, and the fast iterative method. In this paper, we provide a unified view of these algorithms that highlights their similarities and suggests a wider class of Eikonal solvers. We then use this framework to justify applying concurrent priority scheduling techniques to Eikonal solvers. We demonstrate that doing so results in good parallel performance for a problem from seismology. We explain why existing Eikonal solvers may produce different results despite using the same differencing scheme and demonstrate techniques to address these discrepancies. These techniques allow us to obtain deterministic output from our asynchronous fine-grained parallel Eikonal solver.


title: Carrier Tuning of Stoner Ferromagnetism in ThCr$_{\mathbf{2}}$Si$_{\mathbf{2}}$-Structure Cobalt Arsenides

abstract: CaCo$_{2-y}$As$_2$ is an unusual itinerant magnet with signatures of extreme magnetic frustration. The conditions for establishing magnetic order in such itinerant frustrated magnets, either by reducing frustration or increasing electronic correlations, is an open question. Here we use results from inelastic neutron scattering and magnetic susceptibility measurements and density functional theory calculations to show that hole doping in Ca(Co$_{1-x}$Fe$_{x}$)$_{2-y}$As$_{2}$ suppresses magnetic order by quenching the magnetic moment while maintaining the same level of magnetic frustration. The suppression is due to tuning the Fermi energy away from a peak in the electronic density of states originating from a flat conduction band. This results in the complete elimination of the magnetic moment by $x\approx0.25$, providing a clear example of a Stoner-type transition.


title: Chebyshev Polynomials and Inequalities for Kleinian Groups

abstract: The principal character of a representation of the free group of rank two into PSL(2, C) is a triple of complex numbers that determines an irreducible representation uniquely up to conjugacy. It is a central problem in the geometry of discrete groups and low dimensional topology to determine when such a triple represents a discrete group that is not virtually abelian, that is a Kleinian group. A classical necessary condition is J{\o}rgensen's inequality. Here we use certainly shifted Chebyshev polynomials and trace identities to determine new families of such inequalities, some of which are best possible. The use of these polynomials also shows how we can identify the principal character of some important subgroups from that of the group itself.


title: On Security Properties of All-or-nothing Transforms

abstract: All-or-nothing transforms have been defined as bijective mappings on all s-tuples over a specified finite alphabet. These mappings are required to satisfy certain "perfect security" conditions specified using entropies of the probability distribution defined on the input s-tuples. Alternatively, purely combinatorial definitions of AONTs have been given, which involve certain kinds of "unbiased arrays". However, the combinatorial definition makes no reference to probability definitions. In this paper, we examine the security provided by AONTs that satisfy the combinatorial definition. The security of the AONT can depend on the underlying probability distribution of the s-tuples. We show that perfect security is obtained from an AONT if and only if the input s-tuples are equiprobable. However, in the case where the input s-tuples are not equiprobable, we still achieve a weaker security guarantee. We also consider the use of randomized AONTs to provide perfect security for a smaller number of inputs, even when those inputs are not equiprobable.


title: Gemini/GMOS Optical Transmission Spectroscopy of WASP-121b: signs of variability in an ultra-hot Jupiter?

abstract: We present ground-based, spectroscopic observations of two transits of the ultra-hot Jupiter WASP-121b covering the wavelength range $\approx$500 - 950 nm using Gemini/GMOS. We use a Gaussian process framework to model instrumental systematics in the light curves, and also demonstrate the use of the more generalised Student's-T process to verify our results. We find that our measured transmission spectrum, whilst showing overall agreement, is slightly discrepant with results obtained using HST/STIS, particularly for wavelengths shortward of $\approx$650 nm. In contrast to the STIS results, we find evidence for an increasing blueward slope and little evidence for absorption from either TiO or VO in our retrieval, in agreement with a number of recent studies performed at high-resolution. We suggest that this might point to some other absorbers, particularly some combination of recently detected atomic metals, in addition to scattering by hazes, being responsible for the excess optical absorption and observed vertical thermal inversion. Our results are also broadly consistent with previous ground-based photometry and 3D GCM predictions, however, these assumed different chemistry to our retrievals. In addition, we show that the GMOS observations are repeatable over short periods (days), similarly to the HST/STIS observations. Their difference over longer periods (months) could well be the result of temporal variability in the atmospheric properties (i.e. weather) as predicted by theoretical models of ultra-hot Jupiters; however, more mundane explanations such as instrumental systematics and stellar activity cannot be fully ruled out, and we encourage future observations to explore this possibility.


title: Foreword: Manhattan Project Nuclear Science & Technology Developments at Los Alamos

abstract: A collection of articles reviewing Trinity on its 75th anniversary.


title: Fast dissipation of Colliding Alfv\'en Waves in a Magnetically Dominated Plasma

abstract: Magnetic energy around compact objects often dominates over plasma rest mass, and its dissipation can power the object luminosity. We describe a dissipation mechanism which works faster than magnetic reconnection. The mechanism involves two strong Alfv\'en waves with anti-aligned magnetic fields $\boldsymbol{B}_1$ and $\boldsymbol{B}_2$ that propagate in opposite directions along background magnetic field $\boldsymbol{B}_0$ and collide. The collision forms a thin current sheet perpendicular to $\boldsymbol{B}_0$, which absorbs the incoming waves. The current sheet is sustained by electric field $\boldsymbol{E}$ breaking the magnetohydrodynamic condition $E<B$ and accelerating particles to high energies. We demonstrate this mechanism with kinetic plasma simulations using a simple setup of two symmetric plane waves with amplitude $A=B_1/B_0=B_2/B_0$ propagating in a uniform $\boldsymbol{B}_0$. The mechanism is activated when $A>1/2$. It dissipates a large fraction of the wave energy, $f=(2A-1)/A^2$, reaching $100\%$ when $A=1$. The plane geometry allows one to see the dissipation process in a one-dimensional simulation. We also perform two-dimensional simulations, enabling spontaneous breaking of the plane symmetry by the tearing instability of the current sheet. At moderate $A$ of main interest the tearing instability is suppressed. Dissipation transitions to normal, slower, magnetic reconnection at $A\gg 1$. The fast dissipation described in this paper may occur in various objects with perturbed magnetic fields, including magnetars, jets from accreting black holes, and pulsar wind nebulae.


title: Arbitrary order total variation Convergence of Markov semigroups using random grids

abstract: We provide a abstract framework to prove total variation convergence result with arbitrary rate for numerical scheme for SDE. In particular we show that under standard weak approximation properties of scheme such as Euler we can obtain total variation convergence with any desired rate by building a specific approximation for the SDE.


title: Automorphisms and some geodesic properties of ortho-Grassmann graphs

abstract: Let $H$ be a complex Hilbert space. Consider the ortho-Grassmann graph $\Gamma^{\perp}_{k}(H)$ whose vertices are $k$-dimensional subspaces of $H$ (projections of rank $k$) and two subspaces are connected by an edge in this graph if they are compatible and adjacent (the corresponding rank-$k$ projections commute and their difference is an operator of rank $2$). Our main result is the following: if $\dim H\ne 2k$, then every automorphism of $\Gamma^{\perp}_{k}(H)$ is induced by a unitary or anti-unitary operator; if $\dim H=2k\ge 6$, then every automorphism of $\Gamma^{\perp}_{k}(H)$ is induced by a unitary or anti-unitary operator or it is the composition of such an automorphism and the orthocomplementary map. For the case when $\dim H=2k=4$ the statement fails. To prove this statement we compare geodesics of length two in ortho-Grassmann graphs and characterise compatibility (commutativity) in terms of geodesics in Grassmann and ortho-Grassmann graphs. At the end, we extend this result on generalised ortho-Grassmann graphs associated to conjugacy classes of finite-rank self-adjoint operators.


title: Knots are Generic Stable Phases in Semiflexible Polymers

abstract: Semiflexible polymer models are widely used as a paradigm to understand structural phases in biomolecules including folding of proteins. Since stable knots are not so common in real proteins, the existence of stable knots in semiflexible polymers has not been explored much. Here, via extensive replica exchange Monte Carlo simulation we investigate the same for a bead-stick and a bead-spring homopolymer model that covers the whole range from flexible to stiff. We establish the fact that the presence of stable knotted phases in the phase diagram is dependent on the ratio $r_b/r_{\rm{min}}$ where $r_b$ is the equilibrium bond length and $r_{\rm{min}}$ is the distance for the strongest nonbonded contacts. Our results provide evidence for both models that if the ratio $r_b/r_{\rm{min}}$ is outside a small window around unity then depending on the bending stiffness one always encounters stable knotted phases along with the usual frozen and bent-like structures at low temperatures. These findings prompt us to conclude that knots are generic stable phases in semiflexible polymers.


title: Automatic code generation from sketches of mobile applications in end-user development using Deep Learning

abstract: A common need for mobile application development by end-users or in computing education is to transform a sketch of a user interface into wireframe code using App Inventor, a popular block-based programming environment. As this task is challenging and time-consuming, we present the Sketch2aia approach that automates this process. Sketch2aia employs deep learning to detect the most frequent user interface components and their position on a hand-drawn sketch creating an intermediate representation of the user interface and then automatically generates the App Inventor code of the wireframe. The approach achieves an average user interface component classification accuracy of 87,72% and results of a preliminary user evaluation indicate that it generates wireframes that closely mirror the sketches in terms of visual similarity. The approach has been implemented as a web tool and can be used to support the end-user development of mobile applications effectively and efficiently as well as the teaching of user interface design in K-12.


title: The Los Alamos Computing Facility during the Manhattan Project

abstract: This article describes the history of the computing facility at Los Alamos during the Manhattan Project, 1944 to 1946. The hand computations are briefly discussed, but the focus is on the IBM Punch Card Accounting Machines (PCAM). During WWII the Los Alamos facility was one of most advanced PCAM facilities both in the machines and in the problems being solved.


title: A sampling criterion for constrained Bayesian optimization with uncertainties

abstract: We consider the problem of chance constrained optimization where it is sought to optimize a function and satisfy constraints, both of which are affected by uncertainties. The real world declinations of this problem are particularly challenging because of their inherent computational cost. To tackle such problems, we propose a new Bayesian optimization method. It applies to the situation where the uncertainty comes from some of the inputs, so that it becomes possible to define an acquisition criterion in the joint controlled-uncontrolled input space. The main contribution of this work is an acquisition criterion that accounts for both the average improvement in objective function and the constraint reliability. The criterion is derived following the Stepwise Uncertainty Reduction logic and its maximization provides both optimal controlled and uncontrolled parameters. Analytical expressions are given to efficiently calculate the criterion. Numerical studies on test functions are presented. It is found through experimental comparisons with alternative sampling criteria that the adequation between the sampling criterion and the problem contributes to the efficiency of the overall optimization. As a side result, an expression for the variance of the improvement is given.


title: Endurance-Aware Mapping of Spiking Neural Networks to Neuromorphic Hardware

abstract: Neuromorphic computing systems are embracing memristors to implement high density and low power synaptic storage as crossbar arrays in hardware. These systems are energy efficient in executing Spiking Neural Networks (SNNs). We observe that long bitlines and wordlines in a memristive crossbar are a major source of parasitic voltage drops, which create current asymmetry. Through circuit simulations, we show the significant endurance variation that results from this asymmetry. Therefore, if the critical memristors (ones with lower endurance) are overutilized, they may lead to a reduction of the crossbar's lifetime. We propose eSpine, a novel technique to improve lifetime by incorporating the endurance variation within each crossbar in mapping machine learning workloads, ensuring that synapses with higher activation are always implemented on memristors with higher endurance, and vice versa. eSpine works in two steps. First, it uses the Kernighan-Lin Graph Partitioning algorithm to partition a workload into clusters of neurons and synapses, where each cluster can fit in a crossbar. Second, it uses an instance of Particle Swarm Optimization (PSO) to map clusters to tiles, where the placement of synapses of a cluster to memristors of a crossbar is performed by analyzing their activation within the workload. We evaluate eSpine for a state-of-the-art neuromorphic hardware model with phase-change memory (PCM)-based memristors. Using 10 SNN workloads, we demonstrate a significant improvement in the effective lifetime.


title: Machine Learning the period finding algorithm

abstract: We use differentiable programming and gradient descent to find unitary matrices that can be used in the period finding algorithm to extract period information from the state of a quantum computer post application of the oracle. The standard procedure is to use the inverse quantum Fourier transform. Our findings suggest that that this is not the only unitary matrix appropriate for the period finding algorithm, There exist several unitary matrices that can affect out the same transformation and they are significantly different from each other as well. These unitary matrices can be learned by an algorithm. Neural networks can be applied to differentiate such unitary matrices from randomly generated ones indicating that these unitaries do have characteristic features that cannot otherwise be discerned easily.


title: Orbital Order, Superconductivity, Pseudogap and Spectral Weight in High-Tc Cuprates

abstract: After providing a brief genealogy of our recently proposed model for High-Tc cuprates, we investigate the details of the microscopic mechanism that produces an attractive interaction between neighboring holes. We show that a peculiar arrangement of the $p_x$ and $p_y$ oxygen orbitals makes the mutual magnetic interaction of the holes with the localized copper ions to produce a net attractive interaction between themselves, which is responsible for the emergence of a superconducting phase. We also study the connection existing between the proposed pseudogap order parameter and the spectral density. We show that the occurrence of two sharp peaks in the latter, between which the density of states suffers a depletion is a direct consequence of the d-wave character of the pseudogap order parameter dependence on $\mathbf{k}$, which breaks the 90$^\circ$-rotation symmetry of the oxygen lattices. The peak separation in the spectral density works effectively as an overall pseudogap order parameter for the cuprates. We explicitly calculate the spectral density in the strange metal and pseudogap phases of Bi2212, at different temperatures, and show that our results compare very well with the experimental data.


title: On the Symmetry of Blast Waves

abstract: This article presents a brief historical review of G. I. Taylor's solution of the point blast wave problem which was applied to the Trinity test of the first atomic bomb. Lie group symmetry techniques are used to derive Taylor's famous two-fifths law that relates the position of a blast wave to the time after the explosion and the total energy released. The theory of exterior differential systems is combined with the method of characteristics to demonstrate that the solution of the blast wave problem is directly related to the basic relationships that exist between the geometry (or symmetry) and the physics of wave propagation through the equations of motion. The point blast wave model is cast in terms of two exterior differential systems and both systems are shown to be integrable with local solutions for the velocity, pressure, and density along curves in space and time behind the blast wave.


title: Incremental Relaying for Power Line Communication: Performance Analysis and Power Allocation

abstract: In this paper, incremental decode-and-forward (IDF) and incremental selective decode-and-forward (ISDF) relaying are proposed to improve the spectral efficiency of power line communication. Contrary to the traditional decode-and-forward (DF) relaying, IDF and ISDF strategies utilize the relay only if the direct link ceases to attain a certain information rate, thereby improving the spectral efficiency. The path gain through the power line is assumed to be log-normally distributed with high distance-dependent attenuation and the additive noise is from a Bernoulli-Gaussian process. Closed-form expressions for the outage probability, and approximate closed-form expressions for the end-to-end average channel capacity and the average bit error rate for binary phase-shift keying are derived. Furthermore, a closed-form expression for the fraction of times the relay is in use is derived as a measure of the spectral efficiency. Comparative analysis of IDF and ISDF with traditional DF relaying is presented. It is shown that IDF is a specific case of ISDF and can obtain optimal spectral efficiency without compromising the outage performance. By employing power allocation to minimize the outage probability, it is realized that the power should be allocated in accordance with the inter-node distances and channel parameters.


title: Simple Flagellated Soft Robot for Locomotion near Air-Fluid Interface

abstract: A wide range of microorganisms, e.g. bacteria, propel themselves by rotation of soft helical tails, also known as flagella. Due to the small size of these organisms, viscous forces overwhelm inertial effects and the flow is at low Reynolds number. In this fluid-structure problem, a competition between elastic forces and hydrodynamic (viscous) forces leads to a net propulsive force forward. A thorough understanding of this highly coupled fluid-structure interaction problem can not only help us better understand biological propulsion but also help us design bio-inspired functional robots with applications in oil spill cleanup, water quality monitoring, and infrastructure inspection. Here, we introduce arguably the simplest soft robot with a single binary control signal, which is capable of moving along an arbitrary 2D trajectory near air-fluid interface and at the interface between two fluids. The robot exploits the variation in viscosity to move along the prescribed trajectory. Our analysis of this newly introduced soft robot consists of three main components. First, we fabricate this simple robot and use it as an experimental testbed. Second, a discrete differential geometry-based modeling framework is used for simulation of the robot. Upon validation of the simulation tool, the third part of this study employs the simulations to develop a control scheme with a single binary input to make the robot follow any prescribed path.


title: Bayesian constraints on the origin and geology of exo-planetary material using a population of externally polluted white dwarfs

abstract: White dwarfs that have accreted planetary bodies are a powerful probe of the bulk composition of exoplanetary material. In this paper, we present a Bayesian model to explain the abundances observed in the atmospheres of 202 DZ white dwarfs by considering the heating, geochemical differentiation, and collisional processes experienced by the planetary bodies accreted, as well as gravitational sinking. The majority (>60%) of systems are consistent with the accretion of primitive material. We attribute the small spread in refractory abundances observed to a similar spread in the initial planet-forming material, as seen in the compositions of nearby stars. A range in Na abundances in the pollutant material is attributed to a range in formation temperatures from below 1,000K to higher than 1,400K, suggesting that pollutant material arrives in white dwarf atmospheres from a variety of radial locations. We also find that Solar System-like differentiation is common place in exo-planetary systems. Extreme siderophile (Fe, Ni or Cr) abundances in 8 systems require the accretion of a core-rich fragment of a larger differentiated body to at least a 3sigma significance, whilst one system shows evidence that it accreted a crust-rich fragment. In systems where the abundances suggest that accretion has finished (13/202), the total mass accreted can be calculated. The 13 systems are estimated to have accreted masses ranging from the mass of the Moon to half that of Vesta. Our analysis suggests that accretion continues for 11Myrs on average.


title: The Trinity High Explosive Implosion System: The Foundation for Precision Explosive Applications

abstract: This article is set during the 1944 and 1945 final push to complete Project Y -- the Manhattan Project at Los Alamos -- and focuses primarily on overcoming the challenge of creating and demonstrating a successful convergent explosive implosion to turn a subcritical quantity of plutonium into a critical mass. The critical mass would then efficiently yield kilotons of trinitrotoluene (TNT)-equivalent energy in about a microsecond, demonstrating the implosion atomic bomb concept. This work culminated in the Trinity atomic test near Alamogordo on July 16, 1945. This implosion effect demarcated the approach to explosive science and technology the Laboratory has followed ever since, including development of high-explosive synthesis and formulation, small and large test and diagnostic facilities, shock dynamics theory, high-explosive system design engineering, and three-dimensional implosion modeling and simulation using some of the fastest computers in the world. This work also ushered in a period of broader application of precision high explosives in conventional munitions, demolition, mining and oil exploration, and space travel.


title: Pixel-wise Distance Regression for Glacier Calving Front Detection and Segmentation

abstract: Glacier calving front position (CFP) is an important glaciological variable. Traditionally, delineating the CFPs has been carried out manually, which was subjective, tedious and expensive. Automating this process is crucial for continuously monitoring the evolution and status of glaciers. Recently, deep learning approaches have been investigated for this application. However, the current methods get challenged by a severe class-imbalance problem. In this work, we propose to mitigate the class-imbalance between the calving front class and the non-calving front class by reformulating the segmentation problem into a pixel-wise regression task. A Convolutional Neural Network gets optimized to predict the distance values to the glacier front for each pixel in the image. The resulting distance map localizes the CFP and is further post-processed to extract the calving front line. We propose three post-processing methods, one method based on statistical thresholding, a second method based on conditional random fields (CRF), and finally the use of a second U-Net. The experimental results confirm that our approach significantly outperforms the state-of-the-art methods and produces accurate delineation. The Second U-Net obtains the best performance results, resulting in an average improvement of about 21% dice coefficient enhancement.


title: Error-and-Erasure Decoding of Product and Staircase Codes

abstract: High-rate product codes (PCs) and staircase codes (SCs) are ubiquitous codes in high-speed optical communication achieving near-capacity performance on the binary symmetric channel. Their success is mostly due to very efficient iterative decoding algorithms that require very little complexity. In this paper, we extend the density evolution (DE) analysis for PCs and SCs to a channel with ternary output and ternary message passing, where the third symbol marks an erasure. We investigate the performance of a standard error-and-erasure decoder and of its simplification using DE. The proposed analysis can be used to find component code configurations and quantizer levels for the channel output. We also show how the use of even-weight BCH subcodes as component codes can improve the decoding performance at high rates. The DE results are verified by Monte-Carlo simulations, which show that additional coding gains of up to 0.6 dB are possible by ternary decoding, at only a small additional increase in complexity compared to traditional binary message passing.


title: Secrecy Outage Probability of Cognitive Small-Cell Network with Unreliable Backhaul Connections

abstract: In this paper, we investigate the secrecy performance of underlay cognitive small-cell radio network with unreliable backhaul connections. The secondary cognitive small-cell transmitters are connected to macro base station by wireless backhaul links. The small-cell network is sharing the same spectrum with the primary network ensuring that a desired outage probability constraint in the primary network is always satisfied. We propose an optimal transmitter selection (OTS) scheme for small-cell network to transfer information to the destination. The closed-form expression of secrecy outage probability are derived. Our result shows that increasing the primary transmitter's transmit power and the number of small-cell transmitter can improve the system performance. The backhaul reliability of secondary and the desired outage probability of the primary also have significant impact on the system.


title: Data driven reconstruction using frames and Riesz bases

abstract: We study the problem of regularization of inverse problems adopting a purely data driven approach, by using the similarity to the method of regularization by projection. We provide an application of a projection algorithm, utilized and applied in frames theory, as a data driven reconstruction procedure in inverse problems, generalizing the algorithm proposed by the authors in Inverse Problems 36 (2020), n. 12, 125009, based on an orthonormalization procedure for the training pairs. We show some numerical experiments, comparing the different methods.


title: Spheroidal Ambisonics: a Spatial Audio Framework Using Spheroidal Bases

abstract: Ambisonics is an established framework to capture, process, and reproduce spatial sound fields based on its spherical harmonics representation. We propose a generalization of conventional spherical ambisonics to the spheroidal coordinate system and spheroidal microphone arrays, which represent sound fields by means of spheroidal wave functions. This framework is referred to as spheroidal ambisonics and a formulation for the case of prolate spheroidal coordinates is presented. Spheroidal ambisonics allows analytical encoding of sound fields using spheroidal microphone arrays. In addition, an analytical conversion formula from spheroidal ambisonics to spherical ambisonics is derived in order to ensure compatibility with the existing ecosystem of spherical ambisonics. Numerical experiments are performed to verify spheroidal ambisonic encoding and transcoding when used for spatial sound field recording. It is found that the sound field reconstructed from the transcoded coefficients has a zone of accurate reconstruction which is prolonged towards the long axis of a prolate spheroidal microphone array.


title: Power-Balanced Hybrid Optics Boosted Design for Achromatic Extended-Depth-of-Field Imaging via Optimized Mixed OTF

abstract: The power-balanced hybrid optical imaging system is a special design of a diffractive computational camera, introduced in this paper, with image formation by a refractive lens and Multilevel Phase Mask (MPM). This system provides a long focal depth with low chromatic aberrations thanks to MPM and a high energy light concentration due to the refractive lens. We introduce the concept of optical power balance between the lens and MPM which controls the contribution of each element to modulate the incoming light. Additional unique features of our MPM design are the inclusion of quantization of the MPM's shape on the number of levels and the Fresnel order (thickness) using a smoothing function. To optimize optical power-balance as well as the MPM, we build a fully-differentiable image formation model for joint optimization of optical and imaging parameters for the proposed camera using Neural Network techniques. Additionally, we optimize a single Wiener-like optical transfer function (OTF) invariant to depth to reconstruct a sharp image. We numerically and experimentally compare the designed system with its counterparts, lensless and just-lens optical systems, for the visible wavelength interval (400-700)nm and the depth-of-field range (0.5-$\infty$m for numerical and 0.5-2m for experimental). The attained results demonstrate that the proposed system equipped with the optimal OTF overcomes its counterparts (even when they are used with optimized OTF) in terms of reconstruction quality for off-focus distances. The simulation results also reveal that optimizing the optical power-balance, Fresnel order, and the number of levels parameters are essential for system performance attaining an improvement of up to 5dB of PSNR using the optimized OTF compared with its counterpart lensless setup.


title: A finite element analysis model to predict and optimize the mechanical behaviour of bioprinted scaffolds

abstract: Bioprinting is an enabling biofabrication technique to create heterogeneous tissue constructs according to patient-specific geometries and compositions. Optimization of bioinks as per requirements for specific tissue applications is a critical exercise in ensuring clinical translation of the bioprinting technologies. Most notably, optimum hydrogel polymer concentrations are required to ensure adequate mechanical properties of bioprinted constructs without causing significant shear stresses on cells. However, experimental iterations are often tedious for optimizing the bioink properties. In this work, a finite element modelling approach has been undertaken to determine the effect of different bioink parameters like composition, concentration on the range of stresses being experienced by the cells in a bioprinting process. The stress distribution of the cells at different parts of the constructs has also been modelled. It is found that both bioink chemical compositions and stoichiometric concentrations can substantially alter the stress effects experienced by the cells. Similarly, concentrated regions of soft cells near the pore regions were found to increase stress concentrations by almost three times compared to the Von-Mises stress generated around the region of cells away from the pores. The study outlines the importance of finite element models in the rapid development of bioinks.


title: Thirty Minutes Before the Dawn -- Trinity

abstract: The Trinity test of July 16, 1945 marked the scientific apex of the Manhattan Project. Often recognized as the symbolic birth of the nuclear age, Trinity's multifaceted legacy remains just as captivating and complex as it did 75 years ago. This paper examines why the test was necessary from a technical standpoint, shows how Los Alamos scientists planned the event, and explores the physical and emotional aftermaths of Trinity. The author also uses rarely accessed original records to reconstruct the story of Trinity's health hazards, as seen through the eyes of radiation technicians and medical doctors as events unfolded. Trinity was conducted as the Potsdam Conference began, weeks after the collapse of Nazi Germany. It was considered necessary to let President Harry S. Truman know whether the United States possessed a nuclear capability ahead of his negotiations with Joseph Stalin, the Soviet premier. The author examines the competing priorities that drove the timetable for the test: international politics, security and safety. Three weeks after Trinity, a gun-assembled enriched-uranium bomb called Little Boy was used against the Japanese city of Hiroshima. Three days later, Fat Man, a weaponized version of the imploding Trinity device, was dropped on Nagasaki. The author briefly examines these strikes and what impact they may have had on the Japanese surrender. The paper concludes by examining the legacy of the Trinity test 75 years into the age it helped usher in.


title: A Multi-resolution Approach to Expression Recognition in the Wild

abstract: Facial expressions play a fundamental role in human communication. Indeed, they typically reveal the real emotional status of people beyond the spoken language. Moreover, the comprehension of human affect based on visual patterns is a key ingredient for any human-machine interaction system and, for such reasons, the task of Facial Expression Recognition (FER) draws both scientific and industrial interest. In the recent years, Deep Learning techniques reached very high performance on FER by exploiting different architectures and learning paradigms. In such a context, we propose a multi-resolution approach to solve the FER task. We ground our intuition on the observation that often faces images are acquired at different resolutions. Thus, directly considering such property while training a model can help achieve higher performance on recognizing facial expressions. To our aim, we use a ResNet-like architecture, equipped with Squeeze-and-Excitation blocks, trained on the Affect-in-the-Wild 2 dataset. Not being available a test set, we conduct tests and models selection by employing the validation set only on which we achieve more than 90\% accuracy on classifying the seven expressions that the dataset comprises.


title: Cell-Laden alginate biomaterial modelling using three-dimensional (3D) microscale finite element technique

abstract: A novel modelling technique using finite element analysis to mimic the mechanoresponse of cell-laden biomaterial is proposed for the use in bioinks and other tissue engineering applications. Here a hydrogel-based composite biomaterial specimen was used consisting of 5% (V/V) HeLa cells added to alginate solution (4% W/V) and another specimen with no living cell present in alginate solution (4% W/V). Tensile test experiments were performed on both the specimens with a load cell of 25 N. The specimens were bioprinted using an in-house developed three-dimensional (3D) bioprinter. To allow for the nonlinear hyperelastic behavior of the specimen, the specimens were loaded very slowly, at rates of 0.1 mm/min and 0.5 mm/min, during the tensile test. The microscale finite element models developed in Ansys were loaded with similar load rates and their responses were recorded. Both the model results were validated with the experiment results. A very good agreement between the finite element model and the tensile test experiment was observed under the same mechanical stimuli. Hence, the study reveals that bioprinted scaffold can be virtually modeled to obtain its mechanical characteristics beforehand.


title: Covering versus partitioning with the Cantor space

abstract: What topological spaces can be partitioned into copies of the Cantor space $2^\omega$? An obvious necessary condition is that a space can be partitioned into copies of $2^\omega$ only if it can be covered with copies of $2^\omega$. We prove three theorems concerning when this necessary condition is also sufficient. If $X$ is a metrizable space and $|X| \leq \mathfrak{c}^{+\omega}$ (the least limit cardinal $>\!\mathfrak{c}$), then $X$ can be partitioned into copies of $2^\omega$ if and only if $X$ can be covered with copies of $2^\omega$. To show this cardinality bound is sharp, we construct a metrizable space of size $\mathfrak{c}^{+(\omega+1)}$ that can be covered with copies of $2^\omega$, but not partitioned into copies of $2^\omega$. Similarly, if $X$ is first countable and $|X| \leq \mathfrak{c}$, then $X$ can be partitioned into copies of $2^\omega$ if and only if $X$ can be covered with copies of $2^\omega$. On the other hand, there is a first countable space of size $\mathfrak{c}^+$ that can be covered with copies of $2^\omega$, but not partitioned into copies of $2^\omega$. Finally, we show that a completely metrizable space can be partitioned into copies of $2^\omega$ if and only if it can be covered with copies of $2^\omega$ if and only if it has no isolated points.


title: Analysis-suitable unstructured T-splines: Multiple extraordinary points per face

abstract: Analysis-suitable T-splines (AST-splines) are a promising candidate to achieve a seamless integration between the design and the analysis of thin-walled structures in industrial settings. In this work, we generalize AST-splines to allow multiple extraordinary points within the same face. This generalization drastically increases the flexibility to build geometries using AST-splines; e.g., much coarser meshes can be generated to represent a certain geometry. The AST-spline spaces detailed in this work have $C^1$ inter-element continuity near extraordinary points and $C^2$ inter-element continuity elsewhere. We mathematically show that AST-splines with multiple extraordinary points per face are linearly independent and their polynomial basis functions form a non-negative partition of unity. We numerically show that AST-splines with multiple extraordinary points per face lead to optimal convergence rates for second- and fourth-order linear elliptic problems. To illustrate a possible isogeometric framework that is already available, we design the B-pillar and the side outer panel of a car using T-splines with the commercial software Autodesk Fusion360, import the control nets into our in-house code to build AST-splines, and import the B\'ezier extraction information into the commercial software LS-DYNA to solve eigenvalue problems. The results are compared with conventional finite elements. Good agreement is found, but conventional finite elements require significantly more degrees of freedom to reach a converged solution than AST-splines.


title: Nuclear Science for the Manhattan Project and Comparison to Today's ENDF Data

abstract: Nuclear physics advances in the US and Britain, from 1939-1945, are described. The Manhattan Project's work led to an explosion in our knowledge of nuclear science. A conference in April 1943 at Los Alamos provided a simple formula used to compute critical masses, and laid out the research program needed to determine the key nuclear constants. In short order, four university accelerators were disassembled and reassembled at Los Alamos, and methods were established to make measurements on extremely small samples owing to the initial lack of availability of enriched 235 U and plutonium. I trace the program that measured fission cross sections, fission emitted neutron multiplicities and their energy spectra, and transport cross sections, comparing the measurements with our best understanding today as embodied in the Evaluated Nuclear Data File ENDF/B-VIII.0. The large nuclear data uncertainties at the beginning of the project, which often exceeded 25--50%, were reduced by 1945, often to less than 5-10%. 235 U and 239 Pu fission cross section assessments in the fast MeV range were reduced with more accurate measurements, and the neutron multiplicity increased. By a lucky coincidence of cancelling errors, the initial critical mass estimates were close to the final estimated masses. Some images from historical documents from our Los Alamos archives are shown. Many of the original measurements from these early years have not previously been widely available. Through this work, these data have now been archived in the international experimental nuclear reaction data library (EXFOR) in a collaboration with the IAEA and Brookhaven National Laboratory.


title: SRG/ART-XC and NuSTAR observations of the X-ray pulsar GRO J1008-57 in the lowest luminosity state

abstract: We report results of the first broadband observation of the transient X-ray pulsar GRO J1008-57 performed in the quiescent state. Observations were conducted quasi-simultaneously with the Mikhail Pavlinsky ART-XC telescope on board SRG and NuSTAR right before the beginning of a Type I outburst. GRO J1008-57 was detected in the state with the lowest observed luminosity around several $\times 10^{34}$ erg s$^{-1}$ and consequently accreting from the cold disk. Timing analysis allowed to significantly detect pulsations during this state for the first time. The observed pulsed fraction of about 20\% is, however, almost three times lower than in brighter states when the accretion proceeds through the standard disk. We traced the evolution of the broadband spectrum of the source on a scale of three orders of magnitude in luminosity and found that at the lowest luminosities the spectrum transforms into the double-hump structure similarly to other X-ray pulsars accreting at low luminosities (X Persei, GX 304-1, A0535+262) reinforcing conclusion that this spectral shape is typical for these objects.


title: Impact of personal income on mortality by age: biological versus socio-economic effects

abstract: The influence of per capita income on life expectancy is well documented, mostly through studies of multinational samples. However, one expects fairly weak correlations at both ends of the life span, that is to say in early infancy and in age groups of elderly from 85 to 100 years. The reason is that at both ends mortality is largely controled by biological factors rather than by socio-economic conditions. In order to test this conjecture, we explore the influence of income on age groups, separately in France, the United States and South Korea. More precisely in each country we compare income and mortality data in as many regional subunits as possible. One noteworthy constatation is that, contrary to a common view, personal income is only weakly correlated with infant mortality (i.e. mortality under the age of one year). More broadly, we propose as a conjecture that the common pattern revealed by the analysis of the three countries is also valid in other developed countries.


title: Structural Connectome Atlas Construction in the Space of Riemannian Metrics

abstract: The structural connectome is often represented by fiber bundles generated from various types of tractography. We propose a method of analyzing connectomes by representing them as a Riemannian metric, thereby viewing them as points in an infinite-dimensional manifold. After equipping this space with a natural metric structure, the Ebin metric, we apply object-oriented statistical analysis to define an atlas as the Fr\'echet mean of a population of Riemannian metrics. We demonstrate connectome registration and atlas formation using connectomes derived from diffusion tensors estimated from a subset of subjects from the Human Connectome Project.


title: Hasse principle violations in twist families of superelliptic curves

abstract: Conditionally on the $abc$ conjecture, we generalize previous work of Clark and the author to show that a superelliptic curve $C: y^n = f(x)$ of sufficiently high genus has infinitely many twists violating the Hasse Principle if and only if $f(x)$ has no $\mathbb Q$-rational roots. We also show unconditionally that a curve defined by $C: y^{pN}=f(x)$ has infinitely many twists violating the Hasse Principle over any number field $k$ such that $k$ contains the $p$th roots of unity and $f(x)$ has no $k$-rational roots.


title: Picosecond acoustic dynamics in stimulated Brillouin scattering

abstract: Recent experiments demonstrating storage of optical pulses in acoustic phonons based on stimulated Brillouin scattering raise a number of questions about the spectral and temporal capacities of such protocols and the limitations of the theoretical frameworks routinely used to describe them. In this work, we consider the dynamics of photon-phonon scattering induced by optical pulses with temporal widths comparable to the period of acoustic oscillations. We revisit the widely adopted classical formalism of coupled modes and demonstrate its breakdown. We propose a simple extension to generalise the formulation and find potentially measurable consequences in the dynamics of Brillouin experiments involving ultra-short pulses.


title: Collective excitations of superfluid Fermi gases near the transition temperature

abstract: Studying the collective pairing phenomena in a two-component Fermi gas, we predict the appearance near the transition temperature $T_c$ of a well-resolved collective mode of quadratic dispersion. The mode is visible both above and below $T_c$ in the system's response to a driving pairing field. When approaching $T_c$ from below, the phononic and pair-breaking branches, characteristic of the zero temperature behavior, reduce to a very low energy-momentum region when the pair correlation length reaches its critical divergent behavior $\xi_{\rm pair}\propto|T_c-T|^{-1/2}$; elsewhere, they are replaced by the quadratically-dispersed pairing resonance, which thus acts as a precursor of the phase transition. In the strong-coupling and Bose-Einstein Condensate regime, this mode is a weakly-damped propagating mode associated to a Lorentzian resonance. Conversely, in the BCS limit it is a relaxation mode of pure imaginary eigenenergy. At large momenta, the resonance disappears when it is reabsorbed by the lower-edge of the pairing continuum. At intermediate temperatures between 0 and $T_c$, we unify the newly found collective phenomena near $T_c$ with the phononic and pair-breaking branches predicted from previous studies, and we exhaustively classify the roots of the analytically continued dispersion equation, and show that they provided a very good summary of the pair spectral functions.


title: Dependence of integrated, instantaneous, and fluctuating entropy production on the initial state in quantum and classical processes

abstract: We consider the additional entropy production (EP) incurred by a fixed quantum or classical process on some initial state $\rho$, above the minimum EP incurred by the same process on any initial state. We show that this additional EP, which we term the "mismatch cost of $\rho$", has a universal information-theoretic form: it is given by the contraction of the relative entropy between $\rho$ and the least-dissipative initial state $\varphi$ over time. We derive versions of this result for integrated EP incurred over the course of a process, for trajectory-level fluctuating EP, and for instantaneous EP rate. We also show that mismatch cost for fluctuating EP obeys an integral fluctuation theorem. Our results demonstrate a fundamental relationship between "thermodynamic irreversibility" (generation of EP) and "logical irreversibility" (inability to know the initial state corresponding to a given final state). We use this relationship to derive quantitative bounds on the thermodynamics of quantum error correction and to propose a thermodynamically-operationalized measure of the logical irreversibility of a quantum channel. Our results hold for both finite and infinite dimensional systems, and generalize beyond EP to many other thermodynamic costs, including nonadiabatic EP, free energy loss, and entropy gain.

categories: quant-ph cond-mat.stat-mech cs.IT math.IT


title: The variable absorption in the X-ray spectrum of GRB 190114C

abstract: GRB 190114C was a bright burst that occurred in the local Universe (z=0.425). It was the first gamma-ray burst (GRB) ever detected at TeV energies, thanks to MAGIC. We characterize the ambient medium properties of the host galaxy through the study of the absorbing X-ray column density. Joining Swift, XMM-Newton, and NuSTAR observations, we find that the GRB X-ray spectrum is characterized by a high column density that is well in excess of the expected Milky Way value and decreases, by a factor of ~2, around ~$10^5$ s. Such a variability is not common in GRBs. The most straightforward interpretation of the variability in terms of photoionization of the ambient medium is not able to account for the decrease at such late times, when the source flux is less intense. Instead, we interpret the decrease as due to a clumped absorber, denser along the line of sight and surrounded by lower-density gas. After the detection at TeV energies of GRB 190114C, two other GRBs were promptly detected. They share a high value of the intrinsic column density and there are hints for a decrease of the column density, too. We speculate that a high local column density might be a common ingredient for TeV-detected GRBs.


title: Dynamic programming equation for the mean field optimal stopping problem

abstract: We study the optimal stopping problem of McKean-Vlasov diffusions when the criterion is a function of the law of the stopped process. A remarkable new feature in this setting is that the stopping time also impacts the dynamics of the stopped process through the dependence of the coefficients on the law. The mean field stopping problem is introduced in weak formulation in terms of the joint marginal law of the stopped underlying process and the survival process. This specification satisfies a dynamic programming principle. The corresponding dynamic programming equation is an obstacle problem on the Wasserstein space, and is obtained by means of a general It\^o formula for flows of marginal laws of c\`adl\`ag semimartingales. Our verification result characterizes the nature of optimal stopping policies, highlighting the crucial need to randomized stopping. The effectiveness of our dynamic programming equation is illustrated by various examples including the mean-variance optimal stopping problem.


title: The AI Arena: A Framework for Distributed Multi-Agent Reinforcement Learning

abstract: Advances in reinforcement learning (RL) have resulted in recent breakthroughs in the application of artificial intelligence (AI) across many different domains. An emerging landscape of development environments is making powerful RL techniques more accessible for a growing community of researchers. However, most existing frameworks do not directly address the problem of learning in complex operating environments, such as dense urban settings or defense-related scenarios, that incorporate distributed, heterogeneous teams of agents. To help enable AI research for this important class of applications, we introduce the AI Arena: a scalable framework with flexible abstractions for distributed multi-agent reinforcement learning. The AI Arena extends the OpenAI Gym interface to allow greater flexibility in learning control policies across multiple agents with heterogeneous learning strategies and localized views of the environment. To illustrate the utility of our framework, we present experimental results that demonstrate performance gains due to a distributed multi-agent learning approach over commonly-used RL techniques in several different learning environments.


title: Multiple zeros of nonlinear systems

abstract: As an attempt to bridge between numerical analysis and algebraic geometry, this paper formulates the multiplicity for the general nonlinear system at an isolated zero, presents an algorithm for computing the multiplicity structure, proposes a depth-deflation method for accurate computation of multiple zeros, and introduces the basic algebraic theory of the multiplicity. Furthermore, this paper elaborates and proves some fundamental properties of the multiplicity, including local finiteness, consistency, perturbation invariance, and depth-deflatability. As a justification of this formulation, the multiplicity is proved to be consistent with the multiplicity defined in algebraic geometry for the special case of polynomial systems. The proposed algorithms can accurately compute the multiplicity and the multiple zeros using floating point arithmetic even if the nonlinear system is perturbed.


title: A convergence framework for optimal transport on the sphere

abstract: We consider a PDE approach to numerically solving the optimal transportation problem on the sphere. We focus on both the traditional squared geodesic cost and a logarithmic cost, which arises in the reflector antenna design problem. At each point on the sphere, we replace the surface PDE with a generalized Monge-Amp\`ere type equation posed on the tangent plane using normal coordinates. The resulting nonlinear PDE can then be approximated by any consistent, monotone scheme for generalized Monge-Amp\`ere type equations on the plane. Existing techniques for proving convergence do not immediately apply because the PDE lacks both a comparison principle and a unique solution, which makes it difficult to produce a stable, well-posed scheme. By augmenting this discretization with an additional term that constrains the solution gradient, we obtain a strong form of stability. A modification of the Barles-Souganidis convergence framework then establishes convergence to the mean-zero solution of the original PDE.


title: $C^*$-algebraic approach to interacting quantum field theory: Inclusion of Fermi fields

abstract: We extend the $C^*$-algebraic approach to interacting quantum field theory, proposed recently by Detlev Buchholz and one of us (KF) to Fermi fields. The crucial feature of our approach is the use of auxiliary Grassmann variables in a functorial way.


title: Non-asymptotic Confidence Intervals of Off-policy Evaluation: Primal and Dual Bounds

abstract: Off-policy evaluation (OPE) is the task of estimating the expected reward of a given policy based on offline data previously collected under different policies. Therefore, OPE is a key step in applying reinforcement learning to real-world domains such as medical treatment, where interactive data collection is expensive or even unsafe. As the observed data tends to be noisy and limited, it is essential to provide rigorous uncertainty quantification, not just a point estimation, when applying OPE to make high stakes decisions. This work considers the problem of constructing non-asymptotic confidence intervals in infinite-horizon off-policy evaluation, which remains a challenging open question. We develop a practical algorithm through a primal-dual optimization-based approach, which leverages the kernel Bellman loss (KBL) of Feng et al.(2019) and a new martingale concentration inequality of KBL applicable to time-dependent data with unknown mixing conditions. Our algorithm makes minimum assumptions on the data and the function class of the Q-function, and works for the behavior-agnostic settings where the data is collected under a mix of arbitrary unknown behavior policies. We present empirical results that clearly demonstrate the advantages of our approach over existing methods.


title: Remarks on Askey-Wilson polynomials and Meixner polynomials of the second kind

abstract: The purpose of this note is twofold: firstly to characterize all the sequences of orthogonal polynomials $(P_n)_{n\geq 0}$ such that $$ \frac{\triangle}{{\bf \triangle} x(s-1/2)}P_{n+1}(x(s-1/2))=c_n(\triangle +2\,\mathrm{I})P_n(x(s-1/2)), $$ where $\mathrm{I}$ is the identity operator, $x$ defines a class of lattices with, generally, nonuniform step-size, and $\triangle f(s)=f(s+1)-f(s)$; and secondly to present, in a friendly way, a method to deal with these kind of problems.


title: Self-induced spin-orbit torques in metallic ferromagnets

abstract: We present a phenomenological theory of spin-orbit torques in a metallic ferromagnet with spin-relaxing boundaries. The model is rooted in the coupled diffusion of charge and spin in the bulk of the ferromagnet, where we account for the anomalous Hall effects as well as the anisotropic magnetoresistance in the corresponding constitutive relations for both charge and spin sectors. The diffusion equations are supplemented with suitable boundary conditions reflecting the spin-sink capacity of the environment. In inversion-asymmetric heterostructures, the uncompensated spin accumulation exerts a dissipative torque on the order parameter, giving rise to a current-dependent linewidth in the ferromagnetic resonance with a characteristic angular dependence. We compare our model to recent spin-torque ferromagnetic resonance measurements, illustrating how rich self-induced spin-torque phenomenology can arise even in simple magnetic structures.


title: Deep neural network approximation for high-dimensional parabolic Hamilton-Jacobi-Bellman equations

abstract: The approximation of solutions to second order Hamilton--Jacobi--Bellman (HJB) equations by deep neural networks is investigated. It is shown that for HJB equations that arise in the context of the optimal control of certain Markov processes the solution can be approximated by deep neural networks without incurring the curse of dimension. The dynamics is assumed to depend affinely on the controls and the cost depends quadratically on the controls. The admissible controls take values in a bounded set.


title: Content-Preserving Unpaired Translation from Simulated to Realistic Ultrasound Images

abstract: Interactive simulation of ultrasound imaging greatly facilitates sonography training. Although ray-tracing based methods have shown promising results, obtaining realistic images requires substantial modeling effort and manual parameter tuning. In addition, current techniques still result in a significant appearance gap between simulated images and real clinical scans. Herein we introduce a novel content-preserving image translation framework (ConPres) to bridge this appearance gap, while maintaining the simulated anatomical layout. We achieve this goal by leveraging both simulated images with semantic segmentations and unpaired in-vivo ultrasound scans. Our framework is based on recent contrastive unpaired translation techniques and we propose a regularization approach by learning an auxiliary segmentation-to-real image translation task, which encourages the disentanglement of content and style. In addition, we extend the generator to be class-conditional, which enables the incorporation of additional losses, in particular a cyclic consistency loss, to further improve the translation quality. Qualitative and quantitative comparisons against state-of-the-art unpaired translation methods demonstrate the superiority of our proposed framework.


title: Analyzing Human Models that Adapt Online

abstract: Predictive human models often need to adapt their parameters online from human data. This raises previously ignored safety-related questions for robots relying on these models such as what the model could learn online and how quickly could it learn it. For instance, when will the robot have a confident estimate in a nearby human's goal? Or, what parameter initializations guarantee that the robot can learn the human's preferences in a finite number of observations? To answer such analysis questions, our key idea is to model the robot's learning algorithm as a dynamical system where the state is the current model parameter estimate and the control is the human data the robot observes. This enables us to leverage tools from reachability analysis and optimal control to compute the set of hypotheses the robot could learn in finite time, as well as the worst and best-case time it takes to learn them. We demonstrate the utility of our analysis tool in four human-robot domains, including autonomous driving and indoor navigation.


title: Asymptotic posterior normality of the generalized extreme value distribution

abstract: The univariate generalized extreme value (GEV) distribution is the most commonly used tool for analyzing the properties of rare events. The ever greater utilization of Bayesian methods for extreme value analysis warrants detailed theoretical investigation, which has thus far been underdeveloped. Even the most basic asymptotic results are difficult to obtain because the GEV fails to satisfy standard regularity conditions. Here, we prove that the posterior distribution of the GEV parameter vector, given $n$ independent and identically distributed samples, converges in distribution to a trivariate normal distribution. The proof necessitates analyzing integrals of the GEV likelihood function over the entire parameter space, which requires considerable care because the support of the GEV density depends on the parameters in complicated ways.


title: Apprentice for Event Generator Tuning

abstract: Apprentice is a tool developed for event generator tuning. It contains a range of conceptual improvements and extensions over the tuning tool Professor. Its core functionality remains the construction of a multivariate analytic surrogate model to computationally expensive Monte-Carlo event generator predictions. The surrogate model is used for numerical optimization in chi-square minimization and likelihood evaluation. Apprentice also introduces algorithms to automate the selection of observable weights to minimize the effect of mis-modeling in the event generators. We illustrate our improvements for the task of MC-generator tuning and limit setting.


title: Connecting Latent ReLationships over Heterogeneous Attributed Network for Recommendation

abstract: Recently, deep neural network models for graph-structured data have been demonstrating to be influential in recommendation systems. Graph Neural Network (GNN), which can generate high-quality embeddings by capturing graph-structured information, is convenient for the recommendation. However, most existing GNN models mainly focus on the homogeneous graph. They cannot characterize heterogeneous and complex data in the recommendation system. Meanwhile, it is challenging to develop effective methods to mine the heterogeneity and latent correlations in the graph. In this paper, we adopt Heterogeneous Attributed Network (HAN), which involves different node types as well as rich node attributes, to model data in the recommendation system. Furthermore, we propose a novel graph neural network-based model to deal with HAN for Recommendation, called HANRec. In particular, we design a component connecting potential neighbors to explore the influence between neighbors and provide two different strategies with the attention mechanism to aggregate neighbors' information. The experimental results on two real-world datasets prove that HANRec outperforms other state-of-the-art methods.