title: Cosmological signatures of time-asymmetric gravity

abstract: We develop the model proposed by Cort\^es, Gomes & Smolin, to predict cosmological signatures of time-asymmetric extensions of general relativity they proposed recently. Within this class of models the equation of motion of chiral fermions is modified by a torsion term. This term leads to a dispersion law for neutrinos that associates a new time-varying energy with each particle. We find a new neutrino contribution to the Friedmann equation resulting from the torsion term in the Ashtekar connection. In this note we explore the phenomenology of this term and observational consequences for cosmological evolution. We show that constraints on the critical energy density will ordinarily render this term unobservably small, a maximum of order $10^{-25}$ of the neutrino energy density today. However, if the time-asymmetric dark energy is tuned to cancel the cosmological constant, the torsion effect may be a dark matter candidate.


title: Noise Response Data Reveal Novel Controllability Gramian for Nonlinear Network Dynamics

abstract: Control of nonlinear large-scale dynamical networks, e.g., collective behavior of agents interacting via a scale-free connection topology, is a central problem in many scientific and engineering fields. For the linear version of this problem, the so-called controllability Gramian has played an important role to quantify how effectively the dynamical states are reachable by a suitable driving input. In this paper, we first extend the notion of the controllability Gramian to nonlinear dynamics in terms of the Gibbs distribution. Next, we show that, when the networks are open to environmental noise, the newly defined Gramian is equal to the covariance matrix associated with randomly excited, but uncontrolled, dynamical state trajectories. This fact theoretically justifies a simple Monte Carlo simulation that can extract effectively controllable subdynamics in nonlinear complex networks. In addition, the result provides a novel insight into the relationship between controllability and statistical mechanics.


title: Sterile neutrino portal to Dark Matter I: The $U(1)_{B-L}$ case

abstract: In this paper we explore the possibility that the sterile neutrino and Dark Matter sectors in the Universe have a common origin. We study the consequences of this assumption in the simple case of coupling the dark sector to the Standard Model via a global $U(1)_{B-L}$, broken down spontaneously by a dark scalar. This dark scalar provides masses to the dark fermions and communicates with the Higgs via a Higgs portal coupling. We find an interesting interplay between Dark Matter annihilation to dark scalars - the CP-even that mixes with the Higgs and the CP-odd which becomes a Goldstone boson, the Majoron - and heavy neutrinos, as well as collider probes via the coupling to the Higgs. Moreover, Dark Matter annihilation into sterile neutrinos and its subsequent decay to gauge bosons and quarks, charged leptons or neutrinos lead to indirect detection signatures which are close to current bounds on the gamma ray flux from the galactic center and dwarf galaxies.


title: The Relation Between [OIII]/H$\beta$ and Specific Star Formation Rate in Galaxies at $z \sim 2$

abstract: Recent surveys have identified a seemingly ubiquitous population of galaxies with elevated [OIII]/H$\beta$ emission line ratios at $z > 1$, though the nature of this phenomenon continues to be debated. The [OIII]/H$\beta$ line ratio is of interest because it is a main component of the standard diagnostic tools used to differentiate between active galactic nuclei (AGN) and star-forming galaxies, as well as the gas-phase metallicity indicators $O_{23}$ and $R_{23}$. Here, we investigate the primary driver of increased [OIII]/H$\beta$ ratios by median-stacking rest-frame optical spectra for a sample of star-forming galaxies in the 3D-HST survey in the redshift range $z\sim1.4-2.2$. Using $N = 4220$ star-forming galaxies, we stack the data in bins of mass and specific star formation rates (sSFR) respectively. After accounting for stellar Balmer absorption, we measure [OIII]$\lambda5007$\AA/H$\beta$ down to $\mathrm{M} \sim 10^{9.2} \ \mathrm{M_\odot}$ and sSFR $\sim 10^{-9.6} \ \mathrm{yr}^{-1}$, more than an order of magnitude lower than previous work at similar redshifts. We find an offset of $0.59\pm0.05$ dex between the median ratios at $z\sim2$ and $z\sim0$ at fixed stellar mass, in agreement with existing studies. However, with respect to sSFR, the $z \sim 2$ stacks all lie within 1$\sigma$ of the median SDSS ratios, with an average offset of only $-0.06\pm 0.05$. We find that the excitation properties of galaxies are tightly correlated with their sSFR at both $z\sim2$ and $z\sim0$, with a relation that appears to be roughly constant over the last 10 Gyr of cosmic time.


title: Bimodal entanglement entropy distribution in the many-body localization transition

abstract: We introduce the cut averaged entanglement entropy in disordered periodic spin chains and prove it to be a concave function of subsystem size for individual eigenstates. This allows us to identify the entanglement scaling as a function of subsystem size for individual states in inhomogeneous systems. Using this quantity, we probe the critical region between the many-body localized (MBL) and ergodic phases in finite systems. In the middle of the spectrum, we show evidence for bimodality of the entanglement distribution in the MBL critical region, finding both volume law and area law eigenstates over disorder realizations as well as within \emph{single disorder realizations}. The disorder averaged entanglement entropy in this region then scales as a volume law with a coefficient below its thermal value. We discover in the critical region, as we approach the thermodynamic limit, that the cut averaged entanglement entropy density falls on a one-parameter family of curves. Finally, we also show that without averaging over cuts the slope of the entanglement entropy \vs subsystem size can be negative at intermediate and strong disorder, caused by rare localized regions in the system.


title: Minimizing Regret on Reflexive Banach Spaces and Learning Nash Equilibria in Continuous Zero-Sum Games

abstract: We study a general version of the adversarial online learning problem. We are given a decision set $\mathcal{X}$ in a reflexive Banach space $X$ and a sequence of reward vectors in the dual space of $X$. At each iteration, we choose an action from $\mathcal{X}$, based on the observed sequence of previous rewards. Our goal is to minimize regret, defined as the gap between the realized reward and the reward of the best fixed action in hindsight. Using results from infinite dimensional convex analysis, we generalize the method of Dual Averaging (or Follow the Regularized Leader) to our setting and obtain general upper bounds on the worst-case regret that subsume a wide range of results from the literature. Under the assumption of uniformly continuous rewards, we obtain explicit anytime regret bounds in a setting where the decision set is the set of probability distributions on a compact metric space $S$ whose Radon-Nikodym derivatives are elements of $L^p(S)$ for some $p > 1$. Importantly, we make no convexity assumptions on either the set $S$ or the reward functions. We also prove a general lower bound on the worst-case regret for any online algorithm. We then apply these results to the problem of learning in repeated continuous two-player zero-sum games, in which players' strategy sets are compact metric spaces. In doing so, we first prove that if both players play a Hannan-consistent strategy, then with probability 1 the empirical distributions of play weakly converge to the set of Nash equilibria of the game. We then show that, under mild assumptions, Dual Averaging on the (infinite-dimensional) space of probability distributions indeed achieves Hannan-consistency. Finally, we illustrate our results through numerical examples.


title: Directly comparing GW150914 with numerical solutions of Einstein's equations for binary black hole coalescence

abstract: We compare GW150914 directly to simulations of coalescing binary black holes in full general relativity, accounting for all the spin-weighted quadrupolar modes, and separately accounting for all the quadrupolar and octopolar modes. Consistent with the posterior distributions reported in LVC_PE[1] (at 90% confidence), we find the data are compatible with a wide range of nonprecessing and precessing simulations. Followup simulations performed using previously-estimated binary parameters most resemble the data. Comparisons including only the quadrupolar modes constrain the total redshifted mass Mz \in [64 - 82M_\odot], mass ratio q = m2/m1 \in [0.6,1], and effective aligned spin \chi_eff \in [-0.3, 0.2], where \chi_{eff} = (S1/m1 + S2/m2) \cdot\hat{L} /M. Including both quadrupolar and octopolar modes, we find the mass ratio is even more tightly constrained. Simulations with extreme mass ratios and effective spins are highly inconsistent with the data, at any mass. Several nonprecessing and precessing simulations with similar mass ratio and \chi_{eff} are consistent with the data. Though correlated, the components' spins (both in magnitude and directions) are not significantly constrained by the data. For nonprecessing binaries, interpolating between simulations, we reconstruct a posterior distribution consistent with previous results. The final black hole's redshifted mass is consistent with Mf,z between 64.0 - 73.5M_\odot and the final black hole's dimensionless spin parameter is consistent with af = 0.62 - 0.73. As our approach invokes no intermediate approximations to general relativity and can strongly reject binaries whose radiation is inconsistent with the data, our analysis provides a valuable complement to LVC_PE[1].


title: NH$_3$(3,3) and CH$_3$OH near Supernova Remnants: GBT and VLA Observations

abstract: We report on Green Bank Telescope 23.87 GHz NH$_3$(3,3), emission observations in five supernova remnants interacting with molecular clouds (G1.4$-$0.1, IC443, W44, W51C, and G5.7$-$0.0). The observations show a clumpy gas density distribution, and in most cases the narrow line widths of $\sim3-4$\,km\,s$^{-1}$ are suggestive of maser emission. Very Large Array observations reveal 36~GHz and/or 44~GHz CH$_3$OH, maser emission in a majority (72\%) of the NH$_3$, peak positions towards three of these SNRs. This good positional correlation is in agreement with the high densities required for the excitation of each line. Through these observations we have shown that CH$_3$OH, and NH$_3$, maser emission can be used as indicators of high density clumps of gas shocked by supernova remnants, and provide density estimates thereof. Modeling of the optical depth of the NH$_3$(3,3) emission is compared to that of CH$_3$OH, constraining the densities of the clumps to a typical density of the order of $10^5$~cm$^{-3}$ for cospatial masers. Regions of gas with this density are found to exist in the post-shocked gas quite close to the SNR shock front, and may be associated with sites where cosmic rays produce gamma-ray emission via neutral pion decay.


title: Transiting exoplanet candidates from K2 Campaigns 5 and 6

abstract: We introduce a new transit search and vetting pipeline for observations from the K2 mission, and present the candidate transiting planets identified by this pipeline out of the targets in Campaigns 5 and 6. Our pipeline uses the Gaussian Process-based K2SC code to correct for the K2 pointing systematics and simultaneously model stellar variability. The systematics-corrected, variability-detrended light curves are searched for transits with the Box Least Squares method, and a period-dependent detection threshold is used to generate a preliminary candidate list. Two or three individuals vet each candidate manually to produce the final candidate list, using a set of automatically-generated transit fits and assorted diagnostic tests to inform the vetting. We detect 147 single-planet system candidates and 5 multi-planet systems, independently recovering the previously-published hot~Jupiters EPIC 212110888b, WASP-55b (EPIC 212300977b) and Qatar-2b (EPIC 212756297b). We also report the outcome of reconnaissance spectroscopy carried out for all candidates with Kepler magnitude Kp < 13, identifying 12 targets as likely false positives. We compare our results to those of other K2 transit search pipelines, noting that ours performs particularly well for variable and/or active stars, but that the results are very similar overall. All the light curves and code used in the transit search and vetting process are publicly available, as are the follow-up spectra.


title: Gaussian process emulators for computer experiments with inequality constraints

abstract: Physical phenomena are observed in many fields (sciences and engineering) and are often studied by time-consuming computer codes. These codes are analyzed with statistical models, often called emulators. In many situations, the physical system (computer model output) may be known to satisfy inequality constraints with respect to some or all input variables. Our aim is to build a model capable of incorporating both data interpolation and inequality constraints into a Gaussian process emulator. By using a functional decomposition, we propose to approximate the original Gaussian process by a finite-dimensional Gaussian process such that all conditional simulations satisfy the inequality constraints in the whole domain. The mean, mode (maximum a posteriori) and prediction intervals (uncertainty quantification) of the conditional Gaussian process are calculated. To investigate the performance of the proposed model, some conditional simulations with inequality constraints such as boundary, monotonicity or convexity conditions are given. 1. Introduction. In the engineering activity, runs of a computer code can be expensive and time-consuming. One solution is to use a statistical surrogate for conditioning computer model outputs at some input locations (design points). Gaussian process (GP) emulator is one of the most popular choices [23]. The reason comes from the property of the GP that uncertainty quantification can be calculated. Furthermore, it has several nice properties. For example, the conditional GP at observation data (linear equality constraints) is still a GP [5]. Additionally, some inequality constraints (such as monotonicity and convexity) of output computer responses are related to partial derivatives. In such cases, the partial derivatives of the GP are also GPs. Incorporating an infinite number of linear inequality constraints into a GP emulator, the problem becomes more difficult. The reason is that the resulting conditional process is not a GP. In the literature of interpolation with inequality constraints, we find two types of meth-ods. The first one is deterministic and based on splines, which have the advantage that the inequality constraints are satisfied in the whole input domain (see e.g. [16], [24] and [25]). The second one is based on the simulation of the conditional GP by using the subdivision of the input set (see e.g. [1], [6] and [11]). In that case, the inequality constraints are satisfied in a finite number of input locations. Notice that the advantage of such method is that un-certainty quantification can be calculated. In previous work, some methodologies have been based on the knowledge of the derivatives of the GP at some input locations ([11], [21] and [26]). For monotonicity constraints with noisy data, a Bayesian approach was developed in [21]. In [11] the problem is to build a GP emulator by using the prior monotonicity


title: The Vaserstein symbol on real smooth affine threefolds

abstract: We give a necessary and sufficient topological condition for the Vaserstein symbol to be injective on smooth affine real threefolds. More precisely, we show that the Vaserstein symbol is a bijection for such a threefold X if and only if the set of compact connected components of the real points of X (endowed with the Euclidean topology) is empty.


title: Holographic Space-time, Newton's Law and the Dynamics of Black Holes

abstract: We revisit the construction of models of quantum gravity in d dimensional Minkowski space in terms of random tensor models, and correct some mistakes in our previous treatment of the subject. We find a large class of models in which the large impact parameter scattering scales with energy and impact parameter like Newton`s law. These same models also have emergent energy, momentum and angular conservation laws, despite being based on time dependent Hamiltonians. Many of the scattering amplitudes have a Feynman diagram like structure: local interaction vertices connected by propagation of free particles (really Sterman-Weinberg jets of particles). However, there are also amplitudes where jets collide to form large meta-stable objects, with all the scaling properties of black holes: energy, entropy and temperature, as well as the characteristic time scale for the decay of perturbations. We generalize the conjecture of Sekino and Susskind, to claim that all of these models are fast scramblers. The rationale for this claim is that the interactions are invariant under fuzzy subgroups of the group of volume preserving diffeomorphisms, so that they are highly non-local on the holographic screen. We review how this formalism resolves the Firewall Paradox.


title: Obscured AGNs in Bulgeless Hosts discovered by WISE: The Case Study of SDSS J1224+5555

abstract: There is mounting evidence that supermassive black holes form and grow in bulgeless galaxies. However, a robust determination of the fraction of AGNs in bulgeless galaxies, an important constraint to models of supermassive black hole seed formation and merger-free models of AGN fueling, is unknown, since optical studies have been shown to be incomplete for low mass AGNs. In a recent study using the Wide-field Infrared Survey Explorer, we discovered hundreds of bulgeless galaxies that display mid-infrared signatures of extremely hot dust suggestive of powerful accreting massive black holes, despite having no signatures of black hole activity at optical wavelengths. Here we report X-ray follow-up observations of J122434.66+555522.3, a nearby (z=0.052) isolated bulgeless galaxy that contains an unresolved X-ray source detected at the 3 sigma level by XMM-Newton with an observed luminosity uncorrected for intrinsic absorption of L{2-10~keV}=1.1+/-0.4 10^40 ergs/s. Ground-based near-infrared spectroscopy with the Large Binocular Telescope together with multiwavelength observations from ultraviolet to millimeter wavelengths together suggest that J1224+5555 harbors a highly absorbed AGN with an intrinsic absorption of ~NH >10^24 cm^-2. The hard X-ray luminosity of the putative AGN corrected for absorption is L{2-10~keV}~3x10^42 ergs/s, which, depending on the bolometric correction factor, corresponds to a bolometric luminosity of the AGN of 6x10^43ergs/s - 3x10^44 erg/s, and a lower mass limit for the black hole of MBH~2x10^6 Msun, based on the Eddington limit.


title: End-to-end LSTM-based dialog control optimized with supervised and reinforcement learning

abstract: This paper presents a model for end-to-end learning of task-oriented dialog systems. The main component of the model is a recurrent neural network (an LSTM), which maps from raw dialog history directly to a distribution over system actions. The LSTM automatically infers a representation of dialog history, which relieves the system developer of much of the manual feature engineering of dialog state. In addition, the developer can provide software that expresses business rules and provides access to programmatic APIs, enabling the LSTM to take actions in the real world on behalf of the user. The LSTM can be optimized using supervised learning (SL), where a domain expert provides example dialogs which the LSTM should imitate; or using reinforcement learning (RL), where the system improves by interacting directly with end users. Experiments show that SL and RL are complementary: SL alone can derive a reasonable initial policy from a small number of training dialogs; and starting RL optimization with a policy trained with SL substantially accelerates the learning rate of RL.


title: A data driven network approach to rank countries production diversity and food specialization

abstract: The easy access to large data sets has allowed for leveraging methodology in network physics and complexity science to disentangle patterns and processes directly from the data, leading to key insights in the behavior of systems. Here we use to country specific food production data to study binary and weighted topological properties of the bipartite country-food production matrix. This country-food production matrix can be: 1) transformed into overlap matrices which embed information regarding shared production of products among countries, and or shared countries for individual products, 2) identify subsets of countries which produce similar commodities or subsets of commodities shared by a given country allowing for visualization of correlations in large networks, and 3) used to rank country's fitness (the ability to produce a diverse array of products weighted on the type of food commodities) and food specialization (quantified on the number of countries producing that food product weighted on their fitness). Our results show that, on average, countries with high fitness producing highly specialized food commodities also produce low specialization goods, while nations with low fitness producing a small basket of diverse food products, typically produce low specialized food commodities.


title: Fermi-LAT kills dark matter interpretations of AMS-02 data. Or not?

abstract: A number of papers attempt to explain the positron anomaly in cosmic rays, observed by PAMELA and AMS-02, in terms of dark matter (DM) decays or annihilations. However, the recent progress in cosmic gamma-ray studies challenges these attempts. Indeed, as we show, any rational DM model explaining the positron anomaly abundantly produces final state radiation and Inverse Compton gamma rays, which inevitably leads to a contradiction with Fermi-LAT isotropic diffuse gamma-ray background measurements. Furthermore, the Fermi-LAT observation of Milky Way dwarf satellites, supposed to be rich in DM, revealed no significant signal in gamma rays. We propose a generic approach in which the major contribution to cosmic rays comes from the dark matter disc and prove that the tension between the DM origin of the positron anomaly and the cosmic gamma-ray observations can be relieved. We consider both a simple model, in which DM decay/annihilate into charged leptons, and a model-independent minimal case of particle production, and we estimate the optimal thickness of DM disk. Possible mechanisms of formation and its properties are briefly discussed.


title: Bypassing damaged nervous tissue

abstract: It is shown the principal possibility of bypassing damaged demyelinated portions of the nervous tissue, thereby restoring its normal function for the passage of action potentials.


title: Nuclear energy gradients for internally contracted complete active space second-order perturbation theory: Multistate extensions

abstract: We report the development of the theory and computer program for analytical nuclear energy gradients for (extended) multi-state complete active space perturbation theory (CASPT2) with full internal contraction. The vertical shifts are also considered in this work. This is an extension of the fully internally contracted CASPT2 nuclear gradient program, recently developed for a state-specific variant by us [MacLeod and Shiozaki, J. Chem. Phys. 142, 051103 (2015)]; in this extension, the so-called {\lambda} equation is solved to account for the variation of the multi-state CASPT2 energies with respect to the change in the amplitudes obtained in the preceding state-specific CASPT2 calculations, and the Z-vector equations are modified accordingly. The program is parallelized using the MPI3 remote memory access protocol that allows us to perform efficient one-sided communication. The optimized geometries of the ground and excited states of a copper corrole and benzophenone are presented as numerical examples. The code is publicly available under the GNU General Public License.


title: Prefix frequency of lost positions

abstract: The concept of "lost positions" is a recently introduced tool for counting the number of runs in words. We investigate the frequency of lost positions in prefixes of words. This leads to an algorithm that allows to show, using an extensive computer search, that the asymptotic density of runs in binary words is less than $183/193\approx 0.9482$.


title: Predicting with Distributions

abstract: We consider a new learning model in which a joint distribution over vector pairs $(x,y)$ is determined by an unknown function $c(x)$ that maps input vectors $x$ not to individual outputs, but to entire {\em distributions\/} over output vectors $y$. Our main results take the form of rather general reductions from our model to algorithms for PAC learning the function class and the distribution class separately, and show that virtually every such combination yields an efficient algorithm in our model. Our methods include a randomized reduction to classification noise and an application of Le Cam's method to obtain robust learning algorithms.


title: Study of \chi_{bJ}(1P) Properties in the Radiative \Upsilon(2S) Decays

abstract: We report a study of radiative decays of \chi_{bJ}(1P)(J=0,1,2) mesons into 74 hadronic final states comprising charged and neutral pions, kaons, protons; out of these, 41 modes are observed with at least 5 standard deviation significance. Our measurements not only improve the previous measurements by the CLEO Collaboration but also lead to first observations in many new modes. The large sample allows us to probe the total decay width of the \chi_{b0}(1P). In the absence of a statistically significant result, a 90% confidence-level upper limit is set on the width at \Gamma_{total}< 2.4 MeV. Our results are based on 24.7 fb^{-1} of e+e- collision data recorded by the Belle detector at the \Upsilon(2S) resonance, corresponding to (157.8\pm3.6)\times10^6 \Upsilon(2S) decays.


title: Lattice study of large $N_c$ QCD

abstract: We present a lattice simulation study of large $N_c$ regularities of meson and baryon spectroscopy in $SU(N_c)$ gauge theory with two flavors of dynamical fundamental representation fermions. Systems investigated include $N_c=2$, 3, 4, and 5, over a range of fermion masses parametrized by a squared pseudoscalar to vector meson mass ratio between about 0.2 to 0.7. Good agreement with large $N_c$ scaling is observed in the static potential, in meson masses and decay constants, and in baryon spectroscopy. This is an update of the published version of the paper (Phys. Rev. D94 (2016) 034506).


title: Simple D-module components of local cohomology modules

abstract: For a projective variety V in P^n over a field of characteristic zero, with homogeneous ideal I in A = k[x0,x1,...,xn], we consider the local cohomology modules H^i_I(A). These have a structure of holonomic D-module over A, and we investigate their filtration by simple D-modules. In case V is nonsingular, we can describe completely these simple components in terms of the Betti numbers of V.


title: Generation of energy-entangled W states via parametric fluorescence in integrated devices

abstract: Tripartite entangled states, such as GHZ and W states, are typically generated by manipulating two pairs of polarization-entangled photons in bulk optics. Here we propose a scheme to generate W states that are entangled in the energy degree of freedom in an integrated optical circuit. Our approach employs photon pairs generated by spontaneous four-wave mixing (SFWM) in a micro- ring resonator. We also present a feasible procedure for demonstrating the generation of such a state, and we compare polarization-entangled and energy-entangled schemes for the preparation of W states.


title: Dependency Parsing as Head Selection

abstract: Conventional graph-based dependency parsers guarantee a tree structure both during training and inference. Instead, we formalize dependency parsing as the problem of independently selecting the head of each word in a sentence. Our model which we call \textsc{DeNSe} (as shorthand for {\bf De}pendency {\bf N}eural {\bf Se}lection) produces a distribution over possible heads for each word using features obtained from a bidirectional recurrent neural network. Without enforcing structural constraints during training, \textsc{DeNSe} generates (at inference time) trees for the overwhelming majority of sentences, while non-tree outputs can be adjusted with a maximum spanning tree algorithm. We evaluate \textsc{DeNSe} on four languages (English, Chinese, Czech, and German) with varying degrees of non-projectivity. Despite the simplicity of the approach, our parsers are on par with the state of the art.


title: Bicyclic graphs with extremal degree resistance distance

abstract: Let $r(u,v)$ be the resistance distance between two vertices $u, v$ of a simple graph $G$, which is the effective resistance between the vertices in the corresponding electrical network constructed from $G$ by replacing each edge of $G$ with a unit resistor. The degree resistance distance of a simple graph $G$ is defined as ${D_R}(G) = \sum\limits_{\{u,v\} \subseteq V(G)} {[d(u) + d(v)]r(u,v)},$ where $d(u)$ is the degree of the vertex $u$. In this paper, the bicyclic graphs with extremal degree resistance distance are strong-minded. We first determine the $n$-vertex bicyclic graphs having precisely two cycles with minimum and maximum degree resistance distance. We then completely characterize the bicyclic graphs with extremal degree resistance distance.


title: Canonical formulation and path integral for local vacuum energy sequestering

abstract: We establish the Hamiltonian analysis and the canonical path integral for a local formulation of vacuum energy sequestering. In particular, by considering the state of the universe as a superposition of vacuum states corresponding to different values of the cosmological and gravitational constants, the path integral is extended to include integrations over the cosmological and gravitational constants. The result is an extension of the Ng-van Dam form of the path integral of unimodular gravity. It is argued to imply a relation between the fraction of the most likely values of the gravitational and cosmological constants and the average values of the energy density and pressure of matter over spacetime. Finally, we construct and analyze a BRST-exact formulation of the theory, which can be considered as a topological field theory.


title: Enhancing the LexVec Distributed Word Representation Model Using Positional Contexts and External Memory

abstract: In this paper we take a state-of-the-art model for distributed word representation that explicitly factorizes the positive pointwise mutual information (PPMI) matrix using window sampling and negative sampling and address two of its shortcomings. We improve syntactic performance by using positional contexts, and solve the need to store the PPMI matrix in memory by working on aggregate data in external memory. The effectiveness of both modifications is shown using word similarity and analogy tasks.


title: Statistical Pattern Recognition for Driving Styles Based on Bayesian Probability and Kernel Density Estimation

abstract: Driving styles have a great influence on vehicle fuel economy, active safety, and drivability. To recognize driving styles of path-tracking behaviors for different divers, a statistical pattern-recognition method is developed to deal with the uncertainty of driving styles or characteristics based on probability density estimation. First, to describe driver path-tracking styles, vehicle speed and throttle opening are selected as the discriminative parameters, and a conditional kernel density function of vehicle speed and throttle opening is built, respectively, to describe the uncertainty and probability of two representative driving styles, e.g., aggressive and normal. Meanwhile, a posterior probability of each element in feature vector is obtained using full Bayesian theory. Second, a Euclidean distance method is involved to decide to which class the driver should be subject instead of calculating the complex covariance between every two elements of feature vectors. By comparing the Euclidean distance between every elements in feature vector, driving styles are classified into seven levels ranging from low normal to high aggressive. Subsequently, to show benefits of the proposed pattern-recognition method, a cross-validated method is used, compared with a fuzzy logic-based pattern-recognition method. The experiment results show that the proposed statistical pattern-recognition method for driving styles based on kernel density estimation is more efficient and stable than the fuzzy logic-based method.


title: Spread of a Catalytic Branching Random Walk on a Multidimensional Lattice

abstract: For a supercritical catalytic branching random walk on Z^d (d is positive integer) with an arbitrary finite catalysts set we study the spread of particles population as time grows to infinity. Namely, we divide by t the position coordinates of each particle existing at time t and then let t tend to infinity. It is shown that in the limit there are a.s. no particles outside the closed convex surface in R^d which we call the propagation front and, under condition of infinite number of visits of the catalysts set, a.s. there exist particles on the propagation front. We also demonstrate that the propagation front is asymptotically densely populated and derive its alternative representation. Recent strong limit theorems for total and local particles numbers established by the author play an essential role. The results obtained develop ones by Ph.Carmona and Y.Hu (2014) devoted to the spread of catalytic branching random walk on Z. Keywords and phrases: branching random walk, supercritical regime, spread of population, propagation front, many-to-one lemma.


title: Incorporating long-range consistency in CNN-based texture generation

abstract: Gatys et al. (2015) showed that pair-wise products of features in a convolutional network are a very effective representation of image textures. We propose a simple modification to that representation which makes it possible to incorporate long-range structure into image generation, and to render images that satisfy various symmetry constraints. We show how this can greatly improve rendering of regular textures and of images that contain other kinds of symmetric structure. We also present applications to inpainting and season transfer.


title: Nonparametric hypersurfaces moving by powers of Gauss curvature

abstract: We study asymptotic behavior of nonparametric hypersurfaces moving by $\alpha$ powers of Gauss curvature $\alpha > 1/n$. Our work generalizes the results of V. Oliker [Oli91] for $\alpha= 1$.


title: A Survey of Research on Control of Teams of Small Robots in Military Operations

abstract: While a number of excellent review articles on military robots have appeared in existing literature, this paper focuses on a distinct sub-space of related problems: small military robots organized into moderately sized squads, operating in a ground combat environment. Specifically, we consider the following: - Command of practical small robots, comparable to current generation, small unmanned ground vehicles (e.g., PackBots) with limited computing and sensor payload, as opposed to larger vehicle-sized robots or micro-scale robots; - Utilization of moderately sized practical forces of 3-10 robots applicable to currently envisioned military ground operations; - Complex three-dimensional physical environments, such as urban areas or mountainous terrains and the inherent difficulties they impose, including limited and variable fields of observation, difficult navigation, and intermittent communication; - Adversarial environments where the active, intelligent enemy is the key consideration in determining the behavior of the robotic force; and - Purposeful, partly autonomous, coordinated behaviors that are necessary for such a robotic force to survive and complete missions; these are far more complex than, for example, formation control or field coverage behavior.


title: Conforming restricted Delaunay mesh generation for piecewise smooth complexes

abstract: A Frontal-Delaunay refinement algorithm for mesh generation in piecewise smooth domains is described. Built using a restricted Delaunay framework, this new algorithm combines a number of novel features, including: (i) an unweighted, conforming restricted Delaunay representation for domains specified as a (non-manifold) collection of piecewise smooth surface patches and curve segments, (ii) a protection strategy for domains containing curve segments that subtend sharply acute angles, and (iii) a new class of off-centre refinement rules designed to achieve high-quality point-placement along embedded curve features. Experimental comparisons show that the new Frontal-Delaunay algorithm outperforms a classical (statically weighted) restricted Delaunay-refinement technique for a number of three-dimensional benchmark problems.


title: Fluctuations of a surface relaxation model in interacting scale free networks

abstract: Isolated complex networks have been studied deeply in the last decades due to the fact that many real systems can be modeled using these types of structures. However, it is well known that the behavior of a system not only depends on itself, but usually also depends on the dynamics of other structures. For this reason, interacting complex networks and the processes developed on them have been the focus of study of many researches in the last years. One of the most studied subjects in this type of structures is the Synchronization problem, which is important in a wide variety of processes in real systems. In this manuscript we study the synchronization of two interacting scale-free networks, in which each node has $ke$ dependency links with different nodes in the other network. We map the synchronization problem with an interface growth, by studying the fluctuations in the steady state of a scalar field defined in both networks. We find that as $ke$ slightly increases from $ke=0$, there is a really significant decreasing in the fluctuations of the system. However, this considerable improvement takes place mainly for small values of $ke$, when the interaction between networks becomes stronger there is only a slight change in the fluctuations. We characterize how the dispersion of the scalar field depends on the internal degree, and we show that a combination between the decreasing of this dispersion and the integer nature of our growth model are the responsible for the behavior of the fluctuations with $ke$.


title: Minimal prospects for radio detection of extensive air showers in the atmosphere of Jupiter

abstract: One possible approach for detecting ultra-high-energy cosmic rays and neutrinos is to search for radio emission from extensive air showers created when they interact in the atmosphere of Jupiter, effectively utilizing Jupiter as a particle detector. We investigate the potential of this approach. For searches with current or planned radio telescopes we find that the effective area for detection of cosmic rays is substantial (~3*10^7 km^2), but the acceptance angle is so small that the typical geometric aperture (~10^3 km^2 sr) is less than that of existing terrestrial detectors, and cosmic rays also cannot be detected below an extremely high threshold energy (~10^23 eV). The geometric aperture for neutrinos is slightly larger, and greater sensitivity can be achieved with a radio detector on a Jupiter-orbiting satellite, but in neither case is this sufficient to constitute a practical detection technique. Exploitation of the large surface area of Jupiter for detecting ultra-high-energy particles remains a long-term prospect that will require a different technique, such as orbital fluorescence detection.


title: An Attentional Neural Conversation Model with Improved Specificity

abstract: In this paper we propose a neural conversation model for conducting dialogues. We demonstrate the use of this model to generate help desk responses, where users are asking questions about PC applications. Our model is distinguished by two characteristics. First, it models intention across turns with a recurrent network, and incorporates an attention model that is conditioned on the representation of intention. Secondly, it avoids generating non-specific responses by incorporating an IDF term in the objective function. The model is evaluated both as a pure generation model in which a help-desk response is generated from scratch, and as a retrieval model with performance measured using recall rates of the correct response. Experimental results indicate that the model outperforms previously proposed neural conversation architectures, and that using specificity in the objective function significantly improves performances for both generation and retrieval.


title: Robust Inversion Methods for Aerosol Spectroscopy

abstract: The Fast Aerosol Spectrometer (FASP) is a device for spectral aerosol measurements. Its purpose is to safely monitor the atmosphere inside a reactor containment. First we describe the FASP and explain its basic physical laws. Then we introduce our reconstruction methods for aerosol particle size distributions designed for the FASP. We extend known existence results for constrained Tikhonov regularization by uniqueness criteria and use those to generate reasonable models for the size distributions. We apply a Bayesian model-selection framework on these pre-generated models. We compare our algorithm with classical inversion methods using simulated measurements. We then extend our reconstruction algorithm for two-component aerosols, so that we can simultaneously retrieve their particle-size distributions and unknown volume fractions of their two components. Finally we present the results of a numerical study for the extended algorithm.


title: Congruence primes for automorphic forms on unitary groups and applications to the arithmetic of Ikeda lifts

abstract: In this paper we provide a sufficient condition for a prime to be a congruence prime for an automorphic form $f$ on the unitary group $U(n,n)(A_F)$ for a large class of totally real fields $F$ via a divisibility of a special value of the standard $L$-function associated to $f$. We also study $\ell$-adic properties of the Fourier coefficients of an Ikeda lift $I_{\phi}$ (of an elliptic modular form $\phi$) on $U(n,n)(A_{\mathbf{Q}})$ proving that they are $\ell$-adic integers which do not all vanish modulo $\ell$. Finally we combine these results to show that the condition of $\ell$ being a congruence prime for $I_{\phi}$ is controlled by the $\ell$-divisibility of a product of special values of the symmetric square $L$-function of $\phi$. We close the paper by computing an example when our main theorem applies.


title: Weighted $\ell_1$-Minimization for Sparse Recovery under Arbitrary Prior Information

abstract: Weighted $\ell_1$-minimization has been studied as a technique for the reconstruction of a sparse signal from compressively sampled measurements when prior information about the signal, in the form of a support estimate, is available. In this work, we study the recovery conditions and the associated recovery guarantees of weighted $\ell_1$-minimization when arbitrarily many distinct weights are permitted. For example, such a setup might be used when one has multiple estimates for the support of a signal, and these estimates have varying degrees of accuracy. Our analysis yields an extension to existing works that assume only a single support estimate set upon which a constant weight is applied. We include numerical experiments, with both synthetic signals and real video data, that demonstrate the benefits of allowing non-uniform weights in the reconstruction procedure.


title: Discriminants and Automorphism Groups of Veronese subrings of skew polynomial rings

abstract: We study important invariants and properties of the Veronese subalgebras of $q$-skew polynomial rings, including their discriminant, center and automorphism group, as well as cancellation property and the Tits alternative.


title: Coarse indices of twisted operators

abstract: Several formulas for computing coarse indices of twisted Dirac type operators are introduced. One type of such formulas is by composition product in $E$-theory. The other type is by module multiplications in $K$-theory, which also yields an index theoretic interpretation of the duality between Roe algebra and stable Higson corona.


title: Two Higgs doublet models augmented by a scalar color octet

abstract: The LHC is now studying in detail the couplings of the Higgs boson in order to determine if there is new physics. Many recent studies have examined the available fits to Higgs couplings from the perspective of constraining two Higgs doublet models (2HDM). In this paper we extend those studies to include constraints on the one loop couplings of the Higgs to gluons and photons. These couplings are particularly sensitive to the existence of new colored particles that are hard to detect otherwise and we use them to constrain a 2HDM augmented with a color-octet scalar, a possibility motivated by minimal flavor violation. We first study theoretical constraints on this model and then compare them with LHC measurements.


title: RAISR: Rapid and Accurate Image Super Resolution

abstract: Given an image, we wish to produce an image of larger size with significantly more pixels and higher image quality. This is generally known as the Single Image Super-Resolution (SISR) problem. The idea is that with sufficient training data (corresponding pairs of low and high resolution images) we can learn set of filters (i.e. a mapping) that when applied to given image that is not in the training set, will produce a higher resolution version of it, where the learning is preferably low complexity. In our proposed approach, the run-time is more than one to two orders of magnitude faster than the best competing methods currently available, while producing results comparable or better than state-of-the-art. A closely related topic is image sharpening and contrast enhancement, i.e., improving the visual quality of a blurry image by amplifying the underlying details (a wide range of frequencies). Our approach additionally includes an extremely efficient way to produce an image that is significantly sharper than the input blurry one, without introducing artifacts such as halos and noise amplification. We illustrate how this effective sharpening algorithm, in addition to being of independent interest, can be used as a pre-processing step to induce the learning of more effective upscaling filters with built-in sharpening and contrast enhancement effect.


title: Scale-free channeling patterns near the onset of erosion of sheared granular beds

abstract: Erosion shapes our landscape and occurs when a sufficient shear stress is exerted by a fluid on a sedimented layer. What controls erosion at a microscopic level remains debated, especially near the threshold forcing where it stops. Here we study experimentally the collective dynamics of the moving particles, using a set-up where the system spontaneously evolves toward the erosion onset. We find that the spatial organization of the erosion flux is heterogeneous in space, and occurs along channels of local flux $\sigma$ whose distribution displays scaling near threshold and follows $P(\sigma)\sim J/\sigma$, where $J$ is the mean erosion flux. Channels are strongly correlated in the direction of forcing but not in the transverse direction. We show that these results quantitatively agree with a model where the dynamics is governed by the competition of disorder (which channels mobile particles) and particle interactions (which reduces channeling). These observations support that for laminar flows, erosion is a dynamical phase transition which shares similarity with the plastic depinning transition occurring in dirty superconductors. The methodology we introduce here could be applied to probe these systems as well.


title: Vector magnetometry using silicon vacancies in 4H-SiC at ambient conditions

abstract: Point defects in solids promise precise measurements of various quantities. Especially magnetic field sensing using the spin of point defects has been of great interest recently. When optical readout of spin states is used, point defects achieve optical magnetic imaging with high spatial resolution at ambient conditions. Here, we demonstrate that genuine optical vector magnetometry can be realized using the silicon vacancy in SiC, which has an uncommon S=3/2 spin. To this end, we develop and experimentally test sensing protocols based on a reference field approach combined with multi frequency spin excitation. Our works suggest that the silicon vacancy in an industry-friendly platform, SiC, has potential for various magnetometry applications at ambient conditions.


title: RENORM tensor-Pomeron diffractive predictions

abstract: Predictions of the elastic scattering, total-inetastic, and total proton-proton cross sections, based on a Regge theory inspired tensor-Pomeron implementation of the RENORM model for hadronic diffraction, are compared to the latest experimental measurements at the LHC. The measured cross sections are in good agreement within the experimental uncertainties of the data and the theoretical uncertainties of the model, reaching down to the ~1% level.


title: Experimental demonstration of acoustic wave induced magnetization switching of dipole coupled magnetostrictive nanomagnets for ultralow power computing

abstract: We report nanomagnetic switching with Acoustic Waves (AW) launched from interdigitated electrodes that modulate the stress anisotropy of elliptical cobalt nanoscale magnetostrictive magnets (340 nm x 270 nm x 12 nm) delineated on 128 degree Y-cut lithium niobate. The dipole-coupled nanomagnet pairs are in a single-domain state and are initially magnetized along the major axis of the ellipse, with their magnetizations parallel to each other. The magnetizations of nanomagnets having lower shape anisotropy are reversed upon acoustic wave propagation. Thereafter, the magnetization of these nanomagnets remains in the reversed state and demonstrate non-volatility. This executes a 'NOT' operation. This proof of acoustic wave induced magnetic state reversal in dipole-coupled nanomagnets implementing a 'NOT' gate operation could potentially lead to the development of extremely energy-efficient nanomagnetic logic. Furthermore, fabrication complexity is reduced immensely due to the absence of individual contacts to the nanomagnets, leading to lower energy dissipation


title: Lyman-alpha and CIII] Emission in z=7-9 Galaxies: Accelerated Reionization Around Luminous Star Forming Systems?

abstract: We discuss new Keck/MOSFIRE spectroscopic observations of four luminous galaxies at z~7-9 selected to have intense optical line emission by Roberts-Borsani et al. (2016). Previous follow-up has revealed Lyman-alpha in two of the four galaxies. Our new MOSFIRE observations confirm that Lyman-alpha is present in the entire sample. We detect Lyman-alpha emission in COS-zs7-1, confirming its redshift as z=7.154, and we detect Lyman-alpha in EGS-zs8-2 at z=7.477, verifying a tentative detection presented in an earlier study. The ubiquity of Lyman-alpha in this sample is puzzling given that the IGM is likely significantly neutral over 7<z<9. To investigate this result in more detail, we have initiated a campaign to target UV metal emission in the four Lyman-alpha emitters as a probe of both the radiation field and the velocity offset of Lyman-alpha. Here we present the detection of intense CIII] emission in EGS-zs8-1, a galaxy from this sample previously shown to have Lyman-alpha at z=7.73. Photoionization models indicate that an intense radiation field and low metallicity are required to reproduce the intense CIII] and optical line emission. We argue that this extreme radiation field is likely to affect the local environment, increasing the transmission of Lyman-alpha through the galaxy. Moreover, the centroid of CIII] indicates that Lyman-alpha is redshifted from the systemic value by 340 km/s. This velocity offset is larger than that seen in less luminous systems, providing an additional explanation for the transmission of Lyman-alpha emission through the IGM. Since the transmission is further enhanced by the likelihood that such systems are also situated in the densest regions with the largest ionized bubbles, the visibility of Lyman-alpha at z>7 is expected to be strongly luminosity-dependent, with the most effective transmission occurring in systems with intense star formation.


title: Zoneout: Regularizing RNNs by Randomly Preserving Hidden Activations

abstract: We propose zoneout, a novel method for regularizing RNNs. At each timestep, zoneout stochastically forces some hidden units to maintain their previous values. Like dropout, zoneout uses random noise to train a pseudo-ensemble, improving generalization. But by preserving instead of dropping hidden units, gradient information and state information are more readily propagated through time, as in feedforward stochastic depth networks. We perform an empirical investigation of various RNN regularizers, and find that zoneout gives significant performance improvements across tasks. We achieve competitive results with relatively simple models in character- and word-level language modelling on the Penn Treebank and Text8 datasets, and combining with recurrent batch normalization yields state-of-the-art results on permuted sequential MNIST.


title: Geometrical symmetries of nuclear systems: D(3h) and T(d) symmetries in light nuclei

abstract: The role of discrete (or point-group) symmetries in alpha-cluster nuclei is discussed in the framework of the algebraic cluster model which describes the relative motion of the alpha-particles. Particular attention is paid to the discrete symmetry of the geometric arrangement of the alpha-particles, and the consequences for the structure of the corresponding rotational bands. The method is applied to study cluster states in the nuclei 12C and 16O. The observed level sequences can be understood in a simple way as a consequence of the underlying discrete symmetry that characterizes the geometrical configuration of the alpha-particles, i.e. an equilateral triangle with D(3h) symmetry for 12C, and a tetrahedron with T(d) symmetry for 16O. The structure of rotational bands provides a fingerprint of the underlying geometrical configuration of alpha-particles.


title: Scene Grammars, Factor Graphs, and Belief Propagation

abstract: We describe a general framework for probabilistic modeling of complex scenes and inference from ambiguous observations. The approach is motivated by applications in image analysis and is based on the use of priors defined by stochastic grammars. We define a class of grammars that capture relationships between the objects in a scene and provide important contextual cues for statistical inference. The distribution over scenes defined by a probabilistic scene grammar can be represented by a graphical model and this construction can be used for efficient inference with loopy belief propagation. We show experimental results with two different applications. One application involves the reconstruction of binary contour maps. Another application involves detecting and localizing faces in images. In both applications the same framework leads to robust inference algorithms that can effectively combine local information to reason about a scene.


title: Object Kinetic Monte Carlo Simulations of Radiation Damage in Bulk Tungsten Part-II: With a PKA Spectrum Corresponding to 14-MeV Neutrons

abstract: Object kinetic Monte Carlo was employed to study the effect of dose rate on the evolution of vacancy microstructure in polycrystalline tungsten under neutron bombardment. The evolution was followed up to 1.0 displacement per atom (dpa) with point defects generated in accordance with a primary knock-on atom (PKA) spectrum corresponding to 14-MeV neutrons. The present study includes the effect of grain size (2.0 and 4.0 $\mu$m) but excludes the impact of transmutation or pre-existing defects beyond grain boundary sinks. Vacancy cluster density increases with dose rate, while the density of vacancies decreases. Consequently, the average vacancy cluster size and the fraction of vacancies in visible clusters decrease with increasing dose rate. The density of vacancies and vacancy clusters decrease with grain size such that the average size of the clusters remains similar. However, the average size is larger for larger grains at dose rates < 4.5 x 10-7 dpa/s. The trend of vacancy accumulation as a function of dose, dose rate, and grain size is similar to that obtained with the High Flux Isotope Reactor (HFIR) PKA spectrum. However, the amount of vacancy accumulation and the vacancy microstructure are quite different. Compared to the HFIR case, we find that even though the dose rates are 2.5 times higher, the density of vacancies and the average vacancy cluster sizes are lower. In addition, a void lattice forms only for the lowest two dose rates (4.5 x 10-8 and 4.5 x 10-9 dpa/s). In contrast, a void lattice formed at all dose rates studied using the HFIR PKA spectrum. We discuss in detail the factors that lead to these different microstructures.


title: Orbital-Exchange and Fractional Quantum Number Excitations in an f-electron Metal Yb$_2$Pt$_2$Pb

abstract: Exotic quantum states and fractionalized magnetic excitations, such as spinons in one-dimensional chains, are generally viewed as belonging to the domain of 3d transition metal systems with spins 1/2. Our neutron scattering experiments on the 4f-electron metal Yb$_2$Pt$_2$Pb overturn this common wisdom. We observe broad magnetic continuum dispersing in only one direction, which indicates that the underlying elementary excitations are spinons carrying fractional spin-1/2. These spinons are the quantum dynamics of the anisotropic, orbital-dominated Yb moments, and thus these effective quantum spins are emergent variables that encode the electronic orbitals. The unique birthmark of their unusual origin is that only longitudinal spin fluctuations are measurable, while the transverse excitations such as spin waves are virtually invisible to magnetic neutron scattering. The proliferation of these orbital-spinons strips the electrons of their orbital identity, and we thus report here a new electron fractionalization phenomenon, charge-orbital separation.


title: The Unusual Galactic Center Radio Source N3

abstract: Here we report on new, multi-wavelength radio observations of the unusual point source "N3" that appears to be located in the vicinity of the Galactic Center. VLA observations between 2 and 50 GHz reveal that N3 is a compact and bright source (56 mJy at 10 GHz) with a non-thermal spectrum superimposed upon the non-thermal radio filaments (NTFs) of the Radio Arc. Our highest frequency observations place a strict upper limit of 65 x 28 mas on the size of N3. We compare our observations to those of Yusef-Zadeh & Morris (1987) and Lang et al. (1997) and conclude that N3 is variable over long time scales. Additionally, we present the detection of a compact molecular cloud located adjacent to N3 in projection. CH3CN, CH3OH, CS, HC3N, HNCO, SiO, SO, and NH3 are detected in the cloud and most transitions have FWHM line widths of ~20 km/s. The rotational temperature determined from the metastable NH3 transitions ranges from 79 K to 183 K depending on the transitions used. We present evidence that this molecular cloud is interacting with N3. After exploring the relationship between the NTFs, molecular cloud, and N3, we conclude that N3 likely lies within the Galactic Center. We are able to rule out the HII region, young supernova, active star, AGN, and micro-quasar hypotheses for N3. While a micro-blazar may provide a viable explanation for N3, additional observations are needed to determine the physical counterpart of this mysterious source.


title: Validity and regularization of classical half-space equations

abstract: Recent result [Wu and Guo, Comm. Math. Phys., 2015] has shown that over the 2D unit disk, the classical half-space equation (CHS) for the neutron transport does not capture the correct boundary layer behaviour as long believed. In this paper we develop a regularization technique for CHS to any arbitrary order and use its first-order regularization to show that in the case of the 2D unit disk, although CHS misrepresents the boundary layer behaviour, it does give the correct boundary condition for the interior macroscopic (Laplace) equation. Therefore CHS is still a valid equation to recover the correct boundary condition for the interior Laplace equation over the 2D unit disk.


title: Project Icarus: Preliminary Thoughts on the Selection of Probes and Instruments for an Icarus-style Interstellar Mission

abstract: In this paper we outline the range of probes and scientific instruments that will be required in order for Icarus to fulfill its scientific mission of exploring a nearby star, its attendant planetary system, and the intervening interstellar medium. Based on this preliminary analysis, we estimate that the minimum total Icarus scientific payload mass (i.e. the mass of probes and instruments which must be decelerated to rest in the target system to enable a meaningful programme of scientific investigation) will be in the region of 100 tonnes. Of this, approximately 10 tonnes would be allocated for cruise-phase science instruments, and about 35 tonnes (i.e. the average of estimated lower and upper limits of 28 and 41 tonnes) would be contributed by the intra-system science payload itself (i.e. the dry mass of the stellar and planetary probes and their instruments). The remaining ~55 tonnes is allocated for the sub-probe intra-system propulsion requirements (crudely estimated from current Solar System missions; detailed modelling of sub-probe propulsion systems will be needed to refine this figure). The overall mass contributed by the science payload to the total that must be decelerated from the interstellar cruise velocity will be considerably more than 100 tonnes, however, as allowance must be made for the payload structural and infrastructural elements required to support, deploy, and communicate with the science probes and instruments. Based on the earlier Daedalus study, we estimate another factor of two to allow for these components. Pending the outcome of more detailed studies, it therefore appears that an overall science-related payload mass of ~200 tonnes will be required. This paper is a submission of the Project Icarus Study Group.


title: Design of Robust Adaptive Beamforming Algorithms Based on Low-Rank and Cross-Correlation Techniques

abstract: This work presents cost-effective low-rank techniques for designing robust adaptive beamforming (RAB) algorithms. The proposed algorithms are based on the exploitation of the cross-correlation between the array observation data and the output of the beamformer. Firstly, we construct a general linear equation considered in large dimensions whose solution yields the steering vector mismatch. Then, we employ the idea of the full orthogonalization method (FOM), an orthogonal Krylov subspace based method, to iteratively estimate the steering vector mismatch in a reduced-dimensional subspace, resulting in the proposed orthogonal Krylov subspace projection mismatch estimation (OKSPME) method. We also devise adaptive algorithms based on stochastic gradient (SG) and conjugate gradient (CG) techniques to update the beamforming weights with low complexity and avoid any costly matrix inversion. The main advantages of the proposed low-rank and mismatch estimation techniques are their cost-effectiveness when dealing with high dimension subspaces or large sensor arrays. Simulations results show excellent performance in terms of the output signal-to-interference-plus-noise ratio (SINR) of the beamformer among all the compared RAB methods.


title: Optimal Storage Allocation for Wireless Cloud Caching Systems with a Limited Sum Storage Capacity

abstract: In wireless cloud storage systems, the recovery failure probability depends on not only wireless channel conditions but also storage size of each distributed storage node. For an efficient utilization of limited storage capacity and the performance characterization of allocation strategies, we asymptotically analyze the recovery failure probability of a wireless cloud storage system with a sum storage capacity constraint for both high SNR regime and low SNR regime. Then, we find the optimal storage allocation strategy across distributed storage nodes in terms of the asymptotic recovery failure probability. Our analysis reveals that the maximal symmetric allocation is optimal for high SNR regime and the minimal allocation (with $\lfloor T\rfloor$ complete storage nodes and an incomplete storage node) is optimal for low SNR regime, where $T$ is the sum storage capacity. Based on the numerical investigation, we also show that in intermediate SNR regime, a balance allocation between the minimal allocation and the maximal symmetric allocation would not be required if we select one between them according to SNR.


title: Theoretical Modeling Expressions for Networked Enzymatic Signal Processing Steps as Logic Gates Optimized by Filtering

abstract: We describe modeling approaches to a "network" of connected enzyme-catalyzed reactions, with added (bio)chemical processes that introduce biochemical filtering steps into the functioning of such a biocatalytic cascade. Theoretical expressions are derived that allow simple, few-parameter modeling of processes concatenated in such cascades, both with and without filtering. The modeling approach captures and explains features identified in earlier studies of enzymatic processes considered as potential "network components" for multi-step information/signal processing systems.


title: Provable Burer-Monteiro factorization for a class of norm-constrained matrix problems

abstract: We study the projected gradient descent method on low-rank matrix problems with a strongly convex objective. We use the Burer-Monteiro factorization approach to implicitly enforce low-rankness; such factorization introduces non-convexity in the objective. We focus on constraint sets that include both positive semi-definite (PSD) constraints and specific matrix norm-constraints. Such criteria appear in quantum state tomography and phase retrieval applications. We show that non-convex projected gradient descent favors local linear convergence in the factored space. We build our theory on a novel descent lemma, that non-trivially extends recent results on the unconstrained problem. The resulting algorithm is Projected Factored Gradient Descent, abbreviated as ProjFGD, and shows superior performance compared to state of the art on quantum state tomography and sparse phase retrieval applications.

categories: stat.ML cs.DS cs.IT cs.NA math.IT math.OC


title: Allowed patterns of symmetric tent maps via commuter functions

abstract: We introduce a new technique to study pattern avoidance in dynamical systems, namely the use of a commuter function between non-conjugate dynamical systems. We investigate the properties of such a commuter function, specifically $h : [0,1] \to [0,1]$ satisfying $T_1 \circ h = h \circ T_\mu$, where $T_\mu$ denotes a symmetric tent map of height $\mu$. We make use of this commuter function to prove strict inclusion of the set of allowed patterns of $T_\mu$ in the set of allowed patterns of $T_1$.


title: The Origin of Weak Lensing Convergence Peaks

abstract: Weak lensing convergence peaks are a promising tool to probe nonlinear structure evolution at late times, providing additional cosmological information beyond second-order statistics. Previous theoretical and observational studies have shown that the cosmological constraints on $\Omega_m$ and $\sigma_8$ are improved by a factor of up to ~ 2 when peak counts and second-order statistics are combined, compared to using the latter alone. We study the origin of lensing peaks using observational data from the 154 deg$^2$ Canada-France-Hawaii Telescope Lensing Survey. We found that while high peaks (with height $\kappa$ >3.5 $\sigma_\kappa$, where $\sigma_\kappa$ is the r.m.s. of the convergence $\kappa$) are typically due to one single massive halo of ~$10^{15}M_\odot$, low peaks ($\kappa$ <~ $\sigma_\kappa$) are associated with constellations of 2-8 smaller halos (<~$10^{13}M_\odot$). In addition, halos responsible for forming low peaks are found to be significantly offset from the line-of-sight towards the peak center (impact parameter >~ their virial radii), compared with ~0.25 virial radii for halos linked with high peaks, hinting that low peaks are more immune to baryonic processes whose impact is confined to the inner regions of the dark matter halos. Our findings are in good agreement with results from the simulation work by Yang el al. (2011).


title: Modeling and Modifying Response of Biochemical Processes for Biocomputing and Biosensing Signal Processing

abstract: Processes involving multi-input multi-step reaction cascades are used in developing novel biosensing, biocomputing, and decision making systems. In various applications different changes in responses of the constituent processing steps (reactions) in a cascade are desirable in order to allow control of the system's response. Here we consider conversion of convex response to sigmoid by "intensity filtering," as well as "threshold filtering," and we offer a general overview of this field of research. Specifically, we survey rate equation modelling that has been used for enzymatic reactions. This allows us to design modified biochemical processes as "network components" with responses desirable in applications.


title: GMC Collisions as Triggers of Star Formation. II. 3D Turbulent, Magnetized Simulations

abstract: We investigate giant molecular cloud (GMCs) collisions and their ability to induce gravitational instability and thus star formation. This mechanism may be a major driver of star formation activity in galactic disks. We carry out a series of three dimensional, magnetohydrodynamics (MHD), adaptive mesh refinement (AMR) simulations to study how cloud collisions trigger formation of dense filaments and clumps. Heating and cooling functions are implemented based on photo-dissociation region (PDR) models that span the atomic to molecular transition and can return detailed diagnostic information. The clouds are initialized with supersonic turbulence and a range of magnetic field strengths and orientations. Collisions at various velocities and impact parameters are investigated. Comparing and contrasting colliding and non-colliding cases, we characterize morphologies of dense gas, magnetic field structure, cloud kinematic signatures, and cloud dynamics. We present key observational diagnostics of cloud collisions, especially: relative orientations between magnetic fields and density structures, like filaments; $^{13}$CO($J$=2-1), $^{13}$CO($J$=3-2), and $^{12}$CO($J$=8-7) integrated intensity maps and spectra; and cloud virial parameters. We compare these results to observed Galactic clouds.


title: Notes on isocrystals

abstract: For varieties over a perfect field of characteristic p, etale cohomology with Q_l-coefficients is a Weil cohomology theory only when l is not equal to p; the corresponding role for l = p is played by Berthelot's rigid cohomology. In that theory, the coefficient objects analogous to lisse l-adic sheaves are the overconvergent F-isocrystals. This expository article is a brief user's guide for these objects, including some features shared with l-adic cohomology (purity, weights) and some features exclusive to the p-adic case (Newton polygons, convergence and overconvergence). The relationship between the two cases, via the theory of companions, will be treated in a sequel paper.


title: Ferroelectric Domain Patterning Controlled Schottky Junction State in Monolayer MoS$_2$

abstract: We exploit scanning probe controlled domain patterning in a ferroelectric top layer to induce nonvolatile modulation of the conduction characteristic of monolayer MoS$_2$ between a transistor and a junction state. In the presence of a domain wall, MoS$_2$ exhibits rectified I-V that is well described by the thermionic emission model. The induced Schottky barrier height varies from 0.38 eV to 0.57 eV and is tunabe by a SiO$_2$ global back-gate, while the tuning range of depends sensitively on the conduction band tail trapping states. Our work points to a new route to achieving programmable functionalities in van der Waals materials and sheds light on the critical performance limiting factors in these hybrid systems.


title: Improving Coreference Resolution by Learning Entity-Level Distributed Representations

abstract: A long-standing challenge in coreference resolution has been the incorporation of entity-level information - features defined over clusters of mentions instead of mention pairs. We present a neural network based coreference system that produces high-dimensional vector representations for pairs of coreference clusters. Using these representations, our system learns when combining clusters is desirable. We train the system with a learning-to-search algorithm that teaches it which local decisions (cluster merges) will lead to a high-scoring final coreference partition. The system substantially outperforms the current state-of-the-art on the English and Chinese portions of the CoNLL 2012 Shared Task dataset despite using few hand-engineered features.


title: Bayesian index of superiority and the p-value of the conditional test for Poisson parameters

abstract: We consider the problem of comparing two Poisson parameters from the Bayesian perspective. Kawasaki and Miyaoka (2012b) proposed the Bayesian index $P(\lambda_1 < \lambda_2 | X_1,X_2)$ and expressed it using the hypergeometric series. In this paper, under some conditions, we give four other expressions of the Bayesian index in terms of the cumulative distribution functions of beta, $F$, binomial, and negative binomial distribution. Next, we investigate the relationship between the Bayesian index and the $p$-value of the conditional test with the null hypothesis $H_0: \lambda_1 \geq \lambda_2 $ versus an alternative hypothesis $H_1: \lambda_1<\lambda_2 $. Additionally, we investigate the generalized relationship between $P(\lambda_1/\lambda_2 <c | X_1, X_2)$ and the $p$-value of the conditional test with the null hypothesis $H_0: \lambda_1/\lambda_2 \geq c$ versus the alternative $H_1: \lambda_1/\lambda_2 < c$. We illustrate the utility of the Bayesian index using analyses of real data. Our finding suggests that the Bayesian index can potentially be useful in an epidemiology and in a clinical trial.


title: Effect of strong $\bar{\rm p}$-p nuclear forces on the rate of the low-energy three-body protonium formation reaction: $\bar{p} + H_{\mu}(1s) \rightarrow (\bar{p} p)_{\alpha} + \mu^-$

abstract: The effect of the strong $\bar{\rm p}$-p nuclear interaction in a three-charge-particle system with arbitrary masses is investigated. Specifically, the ($\bar{\rm p},\ \mu^-$,\ p) system is considered, where $\bar{p}$ is an antiproton, $\mu^-$ is a muon and p is a proton. A numerical computation in the framework of a detailed few-body approach is carried out for the following protonium (antiprotonic hydrogen) formation three-body reaction: $\bar{p} + H_{\mu}(1s) \rightarrow (\bar{p} p)_{\alpha} + \mu^-$. Here, $H_{\mu}(1s)$ is a ground state muonic hydrogen, i.e. a bound state of p and $\mu^-$. A bound state of $p$ and its counterpart $\bar{p}$ is a protonium atom in a quantum atomic state $\alpha$, i.e. $Pn = (\bar{p}p)_{\alpha}$. The low-energy cross sections and rates of the $Pn$ formation reaction are computed in the framework of a Faddeev-like equation. The strong $\bar{\rm p}$-p interaction is included in these calculations within a first order approximation. It was found, that even in the framework of this approximation the inclusion of the strong interaction results in a quite significant correction to the rate of the three-body reaction. Therefore, the title three-body antiprotonic process with participation of muons should be useful, especially at low-energy collisions, in studying the $\bar{\rm p}$-p nuclear forces and the annihilation channels in $Pn$.


title: Calculations for electron-impact excitation and ionization of beryllium

abstract: The B-spline R-matrix and the convergent close-coupling methods are used to study electron collisions with neutral beryllium over an energy range from threshold to 100 eV. Coupling to the target continuum significantly affects the results for transitions from the ground state, but to a lesser extent the strong transitions between excited states. Cross sections are presented for selected transitions between low-lying physical bound states of beryllium, as well as for elastic scattering, momentum transfer, and ionization. The present cross sections for transitions from the ground state from the two methods are in excellent agreement with each other, and also with other available results based on nonperturbative convergent pseudo-state and time-dependent close-coupling models. The elastic cross section at low energies is dominated by a prominent shape resonance. The ionization from the $(2s2p)^3P$ and $(2s2p)^1P$ states strongly depends on the respective term. The current predictions represent an extensive set of electron scattering data for neutral beryllium, which should be sufficient for most modeling applications.


title: Envelope Functions: Unifications and Further Properties

abstract: Recently, the forward-backward and Douglas-Rachford envelope functions were proposed in the literature. The stationary points of these envelope functions have a close relationship with the solutions of the possibly nonsmooth optimization problem to be solved. The envelopes were shown to be smooth and convex under some additional assumptions. Therefore, these envelope functions create powerful bridges between nonsmooth and smooth optimization. In this paper, we present a general envelope function that unifies and generalizes these envelope functions. We provide properties of the general envelope function that sharpen corresponding known results for the special cases. We also present an envelope function for the generalized alternating projections method (GAP), named the GAP envelope. It enables for convex feasibility problems with two sets, of which one is affine, to be solved by finding any stationary point of the smooth and under some assumptions convex GAP envelope.


title: Improved approximative multicoloring of hexagonal graphs

abstract: In 1999, McDiarmid and Reed conjectued that the approximation ratio $9/8$ of multichromatic number to weighted clique number asymptotically is the best possible for general weighted hexagonal graphs. We prove that there is a proper multicoloring of $G$ that uses at most $15 \lfloor\frac{\omega(G)}{12}\rfloor+ 18 $ colors improving the best previously known asimptotic ratio from 4/3 to 5/4.


title: Transient nutations of an electron spin 1/2 with dipolar coupling to neighbouring nuclear spins

abstract: A model cluster, with a resonant electron spin surrounded with identical nuclear spins occupying the sites of a cubic lattice, and with dipolar coupling between the spins, is used for a theoretical investigation of the transient nutations induced by a strong oscillating field but disturbed by these local dipolar fields. An effective Hamiltonian valid at the time scale of the nutations is first established. An analytical expression for the Dipolar Coherence time tDC is derived by adapting Van Vleck's perturbative approach to the present situation, with weak rather than intermediate dipolar coupling. It is found that the dipolar coupling does not affect tDC at exact ESR resonance (Delta=0), a result reminiscent of a behaviour found in transient NMR in liquids, in an inhomogeneous static field. The expression for tDC also presents a qualitative similarity with a time introduced in that liquid-state NMR context. Introducing a collection of independent clusters, a bridge is built between tDC and the half-width delta of the usual steady state unsaturated ESR line. Assuming a simple cubic lattice, numerical values for (tDC)Delta=delta are given, for both undiluted and dilute samples. It is shown that the TN do not obey the Bloch equations. Precautions to be taken when trying a qualitative comparison between the present results and existing experimental data, and some thermodynamic aspects, are discussed.


abstract: We obtain the Buchdahl compactness limit for a pure Lovelock static fluid star and verify that the limit following from the uniform density Schwarzschild's interior solution, which is universal irrespective of the gravitational theory (Einstein or Lovelock), is true in general. In terms of surface potential $\Phi(r)$, it means at the surface of the star $r=r_{0}$, $\Phi(r_{0}) < 2N(d-N-1)/(d-1)^2$ where $d$, $N$ respectively indicate spacetime dimensions and Lovelock order. For a given $N$, $\Phi(r_{0})$ is maximum for $d=2N+2$ while it is always $4/9$, Buchdahl's limit, for $d=3N+1$. It is also remarkable that for $N=1$ Einstein gravity, or for pure Lovelock in $d=3N+1$, Buchdahl's limit is equivalent to the criteria that gravitational field energy exterior to the star is less than half its gravitational mass, having no reference to the interior at all.


title: Solar Physics Research in the Russian Subcontinent - Current Status and Future

abstract: Modern research in solar physics in Russia is a multifaceted endeavor, which includes multi-wavelength observations from the ground- and space-based instruments, extensive theoretical and numerical modeling studies, new instrument development, and cross-disciplinary and international research. The research is conducted at the research organizations under the auspices of the Russian Academy of Sciences and to a lesser extent, by the research groups at Universities. Here, we review the history of solar physics research in Russia, and provide an update on recent developments.


title: Mild solutions are weak solutions in a class of (non)linear measure-valued evolution equations on a bounded domain

abstract: We study the connection between mild and weak solutions for a class of measure-valued evolution equations on the bounded domain $[0,1]$. Mass moves, driven by a velocity field that is either a function of the spatial variable only, $v=v(x)$, or depends on the solution $\mu$ itself: $v=v[\mu](x)$. The flow is stopped at the boundaries of $[0,1]$, while mass is gated away by a certain right-hand side. In previous works [J. Differential Equations, 259 (2015), pp. 1068-1097] and [SIAM J. Math. Anal., 48 (2016), pp. 1929-1953], we showed the existence and uniqueness of appropriately defined mild solutions for $v=v(x)$ and $v=v[\mu](x)$, respectively. In the current paper we define weak solutions (by specifying the weak formulation and the space of test functions). The main result is that the aforementioned mild solutions are weak solutions, both when $v=v(x)$ and when $v=v[\mu](x)$.


title: Analysis of the multiple system with CP component phi Draconis

abstract: The star phi Dra comprises a spectroscopic binary and a third star that together form a visual triple system. It is one of the brightest chemically peculiar (CP) stars of the upper main sequence. Despite these facts, no comprehensive study of its multiplicity has been performed yet. In this work, we present a detailed analysis of the triple system based on available measurements. We use radial velocities taken from four sources in the literature in a re-analysis of the inner spectroscopic binary (Aab). An incorrect value of the orbital period of the inner system Aab about 27 days was accepted in literature more than forty years. A new solution of orbit with the 128-day period was determined. Relative position measurements of the outer visual binary system (AB) from Washington Double Star Catalog were compared with known orbital models. Furthermore, it was shown that astrometric motion in system AB is well described by the model of Andrade (2005) with a 308-year orbital period. Parameters of A and B components were utilized to estimate individual brightness for all components and their masses from evolutionary tracks. Although we found several facts which support the gravitational bond between them, unbound solution cannot be fully excluded yet.


title: Charmonium dissociation cross sections and charmonium dissociation rates in hadronic matter

abstract: K*-charmonium dissociation reactions in hadronic matter are studied in the Born approximation, in the quark-interchange mechanism, and with a temperature-dependent quark potential. We obtain the temperature dependence of unpolarized cross sections for K*-charmonium dissociation reactions which produce charmed mesons and charmed strange mesons. We use the cross sections for charmonium dissociation in collisions with pion, rho meson, kaon, vector kaon, and eta meson to calculate dissociation rates of charmonium with the five types of mesons. Because of the temperature dependence of the meson masses, dissociation cross sections, and meson distribution functions, the charmonium dissociation rates generally increase with the increase of temperature and decrease with the increase of charmonium momentum from 2.2 GeV/c. We find that the first derivative of the dissociation rate with respect to the charmonium momentum is zero when the charmonium is at rest. While the eta + psi' and eta + chi_c dissociation reactions can be neglected, the J/psi, psi', and chi_c dissociation are caused by the collisions with pion, rho meson, kaon, vector kaon, and eta meson.


title: Estimate of the squeezing function for a class of bounded domains

abstract: We construct a class of bounded domains, on which the squeezing function is not uniformly bounded from below near a smooth and pseudoconvex boundary point.


title: Space-based Microlens Parallax Observation As a Way to Resolve the Severe Degeneracy between Microlens-parallax and Lens-orbital Effect

abstract: In this paper, we demonstrate the severity of the degeneracy between the microlens-parallax and lens-orbital effects by presenting the analysis of the gravitational binary-lens event OGLE-2015-BLG-0768. Despite the obvious deviation from the model based on the the linear observer motion and the static binary, it is found that the residual can be almost equally well explained by either the parallactic motion of the Earth or the rotation of the binary lens axis, resulting in the severe degeneracy between the two effects. We show that the degeneracy can be readily resolved with the additional data provided by space-based microlens parallax observations. Enabling to distinguish between the two higher-order effects, space-based microlens parallax observations will make it possible not only to accurately determine the physical lens parameters but also to further constrain the orbital parameters of binary lenses.


title: Simple fractal calculus from fractal arithmetic

abstract: Non-Newtonian calculus that starts with elementary non-Diophantine arithmetic operations of a Burgin type is applicable to all fractals whose cardinality is continuum. The resulting definitions of derivatives and integrals are simpler from what one finds in the more traditional literature of the subject, and they often work in the cases where the standard methods fail. As an illustration, we perform a Fourier transform of a real-valued function with Sierpi\'nski-set domain. The resulting formalism is as simple as the usual undergraduate calculus.


title: Maximum norm analysis of implicit-explicit backward difference formulas for nonlinear parabolic equations

abstract: We establish optimal order a priori error estimates for implicit-explicit BDF methods for abstract semilinear parabolic equations with time-dependent operators in a complex Banach space settings, under a sharp condition on the non-self-adjointness of the linear operator. Our approach relies on the discrete maximal parabolic regularity of implicit BDF schemes for autonomous linear parabolic equations, recently established in [20], and on ideas from [7]. We illustrate the applicability of our results to four initial and boundary value problems, namely two for second order, one for fractional order, and one for fourth order, namely the Cahn-Hilliard, parabolic equations.


title: A no-go theorem for the $n$-twistor description of a massive particle

abstract: It is proved that the $n$-twistor expression of a particle's four-momentum vector reduces, by a unitary transformation, to the two-twistor expression for a massive particle or the one-twistor expression for a massless particle. Therefore the {\em genuine} $n$-twistor description of a massive particle in four-dimensional Minkowski space fails for the case $n\ge3$.


title: Finding Multiple New Optimal Locations in a Road Network

abstract: We study the problem of optimal location querying for location based services in road networks, which aims to find locations for new servers or facilities. The existing optimal solutions on this problem consider only the cases with one new server. When two or more new servers are to be set up, the problem with minmax cost criteria, MinMax, becomes NP-hard. In this work we identify some useful properties about the potential locations for the new servers, from which we derive a novel algorithm for MinMax, and show that it is efficient when the number of new servers is small. When the number of new servers is large, we propose an efficient 3-approximate algorithm. We verify with experiments on real road networks that our solutions are effective and attains significantly better result quality compared to the existing greedy algorithms.


title: Neural Architectures for Fine-grained Entity Type Classification

abstract: In this work, we investigate several neural network architectures for fine-grained entity type classification. Particularly, we consider extensions to a recently proposed attentive neural architecture and make three key contributions. Previous work on attentive neural architectures do not consider hand-crafted features, we combine learnt and hand-crafted features and observe that they complement each other. Additionally, through quantitative analysis we establish that the attention mechanism is capable of learning to attend over syntactic heads and the phrase containing the mention, where both are known strong hand-crafted features for our task. We enable parameter sharing through a hierarchical label encoding method, that in low-dimensional projections show clear clusters for each type hierarchy. Lastly, despite using the same evaluation dataset, the literature frequently compare models trained using different data. We establish that the choice of training data has a drastic impact on performance, with decreases by as much as 9.85% loose micro F1 score for a previously proposed method. Despite this, our best model achieves state-of-the-art results with 75.36% loose micro F1 score on the well- established FIGER (GOLD) dataset.


title: Recoverable robust spanning tree problem under interval uncertainty representations

abstract: This paper deals with the recoverable robust spanning tree problem under interval uncertainty representations. A polynomial time, combinatorial algorithm for the recoverable spanning tree problem is first constructed. This problem generalizes the incremental spanning tree problem, previously discussed in literature. The algorithm built is then applied to solve the recoverable robust spanning tree problem, under the traditional interval uncertainty representation, in polynomial time. Moreover, the algorithm allows to obtain, under some mild assumptions about the uncertainty intervals,several approximation results for the recoverable robust spanning tree problem under the Bertsimas and Sim interval uncertainty representation and the interval uncertainty representation with a budget constraint.


title: The Zero-Coupon Rate Model for Derivatives Pricing

abstract: The aim of this paper is to present a dual-term structure model of interest rate derivatives in order to solve the two hardest problems in financial modeling: the exact volatility calibration of the entire swaption matrix, and the calculation of bucket vegas for structured products. The model takes a series of long-term zero-coupon rates as basic state variables that are driven directly by one or more Brownian motion. The model volatility is assigned in a matrix form with two terms. A numerical scheme for implementing the model has been developed in the paper. At the end, several examples have been given for the model calibration, the structured products pricing and the calculation of bucket vegas.


title: Proceedings 17th International Workshop on Refinement

abstract: We are proud to present the papers from the 17th Refinement Workshop, co-located with FM 2015 held in Oslo, Norway on June 22nd, 2015. Refinement is one of the cornerstones of a formal approach to software engineering: the process of developing a more detailed design or implementation from an abstract specification through a sequence of mathematically-based steps that maintain correctness with respect to the original specification. This 17th workshop continued a 20+ year tradition under the auspices of the British Computer Society (BCS) FACS special interest group. This is the third volume that has appeared as an EPTCS proceedings, and we would like to thank the editorial board (and in particular Rob van Glabbeek) for their help and cooperation in making this happen. The organisers would like to thank everyone: the authors, BCS-FACS, EPTCS, and the organisers of FM 2015 for their help in organising this workshop, the participants of the workshop, and the reviewers involved in selecting the papers.


title: Building blocks of polarized endomorphisms of normal projective varieties

abstract: An endomorphism $f$ of a projective variety X is polarized (resp. quasi-polarized) if $f^*H$ is linearly equivalent to $qH$ for some ample (resp. nef and big) Cartier divisor $H$ and integer $q > 1$. First, we use cone analysis to show that a quasi-polarized endomorphism is always polarized, and the polarized property descends via any equivariant dominant rational map. Next, we show that a suitable maximal rationally connected fibration (MRC) can be made $f$-equivariant using a construction of N. Nakayama, that $f$ descends to a polarized endomorphism of the base Y of this MRC and that this Y is a Q-abelian variety (quasi-\'etale quotient of an abelian variety). Finally, we show that we can run the minimal model program (MMP) $f$-equivariantly for mildly singular X and reach either a Q-abelian variety or a Fano variety of Picard number one. As a consequence, the building blocks of polarized endomorphisms are those of Q-abelian varieties and those of Fano varieties of Picard number one. Along the way, we show that $f$ always descends to a polarized endomorphism of the Albanese variety Alb(X) of X, and that the pullback of a power of $f$ acts as a scalar multiplication on the Neron-Severi group of X (modulo torsion) when X is smooth and rationally connected. Partial answers about X being of Calabi-Yau type, or Fano type are also given with an extra primitivity assumption on $f$ which seems necessary by an example.


title: On commuting ordinary differential operators with polynomial coefficients corresponding to spectral curves of genus two

abstract: The group of automorphisms of the first Weyl algebra acts on commuting ordinary differential operators with polynomial coefficient. In this paper we prove that for fixed generic spectral curve of genus two the set of orbits is infinite.


title: Rotational control of asymmetric molecules: dipole- vs. polarizability- driven rotational dynamics

abstract: We experimentally study the optical- and terahertz- induced rotational dynamics of asymmetric molecules in the gas phase. Terahertz and optical fields are identified as two distinct control handles over asymmetric molecules, as they couple to the rotational degrees of freedom via the molecular- dipole and polarizability selectively. The distinction between those two rotational handles is highlighted by different types of quantum revivals observed in long duration (>100ps) field-free rotational evolution. The experimental results are in excellent agreement with Random Phase Wave Function simulations [Phys. Rev. A 91, 063420 (2015)] and provide verification of the RPWF as an efficient method for calculating asymmetric molecular dynamics at ambient temperatures, where exact calculation methods are practically not feasible. Our observations and analysis pave the way for orchestrated excitations by both optical and THz fields as complementary rotational handles, that enable a plethora of new possibilities in three-dimensional rotational control of asymmetric molecules.


title: A Conceptual Framework to Assess the Effectiveness of Rubric Tool

abstract: Rubrics are being used in a wide variety of disciplines in higher education to evaluate assessments and provide feedback to students. Rubrics are traditionally implemented as paper-based table format to grade assessments and provide feedback. The advancement of technology has integrated Rubric Tool into Learning Management Systems (LMS) like Blackboard to facilitate online marking. In this study, we aim to evaluate the effectiveness of Rubric Tool. We propose an integrated conceptual framework using various theories from extant literature to assess the effectiveness of Rubric Tool. Further, this framework has been applied in an undergraduate course at the Asian campus of an Australian International University. This is a research in progress paper which adds to the literature on integrating information systems theory, learning theory, political philosophy and communication theory to examine the effectiveness of Rubric Tool in student learning.