title: Prediction of Supernova Rates in Known Galaxy-galaxy Strong-lens Systems

abstract: We propose a new strategy of finding strongly-lensed supernovae (SNe) by monitoring known galaxy-scale strong-lens systems. Strongly lensed SNe are potentially powerful tools for the study of cosmology, galaxy evolution, and stellar populations, but they are extremely rare. By targeting known strongly lensed starforming galaxies, our strategy significantly boosts the detection efficiency for lensed SNe compared to a blind search. As a reference sample, we compile the 128 galaxy-galaxy strong-lens systems from the Sloan Lens ACS Survey (SLACS), the SLACS for the Masses Survey, and the Baryon Oscillation Spectroscopic Survey Emission-Line Lens Survey. Within this sample, we estimate the rates of strongly-lensed Type Ia SN (SNIa) and core-collapse SN (CCSN) to be $1.23 \pm 0.12$ and $10.4 \pm 1.1$ events per year, respectively. The lensed SN images are expected to be widely separated with a median separation of 2 arcsec. Assuming a conservative fiducial lensing magnification factor of 5 for the most highly magnified SN image, we forecast that a monitoring program with a single-visit depth of 24.7 mag (5$\sigma$ point source, $r$ band) and a cadence of 5 days can detect 0.49 strongly-lensed SNIa event and 2.1 strongly-lensed CCSN events per year within this sample. Our proposed targeted-search strategy is particularly useful for prompt and efficient identifications and follow-up observations of strongly-lensed SN candidates. It also allows telescopes with small field of views and limited time to efficiently discover strongly-lensed SNe with a pencil-beam scanning strategy.


title: Optimal Target Stars in the Search for Life

abstract: The selection of optimal targets in the search for life represents a highly important strategic issue. In this Letter, we evaluate the benefits of searching for life around a potentially habitable planet orbiting a star of arbitrary mass relative to a similar planet around a Sun-like star. If recent physical arguments implying that the habitability of planets orbiting low-mass stars is selectively suppressed are correct, we find that planets around solar-type stars may represent the optimal targets.


title: A Potential Cyclotron Resonant Scattering Feature in the ULX Pulsar NGC 300 ULX1 seen by NuSTAR and XMM-Newton

abstract: Based on phase-resolved broadband spectroscopy using $XMM$-$Newton$ and $NuSTAR$, we report on a potential cyclotron resonant scattering feature at $E \sim 13$ keV in the pulsed spectrum of the recently discoverd ULX pulsar NGC 300 ULX1. If this interpretation is correct, the implied magnetic field of the central neutron star is $B \sim 10^{12}$ G (assuming scattering off electrons), similar to that estimated from the observed spin-up of the star, and also similar to known Galactic X-ray pulsars. We discuss the implications of this result for the connection between NGC 300 ULX1 and the other known ULX pulsars, particularly in light of the recent discovery of a likely proton Cyclotron line in another ULX, M51 ULX-8.


title: Co-existence and switching between fast and $\Omega$-slow wind solutions in rapidly rotating massive stars

abstract: Most radiatively-driven massive star winds can be modelled with m-CAK theory resulting in so called fast solution. However, those most rapidly rotating among them, especially when the stellar rotational speed is higher than $\sim 75\%$ of the critical rotational speed, can adopt a different solution called $\Omega$-slow solution characterized by a dense and slow wind. Here, in this work we study the transition region of the solutions where the fast solution changes to the $\Omega$-slow. Using both time-steady and time-dependent numerical codes, we study this transition region for different equatorial models of B-type stars. In all the cases, at certain range of rotational speeds, we found a region where the fast and $\Omega$-slow solution can co-exist. We find that the type of solution obtained in this co-existence region depends heavily on the initial conditions of our models. We also test the stability of the solutions within the co-existence region by performing base density perturbations in the wind. We find that under certain conditions, the fast solution can switch to a $\Omega$-slow solution, or vice versa. Such solution switching may be a possible contributor of material injected into the circumstellar environment of Be stars, without requiring rotational speeds near critical values.


title: Numerical Problems in Coupling Photon Momentum (Radiation Pressure) to Gas

abstract: Radiation pressure (RP; or photon momentum absorbed by gas) is important in a tremendous range of astrophysical systems. But we show the usual method for assigning absorbed photon momentum to gas in numerical radiation-hydrodynamics simulations (integrating over cell volumes or evaluating at cell centers) can severely under-estimate the RP force in the immediate vicinity around un-resolved (point/discrete) sources (and subsequently under-estimate its effects on bulk gas properties), unless photon mean-free-paths are highly-resolved in the fluid grid. The existence of this error is independent of the numerical radiation transfer (RT) method (even in exact ray-tracing/Monte-Carlo methods), because it depends on how the RT solution is interpolated back onto fluid elements. Brute-force convergence (resolving mean-free paths) is impossible in many cases (especially where UV/ionizing photons are involved). Instead, we show a 'face-integrated' method -- integrating and applying the momentum fluxes at interfaces between fluid elements -- better approximates the correct solution at all resolution levels. The 'fix' is simple and we provide example implementations for ray-tracing, Monte-Carlo, and moments RT methods in both grid and mesh-free fluid schemes. We consider an example of star formation in a molecular cloud with UV/ionizing RP. At state-of-the-art resolution, cell-integrated methods under-estimate the net effects of RP by an order of magnitude, leading (incorrectly) to the conclusion that RP is unimportant, while face-integrated methods predict strong self-regulation of star formation and cloud destruction via RP.

categories: astro-ph.GA astro-ph.CO astro-ph.HE astro-ph.IM astro-ph.SR


title: Water Residence Time Estimation by 1D Deconvolution in the Form of a l2-Regularized Inverse Problem With Smoothness, Positivity and Causality Constraints

abstract: The Water Residence Time distribution is the equivalent of the impulse response of a linear system allowing the propagation of water through a medium, e.g. the propagation of rain water from the top of the mountain towards the aquifers. We consider the output aquifer levels as the convolution between the input rain levels and the Water Residence Time, starting with an initial aquifer base level. The estimation of Water Residence Time is important for a better understanding of hydro-bio-geochemical processes and mixing properties of wetlands used as filters in ecological applications, as well as protecting fresh water sources for wells from pollutants. Common methods of estimating the Water Residence Time focus on cross-correlation, parameter fitting and non-parametric deconvolution methods. Here we propose a 1D full-deconvolution, regularized, non-parametric inverse problem algorithm that enforces smoothness and uses constraints of causality and positivity to estimate the Water Residence Time curve. Compared to Bayesian non-parametric deconvolution approaches, it has a fast runtime per test case; compared to the popular and fast cross-correlation method, it produces a more precise Water Residence Time curve even in the case of noisy measurements. The algorithm needs only one regularization parameter to balance between smoothness of the Water Residence Time and accuracy of the reconstruction. We propose an approach on how to automatically find a suitable value of the regularization parameter from the input data only. Tests on real data illustrate the potential of this method to analyze hydrological datasets.


title: Light axions with large hadronic couplings

abstract: We present models in which axions have large couplings to hadrons while remaining naturally light. By assigning $U(1)_{PQ}$ charges to Standard Model quarks such that $U(1)_{PQ}$ is not anomalous under QCD, the derivative couplings naturally arise while no potential is generated for the axion upon QCD confinement. We present simple models that implement this idea both for an axion-like particle and for the QCD axion. We show that both models lead to sizable flavor violations that could be probed by future experiments. Our construction shows that the axion coupling to hadrons can be essentially independent from its mass and thus motivates experimental searches in a wide range of axion parameter space.


title: Quasinormal modes as a distinguisher between general relativity and f(R) gravity: Charged black-holes

abstract: The ring-down phase of black-hole perturbations is governed by the Quasi-Normal modes (QNM) and offer valuable insight into the nature of the objects emitting them. In General relativity, we identify a dimensionless parameter that can estimate the charge of a black-hole from the energetics of quasi-normal modes. We then extend the analysis for f(R) theories of gravity. Using Isaacson's prescription to modified theories of gravity, we obtain an accurate and robust method to estimate the difference in the radiated energies in the scalar and vector modes due to modifications to general relativity. We then discuss how to quantify the deviation from general relativity.


title: Stellar Population Properties of Ultracompact Dwarfs in M87: a Mass-metallicity Correlation Connecting Low-metallicity Globular Clusters and Compact Ellipticals

abstract: We derive stellar population parameters for a representative sample of ultracompact dwarfs (UCDs) and a large sample of massive globular clusters (GCs) with stellar masses $\gtrsim$ 10$^{6}$ $M_{\odot}$ in the central galaxy M87 of the Virgo galaxy cluster, based on model fitting to the Lick-index measurements from both the literature and new observations. After necessary spectral stacking of the relatively faint objects in our initial sample of 40 UCDs and 118 GCs, we obtain 30 sets of Lick-index measurements for UCDs and 80 for GCs. The M87 UCDs have ages $\gtrsim$ 8 Gyr and [$\alpha$/Fe] $\simeq$ 0.4 dex, in agreement with previous studies based on smaller samples. The literature UCDs, located in lower-density environments than M87, extend to younger ages and smaller [$\alpha$/Fe] (at given metallicities) than M87 UCDs, resembling the environmental dependence of the Virgo dE nuclei. The UCDs exhibit a positive mass-metallicity relation (MZR), which flattens and connects compact ellipticals at stellar masses $\gtrsim$ 10$^{8}$ $M_{\odot}$. The Virgo dE nuclei largely follow the average MZR of UCDs, whereas most of the M87 GCs are offset towards higher metallicities for given stellar masses. The difference between the mass-metallicity distributions of UCDs and GCs may be qualitatively understood as a result of their different physical sizes at birth in a self-enrichment scenario or of galactic nuclear cluster star formation efficiency being relatively low in a tidal stripping scenario for UCD formation. The existing observations provide the necessary but not sufficient evidence for tidally stripped dE nuclei being the dominant contributors to the M87 UCDs.


title: Toward a compact fibered squeezing parametric source

abstract: In this work, we investigate three different compact fibered systems generating vacuum squeezing that involve optical cavities limited by the end surface of a fiber and by a curved mirror and containing a thin parametric crystal. These systems have the advantage to couple squeezed states directly to a fiber, allowing the user to benefit from the flexibility of fibers in the use of squeezing. Three types of fibers are investigated: standard single-mode fibers, photonic-crystal large mode area single-mode fibers and short multi-mode fiber taped to a single-mode fiber. The observed squeezing is modest (-0.56 dB, -0.9 dB, -1 dB), but these experiments open the way to miniaturized squeezing devices that could be a very interesting advantage in scaling up quantum systems for quantum processing, opening new perspectives in the domain of integrated quantum optics.


title: Schrodinger-Maxwell systems on compact Riemannian manifolds

abstract: In this paper we are focusing to the following Schr\"odinger-Maxwell system $(\mathcal{SM}_{\Psi(\lambda,\cdot)}^{e})$: \[ \begin{cases} -\Delta_{g}u+\beta(x)u+eu\phi=\Psi(\lambda,x)f(u) & \mathrm{in}\ M -\Delta_{g}\phi+\phi=qu^{2} & \mathrm{\mathrm{in}\ M} \end{cases} \] where $(M,g)$ is a 3-dimensional compact Riemannian manifold without boundary, $e,q>0$ are positive numbers, $f:\mathbb{R}\to\mathbb{R}$ is a continuous function, $\beta\in C^{\infty}(M)$ and $\Psi\in C^{\infty}(\mathbb{R}_{+}\times M)$ are positive functions. By various variational approaches, existence of multiple solutions of the problem $(\mathcal{SM}_{\Psi(\lambda,\cdot)}^{e})$ is established.


title: Resource theory of non-Gaussian operations

abstract: Non-Gaussian states and operations are crucial for various continuous-variable quantum information processing tasks. To quantitatively understand non-Gaussianity beyond states, we establish a resource theory for non-Gaussian operations. In our framework, we consider Gaussian operations as free operations, and non-Gaussian operations as resources. We define entanglement-assisted non-Gaussianity generating power and show that it is a monotone that is non-increasing under the set of free super-operations, i.e., concatenation and tensoring with Gaussian channels. For conditional unitary maps, this monotone can be analytically calculated. As examples, we show that the non-Gaussianity of ideal photon-number subtraction and photon-number addition equal the non-Gaussianity of the single-photon Fock state. Based on our non-Gaussianity monotone, we divide non-Gaussian operations into two classes: (1) the finite non-Gaussianity class, e.g., photon-number subtraction, photon-number addition and all Gaussian-dilatable non-Gaussian channels; and (2) the diverging non-Gaussianity class, e.g., the binary phase-shift channel and the Kerr nonlinearity. This classification also implies that not all non-Gaussian channels are exactly Gaussian-dilatable. Our resource theory enables a quantitative characterization and a first classification of non-Gaussian operations, paving the way towards the full understanding of non-Gaussianity.


title: Packing and covering induced subdivisions

abstract: A class $\mathcal{F}$ of graphs has the induced Erd\H{o}s-P\'osa property if there exists a function $f$ such that for every graph $G$ and every positive integer $k$, $G$ contains either $k$ pairwise vertex-disjoint induced subgraphs that belong to $\mathcal{F}$, or a vertex set of size at most $f(k)$ hitting all induced copies of graphs in $\mathcal{F}$. Kim and Kwon (SODA'18) showed that for a cycle $C_{\ell}$ of length $\ell$, the class of $C_{\ell}$-subdivisions has the induced Erd\H{o}s-P\'osa property if and only if $\ell\le 4$. In this paper, we investigate whether or not the class of $H$-subdivisions has the induced Erd\H{o}s-P\'osa property for other graphs $H$. We completely settle the case when $H$ is a forest or a complete bipartite graph. Regarding the general case, we identify necessary conditions on $H$ for the class of $H$-subdivisions to have the induced Erd\H{o}s-P\'osa property. For this, we provide three basic constructions that are useful to prove that the class of the subdivisions of a graph does not have the induced Erd\H{o}s-P\'osa property. Among remaining graphs, we prove that if $H$ is either the diamond, the $1$-pan, or the $2$-pan, then the class of $H$-subdivisions has the induced Erd\H{o}s-P\'osa property.


title: Indirect Influences, Links Ranking, and Deconstruction of Networks

abstract: The PWP map was introduced by the second author as a tool for ranking nodes in networks. In this work we extend this technique so that it can be used to rank links as well. Applying the Girvan-Newman algorithm a ranking method on links induces a deconstruction method for networks, therefore we obtain new methods for finding clustering and core-periphery structures on networks.


title: Long-range entanglement near a Kondo-destruction quantum critical point

abstract: The numerical renormalization group is used to study quantum entanglement in the Kondo impurity model with a pseudogapped density of states $\rho(\varepsilon)\propto|\varepsilon|^r$ ($r>0$) that vanishes at the Fermi energy $\varepsilon=0$. The model features a Kondo-destruction quantum critical point (QCP) separating a partially screened phase (reached for impurity-band exchange couplings $J>J_c$) from a local-moment phase ($J<J_c$). The impurity contribution $S_e^{imp}$ to the entanglement entropy between a region of radius $R$ around the magnetic impurity and the rest of the host system reveals a characteristic length scale $R^*$ that distinguishes a regime $R\ll R^*$ of maximal critical entanglement from one $R\gg R^*$ of weaker entanglement. Within each phase, $S_e^{imp}$ is a universal function of $R/R^*$ with a power-law decay for $R/R^*\gg 1$. The entanglement length scale $R^*$ diverges on approach to the QCP with a critical exponent that depends only on $r$.


title: Quartz tuning-fork based carbon nanotube transfer into quantum device geometries

abstract: With the objective of integrating single clean, as-grown carbon nanotubes into complex circuits, we have developed a technique to grow nanotubes directly on commercially available quartz tuning forks using a high temperature CVD process. Multiple straight and aligned nanotubes bridge the >100um gap between the two tips. The nanotubes are then lowered onto contact electrodes, electronically characterized in situ, and subsequently cut loose from the tuning fork using a high current. First quantum transport measurements of the resulting devices at cryogenic temperatures display Coulomb blockade characteristics.


title: Homological classification of topological terms in sigma models on homogeneous spaces

abstract: We classify the topological terms (in a sense to be made precise) that may appear in a non-linear sigma model based on maps from an arbitrary worldvolume manifold to a homogeneous space $G/H$ (where $G$ is an arbitrary Lie group and $H \subset G$). We derive a new condition for $G$-invariance of topological terms, which is necessary and sufficient (at least when $G$ is connected), and discuss a variety of examples in quantum mechanics and quantum field theory. In the present work we discuss only terms that may be written in terms of (possibly only locally-defined) differential forms on $G/H$, leading to an action that is manifestly local. Such terms come in one of two types, with prototypical quantum-mechanical examples given by the Aharonov-Bohm effect and the Dirac monopole. The classification is based on the observation that, for topological terms, the maps from the worldvolume to $G/H$ may be replaced by singular homology cycles on $G/H$. In a forthcoming paper we apply the results to phenomenological models in which the Higgs boson is composite.


title: Low-Complexity Distributed Radio Access Network Slicing: Algorithms and Experimental Results

abstract: Radio access network (RAN) slicing is an effective methodology to dynamically allocate networking resources in 5G networks. One of the main challenges of RAN slicing is that it is provably an NP-Hard problem. For this reason, we design near-optimal low-complexity distributed RAN slicing algorithms. First, we model the slicing problem as a congestion game, and demonstrate that such game admits a unique Nash equilibrium (NE). Then, we evaluate the Price of Anarchy (PoA) of the NE, i.e., the efficiency of the NE as compared to the social optimum, and demonstrate that the PoA is upper-bounded by 3/2. Next, we propose two fully-distributed algorithms that provably converge to the unique NE without revealing privacy-sensitive parameters from the slice tenants. Moreover, we introduce an adaptive pricing mechanism of the wireless resources to improve the network owner's profit. We evaluate the performance of our algorithms through simulations and an experimental testbed deployed on the Amazon EC2 cloud, both based on a real-world dataset of base stations from the OpenCellID project. Results conclude that our algorithms converge to the NE rapidly and achieve near-optimal performance, while our pricing mechanism effectively improves the profit of the network owner.


title: HINT: A Hierarchical Independent Component Analysis Toolbox for Investigating Brain Functional Networks using Neuroimaging Data

abstract: Independent component analysis (ICA) is a popular tool for investigating brain organization in neuroscience research. In fMRI studies, an important goal is to study how brain networks are modulated by subjects' clinical and demographic variables. Existing ICA methods and toolboxes don't incorporate subjects' covariates effects in ICA estimation of brain networks, which potentially leads to loss in accuracy and statistical power in detecting brain network differences between subjects' groups. We introduce a Matlab toolbox, HINT (Hierarchical INdependent component analysis Toolbox), that provides a hierarchical covariate-adjusted ICA (hc-ICA) for modeling and testing covariate effects and generates model-based estimates of brain networks on both the population- and individual-level. HINT provides a user-friendly Matlab GUI that allows users to easily load images, specify covariate effects, monitor model estimation via an EM algorithm, specify hypothesis tests, and visualize results. HINT also has a command line interface which allows users to conveniently run and reproduce the analysis with a script. HINT implements a new multi-level probabilistic ICA model for group ICA. It provides a statistically principled ICA modeling framework for investigating covariate effects on brain networks. HINT can also generate and visualize model-based network estimates for user-specified subject groups, which greatly facilitates group comparisons.


title: A Push-Pull Gradient Method for Distributed Optimization in Networks

abstract: In this paper, we focus on solving a distributed convex optimization problem in a network, where each agent has its own convex cost function and the goal is to minimize the sum of the agents' cost functions while obeying the network connectivity structure. In order to minimize the sum of the cost functions, we consider a new distributed gradient-based method where each node maintains two estimates, namely, an estimate of the optimal decision variable and an estimate of the gradient for the average of the agents' objective functions. From the viewpoint of an agent, the information about the decision variable is pushed to the neighbors, while the information about the gradients is pulled from the neighbors (hence giving the name "push-pull gradient method"). The method unifies the algorithms with different types of distributed architecture, including decentralized (peer-to-peer), centralized (master-slave), and semi-centralized (leader-follower) architecture. We show that the algorithm converges linearly for strongly convex and smooth objective functions over a directed static network. In our numerical test, the algorithm performs well even for time-varying directed networks.


title: The Swift/BAT AGN Spectroscopic Survey -- IX. The Clustering Environments of an Unbiased Sample of Local AGN

abstract: We characterize the environments of local accreting supermassive black holes by measuring the clustering of AGN in the Swift/BAT Spectroscopic Survey (BASS). With 548 AGN in the redshift range 0.01<z<0.1 over the full sky from the DR1 catalog, BASS provides the largest, least biased sample of local AGN to date due to its hard X-ray selection (14-195 keV) and rich multiwavelength/ancillary data. By measuring the projected cross-correlation function between the AGN and 2MASS galaxies, and interpreting it via halo occupation distribution (HOD) and subhalo-based models, we constrain the occupation statistics of the full sample, as well as in bins of absorbing column density and black hole mass. We find that AGN tend to reside in galaxy group environments, in agreement with previous studies of AGN throughout a large range of luminosity and redshift, and that on average they occupy their dark matter halos similar to inactive galaxies of comparable stellar mass. We also find evidence that obscured AGN tend to reside in denser environments than unobscured AGN, even when samples were matched in luminosity, redshift, stellar mass, and Eddington ratio. We show that this can be explained either by significantly different halo occupation distributions or statistically different host halo assembly histories. Lastly, we see that massive black holes are slightly more likely to reside in central galaxies than black holes of smaller mass.


title: On the Basel Liquidity Formula for Elliptical Distributions

abstract: A justification of the Basel liquidity formula for risk capital in the trading book is given under the assumption that market risk-factor changes form a Gaussian white noise process over 10-day time steps and changes to P&L are linear in the risk-factor changes. A generalization of the formula is derived under the more general assumption that risk-factor changes are multivariate elliptical. It is shown that the Basel formula tends to be conservative when the elliptical distributions are from the heavier-tailed generalized hyperbolic family. As a by-product of the analysis a Fourier approach to calculating expected shortfall for general symmetric loss distributions is developed.


title: A probabilistic variant of Sperner's theorem and of maximal $r$-cover free families

abstract: A family of sets is called $r$-\emph{cover free} if no set in the family is contained in the union of $r$ (or less) other sets in the family. A $1$-cover free family is simply an antichain with respect to set inclusion. Thus, Sperner's classical result determines the maximal cardinality of a $1$-cover free family of subsets of an $n$-element set. Estimating the maximal cardinality of an $r$-cover free family of subsets of an $n$-element set for $r>1$ was also studied. In this note we are interested in the following probabilistic variant of this problem. Let $S_0,S_1,\ldots, S_r$ be independent and identically distributed random subsets of an $n$-element set. Which distribution minimizes the probability that $S_0\subseteq {\bigcup_{i=1}^r S_i}$? A natural candidate is the uniform distribution on an $r$-cover-free family of maximal cardinality. We show that for $r=1$ such distribution is indeed best possible. In a complete contrast, we also show that this is far from being true for every $r>1$ and $n$ large enough.


title: Critical domains for the first nonzero Neumann eigenvalue in Riemannian manifolds

abstract: The present paper is devoted to geometric optimization problems related to the Neumann eigenvalue problem for the Laplace-Beltrami operator on bounded subdomains $\Omega$ of a Riemannian manifold $(\mathcal{M},g)$. More precisely, we analyze locally extremal domains for the first nontrivial eigenvalue $\mu_2(\Omega)$ with respect to volume preserving domain perturbations, and we show that corresponding notions of criticality arise in the form of overdetermined boundary problems. Our results rely on an extension of Zanger's shape derivative formula which covers the case when $\mu_2(\Omega)$ is not a simple eigenvalue. In the second part of the paper, we focus on product manifolds of the form $\mathcal{M} = \mathbb{R}^k \times \mathcal{N}$, and we classify the subdomains where an associated overdetermined boundary value problem has a solution.


title: Statistical Analysis of Solar Events Associated with Storm Sudden Commencements over One Year of Solar Maximum during Cycle 23: Propagation from the Sun to the Earth and Effects

abstract: Taking the 32 storm sudden commencements (SSCs) listed by ISGI during 2002 as a starting point, we performed a multi-criterion analysis based on observations (propagation time, velocity comparisons, sense of the magnetic field rotation, radio waves) to associate them with solar sources. We identified their effects in the interplanetary medium, and looked at the response of the terrestrial ionized and neutral environment. We find that 28 SSCs can be related to 44 coronal mass ejections (CMEs), 15 with a unique CME and 13 with a series of multiple CMEs, among which 19 involved halo CMEs; 12 of the 19 fastest CMEs with speeds greater than 1000 km/s are halo CMEs. The probability for an SSC to occur is 75% if the CME is a halo CME. The complex interactions between two CMEs and the modification of their trajectories have been examined using joint white-light and multiple-wavelength radio observations. The solar-wind structures at L1 after the shocks leading the 32 SSCs are 12 magnetic clouds (MCs), 6 interplanetary CMEs (ICMEs) without an MC structure, 4 miscellaneous structures, which cannot unambiguously be classified as ICMEs, 5 corotating or stream interaction regions (CIRs/SIRs), and 4 isolated shock events; note than one CIR caused two SSCs. The 11 MCs listed in 3 or more MC catalogs covering the year 2002 are associated with SSCs. For the 3 most intense geomagnetic storms related to MCs, we note 2 sudden increases of the Dst, at the arrival of the sheath and the arrival of the MC itself. The most geoeffective events are MCs, since 92% of them trigger moderate or intense storms, followed by ICMEs (33%). At best, CIRs/SIRs only cause weak storms. We show that these geoeffective events (ICMEs or MCs) trigger or reinforce terrestrial radiowave activity in the magnetosphere, an enhanced convection in the ionosphere, and a stronger response in the thermosphere.


title: Hydrodynamic Flows on Curved Surfaces: Spectral Numerical Methods for Radial Manifold Shapes

abstract: We formulate hydrodynamic equations and spectrally accurate numerical methods for investigating the role of geometry in flows within two-dimensional fluid interfaces. To achieve numerical approximations having high precision and level of symmetry for radial manifold shapes, we develop spectral Galerkin methods based on hyperinterpolation with Lebedev quadratures for $L^2$-projection to spherical harmonics. We demonstrate our methods by investigating hydrodynamic responses as the surface geometry is varied. Relative to the case of a sphere, we find significant changes can occur in the observed hydrodynamic flow responses as exhibited by quantitative and topological transitions in the structure of the flow. We present numerical results based on the Rayleigh-Dissipation principle to gain further insights into these flow responses. We investigate the roles played by the geometry especially concerning the positive and negative Gaussian curvature of the interface. We provide general approaches for taking geometric effects into account for investigations of hydrodynamic phenomena within curved fluid interfaces.


title: From $\gamma$ to Radio - The Electromagnetic Counterpart of GW 170817

abstract: The gravitational waves from the first binary neutron star merger, GW170817, were accompanied by a multi-wavelength electromagnetic counterpart, from $\gamma$-rays to radio. The accompanying gamma-rays, seems at first to confirm the association of mergers with short gamma-ray bursts (sGRBs). The common interpretation was that we see an emission from an sGRB jet seen off-axis. However, a closer examination of the sub-luminous $\gamma$-rays and the peculiar radio afterglow were inconsistent with this simple interpretation. Here we present results of 3D and 2D numerical simulations that follow the hydrodynamics and emission of the outflow from a neutron star merger form its ejection and up to its deceleration by the circum-merger medium. Our results show that the entire set of $\gamma$-rays, X-rays and radio observations can be explained by the emission from a mildly relativistic cocoon material (Lorentz factor $\sim$2-5) that was formed while a jet propagated through the material ejected during the merger. The $\gamma$-rays are generated when the cocoon breaks out from the engulfing ejecta while the afterglow is produced by interaction of the cocoon matter with the interstellar medium. The strong early uv/optical signal may be a Lorentz boosted macronova/kilonova. The fate of the jet itself is currently unknown, but our full-EM models define a path to resolving between successful and choked jet scenarios, outputting coupled predictions for the image size, morphology, observed time-dependent polarization and light curve behavior from radio to X-ray. The predictive power of these models will prove key in interpreting the on-going multi-faceted observations of this unprecedented event.


title: Mumford divisors

abstract: We define the notion of Mumford divisors, argue that they are the natural divisors to study on reduced but non-normal varieties and prove a structure theorem for the Mumford class group. v.2: references added


title: Diffusion-driven transition between two regimes of viscous fingering

abstract: Viscous fingering patterns can form at the interface between two immiscible fluids confined in the gap between a pair of flat plates; whenever the fluid with lower viscosity displaces the one of higher viscosity the interface is unstable. For miscible fluids the situation is more complicated due to the formation of interfacial structure in the thin dimension spanning the gap. Here we study the effect of the inherent diffusion between the two miscible fluids on this structure and on the viscous fingering patterns that emerge. We discover an unexpected transition separating two distinct regimes where the pattern morphologies and mode of onset are different. This transition is marked by a regime of transient stability as the structure of the fingers evolves from having three-dimensional structure to being quasi-two dimensional. The presence of diffusion allows an instability to form where it was otherwise forbidden.


title: Secondary ion mass spectrometry for SRF cavity materials

abstract: Historically, many advances in superconducting radio frequency (SRF) cavities destined for use in advanced particle accelerators have come empirically, through the iterative procedure of modifying processing and then performance testing. However, material structure is directly responsible for performance. Understanding, the link between processing, structure, and performance will streamline and accelerate the research process. In order to connect processing, structure, and performance, accurate and robust materials characterization methods are needed. Here one such method, SIMS, is discussed with focus on analysis of SRF materials. In addition, several examples are presented, showing how SIMS is being used to further understanding of materials based SRF technologies.


title: Thermal to Visible Synthesis of Face Images using Multiple Regions

abstract: Synthesis of visible spectrum faces from thermal facial imagery is a promising approach for heterogeneous face recognition; enabling existing face recognition software trained on visible imagery to be leveraged, and allowing human analysts to verify cross-spectrum matches more effectively. We propose a new synthesis method to enhance the discriminative quality of synthesized visible face imagery by leveraging both global (e.g., entire face) and local regions (e.g., eyes, nose, and mouth). Here, each region provides (1) an independent representation for the corresponding area, and (2) additional regularization terms, which impact the overall quality of synthesized images. We analyze the effects of using multiple regions to synthesize a visible face image from a thermal face. We demonstrate that our approach improves cross-spectrum verification rates over recently published synthesis approaches. Moreover, using our synthesized imagery, we report the results on facial landmark detection-commonly used for image registration-which is a critical part of the face recognition process.


title: Axionic black branes with conformal coupling

abstract: We find neutral and charged black branes solutions with axion fields in the context of a conformally coupled gravitational theory in four dimensions. These solutions describe AdS black branes supported by axion fields that break translational invariance at the boundary, providing for momentum dissipation. The conformally coupled scalar field is regular inside and outside the event horizon and there is no need of any self-interaction, obtaining in this way solutions without fine-tuned parameters. We analyze the thermodynamics of our solutions considering the effects of the axion charges and it is shown that axionic and electric charges must be related such that the conformal scalar field does not contribute to the mass. We compute the holographic DC conductivity and we show how it is affected by the inclusion of the conformal scalar field, which provides a temperature independent behavior. We include a non-linear axionic contribution given by a k-essence term that modifies the DC conductivity providing for more general behaviors. Finally, we endorse our solutions with rotation showing that angular momentum is sustained by the axion charges.


title: Gravitationally Lensed Quasars in Gaia: II. Discovery of 24 Lensed Quasars

abstract: We report the discovery, spectroscopic confirmation and preliminary characterisation of 24 gravitationally lensed quasars identified using Gaia observations. Candidates were selected in the Pan-STARRS footprint with quasar-like WISE colours or as photometric quasars from SDSS, requiring either multiple detections in Gaia or a single Gaia detection near a morphological galaxy. The Pan-STARRS grizY images were modelled for the most promising candidates and 60 candidate systems were followed up with the William Herschel Telescope. 13 of the lenses were discovered as Gaia multiples and 10 as single Gaia detections near galaxies. We also discover 1 lens identified through a quasar emission line in an SDSS galaxy spectrum. The lenses have median image separation 2.13 arcsec and the source redshifts range from 1.06 to 3.36. 4 systems are quadruply-imaged and 20 are doubly-imaged. Deep CFHT data reveal an Einstein ring in one double system. We also report 12 quasar pairs, 10 of which have components at the same redshift and require further follow-up to rule out the lensing hypothesis. We compare the properties of these lenses and other known lenses recovered by our search method to a complete sample of simulated lenses to show the lenses we are missing are mainly those with small separations and higher source redshifts. The initial Gaia data release only catalogues all images of ~ 30% of known bright lensed quasars, however the improved completeness of Gaia data release 2 will help find all bright lensed quasars on the sky.


title: UnibucKernel: A kernel-based learning method for complex word identification

abstract: In this paper, we present a kernel-based learning approach for the 2018 Complex Word Identification (CWI) Shared Task. Our approach is based on combining multiple low-level features, such as character n-grams, with high-level semantic features that are either automatically learned using word embeddings or extracted from a lexical knowledge base, namely WordNet. After feature extraction, we employ a kernel method for the learning phase. The feature matrix is first transformed into a normalized kernel matrix. For the binary classification task (simple versus complex), we employ Support Vector Machines. For the regression task, in which we have to predict the complexity level of a word (a word is more complex if it is labeled as complex by more annotators), we employ v-Support Vector Regression. We applied our approach only on the three English data sets containing documents from Wikipedia, WikiNews and News domains. Our best result during the competition was the third place on the English Wikipedia data set. However, in this paper, we also report better post-competition results.


title: Schopenhauer on Space, Time, Causality and Matter: A Physical Re-examination

abstract: According to Schopenhauer, Kant's arguments about the transcendental ideality of space and time can be extended to matter through the concept of causality and the principle of sufficient reason. In this article, I examine to what extent space, time and causality can be considered a priori concepts in the light of classical and modern physics. The concepts of matter and field, and their possible a priori fundaments, as stated by Schopenhauer, are thus revisited in a modern context.


title: Continuous cohomology of topological quandles

abstract: A continuous cohomology theory for topological quandles is introduced, and compared to the algebraic theories. Extensions of topological quandles are studied with respect to continuous 2-cocycles, and used to show the differences in second cohomology groups for specific topological quandles. A method of computing the cohomology groups of the inverse limit is applied to quandles.


title: The QCD running coupling at all scales and the connection between hadron masses and \Lambda_s

abstract: We report on recent experimental and theoretical developments in our understanding of the QCD running coupling \alpha_s in QCD's nonperturbative regime. They allow us to analytically compute the hadron mass spectrum, with \Lambda_s the only input necessary to this determination. The computed spectra agrees well with experimental data.


title: Couette-Poiseuille flow with partial slip and uniform cross flow for power-law fluids

abstract: Exact solutions are obtained for the steady flow of a power-law fluid between parallel plates with partial slip conditions and uniform cross flow. The problem is properly formulated and similarities are exploited. The exact solutions are obtained in terms of integrals which can be performed, in closed form, in special cases of the power-law index n. Solutions to cases of n=1/2, 1, and 2; representing a pseudo-plastic, a Newtonian, and a dilatant fluid, respectively, are presented. Tendencies to corresponding degenerate cases in the literature are demonstrated. Depending on the strength of the cross flow and the pressure gradient, the flow may be of Couette type with convex, linear, or concave velocity profile; or of Poiseuille type. Borderline cases are identified. Moreover, a case in which the power-law model for the dilatant fluid fails is detected.


title: Surfactant-induced-retardation in lateral migration of droplets in a microfluidic confinement

abstract: In the present study, the role of surfactants on the cross-stream migration of droplets is investigated both experimentally and theoretically. For experimental analysis, sunflower oil is used as the carrier phase and DI water as the dispersed phase which is intermixed with Triton X-100 that acts as the surfactant. A T-junction is used in the microchannel for the purpose of droplet generation. Presence of an imposed pressure driven flow induces droplet deformation and disturbs the equilibrium that results in subsequent surfactant-redistribution along the interface. This further induces a gradient in the surface tension, thus generating a Marangoni stress that significantly alters the droplet dynamics. On subsequent experimental investigation, it is found that presence of surfactants reduces the cross-stream migration velocity of the droplet. Further, it is shown that the effect of surfactants in reducing the cross-stream migration is significantly enhanced for a larger droplet as compared to a smaller one within the same time span, provided the channel height is kept constant. In addition, a larger surfactant concentration is found to induce a greater retardation in cross-stream migration of the droplet, the effect of which is reduced when the initial transverse position of the droplet is shifted closer to the channel centerline. The present analysis can be applied to various droplet based microfluidic as well as medical diagnostic devices where manipulation of droplet trajectory is a major issue. For the theoretical counterpart, an asymptotic approach is adopted in the presence of bulk-insoluble surfactants and under the assumption of negligible fluid inertia. A good match between our theoretical prediction and the experimental results is obtained.


title: A Survey of Deep Learning Techniques for Mobile Robot Applications

abstract: Advancements in deep learning over the years have attracted research into how deep artificial neural networks can be used in robotic systems. This research survey will present a summarization of the current research with a specific focus on the gains and obstacles for deep learning to be applied to mobile robotics.


title: The $\ell^\infty$-Cophenetic Metric for Phylogenetic Trees as an Interleaving Distance

abstract: There are many metrics available to compare phylogenetic trees since this is a fundamental task in computational biology. In this paper, we focus on one such metric, the $\ell^\infty$-cophenetic metric introduced by Cardona et al. This metric works by representing a phylogenetic tree with $n$ labeled leaves as a point in $\mathbb{R}^{n(n+1)/2}$ known as the cophenetic vector, then comparing the two resulting Euclidean points using the $\ell^\infty$ distance. Meanwhile, the interleaving distance is a formal categorical construction generalized from the definition of Chazal et al., originally introduced to compare persistence modules arising from the field of topological data analysis. We show that the $\ell^\infty$-cophenetic metric is an example of an interleaving distance. To do this, we define phylogenetic trees as a category of merge trees with some additional structure; namely labelings on the leaves plus a requirement that morphisms respect these labels. Then we can use the definition of a flow on this category to give an interleaving distance. Finally, we show that, because of the additional structure given by the categories defined, the map sending a labeled merge tree to the cophenetic vector is, in fact, an isometric embedding, thus proving that the $\ell^\infty$-cophenetic metric is, in fact, an interleaving distance.


title: Exploring the use of time-varying graphs for modelling transit networks

abstract: The study of the dynamic relationship between topological structure of a transit network and the mobility patterns of transit vehicles on this network is critical towardsdevising smart and time-aware solutions to transit management and recommendation systems. This paper proposes a time-varying graph (TVG) to model thisrelationship. The effectiveness of this proposed model has been explored by implementing the model in Neo4j graph database using transit feeds generated by bus transit network of the City of Moncton, New Brunswick, Canada. Dynamics in this relationshipalsohave been detected using network metrics such as temporal shortest paths, degree, betweenness and PageRank centralities as well as temporal network diameter and density. Keywords: Transit Networks,Mobility Pattern,Time-Varying Graph model, Graph Databaseand Graph Analytics Keywords: Transit Networks,Mobility Pattern,Time-Varying Graph model, Graph Database and Graph Analytics


title: Compact complex surfaces with no nonconstant meromorphic functions

abstract: In 1949 Siegel gave an example of a complex two-torus with no nonconstant meromorphic functions. In 1964 Kodaira showed that compact complex surfaces with no nonconstant meromorphic must be of the following three types: tori, Hopf type surfaces with first Betti number equal to one, and K3 surfaces. In this paper we show that surfaces of these three types have a dense set of surfaces in their natural moduli spaces with no nonconstant meromorphic functions.


title: Generating Multi-Agent Trajectories using Programmatic Weak Supervision

abstract: We study the problem of training sequential generative models for capturing coordinated multi-agent trajectory behavior, such as offensive basketball gameplay. When modeling such settings, it is often beneficial to design hierarchical models that can capture long-term coordination using intermediate variables. Furthermore, these intermediate variables should capture interesting high-level behavioral semantics in an interpretable and manipulatable way. We present a hierarchical framework that can effectively learn such sequential generative models. Our approach is inspired by recent work on leveraging programmatically produced weak labels, which we extend to the spatiotemporal regime. In addition to synthetic settings, we show how to instantiate our framework to effectively model complex interactions between basketball players and generate realistic multi-agent trajectories of basketball gameplay over long time periods. We validate our approach using both quantitative and qualitative evaluations, including a user study comparison conducted with professional sports analysts.


title: Integrating DRAM Power-Down Modes in gem5 and Quantifying their Impact

abstract: Across applications, DRAM is a significant contributor to the overall system power, with the DRAM access energy per bit up to three orders of magnitude higher compared to on-chip memory accesses. To improve the power efficiency, DRAM technology incorporates multiple power-down modes, each with different trade-offs between achievable power savings and performance impact due to entry and exit delay requirements. Accurate modeling of these low power modes and entry and exit control is crucial to analyze the trade-offs across controller configurations and workloads with varied memory access characteristics. To address this, we integrate the power-down modes into the DRAM controller model in the open-source simulator gem5. This is the first publicly available full-system simulator with DRAM power-down modes, providing the research community a tool for DRAM power analysis for a breadth of use cases. We validate the power-down functionality with sweep tests, which trigger defined memory access characteristics. We further evaluate the model with real HPC workloads, illustrating the value of integrating low power functionality into a full system simulator.


title: Fog Massive MIMO: A User-Centric Seamless Hot-Spot Architecture

abstract: The decoupling of data and control planes, as proposed for 5G networks, will enable the efficient implementation of multitier networks where user equipment (UE) nodes obtain coverage and connectivity through the top-tier macro-cells, and, at the same time, achieve high-throughput low-latency communication through lower tiers in the hierarchy. This paper considers a new architecture for such lower tiers, dubbed fog massive MIMO, where the UEs are able to establish high-throughput low-latency data links in a seamless and opportunistic manner, as they travel through a dense fog of high-capacity wireless infrastructure nodes, referred to as remote radio heads (RRHs). Traditional handover mechanisms in dense multicell networks inherently give rise to frequent handovers and pilot sequence re-assignments, incurring, as a result, excessive protocol overhead and significant latency. In the proposed fog massive MIMO architecture, UEs seamlessly and implicitly associate themselves to the most convenient RRHs in a completely autonomous manner. Each UE makes use of a unique uplink pilot sequence, and pilot contamination is mitigated by a novel coded "on-the-fly" pilot contamination control mechanism. We analyze the spectral efficiency and the outage probability of the proposed architecture via stochastic geometry, using some recent results on unique coverage in Boolean models, and provide a detailed comparison with respect to an idealized baseline massive MIMO cellular system, that neglects protocol overhead and latency due to explicit user-cell association. Our analysis, supported by extensive system simulation, reveals that there exists a "sweet spot" of the per-pilot user load (number of users per pilot), such that the proposed system achieves spectral efficiency close to that of an ideal cellular system with the minimum distance user-base station association and no pilot/handover overhead.


title: Chaos in Continuously Monitored Quantum Systems: An Optimal Path Approach

abstract: We predict that continuously monitored quantum dynamics can be chaotic. The optimal paths between past and future boundary conditions can diverge exponentially in time when there is time-dependent evolution and continuous weak monitoring. Optimal paths are defined by extremizing the global probability density to move between two boundary conditions. We investigate the onset of chaos in pure-state qubit systems with optimal paths generated by a periodic Hamiltonian. Specifically, chaotic quantum dynamics are demonstrated in a scheme where two non-commuting observables of a qubit are continuously monitored, and one measurement strength is periodically modulated. The optimal quantum paths in this example bear similarities to the trajectories of the kicked rotor, or standard map, which is a paradigmatic example of classical chaos. We emphasize connections with the concept of resonance between integrable optimal paths and weak periodic perturbations, as well as our previous work on "multipaths", and connect the optimal path chaos to instabilities in the underlying quantum trajectories.


title: IntPhys: A Framework and Benchmark for Visual Intuitive Physics Reasoning

abstract: In order to reach human performance on complexvisual tasks, artificial systems need to incorporate a sig-nificant amount of understanding of the world in termsof macroscopic objects, movements, forces, etc. Inspiredby work on intuitive physics in infants, we propose anevaluation benchmark which diagnoses how much a givensystem understands about physics by testing whether itcan tell apart well matched videos of possible versusimpossible events constructed with a game engine. Thetest requires systems to compute a physical plausibilityscore over an entire video. It is free of bias and cantest a range of basic physical reasoning concepts. Wethen describe two Deep Neural Networks systems aimedat learning intuitive physics in an unsupervised way,using only physically possible videos. The systems aretrained with a future semantic mask prediction objectiveand tested on the possible versus impossible discrimi-nation task. The analysis of their results compared tohuman data gives novel insights in the potentials andlimitations of next frame prediction architectures.


title: Online Learning: Sufficient Statistics and the Burkholder Method

abstract: We uncover a fairly general principle in online learning: If regret can be (approximately) expressed as a function of certain "sufficient statistics" for the data sequence, then there exists a special Burkholder function that 1) can be used algorithmically to achieve the regret bound and 2) only depends on these sufficient statistics, not the entire data sequence, so that the online strategy is only required to keep the sufficient statistics in memory. This characterization is achieved by bringing the full power of the Burkholder Method --- originally developed for certifying probabilistic martingale inequalities --- to bear on the online learning setting. To demonstrate the scope and effectiveness of the Burkholder method, we develop a novel online strategy for matrix prediction that attains a regret bound corresponding to the variance term in matrix concentration inequalities. We also present a linear-time/space prediction strategy for parameter free supervised learning with linear classes and general smooth norms.


title: Coherent multi-heterodyne spectroscopy using acousto-optic frequency combs

abstract: We propose and characterize experimentally a new source of optical frequency combs for performing multi-heterodyne spectrometry. This comb modality is based on a frequency shifting loop seeded with a CW monochromatic laser. The comb lines are generated by successive passes of the CW laser through an acousto-optic frequency shifter. We report the generation of frequency combs with more than 1500 mutually coherent lines, without resorting to non-linear broadening phenomena or external electronic modulation. The comb line spacing is easily reconfigurable from tens of MHz down to the kHz region. We first use a single acousto-optic frequency comb to conduct self-heterodyne interferometry with a high frequency resolution (500 kHz). By increasing the line spacing to 80 MHz, we demonstrate molecular spectroscopy on the sub-millisecond time scale. In order to reduce the detection bandwidth, we subsequently implement an acousto-optic dual-comb spectrometer with the aid of two mutually coherent frequency shifting loops. In each architecture, the potentiality of acousto-optic frequency combs for spectroscopy is validated by spectral measurements of hydrogen cyanide in the near-infrared region.


title: Geometric decomposition of the conformation tensor in viscoelastic turbulence

abstract: This work introduces a mathematical approach to analysing the polymer dynamics in turbulent viscoelastic flows that uses a new geometric decomposition of the conformation tensor, along with associated scalar measures of the polymer fluctuations. The approach circumvents an inherent difficulty in traditional Reynolds decompositions of the conformation tensor: the fluctuating tensor fields are not positive-definite and so do not retain the physical meaning of the tensor. The geometric decomposition of the conformation tensor yields both mean and fluctuating tensor fields that are positive-definite. The fluctuating tensor in the present decomposition has a clear physical interpretation as a polymer deformation relative to the mean configuration. Scalar measures of this fluctuating conformation tensor are developed based on the non-Euclidean geometry of the set of positive-definite tensors. Drag-reduced viscoelastic turbulent channel flow is then used an example case study. The conformation tensor field, obtained using direct numerical simulations, is analysed using the proposed framework.


title: Selective Area Superconductor Epitaxy to Ballistic Semiconductor Nanowires

abstract: Semiconductor nanowires such as InAs and InSb are promising materials for studying Majorana zero-modes and demonstrating non-Abelian particle exchange relevant for topological quantum computing. While evidence for Majorana bound states in nanowires has been shown, the majority of these experiments are marked by significant disorder. In particular, the interfacial inhomogeneity between the superconductor and nanowire is strongly believed to be the main culprit for disorder and the resulting soft superconducting gap ubiquitous in tunneling studies of hybrid semiconductor-superconductor systems. Additionally, a lack of ballistic transport in nanowire systems can create bound states that mimic Majorana signatures. We resolve these problems through the development of selective-area epitaxy of Al to InSb nanowires, a technique applicable to other nanowires and superconductors. Epitaxial InSb-Al devices generically possess a hard superconducting gap and demonstrate ballistic 1D superconductivity and near perfect transmission of supercurrents in the single mode regime, requisites for engineering and controlling 1D topological superconductivity. Additionally, we demonstrate that epitaxial InSb-Al superconducting island devices, the building blocks for Majorana based quantum computing applications, prepared using selective area epitaxy can achieve micron scale ballistic 1D transport. Our results pave the way for the development of networks of ballistic superconducting electronics for quantum device applications.


title: Measurement driven single temperature engine

abstract: A four stroke quantum engine which alternately interacts with a measurement apparatus and a single heat bath is discussed in detail with respect to the average work and heat as well as to the fluctuations of work and heat. The efficiency and the reliability of such an engine with a harmonic oscillator as working substance are analyzed under different conditions such as different speeds of the work strokes, different temperatures of the heat bath and various strengths of the energy supplying measurement. For imperfect thermalization strokes of finite duration also the power of the engine is analyzed. A comparison with a two-temperature Otto engine is provided in the particular case of adiabatic work and ideal thermalization strokes.


title: Long term availability of raw experimental data in experimental fracture mechanics

abstract: Experimental data availability is a cornerstone for reproducibility in experimental fracture mechanics, which is crucial to the scientific method. This short communication focuses on the accessibility and long term availability of raw experimental data. The corresponding authors of the eleven most cited papers, related to experimental fracture mechanics, for every year from 2000 up to 2016, were kindly asked about the availability of the raw experimental data associated with each publication. For the 187 e-mails sent: 22.46% resulted in outdated contact information, 57.75% of the authors did received our request and did not reply, and 19.79 replied to our request. The availability of data is generally low with only $11$ available data sets (5.9%). The authors identified two main issues for the lacking availability of raw experimental data. First, the ability to retrieve data is strongly attached to the the possibility to contact the corresponding author. This study suggests that institutional e-mail addresses are insufficient means for obtaining experimental data sets. Second, lack of experimental data is also due that submission and publication does not require to make the raw experimental data available. The following solutions are proposed: (1) Requirement of unique identifiers, like ORCID or ResearcherID, to detach the author(s) from their institutional e-mail address, (2) Provide DOIs, like Zenodo or Dataverse, to make raw experimental data citable, and (3) grant providing organizations should ensure that experimental data by public funded projects is available to the public.


title: Resummation for (boosted) top-quark pair production at NNLO+NNLL' in QCD

abstract: We construct predictions for top quark pair differential distributions at hadron colliders that combine state-of-the-art NNLO QCD calculations with double resummation at NNLL' accuracy of threshold logarithms arising from soft gluon emissions and of small mass logarithms. This is the first time a resummed calculation at full NNLO+NNLL' accuracy in QCD for a process with non-trivial color structure has been completed at the differential level. Of main interest to us is the stability of the M_{tt} and top-quark p_T distributions in the boosted regime where fixed order calculations may become strongly dependent on the choice of dynamic scales. With the help of numeric and analytic arguments we confirm that the choice for the factorization and renormalization scales advocated recently by some of the authors is indeed optimal. We further derive a set of optimized kinematics-dependent scales for the matching functions which appear in the resummed calculations. Our NNLO+NNLL' prediction for the top-pair invariant mass is significantly less sensitive to the choice of factorization scale than the fixed order prediction, even at NNLO. Notably, the resummed and fixed order calculations are in nearly perfect agreement with each other in the full M_{tt} range when the optimal dynamic scale is used. For the top-quark p_T distribution the resummation performed here has less of an impact and instead we find that upgrading the matching with fixed-order from NLO+NNLL' to NNLO+NNLL' to be an important effect, a point to be kept in mind when using NLO-based Monte Carlo event generators to calculate this distribution.


title: Dynamic Filtering with Large Sampling Field for ConvNets

abstract: We propose a dynamic filtering strategy with large sampling field for ConvNets (LS-DFN), where the position-specific kernels learn from not only the identical position but also multiple sampled neighbor regions. During sampling, residual learning is introduced to ease training and an attention mechanism is applied to fuse features from different samples. Such multiple samples enlarge the kernels' receptive fields significantly without requiring more parameters. While LS-DFN inherits the advantages of DFN, namely avoiding feature map blurring by position-wise kernels while keeping translation invariance, it also efficiently alleviates the overfitting issue caused by much more parameters than normal CNNs. Our model is efficient and can be trained end-to-end via standard back-propagation. We demonstrate the merits of our LS-DFN on both sparse and dense prediction tasks involving object detection, semantic segmentation, and flow estimation. Our results show LS-DFN enjoys stronger recognition abilities in object detection and semantic segmentation tasks on VOC benchmark and sharper responses in flow estimation on FlyingChairs dataset compared to strong baselines.


title: Efficient treatment of bilinear forms in global optimization

abstract: We efficiently treat bilinear forms in the context of global optimization, by applying McCormick convexification and by extending an approach of Saxena, Bonami and Lee for symmetric quadratic forms to bilinear forms. A key application of our work is in treating "structural convexity" in a symmetric quadratic form.


title: Mechanisms for Primordial Black Hole Production in String Theory

abstract: We consider mechanisms for producing a significant population of primordial black holes (PBHs) within string inspired single field models of inflation. The production of PBHs requires a large amplification in the power spectrum of curvature perturbations between scales associated with CMB and PBH formation. In principle, this can be achieved by temporarily breaking the slow-roll conditions during inflation. In this work, we identify two string setups that can realise this process. In string axion models of inflation, subleading non-perturbative effects can superimpose steep cliffs and gentle plateaus onto the leading axion potential. The cliffs can momentarily violate the slow-roll conditions, and the plateaus can lead to phases of ultra slow-roll inflation. We thus achieve a string motivated model which both matches the Planck observations at CMB scales and produces a population of light PBHs, which can account for an order one fraction of dark matter. In DBI models of inflation, a sharp increase in the speed of sound sourced by a steep downward step in the warp factor can drive the amplification. In this scenario, discovery of PBHs could indicate non-trivial dynamics in the bulk, such as flux-antibrane annihilation at the tip of a warped throat.


title: Diagonal reduction of matrices over commutative semihereditary Bezout rings

abstract: It is proven that every commutative semihereditary Bezout ring in which any regular element is Gelfand (adequate), is an elementary divisor ring.


title: On Bethe vectors in $\mathfrak{gl}_3$-invariant integrable models

abstract: We consider quantum integrable models solvable by the nested algebraic Bethe ansatz and possessing $\mathfrak{gl}_3$-invariant $R$-matrix. We study a new recently proposed approach to construct on-shell Bethe vectors of these models. We prove that the vectors constructed by this method are semi-on-shell Bethe vectors for arbitrary values of Bethe parameters. They thus do become on-shell vectors provided the system of Bethe equations is fulfilled.

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


title: Towards observation of three-nucleon short-range correlations in high Q^2 A(e, e')X reactions

abstract: We present a detailed study of kinematical and dynamical conditions necessary for probing highly elusive three-nucleon short range correlations~(3N-SRCs) in nuclei through inclusive electron scattering. The kinematic requirements that should be satisfied in order to isolate 3N-SRCs in inclusive processes are derived. We demonstrate that a sequence of two short-range NN interactions represents the main mechanism for 3N-SRCs in inclusive processes. Within this mechanism we predict a quadratic dependence of the inclusive cross section ratios of nuclei to $^3$He in the 3N-SRC region to the same ratio measured in 2N-SRC domain. The extended analysis of the available data satisfying the necessary 3N-SRC kinematical conditions is presented. This analysis provides tantalizing signatures of scaling associated with the onset of 3N-SRC dominance. The same data are also consistent with the prediction of the quadratic relation between the ratios measured in the 3N and 2N-SRC regions for nuclei ranging $4 \le A \le 197$. This agreement made it possible to extract $a_3(A)$, the probability of 3N-SRCs relative to the A=3 nucleus. We find $a_3(A)$ to be significantly larger in magnitude than the analogous parameter, $a_2(A)$, for 2N-SRCs.


title: The Effects of Bandpass Variations on Foreground Removal Forecasts for Future CMB Experiments

abstract: Time-dependent and systematic variations in band gain and central frequencies of instruments used to study the Cosmic Microwave Background are important factors in the data-to-map analysis pipeline. If not properly characterized, they could limit the ability of next-generation experiments to remove astrophysical foreground contamination. Uncertainties include the instrument detector band, which could systematically change across the focal plane, as well as the calibration of the instrument used to measure the bands. A potentially major effect is time-dependent gain and band uncertainties caused by atmospheric fluctuations. More specifically, changes in atmospheric conditions lead to frequency-dependent changes in the atmospheric transmission which, in turn, leads to variations in the effective gain and central frequency of the instrument's bandpass. Using atmospheric modeling software and ACTPol bandpasses, we simulate the expected variations in band gain and central frequency for 20, 40, 90, 150, and 240 GHz bands as a function of precipitable water vapor, observing angle, and ground temperature. Combining these effects enables us to set maximum and minimum limits on the expected uncertainties in band gain and central frequency over the course of a full observing season. We then introduce the uncertainties to parametric maximum-likelihood component separation methods on simulated CMB maps to forecast foreground removal performance and likelihoods on the tensor-to-scalar ratio r. We conclude that to measure a $\sigma(r=0)$ ~ $10^{-3}$ with a bias on the recovered $r$ under control, the limit on the uncertainty in the relative gain and central frequency of the bandpass must be <2% and <1%, respectively. We also comment on the possibility of self-calibrating bandpass uncertainties.


title: Smooth deformations of singular contractions of class VII surfaces

abstract: We consider normal compact surfaces $Y$ obtained from a minimal class VII surface $X$ by contraction of a cycle $C$ of $r$ rational curves with $C^2<0$. Our main result states that, if the obtained cusp is smoothable, then $Y$ is globally smoothable. The proof is based on a vanishing theorem for $H^2(\Theta_Y)$. If $r<b_2(X)$ any smooth small deformation of $Y$ is rational, and if $r=b_2(X)$ (i.e. when $X$ is a half-Inoue surface) any smooth small deformation of $Y$ is an Enriques surface. The condition "the cusp is smoothable" in our main theorem can be checked in terms of the intersection numbers of the cycle, using the Looijenga conjecture (which has recently become a theorem). Therefore this is a "decidable" condition. We prove that this condition is always satisfied if $r<b_2(X)\leq 11$. Therefore the singular surface $Y$ obtained by contracting a cycle $C$ of $r$ rational curves in a minimal class VII surface $X$ with $r<b_2(X)\leq 11$ is always smoothable by rational surfaces. The statement holds even for unknown class VII surfaces.


title: Spin Temperature and Density of Cold and Warm HI in the Galactic Disk - Hidden HI -

abstract: We present a method to determine the spin temperature $T_{\rm S}$ and volume density $n$ of HI gas simultaneously along the tangent-point circle of galactic rotation in the Milky Way by using the least-$\chi^2$ method. The best-fit $T_{\rm S}$ is shown to range either in $T_{\rm S}=100-120$ K or in $1000-3000$ K, indicating that the gas is either in cold HI phase with high density and large optical depth, or in warm HI with low density and small optical depth. Averaged values at $3\le R \le 8$ kpc are obtained to be $T_{\rm S}=106.7 \pm 16.0$ K and $n=1.53\pm 0.86$ H cm$^{-3}$ for cold HI, and $1720 \pm 1060$ K and $0.38 \pm 0.10$ H cm$^{-3}$ for warm HI, where $R=8\ |\sin \ l|$ kpc is the galacto-centric distance along the tangent-point circle. The cold HI appears in spiral arms and rings, whereas warm HI in the inter-arm regions. The cold HI is denser by a factor of $\sim 4$ than warm HI. The present analysis has revealed the hidden HI mass in cold and optically thick phase in the galactic disk. The total HI mass inside the solar circle is shown to be greater by a factor of $2 - 2.5$ than the current estimation by optically thin assumption.


title: The Menger and projective Menger properties of function spaces with the set-open topology

abstract: For a Tychonoff space $X$ and a family $\lambda$ of subsets of $X$, we denote by $C_{\lambda}(X)$ the space of all real-valued continuous functions on $X$ with the set-open topology. In this paper, we study the Menger and projective Menger properties of a Hausdorff space $C_{\lambda}(X)$. Our main results state that if $\lambda$ is a $\pi$-network of $X$ then (1) $C_{\lambda}(X)$ is Menger space if and only if it is $\sigma$-compact, if $Y$ is a dense subset of $X$ then (2) $C_{p}(Y\vert X)$ is projective Menger space if and only if it is $\sigma$-pseudocompact.


title: Domain Adaptation with Randomized Expectation Maximization

abstract: Domain adaptation (DA) is the task of classifying an unlabeled dataset (target) using a labeled dataset (source) from a related domain. The majority of successful DA methods try to directly match the distributions of the source and target data by transforming the feature space. Despite their success, state of the art methods based on this approach are either involved or unable to directly scale to data with many features. This article shows that domain adaptation can be successfully performed by using a very simple randomized expectation maximization (EM) method. We consider two instances of the method, which involve logistic regression and support vector machine, respectively. The underlying assumption of the proposed method is the existence of a good single linear classifier for both source and target domain. The potential limitations of this assumption are alleviated by the flexibility of the method, which can directly incorporate deep features extracted from a pre-trained deep neural network. The resulting algorithm is strikingly easy to implement and apply. We test its performance on 36 real-life adaptation tasks over text and image data with diverse characteristics. The method achieves state-of-the-art results, competitive with those of involved end-to-end deep transfer-learning methods.


title: Learning Robotic Assembly from CAD

abstract: In this work, motivated by recent manufacturing trends, we investigate autonomous robotic assembly. Industrial assembly tasks require contact-rich manipulation skills, which are challenging to acquire using classical control and motion planning approaches. Consequently, robot controllers for assembly domains are presently engineered to solve a particular task, and cannot easily handle variations in the product or environment. Reinforcement learning (RL) is a promising approach for autonomously acquiring robot skills that involve contact-rich dynamics. However, RL relies on random exploration for learning a control policy, which requires many robot executions, and often gets trapped in locally suboptimal solutions. Instead, we posit that prior knowledge, when available, can improve RL performance. We exploit the fact that in modern assembly domains, geometric information about the task is readily available via the CAD design files. We propose to leverage this prior knowledge by guiding RL along a geometric motion plan, calculated using the CAD data. We show that our approach effectively improves over traditional control approaches for tracking the motion plan, and can solve assembly tasks that require high precision, even without accurate state estimation. In addition, we propose a neural network architecture that can learn to track the motion plan, and generalize the assembly controller to changes in the object positions.


title: A micro-SQUID with dispersive readout for magnetic scanning microscopy

abstract: We have designed and characterized a micro-SQUID with dispersive readout for use in low temperature scanning probe microscopy systems. The design features a capacitively shunted RF SQUID with a tunable resonance frequency from 5 to 12 GHz, micrometer spatial resolution and integrated superconducting coils for local application of magnetic fields. The SQUID is operated as a nonlinear oscillator with a flux- and power-dependent resonance frequency. Measurements for device characterization and noise benchmarking were carried out at 4 K. The measured flux noise above 10 kHz at 4 K is $80\;\mathrm{n}\Phi_0 \mathrm{Hz}^{-1/2}$ at a bandwidth of 200 MHz. Estimations suggest that one can benefit from parametric gain based on inherent nonlinearity of the Josephson junction and reduce the flux noise to $30\;\mathrm{n}\Phi_0\mathrm{Hz}^{-1/2}$ at 100 mK, which corresponds to $7.5\;\mathrm{\mu_B Hz^{-1/2}}$ for a magnetic moment located at the center of the pickup loop.


title: The Applegate mechanism in Post-Common-Envelope Binaries: Investigating the role of rotation

abstract: Eclipsing time variations (ETVs) are observed in many close binary systems. In particular, for several post-common-envelope binaries (PCEBs) that consist of a white dwarf and a main sequence star, the O-C diagram suggests that real or apparent orbital period variations are driven by Jupiter-mass planets or as a result of magnetic activity, the so-called Applegate mechanism. The latter explains orbital period variations as a result of changes in the stellar quadrupole moment due to magnetic activity. We explore the feasibility of driving ETVs via the Applegate mechanism for a sample of PCEB systems, including a range of different rotations. Using the MESA code we evolve 12 stars with different masses and rotation rates. We apply a simple dynamo model to their radial profiles to investigate on which scale the predicted activity cycle matches the observed modulation period, and quantify the uncertainty, and further calculate the required energies to drive que Applegate mechanism. We show that the Applegate mechanism is energetically feasible in 5 PCEB systems, and note that these are the systems with the highest rotation rate compared to the critical rotation rate of the main-sequence star. The results suggest that the ratio of physical to critical rotation in the main sequence star is an important indicator for the feasibility of Applegate's mechanism, but exploring larger samples will be necessary to probe this hypothesis.


title: Photoproduction of light vector mesons in Xe-Xe ultraperipheral collisions at the LHC and the nuclear density of Xe-129

abstract: We make predictions for cross sections of $\rho$ and $\phi$ vector meson photoproduction in ultraperipheral Xe-Xe collisions at $\sqrt{s_{NN}}=5.44$ TeV. Analyzing the momentum transfer distribution of $\rho$ mesons in this process, we explore the feasibility of extracting the nuclear density of $^{129}$Xe, which is needed in searches for dark matter with Xenon-based detectors.


title: Adaptive Greedy Algorithms for Stochastic Set Cover Problems

abstract: We study adaptive greedy algorithms for the problems of stochastic set cover with perfect and imperfect coverages. In stochastic set cover with perfect coverage, we are given a set of items and a ground set B. Evaluating an item reveals its state which is a random subset of B drawn from the state distribution of the item. Every element in B is assumed to be present in the state of some item with probability 1. For this problem, we show that the adaptive greedy algorithm has an approximation ratio of H(|B|), the |B|th Harmonic number. In stochastic set cover with imperfect coverage, an element in the ground set need not be present in the state of any item. We show a reduction from this problem to the former problem; the adaptive greedy algorithm for the reduced instance has an approxiation ratio of H(|E|), where E is the set of pairs (F, e) such that the state of item F contains e with positive probability.


title: AllenNLP: A Deep Semantic Natural Language Processing Platform

abstract: This paper describes AllenNLP, a platform for research on deep learning methods in natural language understanding. AllenNLP is designed to support researchers who want to build novel language understanding models quickly and easily. It is built on top of PyTorch, allowing for dynamic computation graphs, and provides (1) a flexible data API that handles intelligent batching and padding, (2) high-level abstractions for common operations in working with text, and (3) a modular and extensible experiment framework that makes doing good science easy. It also includes reference implementations of high quality approaches for both core semantic problems (e.g. semantic role labeling (Palmer et al., 2005)) and language understanding applications (e.g. machine comprehension (Rajpurkar et al., 2016)). AllenNLP is an ongoing open-source effort maintained by engineers and researchers at the Allen Institute for Artificial Intelligence.


title: An ADMM-based Coordination and Control Strategy for PV and Storage to Dispatch Stochastic Prosumers: Theory and Experimental Validation

abstract: This paper describes a two-layer control and coordination framework for distributed energy resources. The lower layer is a real-time model predictive control (MPC) executed at 10 s resolution to achieve fine tuning of a given energy set-point. The upper layer is a slower MPC coordination mechanism based on distributed optimization, and solved with the alternating direction method of multipliers (ADMM) at 5 minutes resolution. It is needed to coordinate the power flow among the controllable resources such that enough power is available in real-time to achieve a pre-established energy trajectory in the long term. Although the formulation is generic, it is developed for the case of a battery system and a curtailable PV facility to dispatch stochastic prosumption according to a trajectory at 5 minutes resolution established the day before the operation. The proposed method is experimentally validated in a real-life setup to dispatch the operation of a building with rooftop PV generation (i.e., 101 kW average load, 350 kW peak demand, 82 kW peak PV generation) by controlling a 560 kWh/720 kVA battery and a 13 kW peak curtailable PV facility.


title: Local criteria for triangulation of manifolds

abstract: We present criteria for establishing a triangulation of a manifold. Given a manifold M, a simplicial complex A, and a map H from the underlying space of A to M, our criteria are presented in local coordinate charts for M, and ensure that H is a homeomorphism. These criteria do not require a differentiable structure, or even an explicit metric on M. No Delaunay property of A is assumed. The result provides a triangulation guarantee for algorithms that construct a simplicial complex by working in local coordinate patches. Because the criteria are easily verified in such a setting, they are expected to be of general use.


title: On the Dynamics of Distributed Energy Adoption: Equilibrium, Stability, and Limiting Capacity

abstract: The death spiral hypothesis in electric utility represents a positive feedback phenomenon in which a regulated utility is driven to financial instability by rising prices and declining demand. We establish conditions for the existence of death spiral and conditions of stable adoption of distributed energy resources. We show in particular that linear tariffs always induce death spiral when the fixed operating cost of the utility rises beyond a certain threshold. For two-part tariffs with connection and volumetric charges, the Ramsey pricing that optimizes myopically social welfare subject to the revenue adequacy constraint induces a stable equilibrium. The Ramsey pricing, however, inhibits renewable adoption with a high connection charge. In contrast, a two-part tariff with a small connection charge results in a stable adoption process with a higher level of renewable adoption and greater long-term total consumer surplus. Market data are used to illustrate various solar adoption scenarios.


title: Effective Field Theory description of Phantom Fields

abstract: In this work I show that a simple Field Theory on a non trivial gauge background may behave as a phantom field and contribute to an effective $w<-1$ state equation fluid contribution to cosmology.


title: V-Splines and Bayes Estimate

abstract: Smoothing splines can be thought of as the posterior mean of a Gaussian process regression in a certain limit. By constructing a reproducing kernel Hilbert space with an appropriate inner product, the Bayesian form of the V-spline is derived when the penalty term is a fixed constant instead of a function. An extension to the usual generalized cross-validation formula is utilized to find the optimal V-spline parameters.


title: Optical and Near-Infrared Radial Velocity Content of M Dwarfs: Testing Models with Barnard's Star

abstract: High precision radial velocity (RV) measurements have been central in the study of exoplanets during the last two decades, from the early discovery of hot Jupiters, to the recent mass measurements of Earth-sized planets uncovered by transit surveys. While optical radial-velocity is now a mature field, there is currently a strong effort to push the technique into the near-infrared (nIR) domain (chiefly $Y$, $J$, $H$ and $K$ band passes) to probe planetary systems around late-type stars. The combined lower mass and luminosity of M dwarfs leads to an increased reflex RV signal for planets in the habitable zone compared to Sun-like stars. The estimates on the detectability of planets rely on various instrumental characteristics, but also on a prior knowledge of the stellar spectrum. While the overall properties of M dwarf spectra have been extensively tested against observations, the same is not true for their detailed line profiles, which leads to significant uncertainties when converting a given signal-to-noise ratio to a corresponding RV precision as attainable on a given spectrograph. By combining archival CRIRES and HARPS data with ESPaDOnS data of Barnard's star, we show that state-of-the-art atmosphere models over-predict the $Y$ and $J$-band RV content by more than a factor of $\sim$$2$, while under-predicting the $H$ and $K$-band content by half.


title: The bunkbed conjecture on the complete graph

abstract: The bunkbed conjecture was first posed by Kasteleyn. If $G=(V,E)$ is a finite graph and $H$ some subset of $V$, then the bunkbed of the pair $(G,H)$ is the graph $G\times\{1,2\}$ plus $|H|$ extra edges to connect for every $v\in H$ the vertices $(v,1)$ and $(v,2)$. The conjecture asserts that $(v,1)$ is more likely to connect with $(w,1)$ than with $(w,2)$ in the independent bond percolation model for any $v,w\in V$. This is intuitive because $(v,1)$ is in some sense closer to $(w,1)$ than it is to $(w,2)$. The conjecture has however resisted several attempts of proof. This paper settles the conjecture in the case of a constant percolation parameter and $G$ the complete graph.


title: Identifying Relevant Information Cues for Vulnerability Assessment Using CVSS

abstract: The assessment of new vulnerabilities is an activity that accounts for information from several data sources and produces a `severity' score for the vulnerability. The Common Vulnerability Scoring System (\CVSS) is the reference standard for this assessment. Yet, no guidance currently exists on \emph{which information} aids a correct assessment and should therefore be considered. In this paper we address this problem by evaluating which information cues increase (or decrease) assessment accuracy. We devise a block design experiment with 67 software engineering students with varying vulnerability information and measure scoring accuracy under different information sets. We find that baseline vulnerability descriptions provided by standard vulnerability sources provide only part of the information needed to achieve an accurate vulnerability assessment. Further, we find that additional information on \texttt{assets}, \texttt{attacks}, and \texttt{vulnerability type} contributes in increasing the accuracy of the assessment; conversely, information on \texttt{known threats} misleads the assessor and decreases assessment accuracy and should be avoided when assessing vulnerabilities. These results go in the direction of formalizing the vulnerability communication to, for example, fully automate security assessments.


title: An Online Algorithm for Power-proportional Data Centers with Switching Cost

abstract: Recent studies have shown that power-proportional data centers can save energy cost by dynamically "right-sizing" the data centers based on real-time workload. More servers are activated when the workload increases while some servers can be put into the sleep mode during periods of low load. In this paper, we revisit the dynamic right-sizing problem for heterogeneous data centers with various operational cost and switching cost. We propose a new online algorithm based on a regularization technique, which achieves a better competitive ratio compared to the state-of-the-art greedy algorithm. We further introduce a switching cost offset into the model and extend our algorithm to this new setting. Simulations based on real workload and renewable energy traces show that our algorithms outperform the greedy algorithm in both settings.


title: Quantifying resilience to recurrent ecosystem disturbances using flow-kick dynamics

abstract: Shifting ecosystem disturbance patterns due to climate change (e.g. storms, droughts, wildfires) or direct human interference (e.g. harvests, nutrient loading) highlight the importance of quantifying and strengthening the resilience of desired ecological regimes. Although existing metrics capture resilience to isolated shocks, gradual parameter changes, and continuous noise, quantifying resilience to repeated, discrete disturbances requires novel analytical tools. Here we introduce a flow-kick framework that quantifies resilience to disturbances explicitly in terms of their magnitude and frequency. We present a resilience boundary between disturbances that cause either escape from a basin of attraction or stabilization within it, and use the resilience boundary to build resilience metrics tailored to repeated, discrete disturbances. The flow-kick model suggests that the distance-to-threshold resilience metric overestimates resilience in the context of repeated disturbances. It also reveals counterintuitive triggers for regime shifts, such as increasing recovery times between disturbances, or increasing disturbance magnitude and recovery times proportionately.


title: Distributed Model Predictive Control for Linear Systems with Adaptive Terminal Sets

abstract: In this paper, we propose a distributed model predictive control (DMPC) scheme for linear time-invariant constrained systems which admit a separable structure. To exploit the merits of distributed computation algorithms, the stabilizing terminal controller, value function and invariant terminal set of the DMPC optimization problem need to respect the loosely coupled structure of the system. Although existing methods in the literature address this task, they typically decouple the synthesis of terminal controllers and value functions from the one of terminal sets. In addition, these approaches do not explicitly consider the effect of the current state of the system in the synthesis process. These limitations can lead the resulting DMPC scheme to poor performance since it may admit small or even empty terminal sets. Unlike other approaches, this paper presents a unified framework to encapsulate the synthesis of both the stabilizing terminal controller and invariant terminal set into the DMPC formulation. Conditions for Lyapunov stability and invariance are imposed in the synthesis problem in a way that allows the value function and invariant terminal set to admit the desired distributed structure. We illustrate the effectiveness of the proposed method on several examples including a benchmark spring-mass-damper problem.


title: ALLEGRO Core Degradation Study using MELCOR 2.1

abstract: ALLEGRO is an experimental high-temperature gas (He)-cooled fast reactor (GFR) under development by the European Consortium V4G4 Centre of Excellence on basis of the concept ALLEGRO presented by CEA in 2009. Its main purpose is to demonstrate the viability of GFR in pilot scale and simultaneously serve as a test bed for GFR-related technologies, above all the high-temperature resistant refractory fuel in the conventional oxide fuel driver core. Severe accident studies of the concept with the first core (oxide fuel in stainless steel claddings) are an important part of safety analyses leading to improves in reactor design and safety. This paper extends our previous study using the MELCOR code. Major difference between the previous and this work is decrease of thermal power and power density to 50 MWth and 66,6 MWth/m3 respectively. Analyses are focused on protected scenarios based on total station blackout individually aggravated by loss of primary coolant, water ingress into primary coolant or malfunction of check valves in the decay heat removal system. The results indicate the timing and extent of core degradation and provide valuable data for design of the core catcher in ALLEGRO.


title: The Webster scalar curvature and sharp upper and lower bounds for the first positive eigenvalue of the Kohn-Laplacian on real hypersurfaces

abstract: Let $(M,\theta)$ be a compact strictly pseudoconvex pseudohermitian manifold which is CR embedded into a complex space. In an earlier paper, Lin and the authors gave several sharp upper bounds for the first positive eigenvalue $\lambda_1$ of the Kohn-Laplacian $\Box_b$ on $(M,\theta)$. In the present paper, we give a sharp upper bound for $\lambda_1$, generalizing and extending some previous results. As a corollary, we obtain a Reilly-type estimate when $M$ is embedded into the standard sphere. In another direction, using a Lichnerowicz-type estimate by Chanillo, Chiu, and Yang and an explicit formula for the Webster scalar curvature, we give a lower bound for $\lambda_1$ when the pseudohermitian structure $\theta$ is volume-normalized.


title: Fermi-LAT observations of the 2017 September 10$^{th}$ solar flare

abstract: The Fermi-Large Area Telescope (LAT) detection of the X8.2 GOES class solar flare of 2017 September 10 provides for the first time observations of a long duration high-energy gamma-ray flare associated with a Ground Level Enhancement (GLE). The >100 MeV emission from this flare lasted for more than 12 hours covering both the impulsive and extended phase. We present the localization of the gamma-ray emission and find that it is consistent with the active region (AR) from which the flare occurred over a period lasting more than 6 hours contrary to what was found for the 2012 March 7 flares. The temporal variation of the proton index inferred from the gamma-ray data seems to suggest two phases in acceleration of the proton population. Based on timing arguments we interpret the second phase to be tied to the acceleration mechanism powering the GLE, believed to be particle acceleration at a coronal shock driven by the CME.


title: On Multi-Server Coded Caching in the Low Memory Regime

abstract: In this paper we determine the delivery time for a multi-server coded caching problem when the cache size of each user is small. We propose an achievable scheme based on coded cache content placement, and employ zero-forcing techniques at the content delivery phase. Surprisingly, in contrast to previous multi-server results which were proved to be order-optimal within a multiplicative factor of 2, for the low memory regime we prove that our achievable scheme is optimal. Moreover, we compare the performance of our scheme with the uncoded solution, and show our proposal improvement over the uncoded scheme. Our results also apply to Degrees-of-Freedom (DoF) analysis of Multiple-Input Single-Output Broadcast Channels (MISO-BC) with cache-enabled users, where the multiple-antenna transmitter replaces the role of multiple servers. This shows that interference management in the low memory regime needs different caching techniques compared with medium-high memory regimes discussed in previous works.


title: The puzzling two-proton decay of $^{67}$Kr

abstract: Two-proton (2$p$) radioactivity is a rare decay mode found in a few proton-unbound nuclei. The $2p$-decay lifetime and properties of emitted protons carry invaluable information on nuclear structure in the presence of low-lying proton continuum. The recently measured $2p$ decay of $^{67}$Kr [Phys. Rev. Lett. 117, 162501 (2016)] turned out to be unexpectedly fast. Since $^{67}$Kr is expected to be a deformed system, we investigate the impact of deformation effects on the $2p$ radioactivity. We apply the recently developed Gamow coupled-channel framework, which allows for a precise description of three-body systems in the presence of rotational and vibrational couplings. This is the first application of a three-body approach to a two-nucleon decay from a deformed nucleus. We show that deformation couplings significantly increase the $2p$ decay width of $^{67}$Kr; this finding explains the puzzling experimental data. The calculated angular proton-proton correlations reflect a competition between $1p$ and $2p$ decay modes in this nucleus.