title: Bilinear Assignment Problem: Large Neighborhoods and Experimental Analysis of Algorithms

abstract: The bilinear assignment problem (BAP) is a generalization of the well-known quadratic assignment problem (QAP). In this paper, we study the problem from the computational analysis point of view. Several classes of neigborhood structures are introduced for the problem along with some theoretical analysis. These neighborhoods are then explored within a local search and a variable neighborhood search frameworks with multistart to generate robust heuristic algorithms. Results of systematic experimental analysis have been presented which divulge the effectiveness of our algorithms. In addition, we present several very fast construction heuristics. Our experimental results disclosed some interesting properties of the BAP model, different from those of comparable models. This is the first thorough experimental analysis of algorithms on BAP. We have also introduced benchmark test instances that can be used for future experiments on exact and heuristic algorithms for the problem.


title: A Space-Time Cut Finite Element Method with quadrature in time

abstract: We consider convection-diffusion problems in time-dependent domains and present a space-time finite element method based on quadrature in time which is simple to implement and avoids remeshing procedures as the domain is moving. The evolving domain is embedded in a domain with fixed mesh and a cut finite element method with continuous elements in space and discontinuous elements in time is proposed. The method allows the evolving geometry to cut through the fixed background mesh arbitrarily and thus avoids remeshing procedures. However, the arbitrary cuts may lead to ill-conditioned algebraic systems. A stabilization term is added to the weak form which guarantees well-conditioned linear systems independently of the position of the geometry relative to the fixed mesh and in addition makes it possible to use quadrature rules in time to approximate the space-time integrals. We review here the space-time cut finite element method presented in [13] where linear elements are used in both space and time and extend the method to higher order elements for problems on evolving surfaces (or interfaces). We present a new stabilization term which also when higher order elements are used controls the condition number of the linear systems from cut finite element methods on evolving surfaces. The new stabilization combines the consistent ghost penalty stabilization [1] with a term controlling normal derivatives at the interface.


title: Anti-D3 branes and moduli in non-linear supergravity

abstract: Anti-D3 branes and non-perturbative effects in flux compactifications spontaneously break supersymmetry and stabilise moduli in a metastable de Sitter vacua. The low energy 4D effective field theory description for such models would be a supergravity theory with non-linearly realised supersymmetry. Guided by string theory modular symmetry, we compute this non-linear supergravity theory, including dependence on all bulk moduli. Using either a constrained chiral superfield or a constrained vector field, the uplifting contribution to the scalar potential from the anti-D3 brane can be parameterised either as an F-term or Fayet-Iliopoulos D-term. Using again the modular symmetry, we show that 4D non-linear supergravities that descend from string theory have an enhanced protection from quantum corrections by non-renormalisation theorems. The superpotential giving rise to metastable de Sitter vacua is robust against perturbative string-loop and $\alpha'$ corrections.


title: Transport of hot carriers in plasmonic nanostructures

abstract: Plasmonic hot carrier devices extract excited carriers from metal nanostructures before equilibration, and have the potential to surpass semiconductor light absorbers. However their efficiencies have so far remained well below theoretical limits, which necessitates quantitative prediction of carrier transport and energy loss in plasmonic structures to identify and overcome bottlenecks in carrier harvesting. Here, we present a theoretical and computational framework, Non-Equilibrium Scattering in Space and Energy (NESSE), to predict the spatial evolution of carrier energy distributions that combines the best features of phase-space (Boltzmann) and particle-based (Monte Carlo) methods. Within the NESSE framework, we bridge first-principles electronic structure predictions of plasmon decay and carrier collision integrals at the atomic scale, with electromagnetic field simulations at the nano- to mesoscale. Finally, we apply NESSE to predict spatially-resolved energy distributions of photo-excited carriers that impact the surface of experimentally realizable plasmonic nanostructures at length scales ranging from tens to several hundreds of nanometers, enabling first-principles design of hot carrier devices.

categories: physics.comp-ph cond-mat.mes-hall cond-mat.mtrl-sci physics.optics


title: On semistability of $CAT(0)$ groups

abstract: Does every one-ended $CAT(0)$ group have semistable fundamental group at infinity? As we write, this is an open question. Let $G$ be such a group acting geometrically on the proper $CAT(0)$ space $X$. In this paper we show that in order to establish a positive answer to the question it is only necessary to check that any two geodesic rays in $X$ are properly homotopic. We then show that if the answer to the question is negative, with $(G,X)$ a counter-example, then the boundary of $X$, $\del X$ with the cone topology, must have a weak cut point. This is of interest because a theorem of Papasoglu and the second-named author \cite{PS} has established that there cannot be an example of $(G,X)$ where $\del X$ has a cut point. Thus, the search for a negative answer comes down to the difference between cut points and weak cut points. We also show that the Tits ball of radius $\frac{\pi}{2}$ about that weak cut point is a "cut set" in the sense that it separates $\del X$. Finally, we observe that if a negative example $(G, X)$ exists then $G$ is rank 1.


title: A Pilot Study of Domain Adaptation Effect for Neural Abstractive Summarization

abstract: We study the problem of domain adaptation for neural abstractive summarization. We make initial efforts in investigating what information can be transferred to a new domain. Experimental results on news stories and opinion articles indicate that neural summarization model benefits from pre-training based on extractive summaries. We also find that the combination of in-domain and out-of-domain setup yields better summaries when in-domain data is insufficient. Further analysis shows that, the model is capable to select salient content even trained on out-of-domain data, but requires in-domain data to capture the style for a target domain.


title: Entanglement of a class of non-Gaussian states in disordered harmonic oscillator systems

abstract: For disordered harmonic oscillator systems over the $d$-dimensional lattice, we consider the problem of finding the bipartite entanglement of the uniform ensemble of the energy eigenstates associated with a particular number of modes. Such ensemble define a class of mixed, non-Gaussian entangled states that are labeled, by the energy of the system, in an increasing order. We develop a novel approach to find the exact logarithmic negativity of this class of states. We also prove entanglement bounds and demonstrate that the low energy states follow an area law.


title: Exciton condensate in bilayer transition metal dichalcogenides: strong coupling regime

abstract: Exciton condensation in an electron-hole bilayer system of monolayer transition metal dichalcogenides is analyzed at three different levels of theory to account for screening and quasiparticle renormalization. The large effective masses of the transition metal dichalcogenides place them in a strong coupling regime. In this regime, mean field (MF) theory with either an unscreened or screened interlayer interaction predicts a room temperature condensate. Interlayer and intralayer interactions renormalize the quasiparticle dispersion, and this effect is included in a GW approximation. The renormalization reverses the trends predicted from the unscreened or screened MF theories. In the strong coupling regime, intralayer interactions have a large impact on the magnitude of the order parameter and its functional dependencies on effective mass and carrier density.


title: Modeling alignment enhancement for solid polarized targets

abstract: A model of dynamic orientation using optimized radiofrequency (RF) irradiation produced perpendicular to the holding field is developed for the spin-1 system required for tensor-polarized fixed-target experiments. The derivation applies to RF produced close to the Larmor frequency of the nucleus and requires the electron spin-resonance linewidth to be much smaller than the nuclear magnetic resonance frequency. The rate equations are solved numerically to study a semi-saturated steady-state resulting from the two sources of irradiation: microwave from the DNP process and the additional RF used to manipulate the tensor polarization. The steady-state condition and continuous-wave NMR lineshape are found that optimize the spin-1 alignment in the polycrystalline materials used as solid polarized targets in charged-beam nuclear and particle physics experiments.


title: Ultraslow diffusion in language: Dynamics of appearance of already popular adjectives on Japanese blogs

abstract: What dynamics govern a time series representing the appearance of words in social media data? In this paper, we investigate an elementary dynamics, from which word-dependent special effects are segregated, such as breaking news, increasing (or decreasing) concerns, or seasonality. To elucidate this problem, we investigated approximately three billion Japanese blog articles over a period of six years, and analysed some corresponding solvable mathematical models. From the analysis, we found that a word appearance can be explained by the random diffusion model based on the power-law forgetting process, which is a type of long memory point process related to ARFIMA(0,0.5,0). In particular, we confirmed that ultraslow diffusion (where the mean squared displacement grows logarithmically), which the model predicts in an approximate manner, reproduces the actual data. In addition, we also show that the model can reproduce other statistical properties of a time series: (i) the fluctuation scaling, (ii) spectrum density, and (iii) shapes of the probability density functions.


title: Multipath Multiplexing for Capacity Enhancement in SIMO Wireless Systems

abstract: This paper proposes a novel and simple orthogonal faster than Nyquist (OFTN) data transmission and detection approach for a single input multiple output (SIMO) system. It is assumed that the signal having a bandwidth $B$ is transmitted through a wireless channel with $L$ multipath components. Under this assumption, the current paper provides a novel and simple OFTN transmission and symbol-by-symbol detection approach that exploits the multiplexing gain obtained by the multipath characteristic of wideband wireless channels. It is shown that the proposed design can achieve a higher transmission rate than the existing one (i.e., orthogonal frequency division multiplexing (OFDM)). Furthermore, the achievable rate gap between the proposed approach and that of the OFDM increases as the number of receiver antennas increases for a fixed value of $L$. This implies that the performance gain of the proposed approach can be very significant for a large-scale multi-antenna wireless system. The superiority of the proposed approach is shown theoretically and confirmed via numerical simulations. {Specifically, we have found {upper-bound average} rates of 15 bps/Hz and 28 bps/Hz with the OFDM and proposed approaches, respectively, in a Rayleigh fading channel with 32 receive antennas and signal to noise ratio (SNR) of 15.3 dB. The extension of the proposed approach for different system setups and associated research problems is also discussed.


title: Adaptive Channel Prediction, Beamforming and Scheduling Design for 5G V2I Network

abstract: One of the important use-cases of 5G network is the vehicle to infrastructure (V2I) communication which requires accurate understanding about its dynamic propagation environment. As 5G base stations (BSs) tend to have multiple antennas, they will likely employ beamforming to steer their radiation pattern to the desired vehicle equipment (VE). Furthermore, since most wireless standards employ an OFDM system, each VE may use one or more sub-carriers. To this end, this paper proposes a joint design of adaptive channel prediction, beamforming and scheduling for 5G V2I communications. The channel prediction algorithm is designed without the training signal and channel impulse response (CIR) model. In this regard, first we utilize the well known adaptive recursive least squares (RLS) technique for predicting the next block CIR from the past and current block received signals (a block may have one or more OFDM symbols). Then, we jointly design the beamforming and VE scheduling for each sub-carrier to maximize the uplink channel average sum rate by utilizing the predicted CIR. The beamforming problem is formulated as a Rayleigh quotient optimization where its global optimal solution is guaranteed. And, the VE scheduling design is formulated as an integer programming problem which is solved by employing a greedy search. The superiority of the proposed channel prediction and scheduling algorithms over those of the existing ones is demonstrated via numerical simulations.


title: An explanation of remarkable emission line profiles in post-flare coronal rain

abstract: We study broad red-shifted emission in chromospheric and transition region lines that appears to correspond to a form of post-flare coronal rain. Profiles of Mg II, C II and Si IV lines were obtained using the IRIS instrument before, during and after the X2.1 flare of 11 March 2015 (SOL2015-03- 11T16:22). We analyze the profiles of the five transitions of Mg II (the 3p - 3s h and k transitions, and three lines belonging to the 3d - 3p transitions). We use analytical methods to understand the unusual profiles, together with higher resolution observational data of similar phenomena observed by Jing et al. (2016). The peculiar line ratios indicate anisotropic emission from the strands which have cross-strand line center optical depths (k-line) of between 1 and 10. The lines are broadened by unresolved Alfvenic motions whose energy exceeds the radiation losses in the Mg II lines by an order of magnitude. The decay of the line widths is accompanied by a decay in the brightness, suggesting a causal connection. If the plasma is <~ 99% ionized, ion-neutral collisions can account for the dissipation, otherwise a of dynamical process seems necessary. Our work implies that the motions are initiated during the impulsive phase, to be dissipated as radiation over a period of an hour, predominantly by strong chromospheric lines. The coronal "rain" we observe is far more turbulent that most earlier reports have indicated, with implications for plasma heating mechanisms.


title: Steklov Spectral Geometry for Extrinsic Shape Analysis

abstract: We propose using the Dirichlet-to-Neumann operator as an extrinsic alternative to the Laplacian for spectral geometry processing and shape analysis. Intrinsic approaches, usually based on the Laplace-Beltrami operator, cannot capture the spatial embedding of a shape up to rigid motion, and many previous extrinsic methods lack theoretical justification. Instead, we consider the Steklov eigenvalue problem, computing the spectrum of the Dirichlet-to-Neumann operator of a surface bounding a volume. A remarkable property of this operator is that it completely encodes volumetric geometry. We use the boundary element method (BEM) to discretize the operator, accelerated by hierarchical numerical schemes and preconditioning; this pipeline allows us to solve eigenvalue and linear problems on large-scale meshes despite the density of the Dirichlet-to-Neumann discretization. We further demonstrate that our operators naturally fit into existing frameworks for geometry processing, making a shift from intrinsic to extrinsic geometry as simple as substituting the Laplace-Beltrami operator with the Dirichlet-to-Neumann operator.


title: Complete convergence and records for dynamically generated stochastic processes

abstract: We consider empirical multi-dimensional Rare Events Point Processes that keep track both of the time occurrence of extremal observations and of their severity, for stochastic processes arising from a dynamical system, by evaluating a given potential along its orbits. This is done both in the absence and presence of clustering. A new formula for the piling of points on the vertical direction of bi-dimensional limiting point processes, in the presence of clustering, is given, which is then generalised for higher dimensions. The limiting multi-dimensional processes are computed for systems with sufficiently fast decay of correlations. The complete convergence results are used to study the effect of clustering on the convergence of extremal processes, record time and record values point processes. An example where the clustering prevents the convergence of the record times point process is given.


title: Are all flare ribbons simply connected to the corona?

abstract: We consider the observational basis for the belief that flare ribbons in the chromosphere result from energy transport from the overlying corona. We study ribbons of small flares using magnetic as well as intensity data from the Hinode, SDO and IRIS missions. While most ribbons appear connected to the corona, and they over-lie regions of significant vertical magnetic field, we examine one ribbon with no clear evidence for such connections. Evolving horizontal magnetic fields seen with Hinode suggest that reconnection with pre-existing fields below the corona can explain the data. The identification of just one, albeit small, ribbon, with no apparent connection to the corona, leads us to conclude that at least two mechanisms are responsible for the heating that leads to flare ribbon emission.


title: High-$j$ neutron excitations outside $^{136}$Xe

abstract: The $\nu0h_{9/2}$ and $\nu0i_{13/2}$ strength at $^{137}$Xe, a single neutron outside the $N=82$ shell closure, has been determined using the $^{136}$Xe($\alpha$,$^3$He)$^{137}$Xe reaction carried out at 100 MeV. We confirm the recent observation of the second 13/2$^+$ state and reassess previous data on the 9/2$^-$ states, obtaining spectroscopic factors. These new data provide additional constraints on predictions of the same single-neutron excitations at $^{133}$Sn.


title: What-and-Where to Match: Deep Spatially Multiplicative Integration Networks for Person Re-identification

abstract: Matching pedestrians across disjoint camera views, known as person re-identification (re-id), is a challenging problem that is of importance to visual recognition and surveillance. Most existing methods exploit local regions within spatial manipulation to perform matching in local correspondence. However, they essentially extract \emph{fixed} representations from pre-divided regions for each image and perform matching based on the extracted representation subsequently. For models in this pipeline, local finer patterns that are crucial to distinguish positive pairs from negative ones cannot be captured, and thus making them underperformed. In this paper, we propose a novel deep multiplicative integration gating function, which answers the question of \emph{what-and-where to match} for effective person re-id. To address \emph{what} to match, our deep network emphasizes common local patterns by learning joint representations in a multiplicative way. The network comprises two Convolutional Neural Networks (CNNs) to extract convolutional activations, and generates relevant descriptors for pedestrian matching. This thus, leads to flexible representations for pair-wise images. To address \emph{where} to match, we combat the spatial misalignment by performing spatially recurrent pooling via a four-directional recurrent neural network to impose spatial dependency over all positions with respect to the entire image. The proposed network is designed to be end-to-end trainable to characterize local pairwise feature interactions in a spatially aligned manner. To demonstrate the superiority of our method, extensive experiments are conducted over three benchmark data sets: VIPeR, CUHK03 and Market-1501.


title: Automatic Curation of Golf Highlights using Multimodal Excitement Features

abstract: The production of sports highlight packages summarizing a game's most exciting moments is an essential task for broadcast media. Yet, it requires labor-intensive video editing. We propose a novel approach for auto-curating sports highlights, and use it to create a real-world system for the editorial aid of golf highlight reels. Our method fuses information from the players' reactions (action recognition such as high-fives and fist pumps), spectators (crowd cheering), and commentator (tone of the voice and word analysis) to determine the most interesting moments of a game. We accurately identify the start and end frames of key shot highlights with additional metadata, such as the player's name and the hole number, allowing personalized content summarization and retrieval. In addition, we introduce new techniques for learning our classifiers with reduced manual training data annotation by exploiting the correlation of different modalities. Our work has been demonstrated at a major golf tournament, successfully extracting highlights from live video streams over four consecutive days.


title: Effective one-component model of binary mixture: molecular arrest induced by the spatially correlated stochastic dynamics

abstract: Spatially correlated noise (SCN), i.e. the thermal noise that affects neighbouring particles in a similar manner, is ubiquitous in soft matter systems. In this work, we apply the over-damped SCN-driven Langevin equations as an effective, one-component model of the dynamics in dense binary mixtures. We derive the thermodynamically consistent fluctuation-dissipation relation for SCN to show that it predicts the molecular arrest resembling the glass transition, i.e. the critical slow-down of dynamics in the disordered phases. We show that the mechanism of singular dissipation is embedded in the dissipation matrix, accompanying SCN. We are also able to identify the characteristic length of collective dissipation, which diverges at critical packing. This novel physical quantity conveniently describes the difference between the ergodic and non-ergodic dynamics. The model is fully analytically solvable, one-dimensional and admits arbitrary interactions between the particles. It qualitatively reproduces several different modes of arrested disorder encountered in binary mixtures, including e.g. the re-entrant arrest. The model can be effectively compared to the mode coupling theory.


title: On the Fine Structure Splitting of the 3p43d 4D5/2 and 3p43d 4D7/2 Levels of Fe X

abstract: We study UV spectra obtained with the SO82-B slit spectrograph on board SKYLAB to estimate the fine structure splitting of the Cl-like 3p4 3d 4D J=5/2 and 3p4 3d 4D J=7/2 levels of Fe X. The splitting is of interest because the Zeeman effect mixes these levels, producing a "magnetically induced transition" (MIT) from 3p4 3d 4D J=7/2 to 3p5 2Po J=3/2 for modest magnetic field strengths characteristic of the active solar corona. We estimate the splitting using the Ritz combination formula applied to two lines in the UV region of the spectrum close to 1603.2 Angstrom, which decay from the level 3p4(1D)3d 2G J=7/2 to these two lower levels. The MIT and accompanying spin-forbidden transition lie near 257 Angstrom. By careful inspection of a deep exposure obtained with the S082B instrument we derive a splitting of <~ 7 +/- 3 cm-1. The upper limit arises because of a degeneracy between the effects of non-thermal line broadening and fine-structure splitting for small values of the latter parameter. Although the data were recorded on photographic film, we solved for optimal values of line width and splitting of 8.3 +/- 0.9 and 3.6 +/- 2.7 cm-1.


title: Ergodic Mean Field Games with H\"ormander diffusions

abstract: We prove existence of solutions for a class of systems of subelliptic PDEs arising from Mean Field Game systems with H\"ormander diffusion. These results are motivated by the feedback synthesis Mean Field Game solutions and the Nash equilibria of a large class of $N$-player differential games.


title: Existence and nonexistence of positive solutions to some fully nonlinear equation in one dimension

abstract: In this paper, we consider the existence (and nonexistence) of solutions to \[ -\mathcal{M}_{\lambda,\Lambda}^\pm (u'') + V(x) u = f(u) \quad {\rm in} \ \mathbf{R} \] where $\mathcal{M}_{\lambda,\Lambda}^+$ and $\mathcal{M}_{\lambda,\Lambda}^-$ denote the Pucci operators with $0< \lambda \leq \Lambda < \infty$, $V(x)$ is a bounded function, $f(s)$ is a continuous function and its typical example is a power-type nonlinearity $f(s) =|s|^{p-1}s$ $(p>1)$. In particular, we are interested in positive solutions which decay at infinity, and the existence (and nonexistence) of such solutions is proved.


title: The Interplay of Competition and Cooperation Among Service Providers

abstract: We consider the economics of the interaction between Mobile Virtual Network Operators (MVNOs) and Mobile Network Operators (MNOs). We investigate the incentives of an MNO for offering some of her resources to an MVNO instead of using the resources for her own End-Users (EUs). We consider a market with one MNO and one MVNO, and a continuum of undecided EUs. Two cases for EUs are considered: (i) when EUs need to choose either the MNO or the MVNO, and (ii) when EUs have an outside option. In each of these cases, we consider a non-cooperative framework of sequential game and a cooperative framework of bargaining game. We characterize the Subgame Perfect Nash Equilibria (SPNE) and Nash Bargaining Solution (NBS) of the sequential and bargaining games, respectively. We show that in the non-cooperative framework, SPNE assumes two forms: (i) the MNO invests minimally on her infrastructure and the MVNO leases all the newly invested resources and (ii) the MNO invests more on her infrastructure, the MVNO leases only part of these resources to the MVNO, and the MNO uses the rest herself to attract EUs. Thus, in both, the MNO generates revenue indirectly through the MVNO. In addition to that, in (ii), the MNO generates revenue directly from the EUs. We also prove that in the bargaining framework, the MVNO either reserves all the resources or no resources from the MNO, and the MNO's investments are guided by whether EUs have an outside option. If they don't, then the MNO invests as little as possible on her infrastructure. If they do, then the MNO invests more.


title: Measurement of Atmospheric Neutrino Oscillations at 6-56 GeV with IceCube DeepCore

abstract: We present a measurement of the atmospheric neutrino oscillation parameters using three years of data from the IceCube Neutrino Observatory. The DeepCore infill array in the center of IceCube enables detection and reconstruction of neutrinos produced by the interaction of cosmic rays in the Earth's atmosphere at energies as low as $\sim5$ GeV. That energy threshold permits measurements of muon neutrino disappearance, over a range of baselines up to the diameter of the Earth, probing the same range of $L/E_\nu$ as long-baseline experiments but with substantially higher energy neutrinos. This analysis uses neutrinos from the full sky with reconstructed energies from $5.6$ - $56$ GeV. We measure $\Delta m^2_{32}=2.31^{+0.11}_{-0.13} \times 10^{-3}$ eV$^2$ and $\sin^2 \theta_{23}=0.51^{+0.07}_{-0.09}$, assuming normal neutrino mass ordering. These results are consistent with, and of similar precision to, those from accelerator and reactor-based experiments.


title: Gyroscope Calibration via Magnetometer

abstract: Magnetometers, gyroscopes and accelerometers are commonly used sensors in a variety of applications. The paper proposes a novel gyroscope calibration method in the homogeneous magnetic field by the help of magnetometer. It is shown that, with sufficient rotation excitation, the homogeneous magnetic field vector can be exploited to serve as a good reference for calibrating low-cost gyroscopes. The calibration parameters include the gyroscope scale factor, non-orthogonal coefficient and bias for three axes, as well as its misalignment to the magnetometer frame. Simulation and field test results demonstrate the method's effectiveness.


title: Computing the $k$-resilience of a Synchronized Multi-Robot System

abstract: We study an optimization problem that arises in the design of covering strategies for multi-robot systems. Consider a team of $n$ cooperating robots traveling along predetermined closed and disjoint trajectories. Each robot needs to periodically communicate information to nearby robots. At places where two trajectories are within range of each other, a communication link is established, allowing two robots to exchange information, provided they are "synchronized", i.e., they visit the link at the same time. In this setting a communication graph is defined and a system of robots is called \emph{synchronized} if every pair of neighbors is synchronized. If one or more robots leave the system, then some trajectories are left unattended. To handle such cases in a synchronized system, when a live robot arrives to a communication link and detects the absence of the neighbor, it shifts to the neighboring trajectory to assume the unattended task. If enough robots leave, it may occur that a live robot enters a state of \emph{starvation}, failing to permanently meet other robots during flight. To measure the tolerance of the system under this phenomenon we define the \emph{$k$-resilience} as the minimum number of robots whose removal may cause $k$ surviving robots to enter a state of starvation. We show that the problem of computing the $k$-resilience is NP-hard if $k$ is part of the input, even if the communication graph is a tree. We propose algorithms to compute the $k$-resilience for constant values of $k$ in general communication graphs and show more efficient algorithms for systems whose communication graph is a tree.


title: A Framework for Finding Anomalous Objects at the LHC

abstract: Search for new physics events at the LHC mostly rely on the assumption that the events are characterized in terms of standard-reconstructed objects such as isolated photons, leptons, and jets initiated by QCD-partons. While such strategy works for a vast majority of physics beyond the standard model scenarios, there are examples aplenty where new physics give rise to anomalous objects (such as collimated and equally energetic particles, decays due to long lived particles etc.) in the detectors, which can not be classified as any of the standard-objects. Varied methods and search strategies have been proposed, each of which is trained and optimized for specific models, topologies, and model parameters. Further, as LHC keeps excluding all expected candidates for new physics, the need for a generic method/tool that is capable of finding the unexpected can not be understated. In this paper, we propose one such method that relies on the philosophy that all anomalous objects are $\it{not}$ standard-objects. The anomaly finder, we suggest, simply is a collection of vetoes that eliminate all standard-objects up to a pre-determined acceptance rate. Any event containing at least one anomalous object (that passes all these vetoes), can be identified as a candidate for new physics. Subsequent offline analyses can determine the nature of the anomalous object as well as of the event, paving a robust way to search for these new physics scenarios in a model-independent fashion. Further, since the method relies on learning only the standard-objects, for which control samples are readily available from data, one can build the analysis in an entirely data-driven way.


title: Acceleration of Cosmic Ray Electrons at Weak Shocks in Galaxy Clusters

abstract: According to structure formation simulations, weak shocks with typical Mach number, $M_{\rm s}\lesssim 3$, are expected to form in merging galaxy clusters. The presence of such shocks has been indicated by X-ray and radio observations of many merging clusters. In particular, diffuse radio sources known as radio relics could be explained by synchrotron-emitting electrons accelerated via diffusive shock acceleration (Fermi I) at quasi-perpendicular shocks. Here we also consider possible roles of stochastic acceleration (Fermi II) by compressive MHD turbulence downstream of the shock. Then we explore a puzzling discrepancy that for some radio relics, the shock Mach number inferred from the radio spectral index is substantially larger than that estimated from X-ray observations. This problem could be understood, if shock surfaces associated with radio relics consist of multiple shocks with different strengths.In that case, X-ray observations tend to pick up the part of shocks with lower Mach numbers and higher kinetic energy flux, while radio emissions come preferentially from the part of shocks with higher Mach numbers and higher cosmic ray (CR) production. We also show that the Fermi I reacceleration model with preexisting fossil electrons supplemented by Fermi II acceleration due to postshock turbulence could reproduce observed profiles of radio flux densities and integrated radio spectra of two giant radio relics. This study demonstrates the CR electrons can be accelerated at collisionless shocks in galaxy clusters just like supernova remnant shock in the interstellar medium and interplanetary shocks in the solar wind.


title: Identifying civilians killed by police with distantly supervised entity-event extraction

abstract: We propose a new, socially-impactful task for natural language processing: from a news corpus, extract names of persons who have been killed by police. We present a newly collected police fatality corpus, which we release publicly, and present a model to solve this problem that uses EM-based distant supervision with logistic regression and convolutional neural network classifiers. Our model outperforms two off-the-shelf event extractor systems, and it can suggest candidate victim names in some cases faster than one of the major manually-collected police fatality databases.


title: Modified mean curvature flow of entire locally Lipschitz radial graphs in hyperbolic space

abstract: The asymptotic Plateau problem asks for the existence of smooth complete hypersurfaces of constant mean curvature with prescribed asymptotic boundary at infinity in the hyperbolic space $\mathbb{H}^{n+1}$. The modified mean curvature flow (MMCF) was firstly introduced by Xiao and the second author a few years back, and it provides a tool using geometric flow to find such hypersurfaces with constant mean curvature in $\mathbb{H}^{n+1}$. Similar to the usual mean curvature flow, the MMCF is the natural negative $L^2$-gradient flow of the area-volume functional $\mathcal{I}(\Sigma)=A(\Sigma)+\sigma V(\Sigma)$ associated to a hypersurface $\Sigma$. In this paper, we prove that the MMCF starting from an entire locally Lipschitz continuous radial graph exists and stays radially graphic for all time. In general one cannot expect the convergence of the flow as it can be seen from the flow starting from a horosphere (whose asymptotic boundary is degenerate to a point).


title: On LS-category and topological complexity of connected sum

abstract: The Lusternik-Schnirelmann category and topological complexity are important invariants of manifolds (and more generally, topological spaces). We study the behavior of these invariants under the operation of taking the connected sum of manifolds. We give a complete answer for the LS-categoryof orientable manifolds, $\cat(M\# N)=\max\{\cat M,\cat N\}$. For topological complexity we prove the inequality $\TC (M\# N)\ge\max\{\TC M,\TC N\}$ for simply connected manifolds.


title: Motion Compensated Dynamic MRI Reconstruction with Local Affine Optical Flow Estimation

abstract: This paper proposes a novel framework to reconstruct the dynamic magnetic resonance images (DMRI) with motion compensation (MC). Due to the inherent motion effects during DMRI acquisition, reconstruction of DMRI using motion estimation/compensation (ME/MC) has been studied under a compressed sensing (CS) scheme. In this paper, by embedding the intensity-based optical flow (OF) constraint into the traditional CS scheme, we are able to couple the DMRI reconstruction with motion field estimation. The formulated optimization problem is solved by a primal-dual algorithm with linesearch due to its efficiency when dealing with non-differentiable problems. With the estimated motion field, the DMRI reconstruction is refined through MC. By employing the multi-scale coarse-to-fine strategy, we are able to update the variables(temporal image sequences and motion vectors) and to refine the image reconstruction alternately. Moreover, the proposed framework is capable of handling a wide class of prior information (regularizations) for DMRI reconstruction, such as sparsity, low rank and total variation. Experiments on various DMRI data, ranging from in vivo lung to cardiac dataset, validate the reconstruction quality improvement using the proposed scheme in comparison to several state-of-the-art algorithms.


title: Atomic physics and modern solar spectro-polarimetry

abstract: Observational solar physics is entering a new era with the advent of new 1.5 m class telescopes with adaptive optics, as well as the Daniel K. Inouye 4 m telescope which will become operational in 2019. Major outstanding problems in solar physics all relate to the solar magnetic field. Spectropolarimetry offers the best, and sometimes only, method for accurate measurements of the magnetic field. In this paper we highlight how certain atomic transitions can help us provide both calibration data, as well as diagnostic information on solar magnetic fields, in the presence of residual image distortions through the atmosphere close to, but not at the diffraction limits of large and polarizing telescopes. Particularly useful are spectral lines of neutrals and singly charged ions of iron and other complex atoms. As a proof-of-concept, we explore atomic transitions that might be used to study magnetic fields without the need for an explicit calibration sequence, offering practical solutions to the difficult challenges of calibrating the next generation of solar spectropolarimetric telescopes. Suggestions for additional work on atomic theory and measurements, particularly at infrared wavelengths, are given. There is some promise for continued symbiotic advances between solar physics and atomic physics.


title: Free arrangements with low exponents

abstract: In this article we show that any free hyperplane arrangement with exponents 1's and 2's is a supersolvable arrangement. We conjecture that any free arrangement with exponents 1's, 2's and exactly one 3, is also supersolvable, and we show this conjecture for hyperplane arrangements of ranks 4 and 5, and for inductively free arrangements of any rank.


title: On invariance of plurigenera for foliated surface pairs

abstract: We prove an invariance of plurigenera for some foliated surface pairs of general type.


title: HATS-43b, HATS-44b, HATS-45b, and HATS-46b: Four Short Period Transiting Giant Planets in the Neptune-Jupiter Mass Range

abstract: We report the discovery of four short period extrasolar planets transiting moderately bright stars from photometric measurements of the HATSouth network coupled to additional spectroscopic and photometric follow-up observations. While the planet masses range from 0.26 to 0.90 M$_J$, the radii are all approximately a Jupiter radii, resulting in a wide range of bulk densities. The orbital period of the planets range from 2.7d to 4.7d, with HATS-43b having an orbit that appears to be marginally non-circular (e= 0.173$\pm$0.089). HATS-44 is notable for a high metallicity ([Fe/H]= 0.320$\pm$0.071). The host stars spectral types range from late F to early K, and all of them are moderately bright (13.3<V<14.4), allowing the execution of future detailed follow-up observations. HATS-43b and HATS-46b, with expected transmission signals of 2350 ppm and 1500 ppm, respectively, are particularly well suited targets for atmospheric characterisation via transmission spectroscopy.


title: Hybrid Voltage Control in Distribution Networks Under Limited Communication Rates

abstract: Voltage regulation in distribution networks is challenged by increasing penetration of distributed energy resources (DERs). Thanks to advancement in power electronics, these DERs can be leveraged to regulate the grid voltage by quickly changing the reactive power outputs. This paper develops a hybrid voltage control (HVC) strategy that can seamlessly integrate both local and distributed designs to coordinate the network-wide reactive power resources from DERs. \ws{By minimizing a special voltage mismatch objective, we achieve the proposed HVC architecture using partial primal-dual (PPD) gradient updates that allow for a distributed and online implementation}. The proposed HVC design improves over existing distributed approaches by integrating with local voltage feedback. As a result, it can dynamically adapt to varying system operating conditions while being fully cognizant to the instantaneous availability of communication links. Under the worst-case scenario of a total link outage, the proposed design naturally boils down to a surrogate local control implementation. Numerical tests on realistic feeder cases have been to corroborate our analytical results and demonstrate the algorithmic performance.


title: Counting Conjugacy Classes in $Out(F_N)$

abstract: We show that if a f.g. group $G$ has a non-elementary WPD action on a hyperbolic metric space $X$, then the number of $G$-conjugacy classes of $X$-loxodromic elements of $G$ coming from a ball of radius $R$ in the Cayley graph of $G$ grows exponentially in $R$. As an application we prove that for $N\ge 3$ the number of distinct $Out(F_N)$-conjugacy classes of fully irreducibles $\phi$ from an $R$-ball in the Cayley graph of $Out(F_N)$ with $\log\lambda(\phi)$ on the order of $R$ grows exponentially in $R$.


title: Dynamical phase diagram of parity-time symmetry with competing saturable channels

abstract: Nonlinear channels play a critical role in realizing dynamical functions. Neural ionic channels and non-volatile memristors each derive representative biological and electrical functionalities, such as repetitive firing or pinched hysteresis. In electromagnetics, saturable channels of amplification or absorption provide a large nonlinearity for nonequilibrium wave dynamics, from conventional lasing to mode locking to recent achievements of the non-reciprocity in complex potentials. Here, we investigate the dynamical phase diagram of parity-time symmetric systems, governed by competing nonlinear channels of saturable amplification and absorption. Determined by the relative strength and saturation level of the channels, three distinctive phases of fast- and slow-response equilibriums, and an oscillating nonequilibrium are demonstrated. On phase boundaries, we also reveal the chaotic existence of the strong oscillation state, which allows the non-reciprocal realization of repetitive resonator firing with fully tunable time delays. This work will promote the wave-based realization of nonlinear and chaotic temporal functions, toward light-based neural systems.


title: Towards Optimality in Parallel Scheduling

abstract: To keep pace with Moore's law, chip designers have focused on increasing the number of cores per chip rather than single core performance. In turn, modern jobs are often designed to run on any number of cores. However, to effectively leverage these multi-core chips, one must address the question of how many cores to assign to each job. Given that jobs receive sublinear speedups from additional cores, there is an obvious tradeoff: allocating more cores to an individual job reduces the job's runtime, but in turn decreases the efficiency of the overall system. We ask how the system should schedule jobs across cores so as to minimize the mean response time over a stream of incoming jobs. To answer this question, we develop an analytical model of jobs running on a multi-core machine. We prove that EQUI, a policy which continuously divides cores evenly across jobs, is optimal when all jobs follow a single speedup curve and have exponentially distributed sizes. EQUI requires jobs to change their level of parallelization while they run. Since this is not possible for all workloads, we consider a class of "fixed-width" policies, which choose a single level of parallelization, k, to use for all jobs. We prove that, surprisingly, it is possible to achieve EQUI's performance without requiring jobs to change their levels of parallelization by using the optimal fixed level of parallelization, k*. We also show how to analytically derive the optimal k* as a function of the system load, the speedup curve, and the job size distribution. In the case where jobs may follow different speedup curves, finding a good scheduling policy is even more challenging. We find that policies like EQUI which performed well in the case of a single speedup function now perform poorly. We propose a very simple policy, GREEDY*, which performs near-optimally when compared to the numerically-derived optimal policy.


title: Mixmaster dynamics in the dust time gauge

abstract: We study the Hamiltonian dynamics of the dust-Bianchi IX universe in dust time gauge. This model has three physical metric degrees of freedom, with evolution determined by a time-independent physical Hamiltonian. This approach gives a new physical picture where dust-Bianchi IX dynamics is described by oscillations between dust-Kasner solutions, rather than between vacuum-Kasner solutions. We derive a generalized transition law between these phases, which has a matter component. Sufficiently close to a singularity, we show that this law reduces to the vacuum Belinski-Khalatnikov-Lifshitz map. We include an analysis with dust and a scalar field. Lastly, we describe a path integral quantization using the dust-time physical Hamiltonian and derive an effective action for the dust-Kasner model by integrating out the anisotropy degrees of freedom.


title: Transverse spinning of light with globally unique handedness

abstract: Access to the transverse spin of light has unlocked new regimes in topological photonics and optomechanics. To achieve the transverse spin of nonzero longitudinal fields, various platforms that derive transversely confined waves based on focusing, interference, or evanescent waves have been suggested. Nonetheless, because of the transverse confinement inherently accompanying sign reversal of the field derivative, the resulting transverse spin handedness experiences spatial inversion, which leads to a mismatch between the densities of the wavefunction and its spin component and hinders the global observation of the transverse spin. Here, we reveal a globally pure transverse spin in which the wavefunction density signifies the spin distribution, by employing inverse molding of the eigenmode in the spin basis. Starting from the target spin profile, we analytically obtain the potential landscape and then show that the elliptic-hyperbolic transition around the epsilon-near-zero permittivity allows for the global conservation of transverse spin handedness across the topological interface between anisotropic metamaterials. Extending to the non-Hermitian regime, we also develop annihilated transverse spin modes to cover the entire Poincare sphere of the meridional plane. Our results enable the complete transfer of optical energy to transverse spinning motions and realize the classical analogy of 3-dimensional quantum spin states.


title: Overcoming the Fundamental Barrier Thickness Limits of Ferroelectric Tunnel Junctions through BaTiO3/SrTiO3 Composite Barriers

abstract: Ferroelectric tunnel junctions (FTJs) have attracted increasing research interest as a promising candidate for non-volatile memories. Recently, significant enhancements of tunneling electroresistance (TER) have been realized through modifications of electrode materials. However, direct control of the FTJ performance through modifying the tunneling barrier has not been adequately explored. Here, adding a new direction to FTJ research, we fabricated FTJs with BaTiO3 single barriers (SB-FTJs) and BaTiO3/SrTiO3 composite barriers (CB-FTJs), and reported a systematic study of FTJ performances by varying the barrier thicknesses and compositions. For the SB-FTJs, the TER is limited by pronounced leakage current for ultrathin barriers and extremely small tunneling current for thick barriers. For the CB-FTJs, the extra SrTiO3 barrier provides an additional degree of freedom to modulate the barrier potential and tunneling behavior. The resultant high tunability can be utilized to overcome the barrier thickness limits and enhance the overall CB-FTJ performances beyond those of SB-FTJ. Our results reveal a new paradigm to manipulate the FTJs through designing multilayer tunneling barriers with hybrid functionalities.


title: On the Performance of NOMA-Based Cooperative Relaying Systems over Rician Fading Channels

abstract: Non-orthogonal multiple access (NOMA) is a promising technique for the fifth generation (5G) wireless communications. As users with good channel conditions can serve as relays to enhance the system performance by using successive interference cancellation (SIC), the integration of NOMA and cooperative relaying has recently attracted increasing interests. In this paper, a NOMA-based cooperative relaying system is studied, and an analytical framework is developed to evaluate its performance. Specifically, the performance of NOMA over Rician fading channels is studied, and the exact expression of the average achievable rate is derived. Moreover, we also propose an approximation method to calculate the achievable rate by using the Gauss-Chebyshev Integration. Numerical results confirm that our derived analytical results match well with the Monte Carlo simulations.


title: OBJ2TEXT: Generating Visually Descriptive Language from Object Layouts

abstract: Generating captions for images is a task that has recently received considerable attention. In this work we focus on caption generation for abstract scenes, or object layouts where the only information provided is a set of objects and their locations. We propose OBJ2TEXT, a sequence-to-sequence model that encodes a set of objects and their locations as an input sequence using an LSTM network, and decodes this representation using an LSTM language model. We show that our model, despite encoding object layouts as a sequence, can represent spatial relationships between objects, and generate descriptions that are globally coherent and semantically relevant. We test our approach in a task of object-layout captioning by using only object annotations as inputs. We additionally show that our model, combined with a state-of-the-art object detector, improves an image captioning model from 0.863 to 0.950 (CIDEr score) in the test benchmark of the standard MS-COCO Captioning task.


title: PatchShuffle Regularization

abstract: This paper focuses on regularizing the training of the convolutional neural network (CNN). We propose a new regularization approach named ``PatchShuffle`` that can be adopted in any classification-oriented CNN models. It is easy to implement: in each mini-batch, images or feature maps are randomly chosen to undergo a transformation such that pixels within each local patch are shuffled. Through generating images and feature maps with interior orderless patches, PatchShuffle creates rich local variations, reduces the risk of network overfitting, and can be viewed as a beneficial supplement to various kinds of training regularization techniques, such as weight decay, model ensemble and dropout. Experiments on four representative classification datasets show that PatchShuffle improves the generalization ability of CNN especially when the data is scarce. Moreover, we empirically illustrate that CNN models trained with PatchShuffle are more robust to noise and local changes in an image.


title: The homotopy theory of coalgebras over simplicial comonads

abstract: We apply the Acyclicity Theorem of Hess, Kerdziorek, Riehl, and Shipley (recently corrected by Garner, Kedziorek, and Riehl) to establishing the existence of model category structure on categories of coalgebras over comonads arising from simplicial adjunctions, under mild conditions on the adjunction and the associated comonad. We study three concrete examples of such adjunctions where the left adjoint is comonadic and show that in each case the component of the derived counit of the comparison adjunction at any fibrant object is an isomorphism, while the component of the derived unit at any 1-connected object is a weak equivalence. To prove this last result, we explain how to construct explicit fibrant replacements for 1-connected coalgebras in the image of the canonical comparison functor from the Postnikov decompositions of their underlying simplicial sets. We also show in one case that the derived unit is precisely the Bousfield-Kan completion map.


title: Fast and Robust Determination of Power System Emergency Control Actions

abstract: This paper outlines an optimization framework for choosing fast and reliable control actions in a transmission grid emergency situation. We consider contractual load shedding and generation re-dispatch as exemplary emergency actions. To achieve computational efficiency and scalability, this novel formulation of the robust corrective action optimization is stated in terms of voltages and currents, as opposed to standard formulation in terms of power flows and voltages. The current-voltage formulation is natural for expressing voltage and thermal transmission line constraints, and its major advantage is in the linearity of the power flow equations. The complexity of the current-voltage formulation, which is mainly related to the transformation of voltages and currents into active and reactive power injections, can be largely avoided by stating the cost function directly in terms of currents and using linearized or robust estimates for load and generation constraints. The paper considers five different optimization problem formulations, including the full current-voltage formulation, as well as two convex and two linear approximation to the problem. In a case study on an illustrative case study for the IEEE RTS96 system, we demonstrate pros and cons of the different formulations based on simulations for both normal and contingency operation.


title: Intrapixel effects of CCD and CMOS detectors

abstract: Intrapixel nonuniformity is known to exist in CCD and CMOS image sensors, though the effects in backside illuminated (BSI) CCDs are too small to be a concern for most astronomical observations. However, projects like the Large Synoptic Survey Telescope require precise knowledge of the detector characteristics, and intrapixel effects may need more attention. By scanning CCD and CMOS cameras with a small light spot (unresolved by the optics), we find in the images that the spot's flux, centroid displacement, and ellipticity vary periodically on the pixel scale in most cases. The amplitude of variation depends on not only the detector but also how well the spot is sampled by the pixels. With a spot radius of 2 pixels (encircling 80% energy) as measured, the flux and the ellipticity extracted from the BSI CCD camera vary by 0.2-0.3% (rms) and 0.005 (rms), respectively, while the deviation of the centroid position (rms ~ 0.01 pixel) is not correlated with the pixels. The effects are more pronounced for the BSI CMOS camera and even worse for the frontside illuminated CMOS camera. The results suggest that a closer examination of the intrapixel effects is needed for precision astronomy.


title: Stability of the quantum Sherrington-Kirkpatrick spin glass model

abstract: We study in detail the quantum Sherrington-Kirkpatrick (SK) model, i.e. the infinite-range Ising spin glass in a transverse field, by solving numerically the effective one-dimensional model that the quantum SK model can be mapped to in the thermodynamic limit. We find that the replica symmetric (RS) solution is unstable down to zero temperature, in contrast to some previous claims, and so there is not only a line of transitions in the (longitudinal) field-temperature plane (the de Almeida-Thouless, AT, line) where replica symmetry is broken, but also a quantum de Almeida-Thouless (QuAT) line in the transverse field-longitudinal field plane at $T = 0$. If the QuAT line also occurs in models with short-range interactions its presence might affect the performance of quantum annealers when solving spin glass-type problems with a bias (i.e. magnetic field).


title: Comparing 4-Manifolds in the Pants Complex via Trisections

abstract: Given two smooth, oriented, closed 4-manifolds $M_1$ and $M_2$, we construct two invariants, $D^P(M_1, M_2)$ and $D(M_1, M_2)$, coming from distances in the pants complex and the dual curve complex respectively. To do this, we adapt work of Johnson on Heegaard splittings of 3-manifolds to the trisections of 4-manifolds introduced by Gay and Kirby. Our main results are that the invariants are independent of the choices made throughout the process, as well as interpretations of "nearby" manifolds. This naturally leads to various graphs of 4-manifolds coming from unbalanced trisections, and we briefly explore their properties.


title: Blow-up of solutions for weakly coupled systems of complex Ginzburg-Landau equations

abstract: Blow-up phenomena ofvweakly coupled systems of several evolution equations, especially complex Ginzburg-Landau equationsvis shown by a straightforward ODE approach not so-called test-function method, which gives the natural blow-up rate. The difficulty of the proof is that, unlike the single case, terms which come from the fact that the Laplacian cannot be absorbed into the weakly coupled nonlinearities. A similar ODE approach is applied to heat systems by Mochizuki to obtain the lower estimate of lifespan.


title: The Amplification in Strong Electric and Magnetic Fields

abstract: We consider the amplification in strong electric (magnetic) fields that are uniform along the direction of the electron motion and are not uniform in the transverse direction. It is shown that in such a system the gain is increased compared with that of a magnetic undulator by a value of the order of gamma ( gamma is the relativistic factor). Numerical estimates are presented that determine the possibility of amplification of a test wave in the optical and near infrared bands.


title: Uniqueness for bubbling solutions with collapsing singularities

abstract: The seminal work \cite{bm} by Brezis and Merle showed that the bubbling solutions of the mean field equation have the property of mass concentration. Recently, Lin and Tarantello in \cite{lt} found that the "bubbling implies mass concentration" phenomena might not hold if there is a collapse of singularities. Furthermore, a sharp estimate \cite{llty} for the bubbling solutions has been obtained. In this paper, we prove that there exists at most one sequence of bubbling solutions if the collapsing singularity occurs. The main difficulty comes from that after re-scaling, the difference of two solutions locally converges to an element in the kernel space of the linearized operator. It is well-known that the kernel space is three dimensional. So the main technical ingredient of the proof is to show that the limit after re-scaling is orthogonal to the kernel space.


title: Switching and Data Injection Attacks on Stochastic Cyber-Physical Systems: Modeling, Resilient Estimation and Attack Mitigation

abstract: In this paper, we consider the problem of attack-resilient state estimation, that is to reliably estimate the true system states despite two classes of attacks: (i) attacks on the switching mechanisms and (ii) false data injection attacks on actuator and sensor signals, in the presence of unbounded stochastic process and measurement noise signals. We model the systems under attack as hidden mode stochastic switched linear systems with unknown inputs and propose the use of a multiple-model inference algorithm to tackle these security issues. Moreover, we characterize fundamental limitations to resilient estimation (e.g., upper bound on the number of tolerable signal attacks) and discuss the topics of attack detection, identification and mitigation under this framework. Simulation examples of switching and false data injection attacks on a benchmark system and an IEEE 68-bus test system show the efficacy of our approach to recover resilient (i.e., asymptotically unbiased) state estimates as well as to identify and mitigate the attacks.


title: Adversarial Variational Optimization of Non-Differentiable Simulators

abstract: Complex computer simulators are increasingly used across fields of science as generative models tying parameters of an underlying theory to experimental observations. Inference in this setup is often difficult, as simulators rarely admit a tractable density or likelihood function. We introduce Adversarial Variational Optimization (AVO), a likelihood-free inference algorithm for fitting a non-differentiable generative model incorporating ideas from generative adversarial networks, variational optimization and empirical Bayes. We adapt the training procedure of generative adversarial networks by replacing the differentiable generative network with a domain-specific simulator. We solve the resulting non-differentiable minimax problem by minimizing variational upper bounds of the two adversarial objectives. Effectively, the procedure results in learning a proposal distribution over simulator parameters, such that the JS divergence between the marginal distribution of the synthetic data and the empirical distribution of observed data is minimized. We evaluate and compare the method with simulators producing both discrete and continuous data.


title: Gaussian and plane-wave mixed density fitting for periodic systems

abstract: We introduce a mixed density fitting scheme that uses both a Gaussian and a plane-wave fitting basis to accurately evaluate electron repulsion integrals in crystalline systems. We use this scheme to enable efficient all-electron Gaussian based periodic density functional and Hartree-Fock calculations.


title: On invariant Riemannian metrics on Ledger-Obata spaces

abstract: We study invariant metrics on Ledger-Obata spaces $F^m/\mathrm{diag}(F)$. We give the classification and an explicit construction of all naturally reductive metrics, and also show that in the case $m=3$, any invariant metric is naturally reductive. We prove that a Ledger-Obata space is a geodesic orbit space if and only if the metric is naturally reductive. We then show that a Ledger-Obata space is reducible if and only if it is isometric to the product of Ledger-Obata spaces (and give an effective method of recognising reducible metrics), and that the full connected isometry group of an irreducible Ledger-Obata space $F^m/\mathrm{diag}(F)$ is $F^m$. We deduce that a Ledger-Obata space is a geodesic orbit manifold if and only if it is the product of naturally reductive Ledger-Obata spaces.


title: Risk-Averse Joint Capacity Evaluation of PV Generation and Electric Vehicle Charging Stations in Distribution Networks

abstract: Increasing penetration of distribution generation (DG) and electric vehicles (EVs) calls for an effective way to estimate the achievable capacity connected to the distribution systems, but the exogenous uncertainties of DG outputs and EV charging loads make it challengeable. This study provides a joint capacity evaluation method with a risk threshold setting function for photovoltaic (PV) generation and EV charging stations (EVCSs). The method is mathematically formulated as a distributionally robust joint chance constrained programming model. And the worst-case conditional value at risk (WC-CVaR) approximation and an iterative algorithm based on semidefinite program (SDP) are used to solve the model. Finally, the method test is carried out numerically on IEEE 33-bus radial distribution system.


title: Data-Driven Planning of Plug-in Hybrid Electric Taxi Charging Stations in Urban Environments: A Case in the Central Area of Beijing

abstract: Plug-in electric vehicles (PEVs) can contribute to energy and environmental challenges. Among different types of PEVs, plug-in hybrid electric taxis (PHETs) go in advance. In this study, we provide a spatial and temporal PHET charging demand forecasting method based on one-month global positioning system (GPS)-based taxi travel data in Beijing. Then, using the charging demand forecasting results, a mixed integer linear programming (MILP) model is formulated to plan PHET charging stations in the central area of Beijing. The model minimizes both investment and operation costs of all the PHET charging stations and takes into account the service radius of charging stations, charging demand satisfaction and rational occupation rates of chargers. At last, the test of the planning method is carried out numerically through simulations and the analysis is complemented according to the results.


title: Note on soft theorems and memories in even dimensions

abstract: Recently, it has been shown that the Weinberg's formula for soft graviton production is essentially a Fourier transformation of the formula for gravitational memory which provides an effective way to understand how the classical calculation arises as a limiting case of the quantum result. In this note, we propose a general framework that connects the soft theorems to the radiation fields obtained from classical computation for different theories in even dimensions. We show that the latter is nothing but Fourier transformation of the former. The memory formulas can be derived from radiation fields explicitly.


title: Deep Networks for Compressed Image Sensing

abstract: The compressed sensing (CS) theory has been successfully applied to image compression in the past few years as most image signals are sparse in a certain domain. Several CS reconstruction models have been recently proposed and obtained superior performance. However, there still exist two important challenges within the CS theory. The first one is how to design a sampling mechanism to achieve an optimal sampling efficiency, and the second one is how to perform the reconstruction to get the highest quality to achieve an optimal signal recovery. In this paper, we try to deal with these two problems with a deep network. First of all, we train a sampling matrix via the network training instead of using a traditional manually designed one, which is much appropriate for our deep network based reconstruct process. Then, we propose a deep network to recover the image, which imitates traditional compressed sensing reconstruction processes. Experimental results demonstrate that our deep networks based CS reconstruction method offers a very significant quality improvement compared against state of the art ones.


title: Robust Large Gap Quantum Spin Hall Insulators in Methyl-functionalized III-Bi Buckled Honeycombs

abstract: A large bulk band gap is critical for the applications of quantum spin hall (QSH) insulators in spintronics at room temperature. Based on first-principles calculations, we predict that the methyl-functionalized III-Bi monolayers, namely III-Bi-(CH3)2 (III=Ga, In, Tl) thin films, own QSH states with band gap as large as 0.260, 0.304 and 0.843 eV, respectively, making them suitable for room-temperature applications. The topological characteristics are confirmed by s-px,y band inversion, topological invariant Z2, and the topologically protected edge states. Noticeably, for GaBi/InBi-(CH3)2 films, the s-px,y band inversion occurred in the progress of spin-orbital coupling (SOC), while for TlBi(CH3)2 film, the s-px,y band inversion happened in the progress of chemical bonding. Significantly, the QSH states in III-Bi-(CH3)2 films are robust against the mechanical strains and various methyl coverages, making these films particularly flexible to substrate choice for device applications. Besides, the h-BN substrate is an ideal substrate for III-Bi-(CH3)2 films to realize large gap nontrivial topological states.. These findings demonstrate that the methyl-functionalized III-Bi films may be good QSH effect platforms for topological electronic devices design and fabrication in spintronics.


title: Quantitative $C^1$-estimates by Bismut formulae

abstract: For a $C^2$ function $u$ and an elliptic operator $L$, we prove a quantitative estimate for the derivative $du$ in terms of local bounds on $u$ and $Lu$. An integral version of this estimate is then used to derive a condition for the zero-mean value property of $\Delta u$. An extension to differential forms is also given. Our approach is probabilistic and could easily be adapted to other settings.


title: Enhanced and tunable electric dipole-dipole interactions near a planar metal film

abstract: We investigate the enhanced electric dipole-dipole interaction by surface plasmon polaritons (SPPs) supported by the planar metal film waveguide. By taking two nitrogen-vacancy (NV) center electric-dipoles in diamond as an example, both the coupling strength and collective relaxation of two dipoles are studied with numerical Green Function method. Compared to the two-dipole coupling on planar surface, metal film provides stronger and tunable coupling coefficients. Enhancement of the interaction between coupled NV center dipoles could have applications in both quantum information and energy transfer investigation. Our investigation provides a systematical result for experimental applications based on dipole-dipole interaction mediated with SPPs on planar metal film.


title: The domination number and the least Q-eigenvalue II

abstract: Denote by $L_{g, l}$ the $lollipop$ $graph$ obtained by attaching a pendant path $\mathbb{P}=v_{g}v_{g+1}\cdots v_{g+l}$ ($l\geq 1$) to a cycle $\mathbb{C}=v_{1}v_{2}\cdots v_{g}v_{1}$ ($g\geq 3$). A $\mathcal {F}_{g, l}$-$graph$ of order $n\geq g+1$ is defined to be the graph obtained by attaching $n-g-l$ pendent vertices to some of the nonpendant vertices of $L_{g, l}$ in which each vertex other than $v_{g+l-1}$ is attached at most one pendant vertex. A $\mathcal {F}^{\circ}_{g, l}$-graph is a $\mathcal {F}_{g, l}$-$graph$ in which $v_{g}$ is attached with pendant vertex. Denote by $q_{min}$ the $least$ $Q$-$eigenvalue$ of a graph. In this paper, we proceed on considering the domination number, the least $Q$-eigenvalue of a graph as well as their relation. Further results obtained are as follows: $\mathrm{(i)}$ some results about the changing of the domination number under the structural perturbation of a graph are represented; $\mathrm{(ii)}$ among all nonbipartite unicyclic graphs of order $n$, with both domination number $\gamma$ and girth $g$ ($g\leq n-1$), the minimum $q_{min}$ attains at a $\mathcal {F}_{g, l}$-graph for some $l$; $\mathrm{(iii)}$ among the nonbipartite graphs of order $n$ and with given domination number which contain a $\mathcal {F}^{\circ}_{g, l}$-graph as a subgraph, some lower bounds for $q_{min}$ are represented; $\mathrm{(iv)}$ among the nonbipartite graphs of order $n$ and with given domination number $\frac{n}{2}$, $\frac{n-1}{2}$, the minimum $q_{min}$ is completely determined respectively; $\mathrm{(v)}$ among the nonbipartite graphs of order $n\geq 4$, and with both domination number $\frac{n+1}{3}<\gamma\leq \frac{n}{2}$ and odd-girth (the length of the shortest odd cycle) at most $5$, the minimum $q_{min}$ is completely determined.


title: A signature-based machine learning model for bipolar disorder and borderline personality disorder

abstract: Mobile technologies offer opportunities for higher resolution monitoring of health conditions. This opportunity seems of particular promise in psychiatry where diagnoses often rely on retrospective and subjective recall of mood states. However, getting actionable information from these rather complex time series is challenging, and at present the implications for clinical care are largely hypothetical. This research demonstrates that, with well chosen cohorts (of bipolar disorder, borderline personality disorder, and control) and modern methods, it is possible to objectively learn to identify distinctive behaviour over short periods (20 reports) that effectively separate the cohorts. Participants with bipolar disorder or borderline personality disorder and healthy volunteers completed daily mood ratings using a bespoke smartphone app for up to a year. A signature-based machine learning model was used to classify participants on the basis of the interrelationship between the different mood items assessed and to predict subsequent mood. The signature methodology was significantly superior to earlier statistical approaches applied to this data in distinguishing the participant three groups, clearly placing 75% into their original groups on the basis of their reports. Subsequent mood ratings were correctly predicted with greater than 70% accuracy in all groups. Prediction of mood was most accurate in healthy volunteers (89-98%) compared to bipolar disorder (82-90%) and borderline personality disorder (70-78%).


title: Shadow of Schwarzschild-Tangherlini black holes

abstract: We study the shadow cast by the H$D$ Schwarzschild-Tangherlini black hole, and analytically calculate the influence of extra dimensions on the shadow of a black hole. A black hole casts a shadow as an optical appearance because of its strong gravitational field which is known to be a dark zone covered by a circle for a Schwarzschild black hole. We demonstrate that the null geodesic equation can be integrated by Hamilton-Jacobi approach, which enables us to investigate the shadow cast by the H$D$ Schwarzschild-Tangherlini black holes. Interestingly, it turns out that, for fixed values of the mass parameter, the shadow in H$D$ spacetimes are smaller when compared with 4$D$ Schwarzschild black hole. Further, the shadows of H$D$ Schwarzschild-Tangherlini black holes are concentric circles with a radius of the circle decreases with increase in $D$. We visualize the photon regions and the shadows in various dimensions for different values of the parameters, and the energy emission rates are is also investigated. Our results, in the limit $D=4$, reduced exactly to \emph{vis-$\grave{a}$-vis} Schwarzschild black hole case.


title: Construction and Response of a Highly Granular Scintillator-based Electromagnetic Calorimeter

abstract: A highly granular electromagnetic calorimeter with scintillator strip readout is being developed for future lepton collider experiments. A prototype of 21.5 $X_0$ depth and $180 \times 180 $mm$^2$ transverse dimensions was constructed, consisting of 2160 individually read out $10 \times 45 \times 3$ mm$^3$ scintillator strips. This prototype was tested using electrons of 2--32 GeV at the Fermilab Test Beam Facility in 2009. Deviations from linear energy response were less than 1.1\%, and the intrinsic energy resolution was determined to be $(12.5 \pm 0.1 (\mathrm{stat.}) \pm0.4 (\mathrm{syst.}))\%/\sqrt{E[\mathrm{GeV}]}\oplus (1.2 \pm 0.1(\mathrm{stat.})^{+0.6}_{-0.7}(\mathrm{syst.}))\%$, where the uncertainties correspond to statistical and systematic sources, respectively.


title: Partition-based discrete-time quantum walks

abstract: We introduce a family of discrete-time quantum walks, called two-partition model, based on two equivalence-class partitions of the computational basis, which establish the notion of local dynamics. This family encompasses most versions of unitary discrete-time quantum walks driven by two local operators studied in literature, such as the coined model, Szegedy's model, and the 2-tessellable staggered model. We also analyze the connection of those models with the two-step coined model, which is driven by the square of the evolution operator of the standard discrete-time coined walk. We prove formally that the two-step coined model, an extension of Szegedy model for multigraphs, and the two-tessellable staggered model are unitarily equivalent. Then, selecting one specific model among those families is a matter of taste not generality.


title: Single Image Super-Resolution with Dilated Convolution based Multi-Scale Information Learning Inception Module

abstract: Traditional works have shown that patches in a natural image tend to redundantly recur many times inside the image, both within the same scale, as well as across different scales. Make full use of these multi-scale information can improve the image restoration performance. However, the current proposed deep learning based restoration methods do not take the multi-scale information into account. In this paper, we propose a dilated convolution based inception module to learn multi-scale information and design a deep network for single image super-resolution. Different dilated convolution learns different scale feature, then the inception module concatenates all these features to fuse multi-scale information. In order to increase the reception field of our network to catch more contextual information, we cascade multiple inception modules to constitute a deep network to conduct single image super-resolution. With the novel dilated convolution based inception module, the proposed end-to-end single image super-resolution network can take advantage of multi-scale information to improve image super-resolution performance. Experimental results show that our proposed method outperforms many state-of-the-art single image super-resolution methods.


title: Predicting the Gender of Indonesian Names

abstract: We investigated a way to predict the gender of a name using character-level Long-Short Term Memory (char-LSTM). We compared our method with some conventional machine learning methods, namely Naive Bayes, logistic regression, and XGBoost with n-grams as the features. We evaluated the models on a dataset consisting of the names of Indonesian people. It is not common to use a family name as the surname in Indonesian culture, except in some ethnicities. Therefore, we inferred the gender from both full names and first names. The results show that we can achieve 92.25% accuracy from full names, while using first names only yields 90.65% accuracy. These results are better than the ones from applying the classical machine learning algorithms to n-grams.


title: On locally compact shift-continuous topologies on the $\alpha$-bicyclic monoid

abstract: A topology $\tau$ on a monoid $S$ is called {\em shift-continuous} if for every $a,b\in S$ the two-sided shift $S\to S$, $x\mapsto axb$, is continuous. For every ordinal $\alpha\le \omega$, we describe all shift-continuous locally compact Hausdorff topologies on the $\alpha$-bicyclic monoid $\mathcal{B}_{\alpha}$. More precisely, we prove that the lattice of shift-continuous locally compact Hausdorff topologies on $\mathcal{B}_{\alpha}$ is anti-isomorphic to the segment of $[1,\alpha]$ of ordinals, endowed with the natural well-order. Also we prove that for each ordinal $\alpha$ the $\alpha+1$-bicyclic monoid $\mathcal{B}_{\alpha+1}$ is isomorphic to the Bruck extension of the $\alpha$-bicyclic monoid $\mathcal{B}_{\alpha}$.


title: Scaling relations in large-Prandtl-number natural thermal convection

abstract: In this study we follow Grossmann and Lohse, Phys. Rev. Lett. 86 (2001), who derived various scalings regimes for the dependence of the Nusselt number $Nu$ and the Reynolds number $Re$ on the Rayleigh number $Ra$ and the Prandtl number $Pr$. We focus on theoretical arguments as well as on numerical simulations for the case of large-$Pr$ natural thermal convection. Based on an analysis of self-similarity of the boundary layer equations, we derive that in this case the limiting large-$Pr$ boundary-layer dominated regime is I$_\infty^<$, introduced and defined in [1], with the scaling relations $Nu\sim Pr^0\,Ra^{1/3}$ and $Re\sim Pr^{-1}\,Ra^{2/3}$. Our direct numerical simulations for $Ra$ from $10^4$ to $10^9$ and $Pr$ from 0.1 to 200 show that the regime I$_\infty^<$ is almost indistinguishable from the regime III$_\infty$, where the kinetic dissipation is bulk-dominated. With increasing $Ra$, the scaling relations undergo a transition to those in IV$_u$ of reference [1], where the thermal dissipation is determined by its bulk contribution.


title: Mean curvature flow solitons in the presence of conformal vector fields

abstract: The aim of this paper is to introduce a notion of mean curvature flow soliton general enough to encompass target spaces of constant sectional curvature, Riemannian products or, in increasing generality, warped product spaces.


title: Galois structure of the holomorphic differentials of curves

abstract: Let $X$ be a smooth projective geometrically irreducible curve over a perfect field $k$ of positive characteristic $p$. Suppose $G$ is a finite group acting faithfully on $X$ such that $G$ has non-trivial cyclic Sylow $p$-subgroups. We show that the decomposition of the space of holomorphic differentials of $X$ into a direct sum of indecomposable $k[G]$-modules is uniquely determined by the lower ramification groups and the fundamental characters of closed points of $X$ that are ramified in the cover $X\to X/G$. We apply our method to determine the $\mathrm{PSL}(2,\mathbb{F}_\ell)$-module structure of the space of holomorphic differentials of the reduction of the modular curve $\mathcal{X}(\ell)$ modulo $p$ when $p$ and $\ell$ are distinct odd primes and the action of $\mathrm{PSL}(2,\mathbb{F}_\ell)$ on this reduction is not tamely ramified. This provides some non-trivial congruences modulo appropriate maximal ideals containing $p$ between modular forms arising from isotypic components with respect to the action of $\mathrm{PSL}(2,\mathbb{F}_\ell)$ on $\mathcal{X}(\ell)$.


title: Active Galactic Nuclei: what's in a name?

abstract: Active Galactic Nuclei (AGN) are energetic astrophysical sources powered by accretion onto supermassive black holes in galaxies, and present unique observational signatures that cover the full electromagnetic spectrum over more than twenty orders of magnitude in frequency. The rich phenomenology of AGN has resulted in a large number of different "flavours" in the literature that now comprise a complex and confusing AGN "zoo". It is increasingly clear that these classifications are only partially related to intrinsic differences between AGN, and primarily reflect variations in a relatively small number of astrophysical parameters as well the method by which each class of AGN is selected. Taken together, observations in different electromagnetic bands as well as variations over time provide complementary windows on the physics of different sub-structures in the AGN. In this review, we present an overview of AGN multi-wavelength properties with the aim of painting their "big picture" through observations in each electromagnetic band from radio to gamma-rays as well as AGN variability. We address what we can learn from each observational method, the impact of selection effects, the physics behind the emission at each wavelength, and the potential for future studies. To conclude we use these observations to piece together the basic architecture of AGN, discuss our current understanding of unification models, and highlight some open questions that present opportunities for future observational and theoretical progress.


title: How ants move: individual and collective scaling properties

abstract: The motion of social insects is often used a paradigmatic example of complex adaptive dynamics arising from decentralized individual behavior. In this paper we revisit the topic of the ruling laws behind burst of activity in ants. The analysis, done over previously reported data, reconsider the causation arrows, proposed at individual level, not finding any link between the duration of the ants' activity and its moving speed. Secondly, synthetic trajectories created from steps of different ants, demonstrate that a Markov process can explain the previously reported speed shape profile. Finally we show that as more ants enter the nest, the faster they move, which implies a collective property. Overall these results provides a mechanistic explanation for the reported behavioral laws, and suggest a formal way to further study the collective properties in these scenarios.


title: Contravariant form for reduction algebras and Pieri rule

abstract: We study properties and constructions of contravariant forms on reduction algebras. As an application we compute norms of highest weight vectors in the tensor product of an irreducible finite dimensional representation of the Lie algebra gl(n) with a symmetric or wedge tensor power of its fundamental representation. Their zeroes describe Pieri rules.


title: AutOMP: An Automatic OpenMP Parallelization Generator for Variable-Oriented High-Performance Scientific Codes

abstract: OpenMP is a cross-platform API that extends C, C++ and Fortran and provides shared-memory parallelism platform for those languages. The use of many cores and HPC technologies for scientific computing has been spread since the 1990s, and now takes part in many fields of research. The relative ease of implementing OpenMP, along with the development of multi-core shared memory processors (such as Intel Xeon Phi) makes OpenMP a favorable method for parallelization in the process of modernizing a legacy codes. Legacy scientific codes are usually holding large number of physical arrays which being used and updated by the code routines. In most of the cases the parallelization of such code focuses on loop parallelization. A key step in this parallelization is deciding which of the variables in the parallelized scope should be private (so each thread will hold a copy of them), and which variables should be shared across the threads. Other important step is finding which variables should be synchronized after the loop execution. In this work we present an automatic pre-processor that preforms these stages - AutOMP (Automatic OpenMP). AutOMP recognize all the variables assignments inside a loop. These variables will be private unless the assignment is of an array element which depend on the loop index variable. Afterwards, AutOMP finds the places where threads synchronization is needed, and which reduction operator is to be used. At last, the program provides the parallelization command to be used for parallelizing the loop.


title: On the use of the autonomous Birkhoff equations in Lie series perturbation theory

abstract: In this article, we present the Lie transformation algorithm for autonomous Birkhoff systems. Here, we are referring to Hamiltonian systems that obey a symplectic structure of the general form. Two examples of normalization in the restricted three-body problem are given to illustrate the application of the algorithm in perturbation theory. The efficiency of this algorithm for problems of asymptotic integration in dynamics is discussed for the case where there is a need to use non-canonical variables in the phase space.


title: Clinical Patient Tracking in the Presence of Transient and Permanent Occlusions via Geodesic Feature

abstract: This paper develops a method to use RGB-D cameras to track the motions of a human spinal cord injury patient undergoing spinal stimulation and physical rehabilitation. Because clinicians must remain close to the patient during training sessions, the patient is usually under permanent and transient occlusions due to the training equipment and the movements of the attending clinicians. These occlusions can significantly degrade the accuracy of existing human tracking methods. To improve the data association problem in these circumstances, we present a new global feature based on the geodesic distances of surface mesh points to a set of anchor points. Transient occlusions are handled via a multi-hypothesis tracking framework. To evaluate the method, we simulated different occlusion sizes on a data set captured from a human in varying movement patterns, and compared the proposed feature with other tracking methods. The results show that the proposed method achieves robustness to both surface deformations and transient occlusions.


title: Descriptor System Tools (DSTOOLS) User's Guide

abstract: The Descriptor System Tools (DSTOOLS) is a collection of MATLAB functions for the operation on and manipulation of rational transfer function matrices via their descriptor system realizations. The DSTOOLS collection relies on the Control System Toolbox and several mex-functions based on the Systems and Control Library SLICOT. Many of the implemented functions are based on the computational procedures described in Chapter 10 of the book: "A. Varga, Solving Fault Diagnosis Problems - Linear Synthesis Techniques, Springer, 2017". This document is the User's Guide for the version V0.71 of DSTOOLS. First, we present the mathematical background on rational matrices and descriptor systems. Then, we give in-depth information on the command syntax of the main computational functions. Several examples illustrate the use of the main functions of DSTOOLS.


title: Sur l'op\'erateur $\bar{\partial}$ et les fonctions diff\'erentiables au sens de Whitney dans un domaine $q$-convexe de $\mathbb{C}^n$

abstract: We solve the $\bar{\partial}$-problem for differential forms in the sens of Whitney.


title: Berry Phase and Topological Effects of Phonons

abstract: Phonons as collective excitations of lattice vibrations are the main heat carriers in solids. Tremendous effort has been devoted to investigate phonons and related properties, giving rise to an intriguing field of phononics, which is of great importance to many practical applications, including heat dissipation, thermal barrier coating, thermoelectrics and thermal control devices. Meanwhile, the research of topology-related physics, awarded the 2016 Nobel Prize in Physics, has led to discoveries of various exotic quantum states of matter, including the quantum (anomalous/spin) Hall [Q(A/S)H] effects, topological insulators/semimetals and topological superconductors. An emerging research field is to bring topological concepts for a new paradigm phononics---"topological phononics". In this Perspective, we will briefly introduce this emerging field and discuss the use of novel quantum degrees of freedom like the Berry phase and topology for manipulating phonons in unprecedentedly new ways.


title: A relation on the effective conductivity of composites

abstract: Consider a 2D composites with non-overlapping equal inclusions imbedded in a host material of the normalized unit conductivity. The conductivity of inclusions takes two values $\sigma_1$ and $\sigma_2$ with the probabilities $p$ and $1-p$, respectively. We prove that the effective conductivity tensor of the considered three-phase random composite is equal to the effective conductivity tensor of the two-phase deterministic composite with the same inclusions of the conductivity $\sigma=[p(\sigma_1~-~\sigma_2)+~\sigma_2+\sigma_1\sigma_2] [1+\sigma_1-p(\sigma_1-\sigma_2)]^{-1}$.


title: $Q\bar Q$ ($Q\in \{b, c\}$) spectroscopy using Cornell potential

abstract: The mass spectra and decay properties of heavy quarkonia are computed in nonrelativistic quark-antiquark Cornell potential model. We have employed the numerical solution of Schr\"odinger equation to obtain their mass spectra using only four parameters namely quark mass ($m_c$, $m_b$) and confinement strength ($A_{c\bar c}$, $A_{b\bar b}$). The spin hyperfine, spin-orbit and tensor components of the one gluon exchange interaction are computed perturbatively to determine the mass spectra of excited $S$, $P$, $D$ and $F$ states. Digamma, digluon and dilepton decays of these mesons are computed using the model parameters and numerical wave functions. The predicted spectroscopy and decay properties for quarkonia are found to be consistent with available experimental observations and results from other theoretical models. We also compute mass spectra and life time of the $B_c$ meson without additional parameters. The computed electromagnetic transition widths of heavy quarkonia and $B_c$ mesons are in tune with available experimental data and other theoretical approaches.


title: Reduction of a metapopulation genetic model to an effective one island model

abstract: We explore a model of metapopulation genetics which is based on a more ecologically motivated approach than is frequently used in population genetics. The size of the population is regulated by competition between individuals, rather than by artificially imposing a fixed population size. The increased complexity of the model is managed by employing techniques often used in the physical sciences, namely exploiting time-scale separation to eliminate fast variables and then constructing an effective model from the slow modes. Remarkably, an initial model with 2$\mathcal{D}$ variables, where $\mathcal{D}$ is the number of islands in the metapopulation, can be reduced to a model with a single variable. We analyze this effective model and show that the predictions for the probability of fixation of the alleles and the mean time to fixation agree well with those found from numerical simulations of the original model.


title: Structural Properties of Uncoded Placement Optimization for Coded Delivery

abstract: A centralized coded caching scheme has been proposed by Maddah-Ali and Niesen to reduce the worst-case load of a network consisting of a server with access to N files and connected through a shared link to K users, each equipped with a cache of size M. However, this centralized coded caching scheme is not able to take advantage of a non-uniform, possibly very skewed, file popularity distribution. In this work, we consider the same network setting but aim to reduce the average load under an arbitrary (known) file popularity distribution. First, we consider a class of centralized coded caching schemes utilizing general uncoded placement and a specific coded delivery strategy, which are specified by a general file partition parameter. Then, we formulate the coded caching design optimization problem over the considered class of schemes with 2^K2^N variables to minimize the average load by optimizing the file partition parameter under an arbitrary file popularity. Furthermore, we show that the optimization problem is convex, and the resulting optimal solution generally improves upon known schemes. Next, we analyze structural properties of the optimization problem to obtain design insights and reduce the complexity. Specifically, we obtain an equivalent linear optimization problem with (K+1)N variables under an arbitrary file popularity and an equivalent linear optimization problem with K+1 variables under the uniform file popularity. Under the uniform file popularity, we also obtain the closed form optimal solution, which corresponds to Maddah-Ali-Niesen's centralized coded caching scheme. Finally, we present an information-theoretic converse bound on the average load under an arbitrary file popularity.


title: Reverses of operator Acz\'el inequality

abstract: In this paper we present some inequalities involving operator decreasing functions and operator means. These inequalities provide some reverses of operator Acz\'el inequality dealing with the weighted geometric mean.


title: Planet formation and disk-planet interactions

abstract: This review is based on lectures given at the 45th Saas-Fee Advanced Course 'From Protoplanetary Disks to Planet Formation' held in March 2015 in Les Diablerets, Switzerland. Starting with an overview of the main characterictics of the Solar System and extrasolar planets, we describe the planet formation process in terms of the sequential accretion scenario. First the growth processes of dust particles to planetesimals and subsequently to terrestrial planets or planetary cores are presented. This is followed by the formation process of the giant planets either by core accretion or gravitational instability. Finally, the dynamical evolution of the orbital elements as driven by disk-planet interaction and the overall evolution of multi-object systems is presented.


title: Fitting Prediction Rule Ensembles with R Package pre

abstract: Prediction rule ensembles (PREs) are sparse collections of rules, offering highly interpretable regression and classification models. This paper presents the R package pre, which derives PREs through the methodology of Friedman and Popescu (2008). The implementation and functionality of package pre is described and illustrated through application on a dataset on the prediction of depression. Furthermore, accuracy and sparsity of PREs is compared with that of single trees, random forest and lasso regression in four benchmark datasets. Results indicate that pre derives ensembles with predictive accuracy comparable to that of random forests, while using a smaller number of variables for prediction.