title: Calibration of statistical methods used to constrain pulsar geometries via multiband light curve modelling

abstract: Since its launch in 2008, the Fermi Large Area Telescope (LAT) has detected over 200 -ray pulsars above 100 MeV. This population of pulsars is characterised by a rich diversity of light curve morphologies. Researchers have been using both the radio and -ray light curves to constrain the inclination and observer angles for each of these pulsars. At first, this was done using a by-eye technique and later via statistical approaches. We have also developed two novel statistical approaches that place the radio and -ray data on equal footing despite their disparate relative flux errors. We chose eleven pulsars from the Second Fermi Pulsar Catalog, both old and young, and applied these new techniques as well as the by-eye technique to constrain their geometric parameters using standard pulsar models. We present first results on our comparison of the best-fit parameters yielded by each of the aforementioned techniques. This will assist us in determining the utility of our new statistical approaches, and gauge the overlap of the best-fit parameters (plus errors) from each of the different methods. Such a statistical fitting approach will provide the means for further pulsar magnetospheric model development using light curve data.


title: A closedness theorem over Henselian valued fields with analytic structure

abstract: The main purpose of the paper is to establish a closedness theorem over Henselian valued fields $K$ of equicharacteristic zero (not necessarily algebraically closed) with separated analytic structure. It says that every projection with a projective fiber is a definably closed map. This remains valid also for valued fields with analytic structure induced by a strictly convergent Weierstrass systems, including the classical, complete rank one valued fields with the Tate algebra of strictly convergent power series. As application, we prove two theorems on existence of the limit and on piecewise continuity.


title: Construction and Evaluation of Hierarchical Parcellation of the Brain using fMRI with Prewhitening

abstract: Brain atlases are a ubiquitous tool used for analyzing and interpreting brain imaging datasets. Traditionally, brain atlases divided the brain into regions separated by anatomical landmarks. In the last decade, several attempts have been made to parcellate the brain into regions with distinct functional activity using fMRI. To construct a brain atlas using fMRI, data driven algorithms are used to group voxels with similar functional activity together to form regions. Hierarchical clustering is one parcellation method that has been used for functional parcellation of the brain, resulting in parcellations that align well with cytoarchitectonic divisions of the brain. However, few rigorous data driven evaluations of the method have been performed. Moreover, the effect of removing autocorrelation trends from fMRI time series (prewhitening) on the structure of the resultant atlas has not been previously explored. In this paper, we use hierarchical clustering to produce functional parcellations of the brain using hierarchical clustering. We use both prewhitened and raw fMRI time series to construct the atlas. The resultant atlases were then evaluated for their homogeneity, separation between regions, reproducibility across subjects, and reproducibility across scans.


title: Magnetized Disk-Winds in NGC 3783

abstract: We analyze a 900-ks stacked Chandra/HETG spectrum of NGC 3783 in the context of magnetically-driven accretion-disk wind models in an effort to provide tight constraints on the global conditions of the underlying absorbers. Motivated by the earlier measurements of its absorption measure distribution (AMD) indicating X-ray-absorbing ionic columns that decrease slowly with decreasing ionization parameter, we employ 2D magnetohydrodynamic (MHD) disk-wind models to describe the global outflow. We compute its photoionization structure along with the wind kinematic properties allowing us to further calculate in a self-consistent fashion the shapes of the major X-ray absorption lines. With the wind radial density profile determined by the AMD, the profiles of the ensemble of the observed absorption features are determined by the two global parameters of the MHD wind; i.e. disk inclination \theta_obs and wind density normalization n_o. Considering the most significant absorption features in the (~1.8A-20A) range, we show that the MHD-wind is best described by n(r)~6.9e11(r/ro)^-1.15 [cm^-3] and \theta_obs=44deg. We argue that winds launched by X-ray heating, radiation pressure or even MHD winds but with steeper radial density profiles are strongly disfavored by data. Considering the properties of Fe K band absorption features (i.e. Fe xxv and Fe xxvi), while typically prominent in the AGN X-ray spectra, they appear to be weak in NGC 3783. For the specific parameters of our model obtained by fitting the AMD and the rest of absorption features, these features are found to be weak in agreement with observation.


title: Chromatic splitting for the $K(2)$-local sphere at $p=2$

abstract: We calculate the homotopy type of $L_1L_{K(2)}S^0$ and $L_{K(1)}L_{K(2)}S^0$ at the prime 2, where $L_{K(n)}$ is localization with respect to Morava $K$-theory and $L_1$ localization with respect to $2$-local $K$ theory. In $L_1L_{K(2)}S^0$ we find all the summands predicted by the Chromatic Splitting Conjecture, but we find some extra summands as well. An essential ingredient in our approach is the analysis of the continuous group cohomology $H^\ast_c(\mathbb{G}_2,E_0)$ where $\mathbb{G}_2$ is the Morava stabilizer group and $E_0 = \mathbb{W}[[u_1]]$ is the ring of functions on the height $2$ Lubin-Tate space. We show that the inclusion of the constants $\mathbb{W} \to E_0$ induces an isomorphism on group cohomology, a radical simplification.


title: Possible higher order phase transition in large-$N$ gauge theory at finite temperature

abstract: We analyze the phase structure of $SU(\infty)$ gauge theory at finite temperature using matrix models. Our basic assumption is that the effective potential is dominated by double-trace terms for the Polyakov loops. As a function of the temperature, a background field for the Polyakov loop, and a quartic coupling, it exhibits a universal structure: in the large portion of the parameter space, there is a continuous phase transition analogous to the third-order phase transition of Gross, Witten and Wadia, but the order of phase transition can be higher than third. We show that different confining potentials give rise to drastically different behavior of the eigenvalue density and the free energy. Therefore lattice simulations at large $N$ could probe the order of phase transition and test our results.


title: Brownian motion on Perelman's almost Ricci-flat manifold

abstract: We study Brownian motion and stochastic parallel transport on Perelman's almost Ricci flat manifold $\mathscr M=M\times \mathbb S^N\times I$, whose dimension depends on a parameter $N$ unbounded from above. We construct sequences of projected Brownian motions and stochastic parallel transports which for $N \to \infty$ converge to the corresponding objects for the Ricci flow. In order to make precise this process of passing to the limit, we study the martingale problems for the Laplace operator on $\mathscr M$ and for the horizontal Laplacian on the orthonormal frame bundle $\mathscr{OM}$ . As an application, we see how the characterizations of two-sided bounds on the Ricci curvature established by A. Naber applied to Perelman's manifold lead to the inequalities that characterize solutions of the Ricci flow discovered by Naber and the second author.


title: Universality at large transverse spin in defect CFT

abstract: We study the spectrum of local operators living on a defect in a generic conformal field theory, and their coupling to the local bulk operators. We establish the existence of universal accumulation points in the spectrum at large $s$, $s$ being the charge of the operators under rotations in the space transverse to the defect. Our tools include a formula that inverts the bulk to defect OPE, analogous to the Caron-Huot formula for the four-point function. Analyticity of the formula in $s$ implies that the scaling dimensions of the defect operators are aligned in Regge trajectories $\widehat{\Delta}(s)$. These results require the correlator of two local operators and the defect to be bounded in a certain region, a condition that we do not prove in general. We check our conclusions against examples in perturbation theory and holography, and we make specific predictions concerning the spectrum of defect operators on Wilson lines. We also give an interpretation of the large $s$ spectrum in the spirit of the work of Alday and Maldacena.


title: Mixed weak estimates of Sawyer type for fractional integrals and some related operators

abstract: We prove mixed weak estimates of Sawyer type for fractional operators. More precisely, let $\mathcal{T}$ be either the maximal fractional function $M_\gamma$ or the fractional integral operator $I_\gamma$, $0<\gamma<n$, $1\leq p<n/\gamma$ and $1/q=1/p-\gamma/n$. If $u,v^{q/p}\in A_1$ or if $uv^{-q/{p'}}\in A_1$ and $v^q\in A_\infty(uv^{-q/{p'}})$ then we obtain that the estimate \begin{equation*} uv^{q/p}\left(\left\{x\in \R^n: \frac{|\mathcal{T}(fv)(x)|}{v(x)}>t\right\}\right)^{1/q}\leq \frac{C}{t}\left(\int_{\R^n}|f(x)|^pu(x)^{p/q}v(x)\,dx\right)^{1/p}, \end{equation*} holds for every positive $t$ and every bounded function with compact support. As an important application of the results above we further more exhibe mixed weak estimates for commutators of Calder\'on-Zygmund singular integral and fractional integral operators when the symbol $b$ is in the class Lipschitz-$\delta$, $0<\delta\leq 1$.


title: The Role of Molecular Properties on the Dimerization of Aromatic Compounds

abstract: Recent results have shown the presence and importance of oxygen chemistry during the growth of aromatic compounds, leading to the formation of oxygenated structures that have been identified in various environments. Since the formation of polycyclic aromatic compounds (PAC) bridge the formation of gas-phase species with particle inception, in this work we report a detailed analysis of the effects of molecular characteristics on physical growth of PAC via dimerization. We have included oxygen content, mass, type of bonds (rigid versus rotatable), and shape as main properties of the molecules and studied their effect on the propensity of these structures to form homo-molecular and hetero-molecular dimers. Using enhanced sampling molecular dynamics techniques, we have computed the free energy of dimerization in the temperature range $500-1680$~K. Initial structures used in this study were obtained from experimental data. The results show that although the effects of shape, presence of oxygen, mass, and internal bonds are tightly intertwined, and their relative importance changes with temperature. In general, mass and the presence of rotatable bonds are the best indicators. The results provide knowledge on the inception step and the role that particle characteristics play during inception. In addition, our study highlights the fact that current models that use stabilomers as monomers for physical aggregation are overestimating the importance of this process during particle nucleation.


title: A Holistic Approach to Evaluate EMI Shielding Characteristics of Carbon Nanotube-based Polymer Composites

abstract: The work presents a comprehensive methodology to determine the Shielding Effectiveness (SE) of single-walled carbon nanotubes (SWCNTs)/polymer nanocomposites. Here, an algorithm based on Ant Colony Optimization (ACO) was employed to determine the electrical conductivity ({\sigma}) of these nanocomposites as a function of SWCNT concentration. Then, these {\sigma} values were used to compute the SE values as a function of frequency and concentration. Specifically, a pseudo three-dimensional (3D) percolation model was developed to study the effects of random connectivity of SWCNTs to one another on the {\sigma} values of the nanocomposites. Both intrinsic and tunneling resistances were taken into account. The consequences of the presence of both well-exfoliated and aggregated SWCNTs with varying lengths distributed inhomogeneously on {\sigma} and SE values were investigated.


title: Adjacency and Tensor Representation in General Hypergraphs Part 1: e-adjacency Tensor Uniformisation Using Homogeneous Polynomials

abstract: Adjacency between two vertices in graphs or hypergraphs is a pairwise relationship. It is redefined in this article as 2-adjacency. In general hypergraphs, hyperedges hold for $n$-adic relationship. To keep the $n$-adic relationship the concepts of $k$-adjacency and e-adjacency are defined. In graphs 2-adjacency and e-adjacency concepts match, just as $k$-adjacency and e-adjacency do for $k$-uniform hypergraphs. For general hypergraphs these concepts are different. This paper also contributes in a uniformization process of a general hypergraph to allow the definition of an e-adjacency tensor, viewed as a hypermatrix, reflecting the general hypergraph structure. This symmetric e-adjacency hypermatrix allows to capture not only the degree of the vertices and the cardinality of the hyperedges but also makes a full separation of the different layers of a hypergraph.


title: QCD phase diagram for nonzero isospin-asymmetry

abstract: The QCD phase diagram is studied in the presence of an isospin asymmetry using continuum extrapolated staggered quarks with physical masses. In particular, we investigate the phase boundary between the normal and the pion condensation phases and the chiral/deconfinement transition. The simulations are performed with a small explicit breaking parameter in order to avoid the accumulation of zero modes and thereby stabilize the algorithm. The limit of vanishing explicit breaking is obtained by means of an extrapolation, which is facilitated by a novel improvement program employing the singular value representation of the Dirac operator. Our findings indicate that no pion condensation takes place above $T\approx 160$ MeV and also suggest that the deconfinement crossover continuously connects to the BEC-BCS crossover at high isospin asymmetries. The results may be directly compared to effective theories and model approaches to QCD.


title: Prominence formation and ejection in cool stars

abstract: The observational signatures of prominences have been detected in single and binary G and K type stars for many years now, but recently this has been extended to the M dwarf regime. Prominences carry away both mass and angular momentum when they are ejected and the impact of this mass on any orbiting planets may be important for the evolution of exoplanetary atmospheres. By means of the classification used in the massive star community, that involves knowledge of two parameters (the co-rotation and Alfv\'en radii, $r_K$ and $r_A$), we have determined which cool stars could support prominences. From a model of mechanical support, we have determined that the prominence mass $m_p/M_\star=(E_M/E_G)(r_\star/r_K)^2 F$ where $E_MB_\star^2r_\star^3$ and $E_G = GM_\star^2/r_\star$ are magnetic and gravitational energies and $F$ is a geometric factor. Our calculated masses and ejection frequencies (typically $10^{16}-10^{17}$g and 0.4 d, respectively) are consistent with observations and are sufficient to ensure that an exoplanet orbiting in the habitable zone of an M dwarf could suffer frequent impacts.


title: Structured eigenvalue/eigenvector backward errors of matrix pencils arising in optimal control

abstract: Eigenvalue and eigenpair backward errors are computed for matrix pencils arising in optimal control. In particular, formulas for backward errors are developed that are obtained under block-structure-preserving and symmetry-structure-preserving perturbations. It is shown that these eigenvalue and eigenpair backward errors are sometimes significantly larger than the corresponding backward errors that are obtained under perturbations that ignore the special structure of the pencil.


title: A presentation for the Cox ring of $\overline{M}_{0,6}$

abstract: We compute the relations in the Cox ring of the moduli space $\overline{M}_{0,6}$. This gives a presentation of the Cox ring as a quotient of a polynomial ring with 40 generators by an ideal with 225 generators that come in 5 symmetry classes.


title: Self-Similar $k$-Graph C*-Algebras

abstract: In this paper, we introduce a notion of a self-similar action of a group $G$ on a $k$-graph $\Lambda$, and associate it a universal C*-algebra $\O_{G,\Lambda}$. We prove that $\O_{G,\Lambda}$ can be realized as the Cuntz-Pimsner algebra of a product system. If $G$ is amenable and the action is pseudo free, then $\O_{G,\Lambda}$ is shown to be isomorphic to a "path-like" groupoid C*-algebra. This facilitates studying the properties of $\O_{G,\Lambda}$. We show that $\O_{G,\Lambda}$ is always nuclear and satisfies the Universal Coefficient Theorem; we characterize the simplicity of $\O_{G,\Lambda}$ in terms of the underlying action; and we prove that, whenever $\O_{G,\Lambda}$ is simple, there is a dichotomy: it is either stably finite or purely infinite, depending on whether $\Lambda$ has nonzero graph traces or not. Our main results generalize the recent work of Exel and Pardo on self-similar graphs.


title: Choreographic and Somatic Approaches for the Development of Expressive Robotic Systems

abstract: As robotic systems are moved out of factory work cells into human-facing environments questions of choreography become central to their design, placement, and application. With a human viewer or counterpart present, a system will automatically be interpreted within context, style of movement, and form factor by human beings as animate elements of their environment. The interpretation by this human counterpart is critical to the success of the system's integration: knobs on the system need to make sense to a human counterpart; an artificial agent should have a way of notifying a human counterpart of a change in system state, possibly through motion profiles; and the motion of a human counterpart may have important contextual clues for task completion. Thus, professional choreographers, dance practitioners, and movement analysts are critical to research in robotics. They have design methods for movement that align with human audience perception, can identify simplified features of movement for human-robot interaction goals, and have detailed knowledge of the capacity of human movement. This article provides approaches employed by one research lab, specific impacts on technical and artistic projects within, and principles that may guide future such work. The background section reports on choreography, somatic perspectives, improvisation, the Laban/Bartenieff Movement System, and robotics. From this context methods including embodied exercises, writing prompts, and community building activities have been developed to facilitate interdisciplinary research. The results of this work is presented as an overview of a smattering of projects in areas like high-level motion planning, software development for rapid prototyping of movement, artistic output, and user studies that help understand how people interpret movement. Finally, guiding principles for other groups to adopt are posited.


title: Heat Equation With a Geometric Rough Path Potential in One Space Dimension: Existence and Regularity of Solution

abstract: A solution of the heat equation with a distribution-valued potential is constructed by regularization. When the potential is the generalized derivative of a H\"{o}lder continuous function, regularity of the resulting solution is in line with the standard parabolic theory.


title: Fair Forests: Regularized Tree Induction to Minimize Model Bias

abstract: The potential lack of fairness in the outputs of machine learning algorithms has recently gained attention both within the research community as well as in society more broadly. Surprisingly, there is no prior work developing tree-induction algorithms for building fair decision trees or fair random forests. These methods have widespread popularity as they are one of the few to be simultaneously interpretable, non-linear, and easy-to-use. In this paper we develop, to our knowledge, the first technique for the induction of fair decision trees. We show that our "Fair Forest" retains the benefits of the tree-based approach, while providing both greater accuracy and fairness than other alternatives, for both "group fairness" and "individual fairness.'" We also introduce new measures for fairness which are able to handle multinomial and continues attributes as well as regression problems, as opposed to binary attributes and labels only. Finally, we demonstrate a new, more robust evaluation procedure for algorithms that considers the dataset in its entirety rather than only a specific protected attribute.


title: Recursive Programs for Document Spanners

abstract: A document spanner models a program for Information Extraction (IE) as a function that takes as input a text document (string over a finite alphabet) and produces a relation of spans (intervals in the document) over a predefined schema. A well studied language for expressing spanners is that of the regular spanners: relational algebra over regex formulas, which are obtained by adding capture variables to regular expressions. Equivalently, the regular spanners are the ones expressible in non-recursive Datalog over regex formulas (extracting relations that play the role of EDBs from the input document). In this paper, we investigate the expressive power of recursive Datalog over regex formulas. Our main result is that such programs capture precisely the document spanners computable in polynomial time. Additional results compare recursive programs to known formalisms such as the language of core spanners (that extends regular spanners by allowing to test for string equality) and its closure under difference. Finally, we extend our main result to a recently proposed framework that generalizes both the relational model and document spanners.


title: LHC data and its impact on nCTEQ15 PDFs

abstract: The LHC heavy ion data for W/Z production can provide new incisive information on the PDFs. This data is sensitive to the heavier quark flavors (strange and charm) in a high energy kinematic region; this can facilitate the determination of PDFs in the small x region where previous data was limited. At present, the flavor separation of the proton PDFs is dependent on DIS data from nuclear targets. Therefore, improved nuclear corrections can also yield enhanced flavor determination of both the proton and nuclear PDFs.


title: Repulsive Casimir force in Bose-Einstein Condensate

abstract: We study the Casimir effect for a three dimensional system of ideal free massive Bose gas in a slab geometry with Zaremba and anti-periodic boundary conditions. It is found that for these type of boundary conditions the resulting Casimir force is repulsive in nature, in contrast with usual periodic, Dirichlet or Neumann boundary condition where the Casimir force is attractive (Martin P. A. and Zagrebnov V. A., Europhys. Lett., 73 (2006) 15.). Casimir forces in these boundary conditions also maintain a power law decay function below condensation temperature and exponential decay function above the condensation temperature albeit with a positive sign, identifying the repulsive nature of the force.


title: Multilevel LDPC Lattices with Efficient Encoding and Decoding and a Generalization of Construction D'

abstract: Lattice codes are elegant and powerful structures that not only can achieve the capacity of the AWGN channel but are also a key ingredient to many multiterminal schemes that exploit linearity properties. However, constructing lattice codes that can realize these benefits with low complexity is still a challenging problem. In this paper, efficient encoding and decoding algorithms are proposed for multilevel binary LDPC lattices constructed via Construction D' whose complexity is linear in the total number of coded bits. Moreover, a generalization of Construction D' is proposed that relaxes some of the nesting constraints on the parity-check matrices of the component codes, leading to a simpler and improved design. Based on this construction, low-complexity multilevel LDPC lattices are designed whose performance under multistage decoding is comparable to that of polar lattices and close to that of low-density lattice codes (LDLC) on the power-unconstrained AWGN channel.


title: Excitation of epsilon-near-zero resonance in ultra-thin indium tin oxide shell embedded nanostructured optical fiber

abstract: We report a novel optical waveguide design of a hollow step index fiber modified with a thin layer of indium tin oxide (ITO). We show an excitation of highly confined waveguide mode in the proposed fiber near the wavelength where permittivity of ITO approaches zero. Due to the high field confinement within thin ITO shell inside the fiber, the epsilon-near-zero (ENZ) mode can be characterized by a peak in modal loss of the hybrid waveguide. Our results show that such in-fiber excitation of ENZ mode is due to the coupling of the guided core mode to the thin-film ENZ mode. We also show that the phase matching wavelength, where the coupling takes place, varies depending on the refractive index of the constituents inside the central bore of the fiber. These ENZ nanostructured optical fibers have many potential applications, for example, in ENZ nonlinear and magneto-optics, as in-fiber wavelength-dependent filters, and as subwavelength fluid channel for optical and bio-photonic sensing.


title: Field Theory Done Right

abstract: An effective formalism for white noise analysis, conceptually equivalent to Wilsonian renormalization theory, is introduced. Space-time gets represented by a boolean lattice of coarse regions, energy scales become space-time partitions by lattice regions, and observables are elements of a projective limit with connecting maps given by partial integration of high-energy degrees of freedom. The framework allows for a seamless generalization of the Wick product and the $\mathcal S$-transform to essentially arbitrary L\'evy noises, and we provide a tool to make explicit calculations in several cases of interest, including Gauss, Poisson and Gamma noises (we shall thereby encounter pretty familiar polynomials, like falling factorials and Hermite polynomials). Armed with this, we turn to constructive quantum field theory. We adopt an Euclidean approach and introduce a sufficient condition for reflection positivity, based on our $\mathcal S$-transform, enabling us to construct non-trivial quantum fields by simply specifying compatible families of effective connected $n$-point functions. We exemplify this by producing a field with quartic interaction in dimension $d\leq 8$. Its connected $n$-point functions vanish except for the propagator and the connected $4$-point function, which is that of the $\phi^4$ field up to order $\hbar$. This model satisfies all the physical requirements of a non-trivial quantum field theory.

categories: math-ph cond-mat.stat-mech hep-th math.MP


title: Exchange-torque-induced excitation of perpendicular standing spin waves in nanometer-thick YIG films

abstract: Spin waves in ferrimagnetic yttrium iron garnet (YIG) films with ultralow magnetic damping are relevant for magnon-based spintronics and low-power wave-like computing. The excitation frequency of spin waves in YIG is rather low in weak external magnetic fields because of its small saturation magnetization, which limits the potential of YIG films for high-frequency applications. Here, we demonstrate how exchange-coupling to a CoFeB film enables efficient excitation of high-frequency perpendicular standing spin waves (PSSWs) in nanometer-thick (80 nm and 295 nm) YIG films using uniform microwave magnetic fields. In the 295-nm-thick YIG film, we measure intense PSSW modes up to 10th order. Strong hybridization between the PSSW modes and the ferromagnetic resonance mode of CoFeB leads to characteristic anti-crossing behavior in broadband spin-wave spectra. A dynamic exchange torque at the YIG/CoFeB interface explains the excitation of PSSWs. The localized torque originates from exchange coupling between two dissimilar magnetization precessions in the YIG and CoFeB layers. As a consequence, spin waves are emitted from the YIG/CoFeB interface and PSSWs form when their wave vector matches the perpendicular confinement condition. PSSWs are not excited when the exchange coupling between YIG and CoFeB is suppressed by a Ta spacer layer. Micromagnetic simulations confirm the exchange-torque mechanism.


title: Practically-Self-Stabilizing Vector Clocks in the Absence of Execution Fairness

abstract: Vector clock algorithms are basic wait-free building blocks that facilitate causal ordering of events. As wait-free algorithms, they are guaranteed to complete their operations within a finite number of steps. Stabilizing algorithms allow the system to recover after the occurrence of transient faults, such as soft errors and arbitrary violations of the assumptions according to which the system was designed to behave. We present the first, to the best of our knowledge, stabilizing vector clock algorithm for asynchronous crash-prone message-passing systems that can recover in a wait-free manner after the occurrence of transient faults. In these settings, it is challenging to demonstrate a finite and wait-free recovery from (communication and crash failures as well as) transient faults, bound the message and storage sizes, deal with the removal of all stale information without blocking, and deal with counter overflow events (which occur at different network nodes concurrently). We present an algorithm that never violates safety in the absence of transient faults and provides bounded time recovery during fair executions that follow the last transient fault. The novelty is that in the absence of execution fairness, the algorithm guarantees a bound on the number of times in which the system might violate safety (while existing algorithms might block forever due to the presence of both transient faults and crash failures). Since vector clocks facilitate a number of elementary synchronization building blocks (without requiring remote replica synchronization) in asynchronous systems, we believe that our analytical insights are useful for the design of other systems that cannot guarantee execution fairness.


title: Inelastic Surface Growth

abstract: Inelastic surface growth associated with continuous creation of incompatibility on the boundary of an evolving body is behind a variety of natural and technological processes, including embryonic development and 3D printing. In this paper we extend the recently proposed stress-space-centered theory of surface growth (PRL 119, 048001, 2017) by shifting the focus towards growth induced strains. To illustrate the new development we present several analytically tractable examples.


title: Variational Attention for Sequence-to-Sequence Models

abstract: The variational encoder-decoder (VED) encodes source information as a set of random variables using a neural network, which in turn is decoded into target data using another neural network. In natural language processing, sequence-to-sequence (Seq2Seq) models typically serve as encoder-decoder networks. When combined with a traditional (deterministic) attention mechanism, the variational latent space may be bypassed by the attention model, and thus becomes ineffective. In this paper, we propose a variational attention mechanism for VED, where the attention vector is also modeled as Gaussian distributed random variables. Results on two experiments show that, without loss of quality, our proposed method alleviates the bypassing phenomenon as it increases the diversity of generated sentences.


title: Non-equilibrium effects on thermal ignition using hard sphere molecular dynamics simulations

abstract: The present study addresses the role of molecular non-equilibrium effects in thermal ignition problems. We consider a single binary reaction of the form A+B -> C+C. Molecular dynamics calculations were performed for activation energies ranging between RT and 7.5RT and heat release of 2.5RT and 10RT. The evolution of up to 10,000 particles was calculated as the system undergoes a thermal ignition at constant volume. Ensemble averages of 100 calculations for each parameter set permitted to determine the ignition delay, along with a measure of the stochasticity of the process. A well behaved convergence to large system sizes is also demonstrated. The ignition delay calculations were compared with those obtained at the continuum level using rates derived from kinetic theory: the standard rate assuming that the distribution of the speed of the particles is the Maxwell-Boltzmann distribution, and the perturbed rates by Prigogine and Xhrouet [1] for an isothermal system, and Prigogine and Mahieu [2] for an energy releasing reaction, obtained by the Chapman-Enskog perturbation procedure. The molecular results were found in very good agreement with the latter at low temperatures, confirming that non-equilibrium effects promote the formation of energetic particles, that serve as seeds for subsequent reaction events: i.e., hot spots. This effect was found to lower the ignition delay by up to 30%. At high temperatures, the ignition delay obtained from the standard equilibrium rate was found to be up to 60% longer than the molecular calculations. This effect is due to the rapidity of the reactive collisions that do not allow the system to equilibrate. For this regime, none of the perturbation solutions obtained by the Chapman-Enskog procedure were valid. This study thus shows the importance of non-equilibrium effects in thermal ignition problems, for most temperatures of practical interest.


title: On State Observers for Nonlinear Systems: A New Design and a Unifying Framework

abstract: In this paper we propose a new observer design technique for nonlinear systems. It combines the well-known Kazantzis-Kravaris-Luenberger observer and the recently introduced parameter estimation-based observer, which become special cases of it---extending the realm of applicability of both methods. A second contribution of the paper is the proof that these designs can be recast as particular cases of immersion and invariance observers---providing in this way a unified framework for their analysis and design. Simulation results of a physical system that illustrates the superior performance of the proposed observer compared to other existing observers are presented.


title: A nonamenable "factor" of a Euclidean space

abstract: Answering a question of Benjamini, we present an isometry-invariant random partition of the Euclidean space $\mathbb{R}^d$, $d\geq 3$, into infinite connected indistinguishable pieces, such that the adjacency graph defined on the pieces is the 3-regular infinite tree. Along the way, it is proved that any finitely generated one-ended amenable Cayley graph can be represented in $\mathbb{R}^d$ as an isometry-invariant random partition of $\mathbb{R}^d$ to bounded polyhedra, and also as an isometry-invariant random partition of $\mathbb{R}^d$ to indistinguishable pieces. A new technique is developed to prove indistinguishability for certain constructions, connecting this notion to factor of iid's.


title: Robust Detection of Covariate-Treatment Interactions in Clinical Trials

abstract: Detection of interactions between treatment effects and patient descriptors in clinical trials is critical for optimizing the drug development process. The increasing volume of data accumulated in clinical trials provides a unique opportunity to discover new biomarkers and further the goal of personalized medicine, but it also requires innovative robust biomarker detection methods capable of detecting non-linear, and sometimes weak, signals. We propose a set of novel univariate statistical tests, based on the theory of random walks, which are able to capture non-linear and non-monotonic covariate-treatment interactions. We also propose a novel combined test, which leverages the power of all of our proposed univariate tests into a single general-case tool. We present results for both synthetic trials as well as real-world clinical trials, where we compare our method with state-of-the-art techniques and demonstrate the utility and robustness of our approach.


title: Submodular Optimization for Consensus Networks with Noise-Corrupted Leaders

abstract: We consider the leader selection problem in a network with consensus dynamics where both leader and follower agents are subject to stochastic external disturbances. The performance of the system is quantified by the total steady-state variance of the node states, and the goal is to identify the set of leaders that minimizes this variance. We first show that this performance measure can be expressed as a submodular set function over the nodes in the network. We then use this result to analyze the performance of two greedy, polynomial-time algorithms for leader selection, showing that the leader sets produced by the greedy algorithms are within provable bounds of optimal.


title: The nonlinear heat equation involving highly singular initial values and new blowup and life span results

abstract: In this paper we prove local existence of solutions to the nonlinear heat equation $u_t = \Delta u +a |u|^\alpha u, \; t\in(0,T),\; x=(x_1,\,\cdots,\, x_N)\in {\mathbb R}^N,\; a = \pm 1,\; \alpha>0;$ with initial value $u(0)\in L^1_{\rm{loc}}\left({\mathbb R}^N\setminus\{0\}\right)$, anti-symmetric with respect to $x_1,\; x_2,\; \cdots,\; x_m$ and $|u(0)|\leq C(-1)^m\partial_{1}\partial_{2}\cdot \cdot \cdot \partial_{m}(|x|^{-\gamma})$ for $x_1>0,\; \cdots,\; x_m>0,$ where $C>0$ is a constant, $m\in \{1,\; 2,\; \cdots,\; N\},$ $0<\gamma<N$ and $0<\alpha<2/(\gamma+m).$ This gives a local existence result with highly singular initial values. As an application, for $a=1,$ we establish new blowup criteria for $0<\alpha\leq 2/(\gamma+m)$, including the case $m=0.$ Moreover, if $(N-4)\alpha<2,$ we prove the existence of initial values $u_0 = \lambda f,$ for which the resulting solution blows up in finite time $T_{\max}(\lambda f),$ if $\lambda>0$ is sufficiently small. We also construct blowing up solutions with initial data $\lambda_n f$ such that $\lambda_n^{[({1\over \alpha}-{\gamma+m\over 2})^{-1}]}T_{\max}(\lambda_n f)$ has different finite limits along different sequences $\lambda_n\to 0$. Our result extends the known "small lambda" blow up results for new values of $\alpha$ and a new class of initial data.


title: The length and depth of algebraic groups

abstract: Let $G$ be a connected algebraic group. An unrefinable chain of $G$ is a chain of subgroups $G = G_0 > G_1 > \cdots > G_t = 1$, where each $G_i$ is a maximal connected subgroup of $G_{i-1}$. We introduce the notion of the length (respectively, depth) of $G$, defined as the maximal (respectively, minimal) length of such a chain. Working over an algebraically closed field, we calculate the length of a connected group $G$ in terms of the dimension of its unipotent radical $R_u(G)$ and the dimension of a Borel subgroup $B$ of the reductive quotient $G/R_u(G)$. In particular, a simple algebraic group of rank $r$ has length $\dim B + r$, which gives a natural extension of a theorem of Solomon and Turull on finite quasisimple groups of Lie type. We then deduce that the length of any connected algebraic group $G$ exceeds $\frac{1}{2} \dim G$. We also study the depth of simple algebraic groups. In characteristic zero, we show that the depth of such a group is at most $6$ (this bound is sharp). In the positive characteristic setting, we calculate the exact depth of each exceptional algebraic group and we prove that the depth of a classical group (over a fixed algebraically closed field of positive characteristic) tends to infinity with the rank of the group. Finally we study the chain difference of an algebraic group, which is the difference between its length and its depth. In particular we prove that, for any connected algebraic group $G$, the dimension of $G/R(G)$ is bounded above in terms of the chain difference of $G$.


title: Diverse spatial expression patterns emerge from common transcription bursting kinetics

abstract: In early development, regulation of transcription results in precisely positioned and highly reproducible expression patterns that specify cellular identities. How transcription, a fundamentally noisy molecular process, is regulated to achieve reliable embryonic patterning remains unclear. In particular, it is unknown how gene-specific regulation mechanisms affect kinetic rates of transcription, and whether there are common, global features that govern these rates across a genetic network. Here, we measure nascent transcriptional activity in the gap gene network of early Drosophila embryos and characterize the variability in absolute activity levels across expression boundaries. We demonstrate that boundary formation follows a common transcriptional principle: a single control parameter determines the distribution of transcriptional activity, regardless of gene identity, boundary position, or enhancer-promoter architecture. By employing a minimalist model of transcription, we infer kinetic rates of transcriptional bursting for these patterning genes; we find that the key regulatory parameter is the fraction of time a gene is in an actively transcribing state, while the rate of Pol II loading appears globally conserved. These results point to a universal simplicity underlying the apparently complex transcriptional processes responsible for early embryonic patterning and indicate a path to general rules in transcriptional regulation.


title: Nematic quantum phases in the bilayer honeycomb antiferromagnet

abstract: The spin$-1/2$ Heisenberg antiferromagnet on the honeycomb bilayer lattice is shown to display a rich variety of semiclassical and genuinely quantum phases, controlled by the interplay between intralayer frustration and interlayer exchange. Employing a complementary set of techniques, comprising spin rotationally invariant Schwinger boson mean field theory, bond operators, and series expansions we unveil the quantum phase diagram, analyzing low-energy excitations and order parameters. By virtue of Schwinger bosons we scan the complete range of exchange parameters, covering both long range ordered as well as quantum disordered ground states and reveal the existence of an extended, frustration induced lattice nematic phase in a range of intermediate exchange unexplored so far.


title: On fundamental diffraction limitation of finesse of a Fabry-Perot cavity

abstract: We perform a theoretical study of finesse limitations of a Fabry-Perot (FP) cavity occurring due to finite size, asymmetry, as well as imperfections of the cavity mirrors. A method of numerical simulations of the eigenvalue problem applicable for both the fundamental and high order cavity modes is suggested. Using this technique we find spatial profile of the modes and their round-trip diffraction loss. The results of the numerical simulations and analytical calculations are nearly identical when we consider a conventional FP cavity. The proposed numerical technique has much broader applicability range and is valid for any FP cavity with arbitrary non-spherical mirrors which have cylindrical symmetry but disturbed in an asymmetric way, for example, by tilt or roughness of their mirrors.


title: Well-Balanced Schemes for the Euler Equations with Gravitation: Conservative Formulation Using Global Fluxes

abstract: We develop a second-order well-balanced central-upwind scheme for the compressible Euler equations with gravitational source term. Here, we advocate a new paradigm based on a purely conservative reformulation of the equations using global fluxes. The proposed scheme is capable of exactly preserving steady-state solutions expressed in terms of a nonlocal equilibrium variable. A crucial step in the construction of the second-order scheme is a well-balanced piecewise linear reconstruction of equilibrium variables combined with a well-balanced central-upwind evolution in time, which is adapted to reduce the amount of numerical viscosity when the flow is at (near) steady-state regime. We show the performance of our newly developed central-upwind scheme and demonstrate importance of perfect balance between the fluxes and gravitational forces in a series of one- and two-dimensional examples.


title: Duration of Classicality of Homogeneous Condensates with Attractive Interactions

abstract: Dark matter axions and other highly degenerate bosonic fluids are commonly described by classical field equations. In a recent paper \cite{BECprop} we calculated the duration of classicality of homogeneous condensates with attractive contact interactions and of self-gravitating homogeneous condensates in critical expansion. According to their classical equations of motion, such condensates persist forever. In their quantum evolution parametric resonance causes quanta to jump in pairs out of the condensate into all modes with wavevector less than some critical value. We estimated in each case the time scale over which the condensate is depleted and after which a classical description is invalid.


title: Probability current in zero-spin relativistic quantum mechanics

abstract: We show that the antisymmetric spinor tensor representation of spin-0 relativistic quantum mechanics provides a conserved current with positive definite timelike component, interpretable as probability density. The construction runs in complete analogy to the spin-1/2 case, and provides an analogously natural one-particle Hilbert space description for spin 0. Except for the free particle, the obtained formulation proves to be inequivalent to the one based on the Klein--Gordon equation. The second quantized version may lead to new field theoretical interaction terms for zero-spin particles.


title: Tackling Contention Through Cooperation: A Distributed Federation in LoRaWAN Space

abstract: Low-Power Wide Area Networks (LPWAN) play a key role in the IoT marketplace wherein LoRaWAN is considered a leading solution. Despite the traction of LoRaWAN, research shows that the current contention management mechanisms of LoRaWAN do not scale. This paper tackles contention on LoRaWAN by introducing FLIP, a fully distributed and open architecture for LoRaWAN that fundamentally rethinks how LoRa gateways should be managed and coordinated. FLIP transforms LoRa gateways into a federated network that provides inherent support for roaming while tackling contention using consensus-driven load balancing. FLIP offers identical security guarantees to LoRaWAN, is compatible with existing gateway hardware and requires no updates to end-device hardware or firmware. These features ensure the practicality of FLIP and provide a path to its adoption. We evaluate the performance of FLIP in a large-scale real-world deployment and demonstrate that FLIP delivers scalable roaming and improved contention management in comparison to LoRaWAN. FLIP achieves these benefits within the resource constraints of conventional LoRa gateways and requires no server hardware.


title: An Economic Study of the Effect of Android Platform Fragmentation on Security Updates

abstract: Vendors in the Android ecosystem typically customize their devices by modifying Android Open Source Project (AOSP) code, adding in-house developed proprietary software, and pre-installing third-party applications. However, research has documented how various security problems are associated with this customization process. We develop a model of the Android ecosystem utilizing the concepts of game theory and product differentiation to capture the competition involving two vendors customizing the AOSP platform. We show how the vendors are incentivized to differentiate their products from AOSP and from each other, and how prices are shaped through this differentiation process. We also consider two types of consumers: security-conscious consumers who understand and care about security, and na\"ive consumers who lack the ability to correctly evaluate security properties of vendor-supplied Android products or simply ignore security. It is evident that vendors shirk on security investments in the latter case. Regulators such as the U.S. Federal Trade Commission have sanctioned Android vendors for underinvestment in security, but the exact effects of these sanctions are difficult to disentangle with empirical data. Here, we model the impact of a regulator-imposed fine that incentivizes vendors to match a minimum security standard. Interestingly, we show how product prices will decrease for the same cost of customization in the presence of a fine, or a higher level of regulator-imposed minimum security.


title: The Dirichlet-to-Neumann operator for quantum graphs

abstract: For a compact, connected metric graphs with a boundary that consists of $k$ vertices, we prove that an arbitrary symmetric $k\times k$ matrix with real entries can be realized as the Dirichlet-to-Neumann operator for the Laplacian plus a constant.


title: Non-thermal states in models of filaments: a dynamical study

abstract: We study the origin of the non-thermal profiles observed in filamentary structures in galactic molecular clouds by means of numerical dynamical simulations. We find that such profiles are intrinsic features of the end products of dissipationless collapse in cylindrical symmetry. Moreover, for sufficiently cold initial conditions, we obtain end states characterized by markedly anticorrelated radial density and temperature profiles. Gravitational, dissipationless dynamics alone is thus sufficient to reproduce, at least qualitatively, many of the properties of the observed non-thermal structures.


title: A fast and memory-efficient spectral Galerkin scheme for distributed elliptic optimal control problems

abstract: Many scientific and engineering challenges can be formulated as optimization problems which are constrained by partial differential equations (PDEs). These include inverse problems, control problems, and design problems. As a major challenge, the associated optimization procedures are inherently large-scale. To ensure computational tractability, the design of efficient and robust iterative methods becomes imperative. To meet this challenge, this paper introduces a fast and memory-efficient preconditioned iterative scheme for a class of distributed optimal control problems governed by convection-diffusion-reaction (CDR) equations. As an alternative to low-order discretizations and Schur-complement block preconditioners, the scheme combines a high-order spectral Galerkin method with an efficient preconditioner tailored specifically for the CDR application. The preconditioner is matrix-free and can be applied within linear complexity where the proportionality constant is small. Numerical results demonstrate that the preconditioner is ideal in the sense that appropriate Krylov subspace methods converge within a low number of iterations, independently of the problem size and the Tikhonov regularization parameter.


title: A Primal-Dual Method for Optimal Control and Trajectory Generation in High-Dimensional Systems

abstract: Presented is a method for efficient computation of the Hamilton-Jacobi (HJ) equation for time-optimal control problems using the generalized Hopf formula. Typically, numerical methods to solve the HJ equation rely on a discrete grid of the solution space and exhibit exponential scaling with dimension. The generalized Hopf formula avoids the use of grids and numerical gradients by formulating an unconstrained convex optimization problem. The solution at each point is completely independent, and allows a massively parallel implementation if solutions at multiple points are desired. This work presents a primal-dual method for efficient numeric solution and presents how the resulting optimal trajectory can be generated directly from the solution of the Hopf formula, without further optimization. Examples presented have execution times on the order of milliseconds and experiments show computation scales approximately polynomial in dimension with very small high-order coefficients.


title: Analysis-synthesis model learning with shared features: a new framework for histopathological image classification

abstract: Automated histopathological image analysis offers exciting opportunities for the early diagnosis of several medical conditions including cancer. There are however stiff practical challenges: 1.) discriminative features from such images for separating diseased vs. healthy classes are not readily apparent, and 2.) distinct classes, e.g. healthy vs. stages of disease continue to share several geometric features. We propose a novel Analysis-synthesis model Learning with Shared Features algorithm (ALSF) for classifying such images more effectively. In ALSF, a joint analysis and synthesis learning model is introduced to learn the classifier and the feature extractor at the same time. In this way, the computation load in patch-level based image classification can be much reduced. Crucially, we integrate into this framework the learning of a low rank shared dictionary and a shared analysis operator, which more accurately represents both similarities and differences in histopathological images from distinct classes. ALSF is evaluated on two challenging databases: (1) kidney tissue images provided by the Animal Diagnosis Lab (ADL) at the Pennsylvania State University and (2) brain tumor images from The Cancer Genome Atlas (TCGA) database. Experimental results confirm that ALSF can offer benefits over state of the art alternatives.


title: Automatic Generation of Bounds for Polynomial Systems with Application to the Lorenz System

abstract: This study covers an analytical approach to calculate positively invariant sets of dynamical systems. Using Lyapunov techniques and quantifier elimination methods, an automatic procedure for determining bounds in the state space as an enclosure of attractors is proposed. The available software tools permit an algorithmizable process, which normally requires a good insight into the systems dynamics and experience. As a result we get an estimation of the attractor, whose conservatism only results from the initial choice of the Lyapunov candidate function. The proposed approach is illustrated on the well-known Lorenz system.


title: A compact broadband terahertz range quarter-wave plate

abstract: We detail the design and characterization of a terahertz range achromatic quarter-wave plate based on a stack of aligned variable thickness birefringent sapphire discs. The disc thicknesses and relative rotations of the discs are determined through a basin-hopping Monte Carlo thermal annealing routine. The basin-hopping scheme allows an improved refinement of the required thicknesses and rotations to give a predicted phase error from the ideal $\pi/2$ of only $0.5 \%$, which is a factor of approximately 6 better than previous efforts. Moreover, the large contrast between extraordinary and ordinary axes of sapphire allow us to greatly decrease the overall optical path length of our wave plate design by approximately a factor of 10 over similar designs based on quartz discs. However, this very same contrast requires very precise tolerances in the milled thicknesses of the discs and their assembly. We detail a method to compensate for differences in the thickness from their calculated ideal values. We have constructed one of our designs and found it similar in performance to other configurations, but with our much more compact geometry.

categories: physics.ins-det cond-mat.mes-hall physics.app-ph physics.optics


title: Block-Diagonal and LT Codes for Distributed Computing With Straggling Servers

abstract: We propose two coded schemes for the distributed computing problem of multiplying a matrix by a set of vectors. The first scheme is based on partitioning the matrix into submatrices and applying maximum distance separable (MDS) codes to each submatrix. For this scheme, we prove that up to a given number of partitions the communication load and the computational delay (not including the encoding and decoding delay) are identical to those of the scheme recently proposed by Li et al., based on a single, long MDS code. However, due to the use of shorter MDS codes, our scheme yields a significantly lower overall computational delay when the delay incurred by encoding and decoding is also considered. We further propose a second coded scheme based on Luby Transform (LT) codes under inactivation decoding. Interestingly, LT codes may reduce the delay over the partitioned scheme at the expense of an increased communication load. We also consider distributed computing under a deadline and show numerically that the proposed schemes outperform other schemes in the literature, with the LT code-based scheme yielding the best performance for the scenarios considered.


title: Squares of Hamiltonian cycles in 3-uniform hypergraphs

abstract: We show that every $3$-uniform hypergraph $H=(V,E)$ with $|V(H)|=n$ and minimum pair degree at least $(4/5+o(1))n$ contains a squared Hamiltonian cycle. This may be regarded as a first step towards a hypergraph version of the P\'osa-Seymour conjecture.


title: Smart, Sparse Contours to Represent and Edit Images

abstract: We study the problem of reconstructing an image from information stored at contour locations. We show that high-quality reconstructions with high fidelity to the source image can be obtained from sparse input, e.g., comprising less than $6\%$ of image pixels. This is a significant improvement over existing contour-based reconstruction methods that require much denser input to capture subtle texture information and to ensure image quality. Our model, based on generative adversarial networks, synthesizes texture and details in regions where no input information is provided. The semantic knowledge encoded into our model and the sparsity of the input allows to use contours as an intuitive interface for semantically-aware image manipulation: local edits in contour domain translate to long-range and coherent changes in pixel space. We can perform complex structural changes such as changing facial expression by simple edits of contours. Our experiments demonstrate that humans as well as a face recognition system mostly cannot distinguish between our reconstructions and the source images.


title: Plasmonic enhancement of aligned semiconducting graphene nanoribbons

abstract: We couple photoluminescent semiconducting 7-atom wide armchair edge graphene nanoribbons to plasmonic nanoantenna arrays and demonstrate an enhancement of the photoluminescence and Raman scattering intensity of the nanoribbons by more than one order of magnitude. The increase in signal allows us to study Raman spectra with high signal-to-noise ratio. Using plasmonic enhancement we are able to detect the off-resonant Raman signals from the modified radial breathing-like mode (RBLM) due to physisorbed molecules, the 3rd order RBLM, and C-H vibrations. We find excellent agreement between data and simulations describing the spectral dependence of the enhancement and modifications of the polarization anisotropy. The strong field gradients in the optical near-field further allow us to probe the subwavelength coherence properties of the phonon modes in the nanoribbons. We theoretically model this considering a finite coherence length along the GNR direction. Our results allow estimating the coherence length in graphene nanoribbons.


title: The Elusive Origin of Mercury

abstract: The MESSENGER mission sought to discover what physical processes determined Mercury's high metal to silicate ratio. Instead, the mission has discovered multiple anomalous characteristics about our innermost planet. The lack of FeO and the reduced oxidation state of Mercury's crust and mantle are more extreme than nearly all other known materials in the solar system. In contrast, moderately volatile elements are present in abundances comparable to the other terrestrial planets. No single process during Mercury's formation is able to explain all of these observations. Here, we review the current ideas for the origin of Mercury's unique features. Gaps in understanding the innermost regions of the solar nebula limit testing different hypotheses. Even so, all proposed models are incomplete and need further development in order to unravel Mercury's remaining secrets.


title: Estimating activity cycles with probabilistic methods I. Bayesian Generalised Lomb-Scargle Periodogram with Trend

abstract: Period estimation is one of the central topics in astronomical time series analysis, where data is often unevenly sampled. Especially challenging are studies of stellar magnetic cycles, as there the periods looked for are of the order of the same length than the datasets themselves. The datasets often contain trends, the origin of which is either a real long-term cycle or an instrumental effect, but these effects cannot be reliably separated, while they can lead to erroneous period determinations if not properly handled. In this study we aim at developing a method that can handle the trends properly, and by performing extensive set of testing, we show that this is the optimal procedure when contrasted with methods that do not include the trend directly to the model. The effect of the form of the noise (whether constant or heteroscedastic) on the results is also investigated. We introduce a Bayesian Generalised Lomb-Scargle Periodogram with Trend (BGLST), which is a probabilistic linear regression model using Gaussian priors for the coefficients and uniform prior for the frequency parameter. We show, using synthetic data, that when there is no prior information on whether and to what extent the true model of the data contains a linear trend, the introduced BGLST method is preferable to the methods which either detrend the data or leave the data untrended before fitting the periodic model. Whether to use noise with different than constant variance in the model depends on the density of the data sampling as well as on the true noise type of the process.


title: Converging Shock Flows for a Mie-Gr\"uneisen Equation of State

abstract: Previous work has shown that the one-dimensional (1D) inviscid compressible flow (Euler) equations admit a wide variety of scale-invariant solutions (including the famous Noh, Sedov, and Guderley shock solutions) when the included equation of state (EOS) closure model assumes a certain scale-invariant form. However, this scale-invariant EOS class does not include even simple models used for shock compression of crystalline solids, including many broadly applicable representations of Mie-Gr\"uneisen EOS. Intuitively, this incompatibility naturally arises from the presence of multiple dimensional scales in the Mie-Gr\"uneisen EOS, which are otherwise absent from scale-invariant models that feature only dimensionless parameters (such as the adiabatic index in the ideal gas EOS). The current work extends previous efforts intended to rectify this inconsistency, by using a scale-invariant EOS model to approximate a Mie- Gr\"uneisen EOS form. To this end, the adiabatic bulk modulus for the Mie-Gr\"uneisen EOS is constructed, and its key features are used to motivate the selection of a scale-invariant approximation form. The remaining surrogate model parameters are selected through enforcement of the Rankine-Hugoniot jump conditions for an infinitely strong shock in a Mie-Gr\"uneisen material. Finally, the approximate EOS is used in conjunction with the 1D inviscid Euler equations to calculate a semi-analytical, Guderley-like imploding shock solution in a metal sphere, and to determine if and when the solution may be valid for the underlying Mie-Gr\"uneisen EOS.


title: Pathwise uniqueness of non-uniformly elliptic SDEs with rough coefficients

abstract: In this paper we review and improve pathwise uniqueness results for some types of one-dimensional stochastic differential equations (SDE) involving the local time of the unknown process. The diffusion coefficient of the SDEs we consider is allowed to vanish on a set of positive measure and is not assumed to be smooth. As opposed to various existing results, our arguments are mainly based on the comparison theorem for local time and the occupation time formula. We apply our pathwise uniqueness results to derive strong existence and other properties of solutions for SDEs with rough coefficients.


title: Interventions over Predictions: Reframing the Ethical Debate for Actuarial Risk Assessment

abstract: Actuarial risk assessments might be unduly perceived as a neutral way to counteract implicit bias and increase the fairness of decisions made at almost every juncture of the criminal justice system, from pretrial release to sentencing, parole and probation. In recent times these assessments have come under increased scrutiny, as critics claim that the statistical techniques underlying them might reproduce existing patterns of discrimination and historical biases that are reflected in the data. Much of this debate is centered around competing notions of fairness and predictive accuracy, resting on the contested use of variables that act as "proxies" for characteristics legally protected against discrimination, such as race and gender. We argue that a core ethical debate surrounding the use of regression in risk assessments is not simply one of bias or accuracy. Rather, it's one of purpose. If machine learning is operationalized merely in the service of predicting individual future crime, then it becomes difficult to break cycles of criminalization that are driven by the iatrogenic effects of the criminal justice system itself. We posit that machine learning should not be used for prediction, but rather to surface covariates that are fed into a causal model for understanding the social, structural and psychological drivers of crime. We propose an alternative application of machine learning and causal inference away from predicting risk scores to risk mitigation.


title: Simulation and stability analysis of oblique shock wave/boundary layer interactions at Mach 5.92

abstract: We investigate flow instability created by an oblique shock wave impinging on a Mach 5.92 laminar boundary layer at a transitional Reynolds number. The adverse pressure gradient of the oblique shock causes the boundary layer to separate from the wall, resulting in the formation of a recirculation bubble. For sufficiently large oblique shock angles, the recirculation bubble is unstable to three-dimensional perturbations and the flow bifurcates from its original laminar state. We utilize Direct Numerical Simulation (DNS) and Global Stability Analysis (GSA) to show that this first occurs at a critical shock angle of $\theta = 12.9^o$. At bifurcation, the least stable global mode is non-oscillatory, and it takes place at a spanwise wavenumber $\beta=0.25$, in good agreement with DNS results. Examination of the critical global mode reveals that it originates from an interaction between small spanwise corrugations at the base of the incident shock, streamwise vortices inside the recirculation bubble, and spanwise modulation of the bubble strength. The global mode drives the formation of long streamwise streaks downstream of the bubble. While the streaks may be amplified by either the lift-up effect or by G\"ortler instability, we show that centrifugal instability plays no role in the upstream self-sustaining mechanism of the global mode. We employ an adjoint solver to corroborate our physical interpretation by showing that the critical global mode is most sensitive to base flow modifications that are entirely contained inside the recirculation bubble.


title: Estimating activity cycles with probabilistic methods II. The Mount Wilson Ca H&K data

abstract: Debate over the existence of branches in the stellar activity-rotation diagrams continues. Application of modern time series analysis tools to study the mean cycle periods in chromospheric activity index is lacking. We develop such models, based on Gaussian processes, for one-dimensional time series and apply it to the extended Mount Wilson Ca H&K sample. Our main aim is to study how the previously commonly used assumption of strict harmonicity of the stellar cycles as well as handling of the linear trends affects the results. We introduce three methods of different complexity, starting with the simple Bayesian harmonic model and followed by Gaussian Process models with periodic and quasi-periodic covariance functions. We confirm the existence of two populations in the activity-period diagram. We find only one significant trend in the inactive population, namely that the cycle periods get shorter with increasing rotation. This is in contrast with earlier studies, that postulate the existence of trends in both of the populations. In terms of rotation to cycle period ratio, our data is consistent with only two activity branches such that the active branch merges together with the transitional one. The retrieved stellar cycles are uniformly distributed over the R'HK activity index, indicating that the operation of stellar large-scale dynamos carries smoothly over the Vaughan-Preston gap. At around the solar activity index, however, indications of a disruption in the cyclic dynamo action are seen. Our study shows that stellar cycle estimates depend significantly on the model applied. Such model-dependent aspects include the improper treatment of linear trends, while the assumption of strict harmonicity can result in the appearance of double cyclicities that seem more likely to be explained by the quasi-periodicity of the cycles.


title: Determination of Boolean models by mean values of mixed volumes

abstract: In Weil (2001) formulas were proved for stationary Boolean models $Z$ in $\mathbb{R}^d$ with convex or polyconvex grains, which express the densities of mixed volumes of $Z$ in terms of related mean values of the underlying Poisson particle process $X$. These formulas were then used to show that in dimensions 2 and 3 the mean values of $Z$ determine the intensity $\gamma$ of $X$. For $d=4$ a corresponding result was also stated, but the proof given was incomplete, since in the formula for the mean Euler characteristic $\overline V_0 (Z)$ a term $\overline V^{(0)}_{2,2}(X,X)$ was missing. This was pointed out in Goodey and Weil (2002), where it was also explained that a new decomposition result for mixed volumes and mixed translative functionals would be needed to complete the proof. Such a general decomposition result is now available based on flag measures of the convex bodies involved (Hug, Rataj and Weil (2013, 2017)). Here, we show that such flag representations not only lead to a correct derivation of the 4-dimensional result but even yield a corresponding uniqueness theorem in all dimensions. In the proof of the latter, we make use of Alesker's representation theorem for translation invariant valuations. We also discuss which shape information can be obtained in this way and comment on the situation in the non-stationary case.


title: Exploring the Lyapunov instability properties of high-dimensional atmospheric and climate models

abstract: The stability properties of intermediate-order climate models are investigated by computing their Lyapunov exponents (LEs). The two models considered are PUMA (Portable University Model of the Atmosphere), a primitive-equation simple general circulation model, and MAOOAM (Modular Arbitrary-Order Ocean-Atmosphere Model), a quasi-geostrophic coupled ocean-atmosphere model on a beta-plane. We wish to investigate the effect of the different levels of filtering on the instabilities and dynamics of the atmospheric flows. Moreover, we assess the impact of the oceanic coupling, the dissipation scheme and the resolution on the spectra of LEs. The PUMA Lyapunov spectrum is computed for two different values of the meridional temperature gradient defining the Newtonian forcing. The increase of the gradient gives rise to a higher baroclinicity and stronger instabilities, corresponding to a larger dimension of the unstable manifold and a larger first LE. The convergence rate of the rate functional for the large deviation law of the finite-time Lyapunov exponents (FTLEs) is fast for all exponents, which can be interpreted as resulting from the absence of a clear-cut atmospheric time-scale separation in such a model. The MAOOAM spectra show that the dominant atmospheric instability is correctly represented even at low resolutions. However, the dynamics of the central manifold, which is mostly associated to the ocean dynamics, is not fully resolved because of its associated long time scales, even at intermediate orders. This paper highlights the need to investigate the natural variability of the atmosphere-ocean coupled dynamics by associating rate of growth and decay of perturbations to the physical modes described using the formalism of the covariant Lyapunov vectors and to consider long integrations in order to disentangle the dynamical processes occurring at all time scales.


title: ConvSCCS: convolutional self-controlled case series model for lagged adverse event detection

abstract: With the increased availability of large databases of electronic health records (EHRs) comes the chance of enhancing health risks screening. Most post-marketing detections of adverse drug reaction (ADR) rely on physicians' spontaneous reports, leading to under reporting. To take up this challenge, we develop a scalable model to estimate the effect of multiple longitudinal features (drug exposures) on a rare longitudinal outcome. Our procedure is based on a conditional Poisson model also known as self-controlled case series (SCCS). We model the intensity of outcomes using a convolution between exposures and step functions, that are penalized using a combination of group-Lasso and total-variation. This approach does not require the specification of precise risk periods, and allows to study in the same model several exposures at the same time. We illustrate the fact that this approach improves the state-of-the-art for the estimation of the relative risks both on simulations and on a cohort of diabetic patients, extracted from the large French national health insurance database (SNIIRAM), a SQL database built around medical reimbursements of more than 65 million people. This work has been done in the context of a research partnership between Ecole Polytechnique and CNAMTS (in charge of SNIIRAM).


title: How Well Can Generative Adversarial Networks Learn Densities: A Nonparametric View

abstract: We study in this paper the rate of convergence for learning densities under the Generative Adversarial Networks (GAN) framework, borrowing insights from nonparametric statistics. We introduce an improved GAN estimator that achieves a faster rate, through simultaneously leveraging the level of smoothness in the target density and the evaluation metric, which in theory remedies the mode collapse problem reported in the literature. A minimax lower bound is constructed to show that when the dimension is large, the exponent in the rate for the new GAN estimator is near optimal. One can view our results as answering in a quantitative way how well GAN learns a wide range of densities with different smoothness properties, under a hierarchy of evaluation metrics. As a byproduct, we also obtain improved generalization bounds for GAN with deeper ReLU discriminator network.


title: Superfluid density and compressibility at the superfluid-Mott glass transition

abstract: Systems of disordered interacting bosons with particle-hole symmetry can undergo a quantum phase transition between the superfluid phase and the Mott glass phase which is a gapless incompressible insulator. We employ large-scale Monte Carlo simulations of a two-dimensional site-diluted quantum rotor model to investigate the properties of the superfluid density and the compressibility at this transition. We find that both quantities feature power-law critical behavior with exponents governed by generalized Josephson relations.


title: On the matchings-Jack and hypermap-Jack conjectures for labelled matchings and star hypermaps

abstract: Introduced by Goulden and Jackson in their 1996 paper, the matchings-Jack conjecture and the hypermap-Jack conjecture (also known as the $b$-conjecture) are two major open questions relating Jack symmetric functions, the representation theory of the symmetric groups and combinatorial maps. They show that the coefficients in the power sum expansion of some Cauchy sum for Jack symmetric functions and in the logarithm of the same sum interpolate respectively between the structure constants of the class algebra and the double coset algebra of the symmetric group and between the numbers of orientable and locally orientable hypermaps. They further provide some evidence that these two families of coefficients indexed by three partitions of a given integer $n$ and the Jack parameter $\alpha$ are polynomials in $\beta = \alpha-1$ with non negative integer coefficients of combinatorial significance. This paper is devoted to the case when one of the three partitions is equal to $(n)$. We exhibit some polynomial properties of both families of coefficients and prove a variation of the hypermap-Jack conjecture and the matchings-Jack conjecture involving labelled hypermaps and matchings in some important cases.


title: Pricing double barrier options on homogeneous diffusions: a Neumann series of Bessel functions representation

abstract: This paper develops a novel analytically tractable Neumann series of Bessel functions representation for pricing (and hedging) European-style double barrier knock-out options, which can be applied to the whole class of one-dimensional time-homogeneous diffusions even for the cases where the corresponding transition density is not known. The proposed numerical method is shown to be efficient and simple to implement. To illustrate the flexibility and computational power of the algorithm, we develop an extended jump to default model that is able to capture several empirical regularities commonly observed in the literature.


title: Explicit Reference Governor for the Constrained Control of Time-Delayed Linear Systems

abstract: This paper introduces an explicit reference governor approach for controlling time delay linear systems subject to state and input constraints. The proposed framework relies on suitable invariant sets that can be built using both Lyapunov-Razumikhin and Lyapunov-Krasovskii arguments. The proposed method is validated both numerically and experimentally using several alternative formulations.


title: Graph Clustering with Dynamic Embedding

abstract: Graph clustering (or community detection) has long drawn enormous attention from the research on web mining and information networks. Recent literature on this topic has reached a consensus that node contents and link structures should be integrated for reliable graph clustering, especially in an unsupervised setting. However, existing methods based on shallow models often suffer from content noise and sparsity. In this work, we propose to utilize deep embedding for graph clustering, motivated by the well-recognized power of neural networks in learning intrinsic content representations. Upon that, we capture the dynamic nature of networks through the principle of influence propagation and calculate the dynamic network embedding. Network clusters are then detected based on the stable state of such an embedding. Unlike most existing embedding methods that are task-agnostic, we simultaneously solve for the underlying node representations and the optimal clustering assignments in an end-to-end manner. To provide more insight, we theoretically analyze our interpretation of network clusters and find its underlying connections with two widely applied approaches for network modeling. Extensive experimental results on six real-world datasets including both social networks and citation networks demonstrate the superiority of our proposed model over the state-of-the-art.


title: ReabsNet: Detecting and Revising Adversarial Examples

abstract: Though deep neural network has hit a huge success in recent studies and applica- tions, it still remains vulnerable to adversarial perturbations which are imperceptible to humans. To address this problem, we propose a novel network called ReabsNet to achieve high classification accuracy in the face of various attacks. The approach is to augment an existing classification network with a guardian network to detect if a sample is natural or has been adversarially perturbed. Critically, instead of simply rejecting adversarial examples, we revise them to get their true labels. We exploit the observation that a sample containing adversarial perturbations has a possibility of returning to its true class after revision. We demonstrate that our ReabsNet outperforms the state-of-the-art defense method under various adversarial attacks.


title: A Homological Definition for the Tate--Shafarevich Group of a Pell Conic

abstract: Franz Lemmermeyer's previous work laid the framework for a description of the arithmetic of Pell conics, which is analogous to that of elliptic curves. He describes a group law on conics and conjectures the existence of an analogous Tate--Shafarevich group with order the squared ideals of the narrow class group. In this article, we provide a cohomological definition of the Tate--Shafarevich group and show that its order is as Lemmermeyer conjectured.


title: Entropy solutions for parabolic problems in Musielak framework without using the sign condition of the nonlinearities

abstract: We prove in this paper the existence of solutions of nonlinear parabolic problems in inhomogeneous Musielak Orlicz Sobolev spaces, we assume neither a $\Delta_2$ nor $\nabla_2$ on the Musielak function $\varphi$. The main contribution of our work is to prove the existence of entropy solutions without the sign condition on the nonlinearity. The second term $f$ belongs to $L^1(Q)$.


title: Core-Collapse Supernova Simulations including Neutrino Interactions from the Virial EOS

abstract: Core-collapse supernova explosions are driven by a central engine that converts a small fraction of the gravitational binding energy released during core collapse to outgoing kinetic energy. The suspected mode for this energy conversion is the neutrino mechanism, where a fraction of the neutrinos emitted from the newly formed protoneutron star are absorbed by and heat the matter behind the supernova shock. Accurate neutrino-matter interaction terms are crucial for simulating these explosions. In this proceedings for IAUS 331, SN 1987A, 30 years later, we explore several corrections to the neutrino-nucleon scattering opacity and demonstrate the effect on the dynamics of the core-collapse supernova central engine via two dimensional neutrino-radiation-hydrodynamics simulations. Our results reveal that the explosion properties are sensitive to corrections to the neutral-current scattering cross section at the 10-20% level, but only for densities at or above $\sim 10^{12}$ g cm$^{-3}$


title: T-count and Qubit Optimized Quantum Circuit Design of the Non-Restoring Square Root Algorithm

abstract: Quantum circuits for basic mathematical functions such as the square root are required to implement scientific computing algorithms on quantum computers. Quantum circuits that are based on Clifford+T gates can easily be made fault tolerant but the T gate is very costly to implement. As a result, reducing T-count has become an important optimization goal. Further, quantum circuits with many qubits are difficult to realize, making designs that save qubits and produce no garbage outputs desirable. In this work, we present a T-count optimized quantum square root circuit with only $2 \cdot n +1$ qubits and no garbage output. To make a fair comparison against existing work, the Bennett's garbage removal scheme is used to remove garbage output from existing works. We determined that our proposed design achieves an average T-count savings of $43.44 \%$, $98.95 \%$, $41.06 \%$ and $20.28 \%$ as well as qubit savings of $85.46 \%$, $95.16 \%$, $90.59 \%$ and $86.77 \%$ compared to existing works.


title: On embeddings of locally finite metric spaces into $\ell_p$

abstract: It is known that if finite subsets of a locally finite metric space $M$ admit $C$-bilipschitz embeddings into $\ell_p$ $(1\le p\le \infty)$, then for every $\epsilon>0$, the space $M$ admits a $(C+\epsilon)$-bilipschitz embedding into $\ell_p$. The goal of this paper is to show that for $p\ne 2,\infty$ this result is sharp in the sense that $\epsilon$ cannot be dropped out of its statement.


title: Optimal control for one-qubit quantum sensing

abstract: Quantum systems can be exquisite sensors thanks to their sensitivity to external perturbations. This same characteristic also makes them fragile to external noise. Quantum control can tackle the challenge of protecting quantum sensors from environmental noise, while leaving their strong coupling to the target field to be measured. As the compromise between these two conflicting requirements does not always have an intuitive solution, optimal control based on numerical search could prove very effective. Here we adapt optimal control theory to the quantum sensing scenario, by introducing a cost function that, unlike the usual fidelity of operation, correctly takes into account both the unknown field to be measured and the environmental noise. We experimentally implement this novel control paradigm using a Nitrogen Vacancy center in diamond, finding improved sensitivity to a broad set of time varying fields. The demonstrated robustness and efficiency of the numerical optimization, as well as the sensitivity advantaged it bestows, will prove beneficial to many quantum sensing applications.


title: Multi-color, femtosecond $\gamma$-ray pulse trains driven by comb-like electron beams

abstract: Photon engineering can be exploited to control the nonlinear evolution of the drive pulse in a laser-plasma accelerator (LPA), offering new avenues to tailor electron beam phase space on a femtosecond time scale. One promising option is to drive an LPA with an incoherent stack of two sub-Joule, multi-TW pulses of different colors. Slow self-compression of the bi-color optical driver delays electron dephasing, boosting electron beam energy without accumulation of a massive low-energy tail. The modest energy of the stack affords kHz-scale repetition rate at manageable laser average power. Propagating the stack in a preformed plasma channel induces periodic self-focusing in the trailing pulse, causing oscillations in the size of accelerating bucket. The resulting periodic injection generates, over a mm-scale distance, a train of GeV-scale electron bunches with 5D brightness exceeding $10^{17}$ A/m$^2$. This unconventional comb-like beam, with femtosecond synchronization and controllable energy spacing of components, emits, via Thomson scattering, a train of highly collimated gigawatt $\gamma$-ray pulses. Each pulse, corresponding to a distinct energy band between 2.5 and 25 MeV, contains over $10^6$ photons.


title: Finite collineation groups and birational rigidity

abstract: We classify finite groups $G$ in $\mathrm{PGL}_{4}(\mathbb{C})$ such that $\mathbb{P}^3$ is $G$-birationally rigid.


title: Linear centralization classifier

abstract: A classification algorithm, called the Linear Centralization Classifier (LCC), is introduced. The algorithm seeks to find a transformation that best maps instances from the feature space to a space where they concentrate towards the center of their own classes, while maximimizing the distance between class centers. We formulate the classifier as a quadratic program with quadratic constraints. We then simplify this formulation to a linear program that can be solved effectively using a linear programming solver (e.g., simplex-dual). We extend the formulation for LCC to enable the use of kernel functions for non-linear classification applications. We compare our method with two standard classification methods (support vector machine and linear discriminant analysis) and four state-of-the-art classification methods when they are applied to eight standard classification datasets. Our experimental results show that LCC is able to classify instances more accurately (based on the area under the receiver operating characteristic) in comparison to other tested methods on the chosen datasets. We also report the results for LCC with a particular kernel to solve for synthetic non-linear classification problems.


title: Quantum harmonic oscillator with time dependent mass

abstract: We use the Fourier operator to transform a time dependent mass quantum harmonic oscillator into a frequency dependent one. Then we use Lewis-Ermakov invariants to solve the Schr\"odinger equation by using squeeze operators. Finally we give two examples of time dependencies: quadratically and hyperbolically growing masses.


title: From arteries to boreholes: Transient response of a poroelastic cylinder to fluid injection

abstract: The radially outward flow of fluid through a porous medium occurs in many practical problems, from transport across vascular walls to the pressurisation of boreholes in the subsurface. When the driving pressure is non-negligible relative to the stiffness of the solid structure, the poromechanical coupling between the fluid and the solid can control both the steady-state and the transient mechanics of the system. Very large pressures or very soft materials lead to large deformations of the solid skeleton, which introduce kinematic and constitutive nonlinearity that can have a nontrivial impact on these mechanics. Here, we study the transient response of a poroelastic cylinder to sudden fluid injection. We consider the impacts of kinematic and constitutive nonlinearity, both separately and in combination, and we highlight the central role of driving method in the evolution of the response. We show that the various facets of nonlinearity may either accelerate or decelerate the transient response relative to linear poroelasticity, depending on the boundary conditions and the initial geometry, and that an imposed fluid pressure leads to a much faster response than an imposed fluid flux.


title: Wave kinetics of drift-wave turbulence and zonal flows beyond the ray approximation

abstract: Inhomogeneous drift-wave turbulence can be modeled as an effective plasma where drift waves act as quantumlike particles and the zonal-flow velocity serves as a collective field through which they interact. This effective plasma can be described by a Wigner-Moyal equation (WME), which generalizes the quasilinear wave-kinetic equation (WKE) to the full-wave regime, i.e., resolves the wavelength scale. Unlike waves governed by manifestly quantumlike equations, whose WMEs can be borrowed from quantum mechanics and are commonly known, drift waves have Hamiltonians very different from those of conventional quantum particles. This causes unusual phase-space dynamics that is typically not captured by the WKE. We demonstrate how to correctly model this dynamics with the WME instead. Specifically, we report full-wave phase-space simulations of the zonal-flow formation (zonostrophic instability), deterioration (tertiary instability), and the so-called predator-prey oscillations. We also show how the WME facilitates analysis of these phenomena, namely, (i) we show that full-wave effects critically affect the zonostrophic instability, particularly, its nonlinear stage and saturation; (ii) we derive the tertiary-instability growth rate; and (iii) we demonstrate that, with full-wave effects retained, the predator-prey oscillations do not require zonal-flow collisional damping, contrary to previous studies. We also show how the famous Rayleigh-Kuo criterion, which has been missing in wave-kinetic theories of drift-wave turbulence, emerges from the WME.


title: Using LIP to Gloss Over Faces in Single-Stage Face Detection Networks

abstract: This work shows that it is possible to fool/attack recent state-of-the-art face detectors which are based on the single-stage networks. Successfully attacking face detectors could be a serious malware vulnerability when deploying a smart surveillance system utilizing face detectors. We show that existing adversarial perturbation methods are not effective to perform such an attack, especially when there are multiple faces in the input image. This is because the adversarial perturbation specifically generated for one face may disrupt the adversarial perturbation for another face. In this paper, we call this problem the Instance Perturbation Interference (IPI) problem. This IPI problem is addressed by studying the relationship between the deep neural network receptive field and the adversarial perturbation. As such, we propose the Localized Instance Perturbation (LIP) that uses adversarial perturbation constrained to the Effective Receptive Field (ERF) of a target to perform the attack. Experiment results show the LIP method massively outperforms existing adversarial perturbation generation methods -- often by a factor of 2 to 10.


title: Existence of solutions for critical Choquard equations via the concentration compactness method

abstract: In this paper we consider the nonlinear Choquard equation $$ -\Delta u+V(x)u =\left(\int_{\mathbb{R}^N}\frac{G(y,u)}{|x-y|^{\mu}}dy\right)g(x,u)\hspace{4.14mm}\mbox{in}\hspace{1.14mm} \mathbb{R}^N, $$ where $0<\mu<N$, $N\geq3$, $g(x,u)$ is of critical growth due to the Hardy--Littlewood--Sobolev inequality and $G(x,u)=\displaystyle\int^u_0g(x,s)ds$. Firstly, by assuming that the potential $V(x)$ might be sign-changing, we study the existence of Mountain-Pass solution via a concentration-compactness principle for the Choquard equation. Secondly, under the conditions introduced by Benci and Cerami \cite{BC1}, we also study the existence of high energy solution by using a global compactness lemma for the nonlocal Choquard equation.


title: Computational Modelling of Aquaporin Co-regulation in Cancer

abstract: A computational model of aquaporin regulation in cancer cells has been constructed as a Qualitative Network in the software BioModelAnalyzer (BMA). The model connects some important aquaporins expressed in human cancer to common phenotypes via a number of fundamental, dysregulated signalling pathways. Based on over 60 publications, this model can not only reproduce the results reported in a discrete, qualitative manner, but also reconcile the seemingly incompatible phenotype with research consensus by suggesting molecular mechanisms accountable for it. Novel predictions have also been made by mimicking real-life experiments in the model.


title: Federated Control with Hierarchical Multi-Agent Deep Reinforcement Learning

abstract: We present a framework combining hierarchical and multi-agent deep reinforcement learning approaches to solve coordination problems among a multitude of agents using a semi-decentralized model. The framework extends the multi-agent learning setup by introducing a meta-controller that guides the communication between agent pairs, enabling agents to focus on communicating with only one other agent at any step. This hierarchical decomposition of the task allows for efficient exploration to learn policies that identify globally optimal solutions even as the number of collaborating agents increases. We show promising initial experimental results on a simulated distributed scheduling problem.


title: A New Transferable Interatomic Potential for Molecular Dynamics Simulations of Borosilicate Glasses

abstract: Borosilicate glasses are traditionally challenging to model using atomic scale simulations due to the composition and thermal history dependence of the coordination state of B atoms. Here, we report a new empirical interatomic potential that shows a good transferability over a wide range of borosilicate glasses--ranging from pure silicate to pure borate end members--while relying on a simple formulation and a constant set of energy parameters. In particular, we show that our new potential accurately predicts the compositional dependence of the average coordination number of boron atoms, glass density, overall short-range and medium-range order structure, and shear viscosity values for several borosilicate glasses and liquids. This suggests that our new potential could be used to gain new insights into the structure of a variety of advanced borosilicate glasses to help elucidate composition-structure-property relationships--including in complex nuclear waste immobilization glasses.


title: Beyond saliency: understanding convolutional neural networks from saliency prediction on layer-wise relevance propagation

abstract: Despite the tremendous achievements of deep convolutional neural networks (CNNs) in many computer vision tasks, understanding how they actually work remains a significant challenge. In this paper, we propose a novel two-step understanding method, namely Salient Relevance (SR) map, which aims to shed light on how deep CNNs recognize images and learn features from areas, referred to as attention areas, therein. Our proposed method starts out with a layer-wise relevance propagation (LRP) step which estimates a pixel-wise relevance map over the input image. Following, we construct a context-aware saliency map, SR map, from the LRP-generated map which predicts areas close to the foci of attention instead of isolated pixels that LRP reveals. In human visual system, information of regions is more important than of pixels in recognition. Consequently, our proposed approach closely simulates human recognition. Experimental results using the ILSVRC2012 validation dataset in conjunction with two well-established deep CNN models, AlexNet and VGG-16, clearly demonstrate that our proposed approach concisely identifies not only key pixels but also attention areas that contribute to the underlying neural network's comprehension of the given images. As such, our proposed SR map constitutes a convenient visual interface which unveils the visual attention of the network and reveals which type of objects the model has learned to recognize after training. The source code is available at https://github.com/Hey1Li/Salient-Relevance-Propagation.