title: An alternative definition of cobordism map of ECH

abstract: In this article, we reformulate the cobordism map of embedded contact homology, which is induced by exact symplectic cobordism and defined as direct limit of homomorphisms called filtered ECH cobordism map. The filtered ECH cobordism map is defined by counting embedded holomorphic curves with zero ECH index and we prove that it is independent on almost complex structure by Seiberg Witten theory. Moreover, our definition in fact is equivalent to the existing definition.


title: A Visual Measure of Changes to Weighted Self-Organizing Map Patterns

abstract: Estimating output changes by input changes is the main task in causal analysis. In previous work, input and output Self-Organizing Maps (SOMs) were associated for causal analysis of multivariate and nonlinear data. Based on the association, a weight distribution of the output conditional on a given input was obtained over the output map space. Such a weighted SOM pattern of the output changes when the input changes. In order to analyze the change, it is important to measure the difference of the patterns. Many methods have been proposed for the dissimilarity measure of patterns. However, it remains a major challenge when attempting to measure how the patterns change. In this paper, we propose a visualization approach that simplifies the comparison of the difference in terms of the pattern property. Using this approach, the change can be analyzed by integrating colors and star glyph shapes representing the property dissimilarity. Ecological data is used to demonstrate the usefulness of our approach and the experimental results show that our approach provides the change information effectively.


title: On the geometry of orbifold Gromov-Witten invariants

abstract: We consider the question of how geometric structures of a Deligne-Mumford stack affect its Gromov-Witten invariants. The two geometric structures studied here are {\em gerbes} and {\em root constructions}. In both cases, we explain conjectures on Gromov-Witten theory for these stacks and survey some recent progress on these conjectures.


title: MIHash: Online Hashing with Mutual Information

abstract: Learning-based hashing methods are widely used for nearest neighbor retrieval, and recently, online hashing methods have demonstrated good performance-complexity trade-offs by learning hash functions from streaming data. In this paper, we first address a key challenge for online hashing: the binary codes for indexed data must be recomputed to keep pace with updates to the hash functions. We propose an efficient quality measure for hash functions, based on an information-theoretic quantity, mutual information, and use it successfully as a criterion to eliminate unnecessary hash table updates. Next, we also show how to optimize the mutual information objective using stochastic gradient descent. We thus develop a novel hashing method, MIHash, that can be used in both online and batch settings. Experiments on image retrieval benchmarks (including a 2.5M image dataset) confirm the effectiveness of our formulation, both in reducing hash table recomputations and in learning high-quality hash functions.


title: Modeling high dimensional multichannel brain signals

abstract: Our goal is to model and measure functional and effective (directional) connectivity in multichannel brain physiological signals (e.g., electroencephalograms, local field potentials). The difficulties from analyzing these data mainly come from two aspects: first, there are major statistical and computational challenges for modeling and analyzing high dimensional multichannel brain signals; second, there is no set of universally-agreed measures for characterizing connectivity. To model multichannel brain signals, our approach is to fit a vector autoregressive (VAR) model with potentially high lag order so that complex lead-lag temporal dynamics between the channels can be captured. Estimates of the VAR model will be obtained by our proposed hybrid LASSLE (LASSO+LSE) method which combines regularization (to control for sparsity) and least squares estimation (to improve bias and mean-squared error). Then we employ some measures of connectivity but put an emphasis on partial directed coherence (PDC) which can capture the directional connectivity between channels. PDC is a frequency-specific measure that explains the extent to which the present oscillatory activity in a sender channel influences the future oscillatory activity in a specific receiver channel relative to all possible receivers in the network. The proposed modeling approach provided key insights into potential functional relationships among simultaneously recorded sites during performance of a complex memory task. Specifically, this novel method was successful in quantifying patterns of effective connectivity across electrode locations, and in capturing how these patterns varied across trial epochs and trial types.


title: Nuclear physics insights for new-physics searches using nuclei: Neutrinoless $\beta\beta$ decay and dark matter direct detection

abstract: Experiments using nuclei to probe new physics beyond the Standard Model, such as neutrinoless $\beta\beta$ decay searches testing whether neutrinos are their own antiparticle, and direct detection experiments aiming to identify the nature of dark matter, require accurate nuclear physics input for optimizing their discovery potential and for a correct interpretation of their results. This demands a detailed knowledge of the nuclear structure relevant for these processes. For instance, neutrinoless $\beta\beta$ decay nuclear matrix elements are very sensitive to the nuclear correlations in the initial and final nuclei, and the spin-dependent nuclear structure factors of dark matter scattering depend on the subtle distribution of the nuclear spin among all nucleons. In addition, nucleons are composite and strongly interacting, which implies that many-nucleon processes are necessary for a correct description of nuclei and their interactions. It is thus crucial that theoretical studies and experimental analyses consider $\beta$ decays and dark matter interactions with a coupling to two nucleons, called two-nucleon currents.

categories: nucl-th astro-ph.CO hep-ex hep-ph nucl-ex


title: On Automating the Doctrine of Double Effect

abstract: The doctrine of double effect ($\mathcal{DDE}$) is a long-studied ethical principle that governs when actions that have both positive and negative effects are to be allowed. The goal in this paper is to automate $\mathcal{DDE}$. We briefly present $\mathcal{DDE}$, and use a first-order modal logic, the deontic cognitive event calculus, as our framework to formalize the doctrine. We present formalizations of increasingly stronger versions of the principle, including what is known as the doctrine of triple effect. We then use our framework to simulate successfully scenarios that have been used to test for the presence of the principle in human subjects. Our framework can be used in two different modes: One can use it to build $\mathcal{DDE}$-compliant autonomous systems from scratch, or one can use it to verify that a given AI system is $\mathcal{DDE}$-compliant, by applying a $\mathcal{DDE}$ layer on an existing system or model. For the latter mode, the underlying AI system can be built using any architecture (planners, deep neural networks, bayesian networks, knowledge-representation systems, or a hybrid); as long as the system exposes a few parameters in its model, such verification is possible. The role of the $\mathcal{DDE}$ layer here is akin to a (dynamic or static) software verifier that examines existing software modules. Finally, we end by presenting initial work on how one can apply our $\mathcal{DDE}$ layer to the STRIPS-style planning model, and to a modified POMDP model.This is preliminary work to illustrate the feasibility of the second mode, and we hope that our initial sketches can be useful for other researchers in incorporating DDE in their own frameworks.


title: Pseudocomplementation and Minimal Prime Ideals in Semirings

abstract: In the first section of the present work, we introduce the concept of pseudocomplementation for semirings and show semiring version of some known results in lattice theory. We also introduce semirings with pc-functions and prove some interesting results for minimal prime ideals of such semirings. In the second section, some classical results for minimal prime ideals in ring theory are generalized in the context of semiring theory.


title: Heat Kernels for Non-symmetric Non-local Operators

abstract: We survey the recent progress in the study of heat kernels for a class of non-symmetric non-local operators. We focus on the existence and sharp two-sided estimates of the heat kernels and their connection to jump diffusions.


title: Fringes' Impacts to Astrometry and Photometry of Stars

abstract: Fringes often appear in a CCD frame, especially when a thin CCD chip and a R or I filter is used. 88 CCD frames of the two open clusters NGC 2324 and NGC 1664 with a Johnson I filter taken from the 2.4-m telescope at Yunnan Observatory are used to study the fringes' impacts to the astrometry and photometry of stars. A novel technique proposed by Snodgrass & Carry is applied to remove the fringes in each CCD frame. And an appraisal of this technique is performed to estimate fringes' effects on astrometry and photometry of stars. Our results show that the astrometric and photometric precisions of stars can be improved effectively after the removal of fringes, especially for faint stars.


title: Optical and structural properties of $\mathrm{Eu^{2+}}$ doped BaBrI and BaClI crystals

abstract: The work is necessitated by search for new materials to detect ionizing radiation. The rare-earth ions doped with ternary alkali earth-halide systems are promising scintillators showing high efficiency and energy resolution. Some aspects of crystal growth and data on the structural and luminescence properties of BaBrI and BaClI doped with low concentrations of $\mathrm{Eu^{2+}}$ ions are reported. The crystals are grown by the vertical Bridgman method in sealed quartz ampoule. New crystallography data for BaClI single crystal obtained by single crystal X-ray diffraction method are presented in this paper. Emission, excitation and optical absorption spectra as well as luminescence decay kinetics are studied under excitation by X-ray, vacuum ultraviolet and ultraviolet radiation. The energies of the first 4f-5d transition in $\mathrm{Eu^{2+}}$ and band gap of the crystals have been obtained. We have calculated the electronic band structure of the crystals using density functional theory as implemented in the \latin{Ab Initio}. Calculated band gap energies are in accord with the experimental estimates. The energy of gaps between the occupied Eu$^{2+}$ 4f level and the valence band top are predicted. In addition, positions of lanthanide energy levels in relation to valence band have been constructed using the chemical shift model.


title: Manifesting enhanced cancellations in supergravity: integrands versus integrals

abstract: Examples of "enhanced ultraviolet cancellations" with no known standard-symmetry explanation have been found in a variety of supergravity theories. By examining one- and two-loop examples in four- and five-dimensional half-maximal supergravity, we argue that enhanced cancellations in general cannot be exhibited prior to integration. In light of this, we explore reorganizations of integrands into parts that are manifestly finite and parts that have poor power counting but integrate to zero due to integral identities. At two loops we find that in the large loop-momentum limit the required integral identities follow from Lorentz and SL(2) relabeling symmetry. We carry out a nontrivial check at four loops showing that the identities generated in this way are a complete set. We propose that at $L$ loops the combination of Lorentz and SL($L$) symmetry is sufficient for displaying enhanced cancellations when they happen, whenever the theory is known to be ultraviolet finite up to $(L-1)$ loops.


title: Resource-monotonicity and Population-monotonicity in Connected Cake-cutting

abstract: In the classic cake-cutting problem (Steinhaus, 1948), a heterogeneous resource has to be divided among n agents with different valuations in a proportional way --- giving each agent a piece with a value of at least 1/n of the total. In many applications, such as dividing a land-estate or a time-interval, it is also important that the pieces are connected. We propose two additional requirements: resource-monotonicity (RM) and population-monotonicity (PM). When either the cake or the set of agents changes and the cake is re-divided using the same rule, the utility of all remaining agents must change in the same direction. Classic cake-cutting protocols are neither RM nor PM. Moreover, we prove that no Pareto-optimal proportional division rule can be either RM or PM. Motivated by this negative result, we search for division rules that are weakly-Pareto-optimal --- no other division is strictly better for all agents. We present two such rules. The relative-equitable rule, which assigns the maximum possible relative value equal for all agents, is proportional and PM. The so-called rightmost-mark rule, which is an improved version of the Cut and Choose protocol, is proportional and RM for two agents.


title: A perturbation theory approach to the stability of the Pais-Uhlenbeck oscillator

abstract: We present a detailed analysis of the orbital stability of the Pais-Uhlenbeck oscillator, using Lie-Deprit series and Hamiltonian normal form theories. In particular, we explicitly describe the reduced phase space for this Hamiltonian system and give a proof for the existence of stable orbits for a certain class of self-interaction, found numerically in previous works, by using singular symplectic reduction.


title: Alternative Modes of Interaction in Proximal Human-in-the-Loop Operation of Robots

abstract: Ambiguity and noise in natural language instructions create a significant barrier towards adopting autonomous systems into safety critical workflows involving humans and machines. In this paper, we propose to build on recent advances in electrophysiological monitoring methods and augmented reality technologies, to develop alternative modes of communication between humans and robots involved in large-scale proximal collaborative tasks. We will first introduce augmented reality techniques for projecting a robot's intentions to its human teammate, who can interact with these cues to engage in real-time collaborative plan execution with the robot. We will then look at how electroencephalographic (EEG) feedback can be used to monitor human response to both discrete events, as well as longer term affective states while execution of a plan. These signals can be used by a learning agent, a.k.a an affective robot, to modify its policy. We will present an end-to-end system capable of demonstrating these modalities of interaction. We hope that the proposed system will inspire research in augmenting human-robot interactions with alternative forms of communications in the interests of safety, productivity, and fluency of teaming, particularly in engineered settings such as the factory floor or the assembly line in the manufacturing industry where the use of such wearables can be enforced.


title: Palindromic Decompositions with Gaps and Errors

abstract: Identifying palindromes in sequences has been an interesting line of research in combinatorics on words and also in computational biology, after the discovery of the relation of palindromes in the DNA sequence with the HIV virus. Efficient algorithms for the factorization of sequences into palindromes and maximal palindromes have been devised in recent years. We extend these studies by allowing gaps in decompositions and errors in palindromes, and also imposing a lower bound to the length of acceptable palindromes. We first present an algorithm for obtaining a palindromic decomposition of a string of length n with the minimal total gap length in time O(n log n * g) and space O(n g), where g is the number of allowed gaps in the decomposition. We then consider a decomposition of the string in maximal \delta-palindromes (i.e. palindromes with \delta errors under the edit or Hamming distance) and g allowed gaps. We present an algorithm to obtain such a decomposition with the minimal total gap length in time O(n (g + \delta)) and space O(n g).


title: Radiative-transfer models for explosions from rotating and non-rotating single WC stars. Implications for SN1998bw and LGRB/SNe

abstract: Using 1-D non-Local-Thermodynamic-Equilibrium time-dependent radiative-transfer simulations, we study the ejecta properties required to match the early and late-time photometric and spectroscopic properties of supernovae (SNe) associated with long-duration gamma-ray bursts (LGRBs). To match the short rise time, narrow light curve peak, and extremely broad spectral lines of SN1998bw requires a model with <3Msun ejecta but a high explosion energy of a few 1e52erg and 0.5Msun of Ni56. However the relatively high luminosity, the presence of narrow spectral lines of intermediate mass elements, and the low ionization at the nebular stage are matched with a more standard C-rich Wolf-Rayet (WR) star explosion, with an ejecta of >10Msun, an explosion energy >1e51erg, and only 0.1Msun of Ni56. As the two models are mutually exclusive, the breaking of spherical symmetry is essential to match the early/late photometric/spectroscopic properties of SN1998bw. This conclusion confirms the notion that the ejecta of SN1998bw is aspherical on large scales. More generally, with asphericity, the energetics and Ni56 mass of LGRB/SNe are reduced and their ejecta mass is increased, favoring a massive fast-rotating Wolf-Rayet star progenitor. Contrary to persisting claims in favor of the proto-magnetar model for LGRB/SNe, such progenitor/ejecta properties are compatible with collapsar formation. Ejecta properties of LGRB/SNe inferred from 1D radiative-transfer modeling are fundamentally flawed.


title: Multiple Instance Learning with the Optimal Sub-Pattern Assignment Metric

abstract: Multiple instance data are sets or multi-sets of unordered elements. Using metrics or distances for sets, we propose an approach to several multiple instance learning tasks, such as clustering (unsupervised learning), classification (supervised learning), and novelty detection (semi-supervised learning). In particular, we introduce the Optimal Sub-Pattern Assignment metric to multiple instance learning so as to provide versatile design choices. Numerical experiments on both simulated and real data are presented to illustrate the versatility of the proposed solution.


title: A circuit-preserving mapping from multilevel to Boolean dynamics

abstract: Many discrete models of biological networks rely exclusively on Boolean variables and many tools and theorems are available for analysis of strictly Boolean models. However, multilevel variables are often required to account for threshold effects, in which knowledge of the Boolean case does not generalise straightforwardly. This motivated the development of conversion methods for multilevel to Boolean models. In particular, Van Ham's method has been shown to yield a one-to-one, neighbour and regulation preserving dynamics, making it the de facto standard approach to the problem. However, Van Ham's method has several drawbacks: most notably, it introduces vast regions of "non-admissible" states that have no counterpart in the multilevel, original model. This raises special difficulties for the analysis of interaction between variables and circuit functionality, which is believed to be central to the understanding of dynamic properties of logical models. Here, we propose a new multilevel to Boolean conversion method, with software implementation. Contrary to Van Ham's, our method doesn't yield a one-to-one transposition of multilevel trajectories, however, it maps each and every Boolean state to a specific multilevel state, thus getting rid of the non-admissible regions and, at the expense of (apparently) more complicated, "parallel" trajectories. One of the prominent features of our method is that it preserves dynamics and interaction of variables in a certain manner. As a demonstration of the usability of our method, we apply it to construct a new Boolean counter-example to the well-known conjecture that a local negative circuit is necessary to generate sustained oscillations. This result illustrates the general relevance of our method for the study of multilevel logical models.


title: Security Constrained Multi-Stage Transmission Expansion Planning Considering a Continuously Variable Series Reactor

abstract: This paper introduces a Continuously Variable Series Reactor (CVSR) to the transmission expansion planning (TEP) problem. The CVSR is a FACTS-like device which has the capability of controlling the overall impedance of the transmission line. However, the cost of the CVSR is about one tenth of a similar rated FACTS device which potentially allows large numbers of devices to be installed. The multi-stage TEP with the CVSR considering the $N-1$ security constraints is formulated as a mixed integer linear programming model. The nonlinear part of the power flow introduced by the variable reactance is linearized by a reformulation technique. To reduce the computational burden for a practical large scale system, a decomposition approach is proposed. The detailed simulation results on the IEEE 24-bus and a more practical Polish 2383-bus system demonstrate the effectiveness of the approach. Moreover, the appropriately allocated CVSRs add flexibility to the TEP problem and allow reduced planning costs. Although the proposed decomposition approach cannot guarantee global optimality, a high level picture of how the network can be planned reliably and economically considering CVSR is achieved.


title: Timed Automata with Polynomial Delay and their Expressiveness

abstract: We consider previous models of Timed, Probabilistic and Stochastic Timed Automata, we introduce our model of Timed Automata with Polynomial Delay and we characterize the expressiveness of these models relative to each other.


title: A Scale Free Algorithm for Stochastic Bandits with Bounded Kurtosis

abstract: Existing strategies for finite-armed stochastic bandits mostly depend on a parameter of scale that must be known in advance. Sometimes this is in the form of a bound on the payoffs, or the knowledge of a variance or subgaussian parameter. The notable exceptions are the analysis of Gaussian bandits with unknown mean and variance by Cowan and Katehakis [2015] and of uniform distributions with unknown support [Cowan and Katehakis, 2015]. The results derived in these specialised cases are generalised here to the non-parametric setup, where the learner knows only a bound on the kurtosis of the noise, which is a scale free measure of the extremity of outliers.


title: Edge states in dynamical superlattices

abstract: We address edge states and rich localization regimes available in the one-dimensional (1D) dynamically modulated superlattices, both theoretically and numerically. In contrast to conventional lattices with straight waveguides, the quasi-energy band of infinite modulated superlattice is periodic not only in the transverse Bloch momentum, but it also changes periodically with increase of the coupling strength between waveguides. Due to collapse of quasi-energy bands dynamical superlattices admit known dynamical localization effect. If, however, such a lattice is truncated, periodic longitudinal modulation leads to appearance of specific edge states that exist within certain periodically spaced intervals of coupling constants. We discuss unusual transport properties of such truncated superlattices and illustrate different excitation regimes and enhanced robustness of edge states in them, that is associated with topology of the quasi-energy band.


title: Movement of time-delayed hot spots in Euclidean space for special initial states

abstract: We consider the Cauchy problem for the damped wave equation under the initial state that the sum of an initial position and an initial velocity vanishes. When the initial position is non-zero, non-negative and compactly supported, we study the large time behavior of the spatial null, critical, maximum and minimum sets of the solution. The spatial null set becomes a smooth hyper-surface homeomorphic to a sphere after a large enough time. The spatial critical set has at least three points after a large enough time. The set of spatial maximum points escapes from the convex hull of the support of the initial position. The set of spatial minimum points consists of one point after a large time, and the unique spatial minimum point converges to the centroid of the initial position at time infinity.


title: Tree Edit Distance Cannot be Computed in Strongly Subcubic Time (unless APSP can)

abstract: The edit distance between two rooted ordered trees with $n$ nodes labeled from an alphabet~$\Sigma$ is the minimum cost of transforming one tree into the other by a sequence of elementary operations consisting of deleting and relabeling existing nodes, as well as inserting new nodes. Tree edit distance is a well known generalization of string edit distance. The fastest known algorithm for tree edit distance runs in cubic $O(n^3)$ time and is based on a similar dynamic programming solution as string edit distance. In this paper we show that a truly subcubic $O(n^{3-\varepsilon})$ time algorithm for tree edit distance is unlikely: For $|\Sigma| = \Omega(n)$, a truly subcubic algorithm for tree edit distance implies a truly subcubic algorithm for the all pairs shortest paths problem. For $|\Sigma| = O(1)$, a truly subcubic algorithm for tree edit distance implies an $O(n^{k-\varepsilon})$ algorithm for finding a maximum weight $k$-clique. Thus, while in terms of upper bounds string edit distance and tree edit distance are highly related, in terms of lower bounds string edit distance exhibits the hardness of the strong exponential time hypothesis [Backurs, Indyk STOC'15] whereas tree edit distance exhibits the hardness of all pairs shortest paths. Our result provides a matching conditional lower bound for one of the last remaining classic dynamic programming problems.


title: Physical Layer Security in Wireless Ad Hoc Networks Under A Hybrid Full-/Half-Duplex Receiver Deployment Strategy

abstract: This paper studies physical layer security in a wireless ad hoc network with numerous legitimate transmitter-receiver pairs and eavesdroppers. A hybrid full-/half-duplex receiver deployment strategy is proposed to secure legitimate transmissions, by letting a fraction of legitimate receivers work in the full-duplex (FD) mode sending jamming signals to confuse eavesdroppers upon their information receptions, and letting the other receivers work in the half-duplex mode just receiving their desired signals. The objective of this paper is to choose properly the fraction of FD receivers for achieving the optimal network security performance. Both accurate expressions and tractable approximations for the connection outage probability and the secrecy outage probability of an arbitrary legitimate link are derived, based on which the area secure link number, network-wide secrecy throughput and network-wide secrecy energy efficiency are optimized respectively. Various insights into the optimal fraction are further developed and its closed-form expressions are also derived under perfect self-interference cancellation or in a dense network. It is concluded that the fraction of FD receivers triggers a non-trivial trade-off between reliability and secrecy, and the proposed strategy can significantly enhance the network security performance.


title: Halo-independent determination of the unmodulated WIMP signal in DAMA: the isotropic case

abstract: We present a halo-independent determination of the unmodulated signal corresponding to the DAMA modulation if interpreted as due to dark matter weakly interacting massive particles (WIMPs). First we show how a modulated signal gives information on the WIMP velocity distribution function in the Galactic rest frame, from which the unmodulated signal descends. Then we perform a mathematically-sound profile likelihood analysis in which we profile the likelihood over a continuum of nuisance parameters (namely, the WIMP velocity distribution). As a first application of the method, which is very general and valid for any class of velocity distributions, we restrict the analysis to velocity distributions that are isotropic in the Galactic frame. In this way we obtain halo-independent maximum-likelihood estimates and confidence intervals for the DAMA unmodulated signal. We find that the estimated unmodulated signal is in line with expectations for a WIMP-induced modulation and is compatible with the DAMA background+signal rate. Specifically, for the isotropic case we find that the modulated amplitude ranges between a few percent and about 25% of the unmodulated amplitude, depending on the WIMP mass.


title: Equivalence of recurrence and Liouville property for symmetric Dirichlet forms

abstract: Given a symmetric Dirichlet form $(\mathcal{E},\mathcal{F})$ on a (non-trivial) $\sigma$-finite measure space $(E,\mathcal{B},m)$ with associated Markovian semigroup $\{T_{t}\}_{t\in(0,\infty)}$, we prove that $(\mathcal{E},\mathcal{F})$ is both irreducible and recurrent if and only if there is no non-constant $\mathcal{B}$-measurable function $u:E\to[0,\infty]$ that is \emph{$\mathcal{E}$-excessive}, i.e., such that $T_{t}u\leq u$ $m$-a.e.\ for any $t\in(0,\infty)$. We also prove that these conditions are equivalent to the equality $\{u\in\mathcal{F}_{e}\mid \mathcal{E}(u,u)=0\}=\mathbb{R}\mathbf{1}$, where $\mathcal{F}_{e}$ denotes the extended Dirichlet space associated with $(\mathcal{E},\mathcal{F})$. The proof is based on simple analytic arguments and requires no additional assumption on the state space or on the form. In the course of the proof we also present a characterization of the $\mathcal{E}$-excessiveness in terms of $\mathcal{F}_{e}$ and $\mathcal{E}$, which is valid for any symmetric positivity preserving form.


title: Intelligent bidirectional rapidly-exploring random trees for optimal motion planning in complex cluttered environments

abstract: The sampling based motion planning algorithm known as Rapidly-exploring Random Trees (RRT) has gained the attention of many researchers due to their computational efficiency and effectiveness. Recently, a variant of RRT called RRT* has been proposed that ensures asymptotic optimality. Subsequently its bidirectional version has also been introduced in the literature known as Bidirectional-RRT* (B-RRT*). We introduce a new variant called Intelligent Bidirectional-RRT* (IB-RRT*) which is an improved variant of the optimal RRT* and bidirectional version of RRT* (B-RRT*) algorithms and is specially designed for complex cluttered environments. IB-RRT* utilizes the bidirectional trees approach and introduces intelligent sample insertion heuristic for fast convergence to the optimal path solution using uniform sampling heuristics. The proposed algorithm is evaluated theoretically and experimental results are presented that compares IB-RRT* with RRT* and B-RRT*. Moreover, experimental results demonstrate the superior efficiency of IB-RRT* in comparison with RRT* and B-RRT in complex cluttered environments.


title: Uniform description of the rigged configuration bijection

abstract: We give a uniform description of the bijection $\Phi$ from rigged configurations to tensor products of Kirillov--Reshetikhin crystals of the form $\bigotimes_{i=1}^N B^{r_i,1}$ in dual untwisted types: simply-laced types and types $A_{2n-1}^{(2)}$, $D_{n+1}^{(2)}$, $E_6^{(2)}$, and $D_4^{(3)}$. We give a uniform proof that $\Phi$ is a bijection and preserves statistics. We describe $\Phi$ uniformly using virtual crystals for all remaining types, but our proofs are type-specific. We also give a uniform proof that $\Phi$ is a bijection for $\bigotimes_{i=1}^N B^{r_i,s_i}$ when $r_i$, for all $i$, map to $0$ under an automorphism of the Dynkin diagram. Furthermore, we give a description of the Kirillov--Reshetikhin crystals $B^{r,1}$ using tableaux of a fixed height $k_r$ depending on $r$ in all affine types. Additionally, we are able to describe crystals $B^{r,s}$ using $k_r \times s$ shaped tableaux that are conjecturally the crystal basis for Kirillov--Reshetikhin modules for various nodes $r$.


title: The most effective model for describing the universal behavior of unstable surface growth

abstract: We study a noisy Kuramoto-Sivashinsky (KS) equation which describes unstable surface growth and chemical turbulence. It has been conjectured that the universal long-wavelength behavior of the equation, which is characterized by scale-dependent parameters, is described by a Kardar-Parisi-Zhang (KPZ) equation. We consider this conjecture by analyzing a renormalization-group equation for a class of generalized KPZ equations. We then uniquely determine the parameter values of the KPZ equation that most effectively describes the universal long-wavelength behavior of the noisy KS equation.


title: In Vivo Evaluation of the Secure Opportunistic Schemes Middleware using a Delay Tolerant Social Network

abstract: Over the past decade, online social networks (OSNs) such as Twitter and Facebook have thrived and experienced rapid growth to over 1 billion users. A major evolution would be to leverage the characteristics of OSNs to evaluate the effectiveness of the many routing schemes developed by the research community in real-world scenarios. In this paper, we showcase the Secure Opportunistic Schemes (SOS) middleware which allows different routing schemes to be easily implemented relieving the burden of security and connection establishment. The feasibility of creating a delay tolerant social network is demonstrated by using SOS to power AlleyOop Social, a secure delay tolerant networking research platform that serves as a real-life mobile social networking application for iOS devices. SOS and AlleyOop Social allow users to interact, publish messages, and discover others that share common interests in an intermittent network using Bluetooth, peer-to-peer WiFi, and infrastructure WiFi.


title: Complex curves in pseudoconvex Runge domains containing discrete subsets

abstract: For any pseudoconvex Runge domain $\Omega\subset\mathbb{C}^2$ we prove that every closed discrete subset in $\Omega$ is contained in a properly embedded complex curve in $\Omega$ with any prescribed topology (possibly infinite).


title: Spin-incoherent Luttinger liquid of one-dimensional spin-1 Tonks-Girardeau Bose gas: Spin-dependent properties

abstract: Spin-incoherent Luttinger liquid (SILL) is a different universal class from the Luttinger liquid.\ This difference results from the spin incoherence of the system when the thermal energy of the system is higher than the spin excitation energy.\ We consider one-dimensional spin-$1$ Bose gas in the SILL regime and investigate its spin-dependent many-body properties.\ In Tonks-Girardeau limit, we are able to write down the general wave functions in a harmonic trap.\ We numerically calculate the spin-dependent (spin-plus, minus, and $0$) momentum distributions in the sector of zero magnetization which allows to demonstrate the most significant spin-incoherent feature compared to the spinless or spin-polarized case.\ In contrast to the spinless Bose gas, the momentum distributions are broadened and in the large momentum limit follow the same asymptotic $1/p^4$ dependence but with reduced coefficients.\ While the density matrices and momentum distributions differ between different spin components for small $N$, at large $N$ they approach each other.\ We show these by analytic arguments and numerical calculations up to $N$ $=$ $16$.


title: Gene tree species tree reconciliation with gene conversion

abstract: Gene tree/species tree reconciliation is a recent decisive progress in phylo-genetic methods, accounting for the possible differences between gene histories and species histories. Reconciliation consists in explaining these differences by gene-scale events such as duplication, loss, transfer, which translates mathematically into a mapping between gene tree nodes and species tree nodes or branches. Gene conversion is a very frequent biological event, which results in the replacement of a gene by a copy of another from the same species and in the same gene tree. Including this event in reconciliations has never been attempted because this changes as well the solutions as the methods to construct reconciliations. Standard algorithms based on dynamic programming become ineffective. We propose here a novel mathematical framework including gene conversion as an evolutionary event in gene tree/species tree reconciliation. We describe a randomized algorithm giving in polynomial running time a reconciliation minimizing the number of duplications, losses and conversions. We show that the space of reconciliations includes an analog of the Last Common Ancestor reconciliation, but is not limited to it. Our algorithm outputs any optimal reconciliation with non null probability. We argue that this study opens a wide research avenue on including gene conversion in reconciliation, which can be important for biology.


title: Long-lasting Quantum Memories: Extending the Coherence Time of Superconducting Artificial Atoms in the Ultrastrong-Coupling Regime

abstract: Quantum systems are affected by interactions with their environments, causing decoherence through two processes: pure dephasing and energy relaxation. For quantum information processing it is important to increase the coherence time of Josephson qubits and other artificial two-level atoms. We show theoretically that if the coupling between these qubits and a cavity field is longitudinal and in the ultrastrong-coupling regime, the system is strongly protected against relaxation. Vice versa, if the coupling is transverse and in the ultrastrong-coupling regime, the system is protected against pure dephasing. Taking advantage of the relaxation suppression, we show that it is possible to enhance their coherence time and use these qubits as quantum memories. Indeed, to preserve the coherence from pure dephasing, we prove that it is possible to apply dynamical decoupling. We also use an auxiliary atomic level to store and retrieve quantum information.


title: Rotating hybrid stars with the Dyson-Schwinger quark model

abstract: We study rapidly rotating hybrid stars with the Dyson-Schwinger model for quark matter and the Brueckner-Hartree-Fock many-body theory with realistic two-body and three-body forces for nuclear matter. We determine the maximum gravitational mass, equatorial radius, and rotation frequency of stable stellar configurations by considering the constraints of the Keplerian limit and the secular axisymmetric instability, and compare with observational data. We also discuss the rotational evolution for constant baryonic mass, and find a spinup phenomenon for supramassive stars before they collapse to black holes.


title: Sharp estimates for the anisotropic torsional rigidity and the principal frequency

abstract: In this paper we generalize some classical estimates involving the torsional rigidity and the principal frequency of a convex domain to a class of functionals related to some anisotropic non linear operators.


title: Exact and approximate limit behaviour of the Yule tree's cophenetic index

abstract: In this work we study the limit distribution of an appropriately normalized cophenetic index of the pure-birth tree conditioned on $n$ contemporary tips. We show that this normalized phylogenetic balance index is a submartingale that converges almost surely and in L2. We link our work with studies on trees without branch lengths and show that in this case the limit distribution is a contraction-type distribution, similar to the Quicksort limit distribution. In the continuous branch case we suggest approximations to the limit distribution. We propose heuristic methods of simulating from these distributions and it may be observed that these algorithms result in reasonable tails. Therefore, we propose a way based on the quantiles of the derived distributions for hypothesis testing, whether an observed phylogenetic tree is consistent with the pure-birth process. Simulating a sample by the proposed heuristics is rapid, while exact simulation (simulating the tree and then calculating the index) is a time-consuming procedure. We conduct a power study to investigate how well the cophenetic indices detect deviations from the Yule tree and apply the methodology to empirical phylogenies.


title: Analysis of Different Approaches of Parallel Block Processing for K-Means Clustering Algorithm

abstract: Distributed Computation has been a recent trend in engineering research. Parallel Computation is widely used in different areas of Data Mining, Image Processing, Simulating Models, Aerodynamics and so forth. One of the major usage of Parallel Processing is widely implemented for clustering the satellite images of size more than dimension of 1000x1000 in a legacy system. This paper mainly focuses on the different approaches of parallel block processing such as row-shaped, column-shaped and square-shaped. These approaches are applied for classification problem. These approaches is applied to the K-Means clustering algorithm as this is widely used for the detection of features for high resolution orthoimagery satellite images. The different approaches are analyzed, which lead to reduction in execution time and resulted the influence of improvement in performance measurement compared to sequential K-Means Clustering algorithm.


title: Trace Expressiveness of Timed and Probabilistic Automata

abstract: Automata expressiveness is an essential feature in understanding which of the formalisms available should be chosen for modelling a particular problem. Probabilistic and stochastic automata are suitable for modelling systems exhibiting probabilistic behavior and their expressiveness has been studied relative to non-probabilistic transition systems and Markov chains. In this paper, we consider previous formalisms of Timed, Probabilistic and Stochastic Timed Automata, we present our new model of Timed Automata with Polynomial Delay, we introduce a measure of expressiveness for automata we call trace expressiveness and we characterize the expressiveness of these models relative to each other under this new measure.


title: Vibrations on pulse tube based Dry Dilution Refrigerators for low noise measurements

abstract: Dry Dilution Refrigerators (DDR) based on pulse tube cryo-coolers have started to replace Wet Dilution Refrigerators (WDR) due to the ease and low cost of operation. However these advantages come at the cost of increased vibrations, induced by the pulse tube. In this work, we present the vibration measurements performed on three different commercial DDRs. We describe in detail the vibration measurement system we assembled, based on commercial accelerometers, conditioner and DAQ, and examined the effects of the various damping solutions utilized on three different DDRs, both in the low and high frequency regions. Finally, we ran low temperature, pseudo-massive (30 and 250 g) germanium bolometers in the best vibration-performing system under study and report on the results.


title: Optimal insider control of stochastic Volterra equations

abstract: We study the problem of optimal inside control of a stochastic Volterra equation driven by a Brownian motion and a Poisson random measure. We prove a sufficient and a necessary maximum principle for the optimal control when the trader has only partial information available to her decisions and on the other hand, may have some inside information about the future of the system. The results are applied to the problem of finding the optimal insider portfolio in a financial market where the risky asset price is given by a stochastic Volterra equation.


title: Micromagnetic simulations of spin-torque driven magnetisation dynamics with spatially resolved spin transport and magnetisation texture

abstract: We present a simple and fast method to simulate spin-torque driven magnetisation dynamics in nano-pillar spin-valve structures. The approach is based on the coupling between a spin transport code based on random matrix theory and a micromagnetics finite-elements software. In this way the spatial dependence of both spin transport and magnetisation dynamics is properly taken into account. Our results are compared with experiments. The excitation of the spin-wave modes, in- cluding the threshold current for steady state magnetisation precession and the nonlinear frequency shift of the modes are reproduced correctly. The giant magneto resistance effect and the magnetisa- tion switching also agree with experiment. The similarities with recently described spin-caloritronics devices are also discussed.


title: High-energy exciton transitions in quasi-two-dimensional cadmium chalcogenide nanoplatelets

abstract: Semiconductor nanoparticles of cadmium chalcogenides are known to exhibit pronounced thickness-dependent $E_0$ series of exciton transitions at the $\Gamma$ point of the Brillouin zone (BZ). In this work, we report an experimental evidence for high-energy series of exciton transitions, which originates from BZ points different from the $\Gamma$ point, in the family of cadmium chalcogenide quasi-2D nanoplatelets (NPLs). Intensive UV absorption bands demonstrating a pronounced size effect are observed for CdTe, CdSe, and CdS NPLs in addition to the $E_0$ exciton bands in the visible region. These new bands are attributed to transitions analogous to the $E_1$, $E_1+\Delta_1$, and $E_2$ series observed in bulk crystals. First-principles DFT calculations of the electronic structure and absorption spectra support this explanation and show that the main contribution to these optical transitions comes from $X$ and $M$ points of the 2D BZ, which originate from $L$ and $X$ points of the 3D BZ. At the same time, the $E_0$ series of transitions at the $\Gamma$ point is well described by the multiband effective-mass model. The observation of the UV exciton bands reveals tunable optical properties of cadmium chalcogenide NPLs in UV spectral region, which may be interesting for practical applications.


title: Scaling the Scattering Transform: Deep Hybrid Networks

abstract: We use the scattering network as a generic and fixed ini-tialization of the first layers of a supervised hybrid deep network. We show that early layers do not necessarily need to be learned, providing the best results to-date with pre-defined representations while being competitive with Deep CNNs. Using a shallow cascade of 1 x 1 convolutions, which encodes scattering coefficients that correspond to spatial windows of very small sizes, permits to obtain AlexNet accuracy on the imagenet ILSVRC2012. We demonstrate that this local encoding explicitly learns invariance w.r.t. rotations. Combining scattering networks with a modern ResNet, we achieve a single-crop top 5 error of 11.4% on imagenet ILSVRC2012, comparable to the Resnet-18 architecture, while utilizing only 10 layers. We also find that hybrid architectures can yield excellent performance in the small sample regime, exceeding their end-to-end counterparts, through their ability to incorporate geometrical priors. We demonstrate this on subsets of the CIFAR-10 dataset and on the STL-10 dataset.


title: Entropy spikes as a signature of Lifshitz transitions in the Dirac materials

abstract: We demonstrate theoretically that the characteristic feature of a 2D system undergoing $N$ consequent Lifshitz topological transitions is the occurrence of spikes of entropy per particle $s$ of a magnitude $\pm \ln 2/(J-1/2)$ with $2 \leq J \leq N$ at low temperatures. We derive a general expression for $s$ as a function of chemical potential, temperature and gap magnitude for the gapped Dirac materials. Inside the smallest gap, the dependence of $s$ on the chemical potential exhibits a dip-and-peak structure in the temperature vicinity of the Dirac point. The spikes of the entropy per particles can be considered as a signature of the Dirac materials. These distinctive characteristics of gapped Dirac materials can be detected in transport experiments where the temperature is modulated in gated structures.


title: Numerical Study of Flow Structure and Pedestrian Level Wind Comfort Inside Urban Street Canyons

abstract: In this work we numerically investigate the flow conditions inside uniform and non-uniform street canyons well within the atmospheric boundary layer. The numerical simulations use the steady RANS method with the near-wall modelling approach to simulate wall roughness at the boundary. With the aim of investigating both flow structure in broad terms, and pedestrian comfort in the street canyon between parallel buildings, we test different canyon configurations with varied street width, building width and building height. Turbulent conditions are broadly expected to hold within the physically-realistic range of Reynolds number of order $10^6$ considered here, where we take the building height to be a characteristic length scale, and the free stream velocity as the characteristic velocity. In addition to discussing the features of the canyon and wake flow, we investigate the effects of canyon geometry on pedestrian comfort by using the Extended Land Beaufort Scale for this purpose. We present and compare pedestrian comfort maps for each of our geometries.


title: Variational path-integral approach to back-reactions of composite mesons in the Nambu-Jona-Lasinio model

abstract: For the investigation of back-reactions of composite mesons in the NJL model, a variational path-integral treatment is formulated which yields an effective action $\mathscr{A}_{\rm eff}[D_{\sigma},D_{\pi}; S]$, depending on the propagators $D_\sigma$, $D_\pi$ of $\sigma-$ and $\pi-$mesons and on the full quark propagator $S$. The stationarity conditions $\delta \mathscr{A}_{\rm eff}/ \delta S = 0$, $\delta \mathscr{A}_{\rm eff}/ \delta D_\sigma = 0$, $\delta \mathscr{A}_{\rm eff}/ \delta D_\pi = 0$, then lead to coupled Schwinger-Dyson (SD) equations for the quark self-energy and the meson polarization functions. These results reproduce and extend results of the so-called "$\Phi-$derivable" approach and provide a functional formulation for diagrammatic resummations of $1/N_c-$corrections in the NJL model. Finally, we perform a low-momentum estimate of the quark and meson loop contributions to the polarization function of the pion and on this basis discuss the Goldstone theorem.


title: Theory of non-Markovian dynamics in resonance fluorescence spectrum

abstract: We present a detailed theoretical study of non-Markovian dynamics in the fluorescence spectrum of a driven semiconductor quantum dot (QD), embedded in a cavity and coupled to a three-dimensional (3D) acoustic phonon reservoir. In particular, we investigate the effect of pure dephasing on one of the side-peaks of the Mollow-triplet spectrum, expressed in terms of the off-diagonal element of the reduced system operator. The QD is modeled as a two-level system with an excited state representing a single exciton, and ground state represents the absence of an exciton. Coupling to the radiative modes of the cavity is treated within usual Born-Markov approximation, whereas dot-phonon coupling is discussed within non-Markovian regime beyond Born approximation. Using an equation-of-motion technique, the dot-phonon coupling is solved exactly and found that the exact result coincides with that of obtained within Born approximation. Furthermore, a Markov approximation is carried out with respect to the phonon interaction and compared with the non- Markovian lineshape for different values of the phonon bath temperature. We have found that coupling to the phonons vanishes for a resonant pump laser. For a non-resonant pump, we have characterzied the effect of dot-laser detuning and temperature of the phonon bath on the lineshape. The sideband undergoes a distinct narrowing and aquires an asymmetric shape with increasing phonon bath temperature. We have explained this behavior using a dressed-state picture of the QD levels.


title: Mastering Sketching: Adversarial Augmentation for Structured Prediction

abstract: We present an integral framework for training sketch simplification networks that convert challenging rough sketches into clean line drawings. Our approach augments a simplification network with a discriminator network, training both networks jointly so that the discriminator network discerns whether a line drawing is a real training data or the output of the simplification network, which in turn tries to fool it. This approach has two major advantages. First, because the discriminator network learns the structure in line drawings, it encourages the output sketches of the simplification network to be more similar in appearance to the training sketches. Second, we can also train the simplification network with additional unsupervised data, using the discriminator network as a substitute teacher. Thus, by adding only rough sketches without simplified line drawings, or only line drawings without the original rough sketches, we can improve the quality of the sketch simplification. We show how our framework can be used to train models that significantly outperform the state of the art in the sketch simplification task, despite using the same architecture for inference. We additionally present an approach to optimize for a single image, which improves accuracy at the cost of additional computation time. Finally, we show that, using the same framework, it is possible to train the network to perform the inverse problem, i.e., convert simple line sketches into pencil drawings, which is not possible using the standard mean squared error loss. We validate our framework with two user tests, where our approach is preferred to the state of the art in sketch simplification 92.3% of the time and obtains 1.2 more points on a scale of 1 to 5.


title: Applications of the hyperbolic Ax-Schanuel conjecture

abstract: In 2014, Pila and Tsimerman gave a proof of the Ax-Schanuel conjecture for the $j$-function and, with Mok, have recently announced a proof of its generalization to any (pure) Shimura variety. We refer to this generalization as the hyperbolic Ax-Schanuel conjecture. In this article, we show that the hyperbolic Ax-Schanuel conjecture can be used to reduce the Zilber-Pink conjecture for Shimura varieties to a problem of point counting. We further show that this point counting problem can be tackled in a number of cases using the Pila-Wilkie counting theorem and several arithmetic conjectures. Our methods are inspired by previous applications of the Pila-Zannier method and, in particular, the recent proof by Habegger and Pila of the Zilber-Pink conjecture for curves in abelian varieties.


title: The normal closure of big Dehn twists, and plate spinning with rotating families

abstract: We study the normal closure of a big power of one or several Dehn twists in a Mapping Class Group. We prove that it has a presentation whose relators consists only of commutators between twists of disjoint support, thus answering a question of Ivanov. Our method is to use the theory of projection complexes of Bestvina Bromberg and Fujiwara, together with the theory of rotating families, simultaneously on several spaces.


title: Annexing magic and tune-out wavelengths to the clock transitions of the alkaline-earth metal ions

abstract: We present additional magic wavelengths ($\lambda_{\rm{magic}}$) for the clock transitions in the alkaline-earth metal ions considering circular polarized light aside from our previously reported values in [J. Kaur et al., Phys. Rev. A {\bf 92}, 031402(R) (2015)] for the linearly polarized light. Contributions from the vector component to the dynamic dipole polarizabilities ($\alpha_d(\omega)$) of the atomic states associated with the clock transitions play major roles in the evaluation of these $\lambda_{\rm{magic}}$, hence facilitating in choosing circular polarization of lasers in the experiments. Moreover, the actual clock transitions in these ions are carried out among the hyperfine levels. The $\lambda_{\rm{magic}}$ values in these hyperfine transitions are estimated and found to be different from $\lambda_{\rm{magic}}$ for the atomic transitions due to different contributions coming from the vector and tensor part of $\alpha_d(\omega)$. Importantly, we also present $\lambda_{\rm{magic}}$ values that depend only on the scalar component of $\alpha_d(\omega)$ for their uses in a specially designed trap geometry for these ions so that they can be used unambiguously among any hyperfine levels of the atomic states of the clock transitions. We also present $\alpha_d(\omega)$ values explicitly at the 1064 nm for the atomic states associated with the clock transitions which may be useful for creating "high-field seeking" traps for the above ions using the Nd:YAG laser. The tune out wavelengths at which the states would be free from the Stark shifts are also presented. Accurate values of the electric dipole matrix elements required for these studies are given and trends of electron correlation effects in determining them are also highlighted.


title: Multimodal deep learning approach for joint EEG-EMG data compression and classification

abstract: In this paper, we present a joint compression and classification approach of EEG and EMG signals using a deep learning approach. Specifically, we build our system based on the deep autoencoder architecture which is designed not only to extract discriminant features in the multimodal data representation but also to reconstruct the data from the latent representation using encoder-decoder layers. Since autoencoder can be seen as a compression approach, we extend it to handle multimodal data at the encoder layer, reconstructed and retrieved at the decoder layer. We show through experimental results, that exploiting both multimodal data intercorellation and intracorellation 1) Significantly reduces signal distortion particularly for high compression levels 2) Achieves better accuracy in classifying EEG and EMG signals recorded and labeled according to the sentiments of the volunteer.


title: BPS Alice strings

abstract: When a charged particle encircles around an Alice string, it changes the sign of the electric charge. In this paper we find a BPS-saturated Alice string in $U(1)\times SO(3)$ gauge theory with charged complex scalar fields belonging to the vector representation. After performing BPS completion we solve the BPS equations numerically. We further embed the Alice string into an ${\cal N}=1$ supersymmetric gauge theory to show that it is half BPS.


title: Thompson Sampling for Linear-Quadratic Control Problems

abstract: We consider the exploration-exploitation tradeoff in linear quadratic (LQ) control problems, where the state dynamics is linear and the cost function is quadratic in states and controls. We analyze the regret of Thompson sampling (TS) (a.k.a. posterior-sampling for reinforcement learning) in the frequentist setting, i.e., when the parameters characterizing the LQ dynamics are fixed. Despite the empirical and theoretical success in a wide range of problems from multi-armed bandit to linear bandit, we show that when studying the frequentist regret TS in control problems, we need to trade-off the frequency of sampling optimistic parameters and the frequency of switches in the control policy. This results in an overall regret of $O(T^{2/3})$, which is significantly worse than the regret $O(\sqrt{T})$ achieved by the optimism-in-face-of-uncertainty algorithm in LQ control problems.


title: Distribution of water in the G327.3-0.6 massive star-forming region

abstract: We aim at characterizing the large-scale distribution of H2O in G327.3-0.6, a massive star-forming region made of individual objects in different evolutionary phases. We investigate variations of H2O abundance as function of evolution. We present Herschel continuum maps at 89 and 179 $\mu$m of the whole region and an APEX map at 350 {\mu}m of the IRDC. New spectral HIFI maps toward the IRDC region covering low-energy H2O lines at 987 and 1113 GHz are also presented and combined with HIFI pointed observations of the G327 hot core. We infer the physical properties of the gas through optical depth analysis and radiative transfer modeling. The continuum emission at 89 and 179 {\mu}m follows the thermal continuum emission at longer wavelengths, with a peak at the position of the hot core, a secondary peak in the Hii region, and an arch-like layer of hot gas west of the Hii region. The same morphology is observed in the 1113 GHz line, in absorption toward all dust condensations. Optical depths of ~80 and 15 are estimated and correspond to column densities of 10^15 and 2 10^14 cm-2, for the hot core and IRDC position. These values indicate an H2O to H2 ratio of 3 10^-8 toward the hot core; the abundance of H2O does not change along the IRDC with values of some 10^-8. Infall (over ~ 20") is detected toward the hot core position with a rate of 1-1.3 10^-2 M_sun /yr, high enough to overcome the radiation pressure due to the stellar luminosity. The source structure of the hot core region is complex, with a cold outer gas envelope in expansion, situated between the outflow and the observer, extending over 0.32 pc. The outflow is seen face-on and centered away from the hot core. The distribution of H2O along the IRDC is roughly constant with an abundance peak in the more evolved object. These water abundances are in agreement with previous studies in other massive objects and chemical models.


title: Generalized image deconvolution by exploiting spatially variant point spread functions

abstract: An optical imaging system forms an object image by recollecting light scattered by the object. However, intact optical information of the object delivered through the imaging system is deteriorated by imperfect optical elements and unwanted defects. Image deconvolution, also known as inverse filtering, has been widely exploited as a recovery technique because of its practical feasibility, and operates by assuming the linear shift-invariant property of the imaging system. However, shift invariance is not rigorously hold in all imaging situations and it is not a necessary condition for solving the inverse problem of light propagation. Here, we present a method to solve the linear inverse problem of coherent light propagation without assuming shift invariance. Full characterization of imaging capability of the system is achieved by successively recording optical responses, using various laser illumination angles which are systematically controlled by a digital micro-mirror device. Experimental results show that image distortions caused by optical defocus can be restored by conventional deconvolution, but severe aberrations produced by a tilted lens or an inserted disordered layer can be corrected only by the proposed generalized image deconvolution. This work generalizes the theory of optical imaging and deconvolution, and enables distortion-free imaging under any general imaging condition.


title: Thrust vectoring of an electric solar wind sail with a realistic sail shape

abstract: The shape of a rotating electric solar wind sail under the centrifugal force and solar wind dynamic pressure is modeled to address the sail attitude maintenance and thrust vectoring. The sail rig assumes centrifugally stretched main tethers that extend radially outward from the spacecraft in the sail spin plane. Furthermore, the tips of the main tethers host remote units that are connected by auxiliary tethers at the sail rim. Here, we derive the equation of main tether shape and present both a numerical solution and an analytical approximation for the shape as parametrized both by the ratio of the electric sail force to the centrifugal force and the sail orientation with respect to the solar wind direction. The resulting shape is such that near the spacecraft, the roots of the main tethers form a cone, whereas towards the rim, this coning is flattened by the centrifugal force, and the sail is coplanar with the sail spin plane. Our approximation for the sail shape is parametrized only by the tether root coning angle and the main tether length. Using the approximate shape, we obtain the torque and thrust of the electric sail force applied to the sail. As a result, the amplitude of the tether voltage modulation required for the maintenance of the sail attitude is given as a torque-free solution. The amplitude is smaller than that previously obtained for a rigid single tether resembling a spherical pendulum. This implies that less thrusting margin is required for the maintenance of the sail attitude. For a given voltage modulation, the thrust vectoring is then considered in terms of the radial and transverse thrust components.


title: Hybrid Filtering for a Class of Quantum Systems with Classical Disturbances

abstract: A filtering problem for a class of quantum systems disturbed by a classical stochastic process is investigated in this paper. The classical disturbance process, which is assumed to be described by a linear stochastic differential equation, is modeled by a quantum cavity model. Then the hybrid quantum-classical system is described by a combined quantum system consisting of two quantum cavity subsystems. Quantum filtering theory and a quantum extended Kalman filter method are employed to estimate the states of the combined quantum system. An estimate of the classical stochastic process is derived from the estimate of the combined quantum system. The effectiveness and performance of the proposed methods are illustrated by numerical results.


title: The He atom revisited

abstract: In the Feynman-Kac[1] path integral approach the eigenvalues of a quantum system can be computed using Wiener measure which uses Brownian particle motion. In our previous work[2-3] on such systems we have observed that the Wiener process numerically converges slowly for dimensions greater than two because almost all trajectories will escape to infinity[4]. One can speed up this process by using a Generalized Feynman-Kac (GFK) method[5] in which the new measure associated with the trial function is stationary, so that the convergence rate becomes much faster. We thus achieve an example of Importance Sampling and in the present work we apply it to the Feynman-Kac(FK) path integrals for the ground and first few excited state energies for He to speed up the convergence rate. We calculate the path integrals using space averaging rather than the time averaging as done in the past. The best previous calculations from Variational computations report precisions of Hartrees, whereas in most cases our path integral results obtained for the ground and first excited states of He are lower than these results by about Hartrees or more.


title: Equivalence of Palm measures for determinantal point processes governed by Bergman kernels

abstract: For a determinantal point process induced by the reproducing kernel of the weighted Bergman space $A^2(U, \omega)$ over a domain $U \subset \mathbb{C}^d$, we establish the mutual absolute continuity of reduced Palm measures of any order provided that the domain $U$ contains a non-constant bounded holomorphic function. The result holds in all dimensions. The argument uses the $H^\infty(U)$-module structure of $A^2(U, \omega)$. A corollary is the quasi-invariance of our determinantal point process under the natural action of the group of compactly supported diffeomorphisms of $U$.


title: On the minimum output entropy of random orthogonal quantum channels

abstract: We consider sequences of random quantum channels defined using the Stinespring formula with Haar-distributed random orthogonal matrices. For any fixed sequence of input states, we study the asymptotic eigenvalue distribution of the outputs through tensor powers of random channels. We show that the input states achieving minimum output entropy are tensor products of maximally entangled states (Bell states) when the tensor power is even. This phenomenon is completely different from the one for random quantum channels constructed from Haar-distributed random unitary matrices, which leads us to formulate some conjectures about the regularized minimum output entropy.


title: Freeness of multi-reflection arrangements via primitive vector fields

abstract: In 2002, Terao showed that every reflection multi-arrangement of a real reflection group with constant multiplicity is free by providing a basis of the module of derivations. We first generalize Terao's result to multi-arrangements stemming from well-generated unitary reflection groups, where the multiplicity of a hyperplane depends on the order of its stabilizer. Here the exponents depend on the exponents of the dual reflection representation. We then extend our results further to all imprimitive irreducible unitary reflection groups. In this case the exponents turn out to depend on the exponents of a certain Galois twist of the dual reflection representation that comes from a Beynon-Lusztig type semi-palindromicity of the fake degrees.


title: MURS: Mitigating Memory Pressure in Service-oriented Data Processing Systems

abstract: Although a data processing system often works as a batch processing system, many enterprises deploy such a system as a service, which we call the service-oriented data processing system. It has been shown that in-memory data processing systems suffer from serious memory pressure. The situation becomes even worse for the service-oriented data processing systems due to various reasons. For example, in a service-oriented system, multiple submitted tasks are launched at the same time and executed in the same context in the resources, comparing with the batch processing mode where the tasks are processed one by one. Therefore, the memory pressure will affect all submitted tasks, including the tasks that only incur the light memory pressure when they are run alone. In this paper, we find that the reason why memory pressure arises is because the running tasks produce massive long-living data objects in the limited memory space. Our studies further reveal that the long-living data objects are generated by the API functions that are invoked by the in-memory processing frameworks. Based on these findings, we propose a method to classify the API functions based on the memory usage rate. Further, we design a scheduler called MURS to mitigate the memory pressure. We implement MURS in Spark and conduct the experiments to evaluate the performance of MURS. The results show that when comparing to Spark, MURS can 1) decrease the execution time of the submitted jobs by up to 65.8\%, 2) mitigate the memory pressure in the server by decreasing the garbage collection time by up to 81\%, and 3) reduce the data spilling, and hence disk I/O, by approximately 90\%.


title: Querying Log Data with Metric Temporal Logic (Technical Report)

abstract: We propose a novel framework for ontology-based access to temporal log data using a datalog extension datalogMTL of a Horn fragment of the metric temporal logic MTL. We show that datalogMTL is ExpSpace-complete even with punctual intervals, in which case full MTL is known to be undecidable. We also prove that nonrecursive datalogMTL is PSpace-complete for combined complexity and in AC0 for data complexity. We demonstrate by two real-world use cases that nonrecursive datalogMTL programs can express complex temporal concepts from typical user queries and thereby facilitate access to temporal log data. Our experiments with Siemens turbine data and MesoWest weather data show that datalogMTL ontology-mediated queries are efficient and scale on large datasets.


title: Equation of State Based Slip Spring Model for Entangled Polymer Dynamics

abstract: A mesoscopic, mixed particle- and field-based Brownian dynamics methodology for the simulation of entangled polymer melts has been developed. Polymeric beads consist of several Kuhn segments, and their motion is dictated by the Helmholtz energy of the sample, which is a sum of the entropic elasticity of chain strands between beads, slip springs, and nonbonded interactions. The entanglement effect is introduced by the slip springs, which are springs connecting either nonsuccessive beads on the same chain or beads on different polymer chains. The terminal positions of slip springs are altered during the simulation through a kinetic Monte Carlo hopping scheme, with rate-controlled creation/destruction processes for the slip springs at chain ends. The rate constants are consistent with the free energy function employed and satisfy microscopic reversibility at equilibrium. The free energy of nonbonded interactions is derived from an appropriate equation of state, and it is computed as a functional of the local density by passing an orthogonal grid through the simulation box; accounting for it is necessary for reproducing the correct compressibility of the polymeric material. Parameters invoked by the mesoscopic model are derived from experimental volumetric and viscosity data or from atomistic molecular dynamics simulations, establishing a "bottom-up" predictive framework for conducting slip spring simulations of polymeric systems of specific chemistry. The mesoscopic simulation methodology is implemented for the case of cis-1,4-polyisoprene, whose structure, dynamics, thermodynamics, and linear rheology in the melt state are quantitatively predicted and validated without a posteriori fitting the results to experimental measurements.


title: Thermal Expansion of the Heavy-fermion Superconductor PuCoGa$_{5}$

abstract: We have performed high-resolution powder x-ray diffraction measurements on a sample of $^{242}$PuCoGa$_{5}$, the heavy-fermion superconductor with the highest critical temperature $T_{c}$ = 18.7 K. The results show that the tetragonal symmetry of its crystallographic lattice is preserved down to 2 K. Marginal evidence is obtained for an anomalous behaviour below $T_{c}$ of the $a$ and $c$ lattice parameters. The observed thermal expansion is isotropic down to 150 K, and becomes anisotropic for lower temperatures. This gives a $c/a$ ratio that decreases with increasing temperature to become almost constant above $\sim$150 K. The volume thermal expansion coefficient $\alpha_{V}$ has a jump at $T_{c}$, a factor $\sim$20 larger than the change predicted by the Ehrenfest relation for a second order phase transition. The volume expansion deviates from the curve expected for the conventional anharmonic behaviour described by a simple Gr\"{u}neisen-Einstein model. The observed differences are about ten times larger than the statistical error bars but are too small to be taken as an indication for the proximity of the system to a valence instability that is avoided by the superconducting state.


title: TCP in 5G mmWave Networks: Link Level Retransmissions and MP-TCP

abstract: MmWave communications, one of the cornerstones of future 5G mobile networks, are characterized at the same time by a potential multi-gigabit capacity and by a very dynamic channel, sensitive to blockage, wide fluctuations in the received signal quality, and possibly also sudden link disruption. While the performance of physical and MAC layer schemes that address these issues has been thoroughly investigated in the literature, the complex interactions between mmWave links and transport layer protocols such as TCP are still relatively unexplored. This paper uses the ns-3 mmWave module, with its channel model based on real measurements in New York City, to analyze the performance of the Linux TCP/IP stack (i) with and without link-layer retransmissions, showing that they are fundamental to reach a high TCP throughput on mmWave links and (ii) with Multipath TCP (MP-TCP) over multiple LTE and mmWave links, illustrating which are the throughput-optimal combinations of secondary paths and congestion control algorithms in different conditions.


title: Optical Line Emission from z$\sim$6.8 Sources with Deep Constraints on Ly$\alpha$ Visibility

abstract: We analyze a sample of $z$-dropout galaxies in the CANDELS GOODS South and UDS fields that have been targeted by a dedicated spectroscopic campaign aimed at detecting their Ly$\alpha$ line. Deep IRAC observations at 3.6 and 4.5 $\mu$m are used to determine the strength of optical emission lines affecting these bands at z$\sim$6.5-6.9 in order to i) investigate possible physical differences between Ly$\alpha$ emitting and non-emitting sources; ii) constrain the escape fraction of ionizing photons; iii) provide an estimate of the specific star-formation rate at high redshifts. We find evidence of strong [OIII]+H$\beta$ emission in the average (stacked) SEDs of galaxies both with and without Ly$\alpha$ emission. The blue IRAC [3.6]-[4.5] color of the stack with detected Ly$\alpha$ line can be converted into a rest-frame equivalent width EW([OIII]+H$\beta$)=1500$^{+530}_{-440}\AA$ assuming a flat intrinsic stellar continuum. This strong optical line emission enables a first estimate of f$_{esc}\lesssim$20% on the escape fraction of ionizing photons from Ly$\alpha$ detected objects. The objects with no Ly$\alpha$ line show less extreme EW([OIII]+H$\beta$)=520$^{+170}_{-150}\AA$ suggesting different physical conditions of the HII regions with respect to Ly$\alpha$-emitting ones, or a larger f$_{esc}$. The latter case is consistent with a combined evolution of f$_{esc}$ and the neutral hydrogen fraction as an explanation of the lack of bright Ly$\alpha$ emission at z$>$6. A lower limit on the specific star formation rate, SSFR$>$9.1$Gyr^{-1}$ for $M_{star}=2 \times 10^9 M_{\odot}$ galaxies at these redshifts can be derived from the spectroscopically confirmed sample.


title: LIDAR-based Driving Path Generation Using Fully Convolutional Neural Networks

abstract: In this work, a novel learning-based approach has been developed to generate driving paths by integrating LIDAR point clouds, GPS-IMU information, and Google driving directions. The system is based on a fully convolutional neural network that jointly learns to carry out perception and path generation from real-world driving sequences and that is trained using automatically generated training examples. Several combinations of input data were tested in order to assess the performance gain provided by specific information modalities. The fully convolutional neural network trained using all the available sensors together with driving directions achieved the best MaxF score of 88.13% when considering a region of interest of 60x60 meters. By considering a smaller region of interest, the agreement between predicted paths and ground-truth increased to 92.60%. The positive results obtained in this work indicate that the proposed system may help fill the gap between low-level scene parsing and behavior-reflex approaches by generating outputs that are close to vehicle control and at the same time human-interpretable.


title: How would GW150914 look with future GW detector networks?

abstract: The first detected gravitational wave signal, GW150914, was produced by the coalescence of a stellar-mass binary black hole. Along with the subsequent detection of GW151226, GW170104 and the candidate event LVT151012, this gives us evidence for a population of black hole binaries with component masses in the tens of solar masses. As detector sensitivity improves, this type of source is expected to make a large contribution to the overall number of detections, but has received little attention compared to binary neutron star systems in studies of projected network performance. We simulate the observation of a system like GW150914 with different proposed network configurations, and study the precision of parameter estimates, particularly source location, orientation and masses. We find that the improvements to low frequency sensitivity that are expected with continued commissioning will improve the precision of chirp mass estimates by an order of magnitude, whereas the improvements in sky location and orientation are driven by the expanded network configuration. This demonstrates that both sensitivity and number of detectors will be important factors in the scientific potential of second generation detector networks.


title: Boundedness of Pseudodifferential Operators with symbols in Wiener amalgam spaces on Modulation Spaces

abstract: This paper provides sufficient conditions for the boundedness of Weyl operators on modulation spaces. The Weyl symbols belong to Wiener amalgam spaces, or generalized modulation spaces, as recently renamed by their inventor Hans Feichtinger. This is the first result which relates symbols in Wiener amalgam spaces to operators acting on classical modulation spaces.


title: The role of substrate bias and nitrogen doping on the structural evolution and local elastic modulus of diamond-like carbon films

abstract: Diamond-like carbon (DLC) films are synthesized on Si using plasma enhanced chemical vapor deposition. The role of substrate bias and nitrogen doping on the structural evolution and local elastic modulus of DLC films are systematically investigated. Raman spectroscopic studies reveal that the amount of graphitic C=C sp2 bonding increases with substrate bias and nitrogen doping. The density and hydrogen concentration in the films are found to vary from 0.7 to 2.2 g/cm3 and 16 to 38 atomic %, respectively, depending upon the substrate bias and nitrogen concentration in the DLC films. Atomic force acoustic microscopic (AFAM) analysis shows a direct correlation between local elastic modulus and structural properties estimated by Raman spectroscopy, Rutherford back scattering and elastic recoil detection analysis. AFAM analysis further confirms the evolution of soft second phases at high substrate biases (> -150V) in undoped DLC films. Further, N doping leads to formation of such soft second phases in DLC films even at lower substrate bias of -100 V. The AFAM studies provide a direct microscopic evidence for the "sub-implantation growth model" which predicts the formation of graphitic second phases in DLC matrix at high substrate biases.


title: Multilabel Classification with R Package mlr

abstract: We implemented several multilabel classification algorithms in the machine learning package mlr. The implemented methods are binary relevance, classifier chains, nested stacking, dependent binary relevance and stacking, which can be used with any base learner that is accessible in mlr. Moreover, there is access to the multilabel classification versions of randomForestSRC and rFerns. All these methods can be easily compared by different implemented multilabel performance measures and resampling methods in the standardized mlr framework. In a benchmark experiment with several multilabel datasets, the performance of the different methods is evaluated.


title: A solar blowout jet caused by the eruption of a magnetic flux rope

abstract: We investigate the three-dimensional (3D) magnetic structure of a blowout jet originated in the west edge of NOAA Active Region (AR) 11513 on 02 July 2012 by means of recently developed forced field extrapolation (FFE) model. The results show that the blowout jet was caused by the eruption of the magnetic flux rope (MFR) consisting of twisted field lines. We further calculate the twist number $\mathcal{T}_{w}$ and squashing factor Q of the reconstructed magnetic field and find that (1) the MFR corresponds well to the high $\mathcal{T}_{w}$ region (2) the MFR outer boundary corresponds well to the high Q region, probably interpreting the bright structure at the base of the jet. The twist number of the MFR is estimated to be $\mathcal{T}_{w}=-1.54\pm 0.67$. Thus, the kink instability is regarded as the initiation mechanism of the blowout jet as $\mathcal{T}_{w}$ reaching or even exceeding the threshold value of the kink instability. Our results also indicate that the bright point at the decaying phase is actually comprised of some small loops that are heated by the reconnection occurred above. In summary, the blowout jet is mostly consistent with the scenario proposed by \citet{mcs10} except that the kink instability is found to be a possible trigger.


title: No inflation in type IIA strings on rigid CY spaces

abstract: We investigate whether cosmological inflation is possible in a class of flux compactifications of type IIA strings on rigid Calabi-Yau manifolds, when all perturbative string corrections are taken into account. We confine ourselves to the universal hypermultiplet and an abelian vector multiplet, representing matter in four dimensions. Since all axions can be stabilized by D-instantons, we choose dilaton and a K\"ahler modulus as the only running scalars. Though positivity of their scalar potential can be achieved, we find that there is no slow roll ($\varepsilon > 2$ and $\eta >2$), and no Graceful Exit because the scalar potential has the run-away behaviour resulting in decompactification. We conclude that it is impossible to generate phenomenologically viable inflation in the given class of flux compactifications of type IIA strings, without explicit breaking of $N=2$ local supersymmetry of the low-energy effective action.


title: Value of Information: Sensitivity Analysis and Research Design in Bayesian Evidence Synthesis

abstract: Suppose we have a Bayesian model which combines evidence from several different sources. We want to know which model parameters most affect the estimate or decision from the model, or which of the parameter uncertainties drive the decision uncertainty. Furthermore we want to prioritise what further data should be collected. These questions can be addressed by Value of Information (VoI) analysis, in which we estimate expected reductions in loss from learning specific parameters or collecting data of a given design. We describe the theory and practice of VoI for Bayesian evidence synthesis, using and extending ideas from health economics, computer modelling and Bayesian design. The methods are general to a range of decision problems including point estimation and choices between discrete actions. We apply them to a model for estimating prevalence of HIV infection, combining indirect information from several surveys, registers and expert beliefs. This analysis shows which parameters contribute most of the uncertainty about each prevalence estimate, and provides the expected improvements in precision from collecting specific amounts of additional data.


title: On the Limit Imbalanced Logistic Regression by Binary Predictors

abstract: In this work, we introduce a modified (rescaled) likelihood for imbalanced logistic regression. This new approach makes easier the use of exponential priors and the computation of lasso regularization path. Precisely, we study a limiting behavior for which class imbalance is artificially increased by replication of the majority class observations. If some strong overlap conditions are satisfied, the maximum likelihood estimate converges towards a finite value close to the initial one (intercept excluded) as shown by simulations with binary predictors. This solution corresponds to the extremum of a concave function that we refer to as "rescaled" likelihood. In this context, the use of exponential priors has a clear interpretation as a shift on the predictor means for the minority class. Thanks to the simple binary structure, some random designs give analytic path estimators for the lasso regularization problem. An effective approximate path algorithm by piecewise logarithmic functions based on matrix inversions is also presented. This work was motivated by its potential application to spontaneous reports databases in a pharmacovigilance context.


title: Confirming the oblique rotator model for the extremely slowly rotating O8f?p star HD 108

abstract: The O8f?p star HD 108 is implied to have experienced the most extreme rotational braking of any magnetic, massive star, with a rotational period $P_{\rm rot}$ of at least 55 years, but the upper limit on its spindown timescale is over twice the age estimated from the Hertzsprung-Russell diagram. HD 108's observed X-ray luminosity is also much higher than predicted by the XADM model, a unique discrepancy amongst magnetic O-type stars. Previously reported magnetic data cover only a small fraction ($\sim$3.5\%) of $P_{\rm rot}$, and were furthermore acquired when the star was in a photometric and spectroscopic `low state' at which the longitudinal magnetic field $\langle B_z\rangle$~was likely at a minimum. We have obtained a new ESPaDOnS magnetic measurement of HD 108, 6 years after the last reported measurement. The star is returning to a spectroscopic high state, although its emission lines are still below their maximum observed strength, consistent with the proposed 55-year period. We measured $\langle B_z\rangle=-325 \pm 45$ G, twice the strength of the 2007-2009 observations, raising the lower limit of the dipole surface magnetic field strength to $B_{\rm d} \ge 1$ kG. The simultaneous increase in $\langle B_z\rangle$~and emission strength is consistent with the oblique rotator model. Extrapolation of the $\langle B_z\rangle$~maximum via comparison of HD 108's spectroscopic and magnetic data with the similar Of?p star HD 191612 suggests that $B_{\rm d} > 2$~kG, yielding $t_{\rm S, max}<3$~Myr, compatible with the stellar age. These results also yield a better agreement between the observed X-ray luminosity and that predicted by the XADM model.


title: Probing the dusty stellar populations of the Local Volume Galaxies with JWST/MIRI

abstract: The Mid-Infrared Instrument (MIRI) for the {\em James Webb Space Telescope} (JWST) will revolutionize our understanding of infrared stellar populations in the Local Volume. Using the rich {\em Spitzer}-IRS spectroscopic data-set and spectral classifications from the Surveying the Agents of Galaxy Evolution (SAGE)-Spectroscopic survey of over a thousand objects in the Magellanic Clouds, the Grid of Red supergiant and Asymptotic giant branch star ModelS ({\sc grams}), and the grid of YSO models by Robitaille et al. (2006), we calculate the expected flux-densities and colors in the MIRI broadband filters for prominent infrared stellar populations. We use these fluxes to explore the {\em JWST}/MIRI colours and magnitudes for composite stellar population studies of Local Volume galaxies. MIRI colour classification schemes are presented; these diagrams provide a powerful means of identifying young stellar objects, evolved stars and extragalactic background galaxies in Local Volume galaxies with a high degree of confidence. Finally, we examine which filter combinations are best for selecting populations of sources based on their JWST colours.


title: Arithmetic progressions in middle Nth cantor sets

abstract: We show the middle Nth cantor set contains arithmetic progressions of length at least proportional to N/log_2(N).


title: Coupling of Crop Assignment and Vehicle Routing for Harvest Planning in Agriculture

abstract: A method for harvest planning based on the coupling of crop assignment with vehicle routing is presented. Given a setting with multiple fields, a path network connecting these, multiple depots at which a number of harvesters are initially located, the main question addressed is: Which crop out of a set of different crops to assign to each field when accounting for the given setting? It must be answered by every farm manager at the beginning of every work-cycle starting with plant seeding and ending with harvesting. Rather than solving an assignment problem only, it is here also accounted for the connectivity between fields. In practice, fields are often located distant apart. Traveling costs of machinery and limited harvesting windows demand optimised operation and route planning. Therefore, the proposed method outputs crop assignments to fields and simultaneously determines crop-tours, i.e., optimised sequences in which to service fields of the same crop during harvest. It is of particular relevance for larger farms and groups of farms that collaborate and share machinery. Integer programming based exact algorithms are derived. For large numbers of fields, where these algorithms may not be tractable due to computational constraints, elements of clustering and the solution of local Traveling Salesman Problems are added, thereby on the one hand rendering the method heuristic and in general suboptimal, but on the other hand maintaining large-scale applicability.


title: Water transport to circumprimary habitable zones from icy planetesimal disks in binary star systems

abstract: So far, more than 130 extrasolar planets have been found in multiple stellar systems. Dynamical simulations show that the outcome of the planetary formation process can lead to different planetary architectures (i.e. location, size, mass, and water content) when the star system is single or double. In the late phase of planetary formation, when embryo-sized objects dominate the inner region of the system, asteroids are also present and can provide additional material for objects inside the habitable zone (HZ). In this study, we make a comparison of several binary star systems and aim to show how efficient they are at moving icy asteroids from beyond the snow line into orbits crossing the HZ. We also analyze the influence of secular and mean motion resonances on the water transport towards the HZ. Our study shows that small bodies also participate in bearing a non-negligible amount of water to the HZ. The proximity of a companion moving on an eccentric orbit increases the flux of asteroids to the HZ, which could result in a more efficient water transport on a short timescale, causing a heavy bombardment. In contrast to asteroids moving under the gravitational perturbations of one G-type star and a gas giant, we show that the presence of a companion star not only favors a faster depletion of our disk of planetesimals, but can also bring 4-5 times more water into the whole HZ. However, due to the secular resonance located either inside the HZ or inside the asteroid belt, impacts between icy planetesimals from the disk and big objects in the HZ can occur at high impact speed. Therefore, real collision modeling using a GPU 3D-SPH code show that in reality, the water content of the projectile is greatly reduced and therefore, also the water transported to planets or embryos initially inside the HZ.


title: $p$-adic dynamical systems of $(2,2)$-rational functions with unique fixed point

abstract: We consider a family of $(2,2)$-rational functions given on the set of complex $p$-adic field $\mathbb{C}_p$. Each such function has a unique fixed point. We study $p$-adic dynamical systems generated by the $(2,2)$-rational functions. We show that the fixed point is indifferent and therefore the convergence of the trajectories is not the typical case for the dynamical systems. Siegel disks of these dynamical systems are found. We obtain an upper bound for the set of limit points of each trajectory, i.e., we determine a sufficiently small set containing the set of limit points. For each $(2,2)$-rational function on $\mathbb{C}_p$ there are two points $\hat x_1=\hat x_1(f)$, $\hat x_2=\hat x_2(f)\in \mathbb{C}_p$ which are zeros of its denominator. We give explicit formulas of radiuses of spheres (with the center at the fixed point) containing some points such that the trajectories (under actions of $f$) of the points after a finite step come to $\hat x_1$ or $\hat x_2$. Moreover for a class of $(2,2)$-rational functions we study ergodicity properties of the dynamical systems on the set of $p$-adic numbers $Q_p$. For each such function we describe all possible invariant spheres. We show that the $p$-adic dynamical system reduced on each invariant sphere is not ergodic with respect to Haar measure.


title: On certain orbits of geodesic flow and (a,b)-continued fractions

abstract: In this article, we characterize two kinds of exceptional orbits of the geodesic flow associated with the Modular surface in terms of a two-parameter family of continued fraction expansion of endpoints of the lifts to the hyperbolic plane of the corresponding geodesics. As a consequence, we obtain an extension of Dani correspondence between homogeneous dynamics and Diophantine approximation.


title: Rational ergodicity of Step function Skew Products

abstract: We study rational step function skew products over certain rotations of the circle proving ergodicity and bounded rational ergodicity when rotation number is a quadratic irrational. The latter arises from a consideration of the asymptotic temporal statistics of an orbit as modelled by an associated affine random walk.


title: Gain and time resolution of 45 $\mu$m thin Low Gain Avalanche Detectors before and after irradiation up to a fluence of $10^{15}$ n$_{eq}$/cm$^2$

abstract: Low Gain Avalanche Detectors (LGADs) are silicon sensors with a built-in charge multiplication layer providing a gain of typically 10 to 50. Due to the combination of high signal-to-noise ratio and short rise time, thin LGADs provide good time resolutions. LGADs with an active thickness of about 45 $\mu$m were produced at CNM Barcelona. Their gains and time resolutions were studied in beam tests for two different multiplication layer implantation doses, as well as before and after irradiation with neutrons up to $10^{15}$ n$_{eq}$/cm$^2$. The gain showed the expected decrease at a fixed voltage for a lower initial implantation dose, as well as for a higher fluence due to effective acceptor removal in the multiplication layer. Time resolutions below 30 ps were obtained at the highest applied voltages for both implantation doses before irradiation. Also after an intermediate fluence of $3\times10^{14}$ n$_{eq}$/cm$^2$, similar values were measured since a higher applicable reverse bias voltage could recover most of the pre-irradiation gain. At $10^{15}$ n$_{eq}$/cm$^2$, the time resolution at the maximum applicable voltage of 620 V during the beam test was measured to be 57 ps since the voltage stability was not good enough to compensate for the gain layer loss. The time resolutions were found to follow approximately a universal function of gain for all implantation doses and fluences.