title: Recognition of Ischaemia and Infection in Diabetic Foot Ulcers: Dataset and Techniques

abstract: Recognition and analysis of Diabetic Foot Ulcers (DFU) using computerized methods is an emerging research area with the evolution of image-based machine learning algorithms. Existing research using visual computerized methods mainly focuses on recognition, detection, and segmentation of the visual appearance of the DFU as well as tissue classification. According to DFU medical classification systems, the presence of infection (bacteria in the wound) and ischaemia (inadequate blood supply) has important clinical implications for DFU assessment, which are used to predict the risk of amputation. In this work, we propose a new dataset and computer vision techniques to identify the presence of infection and ischaemia in DFU. This is the first time a DFU dataset with ground truth labels of ischaemia and infection cases is introduced for research purposes. For the handcrafted machine learning approach, we propose a new feature descriptor, namely the Superpixel Color Descriptor. Then we use the Ensemble Convolutional Neural Network (CNN) model for more effective recognition of ischaemia and infection. We propose to use a natural data-augmentation method, which identifies the region of interest on foot images and focuses on finding the salient features existing in this area. Finally, we evaluate the performance of our proposed techniques on binary classification, i.e. ischaemia versus non-ischaemia and infection versus non-infection. Overall, our method performed better in the classification of ischaemia than infection. We found that our proposed Ensemble CNN deep learning algorithms performed better for both classification tasks as compared to handcrafted machine learning algorithms, with 90% accuracy in ischaemia classification and 73% in infection classification.


title: JEDI-net: a jet identification algorithm based on interaction networks

abstract: We investigate the performance of a jet identification algorithm based on interaction networks (JEDI-net) to identify all-hadronic decays of high-momentum heavy particles produced at the LHC and distinguish them from ordinary jets originating from the hadronization of quarks and gluons. The jet dynamics are described as a set of one-to-one interactions between the jet constituents. Based on a representation learned from these interactions, the jet is associated to one of the considered categories. Unlike other architectures, the JEDI-net models achieve their performance without special handling of the sparse input jet representation, extensive pre-processing, particle ordering, or specific assumptions regarding the underlying detector geometry. The presented models give better results with less model parameters, offering interesting prospects for LHC applications.


title: A grouped, selectively weighted false discovery rate procedure

abstract: False discovery rate (FDR) control in structured hypotheses testing is an important topic in simultaneous inference. Most existing methods that aim to utilize group structure among hypotheses either employ the groupwise mixture model or weight all p-values or hypotheses. Thus, their powers can be improved when the groupwise mixture model is inappropriate or when most groups contain only true null hypotheses. Motivated by this, we propose a grouped, selectively weighted FDR procedure, which we refer to as "sGBH". Specifically, without employing the groupwise mixture model, sGBH identifies groups of hypotheses of interest, weights p-values in each such group only, and tests only the selected hypotheses using the weighted p-values. The sGBH subsumes a standard grouped, weighted FDR procedure which we refer to as "GBH". We provide simple conditions to ensure the conservativeness of sGBH, together with empirical evidence on its much improved power over GBH. The new procedure is applied to a gene expression study.


title: $b \bar b$ Kinematic Correlations in Cold Nuclear Matter

abstract: Background: The LHCb Collaboration has studied a number of kinematic correlations between $B$-hadron pairs through their subsequent decays to $J/\psi$ pairs at 7 and 8 TeV for four minimum values of the $J/\psi$ $p_T$. Purpose: In this work, these measurements are compared to calculations of $b \bar b$ pairs and their hadronization and inclusive decays to $J/\psi J/\psi$ are compared to the same observables. Potential cold matter effects on the $b \bar b$ pair observables are discussed to determine which are most likely to provide insights about the system and why. Methods: The calculations, employing the exclusive HVQMNR code, assume the same intrinsic $k_T$-broadening and fragmentation as in [R. Vogt, Phys. Rev. C {\bf 98} (2018) 034907]. The pair distributions presented by LHCb are calculated in this approach, both for the parent $b \bar b$ and the $J/\psi J/\psi$ pairs produced in their decay. The sensitivity of the results to the intrinsic $k_T$ broadening is shown. The theoretical uncertainties due to the $b$ quark mass and scale variations on both the initial $b \bar b$ pairs and the resulting $J/\psi$ pairs are also shown. Possible effects due to the presence of the nucleus are studied by increasing the size of the $k_T$ broadening and modification of the fragmentation parameter. Results: Good agreement with the LHCb data is found for all observables. The parent $b \bar b$ distributions are more sensitive to the $k_T$ broadening than are the final-state $J/\psi$ pairs. Conclusions: Next-to-leading order calculations with $k_T$ broadening, as in [R. Vogt, Phys. Rev. C {\bf 98} (2018) 034907], can describe all correlated observables. Multiple measurements of correlated observables are sensitive to different nuclear effects which can help distinguish between them.


title: The conjugacy growth of the soluble Baumslag-Solitar groups

abstract: In this paper we give asymptotics for the conjugacy growth of the soluble Baumslag-Solitar groups $BS(1,k)$, $k\geq 2$, with respect to the standard generating set, by providing a complete description of geodesic conjugacy representatives. We show that the conjugacy growth series for these groups are transcendental, and give formulas for the series. As a result of our computation we also establish that in each $BS(1,k)$ the conjugacy and standard growth rates are equal.


title: On Strominger K\"ahler-like manifolds with degenerate torsion

abstract: In this paper, we study a special type of compact Hermitian manifolds that are Strominger K\"ahler-like, or SKL for short. This condition means that the Strominger connection (also known as Bismut connection) is K\"ahler-like, in the sense that its curvature tensor obeys all the symmetries of the curvature of a K\"ahler manifold. Previously, we have shown that any SKL manifold $(M^n,g)$ is always pluriclosed, and when the manifold is compact and $g$ is not K\"ahler, it can not admit any balanced or strongly Gauduchon (in the sense of Popovici) metric. Also, when $n=2$, the SKL condition is equivalent to the Vaisman condition. In this paper, we give a classification for compact non-K\"ahler SKL manifolds in dimension $3$ and those with degenerate torsion in higher dimensions. We also present some properties about SKL manifolds in general dimensions, for instance, given any compact non-K\"ahler SKL manifold, its K\"ahler form represents a non-trivial Aeppli cohomology class, the metric can never be locally conformal K\"ahler when $n\geq 3$, and the manifold does not admit any Hermitian symplectic metric.


title: On Separating Points for Ensemble Controllability

abstract: Recent years have witnessed a wave of research activities in systems science toward the study of population systems. The driving force behind this shift was geared by numerous emerging and ever-changing technologies in life and physical sciences and engineering, from neuroscience, biology, and quantum physics to robotics, where many control-enabled applications involve manipulating a large ensemble of structurally identical dynamic units, or agents. Analyzing fundamental properties of ensemble control systems in turn plays a foundational and critical role in enabling and, further, advancing these applications, and the analysis is largely beyond the capability of classical control techniques. In this paper, we consider an ensemble of time-invariant linear systems evolving on an infinite-dimensional space of continuous functions. We exploit the notion of separating points and techniques of polynomial approximation to develop necessary and sufficient ensemble controllability conditions. In particular, we introduce an extended notion of controllability matrix, called Ensemble Controllability Gramian. This means enables the characterization of ensemble controllability through evaluating controllability of each individual system in the ensemble. As a result, the work provides a unified framework with a systematic procedure for analyzing control systems defined on an infinite-dimensional space by a finite-dimensional approach.


title: Dual Adversarial Inference for Text-to-Image Synthesis

abstract: Synthesizing images from a given text description involves engaging two types of information: the content, which includes information explicitly described in the text (e.g., color, composition, etc.), and the style, which is usually not well described in the text (e.g., location, quantity, size, etc.). However, in previous works, it is typically treated as a process of generating images only from the content, i.e., without considering learning meaningful style representations. In this paper, we aim to learn two variables that are disentangled in the latent space, representing content and style respectively. We achieve this by augmenting current text-to-image synthesis frameworks with a dual adversarial inference mechanism. Through extensive experiments, we show that our model learns, in an unsupervised manner, style representations corresponding to certain meaningful information present in the image that are not well described in the text. The new framework also improves the quality of synthesized images when evaluated on Oxford-102, CUB and COCO datasets.


title: Highlights of Galactic Physics with VERITAS

abstract: VERITAS (Very Energetic Radiation Imaging Telescope Array System) is one of the most sensitive currently operating arrays of imaging atmospheric Cherenkov telescopes, which detect very high-energy (VHE; E > 100 GeV) gamma rays. VERITAS is currently in its 11th year of full-array operations with four 12m-diameter telescopes. Many Galactic sources of VHE gamma rays have been detected by VERITAS, such as pulsar wind nebulae, binary systems, and supernova remnants, and the study of VHE emission from these objects has enabled a deeper understanding of the underlying physical processes responsible for the observed gamma rays. Recent highlights from the VERITAS Galactic science program will be presented, including results on pulsar searches, follow-up of sources detected by HAWC, and the 50-year-period binary PSR J2032+4127.


title: Some Gruss type inequalities for Frechet differentiable mappings

abstract: Let X be a Hilbert C^*-module on C^*-algebra A and p in A. We denote by Dp(A;X) the set of all continuous functions f on A, which are Frechet differentiable on a open neighborhood U of p. Then, we introduce some generalized semi-inner products on Dp(A;X), and using them some Gruss type inequalities in semi-inner product C^*-module Dp(A;X) and Dp(A;X^n) are established.


title: Density Effects on the Post-shock Turbulence Structure and Dynamics

abstract: Turbulence structure resulting from multi-fluid or multi-species, variable-density isotropic turbulence interaction with a Mach 2 shock is studied using turbulence-resolving shock-capturing simulations and Eulerian (grid) and Lagrangian (particle) methods. The complex roles density play in the modification of turbulence by the shock wave are identified. Statistical analyses of the velocity gradient tensor (VGT) show that the density variations significantly change the turbulence structure and flow topology. Specifically, a stronger symmetrization of the joint probability density function (PDF) of second and third invariants of the anisotropic velocity gradient tensor, PDF$(Q^\ast, R^\ast)$, as well as the PDF of the vortex stretching contribution to the enstrophy equation, are observed in the multi-species case. Furthermore, subsequent to the interaction with the shock, turbulent statistics also acquire a differential distribution in regions having different densities. This results in a nearly symmetrical PDF$(Q^\ast, R^\ast)$ in heavy fluid regions, while the light fluid regions retain the characteristic tear-drop shape. To understand this behavior and the return to "standard" turbulence structure as the flow evolves away from the shock, Lagrangian dynamics of the VGT and its invariants are studied by considering particle residence times and conditional particle variables in different flow regions. The pressure Hessian contributions to the VGT invariants transport equations are shown to be not only affected by the shock wave, but also by the density in the multi-fluid case, making them critically important to the flow dynamics and turbulence structure.


title: Quenched Dynamics of Artificial Spin Ice: Coarsening versus Kibble-Zurek

abstract: Artificial spin ices are ideal frustrated model systems in which to explore or design emergent phenomena with unprecedented characterization of the constituent degrees of freedom. In square spin ice, violations of the ice rule are topological excitations essential to the kinetics of the system, providing an ideal testbed for studying the dynamics of such defects under varied quench rates. In this work we describe the first test of the Kibble-Zurek mechanism and critical coarsening in colloidal square and colloidal hexagonal ice under quenches from a weakly interacting liquid state into a strongly interacting regime. As expected, for infinitely slow quenches, the system is defect free, while for increasing quench rate, an increasing number of defects remain in the sample. For square ice, we find regimes in which the defect population decreases as a power law with decreasing quench rate. A detailed scaling analysis shows that for a wide range of parameters, including quench rates that are accessible by experiments, the behavior is described by critical coarsening rather than by the Kibble-Zurek mechanism, since the defect-defect interactions are long ranged. For quenches closer to the critical point, however, there can be a competition between the two mechanisms.


title: Modelling Mg II During Solar Flares. I. Partial Frequency Redistribution, Opacity, and Coronal Irradiation

abstract: The Interface Region Imaging Spectrograph (IRIS) has routinely observed the flaring Mg II NUV spectrum, offering excellent diagnostic potential and a window into the location of energy deposition. A number of studies have forward modelled both the general properties of these lines and specific flare observations. Generally these have forward modelled radiation via post-processing of snapshots from hydrodynamic flare simulations through radiation transfer codes. There has, however, not been a study of how the physics included in these radiation transport codes affects the solution. A baseline setup for forward modelling MgII in flares is presented and contrasted with approaches that add or remove complexity. It is shown for Mg II: (1) PRD is still required during flare simulations despite the increased densities, (2) using full angle-dependent PRD affects the solution but takes significantly longer to process a snapshot, (3) including Mg I in NLTE results in negligible differences to the Mg II lines but does affect the NUV quasi-continuum, (4) only hydrogen and Mg II need to be included in NLTE, (5) ideally the non-equilibrium hydrogen populations, with non-thermal collisional rates, should be used rather than the statistical equilibrium populations, (6) an atom consisting of only the ground state, h & k upper levels, and continuum level is insufficient to model the resonance lines, and (7) irradiation from a hot, dense flaring transition region can affect the formation of Mg II. We discuss modifications to the RH code allowing straightforward inclusion of transition region and coronal irradiation in flares.


title: Numerical solution for tachyon vacuum in the Schnabl gauge

abstract: Based on the level truncation scheme, we develop a new numerical method to evaluate the tachyon vacuum solution in the Schnabl gauge up to level $L=24$. We confirm the prediction that the energy associated to this numerical solution has a local minimum at level $L=12$. Extrapolating the energy data of $L \leq 24$ to infinite level, we observe that the energy goes towards the analytical value $-1$, nevertheless the precision of the extrapolation is lower than in the Siegel gauge. Furthermore, we analyze the Ellwood invariant and show that its value converges monotonically towards the expected analytical result. We also study the tachyon vacuum expectation value (vev) and some other coefficients of the solution. Finally, some consistency checks of the solution are performed, and we briefly discuss the search for other Schnabl gauge numerical solutions.


title: Dimesoatom breakup in the Coulomb field

abstract: Momentum and angular distributions of charged meson pairs $h^+h^-$ ($h=\pi, K$) from elementary atoms (EA) breakup (ionization) in the Coulomb field of a target atom is considered in the Born and Glauber approximations. Exploiting the fact that the atomic screening of the target Coulomb potential is important at small transfer momenta, while multi-photon exchanges are essential at large transfer momenta we express the cross sections of EA breakup as a sum of two terms. In the region of modest transfer momenta the cross section is determined by the single-photon exchange (first Born approximation) accounting for the target atoms screening, whereas at large transfer momenta using the unscreened potential allows to take into account all multi-photon exchanges and obtain the cross section of EA breakup in the close analytical form.


title: Interacting Agegraphic Dark Energy Model in DGP Braneworld Cosmology: Dynamical System Approach

abstract: A proposal to study the effect of interaction in an agegraphic dark energy model in DGP brane-world cosmology is presented in this manuscript. After explaining the details, we proceed to apply the dynamical system approach to the model to analyze its stability. We first, constrain model parameters with a variety of independent observational data such as cosmic microwave background anisotropies, baryon acoustic oscillation peaks and observational Hubble data. Then, we obtain the critical points related to different cosmological epochs. In particular, we conclude that in the presence of interaction, dark energy dominated era could be a stable point if model parameters $n$ and $\beta$, obey a given constraint. Also, big rip singularity is avoidable in this model.


title: Families of stable bundles on the fibres of the hyperk\"ahler twistor projection

abstract: Given a holomorphic vector bundle $E$ on the twistor space $\mathrm{Tw}(M)$ of a simple hyperk\"ahler manifold $M$, we view it as a family of bundles $\left\{E_I\right\}$ on the fibres $\pi^{-1}(I)$ of the twistor projection $\pi : \mathrm{Tw}(M) \to \mathbb{CP}^1$, and study the relationship between stability of $E$ and its fibrewise stability. We verify that the argument of Teleman establishing the Zariski openness of stability and semi-stability in families of bundles applies in the case of the family $\left\{E_I\right\}$. We prove a partial converse to a result of Kaledin and Verbitsky, showing that an irreducible bundle $E$ on $\mathrm{Tw}(M)$ is generically fibrewise stable if the rank of $E$ is 2 or 3, or at least one element of the family $\left\{E_I\right\}$ is a simple bundle, in the sense that $\mathrm{Hom}(E_I, E_I) = \mathbb{C}$.


title: Adjoining only the things you want: a survey of Strong Chang's Conjecture and related topics

abstract: We survey some old and new results on strong variants of Chang's Conjecture and related topics.


title: First resolved observations of a highly asymmetric debris disc around HD 160305 with VLT/SPHERE

abstract: Context. Direct imaging of debris discs gives important information about their nature, their global morphology, and allows us to identify specific structures possibly in connection with the presence of gravitational perturbers. It is the most straightforward technique to observe planetary systems as a whole. Aims. We present the first resolved images of the debris disc around the young F-type star HD 160305, detected in scattered light using the VLT/SPHERE instrument in the near infrared. Methods. We used a post-processing method based on angular differential imaging and synthetic images of debris discs produced with a disc modelling code (GRaTer) to constrain the main characteristics of the disc around HD 160305. All of the point sources in the field of the IRDIS camera were analysed with an astrometric tool to determine whether they are bound objects or background stars. Results. We detect a very inclined (~ 82{\deg}) ring-like debris disc located at a stellocentric distance of about 86au (deprojected width ~27 au). The disc displays a brightness asymmetry between the two sides of the major axis, as can be expected from scattering properties of dust grains. We derive an anisotropic scattering factor g>0.5. A second right-left asymmetry is also observed with respect to the minor axis. We measure a surface brightness ratio of 0.73 $\pm$ 0.18 between the bright and the faint sides. Because of the low signal-to-noise ratio (S/N) of the images we cannot easily discriminate between several possible explanations for this left-right asymmetry, such as perturbations by an unseen planet, the aftermath of the breakup of a massive planetesimal, or the pericenter glow effect due to an eccentric ring. Two epochs of observations allow us to reject the companionship hypothesis for the 15 point sources present in the field.


title: Towers and clubs

abstract: We revisit several results concerning club principles and nonsaturation of the nonstationary ideal, attempting to improve them in various ways. So we typically deal with a (non necessarily normal) ideal $J$ extending the nonstationary ideal on a regular uncountable (non necessarily successor) cardinal $\kappa$, our goal being to witness the nonsaturation of $J$ by the existence of towers (of length possibly greater than $\kappa^+$).


title: Quantum Cosmology with many fluids and the choice of cosmological time

abstract: In this work we propose the quantization of a cosmological model describing the primordial universe filled with five barotropic fluids, namely: radiation, dust, vacuum, cosmic strings and domain walls. We intend to identify which fluid is best suited to provide phenomenologically the temporal variable in accordance with the observable universe. Through the Galerkin spectral method and the finite difference method in the Crank-Nicolson scheme (vacuum case), the quantum cosmological solutions are obtained and compared. We, also, compare the quantum cosmological solutions with the corresponding classical ones. The vacuum case is especially interesting because it provides a tunneling transition mechanism from the quantum to the classical phase and the possibility of calculating quantum tunneling probabilities.


title: Robust parametric modeling of Alzheimer's disease progression

abstract: Quantitative characterization of disease progression using longitudinal data can provide long-term predictions for the pathological stages of individuals. This work studies the robust modeling of Alzheimer's disease progression using parametric methods. The proposed method linearly maps the individual's age to a disease progression score (DPS) and jointly fits constrained generalized logistic functions to the longitudinal dynamics of biomarkers as functions of the DPS using M-estimation. Robustness of the estimates is quantified using bootstrapping via Monte Carlo resampling, and the estimated inflection points of the fitted functions are used to temporally order the modeled biomarkers in the disease course. Kernel density estimation is applied to the obtained DPSs for clinical status classification using a Bayesian classifier. Different M-estimators and logistic functions, including a novel type proposed in this study, called modified Stannard, are evaluated on the data from the Alzheimer's Disease Neuroimaging Initiative (ADNI) for robust modeling of volumetric MRI and PET biomarkers, CSF measurements, as well as cognitive tests. The results show that the modified Stannard function fitted using the logistic loss achieves the best modeling performance with an average normalized MAE of 0.991 across all biomarkers and bootstraps. Applied to the ADNI test set, this model achieves a multiclass AUC of 0.934 in clinical status classification. The obtained results for the proposed model outperform almost all state-of-the-art results in predicting biomarker values and classifying clinical status. Finally, the experiments show that the proposed model, trained using abundant ADNI data, generalizes well to data from the National Alzheimer's Coordinating Center (NACC) with an average normalized MAE of 1.182 and a multiclass AUC of 0.929.


title: Mixed pooling of seasonality for time series forecasting: An application to pallet transport data

abstract: Multiple seasonal patterns play a key role in time series forecasting, especially for business time series where seasonal effects are often dramatic. Previous approaches including Fourier decomposition, exponential smoothing, and seasonal autoregressive integrated moving average (SARIMA) models do not reflect the distinct characteristics of each period in seasonal patterns. We propose a mixed hierarchical seasonality (MHS) model. Intermediate parameters for each seasonal period are first estimated, and a mixture of intermediate parameters is taken. This results in a model that automatically learns the relative importance of each seasonality and addresses the interactions between them. The model is implemented with Stan, a probabilistic language, and was compared with three existing models on a real-world dataset of pallet transport from a logistic network. Our new model achieved considerable improvements in terms of out of sample prediction error (MAPE) and predictive density (ELPD) compared to complete pooling, Fourier decomposition, and SARIMA model.


title: A hierarchical model for estimating exposure-response curves from multiple studies

abstract: Cookstove replacement trials have found mixed results on their impact on respiratory health. The limited range of concentrations and small sample sizes of individual studies are important factors that may be limiting their statistical power. We present a hierarchical approach to modeling exposure concentrations and pooling data from multiple studies in order to estimate a common exposure-response curve. The exposure concentration model accommodates temporally sparse, clustered longitudinal observations. The exposure-response curve model provides a flexible, semi-parametric estimate of the exposure-response relationship while accommodating heterogeneous clustered data. We apply this model to data from three studies of cookstoves and respiratory infections in children in Nepal, which represent three study types: crossover trial, parallel trial, and case-control study. We find evidence of increased odds of disease for particulate matter concentrations between 50 and 200 $\mu$g/m$^3$ and a flattening of the exposure-response curve for higher exposure concentrations. The model we present can incorporate additional studies and be applied to other settings.


title: Evaluating Empathy in Artificial Agents

abstract: The novel research area of computational empathy is in its infancy and moving towards developing methods and standards. One major problem is the lack of agreement on the evaluation of empathy in artificial interactive systems. Even though the existence of well-established methods from psychology, psychiatry and neuroscience, the translation between these methods and computational empathy is not straightforward. It requires a collective effort to develop metrics that are more suitable for interactive artificial agents. This paper is aimed as an attempt to initiate the dialogue on this important problem. We examine the evaluation methods for empathy in humans and provide suggestions for the development of better metrics to evaluate empathy in artificial agents. We acknowledge the difficulty of arriving at a single solution in a vast variety of interactive systems and propose a set of systematic approaches that can be used with a variety of applications and systems.


title: Uncharacteristic second order martensitic transformation in metals via epitaxial stress fields

abstract: While most phase transformations, e.g. ferroelectric or ferromagnetic, can be first or second order depending on external applied fields, martensitic transformations in metallic alloys are nearly universally first order. We demonstrate that epitaxial stress originating from the incorporation of a tailored second phase can modify the free energy landscape that governs the phase transition and change its order from first to second. High-fidelity molecular dynamics simulations show a remarkable change in the character of the martensitic transformation in Ni-Al alloys near the critical point. We observe the continuous evolution of the transformation order parameter and scaling with power-law exponents comparable to those in other ferroic transitions exhibiting critical behavior. Our theoretical work provides a foundation to recent experimental and computational results on martensites near critical points.


title: Graph Convolutional Networks for Coronary Artery Segmentation in Cardiac CT Angiography

abstract: Detection of coronary artery stenosis in coronary CT angiography (CCTA) requires highly personalized surface meshes enclosing the coronary lumen. In this work, we propose to use graph convolutional networks (GCNs) to predict the spatial location of vertices in a tubular surface mesh that segments the coronary artery lumen. Predictions for individual vertex locations are based on local image features as well as on features of neighboring vertices in the mesh graph. The method was trained and evaluated using the publicly available Coronary Artery Stenoses Detection and Quantification Evaluation Framework. Surface meshes enclosing the full coronary artery tree were automatically extracted. A quantitative evaluation on 78 coronary artery segments showed that these meshes corresponded closely to reference annotations, with a Dice similarity coefficient of 0.75/0.73, a mean surface distance of 0.25/0.28 mm, and a Hausdorff distance of 1.53/1.86 mm in healthy/diseased vessel segments. The results showed that inclusion of mesh information in a GCN improves segmentation overlap and accuracy over a baseline model without interaction on the mesh. The results indicate that GCNs allow efficient extraction of coronary artery surface meshes and that the use of GCNs leads to regular and more accurate meshes.


title: Raw-to-End Name Entity Recognition in Social Media

abstract: Taking word sequences as the input, typical named entity recognition (NER) models neglect errors from pre-processing (e.g., tokenization). However, these errors can influence the model performance greatly, especially for noisy texts like tweets. Here, we introduce Neural-Char-CRF, a raw-to-end framework that is more robust to pre-processing errors. It takes raw character sequences as inputs and makes end-to-end predictions. Word embedding and contextualized representation models are further tailored to capture textual signals for each character instead of each word. Our model neither requires the conversion from character sequences to word sequences, nor assumes tokenizer can correctly detect all word boundaries. Moreover, we observe our model performance remains unchanged after replacing tokenization with string matching, which demonstrates its potential to be tokenization-free. Extensive experimental results on two public datasets demonstrate the superiority of our proposed method over the state of the art. The implementations and datasets are made available at: https://github.com/LiyuanLucasLiu/Raw-to-End.


title: On Strict (Outer-)Confluent Graphs

abstract: A strict confluent (SC) graph drawing is a drawing of a graph with vertices as points in the plane, where vertex adjacencies are represented not by individual curves but rather by unique smooth paths through a planar system of junctions and arcs. If all vertices of the graph lie in the outer face of the drawing, the drawing is called a strict outerconfluent (SOC) drawing. SC and SOC graphs were first considered by Eppstein et al. in Graph Drawing 2013. Here, we establish several new relationships between the class of SC graphs and other graph classes, in particular string graphs and unit-interval graphs. Further, we extend earlier results about special bipartite graph classes to the notion of strict outerconfluency, show that SOC graphs have cop number two, and establish that tree-like ($\Delta$-)SOC graphs have bounded cliquewidth.


title: Measuring the ionisation fraction in a jet from a massive protostar

abstract: It is important to determine if massive stars form via disc accretion, like their low-mass counterparts. Theory and observation indicate that protostellar jets are a natural consequence of accretion discs and are likely to be crucial for removing angular momentum during the collapse. However, massive protostars are typically rarer, more distant and more dust enshrouded, making observational studies of their jets more challenging. A fundamental question is whether the degree of ionisation in jets is similar across the mass spectrum. Here we determine an ionisation fraction of $\sim5-12\%$ in the jet from the massive protostar G35.20-0.74N, based on spatially coincident infrared and radio emission. This is similar to the values found in jets from lower-mass young stars, implying a unified mechanism of shock ionisation applies in jets across most of the protostellar mass spectrum, up to at least $\sim10$ solar masses.


title: UAV Surveillance Under Visibility and Dwell-Time Constraints: A Sampling-Based Approach

abstract: A framework is introduced for planning unmanned aerial vehicle flight paths for visual surveillance of ground targets, each having particular viewing requirements. Specifically, each target is associated with a set of imaging parameters, including a desired (i) tilt angle, (ii) azimuth, with the option of a 360-degree view, and (iii) dwell-time. Tours are sought to image the targets, while minimizing both the total mission time and the time required to reach the initial target. An epsilon-constraint scalarization is used to pose the multi-objective problem as a constrained optimization, which, through careful discretization, can be approximated as a discrete graph-search. It is shown that, in many cases, this approximation is equivalent to a generalized traveling salesperson problem. A heuristic procedure for solving the discrete approximation and recovering solutions to the full routing problem is presented, and is shown to have resolution completeness properties. Algorithms are illustrated through numerical studies.


title: From Crystallized Adaptivity to Fluid Adaptivity in Deep Reinforcement Learning -- Insights from Biological Systems on Adaptive Flexibility

abstract: Recent developments in machine-learning algorithms have led to impressive performance increases in many traditional application scenarios of artificial intelligence research. In the area of deep reinforcement learning, deep learning functional architectures are combined with incremental learning schemes for sequential tasks that include interaction-based, but often delayed feedback. Despite their impressive successes, modern machine-learning approaches, including deep reinforcement learning, still perform weakly when compared to flexibly adaptive biological systems in certain naturally occurring scenarios. Such scenarios include transfers to environments different than the ones in which the training took place or environments that dynamically change, both of which are often mastered by biological systems through a capability that we here term "fluid adaptivity" to contrast it from the much slower adaptivity ("crystallized adaptivity") of the prior learning from which the behavior emerged. In this article, we derive and discuss research strategies, based on analyzes of fluid adaptivity in biological systems and its neuronal modeling, that might aid in equipping future artificially intelligent systems with capabilities of fluid adaptivity more similar to those seen in some biologically intelligent systems. A key component of this research strategy is the dynamization of the problem space itself and the implementation of this dynamization by suitably designed flexibly interacting modules.


title: Multimodal Emotion Recognition Using Deep Canonical Correlation Analysis

abstract: Multimodal signals are more powerful than unimodal data for emotion recognition since they can represent emotions more comprehensively. In this paper, we introduce deep canonical correlation analysis (DCCA) to multimodal emotion recognition. The basic idea behind DCCA is to transform each modality separately and coordinate different modalities into a hyperspace by using specified canonical correlation analysis constraints. We evaluate the performance of DCCA on five multimodal datasets: the SEED, SEED-IV, SEED-V, DEAP, and DREAMER datasets. Our experimental results demonstrate that DCCA achieves state-of-the-art recognition accuracy rates on all five datasets: 94.58% on the SEED dataset, 87.45% on the SEED-IV dataset, 84.33% and 85.62% for two binary classification tasks and 88.51% for a four-category classification task on the DEAP dataset, 83.08% on the SEED-V dataset, and 88.99%, 90.57%, and 90.67% for three binary classification tasks on the DREAMER dataset. We also compare the noise robustness of DCCA with that of existing methods when adding various amounts of noise to the SEED-V dataset. The experimental results indicate that DCCA has greater robustness. By visualizing feature distributions with t-SNE and calculating the mutual information between different modalities before and after using DCCA, we find that the features transformed by DCCA from different modalities are more homogeneous and discriminative across emotions.


title: Random sequential adsorption of spheres on a cylinder

abstract: Inspired by observations of beads packed on a thin string in such systems as sea-grapes and dental plaque, we study the random sequential adsorption of spheres on a cylinder. We determine the asymptotic fractional coverage of the cylinder as a function of the sole parameter in the problem, the ratio of the sphere radius to the cylinder radius (for a very long cylinder) using a combination of analysis and numerical simulations. Examining the asymptotic structures, we find weak chiral ordering on sufficiently small spatial scales. Experiments involving colloidal microspheres that can attach irreversibly to a silica wire via electrostatic forces or DNA hybridization allow us to verify our predictions for the asymptotic coverage.


title: Experimental quantum repeater without quantum memory

abstract: Quantum repeaters -- important components of a scalable quantum internet -- enable the entanglement to be distributed over long distances. The standard paradigm for a quantum repeater relies on a necessary demanding requirement of quantum memory. Despite significant progress, the limited performance of quantum memory makes practical quantum repeaters still a great challenge. Remarkably, a proposed all-photonic quantum repeater avoids the need for quantum memory by harnessing the graph states in the repeater nodes. Here we perform an experimental demonstration of an all-photonic quantum repeater using linear optics. By manipulating a 12-photon interferometer, we implement a 2-by-2 parallel all-photonic quantum repeater, and observe an 89% enhancement of entanglement-generation rate over the standard parallel entanglement swapping. These results open a new way towards designing repeaters with efficient single-photon sources and photonic graph states, and suggest that the all-photonic scheme represents an alternative path -- parallel to that of matter-memory-based schemes -- towards realizing practical quantum repeaters.


title: A 3-Dimensional Likelihood analysis method for detecting extended sources in VERITAS

abstract: Gamma ray observations from a few hundred MeV up to tens of TeV are a valuable tool for studying particle acceleration and diffusion within our galaxy. Constructing a coherent physical picture of particle accelerators such as supernova remnants, pulsar wind nebulae, and star-forming regions requires the ability to detect extended regions of gamma ray emission, to analyze small-scale spatial variation within these regions, and to synthesize data from multiple observatories across multiple wavebands. Imaging atmospheric Cherenkov telescopes (IACTs) provide fine angular resolution (<0.1$^\circ$) for gamma rays above 100 GeV. However, their limited fields of view typically make detection of extended sources challenging. Maximum likelihood methods are well-suited to simultaneous analysis of multiple fields with overlapping sources and to combining data from multiple gamma ray observatories. Such methods also offer an alternative approach to estimating the IACT cosmic ray background and consequently an enhanced sensitivity to sources that may be as large as the telescope field of view. We report here on the current status and performance of a maximum likelihood technique for the IACT VERITAS.


title: Epsilon factors of symplectic type characters in the wild case

abstract: By work of John Tate we can associate an epsilon factor with every multiplicative character of a local field. In this paper we determine the explicit signs of the epsilon factors for symplectic type characters of $K^\times$, where $K/F$ is a wildly ramified quadratic extension of a non-Archimedean local field $F$ of characteristic zero.


title: Large-Scale-Exploit of GitHub Repository Metadata and Preventive Measures

abstract: When working with Git, a popular version-control system, email addresses are part of the metadata for each individual commit. When those commits are pushed to remote hosting services like GitHub, those email addresses become visible not only to fellow developers, but also to malicious actors aiming to exploit them. As a part of our research we created a tool that leverages the publicly available GitHub API to collect user data. Analysis of this data not only gives access to millions of email addresses in very little time, but is also powerful and dense enough to create targeted phishing attacks posing a great threat to all GitHub users and their private, potentially sensitive data. Even worse, existing countermeasures fail to effectively protect against such exploits. As a consequence and main conclusion of this paper, we suggest multiple preventive measures that should be implemented as soon as possible. We also consider it the duty of both companies like GitHub and well informed software engineers to inform fellow developers about the risk of exposing private email addresses in Git commits published publicly.


title: The generalization error of random features regression: Precise asymptotics and double descent curve

abstract: Deep learning methods operate in regimes that defy the traditional statistical mindset. Neural network architectures often contain more parameters than training samples, and are so rich that they can interpolate the observed labels, even if the latter are replaced by pure noise. Despite their huge complexity, the same architectures achieve small generalization error on real data. This phenomenon has been rationalized in terms of a so-called `double descent' curve. As the model complexity increases, the test error follows the usual U-shaped curve at the beginning, first decreasing and then peaking around the interpolation threshold (when the model achieves vanishing training error). However, it descends again as model complexity exceeds this threshold. The global minimum of the test error is found above the interpolation threshold, often in the extreme overparametrization regime in which the number of parameters is much larger than the number of samples. Far from being a peculiar property of deep neural networks, elements of this behavior have been demonstrated in much simpler settings, including linear regression with random covariates. In this paper we consider the problem of learning an unknown function over the $d$-dimensional sphere $\mathbb S^{d-1}$, from $n$ i.i.d. samples $(\boldsymbol x_i, y_i)\in \mathbb S^{d-1} \times \mathbb R$, $i\le n$. We perform ridge regression on $N$ random features of the form $\sigma(\boldsymbol w_a^{\mathsf T} \boldsymbol x)$, $a\le N$. This can be equivalently described as a two-layers neural network with random first-layer weights. We compute the precise asymptotics of the test error, in the limit $N,n,d\to \infty$ with $N/d$ and $n/d$ fixed. This provides the first analytically tractable model that captures all the features of the double descent phenomenon without assuming ad hoc misspecification structures.


title: The NANOGrav Program for Gravitational Waves and Fundamental Physics

abstract: We describe the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) and its efforts to directly detect and study gravitational waves and other synergistic physics and astrophysics using radio timing observations of millisecond pulsars.


title: Sequential Computer Experimental Design for Estimating an Extreme Probability or Quantile

abstract: A computer code can simulate a system's propagation of variation from random inputs to output measures of quality. Our aim here is to estimate a critical output tail probability or quantile without a large Monte Carlo experiment. Instead, we build a statistical surrogate for the input-output relationship with a modest number of evaluations and then sequentially add further runs, guided by a criterion to improve the estimate. We compare two criteria in the literature. Moreover, we investigate two practical questions: how to design the initial code runs and how to model the input distribution. Hence, we close the gap between the theory of sequential design and its application.


title: Moments of the weighted Cantor measures

abstract: Based on the seminal work of Hutchinson, we investigate properties of {\em $\alpha$-weighted Cantor measures} whose support is a fractal contained in the unit interval. Here, $\alpha$ is a vector of nonnegative weights summing to $1$, and the corresponding weighted Cantor measure $\mu^\alpha$ is the unique Borel probability measure on $[0,1]$ satisfying $ \mu^\alpha(E) = \sum_{ n=0 }^{N-1} \alpha_n\mu^\alpha( \varphi_n^{-1}(E) )$ where $\varphi_n: x\mapsto (x+n)/N$. In Sections 1 and 2 we examine several general properties of the measure $\mu^\alpha$ and the associated Legendre polynomials in $L_{\mu^\alpha}^2[0,1]$. In Section 3, we (1) compute the Laplacian and moment generating function of $\mu^\alpha$, (2) characterize precisely when the moments $I_m = \int_{[0,1]}x^m\,d\mu^\alpha$ exhibit either polynomial or exponential decay, and (3) describe an algorithm which estimates the first $m$ moments within uniform error $\varepsilon$ in $O( (\log\log(1/\varepsilon))\cdot m\log m )$. We also state analogous results in the natural case where $\alpha$ is {\em palindromic} for the measure $\nu^{\alpha}$ attained by shifting $\mu^{\alpha}$ to $[-1/2,1/2]$.


title: Desalination due to Electrical Image Forces

abstract: It will be shown that for a solution of salt dissolved in water in contact with a metallic wall, the concentration of salt ions (both positive and negative) within a few Angstroms of the wall can be large enough to exceed the solubility limit of the salt, as a result of electrical image charge forces. In addition, since the dielectric constant of water increases from 2.1 at the wall to 81 at about a nanometer from a solid wall, there will be an attractive image potential near the plane on which this increase of the dielectric constant occurs. The possible existence of these image potentials suggests that the salt can be removed from the water by making salt water flow between an array of parallel solid plates..


title: An Introduction to Stellarators: From magnetic fields to symmetries and optimization

abstract: In this self-contained document, we aim to present the basic theoretical building blocks to understand modeling of stellarator magnetic fields, some of the challenges associated with modeling, and optimization for designing stellarators. As often as possible, the ideas will be presented using equations and pictures, and references to other relevant introductory material will be included. This document is accessible to those who may not have a physics background but are interested in applications of mathematical and computational tools to stellarator research.


title: Placing quantified variants of 3-SAT and Not-All-Equal 3-SAT in the polynomial hierarchy

abstract: The complexity of variants of 3-SAT and Not-All-Equal 3-SAT is well studied. However, in contrast, very little is known about the complexity of the problems' quantified counterparts. In the first part of this paper, we show that $\forall \exists$ 3-SAT is $\Pi_2^P$-complete even if (1) each variable appears exactly twice unnegated and exactly twice negated, (2) each clause is a disjunction of exactly three distinct variables, and (3) the number of universal variables is equal to the number of existential variables. Furthermore, we show that the problem remains $\Pi_2^P$-complete if (1a) each universal variable appears exactly once unnegated and exactly once negated, (1b) each existential variable appears exactly twice unnegated and exactly twice negated, and (2) and (3) remain unchanged. On the other hand, the problem becomes NP-complete for certain variants in which each universal variable appears exactly once. In the second part of the paper, we establish $\Pi_2^P$-completeness for $\forall \exists$ Not-All-Equal 3-SAT even if (1') the Boolean formula is linear and monotone, (2') each universal variable appears exactly once and each existential variable appears exactly three times, and (3') each clause is a disjunction of exactly three distinct variables that contains at most one universal variable. On the positive side, we uncover variants of $\forall \exists$ Not-All-Equal 3-SAT that are co-NP-complete or solvable in polynomial time.


title: Metallicity has followed local gravitational potential of galaxies since z=3

abstract: The MZ relation between stellar mass (M*) and metallicity (Z) of nearby galaxies has been described as both a global and local property, i.e. valid also on sub-galaxy scales. Here we show that Z has remained a local property, following the gravitational potential, since z=3. In absorption the MZ relation has been well studied, and was in place already at z=5.1. A recent absorption study of GRB galaxies revealed a close match to Damped Ly{\alpha} (DLA) galaxies, surprising due to their vastly different impact parameters and leading the authors to suggest that local metallicity follows the local gravitational potential. In this paper we formulate an observational test of this hypothesis. The test, in essence, forms a prediction that the velocity dispersion of the absorbing gas in galaxy halos, normalized by the central velocity dispersion, must follow a steep log scale slope of -0.015 dex/kpc as a function of impact parameter out to at least 20-30 kpc. We then compile an archival data and literature based sample of galaxies seen in both emission and absorption suitable for the test, and find that current data confirm the hypothesis out to 40-60 kpc. In addition we show that the distribution of the velocity offsets between z em and z abs favours a model where DLA systems are composed of individual sub-clouds distributed along the entire path through the halo, and disfavours a model where they are one single cloud with a bulk motion and internal sub-structure.


title: Exploring galaxy dark matter halos across redshifts with strong quasar absorbers

abstract: Quasar lines of sight intersect intervening galaxy discs or circum-galactic environments at random impact parameters and potential well depths. Absorption line velocity widths ($\Delta v_{90}$) are known to scale with host galaxy stellar masses, and inversely with the projected separation from the quasar line of sight. Its dependence on stellar mass can be eliminated by normalising with the emission-line widths of the host galaxies, $\sigma_{em}$, so that absorbers with a range of $\Delta v_{90}$ values can be compared directly. Using a sample of DLA systems at 0.2 < z < 3.2 with spectroscopically confirmed host galaxies, we find that the velocity ratio $\Delta v_{90}/\sigma_{em}$ decreases with projected distances from the hosts. We compare the data with expectations of line-of-sight velocity dispersions derived for different dark matter halo mass distributions, and find that models with steeper radial dark matter profiles provide a better fit to the observations, although the scatter remains large. Gas outflows from the galaxies may cause an increased scatter, or scale radii of dark matter halo models may not be representative for the galaxies. We demonstrate by computing virial velocities, that metal-rich DLAs that belong to massive galaxy halos (M$_{halo} \approx 10^{12}$ M$_{\odot}$) mostly remain gravitationally bound to the halos.


title: Observation of Dynamical Fermionization

abstract: We observe dynamical fermionization, where the momentum distribution of a Tonks-Girardeau (T-G) gas of strongly interacting bosons in 1D evolves from bosonic to fermionic after its axial confinement is removed. The asymptotic momentum distribution after expansion in 1D is the distribution of rapidities, which are the conserved quantities associated with many-body integrable systems. Rapidities have not previously been measured in any interacting many-body quantum system. Our measurements agree well with T-G gas theory. We also study momentum evolution after the trap depth is suddenly changed to a new non-zero value. We observe the predicted bosonic-fermionic oscillations and see deviations from the theory outside of the T-G gas limit.


title: End-to-End Learning from Complex Multigraphs with Latent-Graph Convolutional Networks

abstract: We study the problem of end-to-end learning from complex multigraphs with potentially very large numbers of edges between two vertices, each edge labeled with rich information. Examples range from communication networks to flights between airports or financial transaction graphs. We propose Latent-Graph Convolutional Networks (L-GCNs), which propagate information from these complex edges to a latent adjacency tensor, after which further downstream tasks can be performed, such as node classification. We evaluate the performance of several variations of the model on two synthetic datasets simulating fraud in financial transaction networks, ensuring the model must make use of edge labels in order to achieve good classification performance. We find that allowing for nonlinear interactions on a per-neighbor basis boosts performance significantly, while showing promising results in an inductive setting. Finally, we demonstrate the use of L-GCNs on real-world data in the form of an urban transportation network.


title: Systematization of Knowledge and Implementation: Short Identity-Based Signatures

abstract: Identity-Based signature schemes are gaining a lot of popularity every day. Over the last decade, there has been a lot of schemes that have been proposed. Several libraries are there that implement identity-based cryptosystems that include identity-based signature schemes like the JPBC library which is written in Java and the charm-crypto library written in python. However, these libraries do not contain all of the popular schemes, rather the JPBC library contains only one identity-based signature scheme and the charm-crypto contains three. Furthermore, the implemented schemes are designed to work on one particular pairing curve. In pairing-based cryptosystems, even for a given signature scheme, the size of the signature and the performance i.e. the time to sign and verify depends on the chosen pairing curve. There are many applications in which the signature size is of more importance than the performance and similarly other applications where the performance is of more importance than signature size. In this work, we describe the popular signature schemes and their implementation using the JPBC library and describe how different pairing curves affect the signature size and performance. We also provide two methods to further shorten the signature size which is not present in the libraries by default.


title: Non-reciprocity using quadrature-phase time-varying slab resonators

abstract: In this paper, it is shown that non-reciprocity can be observed in time-varying media without employing spatio-temporal modulated permittivities. We show that by using only two one dimensional Fabry-Perot slabs with time-periodic permittivities having quadrature phase difference, it is possible to achieve considerable non-reciprocity in transmission at the incidence frequency. To analyze such scenario,generalized transfer matrices are introduced to find the wave amplitudes of all harmonics in all space. The results are verified by in-house FDTD simulations. Moreover, in order to have a simple model of such time-varying slab resonators, a time-perturbed coupled mode theory is developed for multiple resonances, and it is shown that the results obtained by this method and the analytical method are in excellent agreement.


title: Robust One-Bit Recovery via ReLU Generative Networks: Near-Optimal Statistical Rate and Global Landscape Analysis

abstract: We study the robust one-bit compressed sensing problem whose goal is to design an algorithm that faithfully recovers any sparse target vector $\theta_0\in\mathbb{R}^d$ \textit{uniformly} via $m$ quantized noisy measurements. Specifically, we consider a new framework for this problem where the sparsity is implicitly enforced via mapping a low dimensional representation $x_0 \in \mathbb{R}^k$ through a known $n$-layer ReLU generative network $G:\mathbb{R}^k\rightarrow\mathbb{R}^d$ such that $\theta_0 = G(x_0)$. Such a framework poses low-dimensional priors on $\theta_0$ without a known sparsity basis. We propose to recover the target $G(x_0)$ solving an unconstrained empirical risk minimization (ERM). Under a weak \textit{sub-exponential measurement assumption}, we establish a joint statistical and computational analysis. In particular, we prove that the ERM estimator in this new framework achieves a statistical rate of $m=\widetilde{\mathcal{O}}(kn \log d /\varepsilon^2)$ recovering any $G(x_0)$ uniformly up to an error $\varepsilon$. When the network is shallow (i.e., $n$ is small), we show this rate matches the information-theoretic lower bound up to logarithm factors of $\varepsilon^{-1}$. From the lens of computation, we prove that under proper conditions on the network weights, our proposed empirical risk, despite non-convexity, has no stationary point outside of small neighborhoods around the true representation $x_0$ and its negative multiple; furthermore, we show that the global minimizer of the empirical risk stays within the neighborhood around $x_0$ rather than its negative multiple under further assumptions on the network weights.


title: Search for Very-High-Energy (E $>$ 100 GeV) Emission from Geminga Supernova by VERITAS

abstract: Geminga is a nearby (250 pc) middle-aged (spin-down time scale ~12,000 years) pulsar associated with a supernova remnant. Geminga has been a prime candidate for the origin of the unexpectedly high flux of cosmic-ray positrons above 10 GeV detected at Earth. Extended TeV gamma-ray emission from a 2-degree region around the Geminga pulsar was detected by the HAWC observatory, thus suggesting efficient, high-energy leptonic acceleration. Fermi-LAT observations show that the density of GeV leptons in the TeV nebula is lower than predicted by single zone and two zone diffusion models constrained with the HAWC measurements. However, the energy gap between Fermi-LAT and HAWC (~500 GeV to ~1 TeV) remains under-examined. The VERITAS gamma-ray observatory is sensitive in the energy range from 100 GeV to greater than 30 TeV, filling the gap between Fermi-LAT and HAWC. Therefore, VERITAS measurements potentially provide missing information. VERITAS has observed Geminga for 93 hours since 2009 including 28 hours in the 2018/2019 season. However, the standard VERITAS data analysis techniques have insufficient sensitivity to sources extended at the scale of the HAWC detection, due to difficulties with background estimation. We developed the Matched Runs Method (MRM) for VERITAS analysis of spatially extended sources. MRM has been demonstrated to be an effective technique by applying it to archival VERITAS data, and we are currently applying it to the Geminga observations. Here we present the summary of the MRM.


title: Laboratory biases hinder Eco-Evo-Devo integration: hints from the microworld

abstract: How specific environmental contexts contribute to the robustness and variation of developmental trajectories and evolutionary transitions is a central point in Eco-Evo-Devo. However, the articulation of ecological, evolutionary and developmental processes into integrative frameworks has been elusive, partly because standard experimental designs neglect or oversimplify ecologically meaningful contexts. Microbial models are useful to expose and discuss two possible sources of bias associated with gene-centered experimental designs: the use of laboratory strains and laboratory environmental conditions. We illustrate our point by showing how contrasting developmental phenotypes in Myxococcus xanthus depend on the joint variation of temperature and substrate stiffness. Microorganismal development can provide key information for better understanding the role of environmental conditions in the evolution of developmental variation, and to overcome some of the limitations associated with current experimental approaches.


title: The relative $\mathcal{L}$-invariant of a compact $4$-manifold

abstract: In this paper, we introduce the relative $\mathcal{L}$-invariant $r\mathcal{L}(X)$ of a smooth, orientable, compact 4-manifold $X$ with boundary. This invariant is defined by measuring the lengths of certain paths in the cut complex of a trisection surface for $X$. This is motivated by the definition of the $\mathcal{L}$-invariant for smooth, orientable, closed 4-manifolds by Kirby and Thompson. We show that if $X$ is a rational homology ball, then $r\mathcal{L}(X)=0$ if and only if $X\cong B^4$. In order to better understand relative trisections, we also produce an algorithm to glue two relatively trisected 4-manifold by any Murasugi sum or plumbing in the boundary, and also prove that any two relative trisections of a given 4-manifold $X$ are related by interior stabilization, relative stabilization, and the relative double twist, which we introduce in this paper as a trisection version of one of Piergallini and Zuddas's moves on open book decompositions. Previously, it was only known (by Gay and Kirby) that relative trisections inducing equivalent open books on $X$ are related by interior stabilizations.


title: Uplift Modeling for Multiple Treatments with Cost Optimization

abstract: Uplift modeling is an emerging machine learning approach for estimating the treatment effect at an individual or subgroup level. It can be used for optimizing the performance of interventions such as marketing campaigns and product designs. Uplift modeling can be used to estimate which users are likely to benefit from a treatment and then prioritize delivering or promoting the preferred experience to those users. An important but so far neglected use case for uplift modeling is an experiment with multiple treatment groups that have different costs, such as for example when different communication channels and promotion types are tested simultaneously. In this paper, we extend standard uplift models to support multiple treatment groups with different costs. We evaluate the performance of the proposed models using both synthetic and real data. We also describe a production implementation of the approach.


title: Phenomenological Arrhenius Analyses in Plasmon-Enhanced Catalysis

abstract: A range of chemical reactions occurring on the surfaces of metal nanoparticles exhibit enhanced rates under plasmonic excitation. Recent analyses based on Arrhenius law fitting have argued in favor of a purely photothermal mechanism of enhancement and suggested the lack of an involvement of hot electrons. However, there is a caveat as shown here: under certain scenarios, it is practically impossible to distinguish between a photochemical (non-thermal) effect of plasmonic excitation and a purely photothermal one using a phenomenological Arrhenius fitting of the reaction rates alone.


title: Stability of explicit Runge-Kutta methods for high order finite element approximation of linear parabolic equations

abstract: We study the stability of explicit Runge-Kutta methods for high order Lagrangian finite element approximation of linear parabolic equations and establish bounds on the largest eigenvalue of the system matrix which determines the largest permissible time step. A bound expressed in terms of the ratio of the diagonal entries of the stiffness and mass matrices is shown to be tight within a small factor which depends only on the dimension and the choice of the reference element and basis functions but is independent of the mesh or the coefficients of the initial-boundary value problem under consideration. Another bound, which is less tight and expressed in terms of mesh geometry, depends only on the number of mesh elements and the alignment of the mesh with the diffusion matrix. The results provide an insight into how the interplay between the mesh geometry and the diffusion matrix affects the stability of explicit integration schemes when applied to a high order finite element approximation of linear parabolic equations on general nonuniform meshes.


title: Properties of the Interstellar Medium along Sight Lines to Nearby Planet Hosting Stars

abstract: We analyze high-resolution ultraviolet spectra of three nearby exoplanet host stars (HD 192310, HD 9826, and HD 206860) to study interstellar properties along their lines of sight and to search for the presence of astrospheric absorption. Using HST/STIS spectra of the Lyman-alpha, Mg II, and Fe II lines, we identify three interstellar velocity components in the lines of sight to each star. We can reliably assign eight of the nine components to partially ionized clouds found by Redfield & Linsky (2008) on the basis of the star's location in Galactic coordinates and agreement of measured radial velocities with velocities predicted from the cloud velocity vectors. None of the stars show blue-shifted absorption indicative of an astrosphere, implying that the stars are in regions of ionized interstellar gas. Coupling astrospheric and local interstellar medium measurements is necessary to evaluate the host star electromagnetic and particle flux, which have profound impacts on the atmospheres of their orbiting planets. We present a table of all known exoplanets located within 20 pc of the Sun listing their interstellar properties and velocities predicted from the local cloud velocity vectors.


title: Maximum Relevance and Minimum Redundancy Feature Selection Methods for a Marketing Machine Learning Platform

abstract: In machine learning applications for online product offerings and marketing strategies, there are often hundreds or thousands of features available to build such models. Feature selection is one essential method in such applications for multiple objectives: improving the prediction accuracy by eliminating irrelevant features, accelerating the model training and prediction speed, reducing the monitoring and maintenance workload for feature data pipeline, and providing better model interpretation and diagnosis capability. However, selecting an optimal feature subset from a large feature space is considered as an NP-complete problem. The mRMR (Minimum Redundancy and Maximum Relevance) feature selection framework solves this problem by selecting the relevant features while controlling for the redundancy within the selected features. This paper describes the approach to extend, evaluate, and implement the mRMR feature selection methods for classification problem in a marketing machine learning platform at Uber that automates creation and deployment of targeting and personalization models at scale. This study first extends the existing mRMR methods by introducing a non-linear feature redundancy measure and a model-based feature relevance measure. Then an extensive empirical evaluation is performed for eight different feature selection methods, using one synthetic dataset and three real-world marketing datasets at Uber to cover different use cases. Based on the empirical results, the selected mRMR method is implemented in production for the marketing machine learning platform. A description of the production implementation is provided and an online experiment deployed through the platform is discussed.


title: Resonant Machine Learning Based on Complex Growth Transform Dynamical Systems

abstract: Traditional energy-based learning models associate a single energy metric to each configuration of variables involved in the underlying optimization process. Such models associate the lowest energy state to the optimal configuration of variables under consideration, and are thus inherently dissipative. In this paper we propose an energy-efficient learning framework that exploits structural and functional similarities between a machine learning network and a general electrical network satisfying the Tellegen's theorem. In contrast to the standard energy-based models, the proposed formulation associates two energy components, namely, active and reactive energy to the network. This ensures that the network's active-power is dissipated only during the process of learning, whereas the reactive-power is maintained to be zero at all times. As a result, in steady-state, the learned parameters are stored and self-sustained by electrical resonance determined by the network's nodal inductances and capacitances. Based on this approach, this paper introduces three novel concepts: (a) A learning framework where the network's active-power dissipation is used as a regularization for a learning objective function that is subjected to zero total reactive-power constraint; (b) A dynamical system based on complex-domain, continuous-time growth transforms which optimizes the learning objective function and drives the network towards electrical resonance under steady-state operation; and (c) An annealing procedure that controls the trade-off between active-power dissipation and the speed of convergence. As a representative example, we show how the proposed framework can be used for designing resonant support vector machines (SVMs), where we show that the support-vectors correspond to an LC network with self-sustained oscillations.


title: Multi-Task Self-Supervised Learning for Disfluency Detection

abstract: Most existing approaches to disfluency detection heavily rely on human-annotated data, which is expensive to obtain in practice. To tackle the training data bottleneck, we investigate methods for combining multiple self-supervised tasks-i.e., supervised tasks where data can be collected without manual labeling. First, we construct large-scale pseudo training data by randomly adding or deleting words from unlabeled news data, and propose two self-supervised pre-training tasks: (i) tagging task to detect the added noisy words. (ii) sentence classification to distinguish original sentences from grammatically-incorrect sentences. We then combine these two tasks to jointly train a network. The pre-trained network is then fine-tuned using human-annotated disfluency detection training data. Experimental results on the commonly used English Switchboard test set show that our approach can achieve competitive performance compared to the previous systems (trained using the full dataset) by using less than 1% (1000 sentences) of the training data. Our method trained on the full dataset significantly outperforms previous methods, reducing the error by 21% on English Switchboard.


title: Macroscopic electron-hole distribution in silicon and cubic silicon carbide

abstract: Electron excitations at silicon and 3C-SiC surfaces caused by an intense femtosecond laser pulse can be calculated by solving the time-dependent density functional theory and the Maxwell's equation simultaneously. The energy absorption, carrier density, and electron-hole quasi-temperatures decrease exponentially in 100 nm from the surface. The electron and hole quasi-temperatures have finite values even at large distances from the surface because of a specific photo-absorption channel. Although the quasi-temperature in the silicone shows smooth exponential descrease, 3C-SiC shows stepwise decrease because of the change of concerning bands. The quasi-temperature depends not only on the excitation process, i.e., tunnel and multi-photon absorption, but also on the band structure significantly.


title: Residual Force Polytope: Admissible Task-Space Forces of Dynamic Trajectories

abstract: We propose a representation for the set of forces a robot can counteract using full system dynamics: the residual force polytope. Given the nominal torques required by a dynamic motion, this representation models the forces which can be sustained without interfering with that motion. The residual force polytope can be used to analyze and compare the set of admissible forces of different trajectories, but it can also be used to define metrics for solving optimization problems, such as in trajectory optimization or system design. We demonstrate how such a metric can be applied to trajectory optimization and compare it against other objective functions typically used. Our results show that the trajectories computed by optimizing objectives defined as functions of the residual force polytope are more robust to unknown external disturbances. The computational cost of these metrics is relatively high and not compatible with the short planning times required by online methods, but they are acceptable for planning motions offline.


title: A tractable case of the Turing automorphism problem: bi-uniformly $E_0$-invariant Cantor homeomorphisms

abstract: A function $F:2^\omega\to 2^\omega$ is an $E_0$-isomorphism if for all $x,y\in 2^\omega$, we have $xE_0y\iff f(x)E_0 f(y)$, where $xE_0y\iff(\exists a)(\forall n\ge b) x(n)=y(n)$. If such witnesses $a$ for $xE_0 y$ and for $f(x)E_0 f(y)$ depend on each other but not on $x$, $y$, then $F$ is called bi-uniform. It is shown that a homeomorphism of Cantor space which is a bi-uniform $E_0$-isomorphism can induce only the trivial automorphism of the Turing degrees.


title: Gordian complexes of knots and virtual knots given by region crossing changes and arc shift moves

abstract: Gordian complex of knots was defined by Hirasawa and Uchida as the simplicial complex whose vertices are knot isotopy classes in $\mathbb{S}^3$. Later Horiuchi and Ohyama defined Gordian complex of virtual knots using $v$-move and forbidden moves. In this paper we discuss Gordian complex of knots by region crossing change and Gordian complex of virtual knots by arc shift move. Arc shift move is a local move in the virtual knot diagram which results in reversing orientation locally between two consecutive crossings. We show the existence of an arbitrarily high dimensional simplex in both the Gordian complexes, i.e., by region crossing change and by the arc shift move. For any given knot (respectively, virtual knot) diagram we construct an infinite family of knots (respectively, virtual knots) such that any two distinct members of the family have distance one by region crossing change (respectively, arc shift move). We show that that the constructed virtual knots have the same affine index polynomial.


abstract: When working with decomposition-based algorithms, an appropriate set of weights might improve quality of the final solution. A set of uniformly distributed weights usually leads to well-distributed solutions on a Pareto front. However, there are two main difficulties with this approach. Firstly, it may fail depending on the problem geometry. Secondly, the population size becomes not flexible as the number of objectives increases. In this paper, we propose the MOEA/D with Uniformly Randomly Adaptive Weights (MOEA/DURAW) which uses the Uniformly Randomly method as an approach to subproblems generation, allowing a flexible population size even when working with many objective problems. During the evolutionary process, MOEA/D-URAW adds and removes subproblems as a function of the sparsity level of the population. Moreover, instead of requiring assumptions about the Pareto front shape, our method adapts its weights to the shape of the problem during the evolutionary process. Experimental results using WFG41-48 problem classes, with different Pareto front shapes, shows that the present method presents better or equal results in 77.5% of the problems evaluated from 2 to 6 objectives when compared with state-of-the-art methods in the literature.


title: Toward an accurate strongly-coupled many-body theory within the equation of motion framework

abstract: Non-perturbative aspects of the quantum many-body problem are revisited, discussed and advanced in the equation of motion framework. We compare the approach to the two-fermion response function truncated on the two-body level by the cluster expansion of the dynamical interaction kernel to the approach known as time blocking approximation. Such a comparison leads to an extended many-body theory with non-perturbative treatment of high-order configurations. The present implementation of the advanced theory introduces a new class of solutions for the response functions, which include explicitly beyond-mean-field correlations between up to six fermions. The novel approach, which includes configurations with two quasiparticles coupled to two phonons (2q$\otimes$2phonon), is discussed in detail for the particle-hole nuclear response and applied to medium-mass nuclei. The proposed developments are implemented numerically on the basis of the relativistic effective meson-nucleon Lagrangian and compared to the models confined by two-fermion and four-fermion configurations, which are considered as state-of-the-art for the response theory in nuclear structure calculations. The results obtained for the dipole response of $^{42,48}$Ca and $^{68}$Ni nuclei in comparison to available experimental data show that the higher configurations are necessary for a successful description of both gross and fine details of the spectra in both high-energy and low-energy sectors.


title: Secure Coded Cooperative Computation at the Heterogeneous Edge against Byzantine Attacks

abstract: Edge computing is emerging as a new paradigm to allow processing data at the edge of the network, where data is typically generated and collected, by exploiting multiple devices at the edge collectively. However, offloading tasks to other devices leaves the edge computing applications at the complete mercy of an attacker. One of the attacks, which is also the focus of this work, is Byzantine attacks, where one or more devices can corrupt the offloaded tasks. Furthermore, exploiting the potential of edge computing is challenging mainly due to the heterogeneous and time-varying nature of the devices at the edge. In this paper, we develop a secure coded cooperative computation mechanism (SC3) that provides both security and computation efficiency guarantees by gracefully combining homomorphic hash functions and coded cooperative computation. Homomorphic hash functions are used against Byzantine attacks and coded cooperative computation is used to improve computation efficiency when edge resources are heterogeneous and time-varying. Simulations results show that SC3 improves task completion delay significantly.


title: Robust estimation of the mean with bounded relative standard deviation

abstract: Many randomized approximation algorithms operate by giving a procedure for simulating a random variable $X$ which has mean $\mu$ equal to the target answer, and a relative standard deviation bounded above by a known constant $c$. Examples of this type of algorithm includes methods for approximating the number of satisfying assignments to 2-SAT or DNF, the volume of a convex body, and the partition function of a Gibbs distribution. Because the answer is usually exponentially large in the problem input size, it is typical to require an estimate $\hat \mu$ satisfy $\mathbb{P}(|\hat \mu/\mu - 1| > \epsilon) \leq \delta$, where $\epsilon$ and $\delta$ are user specified nonnegative parameters. The current best algorithm uses $\lceil 2c^2\epsilon^{-2}(1+\epsilon)^2 \ln(2/\delta) \rceil$ samples to achieve such an estimate. By modifying the algorithm in order to balance the tails, it is possible to improve this result to $\lceil 2(c^2\epsilon^{-2} + 1)/(1-\epsilon^2)\ln(2/\delta) \rceil$ samples. Aside from the theoretical improvement, we also consider how to best implement this algorithm in practice. Numerical experiments show the behavior of the estimator on distributions where the relative standard deviation is unknown or infinite.


title: HONEM: Learning Embedding for Higher Order Networks

abstract: Representation learning on networks offers a powerful alternative to the oft painstaking process of manual feature engineering, and as a result, has enjoyed considerable success in recent years. However, all the existing representation learning methods are based on the first-order network (FON), that is, the network that only captures the pairwise interactions between the nodes. As a result, these methods may fail to incorporate non-Markovian higher-order dependencies in the network. Thus, the embeddings that are generated may not accurately represent of the underlying phenomena in a network, resulting in inferior performance in different inductive or transductive learning tasks. To address this challenge, this paper presents HONEM, a higher-order network embedding method that captures the non-Markovian higher-order dependencies in a network. HONEM is specifically designed for the higher-order network structure (HON) and outperforms other state-of-the-art methods in node classification, network re-construction, link prediction, and visualization for networks that contain non-Markovian higher-order dependencies.


title: Analysis of Various Transformer Structures for High Frequency Isolation Applications

abstract: High frequency transformers are an integral part of power electronics devices and their parasitic parameters influence the performance and efficiency of the overall system. In this paper, transformer leakage inductances and parasitic capacitances are analyzed using finite element method (FEM) for different structures and windings arrangements of high frequency transformers. Also, magnetic field, electric field, and voltage distribution within the transformer is simulated and analyzed. Six different high frequency transformers with toroidal, EE, and UU cores with different windings are investigated for a 400(V)/400(V), 8 kVA transformer operating at 10 kHz. Additionally, interleaved windings for EE core are simulated and results compared with previous outcomes. Analysis results will help categorize each structure, based on its balance between leakage inductances and series parasitic capacitance. This information can later be used for optimal selection of transformers as a function of their operating frequency and enable designers to compromise between various parameters in different applications, especially new fast switches such as SiC and GaN.


title: SFSegNet: Parse Freehand Sketches using Deep Fully Convolutional Networks

abstract: Parsing sketches via semantic segmentation is attractive but challenging, because (i) free-hand drawings are abstract with large variances in depicting objects due to different drawing styles and skills; (ii) distorting lines drawn on the touchpad make sketches more difficult to be recognized; (iii) the high-performance image segmentation via deep learning technologies needs enormous annotated sketch datasets during the training stage. In this paper, we propose a Sketch-target deep FCN Segmentation Network(SFSegNet) for automatic free-hand sketch segmentation, labeling each sketch in a single object with multiple parts. SFSegNet has an end-to-end network process between the input sketches and the segmentation results, composed of 2 parts: (i) a modified deep Fully Convolutional Network(FCN) using a reweighting strategy to ignore background pixels and classify which part each pixel belongs to; (ii) affine transform encoders that attempt to canonicalize the shaking strokes. We train our network with the dataset that consists of 10,000 annotated sketches, to find an extensively applicable model to segment stokes semantically in one ground truth. Extensive experiments are carried out and segmentation results show that our method outperforms other state-of-the-art networks.


title: 15-nodal quartic surfaces. Part II: The automorphism group

abstract: We describe a set of generators and defining relations for the group of birational automorphisms of a general 15-nodal quartic surface in the complex projective 3-dimensional space.


title: Towards Knowledge-Based Recommender Dialog System

abstract: In this paper, we propose a novel end-to-end framework called KBRD, which stands for Knowledge-Based Recommender Dialog System. It integrates the recommender system and the dialog generation system. The dialog system can enhance the performance of the recommendation system by introducing knowledge-grounded information about users' preferences, and the recommender system can improve that of the dialog generation system by providing recommendation-aware vocabulary bias. Experimental results demonstrate that our proposed model has significant advantages over the baselines in both the evaluation of dialog generation and recommendation. A series of analyses show that the two systems can bring mutual benefits to each other, and the introduced knowledge contributes to both their performances.


title: Explicit Krein Resolvent Identities for Singular Sturm-Liouville Operators with Applications to Bessel Operators

abstract: We derive explicit Krein resolvent identities for generally singular Sturm-Liouville operators in terms of boundary condition bases and the Lagrange bracket. As an application of the resolvent identities obtained, we compute the trace of the resolvent difference of a pair of self-adjoint realizations of the Bessel expression $-d^2/dx^2+(\nu^2-(1/4))x^{-2}$ on $(0,\infty)$ for values of the parameter $\nu\in[0,1)$ and use the resulting trace formula to explicitly determine the spectral shift function for the pair.


title: Lifshitz Transition and Nontrivial H-Doping Effect in Cr-based Superconductor KCr$_3$As$_3$H$_x$

abstract: We report the first-principles study on the H-intercalated Cr-based superconductor KCr$_3$As$_3$H$_x$. Our results show a paramagnetic ground state for KCr$_3$As$_3$H. The electronic structure consists of two quasi-one-dimensional (Q1D) Fermi-surfaces and one 3D Fermi-surface which are mainly contributed by Cr-d$_{z^2}$, d$_{x^2-y^2}$ and d$_{xy}$ orbitals. The bare electron susceptibility shows a $\Gamma$-centered imaginary peak, indicating possible ferromagnetic spin fluctuations. Upon moderate hole doping, the system undergoes a Lifshitz transition, which may enhance the Q1D feature of the system. The Bader charge analysis and electron localization functions reveal a strong bonding nature of hydrogen in KCr$_3$As$_3$H, which results in a nontrivial electron doping in KCr$_3$As$_3$H.


title: Molecular envelope around the HII region RCW 120

abstract: The H II region RCW120 is a well-known object, which is often considered as a target to verify theoretical models of gas and dust dynamics in the interstellar medium. However, the exact geometry of RCW120 is still a matter of debate. In this work, we analyse observational data on molecular emission in RCW120 and show that 13CO(2-1) and C18O(2-1) lines are fitted by a 2D model representing a ring-like face-on structure. The changing of the C18O(3-2) line profile from double-peaked to single-peaked from the dense molecular Condensation 1 might be a signature of stalled expansion in this direction. In order to explain a self-absorption dip of the 13CO(2-1) and 13CO(3-2) lines, we suggest that RCW120 is surrounded by a diffuse molecular cloud, and find confirmation of this cloud on a map of interstellar extinction. Optically thick 13CO(2-1) emission and the infrared 8 um PAH band form a neutral envelope of the H II region resembling a ring, while the envelope breaks into separate clumps on images made with optically thin C18O(2-1) line and far-infrared dust emission.


title: Envelope Polyhedra

abstract: This paper presents an additional class of regular polyhedra--envelope polyhedra--made of regular polygons, where the arrangement of polygons (creating a single surface) around each vertex is identical; but dihedral angles between faces need not be identical, and some of the dihedral angles are 0 degrees (i.e., some polygons are placed back to back). For example, squares, 6 around a point, is produced by deleting the triangles from the rhombicuboctahedron, creating a hollow polyhedron of genus 7 with triangular holes connecting 18 interior and 18 exterior square faces. An empty cube missing its top and bottom faces becomes an envelope polyhedron, squares, 4 around a point, with a toroidal topology. This definition leads to many interesting finite and infinite multiply connected regular polygon networks, including one infinite network with squares, 14 around a point, and another with triangles, 18 around a point. These are introduced just over 50 years after my related paper on infinite spongelike pseudopolyhedra in American Mathematical Monthly (Gott, 1967).


title: Cascaded Lattice Boltzmann Modeling and Simulations of Three-Dimensional Non-Newtonian Fluid Flows

abstract: Non-Newtonian fluid flows, especially in three dimensions (3D), arise in numerous settings of interest to physics. Prior studies using the lattice Boltzmann method (LBM) of such flows have so far been limited to mainly to two dimensions and used less robust collision models. In this paper, we develop a new 3D cascaded LBM based on central moments and multiple relaxation times on a three-dimensional, nineteen velocity (D3Q19) lattice for simulation of generalized Newtonian (power law) fluid flows. The relaxation times of the second order moments are varied locally based on the local shear rate and parameterized by the consistency coefficient and the power law index of the nonlinear constitutive relation of the power law fluid. Numerical validation study of the 3D cascaded LBM for various benchmark problems, including the complex 3D non-Newtonian flow in a cubic cavity at different Reynolds numbers and power law index magnitudes encompassing shear thinning and shear thickening fluids, are presented. Furthermore, numerical stability comparisons of the proposed advanced LBM scheme against the LBM based on other collision models, such as the SRT model and MRT model based on raw moments, are made. Numerical results demonstrate the accuracy, second order grid convergence and significant improvements in stability of the 3D cascaded LBM for simulation of 3D non-Newtonian flows of power law fluids.


title: Multimode Lasing in Wave-Chaotic Semiconductor Microlasers

abstract: We investigate experimentally and theoretically the lasing behavior of dielectric microcavity lasers with chaotic ray dynamics. Experiments show multimode lasing for both D-shaped and stadium-shaped wave-chaotic cavities. Theoretical calculations also find multimode lasing for different shapes, sizes and refractive indices. While there are quantitative differences between the theoretical lasing spectra of the stadium and D-cavity, due to the presence of scarred modes with anomalously high quality factors, these differences decrease as the system size increases, and are also substantially reduced when the effects of surface roughness are taken into account. Lasing spectra calculations are based on Steady-State Ab Initio Laser Theory, and indicate that gain competition is not sufficient to result in single-mode lasing in these systems.


title: Treatment of Advanced Divertor Configurations in the Flux-Coordinate Independent turbulence code GRILLIX

abstract: Advanced divertor configurations modify the magnetic geometry of the diverter to achieve a combination of strong magnetic flux expansion, increased connection length and higher divertor volume - to improve detachment stability, neutral/impurity confinement and heat-channel broadening. In this paper, we discuss the modification of the Flux-Coordinate Independent (FCI) turbulence code GRILLIX to treat generalised magnetic geometry, to allow for the investigation of the effect of magnetic geometry on turbulent structures in the edge and SOL. The development of grids and parallel operators from numerically-defined magnetic equilibria is discussed, as is the application of boundary conditions via penalisation, with the finite-width method generalised to treat complex non-conformal boundaries. Initial testing of hyperbolic (advection) and parabolic (diffusion) test cases is presented for the Snowflake scenario.


title: How does Object-Oriented Code Refactoring Influence Software Quality? Research Landscape and Challenges

abstract: Context: Software refactoring aims to improve software quality and developer productivity. Numerous empirical studies investigating the impact of refactoring activities on software quality have been conducted over the last two decades. Objective: This study aims to perform a comprehensive systematic mapping study of existing empirical studies on evaluation of the effect of object-oriented code refactoring activities on software quality attributes. Method: We followed a multi-stage scrutinizing process to select 142 primary studies published till December 2017. The selected primary studies were further classified based on several aspects to answer the research questions defined for this work. In addition, we applied vote-counting approach to combine the empirical results and their analysis reported in primary studies. Results: The findings indicate that studies conducted in academic settings found more positive impact of refactoring on software quality than studies performed in industries. In general, refactoring activities caused all quality attributes to improve or degrade except for cohesion, complexity, inheritance, fault-proneness and power consumption attributes. Furthermore, individual refactoring activities have variable effects on most quality attributes explored in primary studies, indicating that refactoring does not always improve all quality attributes. Conclusions: This study points out several open issues which require further investigation, e.g., lack of industrial validation, lesser coverage of refactoring activities, limited tool support, etc.


title: Four-phonon scattering diminishes the optical phonon contribution and isotope effect to thermal conductivity of III-V semiconductors

abstract: Recent studies reveal that four-phonon scattering is generally important in determining thermal conductivities of solids. However, these studies have been focused on materials where thermal conductivity $\kappa$ is dominated by acoustic phonons, and the impact of four phonon scattering, although significant, is still generally smaller than three-phonon scattering. In this work, taking AlSb as example, we demonstrated that four-phonon scattering is even more critical to three-phonon scattering as it diminishes optical phonon thermal transport, and therefore significantly reduces the thermal conductivities of materials in which optical branches have long three-phonon lifetimes. Also, our calculations show that four-phonon scattering can play an extremely important role in weakening the isotope effect on $\kappa$. Specifically, four-phonon scattering reduces the room-temperature $\kappa$ of the isotopically pure and natural-occurring AlSb by 70$\%$ and 50$\%$, respectively. The reduction for isotopically pure and natural-occurring c-GaN is about 34$\%$ and 27$\%$, respectively. For isotopically-pure w-GaN, the reduction is about 13$\%$ at room temperature and 25$\%$ at 400 K. These results provided important guidance for experimentalists for achieving high thermal conductivities in III-V compounds for applications in semiconductor industry.


title: Improved measurements of the absolute branching fractions of the inclusive decays $D^{+(0)}\to\phi X$

abstract: By analyzing 2.93 fb$^{-1}$ of $e^+e^-$ annihilation data taken at the center-of-mass energy $\sqrt s=$ 3.773 GeV with the BESIII detector, we determine the branching fractions of the inclusive decays $D^+\to\phi X$ and $D^0\to\phi X$ to be $(1.135\pm0.034\pm0.031)\%$ and $(1.091\pm0.027\pm0.035)\%$, respectively, where $X$ denotes any possible particle combination. The first uncertainties are statistical and the second systematic. We also determine the branching fractions of the decays $D\to\phi X$ and their charge conjugate modes $\bar{D}\to\phi \bar{X}$ separately for the first time, and no significant CP asymmetry is observed.