title: Fractal Universality in Near-Threshold Magnetic Lanthanide Dimers

abstract: Ergodic quantum systems are often quite alike, whereas nonergodic, fractal systems are unique and display characteristic properties. We explore one of these fractal systems, weakly bound dysprosium lanthanide molecules, in an external magnetic field. As recently shown, colliding ultracold magnetic dysprosium atoms display a soft chaotic behavior with a small degree of disorder. We broaden this classification by investigating the generalized inverse participation ratio and fractal dimensions for large sets of molecular wave functions. Our exact close-coupling simulations reveal a dynamic phase transition from partially localized states to totally delocalized states and universality in its distribution by increasing the magnetic field strength to only a hundred Gauss (or 10 mT). Finally, we prove the existence of nonergodic delocalized phase in the system and explain the violation of ergodicity by strong coupling between near-threshold molecular states and the nearby continuum.


title: Segmentation of the prostate and organs at risk in male pelvic CT images using deep learning

abstract: Inter-and intra-observer variation in delineating regions of interest (ROIs) occurs because of differences in expertise level and preferences of the radiation oncologists. We evaluated the accuracy of a segmentation model using the U-Net structure to delineate the prostate, bladder, and rectum in male pelvic CT images. The dataset used for training and testing the model consisted of raw CT scan images of 85 prostate cancer patients. We designed a 2D U-Net model to directly learn a mapping function that converts a 2D CT grayscale image to its corresponding 2D OAR segmented image. Our network contains blocks of convolution 2D layers with variable kernel sizes, channel number, and activation functions. On the left side of the U-Net model, we used three 3x3 convolutions, each followed by a rectified linear unit (ReLu) (activation function), and one max pooling operation. On the right side of the U-Net model, we used a 2x2 transposed convolution and two 3x3 convolution networks followed by a ReLu activation function. The automatic segmentation using the U-Net generated an average dice similarity coefficient (DC) and standard deviation (SD) of the following: DC +- SD (0.88 +- 0.12), (0.95 +- 0.04), and (0.92 +- 0.06) for the prostate, bladder, and rectum, respectively. Furthermore, the mean of average surface Hausdorff distance (ASHD) and SD were 1.2 +- 0.9 mm, 1.08 +- 0.8 mm, and 0.8 +- 0.6 mm for the prostate, bladder, and rectum, respectively. Our proposed method, which employs the U-Net structure, is highly accurate and reproducible for automated ROI segmentation. This provides a foundation to improve automatic delineation of the boundaries between the target and surrounding normal soft tissues on a standard radiation therapy planning CT scan.


title: Bounding Multivariate Trigonometric Polynomials with Applications to Filter Bank Design

abstract: The extremal values of multivariate trigonometric polynomials are of interest in fields ranging from control theory to filter design, but finding the extremal values of such a polynomial is generally NP-Hard. In this paper, we develop simple and efficiently computable estimates of the extremal values of a multivariate trigonometric polynomial directly from its samples. We provide an upper bound on the modulus of a complex trigonometric polynomial, and develop upper and lower bounds for real trigonometric polynomials. For a univarite polynomial, these bounds are tighter than existing bounds, and the extension to multivariate polynomials is new. As an application, the lower bound provides a sufficient condition to certify global positivity of a real trigonometric polynomial. We use this condition to motivate a new algorithm for multi-dimensional, multirate, perfect reconstruction filter bank design. We demonstrate our algorithm by designing a 2D perfect reconstruction filter bank.


title: Volatility estimation in fractional Ornstein-Uhlenbeck models

abstract: In this article we study the asymptotic behaviour of the realized quadratic variation of a process $\int_{0}^{t}u_{s}dY_{s}^{(1)}$% , where $u$ is a $\beta$-H\"older continuous process with $\beta > 1-H$ and $Y_{t}^{(1)}=\int_{0}^{t}e^{-s}dB^{H}_{a_s}$, where $a_{t}=He^{\frac{t% }{H}} $ and $B^H$ is a fractional Brownian motion, is connected to the fractional Ornstein-Uhlenbeck process of the second kind. We prove almost sure convergence uniformly in time, and a stable weak convergence for the realized quadratic variation. As an application, we construct strongly consistent estimator for the integrated volatility parameter in a model driven by $Y^{(1)}$.


title: Revisiting Totally Positive Differential Systems: A Tutorial and New Results

abstract: A matrix is called totally nonnegative (TN) if all its minors are nonnegative, and totally positive (TP) if all its minors are positive. Multiplying a vector by a TN matrix does not increase the number of sign variations in the vector. In a largely forgotten paper, Schwarz (1970) considered matrices whose exponentials are TN or TP. He also analyzed the evolution of the number of sign changes in the vector solutions of the corresponding linear system. In a seemingly different line of research, Smillie (1984), Smith (1991), and others analyzed the stability of nonlinear tridiagonal cooperative systems by using the number of sign variations in the derivative vector as an integer-valued Lyapunov function. We provide a tutorial on these fascinating research topics and show that they are intimately related. This allows to derive generalizations of the results by Smillie (1984) and Smith (1991) while simplifying the proofs. It also opens the door to many new and interesting research directions.


title: A Learning Approach for Low-Complexity Optimization of Energy Efficiency in Multi-Carrier Wireless Networks

abstract: This paper proposes computationally efficient algorithms to maximize the energy efficiency in multi-carrier wireless interference networks, by a suitable allocation of the system radio resources, namely the transmit powers and subcarrier assignment. The problem is formulated as the maximization of the system Global Energy-Efficiency subject to both maximum power and minimum rate constraints. This leads to a challenging non-convex fractional problem, which is tackled through an interplay of fractional programming, learning, and game theory. The proposed algorithmic framework is provably convergent and has a complexity linear in both the number of users and subcarriers, whereas other available solutions can only guarantee a polynomial complexity in the number of users and subcarriers. Numerical results show that the proposed method performs similarly as other, more complex, algorithms.


title: ADMM for Multiaffine Constrained Optimization

abstract: We expand the scope of the alternating direction method of multipliers (ADMM). Specifically, we show that ADMM, when employed to solve problems with multiaffine constraints that satisfy certain verifiable assumptions, converges to the set of constrained stationary points if the penalty parameter in the augmented Lagrangian is sufficiently large. When the Kurdyka-\L{}ojasiewicz (K-\L{}) property holds, this is strengthened to convergence to a single constrained stationary point. Our analysis applies under assumptions that we have endeavored to make as weak as possible. It applies to problems that involve nonconvex and/or nonsmooth objective terms, in addition to the multiaffine constraints that can involve multiple (three or more) blocks of variables. To illustrate the applicability of our results, we describe examples including nonnegative matrix factorization, sparse learning, risk parity portfolio selection, nonconvex formulations of convex problems, and neural network training. In each case, our ADMM approach encounters only subproblems that have closed-form solutions.


title: Spin Seebeck effect and ballistic transport of quasi-acoustic magnons in room-temperature yttrium iron garnet films

abstract: We studied the transient behavior of the spin current generated by the longitudinal spin Seebeck effect (LSSE) in a set of platinum-coated yttrium iron garnet (YIG) films of different thicknesses. The LSSE was induced by means of pulsed microwave heating of the Pt layer and the spin currents were measured electrically using the inverse spin Hall effect in the same layer. We demonstrate that the time evolution of the LSSE is determined by the evolution of the thermal gradient triggering the flux of thermal magnons in the vicinity of the YIG/Pt interface. These magnons move ballistically within the YIG film with a constant group velocity, while their number decays exponentially within an effective propagation length. The ballistic flight of the magnons with energies above 20K is a result of their almost linear dispersion law, similar to that of acoustic phonons. By fitting the time-dependent LSSE signal for different film thicknesses varying by almost an order of magnitude, we found that the effective propagation length is practically independent of the YIG film thickness. We consider this fact as strong support of a ballistic transport scenario - the ballistic propagation of quasi-acoustic magnons in room temperature YIG.


title: All nearest neighbor calculation based on Delaunay graphs

abstract: When we have two data sets and want to find the nearest neighbour of each point in the first dataset among points in the second one, we need the all nearest neighbour operator. This is an operator in spatial databases that has many application in different fields such as GIS and VLSI circuit design. Existing algorithms for calculating this operator assume that there is no pre computation on these data sets. These algorithms has o(n*m*d) time complexity where n and m are the number of points in two data sets and d is the dimension of data points. With assumption of some pre computation on data sets algorithms with lower time complexity can be obtained. One of the most common pre computation on spatial data is Delaunay graphs. In the Delaunay graph of a data set each point is linked to its nearest neighbours. In this paper, we introduce an algorithm for computing the all nearest neighbour operator on spatial data sets based on their Delaunay graphs. The performance of this algorithm is compared with one of the best existing algorithms for computing ANN operator in terms of CPU time and the number of IOs. The experimental results show that this algorithm has better performance than the other.


title: Elliptic flow from Coulomb interaction and low density elastic scattering

abstract: In high energy heavy ion collisions and interacting cold atom systems, large elliptic flow anisotropies have been observed. For the large opacity ($\rho\sigma L\sim 10^{3}$) of the latter hydrodynamics is a natural consequence, but for the small opacity ($\rho\sigma L\sim 1$) of the former hydrodynamic description is questionable. To shed light onto the situation, we simulate the expansion of a low density Argon ion (or atom) system, initially trapped in an elliptical region, under the Coulomb interaction (or elastic scattering). Significant elliptic anisotropy is found in both cases, and the anisotropy depends on the initial spatial eccentricity and the density of the system. The results may provide insights into the physics of anisotropic flow in high energy heavy ion collisions and its role in the study of quantum chromodynamics.


title: Tunability: Importance of Hyperparameters of Machine Learning Algorithms

abstract: Modern supervised machine learning algorithms involve hyperparameters that have to be set before running them. Options for setting hyperparameters are default values from the software package, manual configuration by the user or configuring them for optimal predictive performance by a tuning procedure. The goal of this paper is two-fold. Firstly, we formalize the problem of tuning from a statistical point of view, define data-based defaults and suggest general measures quantifying the tunability of hyperparameters of algorithms. Secondly, we conduct a large-scale benchmarking study based on 38 datasets from the OpenML platform and six common machine learning algorithms. We apply our measures to assess the tunability of their parameters. Our results yield default values for hyperparameters and enable users to decide whether it is worth conducting a possibly time consuming tuning strategy, to focus on the most important hyperparameters and to chose adequate hyperparameter spaces for tuning.


title: Mapping the Invocation Structure of Online Political Interaction

abstract: The surge in political information, discourse, and interaction has been one of the most important developments in social media over the past several years. There is rich structure in the interaction among different viewpoints on the ideological spectrum. However, we still have only a limited analytical vocabulary for expressing the ways in which these viewpoints interact. In this paper, we develop network-based methods that operate on the ways in which users share content; we construct \emph{invocation graphs} on Web domains showing the extent to which pages from one domain are invoked by users to reply to posts containing pages from other domains. When we locate the domains on a political spectrum induced from the data, we obtain an embedded graph showing how these interaction links span different distances on the spectrum. The structure of this embedded network, and its evolution over time, helps us derive macro-level insights about how political interaction unfolded through 2016, leading up to the US Presidential election. In particular, we find that the domains invoked in replies spanned increasing distances on the spectrum over the months approaching the election, and that there was clear asymmetry between the left-to-right and right-to-left patterns of linkage.


title: The Beta-Bernoulli process and algebraic effects

abstract: In this paper we use the framework of algebraic effects from programming language theory to analyze the Beta-Bernoulli process, a standard building block in Bayesian models. Our analysis reveals the importance of abstract data types, and two types of program equations, called commutativity and discardability. We develop an equational theory of terms that use the Beta-Bernoulli process, and show that the theory is complete with respect to the measure-theoretic semantics, and also in the syntactic sense of Post. Our analysis has a potential for being generalized to other stochastic processes relevant to Bayesian modelling, yielding new understanding of these processes from the perspective of programming.


title: Two Families of Monogenic $S_4$ Quartic Number Fields

abstract: Consider the integral polynomials $f_{a,b}(x)=x^4+ax+b$ and $g_{c,d}(x)=x^4+cx^3+d$. Suppose $f_{a,b}(x)$ and $g_{c,d}(x)$ are irreducible, $b\mid a$, and the integers $b$, $d$, $256d-27c^4$, and $\dfrac{256b^3-27a^4}{\gcd(256b^3,27a^4)}$ are all square-free. Using the Montes algorithm, we show that a root of $f_{a,b}(x)$ or $g_{c,d}(x)$ defines a monogenic extension of $\mathbb{Q}$ and serves as a generator for a power integral basis of the ring of integers. In fact, we show monogeneity for slightly more general families. Further, we obtain lower bounds on the density of polynomials generating monogenic $S_4$ fields within the families $f_{b,b}(x)$ and $g_{1,d}(x)$.


title: Optimal Control of Autonomous Vehicles Approaching A Traffic Light

abstract: This paper devotes to the development of an optimal acceleration/speed profile for autonomous vehicles approaching a traffic light. The design objective is to achieve both short travel time and low energy consumption as well as avoid idling at a red light. This is achieved by taking full advantage of the traffic light information based on vehicle-to-infrastructure communication. The problem is modeled as a mixed integer programming, which is equivalently transformed into optimal control problems by relaxing the integer constraint. Then the direct adjoining approach is used to solve both free and fixed terminal time optimal control problems subject to state constraints. By an elaborate analysis, we are able to produce a real-time online analytical solution, distinguishing our method from most existing approaches based on numerical calculations. Extensive simulations are executed to compare the performance of autonomous vehicles under the proposed speed profile and human driving vehicles. The results show quantitatively the advantages of the proposed algorithm in terms of energy consumption and travel time.


title: Subsequential tightness of the maximum of two dimensional Ginzburg-Landau fields

abstract: We prove subsequential tightness of centered maxima of two-dimensional Ginzburg-Landau fields with bounded elliptic contrast.


title: Exoplanet Classification and Yield Estimates for Direct Imaging Missions

abstract: Future NASA concept missions that are currently under study, like Habitable Exoplanet Imaging Mission (HabEx) & Large Ultra-Violet Optical Infra Red (LUVOIR) Surveyor, would discover a large diversity of exoplanets. We propose here a classification scheme that distinguishes exoplanets into different categories based on their size and incident stellar flux, for the purpose of providing the expected number of exoplanets observed (yield) with direct imaging missions. The boundaries of this classification can be computed using the known chemical behavior of gases and condensates at different pressures and temperatures in a planetary atmosphere. In this study, we initially focus on condensation curves for sphalerite ZnS, H2O, CO2 and CH4. The order in which these species condense in a planetary atmosphere define the boundaries between different classes of planets. Broadly, the planets are divided into rocky (0.5 - 1.0RE), super-Earths (1.0- 1.75RE), sub-Neptunes (1.75-3.5RE), sub-Jovians (3.5 - 6.0RE) and Jovians (6-14.3RE) based on their planet sizes, and 'hot', 'warm' and 'cold' based on the incident stellar flux. We then calculate planet occurrence rates within these boundaries for different kinds of exoplanets, \eta_{planet}, using the community co-ordinated results of NASA's Exoplanet Program Analysis Group's Science Analysis Group-13 (SAG-13). These occurrence rate estimates are in turn used to estimate the expected exoplanet yields for direct imaging missions of different telescope diameter.


title: Points on nodal lines with given direction

abstract: We study of the directional distribution function of nodal lines for eigenfunctions of the Laplacian on a planar domain. This quantity counts the number of points where the normal to the nodal line points in a given direction. We give upper bounds for the flat torus, and compute the expected number for arithmetic random waves.


title: PySE: Software for Extracting Sources from Radio Images

abstract: PySE is a Python software package for finding and measuring sources in radio telescope images. The software was designed to detect sources in the LOFAR telescope images, but can be used with images from other radio telescopes as well. We introduce the LOFAR Telescope, the context within which PySE was developed, the design of PySE, and describe how it is used. Detailed experiments on the validation and testing of PySE are then presented, along with results of performance testing. We discuss some of the current issues with the algorithms implemented in PySE and their inter- action with LOFAR images, concluding with the current status of PySE and its future development.


title: Z\'ero-cycles sur les espaces homog\`enes et probl\`eme de Galois inverse

abstract: Let X be a smooth compactification of a homogeneous space of a linear algebraic group G over a number field k. We establish the conjecture of Colliot-Th\'el\`ene, Sansuc, Kato and Saito on the image of the Chow group of zero-cycles of X in the product of the same groups over all the completions of k. When G is semisimple and simply connected and the geometric stabiliser is finite and supersolvable, we show that rational points of X are dense in the Brauer-Manin set. For finite supersolvable groups, in particular for finite nilpotent groups, this yields a new proof of Shafarevich's theorem on the inverse Galois problem, and solves, at the same time, Grunwald's problem, for these groups. ----- Soit X une compactification lisse d'un espace homog\`ene d'un groupe alg\'ebrique lin\'eaire G sur un corps de nombres k. Nous \'etablissons la conjecture de Colliot-Th\'el\`ene, Sansuc, Kato et Saito sur l'image du groupe de Chow des z\'ero-cycles de X dans le produit des m\^emes groupes sur tous les compl\'et\'es de k. Lorsque G est semi-simple et simplement connexe et que le stabilisateur g\'eom\'etrique est fini et hyper-r\'esoluble, nous montrons que les points rationnels de X sont denses dans l'ensemble de Brauer-Manin. Pour les groupes finis hyper-r\'esolubles, en particulier pour les groupes finis nilpotents, cela donne une nouvelle preuve du th\'eor\`eme de Shafarevich sur le probl\`eme de Galois inverse et r\'esout en m\^eme temps, pour ces groupes, le probl\`eme de Grunwald.


title: Evolution of passive film behavior on mechanically polished NiTi shape memory alloy during forward martensitic transformation

abstract: The corrosion property of Nickel-titanium (NiTi) shape memory alloy is investigated during forward martensitic transformation between 40C and 0C. The Differential Scanning Calorimetry technique is used to find the forward and reverse martensitic transformation temperatures. NiTi shows the present of martensitic phase at near room temperature by cooling from 80C. The change in corrosion behavior is monitored by electrochemical Open Circuit Potential (OCP) during phase transformation. Cyclic potentiodynamic polarization method is used to identify the corrosion compound of NiTi. Corrosion rate at each temperature extracted from polarization data and activation energy for corrosion reaction in martensite and austenite phase are calculated based on an Arrhenius relation between corrosion current density and reciprocal temperature.


title: Small-scale dynamos in simulations of stratified turbulent convection

abstract: Small-scale dynamo action is often held responsible for the generation of quiet-Sun magnetic fields. We aim to determine the excitation conditions and saturation level of small-scale dynamos in non-rotating turbulent convection at low magnetic Prandtl numbers. We use high resolution direct numerical simulations of weakly stratified turbulent convection. We find that the critical magnetic Reynolds number for dynamo excitation increases as the magnetic Prandtl number is decreased, which might suggest that small-scale dynamo action is not automatically evident in bodies with small magnetic Prandtl numbers as the Sun. As a function of the magnetic Reynolds number (${\rm Rm}$), the growth rate of the dynamo is consistent with an ${\rm Rm}^{1/2}$ scaling. No evidence for a logarithmic increase of the growth rate with ${\rm Rm}$ is found.


title: Local Coupling Property for Markov Processes with Applications to L\'evy Processes

abstract: In this article, we define the new concept of local coupling property for Markov processes and study its relationship with distributional properties of the transition probability. In the special case of L\'evy processes we show that this property is equivalent to the absolute continuity of the transition probability and also provide a sufficient condition for it in terms of the L\'evy measure. Our result is stronger than existing results for absolute continuity of L\'evy distributions.


title: Artificial Noise Aided Secure Cognitive Beamforming for Cooperative MISO-NOMA Using SWIPT

abstract: Cognitive radio (CR) and non-orthogonal multiple access (NOMA) have been deemed two promising technologies due to their potential to achieve high spectral efficiency and massive connectivity. This paper studies a multiple-input single-output NOMA CR network relying on simultaneous wireless information and power transfer (SWIPT) conceived for supporting a massive population of power limited battery-driven devices. In contrast to most of the existing works, which use an ideally linear energy harvesting model, this study applies a more practical non-linear energy harvesting model. In order to improve the security of the primary network, an artificial-noise-aided cooperative jamming scheme is proposed. The artificial-noise-aided beamforming design problems are investigated subject to the practical secrecy rate and energy harvesting constraints. Specifically, the transmission power minimization problems are formulated under both perfect channel state information (CSI) and the bounded CSI error model. The problems formulated are non-convex, hence they are challenging to solve. A pair of algorithms either using semidefinite relaxation (SDR) or a cost function are proposed for solving these problems. Our simulation results show that the proposed cooperative jamming scheme succeeds in establishing secure communications and NOMA is capable of outperforming the conventional orthogonal multiple access in terms of its power efficiency. Finally, we demonstrate that the cost function algorithm outperforms the SDR-based algorithm.


title: Aggregative Coarsening for Multilevel Hypergraph Partitioning

abstract: Algorithms for many hypergraph problems, including partitioning, utilize multilevel frameworks to achieve a good trade-off between the performance and the quality of results. In this paper we introduce two novel aggregative coarsening schemes and incorporate them within state-of-the-art hypergraph partitioner Zoltan. Our coarsening schemes are inspired by the algebraic multigrid and stable matching approaches. We demonstrate the effectiveness of the developed schemes as a part of multilevel hypergraph partitioning framework on a wide range of problems.


title: An Expanded Local Variance Gamma model

abstract: The paper proposes an expanded version of the Local Variance Gamma model of Carr and Nadtochiy by adding drift to the governing underlying process. Still in this new model it is possible to derive an ordinary differential equation for the option price which plays a role of Dupire's equation for the standard local volatility model. It is shown how calibration of multiple smiles (the whole local volatility surface) can be done in such a case. Further, assuming the local variance to be a piecewise linear function of strike and piecewise constant function of time this ODE is solved in closed form in terms of Confluent hypergeometric functions. Calibration of the model to market smiles does not require solving any optimization problem and, in contrast, can be done term-by-term by solving a system of non-linear algebraic equations for each maturity, which is fast.


title: MILE: A Multi-Level Framework for Scalable Graph Embedding

abstract: Recently there has been a surge of interest in designing graph embedding methods. Few, if any, can scale to a large-sized graph with millions of nodes due to both computational complexity and memory requirements. In this paper, we relax this limitation by introducing the MultI-Level Embedding (MILE) framework -- a generic methodology allowing contemporary graph embedding methods to scale to large graphs. MILE repeatedly coarsens the graph into smaller ones using a hybrid matching technique to maintain the backbone structure of the graph. It then applies existing embedding methods on the coarsest graph and refines the embeddings to the original graph through a graph convolution neural network that it learns. The proposed MILE framework is agnostic to the underlying graph embedding techniques and can be applied to many existing graph embedding methods without modifying them. We employ our framework on several popular graph embedding techniques and conduct embedding for real-world graphs. Experimental results on five large-scale datasets demonstrate that MILE significantly boosts the speed (order of magnitude) of graph embedding while generating embeddings of better quality, for the task of node classification. MILE can comfortably scale to a graph with 9 million nodes and 40 million edges, on which existing methods run out of memory or take too long to compute on a modern workstation. Our code and data are publicly available with detailed instructions for adding new base embedding methods: \url{https://github.com/jiongqian/MILE}.


title: Status report of the ESCULAP project at Orsay: External injection of low energy electrons in a Plasma

abstract: The ESCULAP project aims at studying external injection of low energy (\SI{10}{MeV}) electrons in a plasma in the quasilinear regime. This facility will use the photo injector PHIL and the high power laser LASERIX. We will give a status report of the preliminary work on the facility and the status of the two machines. We will also present the results of simulations showing the expected performances of the facility.


title: Evolution of ferromagnetism in two-dimensional electron gas of LaTiO3/SrTiO3

abstract: Understanding, creating, and manipulating spin polarization of two-dimensional electron gases at complex oxide interfaces presents an experimental challenge. For example, despite almost a decade long research effort, the microscopic origin of ferromagnetism in LaAlO3/SrTiO3 heterojunction is still an open question. Here, by using a prototypical two-dimensional electron gas (2DEG) which emerges at the interface between band insulator SrTiO3 and antiferromagnetic Mott insulator LaTiO3 , the experiment reveals the evidence for magnetic phase separation in hole-doped Ti d1 t2g system resulting in spin-polarized 2DEG. The details of electronic and magnetic properties of the 2DEG were investigated by temperature-dependent d.c. transport, angle-dependent X-ray photoemission spectroscopy, and temperature-dependent magnetoresistance. The observation of clear hysteresis in magnetotransport at low magnetic fields implies spin-polarization from magnetic islands in the hole rich LaTiO3 near the interface. These findings emphasize the role of magnetic instabilities in doped Mott insulators thus providing another path for designing all-oxide structures relevant to spintronics applications.


title: Internal and surface waves in vibrofluidized granular materials: Role of cohesion

abstract: Wave phenomena in vibrofluidized dry and partially wet granular materials confined in a quasi-two-dimensional geometry are investigated with numerical simulations considering individual particles as hard spheres. Short ranged cohesive interactions arising from the formation of liquid bridges between adjacent particles are modeled by changing the velocity dependent coefficient of restitution. Such a change effectively suppresses the formation of surface waves, in agreement with previous experimental observations. The difference in pattern creation arises from the suppressed momentum transfer due to wetting and it can be quantitatively understood from an analysis of binary impacts.


title: Hall viscosity and nonlocal conductivity of "gapped graphene"

abstract: We calculate the Hall viscosity and the nonlocal (i.e., dependent on wave vector $\bm q$) Hall conductivity of "gapped graphene" (a non-topological insulator with two valleys) in the presence of a strong perpendicular magnetic field. Using the linear-response theory at zero temperature within the Dirac approximation for the Landau levels, we present analytical expressions for both valley and total Hall viscosity and conductivity up to $\bm q^2$ at all frequencies. Although the final formulas for total Hall viscosity and conductivity are similar to the ones previously obtained for gapless graphene, the derivation reveals a significant difference between the two systems. First of all, both the Hall viscosity and the Hall conductivity vanish when the Fermi level lies in the gap that separates the lowest Landau level in the conduction band from the highest Landau level in the valence band. It is only when the Fermi level {\it is not} in the gap that the familiar formulas of gapless graphene are recovered. Second, in the case of gapped graphene, it is not possible (at least within our present approach) to define a single-valley Hall viscosity: this quantity diverges with a strength proportional to the magnitude of the gap. It is only when both valleys are included that the diverging terms, having opposite signs in the two valleys, cancel out and the familiar result is recovered. In contrast to this, the nonlocal Hall conductivity is finite in each valley. These results indicate that the Hoyos-Son formula connecting the Hall viscosity to the coefficient of $q^2$ in the small-$q$ expansion of the $q$-dependent Hall conductivity cannot be applied to each valley, but only to the system as a whole. The problem of defining a "valley Hall viscosity" remains open.


title: Multiscale Planar Graph Generation

abstract: The study of network representations of physical, biological, and social phenomena can help us better understand the structural and functional dynamics of their networks and formulate predictive models of these phenomena. However, due to the scarcity of real-world network data owing to factors such as cost and effort required in collection of network data and the sensitivity of this data towards theft and misuse, engineers and researchers often rely on synthetic data for simulations, hypothesis testing, decision making, and algorithm engineering. An important characteristic of infrastructure networks such as roads, water distribution and other utility systems is that they can be embedded in a plane, therefore to simulate these system we need realistic networks which are also planar. While the currently-available synthetic network generators can model networks that exhibit realism, they do not guarantee or achieve planarity. Therefore, in this paper we present a flexible algorithm that can synthesize realistic networks that are planar. The method follows a multi-scale randomized editing approach generating a hierarchy of coarsened networks of a given planar graph and introducing edits at various levels in the hierarchy. The method preserves the structural properties with minimal bias including the planarity of the network, while introducing realistic variability at multiple scales.


title: Community detection using boundary nodes in complex networks

abstract: We propose a new local community detection algorithm that finds communities by identifying borderlines between them using boundary nodes. Our method performs label propagation for community detection, where nodes decide their labels based on the largest "benefit score" exhibited by their immediate neighbors as an attractor to their communities. We try different metrics and find that using the number of common neighbors as benefit scores leads to better decisions for community structure. The proposed algorithm has a local approach and focuses only on boundary nodes during iterations of label propagation, which eliminates unnecessary steps and shortens the overall execution time. It preserves small communities as well as big ones and can outperform other algorithms in terms of the quality of the identified communities, especially when the community structure is subtle. The algorithm has a distributed nature and can be used on large networks in a parallel fashion.


title: Viscous electron flow in mesoscopic two-dimensional electron gas

abstract: We report electrical and magneto transport measurements in mesoscopic size, two-dimensional (2D) electron gas in a GaAs quantum well. Remarkably, we find that the probe configuration and sample geometry strongly affects the temperature evolution of local resistance. We attribute all transport properties to the presence of hydrodynamic effects. Experimental results confirm the theoretically predicted significance of viscous flow in mesoscopic devices.


title: On tradeoffs between width- and fill-like graph parameters

abstract: In this work we consider two two-criteria optimization problems: given an input graph, the goal is to find its interval (or chordal) supergraph that minimizes the number of edges and its clique number simultaneously. For the interval supergraph, the problem can be restated as simultaneous minimization of the pathwidth $pw(G)$ and the profile $p(G)$ of the input graph $G$. We prove that for an arbitrary graph $G$ and an integer $t\in\{1,\ldots,pw(G)+1\}$, there exists an interval supergraph $G'$ of $G$ such that for its clique number it holds $\omega(G')\leq(1+\frac{2}{t})(pw(G)+1)$ and the number of its edges is bounded by $|E(G')|\leq(t+2)p(G)$. In other words, the pathwidth and the profile of a graph can be simultaneously minimized within the factors of $1+\frac{2}{t}$ (plus a small constant) and $t+2$, respectively. Note that for a fixed $t$, both upper bounds provide constant factor approximations. On the negative side, we show an example that proves that, for some graphs, there is no solution in which both parameters are optimal. In case of finding a chordal supergraph, the two corresponding graph parameters that reflect its clique size and number of edges are the treewidth and fill-in. We obtain that the treewidth and the fill-in problems are also `orthogonal' in the sense that for some graphs, a solution that minimizes one of those parameters cannot minimize the other. As a motivating example, we recall graph searching games which illustrates a need of simultaneous minimization of these pairs of graph parameters.


title: Simultaneous cores with restrictions and a question of Zaleski and Zeilberger

abstract: IMPORTANT NOTE: This paper is much rougher than I'd usually submit, and not entirely complete, though the main theorems and proofs should not be hard to follow. Given the ongoing strike at UK Universities it may be some time before I get to complete it to my satisfaction, and in the meantime people I've shared the preliminary draft with would like to be able to reference it. Hence I'm uploading it in its current form, and will update it later. The main new result of this paper is to count the number of (n,n+1)-core partitions with odd parts, answering a question of Zaleski and Zeilberger with bounty a charitable contribution to the OEIS. Along the way, we prove a general theorem giving a recurrence for (n,n+1)-core parts whose smallest part and consecutive part differences are restricted to lie in an arbitrary set M. This theorem unifies many known results about (n,n+1)-core partitions with restrictions. We end with discussions of extensions of the general theorem that keep track of the largest part, number of parts, and size of the partition, and about a few cases where the same methods work on more general simultaneous cores.


title: Partition-crossing hypergraphs

abstract: For a finite set $X$, we say that a set $H\subseteq X$ crosses a partition ${\cal P}=(X_1,\dots,X_k)$ of $X$ if $H$ intersects $\min (|H|,k)$ partition classes. If $|H|\geq k$, this means that $H$ meets all classes $X_i$, whilst for $|H|\leq k$ the elements of the crossing set $H$ belong to mutually distinct classes. A set system ${\cal H}$ crosses ${\cal P}$, if so does some $H\in {\cal H}$. The minimum number of $r$-element subsets, such that every $k$-partition of an $n$-element set $X$ is crossed by at least one of them, is denoted by $f(n,k,r)$. The problem of determining these minimum values for $k=r$ was raised and studied by several authors, first by Sterboul in 1973 [Proc. Colloq. Math. Soc. J. Bolyai, Vol. 10, Keszthely 1973, North-Holland/American Elsevier, 1975, pp. 1387--1404]. The present authors determined asymptotically tight estimates on $f(n,k,k)$ for every fixed $k$ as $n\to \infty$ [Graphs Combin., 25 (2009), 807--816]. Here we consider the more general problem for two parameters $k$ and $r$, and establish lower and upper bounds for $f(n,k,r)$. For various combinations of the three values $n,k,r$ we obtain asymptotically tight estimates, and also point out close connections of the function $f(n,k,r)$ to Tur\'an-type extremal problems on graphs and hypergraphs, or to balanced incomplete block designs.


title: A Resilient Image Matching Method with an Affine Invariant Feature Detector and Descriptor

abstract: Image feature matching is to seek, localize and identify the similarities across the images. The matched local features between different images can indicate the similarities of their content. Resilience of image feature matching to large view point changes is challenging for a lot of applications such as 3D object reconstruction, object recognition and navigation, etc, which need accurate and robust feature matching from quite different view points. In this paper we propose a novel image feature matching algorithm, integrating our previous proposed Affine Invariant Feature Detector (AIFD) and new proposed Affine Invariant Feature Descriptor (AIFDd). Both stages of this new proposed algorithm can provide sufficient resilience to view point changes. With systematic experiments, we can prove that the proposed method of feature detector and descriptor outperforms other state-of-the-art feature matching algorithms especially on view points robustness. It also performs well under other conditions such as the change of illumination, rotation and compression, etc.


title: Fluid-supported elastic sheet under compression: Multifold solutions

abstract: The properties of a hinged floating elastic sheet of finite length under compression are considered. Numerical continuation is used to compute spatially localized buckled states with many spatially localized folds. Both symmetric and antisymmetric states are computed and the corresponding bifurcation diagrams determined. Weakly nonlinear analysis is used to analyze the transition from periodic wrinkles to single fold and multifold states and to compute their energy. States with the same number of folds have energies that barely differ from each other and the energy gap decreases exponentially as localization increases. The stability properties of the different competing states are established and provide a basis for the study of fold interactions.


title: Numerical relativity in spherical coordinates with the Einstein Toolkit

abstract: Numerical relativity codes that do not make assumptions on spatial symmetries most commonly adopt Cartesian coordinates. While these coordinates have many attractive features, spherical coordinates are much better suited to take advantage of approximate symmetries in a number of astrophysical objects, including single stars, black holes and accretion disks. While the appearance of coordinate singularities often spoils numerical relativity simulations in spherical coordinates, especially in the absence of any symmetry assumptions, it has recently been demonstrated that these problems can be avoided if the coordinate singularities are handled analytically. This is possible with the help of a reference-metric version of the Baumgarte-Shapiro-Shibata-Nakamura formulation together with a proper rescaling of tensorial quantities. In this paper we report on an implementation of this formalism in the Einstein Toolkit. We adapt the Einstein Toolkit infrastructure, originally designed for Cartesian coordinates, to handle spherical coordinates, by providing appropriate boundary conditions at both inner and outer boundaries. We perform numerical simulations for a disturbed Kerr black hole, extract the gravitational wave signal, and demonstrate that the noise in these signals is orders of magnitude smaller when computed on spherical grids rather than Cartesian grids. With the public release of our new Einstein Toolkit thorns, our methods for numerical relativity in spherical coordinates will become available to the entire numerical relativity community.


title: Tinkertoys for the $E_8$ Theory

abstract: We construct the 4D N=2 SCFTs of class-S, which stem from the $E_8$ (2,0) theory. There are 49,836 isolated SCFTs which arise as 3-punctured spheres. Of these, 149 are "mixed" (contain free hypermultiplets accompanying the interacting SCFT) and 775 have enhanced global symmetries (beyond the manifest global symmetry associated to the punctures). Among the 49,836 3-punctured spheres we find (after removing any free hypermultiplets which may be present) 29 that are product SCFTs. Turning to 4-punctured spheres, we find 1,025,438 4D SCFTs arising as a gauging (with a simple gauge group) of a pair of 3-punctured spheres. We discuss a number of applications, including recovering several known 4D SCFTs. Our full set of results can be accessed on the Web at https://golem.ph.utexas.edu/class-S/E8/ .


title: Cognition and Reality

abstract: We discuss the two moments of human cognition, namely, apprehension (A), whereby a coherent perception emerges from the recruitment of neuronal groups, and judgment(B),that entails the comparison of two apprehensions acquired at different times, coded in a suitable language and retrieved by memory. (B) entails self-consciousness, in so far as the agent who expresses the judgment must be aware that the two apprehensions are submitted to his/her own scrutiny and that it is his/her task to extract a mutual relation. Since (B) lasts around 3 seconds, the semantic value of the pieces under comparison must be decided within that time. This implies a fast search of the memory contents. As a fact, exploring human subjects with sequences of simple words, we find evidence of a limited time window , corresponding to the memory retrieval of a linguistic item in order to match it with the next one in a text flow (be it literary, or musical,or figurative). While apprehension is globally explained as a Bayes inference, judgment tresults from an inverse Bayes inference. As a consequence, two hermeneutics emerge (called respectively circle and coil). The first one acts in a pre-assigned space of features. The second one provides the discovery of novel features, thus unveiling previously unknown aspects and hence representing the road to reality.


title: Gravitational Wave Searches with Pulsar Timing Arrays. I: Cancellation of Clock and Ephemeris Noises

abstract: We propose a data processing technique to cancel monopole and dipole noise sources (such as clock and ephemeris noises respectively) in pulsar timing array searches for gravitational radiation. These noises are the dominant sources of correlated timing fluctuations in the lower-part ($ \approx 10^{-9} - 10^{-8}$ Hz) of the gravitational wave band accessible by pulsar timing experiments. After deriving the expressions that reconstruct these noises from the timing data, we estimate the gravitational wave sensitivity of our proposed processing technique to single-source signals to be at least one order of magnitude higher than that achievable by directly processing the timing data from an equal-size array. Since arrays can generate pairs of clock and ephemeris-free timing combinations that are no longer affected by correlated noises, we implement with them the cross-correlation statistic to search for an isotropic stochastic gravitational wave background. We find the resulting optimal signal-to-noise ratio to be more than one order of magnitude larger than that obtainable by correlating pairs of timing data from arrays of equal size.


title: The impact of domain walls on the chiral magnetic effect in hot QCD matter

abstract: The Chiral Magnetic Effect (CME) -- the separation of positive and negative electric charges along the direction of the external magnetic field in quark-gluon plasma and other topologically non-trivial media -- is a consequence of the coupling of electrodynamics to the topological gluon field fluctuations that form metastable $CP$-odd domains. In phenomenological models it is usually assumed that the domains are uniform and the influence of the domain walls on the electric current flow is not essential. This paper challenges the latter assumption. A simple model consisting of a uniform spherical domain in a uniform time-dependent magnetic field is introduced and analytically solved. It is shown that (i) no electric current flows into or out of the domain, (ii) the charge separation current, viz. the total electric current flowing inside the domain in the external field direction, is a dissipative Ohm current, (iii) the CME effect can be produced either by the anomalous current or by the boundary conditions on the domain wall and (iv) the charge separation current oscillates in plasma long after the external field decays. These properties are qualitatively different from the CME in an infinite medium.


title: Searching for, and quantifying, non-convexity of functions

abstract: Convexity plays a prominent role in a number of problems, but practical considerations frequently give rise to non-convex functions. We suggest a method for determining convex regions, and also for assessing the lack of convexity in the other regions. The method relies on a specially constructed decomposition of symmetric matrices, such as the Hessian. We illustrate theoretical results using several examples, one of which analyses a problem arising in risk measurement and management in insurance and finance.


title: Bayesian shape modelling of cross-sectional geological data

abstract: Shape information is of great importance in many applications. For example, the oil-bearing capacity of sand bodies, the subterranean remnants of ancient rivers, is related to their cross-sectional shapes. The analysis of these shapes is therefore of some interest, but current classifications are simplistic and ad hoc. In this paper, we describe the first steps towards a coherent statistical analysis of these shapes by deriving the integrated likelihood for data shapes given class parameters. The result is of interest beyond this particular application.


title: Low energy paths for octahedral tilting in inorganic halide perovskites

abstract: Instabilities relating to cooperative octahedral tilting is common in materials with perovskite structures, and in particular in the sub class of halide perovskites. In this work, the energetics of octahedral tilting in the inorganic metal halide perovskites CsPbI$_3$ and CsSnI$_3$ are investigated using first-principles density functional theory calculations. Several low energy paths between symmetry equivalent variants of the stable orthorhombic (\textit{Pnma}) perovskite variant are identified and investigated. The results are in favor of the presence of dynamic disorder in the octahedral tilting phase transitions of inorganic halide perovskites. In particular, one specific type of path, corresponding to an out-of-phase "tilt switch", is found to have significantly lower energy barrier than the others, which indicates the existence of a temperature range where the dynamic fluctuations of the octahedra follow only this type of path. This could produce a time averaged structure corresponding to the intermediate tetragonal (\textit{P4/mbm}) structure observed in experiments. Deficiencies of the commonly employed simple one-dimensional "double well" potentials in describing the dynamics of the octahedra are pointed out and discussed.


title: Magnetic order and disorder in nanomagnets probed by superconducting vortices

abstract: We have studied two nanomagnet systems with strong (Co/Pd multilayers) and weak (NdCo alloy films) stray magnetic fields by probing the out-of-plane magnetic states with superconducting vortices. The hybrid samples are made of array of nanomagnets embedded in superconducting Nb thin films. The vortex motion detects relevant magnetic state features, since superconducting vortices are able to discriminate between different magnetic stray field strengths and directions. The usual matching effect between the superconducting vortex lattice and the periodic pinning array can be quenched by means of disorder magnetic potentials with strong stray fields at random. Ordered stray fields retrieve the matching effect and yield asymmetry and shift in the vortex dissipation signal. Furthermore vortices can discriminate the sizes of the nanomagnet magnetic domains, detecting magnetic domain sizes as small as 70 nm. In addition, we observe that the vortex cores play the crucial role instead of the supercurrents around the vortex.


title: Model-Based Identification and Control of a One-Legged Hopping Robot

abstract: Spring-mass models are well established tools for the analysis and control of legged locomotion. Among the alternatives, spring-loaded inverted pendulum (SLIP) model has shown to be a very accurate descriptor of animal locomotion. Despite its wide use, the SLIP model includes non-integrable stance dynamics that prevent analytical solutions for its equations of motion. Fortunately, there are approximate analytical solutions for different SLIP variants. However, the practicality of such approximations are mostly tested on simulation studies with a few notable exceptions. This thesis extends upon a recent approximation to a hip torque actuated dissipative SLIP (TD-SLIP) model that uses torque actuation to compensate for energy losses. Systematic experiments for careful assessment of the predictive performance of the approximate analytical solution is presented on a well-instrumented one-legged hopping robot which is revised to enhance compatibility and accuracy of the system. Electronic structure of the robot is modified according to TD-SLIP model such that robot uses a real-time operating system to increase processing speed. Using the parameters and results generated by the predictive performance of the approximate analytical solution, a model-based controller is designed and implemented on the robot platform to generate a stable closed-loop running behaviour on the one legged hoping robot platform. In addition, ground reaction forces during the stance phase on the experimental platform is investigated and compared with the human running and the traditional SLIP model data to understand if torque-actuated models approximate natural locomotion better than traditional model.


title: Broadband EPR Spectroscopy in Diverse Field Conditions Using Optically Detected Nitrogen-Vacancy Centers in Diamond

abstract: Paramagnetic magnetic resonance, a powerful technique for characterizing and identifying chemical targets, is increasingly used for imaging; however, low spin polarization at room temperature and moderate magnetic fields poses challenges for detecting small numbers of spins. In this work, we use fluorescence from nitrogen-vacancy (NV) centers in diamond to detect the electron paramagnetic resonance (EPR) spectrum of optically inactive target spins under various conditions of field magnitude and orientation. The protocol requires neither direct microwave manipulation of the NV spins nor spectral overlap between NV and target spin resonances, thus enabling broadband detection. This unexpected non-resonant coupling is attributable to a two-phonon process that relaxes NV spins proximate to the fluctuating dipole moment of the target spin, suggesting that the sensitivity is determined by the dipole-dipole coupling strength. This approach holds promise for sensitive EPR detection, particularly in settings where control over the diamond-crystal orientation is difficult. This is notably the case for biological sensing applications, where nanodiamonds are being pursued for their bright, stable fluorescence and biocompatibility.


title: On the Boundary Point Principle for divergence-type equations

abstract: We provide some versions of the Zaremba-Hopf-Oleinik boundary point lemma for general elliptic and parabolic equations in divergence form under the sharp requirements on the coefficients of equations and on the boundaries of domains.


title: Metric and geometric relaxations of self-contracted curves

abstract: Self-contractedness (or self-expandedness, depending on the orientation) is hereby extended in two natural ways giving rise, for any $\lambda\in\lbrack-1,1)$, to the metric notion of $\lambda $-curve and the (weaker) geometric notion of $\lambda$-cone property ($\lambda$-eel). In the Euclidean space $\mathbb{R}^{d}$ it is established that for $\lambda\in\lbrack-1,1/d)$ bounded $\lambda$-curves have finite length. For $\lambda\geq 1/\sqrt{5}$ it is always possible to construct bounded curves of infinite length in ${\mathbb{R}}^{3}$ which do satisfy the $\lambda $-cone property. This can never happen in ${\mathbb{R}}^{2}$ though: it is shown that all bounded planar curves with the $\lambda$-cone property have finite length.


title: Local and global methods in representations of Hecke algebras

abstract: This paper aims at developing a "local--global" approach for various types of finite dimensional algebras, especially those related to Hecke algebras. The eventual intention is to apply the methods and applications developed here to the cross-characteristic representation theory of finite groups of Lie type. The authors first review the notions of quasi-hereditary and stratified algebras over a Noetherian commutative ring. They prove that many global properties of these algebras hold if and only if they hold locally at every prime ideal. When the commutative ring is sufficiently good, it is often sufficient to check just the prime ideals of height at most one. These methods are applied to construct certain generalized q-Schur algebras, proving they are often quasi-hereditary (the "good" prime case) but always stratified. Finally, these results are used to prove a triangular decomposition matrix theorem for the modular representations of Hecke algebras at good primes. In the bad prime case, the generalized q-Schur algebras are at least stratified, and a block triangular analogue of the good prime case is proved, where the blocks correspond to Kazhdan-Lusztig cells.


title: Learning for Constrained Optimization: Identifying Optimal Active Constraint Sets

abstract: In many engineered systems, optimization is used for decision making at time-scales ranging from real-time operation to long-term planning. This process often involves solving similar optimization problems over and over again with slightly modified input parameters, often under tight latency requirements. We consider the problem of using the information available through this repeated solution process to directly learn a model of the optimal solution as a function of the input parameters, thus reducing the need to solve computationally expensive large-scale parametric programs in real time. Our proposed method is based on learning relevant sets of active constraints, from which the optimal solution can be obtained efficiently. Using active sets as features preserves information about the physics of the system, enables interpretable models, accounts for relevant safety constraints, and is easy to represent and encode. However, the total number of active sets is also very large, as it grows exponentially with system size. The key contribution of this paper is a streaming algorithm that learns the relevant active sets from training samples consisting of the input parameters and the corresponding optimal solution, without any assumptions on the problem structure. The algorithm comes with theoretical performance guarantees, and is known to converge fast for problem instances with a small number of relevant active sets. It can thus be used to establish the practicability of the learning method. Through extensive experiments on the Optimal Power Flow problem, we observe that often only a few active sets are relevant in practice, suggesting that the active sets is the appropriate level of abstraction for a learning algorithm to target.


title: Modeling Others using Oneself in Multi-Agent Reinforcement Learning

abstract: We consider the multi-agent reinforcement learning setting with imperfect information in which each agent is trying to maximize its own utility. The reward function depends on the hidden state (or goal) of both agents, so the agents must infer the other players' hidden goals from their observed behavior in order to solve the tasks. We propose a new approach for learning in these domains: Self Other-Modeling (SOM), in which an agent uses its own policy to predict the other agent's actions and update its belief of their hidden state in an online manner. We evaluate this approach on three different tasks and show that the agents are able to learn better policies using their estimate of the other players' hidden states, in both cooperative and adversarial settings.


title: Protocol-Dependence and State Variables in the Force-Moment Ensemble

abstract: Stress-based ensembles incorporating temperature-like variables have been proposed as a route to an equation of state for granular materials. To test the efficacy of this approach, we perform experiments on a two-dimensional photoelastic granular system under three loading conditions: uniaxial compression, biaxial compression, and simple shear. From the interparticle forces, we find that the distributions of the normal component of the coarse-grained force-moment tensor are exponential-tailed, while the deviatoric component is Gaussian-distributed. This implies that the correct stress-based statistical mechanics conserves both the force-moment tensor and the Maxwell-Cremona force-tiling area. As such, two variables of state arise: the tensorial angoricity ($\hat{\alpha}$) and a new temperature-like quantity associated with the force-tile area which we name {\it keramicity} ($\kappa$). Each quantity is observed to be inversely proportional to the global confining pressure; however only $\kappa$ exhibits the protocol-independence expected of a state variable, while $\hat{\alpha}$ behaves as a variable of process.


title: Selecting optimal subgroups for treatment using many covariates

abstract: We consider the problem of selecting the optimal subgroup to treat when data on covariates is available from a randomized trial or observational study. We distinguish between four different settings including (i) treatment selection when resources are constrained, (ii) treatment selection when resources are not constrained, (iii) treatment selection in the presence of side effects and costs, and (iv) treatment selection to maximize effect heterogeneity. We show that, in each of these cases, the optimal treatment selection rule involves treating those for whom the predicted mean difference in outcomes comparing those with versus without treatment, conditional on covariates, exceeds a certain threshold. The threshold varies across these four scenarios but the form of the optimal treatment selection rule does not. The results suggest a move away from traditional subgroup analysis for personalized medicine. New randomized trial designs are proposed so as to implement and make use of optimal treatment selection rules in health care practice.


title: Derived stacks in symplectic geometry

abstract: This is a survey paper on derived symplectic geometry, that will appear as a chapter contribution to the book "New Spaces for Mathematics and Physics", edited by Mathieu Anel and Gabriel Catren. Our goal is to explain how derived stacks can be useful for ordinary symplectic geometry, with an emphasis on examples coming from classical topological field theories. More precisely, we use classical Chern-Simons theory and moduli spaces of flat $G$-bundles and $G$-local systems as leading examples in our journey. We start in the introduction by reviewing various point-of-views on classical Chern--Simons theory and moduli of flat connections. In the main body of the Chapter we try to convince the reader how derived symplectic geometry (after Pantev-To\"en-Vaqui\'e-Vezzosi somehow reconciles all these different point-of-views.


title: Finite-temperature phase structure of SU(4) gauge theory with multiple fermion representations

abstract: We investigate the phase structure of SU(4) gauge theory with the gauge field simultaneously coupled to two flavors of fermion in the fundamental representation and two flavors of fermion in the two-index antisymmetric representation. We find that the theory has only two phases, a low-temperature phase with both species of fermion confined and chirally broken, and a high-temperature phase with both species of fermion deconfined and chirally restored. The single phase transition in the theory appears to be first order, in agreement with theoretical predictions.


title: The second moment of the Siegel transform in the space of symplectic lattices

abstract: Using results from spectral theory of Eisenstein series, we prove a formula for the second moment of the Siegel transform when averaged over the subspace of symplectic lattices. This generalizes the classical formula of Rogers for the second moment in the full space of unimodular lattices. Using this new formula we give very strong bounds for the discrepancy of the number of lattice points in an Borel set, which hold for generic symplectic lattices.


title: Optimizing over a Restricted Policy Class in Markov Decision Processes

abstract: We address the problem of finding an optimal policy in a Markov decision process under a restricted policy class defined by the convex hull of a set of base policies. This problem is of great interest in applications in which a number of reasonably good (or safe) policies are already known and we are only interested in optimizing in their convex hull. We show that this problem is NP-hard to solve exactly as well as to approximate to arbitrary accuracy. However, under a condition that is akin to the occupancy measures of the base policies having large overlap, we show that there exists an efficient algorithm that finds a policy that is almost as good as the best convex combination of the base policies. The running time of the proposed algorithm is linear in the number of states and polynomial in the number of base policies. In practice, we demonstrate an efficient implementation for large state problems. Compared to traditional policy gradient methods, the proposed approach has the advantage that, apart from the computation of occupancy measures of some base policies, the iterative method need not interact with the environment during the optimization process. This is especially important in complex systems where estimating the value of a policy can be a time consuming process.


title: Shaping Influence and Influencing Shaping: A Computational Red Teaming Trust-based Swarm Intelligence Model

abstract: Sociotechnical systems are complex systems, where nonlinear interaction among different players can obscure causal relationships. The absence of mechanisms to help us understand how to create a change in the system makes it hard to manage these systems. Influencing and shaping are social operators acting on sociotechnical systems to design a change. However, the two operators are usually discussed in an ad-hoc manner, without proper guiding models and metrics which assist in adopting these models successfully. Moreover, both social operators rely on accurate understanding of the concept of trust. Without such understanding, neither of these operators can create the required level to create a change in a desirable direction. In this paper, we define these concepts in a concise manner suitable for modelling the concepts and understanding their dynamics. We then introduce a model for influencing and shaping and use Computational Red Teaming principles to design and demonstrate how this model operates. We validate the results computationally through a simulation environment to show social influencing and shaping in an artificial society.


title: Square function estimates, BMO Dirichlet problem, and absolute continuity of harmonic measure on lower-dimensional sets

abstract: In the recent work [DFM1, DFM2] G. David, J. Feneuil, and the first author have launched a program devoted to an analogue of harmonic measure for lower-dimensional sets. A relevant class of partial differential equations, analogous to the class of elliptic PDEs in the classical context, is given by linear degenerate equations with the degeneracy suitably depending on the distance to the boundary. The present paper continues this line of research and focuses on the criteria of quantitative absolute continuity of the newly defined harmonic measure with respect to the Hausdorff measure, $\omega\in A_\infty(\sigma)$, in terms of solvability of boundary value problems. The authors establish, in particular, square function estimates and solvability of the Dirichlet problem in BMO for domains with lower-dimensional boundaries under the underlying assumption $\omega\in A_\infty(\sigma)$. More generally, it is proved that in all domains with Ahlfors regular boundaries the BMO solvability of the Dirichlet problem is necessary and sufficient for the absolute continuity of the harmonic measure.


title: Topological waves in fluids with odd viscosity

abstract: Fluids in which both time-reversal and parity are broken can display a dissipationless viscosity that is odd under each of these symmetries. Here, we show how this odd viscosity has a dramatic effect on topological sound waves in fluids, including the number and spatial profile of topological edge modes. Odd viscosity provides a short-distance cutoff that allows us to define a bulk topological invariant on a compact momentum space. As the sign of odd viscosity changes, a topological phase transition occurs without closing the bulk gap. Instead, at the transition point, the topological invariant becomes ill-defined because momentum space cannot be compactified. This mechanism is unique to continuum models and can describe fluids ranging from electronic to chiral active systems.


title: ABC Samplers

abstract: This Chapter, "ABC Samplers", is to appear in the forthcoming Handbook of Approximate Bayesian Computation (2018). It details the main ideas and algorithms used to sample from the ABC approximation to the posterior distribution, including methods based on rejection/importance sampling, MCMC and sequential Monte Carlo.


title: Byzantine-Resilient Distributed Observers for LTI Systems

abstract: Consider a linear time-invariant (LTI) dynamical system monitored by a network of sensors, modeled as nodes of an underlying directed communication graph. We study the problem of collaboratively estimating the state of the system when certain nodes are compromised by adversaries. Specifically, we consider a Byzantine adversary model, where a compromised node possesses complete knowledge of the system dynamics and the network, and can deviate arbitrarily from the rules of any prescribed algorithm. We first characterize certain fundamental limitations of any distributed state estimation algorithm in terms of the measurement and communication structure of the nodes. We then develop an attack-resilient, provably correct state estimation algorithm that admits a fully distributed implementation. To characterize feasible network topologies that guarantee success of our proposed technique, we introduce a notion of `strong-robustness' that captures both measurement and communication redundancy. Finally, by drawing connections to bootstrap percolation theory, we argue that given an LTI system and an associated sensor network, the `strong-robustness' property can be checked in polynomial time.


title: Hadronic Superpartners from Superconformal and Supersymmetric Algebra

abstract: Through the embedding of superconformal quantum mechanics into AdS space, it is possible to construct an effective supersymmetric QCD light-front Hamiltonian for hadrons, which includes a spin-spin interaction between the hadronic constituents. A specific breaking of conformal symmetry determines a unique effective quark-confining potential for light hadrons, as well as remarkable connections between the meson, baryon, and tetraquark spectra. The pion is massless in the chiral limit and has no supersymmetric partner. The excitation spectra of relativistic light-quark meson, baryon and tetraquark bound states lie on linear Regge trajectories with identical slopes in the radial and orbital quantum numbers. Although conformal symmetry is strongly broken by the heavy quark mass, the basic underlying supersymmetric mechanism, which transforms mesons to baryons (and baryons to tetraquarks) into each other, still holds and gives remarkable connections across the entire spectrum of light, heavy-light and double-heavy hadrons. Here we show that all the observed hadrons can be related through this effective supersymmetric QCD, and that it can be used to identify the structure of the new charmonium states.


title: On the Suitability of $L_p$-norms for Creating and Preventing Adversarial Examples

abstract: Much research effort has been devoted to better understanding adversarial examples, which are specially crafted inputs to machine-learning models that are perceptually similar to benign inputs, but are classified differently (i.e., misclassified). Both algorithms that create adversarial examples and strategies for defending against them typically use $L_p$-norms to measure the perceptual similarity between an adversarial input and its benign original. Prior work has already shown, however, that two images need not be close to each other as measured by an $L_p$-norm to be perceptually similar. In this work, we show that nearness according to an $L_p$-norm is not just unnecessary for perceptual similarity, but is also insufficient. Specifically, focusing on datasets (CIFAR10 and MNIST), $L_p$-norms, and thresholds used in prior work, we show through online user studies that "adversarial examples" that are closer to their benign counterparts than required by commonly used $L_p$-norm thresholds can nevertheless be perceptually different to humans from the corresponding benign examples. Namely, the perceptual distance between two images that are "near" each other according to an $L_p$-norm can be high enough that participants frequently classify the two images as representing different objects or digits. Combined with prior work, we thus demonstrate that nearness of inputs as measured by $L_p$-norms is neither necessary nor sufficient for perceptual similarity, which has implications for both creating and defending against adversarial examples. We propose and discuss alternative similarity metrics to stimulate future research in the area.


title: Resilient Leader-Follower Consensus to Arbitrary Reference Values

abstract: The problem of consensus in the presence of misbehaving agents has increasingly attracted attention in the literature. Prior results have established algorithms and graph structures for multi-agent networks which guarantee the consensus of normally behaving agents in the presence of a bounded number of misbehaving agents. The final consensus value is guaranteed to fall within the convex hull of initial agent states. However, the problem of consensus tracking considers consensus to arbitrary reference values which may not lie within such bounds. Conditions for consensus tracking in the presence of misbehaving agents has not been fully studied. This paper presents conditions for a network of agents using the W-MSR algorithm to achieve this objective.


title: Translating and Segmenting Multimodal Medical Volumes with Cycle- and Shape-Consistency Generative Adversarial Network

abstract: Synthesized medical images have several important applications, e.g., as an intermedium in cross-modality image registration and as supplementary training samples to boost the generalization capability of a classifier. Especially, synthesized computed tomography (CT) data can provide X-ray attenuation map for radiation therapy planning. In this work, we propose a generic cross-modality synthesis approach with the following targets: 1) synthesizing realistic looking 3D images using unpaired training data, 2) ensuring consistent anatomical structures, which could be changed by geometric distortion in cross-modality synthesis and 3) improving volume segmentation by using synthetic data for modalities with limited training samples. We show that these goals can be achieved with an end-to-end 3D convolutional neural network (CNN) composed of mutually-beneficial generators and segmentors for image synthesis and segmentation tasks. The generators are trained with an adversarial loss, a cycle-consistency loss, and also a shape-consistency loss, which is supervised by segmentors, to reduce the geometric distortion. From the segmentation view, the segmentors are boosted by synthetic data from generators in an online manner. Generators and segmentors prompt each other alternatively in an end-to-end training fashion. With extensive experiments on a dataset including a total of 4,496 CT and magnetic resonance imaging (MRI) cardiovascular volumes, we show both tasks are beneficial to each other and coupling these two tasks results in better performance than solving them exclusively.


title: Learning Binary Latent Variable Models: A Tensor Eigenpair Approach

abstract: Latent variable models with hidden binary units appear in various applications. Learning such models, in particular in the presence of noise, is a challenging computational problem. In this paper we propose a novel spectral approach to this problem, based on the eigenvectors of both the second order moment matrix and third order moment tensor of the observed data. We prove that under mild non-degeneracy conditions, our method consistently estimates the model parameters at the optimal parametric rate. Our tensor-based method generalizes previous orthogonal tensor decomposition approaches, where the hidden units were assumed to be either statistically independent or mutually exclusive. We illustrate the consistency of our method on simulated data and demonstrate its usefulness in learning a common model for population mixtures in genetics.


title: Event-Triggered Controller Synthesis for Dynamical Systems with Temporal Logic Constraints

abstract: In this work, we propose an event-triggered con- trol framework for dynamical systems with temporal logical constraints. Event-triggered control methodologies have proven to be very efficient in reducing sensing, communication and computation costs. When a continuous feedback control is re- placed with an event-triggered strategy, the corresponding state trajectories also differ. In a system with logical constraints, such small deviation in the trajectory might lead to unsatisfiability of the logical constraints. In this work, we develop an approach where we ensure that the event-triggered state trajectory is confined within an tube of the ideal trajectory associated with the continuous state feedback. At the same time, we will ensure satisfiability of the logical constraints as well. Furthermore, we show that the proposed method works for delayed systems as long as the delay is bounded by a certain quantity.


title: The normalized volume of a singularity is lower semicontinuous

abstract: We show that in any $\mathbb{Q}$-Gorenstein flat family of klt singularities, normalized volumes are lower semicontinuous with respect to the Zariski topology. A quick consequence is that smooth points have the largest normalized volume among all klt singularities. Using an alternative characterization of K-semistability developed by Li, Liu and Xu, we show that K-semistability is a very generic or empty condition in any $\mathbb{Q}$-Gorenstein flat family of log Fano pairs.


title: Semantic segmentation of trajectories with agent models

abstract: In many cases, such as trajectories clustering and classification, we often divide a trajectory into segments as preprocessing. In this paper, we propose a trajectory semantic segmentation method based on learned behavior models. In the proposed method, we learn some behavior models from video sequences. Next, using learned behavior models and a hidden Markov model, we segment a trajectory into semantic segments. Comparing with the Ramer-Douglas-Peucker algorithm, we show the effectiveness of the proposed method.


title: Computational Red Teaming in a Sudoku Solving Context: Neural Network Based Skill Representation and Acquisition

abstract: In this paper we provide an insight into the skill representation, where skill representation is seen as an essential part of the skill assessment stage in the Computational Red Teaming process. Skill representation is demonstrated in the context of Sudoku puzzle, for which the real human skills used in Sudoku solving, along with their acquisition, are represented computationally in a cognitively plausible manner, by using feed-forward neural networks with back-propagation, and supervised learning. The neural network based skills are then coupled with a hard-coded constraint propagation computational Sudoku solver, in which the solving sequence is kept hard-coded, and the skills are represented through neural networks. The paper demonstrates that the modified solver can achieve different levels of proficiency, depending on the amount of skills acquired through the neural networks. Results are encouraging for developing more complex skill and skill acquisition models usable in general frameworks related to the skill assessment aspect of Computational Red Teaming.


title: Cloth Manipulation Using Random-Forest-Based Imitation Learning

abstract: We present a novel approach for robust manipulation of high-DOF deformable objects such as cloth. Our approach uses a random forest-based controller that maps the observed visual features of the cloth to an optimal control action of the manipulator. The topological structure of this random forest-based controller is determined automatically based on the training data consisting visual features and optimal control actions. This enables us to integrate the overall process of training data classification and controller optimization into an imitation learning (IL) approach. Our approach enables learning of robust control policy for cloth manipulation with guarantees on convergence.We have evaluated our approach on different multi-task cloth manipulation benchmarks such as flattening, folding and twisting. In practice, our approach works well with different deformable features learned based on the specific task or deep learning. Moreover, our controller outperforms a simple or piecewise linear controller in terms of robustness to noise. In addition, our approach is easy to implement and does not require much parameter tuning.


title: Directional Statistics-based Deep Metric Learning for Image Classification and Retrieval

abstract: Deep distance metric learning (DDML), which is proposed to learn image similarity metrics in an end-to-end manner based on the convolution neural network, has achieved encouraging results in many computer vision tasks.$L2$-normalization in the embedding space has been used to improve the performance of several DDML methods. However, the commonly used Euclidean distance is no longer an accurate metric for $L2$-normalized embedding space, i.e., a hyper-sphere. Another challenge of current DDML methods is that their loss functions are usually based on rigid data formats, such as the triplet tuple. Thus, an extra process is needed to prepare data in specific formats. In addition, their losses are obtained from a limited number of samples, which leads to a lack of the global view of the embedding space. In this paper, we replace the Euclidean distance with the cosine similarity to better utilize the $L2$-normalization, which is able to attenuate the curse of dimensionality. More specifically, a novel loss function based on the von Mises-Fisher distribution is proposed to learn a compact hyper-spherical embedding space. Moreover, a new efficient learning algorithm is developed to better capture the global structure of the embedding space. Experiments for both classification and retrieval tasks on several standard datasets show that our method achieves state-of-the-art performance with a simpler training procedure. Furthermore, we demonstrate that, even with a small number of convolutional layers, our model can still obtain significantly better classification performance than the widely used softmax loss.


title: Killing boundary data for anti-de Sitter-like spacetimes

abstract: Given an initial-boundary value problem for an anti-de Sitter-like spacetime, we analyse conditions on the conformal boundary ensuring the existence of Killing vectors in the spacetime arising from this problem. This analysis makes use of a system of conformal wave equations describing the propagation of the Killing equation first considered by Paetz. We identify an obstruction tensor constructed from Killing vector candidate and the Cotton tensor of the conformal boundary whose vanishing is a necessary condition for the existence of Killing vectors in the spacetime. This obstruction tensor vanishes if the conformal boundary is conformally flat.


title: Thermal tuning capabilities of semiconductor metasurface resonators

abstract: Metasurfaces exploit the ability to engineer the optical phase, amplitude and polarization at subwavelength dimensions providing unprecedented control of light. The realization of the all dielectric approach to metasurfaces has led to the demonstration of extensive flat optical elements and functionalities with low losses. However, to reach their ultimate potential, metasurfaces must move beyond static operation and incorporate active tunability and reconfigurable functions. The central challenge is achieving large tunability in subwavelength resonator elements which require large optical effects in response to external stimuli. Here we study the thermal tunability of high-index silicon and germanium semiconductor resonators over a large temperature range. We demonstrate thermal tuning of Mie resonances due to the normal positive thermo-optic effect (dn/dT >0) over a wide infrared range. We show that at higher temperatures and long wavelengths the sign of the thermo-optic coefficient is reversed (dn/dT<0) culminating in a negative induced index due to thermal excitation of free carriers. We also demonstrate the tuning of high order Mie resonances by several linewidths with a temperature swing of {\Delta}T<100K. Finally, we exploit the larger thermo-optic coefficient at NIR wavelengths in Si metasurfaces to realize optical switching and tunable metafilters.


title: Polynomial Treedepth Bounds in Linear Colorings

abstract: Low-treedepth colorings are an important tool for algorithms that exploit structure in classes of bounded expansion; they guarantee subgraphs that use few colors have bounded treedepth. These colorings have an implicit tradeoff between the total number of colors used and the treedepth bound, and prior empirical work suggests that the former dominates the run time of existing algorithms in practice. We introduce $p$-linear colorings as an alternative to the commonly used $p$-centered colorings. They can be efficiently computed in bounded expansion classes and use at most as many colors as $p$-centered colorings. Although a set of $k<p$ colors from a $p$-centered coloring induces a subgraph of treedepth at most $k$, the same number of colors from a $p$-linear coloring may induce subgraphs of larger treedepth. We establish a polynomial upper bound on the treedepth in general graphs, and give tighter bounds in trees and interval graphs via constructive coloring algorithms. We also give a co-NP-completeness reduction for recognizing $p$-linear colorings and discuss ways to overcome this limitation in practice. This preprint extends results that appeared in [9]; for full proofs omitted from [9], see previous versions of this preprint.


title: High-resolution Observations of Low-luminosity Gigahertz-Peaked Spectrum and Compact Steep Spectrum Sources

abstract: We present Very Long Baseline Interferometry observations of a faint and low-luminosity ($L_{\rm 1.4 GHz} < 10^{27}~\mbox{W Hz}^{-1}$) Gigahertz-Peaked Spectrum (GPS) and Compact Steep Spectrum (CSS) sample. We select eight sources from deep radio observations that have radio spectra characteristic of a GPS or CSS source and an angular size of $\theta \lesssim 2$ arcsec, and detect six of them with the Australian Long Baseline Array. We determine their linear sizes, and model their radio spectra using Synchrotron Self Absorption (SSA) and Free Free Absorption (FFA) models. We derive statistical model ages, based on a fitted scaling relation, and spectral ages, based on the radio spectrum, which are generally consistent with the hypothesis that GPS and CSS sources are young and evolving. We resolve the morphology of one CSS source with a radio luminosity of $10^{25}~\mbox{W Hz}^{-1}$, and find what appear to be two hotspots spanning 1.7 kpc. We find that our sources follow the turnover-linear size relation, and that both homogenous SSA and an inhomogeneous FFA model can account for the spectra with observable turnovers. All but one of the FFA models do not require a spectral break to account for the radio spectrum, while all but one of the alternative SSA and power law models do require a spectral break to account for the radio spectrum. We conclude that our low-luminosity sample is similar to brighter samples in terms of their spectral shape, turnover frequencies, linear sizes, and ages, but cannot test for a difference in morphology.


title: Sufficient variable screening via directional regression with censored response

abstract: We in this paper propose a directional regression based approach for ultrahigh dimensional sufficient variable screening with censored responses. The new method is designed in a model-free manner and thus can be adapted to various complex model structures. Under some commonly used assumptions, we show that the proposed method enjoys the sure screening property when the dimension p diverges at an exponential rate of the sample size n. To improve the marginal screening method, the corresponding iterative screening algorithm and stability screening algorithm are further equipped. We demonstrate the effectiveness of the proposed method through simulation studies and a real data analysis.


title: Propagation of chaos for the Keller-Segel equation over bounded domains

abstract: In this paper we rigorously justify the propagation of chaos for the parabolic-elliptic Keller-Segel equation over bounded convex domains. The boundary condition under consideration is the no-flux condition. As intermediate steps, we establish the well-posedness of the associated stochastic equation as well as the well-posedness of the Keller-Segel equation for bounded weak solutions.


title: A Multi-Disciplinary Review of Knowledge Acquisition Methods: From Human to Autonomous Eliciting Agents

abstract: This paper offers a multi-disciplinary review of knowledge acquisition methods in human activity systems. The review captures the degree of involvement of various types of agencies in the knowledge acquisition process, and proposes a classification with three categories of methods: the human agent, the human-inspired agent, and the autonomous machine agent methods. In the first two categories, the acquisition of knowledge is seen as a cognitive task analysis exercise, while in the third category knowledge acquisition is treated as an autonomous knowledge-discovery endeavour. The motivation for this classification stems from the continuous change over time of the structure, meaning and purpose of human activity systems, which are seen as the factor that fuelled researchers' and practitioners' efforts in knowledge acquisition for more than a century. We show through this review that the KA field is increasingly active due to the higher and higher pace of change in human activity, and conclude by discussing the emergence of a fourth category of knowledge acquisition methods, which are based on red-teaming and co-evolution.


title: Single-View Food Portion Estimation: Learning Image-to-Energy Mappings Using Generative Adversarial Networks

abstract: Due to the growing concern of chronic diseases and other health problems related to diet, there is a need to develop accurate methods to estimate an individual's food and energy intake. Measuring accurate dietary intake is an open research problem. In particular, accurate food portion estimation is challenging since the process of food preparation and consumption impose large variations on food shapes and appearances. In this paper, we present a food portion estimation method to estimate food energy (kilocalories) from food images using Generative Adversarial Networks (GAN). We introduce the concept of an "energy distribution" for each food image. To train the GAN, we design a food image dataset based on ground truth food labels and segmentation masks for each food image as well as energy information associated with the food image. Our goal is to learn the mapping of the food image to the food energy. We can then estimate food energy based on the energy distribution. We show that an average energy estimation error rate of 10.89% can be obtained by learning the image-to-energy mapping.


title: Energy Levels, Lifetimes and Transition rates for P-like ions from Cr X to Zn XVI from large-scale Relativistic Multiconfiguration Calculations

abstract: The fully relativistic multiconfiguration Dirac--Hartree--Fock method is used to compute excitation energies and lifetimes for the 143 lowest states of the $3s^23p^3$, $3s3p^4$, $3s^23p^23d$, $3s3p^33d$, $3p^5$, $3s^23p3d^2$ configurations in P-like ions from Cr X to Zn XVI. Multipole (E1, M1, E2, M2) transition rates, line strengths, oscillator strengths, and branching fractions among these states are also given. Valence-valence and core-valence electron correlation effects are systematically accounted for using large basis function expansions. Computed excitation energies are compared with the NIST ASD and CHIANTI compiled values and previous calculations. The mean average absolute difference, removing obvious outliers, between computed and observed energies for the 41 lowest identified levels in Fe XII is only 0.057 \%, implying that the computed energies are accurate enough to aid identification of new emission lines from the sun and other astrophysical sources. The amount of energy and transition data of high accuracy is significantly increased for several P-like ions of astrophysics interest, where experimental data are still very scarce.


title: Phenomenological approach to study the degree of the itinerancy of the $5f$ electrons in actinide ferromagnets with spin fluctuation theory

abstract: Actinide compounds with 5f electrons have been attracting much attention because of their interesting magnetic and electronic properties such as heavy fermion state, unconventional superconductivity, co-existence of the superconductivity and magnetism. Recently, we have reported a phenomenological analysis on 80 actinide ferromagnets with the spin fluctuation theory originally developed to explain the ferromagnetic properties of itinerant ferromagnets in the 3d transition metals and their intermetallics (N. Tateiwa et al., Phys. Rev. B 96, 035125 (2017)). Our study suggests the itinerancy of the $5f$ electrons in most of the actinide ferromagnets and the applicability of the spin fluctuation theory to actinide 5f system. In this paper, we present a new analysis for the spin fluctuation parameter obtained with a different theoretical formula not used in the reference. We also discuss the results of the analysis from different points of views.


title: The maximum negative hypergeometric distribution

abstract: An urn contains a known number of balls of two different colors. We describe the random variable counting the smallest number of draws needed in order to observe at least $\,c\,$ of both colors when sampling without replacement for a pre-specified value of $\,c=1,2,\ldots\,$. This distribution is the finite sample analogy to the maximum negative binomial distribution described by Zhang, Burtness, and Zelterman (2000). We describe the modes, approximating distributions, and estimation of the contents of the urn.


title: Hydrodynamic limits for long-range asymmetric interacting particle systems

abstract: We consider the hydrodynamic scaling behavior of the mass density with respect to a general class of mass conservative interacting particle systems on ${\mathbb Z}^n$, where the jump rates are asymmetric and long-range of order $\|x\|^{-(n+\alpha)}$ for a particle displacement of order $\|x\|$. Two types of evolution equations are identified depending on the strength of the long-range asymmetry. When $0<\alpha<1$, we find a new integro-partial differential hydrodynamic equation, in an anomalous space-time scale. On the other hand, when $\alpha\geq 1$, we derive a Burgers hydrodynamic equation, as in the finite-range setting, in Euler scale.


title: Generic HKT geometries in the harmonic superspace approach

abstract: We explain how a generic HKT geometry can be derived using the language of N = 4 supersymmetric quantum mechanics. To this end, one should consider a Lagrangian involving several (4,4,0) multiplets defined in harmonic superspace and subject to nontrivial harmonic constraints. Conjecturally, this general construction worked out earlier by Delduc and Ivanov gives a complete classification of all HKT geometries. Each such geometry is generated by two different functions (potentials) of a special type that depend on harmonic superfields and on harmonics. Given these two potentials, one can derive the vielbeins, metric, connections and curvatures, but this is not so simple: one should solve rather complicated differential equations. We illustrate the general construction by giving a detailed derivation of the metric for the hyper-Kaehler Taub-NUT manifold. In the generic case, we arrive at an HKT geometry. In this paper, we give a simple proof of this assertion.


title: Variational Integrators for Inertial Magnetohydrodynamics

abstract: Recently, an extended version of magnetohydrodynamics that incorporates electron inertia, dubbed inertial magnetohydrodynamics, has been proposed. This model features a noncanonical Hamiltonian formulation with a number of conserved quantities, including the total energy and modified versions of magnetic and cross helicity. In this work, a variational integrator is presented which preserves these conservation laws to machine accuracy. As long as effects due to finite electron mass are neglected, the scheme preserves the magnetic field line topology so that unphysical reconnection is absent. Only when effects of finite electron mass are added, magnetic reconnection takes place. The excellent conservation properties of the method are illustrated by numerical examples in 2D.