title: Jetting Through The Primordial Universe

abstract: Collisions of heavy ion nuclei at relativistic speeds (close to the speed of light) creates a high temperature and very dense form of matter, now known to consist of de-confined quarks and gluons, named the quark gluon plasma (QGP). In this thesis, Run1 experimental data from pp and heavy ion collisions at the CERN LHC is analyzed with the CMS detector. The pp jet cross section is compared with next to leading order theoretical calculations supplemented with non perturbative corrections for three different jet radii highlighting better comparisons for larger radii jets. Measurement of the jet yield followed by the nuclear modification factors in proton-lead at 5.02 TeV and lead-lead collisions at 2.76 TeV are presented. A new data driven technique is introduced to estimate and correct for the fake jet contribution in PbPb for low transverse momenta jets. The nuclear modification factors studied in this thesis show jet quenching to be attributed to final state effects, have a strong correlation to the event centrality, a weak inverse correlation to the jet transverse momenta and an apparent independence on the jet radii in the kinematic range studied. These measurements are compared with leading theoretical model calculations and other experimental results at the LHC leading to unanimous agreement on the qualitative nature of jet quenching. This thesis also features novel updates to the Monte Carlo heavy ion event generator JEWEL (Jet Evolution With Energy Loss) including the boson-jet production channels and also background subtraction techniques to reduce the effect of the thermal background. Keeping track of these jet-medium recoils in JEWEL due to the background subtraction techniques significantly improves its descriptions of several jet structure and sub-structure measurements at the LHC. [Shortened abstract]


title: Data Driven Stability Analysis of Black-box Switched Linear Systems

abstract: Can we conclude the stability of an unknown dynamical system from the knowledge of a finite number of snapshots of trajectories? We tackle this black-box problem for switched linear systems. We show that, for any given random set of observations, one can give probabilistic stability guarantees. The probabilistic nature of these guarantees implies a trade-off between their quality and the desired level of confidence. We provide an explicit way of computing the best stability-like guarantee, as a function of both the number of observations and the required level of confidence. Our proof techniques rely on geometrical analysis, chance-constrained optimization, and stability analysis tools for switched systems, including the joint spectral radius.


title: Bit-Tactical: Exploiting Ineffectual Computations in Convolutional Neural Networks: Which, Why, and How

abstract: We show that, during inference with Convolutional Neural Networks (CNNs), more than 2x to $8x ineffectual work can be exposed if instead of targeting those weights and activations that are zero, we target different combinations of value stream properties. We demonstrate a practical application with Bit-Tactical (TCL), a hardware accelerator which exploits weight sparsity, per layer precision variability and dynamic fine-grain precision reduction for activations, and optionally the naturally occurring sparse effectual bit content of activations to improve performance and energy efficiency. TCL benefits both sparse and dense CNNs, natively supports both convolutional and fully-connected layers, and exploits properties of all activations to reduce storage, communication, and computation demands. While TCL does not require changes to the CNN to deliver benefits, it does reward any technique that would amplify any of the aforementioned weight and activation value properties. Compared to an equivalent data-parallel accelerator for dense CNNs, TCLp, a variant of TCL improves performance by 5.05x and is 2.98x more energy efficient while requiring 22% more area.


title: Three colour bipartite Ramsey number of cycles and paths

abstract: The $k$-colour bipartite Ramsey number of a bipartite graph $H$ is the least integer $n$ for which every $k$-edge-coloured complete bipartite graph $K_{n,n}$ contains a monochromatic copy of $H$. The study of bipartite Ramsey numbers was initiated, over 40 years ago, by Faudree and Schelp and, independently, by Gy\'arf\'as and Lehel, who determined the $2$-colour Ramsey number of paths. In this paper we determine asymptotically the $3$-colour bipartite Ramsey number of paths and (even) cycles.


title: Results from phase 1 of the HAYSTAC microwave cavity axion experiment

abstract: We report on the results from a search for dark matter axions with the HAYSTAC experiment using a microwave cavity detector at frequencies between 5.6-5.8$\, \rm Ghz$. We exclude axion models with two photon coupling $g_{a\gamma\gamma}\,\gtrsim\,2\times10^{-14}\,\rm GeV^{-1}$, a factor of 2.7 above the benchmark KSVZ model over the mass range 23.15$\,<\,$$m_a \,$<$\,$24.0$\,\mu\rm eV$. This doubles the range reported in our previous paper. We achieve a near-quantum-limited sensitivity by operating at a temperature $T<h\nu/2k_B$ and incorporating a Josephson parametric amplifier (JPA), with improvements in the cooling of the cavity further reducing the experiment's system noise temperature to only twice the Standard Quantum Limit at its operational frequency, an order of magnitude better than any other dark matter microwave cavity experiment to date. This result concludes the first phase of the HAYSTAC program utilizing a conventional copper cavity and a single JPA.


title: Cascades and Dissipative Anomalies in Nearly Collisionless Plasma Turbulence

abstract: We develop first-principles theory of kinetic plasma turbulence governed by the Vlasov-Maxwell-Landau equations in the limit of vanishing collision rates. Following an exact renormalization-group approach pioneered by Onsager, we demonstrate the existence of a "collisionless range" of scales (lengths and velocities) in 1-particle phase space where the ideal Vlasov-Maxwell equations are satisfied in a "coarse-grained sense". Entropy conservation may nevertheless be violated in that range by a "dissipative anomaly" due to nonlinear entropy cascade. We derive "4/5th-law" type expressions for the entropy flux, which allow us to characterize the singularities (structure-function scaling exponents) required for its non-vanishing. Conservation laws of mass, momentum and energy are not afflicted with anomalous transfers in the collisionless limit. In a subsequent limit of small gyroradii, however, anomalous contributions to inertial-range energy balance may appear due both to cascade of bulk energy and to turbulent redistribution of internal energy in phase space. In that same limit the "generalized Ohm's law" derived from the particle momentum balances reduces to an "ideal Ohm's law", but only in a coarse-grained sense that does not imply magnetic flux-freezing and that permits magnetic reconnection at all inertial-range scales. We compare our results with prior theory based on the gyrokinetic (high gyro-frequency) limit, with numerical simulations, and with spacecraft measurements of the solar wind and terrestrial magnetosphere.


title: On the information in spike timing: neural codes derived from polychronous groups

abstract: There is growing evidence regarding the importance of spike timing in neural information processing, with even a small number of spikes carrying information, but computational models lag significantly behind those for rate coding. Experimental evidence on neuronal behavior is consistent with the dynamical and state dependent behavior provided by recurrent connections. This motivates the minimalistic abstraction investigated in this paper, aimed at providing insight into information encoding in spike timing via recurrent connections. We employ information-theoretic techniques for a simple reservoir model which encodes input spatiotemporal patterns into a sparse neural code, translating the polychronous groups introduced by Izhikevich into codewords on which we can perform standard vector operations. We show that the distance properties of the code are similar to those for (optimal) random codes. In particular, the code meets benchmarks associated with both linear classification and capacity, with the latter scaling exponentially with reservoir size.


title: Laboratory Measurements of X-Ray Emission from Highly Charged Argon Ions

abstract: Uncertainties in atomic models will introduce noticeable additional systematics in calculating the flux of weak dielectronic recombination (DR) satellite lines, affecting the detection and flux measurements of other weak spectral lines. One important example is the Ar XVII He-beta DR, which is expected to be present in emission from the hot intracluster medium (ICM) of galaxy clusters and could impact measurements of the flux of the 3.5 keV line that has been suggested as a secondary emission from a dark matter interaction. We perform a set of experiments using the Lawrence Livermore National Laboratory's electron beam ion trap (EBIT-I) and the X-Ray Spectrometer quantum calorimeter (XRS/EBIT), to test the Ar XVII He-beta DR origin of the 3.5 keV line. We measured the X-ray emission following resonant DR onto helium-like and lithium-like Argon using EBIT-I's Maxwellian simulator mode at a simulated electron temperature of Te=1.74 keV. The measured flux of the Ar XVII He-beta DR lined is too weak to account for the flux in the 3.5 keV line assuming reasonable plasma parameters. We, therefore, rule out Ar XVII He-beta DR as a significant contributor to the 3.5 keV line. A comprehensive comparison between the atomic theory and the EBIT experiment results is also provided.


title: Set Theory-Based Safety Supervisory Control for Wind Turbines to Ensure Adequate Frequency Response

abstract: Inadequate frequency response can arise due to a high penetration of wind turbine generators (WTGs) and requires a frequency support function to be integrated in the WTG. The appropriate design for these controllers to ensure adequate response has not been investigated thoroughly. In this paper, a safety supervisory control (SSC) is proposed to synthesize the supportive modes in WTGs to guarantee performance. The concept, region of safety (ROS), is stated for safe switching synthesis. An optimization formula is proposed to calculate the largest ROS. By assuming a polynomial structure, the problem can be solved by a sum of squares program. A feasible result will generate a polynomial, the zero sublevel set of which represents the ROS and is employed as the safety supervisor. A decentralized communication architecture is proposed for small-scale systems. Moreover, a scheduling loop is suggested so that the supervisor updates its boundary with respect to the renewable penetration level to be robust with respect to variations in system inertia. The proposed controller is first verified on a single-machine three-phase nonlinear microgrid, and then implemented on the IEEE 39-bus system. Both results indicate that the proposed framework and control configuration can guarantee adequate response without excessive conservativeness.


title: Markov chains under nonlinear expectation

abstract: In this paper, we consider continuous-time Markov chains with a finite state space under nonlinear expectations. We define so-called Q-operators as an extension of Q-matrices or rate matrices to a nonlinear setup, where the nonlinearity is due to model uncertainty. The main result gives a full characterization of convex Q-operators in terms of a positive maximum principle, a dual representation by means of Q-matrices, continuous-time Markov chains under convex expectations and nonlinear ordinary differential equations. This extends a classical characterization of generators of Markov chains to the case of model uncertainty in the generator. We further derive a primal and dual representation of the convex semigroup arising from a Markov chain under a convex expectation via the Fenchel-Legendre transformation of its generator. We illustrate the results with several numerical examples, where we compute price bounds for European contingent claims under model uncertainty in terms of the rate matrix.


title: Minimum $T$-Joins and Signed-Circuit Covering

abstract: Let $G$ be a graph and $T$ be a vertex subset of $G$ with even cardinality. A $T$-join of $G$ is a subset $J$ of edges such that a vertex of $G$ is incident with an odd number of edges in $J$ if and only if the vertex belongs to $T$. Minimum $T$-joins have many applications in combinatorial optimizations. In this paper, we show that a minimum $T$-join of a connected graph $G$ has at most $|E(G)|-\frac 1 2 |E(\widehat{\, G\,})|$ edges where $\widehat{\,G\,}$ is the maximum bidegeless subgraph of $G$. Further, we are able to use this result to show that every flow-admissible signed graph $(G,\sigma)$ has a signed-circuit cover with length at most $\frac{19} 6 |E(G)|$. Particularly, a 2-edge-connected signed graph $(G,\sigma)$ with even negativeness has a signed-circuit cover with length at most $\frac 8 3 |E(G)|$.


title: Community Interaction and Conflict on the Web

abstract: Users organize themselves into communities on web platforms. These communities can interact with one another, often leading to conflicts and toxic interactions. However, little is known about the mechanisms of interactions between communities and how they impact users. Here we study intercommunity interactions across 36,000 communities on Reddit, examining cases where users of one community are mobilized by negative sentiment to comment in another community. We show that such conflicts tend to be initiated by a handful of communities---less than 1% of communities start 74% of conflicts. While conflicts tend to be initiated by highly active community members, they are carried out by significantly less active members. We find that conflicts are marked by formation of echo chambers, where users primarily talk to other users from their own community. In the long-term, conflicts have adverse effects and reduce the overall activity of users in the targeted communities. Our analysis of user interactions also suggests strategies for mitigating the negative impact of conflicts---such as increasing direct engagement between attackers and defenders. Further, we accurately predict whether a conflict will occur by creating a novel LSTM model that combines graph embeddings, user, community, and text features. This model can be used toreate early-warning systems for community moderators to prevent conflicts. Altogether, this work presents a data-driven view of community interactions and conflict, and paves the way towards healthier online communities.


title: Almost prime values of the order of abelian varieties over finite fields

abstract: Let $E/\mathbb Q$ be an elliptic curve, and denote by $N(p)$ the number of $\mathbb{F}_p$-points of the reduction modulo $p$ of $E$. A conjecture of Koblitz, refined by Zywina, states that the number of primes $p \leq X$ at which $N(p)$ is also prime is asymptotic to $C_E \cdot X / \log(X)^2$, where $C_E$ is an arithmetically-defined non-negative constant. Following Miri-Murty (2001) and others, Y.R. Liu (2006) and David-Wu (2012) study the number of prime factors of $N(p)$. We generalize their arguments to abelian varieties $A / \mathbb Q$ whose adelic Galois representation has open image in $\textrm{GSp}_{2g} \widehat{\mathbb Z}$. Our main result, after David-Wu, finds a conditional lower bound on the number of primes at which $\# A_p ( \mathbb{F}_p)$ has few prime factors. We also present some experimental evidence in favor of a generalization of Koblitz's conjecture to this context.


title: Ideal convergent subseries in Banach spaces

abstract: Assume that $\mathcal{I}$ is an ideal on $\mathbb{N}$, and $\sum_n x_n$ is a divergent series in a Banach space $X$. We study the Baire category, and the measure of the set $A(\mathcal{I}):=\left\{t \in \{0,1\}^{\mathbb{N}} \colon \sum_n t(n)x_n \textrm{ is } \mathcal{I}\textrm{-convergent}\right\}$. In the category case, we assume that $\mathcal{I}$ has the Baire property and $\sum_n x_n$ is not unconditionally convergent, and we deduce that $A(\mathcal{I})$ is meager. We also study the smallness of $A(\mathcal{I})$ in the measure case when the Haar probability measure $\lambda$ on $\{0,1\}^{\mathbb{N}}$ is considered. If $\mathcal{I}$ is analytic or coanalytic, and $\sum_n x_n$ is $\mathcal{I}$-divergent, then $\lambda(A(\mathcal{I}))=0$ which extends the theorem of Dindo\v{s}, \v{S}al\'at and Toma. Generalizing one of their examples, we show that, for every ideal $\mathcal{I}$ on $\mathbb{N}$, with the property of long intervals, there is a divergent series of reals such that $\lambda(A(Fin))=0$ and $\lambda(A(\mathcal{I}))=1$.


title: Semimartingales and Shrinkage of Filtration

abstract: We consider a complete probability space $(\Omega,\mathcal{F},\mathbb{P})$, which is endowed with two filtrations, $\mathbb{G}$ and $\mathbb{F}$, assumed to satisfy the usual conditions and such that $\mathbb{F} \subset \mathbb{G}$. On this probability space we consider a real valued special $\mathbb{G}$-semimartingale $X$. The purpose of this work is to study the following two problems: A. If $X$ is $\mathbb{F}$-adapted, compute the $\mathbb{F}$-semimartingale characteristics of $X$ in terms of the $\mathbb{G}$-semimartingale characteristics of $X$. B. If $X$ is not $\mathbb{F}$-adapted, given that the $\mathbb{F}$-optional projection of $X$ is a special semimartingale, compute the $\mathbb{F}$-semimartingale characteristics of $\mathbb{F}$-optional projection of $X$ in terms of the $\mathbb{G}$-canonical decomposition and $\mathbb{G}$-semimartingale characteristics of $X$.


title: Biharmonic constant mean curvature surfaces in Killing submersions

abstract: A $3$-dimensional Riemannian manifold is called Killing submersion if it admits a Riemannian submersion over a surface such that its fibers are the trajectories of a complete unit Killing vector field. In this paper, we give a characterization of proper biharmonic CMC surfaces in a Killing submersion. In the last part, we also classify the proper biharmonic Hopf cylinders in a Killing submersion.


title: Orbifold Vertex Operator Algebras and the Positivity Condition

abstract: In this note we show that the irreducible twisted modules of a holomorphic, $C_2$-cofinite vertex operator algebra $V$ have $L_0$-weights at least as large as the smallest $L_0$-weight of $V$. Hence, if $V$ is of CFT-type, then the twisted $V$-modules are almost strictly positively graded. This in turn implies that the fixed-point vertex operator subalgebra $V^G$ for a finite, solvable group of automorphisms of $V$ almost satisfies the positivity condition. These and some further results are obtained by a careful analysis of Dong, Li and Mason's twisted modular invariance.


title: Gravitational Coupling and the Cosmological Constant

abstract: We deal with a dynamical mechanism in which a large cosmological constant, as suggested by inflationary scenarios, decays due to expansion of the universe. This mechanism has its origin in the gravitational coupling of the vacuum density. We assume that the vacuum couples anomalously to gravity that is the metric tensor that appears the gravitational part is not the same as that appears the matter part as suggested by weak equivalence principle. Instead, the two metric tensors are taken to be conformally related. We show that this provides a dynamical mechanism which works during expansion of the universe. We also consider some observational consequences of such a gravitational model.


title: Edge-decomposing graphs into coprime forests

abstract: The Barat-Thomassen conjecture, recently proved in [Bensmail et al.: A proof of the Barat-Thomassen conjecture. J. Combin. Theory Ser. B, 124:39-55, 2017.], asserts that for every tree T, there is a constant $c_T$ such that every $c_T$-edge connected graph G with number of edges (size) divisible by the size of T admits an edge partition into copies of T (a T-decomposition). In this paper, we investigate in which case the connectivity requirement can be dropped to a minimum degree condition. For instance, it was shown in [Bensmail et al.: Edge-partitioning a graph into paths: beyond the Barat-Thomassen conjecture. arXiv:1507.08208] that when T is a path with k edges, there is a constant $d_k$ such that every 24-edge connected graph G with size divisible by k and minimum degree $d_k$ has a T-decomposition. We show in this paper that when F is a coprime forest (the sizes of its components being a coprime set of integers), any graph G with sufficiently large minimum degree has an F-decomposition provided that the size of F divides the size of G (no connectivity is required). A natural conjecture asked in [Bensmail et al.: Edge-partitioning a graph into paths: beyond the Barat-Thomassen conjecture. arXiv:1507.08208] asserts that for a fixed tree T, any graph G of size divisible by the size of T with sufficiently high minimum degree has a T-decomposition, provided that G is sufficiently highly connected in terms of the maximal degree of T. The case of maximum degree 2 is answered by paths. We provide a counterexample to this conjecture in the case of maximum degree 3.


title: Geodesic Obstacle Representation of Graphs

abstract: An obstacle representation of a graph is a mapping of the vertices onto points in the plane and a set of connected regions of the plane (called obstacles) such that the straight-line segment connecting the points corresponding to two vertices does not intersect any obstacles if and only if the vertices are adjacent in the graph. The obstacle representation and its plane variant (in which the resulting representation is a plane straight-line embedding of the graph) have been extensively studied with the main objective of minimizing the number of obstacles. Recently, Biedl and Mehrabi (GD 2017) studied grid obstacle representations of graphs in which the vertices of the graph are mapped onto the points in the plane while the straight-line segments representing the adjacency between the vertices is replaced by the $L_1$ (Manhattan) shortest paths in the plane that avoid obstacles. In this paper, we introduce the notion of geodesic obstacle representations of graphs with the main goal of providing a generalized model, which comes naturally when viewing line segments as shortest paths in the Euclidean plane. To this end, we extend the definition of obstacle representation by allowing some obstacles-avoiding shortest path between the corresponding points in the underlying metric space whenever the vertices are adjacent in the graph. We consider both general and plane variants of geodesic obstacle representations (in a similar sense to obstacle representations) under any polyhedral distance function in $\mathbb{R}^d$ as well as shortest path distances in graphs. Our results generalize and unify the notions of obstacle representations, plane obstacle representations and grid obstacle representations, leading to a number of questions on such embeddings.


title: Electrode configuration and electrical dissipation of mechanical energy in quartz crystal resonators

abstract: Mechanical resonators made with monolithic piezoelectric quartz crystals are promising for studying new physical phenomena. High mechanical quality factors ($Q$) exhibited by the mm-sized quartz resonators make them ideal for studying weak couplings or long timescales in the quantum regime. However, energy losses through mechanical supports pose a serious limiting factor for obtaining high quality factors. Here we investigate how the $Q$ of quartz resonators at deep cryogenic temperatures can be limited by several types of losses related to anchoring. We first introduce means to reduce the mechanical losses by more than an order of magnitude in a no-clamping scheme, obtaining $Q$-factors of $10^8$ of the lowest shear mode. We can exclude a wide coverage of aluminum metallization on the disk or bond wires as sources of dissipation. However, we find a dramatic reduction of the $Q$-factor accompanying an electrode configuration that involves strong focusing of the vibrations in the disk center. We propose a circuit model that accounts for the reduced mechanical $Q$-factor in terms of electrical losses. In particular, we show how the limiting factor for losses can be small ohmic dissipation in a grounding connection, which can be interpreted as electrical anchor losses of the mechanical device.


title: Quasi-Equilibrium Problems with Non-self Constraint Map

abstract: In 2016 Aussel, Sultana and Vetrivel developed the concept of projected solution for quasi-variational inequality problems and projected Nash equilibrium. We introduce a new concept of solution for quasi-equilibrium problems and we study the existence of such solutions. Additionally, as a consequence of our results, we give existence results of projected solutions for quasi-optimization problems, quasi-variational inequalities problems and generalized Nash equilibrium problems.


title: Computational Complexity of Generalized Push Fight

abstract: We analyze the computational complexity of optimally playing the two-player board game Push Fight, generalized to an arbitrary board and number of pieces. We prove that the game is PSPACE-hard to decide who will win from a given position, even for simple (almost rectangular) hole-free boards. We also analyze the mate-in-1 problem: can the player win in a single turn? One turn in Push Fight consists of up to two "moves" followed by a mandatory "push". With these rules, or generalizing the number of allowed moves to any constant, we show mate-in-1 can be solved in polynomial time. If, however, the number of moves per turn is part of the input, the problem becomes NP-complete. On the other hand, without any limit on the number of moves per turn, the problem becomes polynomially solvable again.


title: Nonadiabatic corrections to electric current in molecular junction due to nuclear motion at the molecule-electrode interfaces

abstract: We present quantum electron transport theory that incorporates dynamical effects of motion of atoms on electrode-molecule interfaces in the calculations of the electric current. The theory is based on non-equilibrium Green's functions. We separate time scales in the Green's functions on fast relative time and slow central time. The derivative with respect to the central time serves as a small parameter in the theory. We solve the real-time Kadanoff-Baym equations for molecular Green's functions using Wigner representation and keep terms up to the second order with respect to the central time derivatives. Molecular Green's functions and consequently the electric current are expressed as functions of molecular junction coordinates as well as velocities and accelerations of molecule-electrode interface nuclei. We apply the theory to model a molecular system and study the effects of non-adiabatic nuclear motion on molecular junction conductivity.


title: Label-free imaging of cholesterol and lipid distributions in model membranes

abstract: Over recent decades, lipid membranes have become standard models for examining the biophysics and biochemistry of cell membranes. Interrogation of lipid domains within biomembranes is generally done with fluorescence microscopy via exogenous chemical probes. However, most fluorophores have limited partitioning tunability, with the majority segregating in the least biologically relevant domains (i.e., low-density liquid domains). Therefore, a molecular-level picture of the majority of non-labeled lipids forming the membrane is still elusive. Here, we present simple, label-free imaging of domain formation in lipid monolayers, with chemical selectivity in unraveling lipid and cholesterol composition in all domain types. Exploiting conventional vibrational contrast in spontaneous Raman imaging, combined with chemometrics analysis, allows for examination of ternary systems containing saturated lipids, unsaturated lipids, and cholesterol. We confirm features commonly observed by fluorescence microscopy, and provide an unprecedented analysis of cholesterol distribution at the single-membrane level.


title: Local Kernels that Approximate Bayesian Regularization and Proximal Operators

abstract: In this work, we broadly connect kernel-based filtering (e.g. approaches such as the bilateral filters and nonlocal means, but also many more) with general variational formulations of Bayesian regularized least squares, and the related concept of proximal operators. The latter set of variational/Bayesian/proximal formulations often result in optimization problems that do not have closed-form solutions, and therefore typically require global iterative solutions. Our main contribution here is to establish how one can approximate the solution of the resulting global optimization problems with use of locally adaptive filters with specific kernels. Our results are valid for small regularization strength but the approach is powerful enough to be useful for a wide range of applications because we expose how to derive a "kernelized" solution to these problems that approximates the global solution in one-shot, using only local operations. As another side benefit in the reverse direction, given a local data-adaptive filter constructed with a particular choice of kernel, we enable the interpretation of such filters in the variational/Bayesian/proximal framework.


title: Response to the authors of 'On the (un)effectiveness of Proton Boron Capture in Proton Therapy'

abstract: This manuscript provides a response to a recent report by Mazzone et al. available online on arXiv that, in turn, tentatively aims at demonstrating the inefficacy of proton boron capture in hadrotherapy. We clarify that Mazzone et al. do not add any scientific or technical insights to the points extensively discussed in the original manuscript by Cirrone et al., and/or in the series of iterations had with the Referee, which ultimately lead to the publication of our original and pioneering experimental work. Here we summarize some of the key points of the long scientific debate we had during the review process of paper by Cirrone et al., which are very similar to the considerations presented by Mazzone et al.. In conclusion, no quantitative explanation of our robust experimental achievements presented in Cirrone et al. is provided in Mazzone et al.


title: Current fluctuations in quantum absorption refrigerators

abstract: Absorption refrigerators transfer thermal energy from a cold bath to a hot bath without input power by utilizing heat from an additional "work" reservoir. Particularly interesting is a three-level design for a quantum absorption refrigerator, which can be optimized to reach the maximal (Carnot) cooling efficiency. Previous studies of three-level chillers focused on the behavior of the averaged cooling current. Here, we go beyond that and study the full counting statistics of heat exchange in a three-level chiller model. We explain how to obtain the complete cumulant generating function of the refrigerator in steady state, then derive a partial cumulant generating function, which yields closed-form expressions for both the averaged cooling current and its noise. Our analytical results and simulations are beneficial for the design of nanoscale engines and cooling systems far from equilibrium, with their performance optimized according to different criteria, efficiency, power, fluctuations and dissipation.


title: Accelerated Wirtinger Flow for Multiplexed Fourier Ptychographic Microscopy

abstract: Fourier ptychographic microscopy enables gigapixel-scale imaging, with both large field-of-view and high resolution. Using a set of low-resolution images that are recorded under varying illumination angles, the goal is to computationally reconstruct high-resolution phase and amplitude images. To increase temporal resolution, one may use multiplexed measurements where the sample is illuminated simultaneously from a subset of the angles. In this paper, we develop an algorithm for Fourier ptychographic microscopy with such multiplexed illumination. Specifically, we consider gradient descent type updates and propose an analytical step size that ensures the convergence of the iterates to a stationary point. Furthermore, we propose an accelerated version of our algorithm (with the same step size) which significantly improves the convergence speed. We demonstrate that the practical performance of our algorithm is identical to the case where the step size is manually tuned. Finally, we apply our parameter-free approach to real data and validate its applicability.


title: The dynamics of the angular and radial density correlation scaling exponents in fractal to non-fractal morphodynamics

abstract: Fractal/non-fractal morphological transitions allow for the systematic study of the physics behind fractal morphogenesis in nature. In these systems, the fractal dimension is considered a non-thermal order parameter, commonly and equivalently computed from the scaling of the two-point radial- or angular-density correlations. However, these two quantities lead to discrepancies during the analysis of basic systems, such as in the diffusion-limited aggregation fractal. Hence, the corresponding clarification regarding the limits of the radial/angular scaling equivalence is needed. In this work, considering three fundamental fractal/non-fractal transitions in two dimensions, we show that the unavoidable emergence of growth anisotropies is responsible for the breaking-down of the radial/angular equivalence. Specifically, we show that the angular scaling behaves as a critical power-law, whereas the radial scaling as an exponential that, under the fractal dimension interpretation, resemble first- and second-order transitions, respectively. Remarkably, these and previous results can be unified under a single fractal dimensionality equation.


title: TRAJEDI: Trajectory Dissimilarity

abstract: The vast increase in our ability to obtain and store trajectory data necessitates trajectory analytics techniques to extract useful information from this data. Pair-wise distance functions are a foundation building block for common operations on trajectory datasets including constrained SELECT queries, k-nearest neighbors, and similarity and diversity algorithms. The accuracy and performance of these operations depend heavily on the speed and accuracy of the underlying trajectory distance function, which is in turn affected by trajectory calibration. Current methods either require calibrated data, or perform calibration of the entire relevant dataset first, which is expensive and time consuming for large datasets. We present TRAJEDI, a calibrationaware pair-wise distance calculation scheme that outperforms naive approaches while preserving accuracy. We also provide analyses of parameter tuning to trade-off between speed and accuracy. Our scheme is usable with any diversity, similarity or k-nearest neighbor algorithm.


title: Low-Rank Solution Methods for Stochastic Eigenvalue Problems

abstract: We study efficient solution methods for stochastic eigenvalue problems arising from discretization of self-adjoint partial differential equations with random data. With the stochastic Galerkin approach, the solutions are represented as generalized polynomial chaos expansions. A low-rank variant of the inverse subspace iteration algorithm is presented for computing one or several minimal eigenvalues and corresponding eigenvectors of parameter-dependent matrices. In the algorithm, the iterates are approximated by low-rank matrices, which leads to significant cost savings. The algorithm is tested on two benchmark problems, a stochastic diffusion problem with some poorly separated eigenvalues, and an operator derived from a discrete stochastic Stokes problem whose minimal eigenvalue is related to the inf-sup stability constant. Numerical experiments show that the low-rank algorithm produces accurate solutions compared to the Monte Carlo method, and it uses much less computational time than the original algorithm without low-rank approximation.


title: Simultaneous Broadband Vector Magnetometry Using Solid-State Spins

abstract: We demonstrate a vector magnetometer that simultaneously measures all Cartesian components of a dynamic magnetic field using an ensemble of nitrogen-vacancy (NV) centers in a single-crystal diamond. Optical NV-diamond measurements provide high-sensitivity, broadband magnetometry under ambient or extreme physical conditions; and the fixed crystallographic axes inherent to this solid-state system enable vector sensing free from heading errors. In the present device, multi-channel lock-in detection extracts the magnetic-field-dependent spin resonance shifts of NVs oriented along all four tetrahedral diamond axes from the optical signal measured on a single detector. The sensor operates from near DC up to a $12.5$ kHz measurement bandwidth; and simultaneously achieves $\sim\!50$ pT/$\sqrt{\text{Hz}}$ magnetic field sensitivity for each Cartesian component, which is to date the highest demonstrated sensitivity of a full vector magnetometer employing solid-state spins. Compared to optimized devices interrogating the four NV orientations sequentially, the simultaneous vector magnetometer enables a $4\times$ measurement speedup. This technique can be extended to pulsed-type sensing protocols and parallel wide-field magnetic imaging.

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


title: DeepMoTIon: Learning to Navigate Like Humans

abstract: We present a novel human-aware navigation approach, where the robot learns to mimic humans to navigate safely in crowds. The presented model, referred to as DeepMoTIon, is trained with pedestrian surveillance data to predict human velocity in the environment. The robot processes LiDAR scans via the trained network to navigate to the target location. We conduct extensive experiments to assess the components of our network and prove their necessity to imitate humans. Our experiments show that DeepMoTIion outperforms all the benchmarks in terms of human imitation, achieving a 24% reduction in time series-based path deviation over the next best approach. In addition, while many other approaches often failed to reach the target, our method reached the target in 100% of the test cases while complying with social norms and ensuring human safety.


title: Reduced modeling of porous media convection in a minimal flow unit at large Rayleigh number

abstract: Direct numerical simulations (DNS) indicate that at large values of the Rayleigh number ($Ra$) convection in porous media self-organizes into narrowly-spaced columnar flows, with more complex spatiotemporal features being confined to boundary layers near the top and bottom walls. In this investigation of high-$Ra$ porous media convection in a minimal flow unit, two reduced modeling strategies are proposed that exploit these specific flow characteristics. Both approaches utilize the idea of decomposition since the flow exhibits different dynamics in different regions of the domain: small-scale cellular motions generally are localized within the thermal and vorticity boundary layers near the upper and lower walls, while in the interior, the flow exhibits persistent large-scale structures and only a few low (horizontal) wavenumber Fourier modes are active. Accordingly, in the first strategy, the domain is decomposed into two near-wall regions and one interior region. Our results confirm that suppressing the interior high-wavenumber modes has negligible impact on the essential structural features and transport properties of the flow. In the second strategy, a hybrid reduced model is constructed by using Galerkin projection onto a fully \emph{a priori} eigenbasis drawn from energy stability and upper bound theory, thereby extending the model reduction strategy developed by Chini \emph{et al.} (\emph{Physica~D}, vol. 240, 2011, pp. 241--248) to large $Ra$. The results indicate that the near-wall upper-bound eigenmodes can economically represent the small-scale rolls within the exquisitely-thin thermal boundary layers. Relative to DNS, the hybrid algorithm enables over an order-of-magnitude increase in computational efficiency with only a modest loss of accuracy.


title: Multiscale Methods for Model Order Reduction of Non Linear Multiphase Flow Problems

abstract: Numerical simulations for flow and transport in subsurface porous media often prove computationally prohibitive due to property data availability at multiple spatial scales that can vary by orders of magnitude. A number of model order reduction approaches are available in the existing literature that alleviate this issue by approximating the solution at a coarse scale. We attempt to present a comparison between two such model order reduction techniques, namely: (1) adaptive numerical homogenization and (2) generalized multiscale basis functions. We rely upon a non-linear, multi-phase, black-oil model formulation, commonly encountered in the oil and gas industry, as the basis for comparing the aforementioned two approaches. An expanded mixed finite element formulation is used to separate the spatial scales between non-linear, flow and transport problems. To the author's knowledge this is the first time these approaches have been described for a practical non-linear, multiphase flow problem of interest. A numerical benchmark is setup using fine scale property information from the 10$^{th}$ SPE comparative project dataset for the purpose of comparing accuracies of these two schemes. An adaptive criterion is employed in by both the schemes for local enrichment that allows us to preserve solution accuracy compared to the fine scale benchmark problem. The numerical results indicate that both schemes are able to adequately capture the fine scale features of the model problem at hand.


title: Random Partitions and Cohen-Lenstra Heuristics

abstract: We investigate combinatorial properties of a family of probability distributions on finite abelian p-groups. This family includes several well-known distributions as specializations. These specializations have been studied in the context of Cohen-Lenstra heuristics and cokernels of families of random p-adic matrices.


title: Design, construction, and characterization of a compact DD neutron generator designed for 40Ar/39Ar geochronology

abstract: A next-generation, high-flux DD neutron generator has been designed, commissioned, and characterized, and is now operational in a new facility at the University of California Berkeley. The generator, originally designed for 40Ar/39Ar dating of geological materials, has since served numerous additional applications, including medical isotope production studies, with others planned for the near future. In this work, we present an overview of the High Flux Neutron Generator (HFNG) which includes a variety of simulations, analytical models, and experimental validation of results. Extensive analysis was performed in order to characterize the neutron yield, flux, and energy distribution at specific locations where samples may be loaded for irradiation. A notable design feature of the HFNG is the possibility for sample irradiation internal to the cathode, just 8 mm away from the neutron production site, thus maximizing the neutron flux (n/cm2/s). The generator's maximum neutron flux at this irradiation position is 2.58e7 n/cm2/s +/- 5% (approximately 3e8 n/s total yield) as measured via activation of small natural indium foils. However, future development is aimed at achieving an order of magnitude increase in flux. Additionally, the deuterium ion beam optics were optimized by simulations for various extraction configurations in order to achieve a uniform neutron flux distribution and an acceptable heat load. Finally, experiments were performed in order to benchmark the modeling and characterization of the HFNG.


title: Contour Parametrization via Anisotropic Mean Curvature Flows

abstract: We present a new implementation of anisotropic mean curvature flow for contour recognition. Our procedure couples the mean curvature flow of planar closed smooth curves, with an external field from a potential of point-wise charges. This coupling constrains the motion when the curve matches a picture placed as background. We include a stability criteria for our numerical approximation.


title: Dynamically Efficient Kinematics for Hyper-Redundant Manipulators

abstract: A hyper-redundant robotic arm is a manipulator with many degrees of freedom, capable of executing tasks in cluttered environments where robotic arms with fewer degrees of freedom are unable to operate. This paper introduces a new method for modeling those manipulators in a completely dynamic way. The proposed method enables online changes of the kinematic structure with the use of a special function; termed "meta-controlling function". This function can be used to develop policies to reduce drastically the computational cost for a single task, and to robustly control the robotic arm, even in the event of partial damage. The direct and inverse kinematics are solved for a generic three-dimensional articulated hyper-redundant arm, that can be used as a proof of concept for more specific structures. To demonstrate the robustness of our method, experimental simulation results, for a basic "meta-controlling" function, are presented.


title: A new route to finding bounds on the generalized spectrum of many physical operators

abstract: Here we obtain bounds on the spectrum of that operator whose inverse, when it exists, gives the Green's function. We consider the wide of physical problems that can be cast in a form where a constitutive equation ${\bf J}({\bf x})={\bf L}({\bf x}){\bf E}({\bf x})-{\bf h}({\bf x})$ with a source term ${\bf h}({\bf x})$ holds for all ${\bf x}$ in some domain $\Omega$, and relates fields ${\bf E}$ and ${\bf J}$ that satisfy appropriate differential constraints, symbolized by ${\bf E}\in\cal{E}_\Omega$ and ${\bf J}\in\cal{J}_\Omega$ where $\cal{E}_\Omega$ and $\cal{J}_\Omega$ are orthogonal spaces that span the space $\cal{H}_\Omega$ of square-integrable fields in which ${\bf h}$ lies. Boundedness and coercivity conditions on the moduli ${\bf L}({\bf x})$ ensure there exists a unique ${\bf E}$ for any given ${\bf h}$, i.e., ${\bf E}={\bf G}_\Omega{\bf h}$ which then establishes the existence of the Green's function ${\bf G}_\Omega$. We show that the coercivity condition is guaranteed to hold if weaker conditions, involving generalized quasiconvex functions, are satisfied. The advantage is that these weaker conditions are easier to verify, and for multiphase materials they can be independent of the geometry of the phases. For ${\bf L}({\bf x} )$ depending linearly on a vector of parameters ${\bf z}=(z_1, z_2,\ldots, z_n)$, we obtain constraints on ${\bf z}$ that ensure the Green's function exists, and hence which provide bounds on the spectrum.


title: Thermal Management in Fine-Grained 3-D Integrated Circuits

abstract: For beyond 2-D CMOS logic, various 3-D integration approaches specially transistor based 3-D integrations such as monolithic 3-D [1], Skybridge [2], SN3D [3] holds most promise. However, such 3D architectures within small form factor increase hotspots and demand careful consideration of thermal management at all levels of integration [4] as stacked transistors are detached from the substrate (i.e., heat sink). Traditional system level approaches such as liquid cooling [5], heat spreader [6], etc. are inadequate for transistor level 3-D integration and have huge cost overhead [7]. In this paper, we investigate the thermal profile for transistor level 3-D integration approaches through finite element based modeling. Additionally, we propose generic physical level heat management features for such transistor level 3-D integration and show their application through detailed thermal modeling and simulations. These features include a thermal junction and heat conducting nano pillar. The heat junction is a specialized junction to extract heat from a selected region in 3-D; it allows heat conduction without interference with the electrical activities of the circuit. In conjunction with the junction, our proposed thermal pillars enable heat dissipation through the substrate; these pillars are analogous to TSVs/Vias, but carry only heat. Such structures are generic and is applicable to any transistor level 3-D integration approaches. We perform 3-D finite element based analysis to capture both static and transient thermal behaviors of 3-D circuits, and show the effectiveness of heat management features. Our simulation results show that without any heat extraction feature, temperature for 3-D integrated circuits increased by almost 100K-200K. However, proposed heat extraction feature is very effective in heat management, reducing temperature from heated area by up to 53%.


title: Geodesic nets with three boundary vertices

abstract: We prove that a geodesic net with three boundary (= unbalanced) vertices on a non-positively curved plane has at most one balanced vertex. We do not assume any a priori bound for the degrees of unbalanced vertices. The result seems to be new even in the Euclidean case. We demonstrate by examples that the result is not true for metrics of positive curvature on the plane, and that there are no immediate generalizations of this result for geodesic nets with four unbalanced vertices.


title: A Large-Scale Multi-Institutional Evaluation of Advanced Discrimination Algorithms for Buried Threat Detection in Ground Penetrating Radar

abstract: In this paper we consider the development of algorithms for the automatic detection of buried threats using ground penetrating radar (GPR) measurements. GPR is one of the most studied and successful modalities for automatic buried threat detection (BTD), and a large variety of BTD algorithms have been proposed for it. Despite this, large-scale comparisons of GPR-based BTD algorithms are rare in the literature. In this work we report the results of a multi-institutional effort to develop advanced buried threat detection algorithms for a real-world GPR BTD system. The effort involved five institutions with substantial experience with the development of GPR-based BTD algorithms. In this paper we report the technical details of the advanced algorithms submitted by each institution, representing their latest technical advances, and many state-of-the-art GPR-based BTD algorithms. We also report the results of evaluating the algorithms from each institution on the large experimental dataset used for development. The experimental dataset comprised 120,000 m^2 of GPR data using surface area, from 13 different lanes across two US test sites. The data was collected using a vehicle-mounted GPR system, the variants of which have supplied data for numerous publications. Using these results, we identify the most successful and common processing strategies among the submitted algorithms, and make recommendations for GPR-based BTD algorithm design.


title: A White Paper Submitted to The National Academy of Science's Committee on Exoplanet Science Strategy: Observing Exoplanets with the James Webb Space Telescope

abstract: The James Webb Space Telescope (JWST) will revolutionize our understanding of exoplanets with transit spectroscopy of a wide range of mature planets close to their host stars ($<$2 AU) and with coronagraphic imaging and spectroscopy of young objects located further out ($>$10 AU). The census of exoplanets has revealed an enormous variety of planets orbiting stars of all ages and spectral types. With TESS adding to this census with its all-sky survey of the closest, brightest stars, the challenge of the coming decade will be to move from demography to physical characterization. This white paper discusses the wide variety of exoplanet opportunities enabled by JWST's sensitivity and stability, its high angular resolution, and its suite of powerful instruments. JWST observations will advance our understanding of the atmospheres of young to mature planets and will provide new insights into planet formation.


title: Model Structural Inference using Local Dynamic Operators

abstract: This paper focuses on the problem of quantifying the effects of model-structure uncertainty in the context of time-evolving dynamical systems. This is motivated by multi-model uncertainty in computer physics simulations: developers often make different modeling choices in numerical approximations and process simplifications, leading to different numerical codes that ostensibly represent the same underlying dynamics. We consider model-structure inference as a two-step methodology: the first step is to perform system identification on numerical codes for which it is possible to observe the full state; the second step is structural uncertainty quantification (UQ), in which the goal is to search candidate models "close" to the numerical code surrogates for those that best match a quantity-of-interest (QOI) from some empirical dataset. Specifically, we: (1) define a discrete, local representation of the structure of a partial differential equation, which we refer to as the "local dynamical operator" (LDO); (2) identify model structure non-intrusively from numerical code output; (3) non-intrusively construct a reduced order model (ROM) of the numerical model through POD-DEIM-Galerkin projection; (4) perturb the ROM dynamics to approximate the behavior of alternate model structures; and (5) apply Bayesian inference and energy conservation laws to calibrate a LDO to a given QOI. We demonstrate these techniques using the two-dimensional rotating shallow water (RSW) equations as an example system.


title: Precision Space Astrometry as a Tool to Find Earth-like Exoplanets

abstract: Because of the recent technological advances, the key technologies needed for precision space optical astrometry are now in hand. The Microarcsecond Astrometry Probe (MAP) mission concept is designed to find 1 Earth mass planets at 1AU orbit (scaled to solar luminosity) around the nearest ~90 FGK stars. The MAP payload includes i) a single three-mirror anastigmatic telescope with a 1-m primary mirror and metrology subsystems, and ii) a camera. The camera focal plane consists of 42 detectors, providing a Nyquist sampled FOV of 0.4-deg. Its metrology subsystems ensure that MAP can achieve the 0.8 uas astrometric precision in 1 hr, which is required to detect Earth-like exoplanets in our stellar neighborhood. MAP mission could provide ~10 specific targets for a much larger coronagraphic mission that would measure its spectra. We argue for the development of the space astrometric missions capable of finding Earth-2.0. Given the current technology readiness such missions relying on precision astrometry could be flown in the next decade, perhaps in collaboration with other national space agencies.


title: Mobile Edge Computing for Cellular-Connected UAV: Computation Offloading and Trajectory Optimization

abstract: This paper studies a new mobile edge computing (MEC) setup where an unmanned aerial vehicle (UAV) is served by cellular ground base stations (GBSs) for computation offloading. The UAV flies between a give pair of initial and final locations, during which it needs to accomplish certain computation tasks by offloading them to some selected GBSs along its trajectory for parallel execution. Under this setup, we aim to minimize the UAV's mission completion time by optimizing its trajectory jointly with the computation offloading scheduling, subject to the maximum speed constraint of the UAV, and the computation capacity constraints at GBSs. The joint UAV trajectory and computation offloading optimization problem is, however, non-convex and thus difficult to be solved optimally. To tackle this problem, we propose an efficient algorithm to obtain a high-quality suboptimal solution. Numerical results show that the proposed design significantly reduces the UAV's mission completion time, as compared to benchmark schemes.


title: Quantum heat engine based on trapped Bose gases: Its maximum efficiency can approach the Carnot value at finite power

abstract: It was reported that, if and only if the specific heat, correlation length, and dynamical exponents $\alpha, \nu$ and $z$, fulfill the condition $\alpha-z\nu>0$, the phase transitions can enable a quantum heat engine to approach Carnot efficiency at finite power. We start our analysis via a different approach in which the effects of interaction and fluctuations on the Hamiltonian of a trapped dilute Bose gas belonging to the same universality as $XY$ model. Based on models of quantum Otto heat engines, we find the general expression of the efficiency which includes the correction due to interaction and fluctuations at the critical point, and show that, near the Bose-Einstein-condensation point with $\alpha-z\nu<0$, energy fluctuations could enable attaintment of the Carnot efficiency with nonvanishing power. Such quantum heat engines can also be realized by changing the shape of the trap confining the ideal and weakly interacting Bose gas during the adiabatic processes of the cycle. These quantum heat engines working with the trapped Bose gases, which are based on techniques of cooling Bose condensates and could be realizable at present technology.


title: Attention-based Graph Neural Network for Semi-supervised Learning

abstract: Recently popularized graph neural networks achieve the state-of-the-art accuracy on a number of standard benchmark datasets for graph-based semi-supervised learning, improving significantly over existing approaches. These architectures alternate between a propagation layer that aggregates the hidden states of the local neighborhood and a fully-connected layer. Perhaps surprisingly, we show that a linear model, that removes all the intermediate fully-connected layers, is still able to achieve a performance comparable to the state-of-the-art models. This significantly reduces the number of parameters, which is critical for semi-supervised learning where number of labeled examples are small. This in turn allows a room for designing more innovative propagation layers. Based on this insight, we propose a novel graph neural network that removes all the intermediate fully-connected layers, and replaces the propagation layers with attention mechanisms that respect the structure of the graph. The attention mechanism allows us to learn a dynamic and adaptive local summary of the neighborhood to achieve more accurate predictions. In a number of experiments on benchmark citation networks datasets, we demonstrate that our approach outperforms competing methods. By examining the attention weights among neighbors, we show that our model provides some interesting insights on how neighbors influence each other.


title: Joint Optimization of Scheduling and Routing in Multicast Wireless Ad-Hoc Network Using Soft Graph Coloring and Non-linear Cubic Games

abstract: In this paper we present matrix game-theoretic models for joint routing, network coding, and scheduling problem. First routing and network coding are modeled by using a new approach based on compressed topology matrix that takes into account the inherent multicast gain of the network. The scheduling is optimized by a new approach called network graph soft coloring. Soft graph coloring is designed by switching between different components of a wireless network graph, which we refer to as graph fractals, with appropriate usage rates. The network components, represented by graph fractals, are a new paradigm in network graph partitioning that enables modeling of the network optimization problem by using the matrix game framework. In the proposed game which is a nonlinear cubic game, the strategy sets of the players are links, path, and network components. The outputs of this game model are mixed strategy vectors of the second and the third players at equilibrium. Strategy vector of the second player specifies optimum multi-path routing and network coding solution while mixed strategy vector of the third players indicates optimum switching rate among different network components or membership probabilities for optimal soft scheduling approach. Optimum throughput is the value of the proposed nonlinear cubic game at equilibrium. The proposed nonlinear cubic game is solved by extending fictitious playing method. Numerical and simulation results prove the superior performance of the proposed techniques compared to the conventional schemes using hard graph coloring.


title: Enhancing Evolutionary Conversion Rate Optimization via Multi-armed Bandit Algorithms

abstract: Conversion rate optimization means designing web interfaces such that more visitors perform a desired action (such as register or purchase) on the site. One promising approach, implemented in Sentient Ascend, is to optimize the design using evolutionary algorithms, evaluating each candidate design online with actual visitors. Because such evaluations are costly and noisy, several challenges emerge: How can available visitor traffic be used most efficiently? How can good solutions be identified most reliably? How can a high conversion rate be maintained during optimization? This paper proposes a new technique to address these issues. Traffic is allocated to candidate solutions using a multi-armed bandit algorithm, using more traffic on those evaluations that are most useful. In a best-arm identification mode, the best candidate can be identified reliably at the end of evolution, and in a campaign mode, the overall conversion rate can be optimized throughout the entire evolution process. Multi-armed bandit algorithms thus improve performance and reliability of machine discovery in noisy real-world environments.


title: Stochastic Models of Coalition Games for Spectrum Sharing in Large Scale Interference Channels

abstract: In this paper, we present a framework for the analysis of self-organized distributed coalition formation process for spectrum sharing in interference channel for large-scale ad hoc networks. In this approach, we use the concept of coalition clusters within the network where mutual interdependency between different clusters is characterized by the concept of spatial network correlation. Then by using stochastic models of the process we give up some details characteristic for coalition game theory in order to be able to include some additional parameters for network scaling. Applications of this model are a) Estimation of average time to reach grand coalition and its variance through closed-form equations. These parameters are important in designing the process in a dynamic environment. b) Dimensioning the coalition cluster within the network c) Modelling the network spatial correlation characterizing mutual visibility of the interfering links. d) Modeling of the effect of the new link activation/inactivation on the coalition forming process. e) Modeling the effect of link mobility on the coalition-forming process.


title: Disentangled cooperative orderings in artificial rare-earth nickelates

abstract: Coupled transitions between distinct ordered phases are important aspects behind the rich phase complexity of correlated oxides that hinders our understanding of the underlying phenomena. For this reason, fundamental control over complex transitions has become a leading motivation of the designer approach to materials. We have devised a series of new superlattices by combining a Mott insulator and a correlated metal to form ultra-short period superlattices, which allow one to disentangle the simultaneous orderings in $RE$NiO$_3$. Tailoring an incommensurate heterostructure period relative to the bulk charge ordering pattern suppresses the charge order transition while preserving metal-insulator and antiferromagnetic transitions. Such selective decoupling of the entangled phases resolves the long-standing puzzle about the driving force behind the metal-insulator transition and points to the site selective Mott transition as the operative mechanism. This designer approach emphasizes the potential of heterointerfaces for selective control of simultaneous transitions in complex materials with entwined broken symmetries.


title: Cluster Size Optimization in Cooperative Spectrum Sensing

abstract: In this paper, we study and optimize the cooperation cluster size in cooperative spectrum sensing to maximize the throughput of secondary users (SUs). To calculate the effective throughput, we assume each SU spends just 1 symbol to negotiate with the other SUs in its transmission range. This is the minimum overhead required for each SU to broadcast its sensing decision to the other members of the cluster. When the number of SUs is large, the throughput spent for the negotiation is noticeable and therefore increasing the cooperation cluster size does not improve the effective throughput anymore. In this paper, we calculate the effective throughput as a function of the cooperation cluster size, and then we maximize the throughput by finding the optimal cluster size. Various numerical results show that when decisions are combined by the OR-rule, the optimum cooperation cluster size is less than when the AND-rule is used. On the other hand, the optimum cluster size monotonically decreases with the increase in the average SNR of the SUs. Another interesting result is that when the cluster size is optimized the OR-rule always outperforms the AND-rule.


title: Tokunaga self-similarity arises naturally from time invariance

abstract: The Tokunaga condition is an algebraic rule that provides a detailed description of the branching structure in a self-similar tree. Despite a solid empirical validation and practical convenience, the Tokunaga condition lacks a theoretical justification. Such a justification is suggested in this work. We define a geometric branching processes $\mathcal{G}(s)$ that generates self-similar rooted trees. The main result establishes the equivalence between the invariance of $\mathcal{G}(s)$ with respect to a time shift and a one-parametric version of the Tokunaga condition. In the parameter region where the process satisfies the Tokunaga condition (and hence is time invariant), $\mathcal{G}(s)$ enjoys many of the symmetries observed in a critical binary Galton-Watson branching process and reproduce the latter for a particular parameter value.


title: K-shell decomposition reveals hierarchical cortical organization of the human brain

abstract: In recent years numerous attempts to understand the human brain were undertaken from a network point of view. A network framework takes into account the relationships between the different parts of the system and enables to examine how global and complex functions might emerge from network topology. Previous work revealed that the human brain features 'small world' characteristics and that cortical hubs tend to interconnect among themselves. However, in order to fully understand the topological structure of hubs one needs to go beyond the properties of a specific hub and examine the various structural layers of the network. To address this topic further, we applied an analysis known in statistical physics and network theory as k-shell decomposition analysis. The analysis was applied on a human cortical network, derived from MRI\DSI data of six participants. Such analysis enables us to portray a detailed account of cortical connectivity focusing on different neighborhoods of interconnected layers across the cortex. Our findings reveal that the human cortex is highly connected and efficient, and unlike the internet network contains no isolated nodes. The cortical network is comprised of a nucleus alongside shells of increasing connectivity that formed one connected giant component. All these components were further categorized into three hierarchies in accordance with their connectivity profile, with each hierarchy reflecting different functional roles. Such a model may explain an efficient flow of information from the lowest hierarchy to the highest one, with each step enabling increased data integration. At the top, the highest hierarchy (the nucleus) serves as a global interconnected collective and demonstrates high correlation with consciousness related regions, suggesting that the nucleus might serve as a platform for consciousness to emerge.


title: Classical and quantum dissipative dynamics in Josephson junctions: an Arnold problem, bifurcation and capture into resonance

abstract: We theoretically study the phase dynamics in Josephson junctions, which maps onto the oscillatory motion of a point-like particle in the washboard potential. Under appropriate driving and damping conditions, the Josephson phase undergoes intriguing bistable dynamics near a saddle point in the quasienergy landscape. The bifurcation mechanism plays a critical role in superconducting quantum circuits with relevance to non-demolition measurements such as high-fidelity readout of qubit states. We address the question `what is the probability of capture into either basin of attraction' and answer it concerning both classical and quantum dynamics. Consequently, we derive the Arnold probability and numerically analyze its implementation of the controlled dynamical switching between two steady states under the various nonequilibrium conditions.


title: Enhanced Optimization with Composite Objectives and Novelty Selection

abstract: An important benefit of multi-objective search is that it maintains a diverse population of candidates, which helps in deceptive problems in particular. Not all diversity is useful, however: candidates that optimize only one objective while ignoring others are rarely helpful. This paper proposes a solution: The original objectives are replaced by their linear combinations, thus focusing the search on the most useful tradeoffs between objectives. To compensate for the loss of diversity, this transformation is accompanied by a selection mechanism that favors novelty. In the highly deceptive problem of discovering minimal sorting networks, this approach finds better solutions, and finds them faster and more consistently than standard methods. It is therefore a promising approach to solving deceptive problems through multi-objective optimization.


title: Evolutionary Architecture Search For Deep Multitask Networks

abstract: Multitask learning, i.e. learning several tasks at once with the same neural network, can improve performance in each of the tasks. Designing deep neural network architectures for multitask learning is a challenge: There are many ways to tie the tasks together, and the design choices matter. The size and complexity of this problem exceeds human design ability, making it a compelling domain for evolutionary optimization. Using the existing state of the art soft ordering architecture as the starting point, methods for evolving the modules of this architecture and for evolving the overall topology or routing between modules are evaluated in this paper. A synergetic approach of evolving custom routings with evolved, shared modules for each task is found to be very powerful, significantly improving the state of the art in the Omniglot multitask, multialphabet character recognition domain. This result demonstrates how evolution can be instrumental in advancing deep neural network and complex system design in general.


title: Signal Optimization with HV divider of MCP-PMT for JUNO

abstract: The Jiangmen Underground Neutrino Observatory (JUNO) is proposed to determine the neutrino mass hierarchy using a 20 kiloton underground liquid scintillator detector (CD). One of the keys is the energy resolution of the CD to reach <3% at 1 MeV, where totally 15,000 MCP-PMT will be used. The optimization of the 20-inch MCP-PMT is very important for better detection efficiency and stable performance. In this work, we will show the study to optimize the MCP-PMT working configuration for charge measurement. Particularly, the quality of PMT signal is another key for high-precision neutrino experiments while most of these experiments are affected by the overshoot of PMT signal from the positive HV scheme. The overshoot coupled with positive HV which is troubling trigger, dead time and precise charge measurement, we have studied to control it to less than 1% of signal amplitude for a better physics measurement. In this article, on the one hand, the optimized HV divider ratio will be presented here to improve its collection efficiency; on the other hand, we will introduce the method to reduce the ratio of overshoot from 10% to 1%.


title: Emergence of spiral dark solitons in the merging of rotating Bose-Einstein condensates

abstract: Merging of isolated Bose-Einstein condensates (BECs) is an important topic due to its relevance to matter-wave interferometry and the Kibble-Zurek mechanism. Many past research focused on merging of BECs with uniform initial phases. In our recent numerical study (Phys. Rev. A 97, 013612 (2018)), we revealed that upon merging of rotating BECs with non-uniform initial phases, spiral-shaped dark solitons can emerge. These solitons facilitate angular momentum transfer and allow the merged condensate to rotate even in the absence of quantized vortices. More strikingly, the sharp endpoints of these spiral solitons can induce rotational motion in the BECs like vortices but with effectively a fraction of a quantized circulation. This paper reports our systematic study on the merging dynamics of rotating BECs. We discuss how the potential barrier that initially separates the BECs can affect the profile of the spiral solitons. We also show that the number of spiral solitons created in the BECs matches the relative winding number of the rotating BECs. The underlying mechanism of the observed soliton dynamics is explained.


title: Integrated Optimization of Partitioning, Scheduling and Floorplanning for Partially Dynamically Reconfigurable Systems

abstract: Confronted with the challenge of high performance for applications and the restriction of hardware resources for field-programmable gate arrays (FPGAs), partial dynamic reconfiguration (PDR) technology is anticipated to accelerate the reconfiguration process and alleviate the device shortage. In this paper, we propose an integrated optimization framework for task partitioning, scheduling and floorplanning on partially dynamically reconfigurable FPGAs. The partitions, schedule, and floorplan of the tasks are represented by the partitioned sequence triple P-ST (PS,QS,RS), where (PS,QS) is a hybrid nested sequence pair (HNSP) for representing the spatial and temporal partitions, as well as the floorplan, and RS is the partitioned dynamic configuration order of the tasks. The floorplanning and scheduling of task modules can be computed from the partitioned sequence triple P-ST in O(n^2) time. To integrate the exploration of the scheduling and floorplanning design space, we use a simulated annealing-based search engine and elaborate a perturbation method, where a randomly chosen task module is removed from the partition sequence triple and then inserted back into a proper position selected from all the (n+1)^3 possible combinations of partitions, schedule and floorplan. The experimental results demonstrate the efficiency and effectiveness of the proposed framework.


title: Counting trees in a graph

abstract: We discuss a recursive formula for number of spanning trees in a graph. The paper is written primary for school students.


title: Charge Transfer Effects in Naturally Occurring van der Waals Heterostructures (PbSe)1.16(TiSe2)m (m=1, 2)

abstract: Van der Waals heterostructures (VDWHs) exhibit rich properties and thus has potential for applications, and charge transfer between different layers in a heterostructure often dominates its properties and device performance. It is thus critical to reveal and understand the charge transfer effects in VDWHs, for which electronic structure measurements have proven to be effective. Using angle-resolved photoemission spectroscopy, we studied the electronic structures of (PbSe)1.16(TiSe2)m(m=1, 2), which are naturally occurring VDWHs, and discovered several striking charge transfer effects. When the thickness of the TiSe2 layers is halved from m=2 to m=1, the amount of charge transferred increases unexpectedly by more than 250%. This is accompanied by a dramatic drop in the electron-phonon interaction strength far beyond the prediction by first-principles calculations and, consequently, superconductivity only exists in the m=2 compound with strong electron-phonon interaction, albeit with lower carrier density. Furthermore, we found that the amount of charge transferred in both compounds is nearly halved when warmed from below 10 K to room temperature, due to the different thermal expansion coefficients of the constituent layers of these misfit compounds. These unprecedentedly large charge transfer effects might widely exist in VDWHs composed of metal-semiconductor contacts; thus, our results provide important insights for further understanding and applications of VDWHs.


title: Exploring Extreme Space Weather Factors of Exoplanetary Habitability

abstract: It is currently unknown how common life is on exoplanets, or how long planets can remain viable for life. To date, we have a superficial notion of habitability, a necessary first step, but so far lacking an understanding of the detailed interaction between stars and planets over geological timescales, dynamical evolution of planetary systems, and atmospheric evolution on planets in other systems. A planet mass, net insolation, and atmospheric composition alone are insufficient to determine the probability that life on a planet could arise or be detected. The latter set of planetary considerations, among others, underpin the concept of the habitable zone (HZ), defined as the circumstellar region where standing bodies of liquid water could be supported on the surface of a rocky planet. However, stars within the same spectral class are often treated in the same way in HZ studies, without any regard for variations in activity among individual stars. Such formulations ignore differences in how nonthermal emission and magnetic energy of transient events in different stars affect the ability of an exoplanet to retain its atmosphere.In the last few years there has been a growing appreciation that the atmospheric chemistry, and even retention of an atmosphere in many cases, depends critically on the high-energy radiation and particle environments around these stars. Indeed, recent studies have shown stellar activity and the extreme space weather, such as that created by the frequent flares and coronal mass ejections (CMEs) from the active stars and young Sun, may have profoundly affected the chemistry and climate and thus habitability of the early Earth and terrestrial type exoplanets. The goal of this white paper is to identify and describe promising key research goals to aid the field of the exoplanetary habitability for the next 20 years.


title: Optimum Linear Codes with Support Constraints over Small Fields

abstract: We consider the problem of designing optimal linear codes (in terms of having the largest minimum distance) subject to a support constraint on the generator matrix. We show that the largest minimum distance can be achieved by a subcode of a Reed-Solomon code of small field size. As a by-product of this result, we settle the GM-MDS conjecture of Dau et. al. in the affirmative.


title: On a metric generalization of the $tt$-degrees and effective dimension theory

abstract: In this article, we study an analogue of $tt$-reducibility for points in computable metric spaces. We characterize the notion of the metric $tt$-degree in the context of first-level Borel isomorphism. Then, we study this concept from the perspectives of effective topological dimension theory and of effective fractal dimension theory.


title: Extensions of $\overline\partial$-closed forms on compact generalized Hermitian manifolds

abstract: In this paper, we first get a criterion formula for whether a differential form is holomorphic with respect to the generalized complex structure induced by $\epsilon$. Next, we get the local extensions of $\overline\partial$-closed forms on a smooth family of compact generalized Hermitian manifolds by using this criterion. Finally, as an application, we use this extension to get the invariance of the generalized Hodge number of the deformations of compact generalized Hermitian manifolds with $\partial\overline\partial$-lemma holds.


title: Global stability of a network-based SIRS epidemic model with nonmonotone incidence rate

abstract: This paper studies the dynamics of a network-based SIRS epidemic model with vaccination and a nonmonotone incidence rate. This type of nonlinear incidence can be used to describe the psychological or inhibitory effect from the behavioral change of the susceptible individuals when the number of infective individuals on heterogeneous networks is getting larger. Using the analytical method, epidemic threshold $R_0$ is obtained. When $R_0$ is less than one, we prove the disease-free equilibrium is globally asymptotically stable and the disease dies out, while $R_0$ is greater than one, there exists a unique endemic equilibrium. By constructing a suitable Lyapunov function, we also prove the endemic equilibrium is globally asymptotically stable if the inhibitory factor $\alpha$ is sufficiently large. Numerical experiments are also given to support the theoretical results. It is shown both theoretically and numerically a larger $\alpha$ can accelerate the extinction of the disease and reduce the level of disease.


title: Influence of the Event Rate on Discrimination Abilities of Bankruptcy Prediction Models

abstract: In bankruptcy prediction, the proportion of events is very low, which is often oversampled to eliminate this bias. In this paper, we study the influence of the event rate on discrimination abilities of bankruptcy prediction models. First the statistical association and significance of public records and firmographics indicators with the bankruptcy were explored. Then the event rate was oversampled from 0.12% to 10%, 20%, 30%, 40%, and 50%, respectively. Seven models were developed, including Logistic Regression, Decision Tree, Random Forest, Gradient Boosting, Support Vector Machine, Bayesian Network, and Neural Network. Under different event rates, models were comprehensively evaluated and compared based on Kolmogorov-Smirnov Statistic, accuracy, F1 score, Type I error, Type II error, and ROC curve on the hold-out dataset with their best probability cut-offs. Results show that Bayesian Network is the most insensitive to the event rate, while Support Vector Machine is the most sensitive.


title: Surface wrinkling of an elastic block subject to biaxial loading by an energy method

abstract: Wrinkles are often observed on the surfaces of compressed soft materials in nature. In the past few decades, the fascinating surface patterns have been studied extensively by using the linear bifurcation analysis under plane strain. The bifurcation concerns the non-uniqueness solutions, however, it delivers little information about the surface instability before and after the threshold. In this paper, we study surface wrinkling of a finite elastic block of general elastic materials subject to biaxial loading by an energy method. The first and second variations of the strain energy functional are systematically studied, and an eigenvalue problem is proposed whether the second variation is positive definite. We illustrate our analysis by using neo-Hookean materials as an example. Accordingly, we show that the initially flat state has the lowest energy and is stable before the stretches reach the threshold at which the surface wrinkling occurs. We also find that the threshold is independent of the size of the block and coincides with that of the surface instability of an elastic half-space studied by Biot (1963) with the linear bifurcation analysis. However, the stability region cannot be obtained by using the linear stability analysis. In contrast to the size-independent threshold, the wavelength of surface wrinkling depends on the size of the block. We first show that a two-dimensional rather than a three-dimensional perturbation has lower energy and is more likely to trigger the surface wrinkling in the instability region. The same stretch threshold of a finite block and a half-space could shed light on the relation of surface instabilities between finite and infinite bodies.


title: A Tutorial On Autonomous Vehicle Steering Controller Design, Simulation and Implementation

abstract: In this tutorial, we detailed simple controllers for autonomous parking and path following for self-driving cars providing practical methods for curvature computation.


title: Speech Recognition: Keyword Spotting Through Image Recognition

abstract: The problem of identifying voice commands has always been a challenge due to the presence of noise and variability in speed, pitch, etc. We will compare the efficacies of several neural network architectures for the speech recognition problem. In particular, we will build a model to determine whether a one second audio clip contains a particular word (out of a set of 10), an unknown word, or silence. The models to be implemented are a CNN recommended by the Tensorflow Speech Recognition tutorial, a low-latency CNN, and an adversarially trained CNN. The result is a demonstration of how to convert a problem in audio recognition to the better-studied domain of image classification, where the powerful techniques of convolutional neural networks are fully developed. Additionally, we demonstrate the applicability of the technique of Virtual Adversarial Training (VAT) to this problem domain, functioning as a powerful regularizer with promising potential future applications.


title: Efficient Determination of Equivalence for Encrypted Data

abstract: Secure computation of equivalence has fundamental application in many different areas, including healthcare. We study this problem in the context of matching an individual identity to link medical records across systems. We develop an efficient solution for equivalence based on existing work that can evaluate the greater than relation. We implement the approach and demonstrate its effectiveness on data, as well as demonstrate how it meets regulatory criteria for risk.


title: Quantum-secure message authentication via blind-unforgeability

abstract: Formulating and designing authentication of classical messages in the presence of adversaries with quantum query access has been a longstanding challenge, as the familiar classical notions of unforgeability do not directly translate into meaningful notions in the quantum setting. A particular difficulty is how to fairly capture the notion of "predicting an unqueried value" when the adversary can query in quantum superposition. We propose a natural definition of unforgeability against quantum adversaries called blind unforgeability. This notion defines a function to be predictable if there exists an adversary who can use "partially blinded" oracle access to predict values in the blinded region. We support the proposal with a number of technical results. We begin by establishing that the notion coincides with EUF-CMA in the classical setting and go on to demonstrate that the notion is satisfied by a number of simple guiding examples, such as random functions and quantum-query-secure pseudorandom functions. We then show the suitability of blind unforgeability for supporting canonical constructions and reductions. We prove that the "hash-and-MAC" paradigm and the Lamport one-time digital signature scheme are indeed unforgeable according to the definition. To support our analysis, we additionally define and study a new variety of quantum-secure hash functions called Bernoulli-preserving. Finally, we demonstrate that blind unforgeability is stronger than a previous definition of Boneh and Zhandry [EUROCRYPT '13, CRYPTO '13] in the sense that we can construct an explicit function family which is forgeable by an attack that is recognized by blind-unforgeability, yet satisfies the definition by Boneh and Zhandry.


title: The Dirichlet Casimir Energy for $\phi^4$ Theory in a Rectangle

abstract: In this article, we present the zero and first-order radiative correction to the Dirichlet Casimir energy for massive and massless scalar field confined in a rectangle. This calculation procedure was conducted in two spatial dimensions and for the case of the first-order correction term is new. The renormalization program that we have used in this work, allows all influences from the dominant boundary conditions (e.g. the Dirichlet boundary condition) be automatically reflected in the counterterms. This permission usually makes the counterterms position-dependent. Along with the renormalization program, a supplementary regularization technique was performed in this work. In this regularization technique, that we have named Box Subtraction Scheme (BSS), two similar configurations were introduced and the zero point energies of these two configurations were subtracted from each other using appropriate limits. This regularization procedure makes the usage of any analytic continuation techniques unnecessary. In the present work, first, we briefly present calculation of the leading order Casimir energy for the massive scalar field in a rectangle via BSS. Next, the first order correction to the Casimir energy is calculated by applying the mentioned renormalization and regularization procedures. Finally, all the necessary limits of obtained answers for both massive and massless cases are discussed.


title: \bf $\delta M$ Formalism and Anisotropic Chaotic Inflation Power Spectrum

abstract: A new analytical approach to linear perturbations in anisotropic inflation has been introduced in [A. Talebian-Ashkezari, N. Ahmadi and A.A. Abolhasanib, JCAP 03(2018)001] under the name of $\delta M$ formalism. In this paper we apply the mentioned approach to a model of anisotropic inflation driven by a scalar field, coupled to the kinetic term of a vector field with a $U(1)$ symmetry. The $\delta M$ formalism provides an efficient way of computing tensor- tensor, tensor- scalar as well as scalar- scalar 2-point correlations that are needed for the analysis of the observational features of an anisotropic model on the CMB. A comparison between $\delta M$ results and the tedious calculations using in-in formalism shows the aptitude of the $\delta M$ formalism in calculating accurate two point correlation functions between physical modes of the system.


title: Variance Networks: When Expectation Does Not Meet Your Expectations

abstract: Ordinary stochastic neural networks mostly rely on the expected values of their weights to make predictions, whereas the induced noise is mostly used to capture the uncertainty, prevent overfitting and slightly boost the performance through test-time averaging. In this paper, we introduce variance layers, a different kind of stochastic layers. Each weight of a variance layer follows a zero-mean distribution and is only parameterized by its variance. We show that such layers can learn surprisingly well, can serve as an efficient exploration tool in reinforcement learning tasks and provide a decent defense against adversarial attacks. We also show that a number of conventional Bayesian neural networks naturally converge to such zero-mean posteriors. We observe that in these cases such zero-mean parameterization leads to a much better training objective than conventional parameterizations where the mean is being learned.


title: How to solve the stochastic partial differential equation that gives a Mat\'ern random field using the finite element method

abstract: This tutorial teaches parts of the finite element method (FEM), and solves a stochastic partial differential equation (SPDE). The contents herein are considered "known" in the numerics literature, but for statisticians it is very difficult to find a resource for learning these ideas in a timely manner (without doing a year's worth of courses in numerics). The goal of this tutorial is to be pedagogical and explain the computations/theory to a statistician. This is not a practical tutorial, there is little computer code, and no data analysis.


title: Bright luminescence from indirect and strongly bound excitons in hBN

abstract: A quantitative analysis of the excitonic luminescence efficiency in hexagonal boron nitride (hBN) is carried out by cathodoluminescence in the ultraviolet range and compared with zinc oxide and diamond single crystals. A high quantum yield value of ~50% is found for hBN at 10 K comparable to that of direct bandgap semiconductors. This bright luminescence at 215 nm remains stable up to room temperature, evidencing the strongly bound character of excitons in bulk hBN. Ab initio calculations of the exciton dispersion confirm the indirect nature of the lowest-energy exciton whose binding energy is found equal to 300 meV, in agreement with the thermal stability observed in luminescence. The direct exciton is found at a higher energy but very close to the indirect one, which solves the long debated Stokes shift in bulk hBN.


title: Multi-Agent Submodular Optimization

abstract: Recent years have seen many algorithmic advances in the area of submodular optimization: (SO) $\min/\max~f(S): S \in \mathcal{F}$, where $\mathcal{F}$ is a given family of feasible sets over a ground set $V$ and $f:2^V \rightarrow \mathbb{R}$ is submodular. This progress has been coupled with a wealth of new applications for these models. Our focus is on a more general class of \emph{multi-agent submodular optimization} (MASO) which was introduced by Goel et al. in the minimization setting: $\min \sum_i f_i(S_i): S_1 \uplus S_2 \uplus \cdots \uplus S_k \in \mathcal{F}$. Here we use $\uplus$ to denote disjoint union and hence this model is attractive where resources are being allocated across $k$ agents, each with its own submodular cost function $f_i()$. In this paper we explore the extent to which the approximability of the multi-agent problems are linked to their single-agent {\em primitives}, referred to informally as the {\em multi-agent gap}. We present different reductions that transform a multi-agent problem into a single-agent one. For maximization we show that (MASO) admits an $O(\alpha)$-approximation whenever (SO) admits an $\alpha$-approximation over the multilinear formulation, and thus substantially expanding the family of tractable models. We also discuss several family classes (such as spanning trees, matroids, and $p$-systems) that have a provable multi-agent gap of 1. In the minimization setting we show that (MASO) has an $O(\alpha \cdot \min \{k, \log^2 (n)\})$-approximation whenever (SO) admits an $\alpha$-approximation over the convex formulation. In addition, we discuss the class of "bounded blocker" families where there is a provably tight O$(\log n)$ gap between (MASO) and (SO).


title: Visco-Energetic solutions to some rate-independent systems in damage, delamination, and plasticity

abstract: This paper revolves around a newly introduced weak solvability concept for rate-independent systems, alternative to the notions of Energetic and Balanced Viscosity solutions. Visco-Energetic solutions have been recently obtained by passing to the time-continuous limit in a time-incremental scheme, akin to that for Energetic solutions, but perturbed by a `viscous' correction term, as in the case of Balanced Viscosity solutions. However, for Visco-Energetic solutions this viscous correction is tuned by a fixed parameter. The resulting solution notion turns out to describe a kind of evolution in between Energetic and Balanced Viscosity evolution. In this paper we aim to investigate the application of Visco-Energetic solutions to the paradigmatic example of perfect plasticity, and to nonsmooth rate-independent processes in solid mechanics such as damage and plasticity at finite strains. With the limit passage from adhesive contact to brittle delamination, we also provide a first result of Evolutionary Gamma-convergence for Visco-Energetic solutions. The analysis of these applications reveals the wide applicability of this solution concept and confirms its intermediate character.


title: A Minimax Surrogate Loss Approach to Conditional Difference Estimation

abstract: We present a new machine learning approach to estimate personalized treatment effects in the classical potential outcomes framework with binary outcomes. To overcome the problem that both treatment and control outcomes for the same unit are required for supervised learning, we propose surrogate loss functions that incorporate both treatment and control data. The new surrogates yield tighter bounds than the sum of losses for treatment and control groups. A specific choice of loss function, namely a type of hinge loss, yields a minimax support vector machine formulation. The resulting optimization problem requires the solution to only a single convex optimization problem, incorporating both treatment and control units, and it enables the kernel trick to be used to handle nonlinear (also non-parametric) estimation. Statistical learning bounds are also presented for the framework, and experimental results.


title: Continuous solutions of a second order iterative equation

abstract: In this paper we study the existence of continuous solutions and their constructions for a second order iterative functional equation, which involves iterate of the unknown function and a nonlinear term. Imposing Lipschitz conditions to those given functions, we prove the existence of continuous solutions on the whole $\mathbb{R}$ by applying the contraction principle. In the case without Lipschitz conditions we hardly use the contraction principle, but we construct continuous solutions on $\mathbb{R}$ recursively with a partition of $\mathbb{R}$.