title: Extraction of neutrino mixing parameters from experiments with multiple identical detectors

abstract: The statistical method used in the analyses of measurements of neutrino oscillation mixing angle $\theta_{13}$ by the Daya Bay Collaboration is based on variational minimization of a $\chi^2$ function defined in terms of quantities of interest and pull factors which are introduced to deal with effects of systematic uncertainties. For both experiments, the number of parameters that need to be determined is great than the number of available data points (20 vs 6 for the Daya Bay and 12 vs 2 for the RENO). While the results for the mixing angle and the normalization factor were reported, results for the other parameters (pull factors) were omitted in their publications. There exist multiple sets of parameters from the minimization of the $\chi^2$ function. We investigate the sensitivity of the extracted mixing angle on this non-uniqueness of minimization results for the Daya Bay data using two methods of minimization. We report results for all parameters, including those of physics interest and pull factors. The obtained results for the mixing angle and the normalization factor are in agreement with those reported by the Daya Bay Collaboration. Furthermore, we present plots of confidence level contours in the space of the mixing angle and normalization factor. We also present results from fittings using a reduced $\chi^2$ function with fewer parameters than the one employed by the Daya Bay Collaboration.


title: Speech Dereverberation Based on Integrated Deep and Ensemble Learning Algorithm

abstract: Reverberation, which is generally caused by sound reflections from walls, ceilings, and floors, can result in severe performance degradation of acoustic applications. Due to a complicated combination of attenuation and time-delay effects, the reverberation property is difficult to characterize, and it remains a challenging task to effectively retrieve the anechoic speech signals from reverberation ones. In the present study, we proposed a novel integrated deep and ensemble learning algorithm (IDEA) for speech dereverberation. The IDEA consists of offline and online phases. In the offline phase, we train multiple dereverberation models, each aiming to precisely dereverb speech signals in a particular acoustic environment; then a unified fusion function is estimated that aims to integrate the information of multiple dereverberation models. In the online phase, an input utterance is first processed by each of the dereverberation models. The outputs of all models are integrated accordingly to generate the final anechoic signal. We evaluated the IDEA on designed acoustic environments, including both matched and mismatched conditions of the training and testing data. Experimental results confirm that the proposed IDEA outperforms single deep-neural-network-based dereverberation model with the same model architecture and training data.


title: Noisy Expectation-Maximization: Applications and Generalizations

abstract: We present a noise-injected version of the Expectation-Maximization (EM) algorithm: the Noisy Expectation Maximization (NEM) algorithm. The NEM algorithm uses noise to speed up the convergence of the EM algorithm. The NEM theorem shows that injected noise speeds up the average convergence of the EM algorithm to a local maximum of the likelihood surface if a positivity condition holds. The generalized form of the noisy expectation-maximization (NEM) algorithm allow for arbitrary modes of noise injection including adding and multiplying noise to the data. We demonstrate these noise benefits on EM algorithms for the Gaussian mixture model (GMM) with both additive and multiplicative NEM noise injection. A separate theorem (not presented here) shows that the noise benefit for independent identically distributed additive noise decreases with sample size in mixture models. This theorem implies that the noise benefit is most pronounced if the data is sparse. Injecting blind noise only slowed convergence.


title: A novel concept of pseudo ternary diffusion couple for the estimation of diffusion coefficients in multicomponent systems

abstract: A pseudo ternary diffusion couple technique in a multicomponent system by simplifying the mathematical complications of Onsager formalism is proposed for the estimation of composition dependent values of the interdiffusion coefficients. This is otherwise impossible following the conventional method in a system with more than three components. Other alternative methods estimate the average diffusion coefficients over a composition range of random choice and lack physical significance. This method can be followed in a multicomponent system with any number of components on the condition that only three components develop diffusion profiles keeping others as constant.


title: A3T: Adversarially Augmented Adversarial Training

abstract: Recent research showed that deep neural networks are highly sensitive to so-called adversarial perturbations, which are tiny perturbations of the input data purposely designed to fool a machine learning classifier. Most classification models, including deep learning models, are highly vulnerable to adversarial attacks. In this work, we investigate a procedure to improve adversarial robustness of deep neural networks through enforcing representation invariance. The idea is to train the classifier jointly with a discriminator attached to one of its hidden layer and trained to filter the adversarial noise. We perform preliminary experiments to test the viability of the approach and to compare it to other standard adversarial training methods.


title: Early Universe with modified scalar-tensor theory of gravity

abstract: Scalar-tensor theory of gravity with non-minimal coupling is a fairly good candidate for dark energy, required to explain late-time cosmic evolution. Here we study the very early stage of evolution of the universe with a modified version of the theory, which includes scalar curvature squared term. One of the key aspects of the present study is that, the quantum dynamics of the action under consideration ends up generically with de-Sitter expansion under semiclassical approximation, rather than power-law. This justifies the analysis of inflationary regime with de-Sitter expansion. The other key aspect is that, while studying gravitational perturbation, the perturbed generalized scalar field equation obtained from the perturbed action, when matched with the perturbed form of the background scalar field equation, relates the coupling parameter and the potential exactly in the same manner as the solution of classical field equations does, assuming de-Sitter expansion. The study also reveals that the quantum theory is well behaved, inflationary parameters fall well within the observational limit and quantum perturbation analysis shows that the power-spectrum does not deviate considerably from the standard one obtained from minimally coupled theory.


title: Emergent memory in cell signaling: Persistent adaptive dynamics in cascades can arise from the diversity of relaxation time-scales

abstract: The mitogen-activated protein kinase (MAPK) signaling cascade, an evolutionarily conserved motif present in all eukaryotic cells, is involved in coordinating critical cell-fate decisions, regulating protein synthesis, and mediating learning and memory. While the steady-state behavior of the pathway stimulated by a time-invariant signal is relatively well-understood, we show using a computational model that it exhibits a rich repertoire of transient adaptive responses to changes in stimuli. When the signal is switched on, the response is characterized by long-lived modulations in frequency as well as amplitude. On withdrawing the stimulus, the activity decays over timescales much longer than that of phosphorylation-dephosphorylation processes, exhibiting reverberations characterized by repeated spiking in the activated MAPK concentration. The long-term persistence of such post-stimulus activity suggests that the cascade retains memory of the signal for a significant duration following its removal, even in the absence of any explicit feedback or cross-talk with other pathways. We find that the molecular mechanism underlying this behavior is related to the existence of distinct relaxation rates for the different cascade components. This results in the imbalance of fluxes between different layers of the cascade, with the repeated reuse of activated kinases as enzymes when they are released from sequestration in complexes leading to one or more spike events following the removal of the stimulus. The persistent adaptive response reported here, indicative of a cellular "short-term" memory, suggests that this ubiquitous signaling pathway plays an even more central role in information processing by eukaryotic cells.


title: Multiple Imputation: A Review of Practical and Theoretical Findings

abstract: Multiple imputation is a straightforward method for handling missing data in a principled fashion. This paper presents an overview of multiple imputation, including important theoretical results and their practical implications for generating and using multiple imputations. A review of strategies for generating imputations follows, including recent developments in flexible joint modeling and sequential regression/chained equations/fully conditional specification approaches. Finally, we compare and contrast different methods for generating imputations on a range of criteria before identifying promising avenues for future research.


title: HelMod in the works: from direct observations to the local interstellar spectrum of cosmic-ray electrons

abstract: The local interstellar spectrum (LIS) of cosmic-ray (CR) electrons for the energy range 1 MeV to 1 TeV is derived using the most recent experimental results combined with the state-of-the-art models for CR propagation in the Galaxy and in the heliosphere. Two propagation packages, GALPROP and HelMod, are combined to provide a single framework that is run to reproduce direct measurements of CR species at different modulation levels, and at both polarities of the solar magnetic field. An iterative maximum-likelihood method is developed that uses GALPROP-predicted LIS as input to HelMod, which provides the modulated spectra for specific time periods of the selected experiments for model-data comparison. The optimized HelMod parameters are then used to adjust GALPROP parameters to predict a refined LIS with the procedure repeated subject to a convergence criterion. The parameter optimization uses an extensive data set of proton spectra from 1997-2015. The proposed CR electron LIS accommodates both the low-energy interstellar spectra measured by Voyager 1 as well as the high-energy observations by PAMELA and AMS-02 that are made deep in the heliosphere; it also accounts for Ulysses counting rate features measured out of the ecliptic plane. The interstellar and heliospheric propagation parameters derived in this study agree well with our earlier results for CR protons, helium nuclei, and anti-protons propagation and LIS obtained in the same framework.


title: Leakage suppression by ultrafast pulse shaping

abstract: We consider the leakage suppression problem of a three-level system in which the first two levels are the qubit system and the third, weakly-coupled to the second, is the leakage state. We show that phase- and amplitude-controlled two (three) pulses are sufficient for arbitrary qubit controls from the ground (an arbitrary) initial state, with leakage suppressed up to the first order of perturbation without additional pulse-area cost. A proof-of-principle experiment was performed with shaped femtosecond optical pulses and atomic rubidium showing a good agreement with the theory.


title: Gamma-ray Bursts Induced by Turbulent Reconnection

abstract: We consider a simple model for gamma-ray bursts induced by magnetic reconnection in turbulent media. The magnetic field in a jet is subject to kink instabilities, which distort the regular structure of the spiral magnetic field, drive turbulence, and trigger reconnection. The resulting reconnection takes place in a high Reynolds number medium, where turbulence is further enhanced and in turn accelerates the reconnection process. This boot-strap reconnection gives rise to bursts of reconnection events, through which the free energy of magnetic field is transformed into a gamma-ray burst. The efficiency of magnetic reconnection and magnetic energy dissipation is not constrained by microphysical properties of plasmas.


title: MINE: Mutual Information Neural Estimation

abstract: We argue that the estimation of mutual information between high dimensional continuous random variables can be achieved by gradient descent over neural networks. We present a Mutual Information Neural Estimator (MINE) that is linearly scalable in dimensionality as well as in sample size, trainable through back-prop, and strongly consistent. We present a handful of applications on which MINE can be used to minimize or maximize mutual information. We apply MINE to improve adversarially trained generative models. We also use MINE to implement Information Bottleneck, applying it to supervised classification; our results demonstrate substantial improvement in flexibility and performance in these settings.


title: How Many Samples Required in Big Data Collection: A Differential Message Importance Measure

abstract: Information collection is a fundamental problem in big data, where the size of sampling sets plays a very important role. This work considers the information collection process by taking message importance into account. Similar to differential entropy, we define differential message importance measure (DMIM) as a measure of message importance for continuous random variable. It is proved that the change of DMIM can describe the gap between the distribution of a set of sample values and a theoretical distribution. In fact, the deviation of DMIM is equivalent to Kolmogorov-Smirnov statistic, but it offers a new way to characterize the distribution goodness-of-fit. Numerical results show some basic properties of DMIM and the accuracy of the proposed approximate values. Furthermore, it is also obtained that the empirical distribution approaches the real distribution with decreasing of the DMIM deviation, which contributes to the selection of suitable sampling points in actual system.


title: State Variation Mining: On Information Divergence with Message Importance in Big Data

abstract: Information transfer which reveals the state variation of variables usually plays a vital role in big data analytics and processing. In fact, the measures for information transfer could reflect the system change by use of the variable distributions, similar to KL divergence and Renyi divergence. Furthermore, in terms of the information transfer in big data, small probability events usually dominate the importance of the total message to some degree. Therefore, it is significant to design an information transfer measure based on the message importance which emphasizes the small probability events. In this paper, we propose a message importance transfer measure (MITM) and investigate its characteristics and applications on three aspects. First, the message importance transfer capacity based on MITM is presented to offer an upper bound for the information transfer process with disturbance. Then, we extend the MITM to the continuous case and discuss the robustness by using it to measuring information distance. Finally, we utilize the MITM to guide the queue length selection in the caching operation of mobile edge computing.


title: Deep Stereo Matching with Explicit Cost Aggregation Sub-Architecture

abstract: Deep neural networks have shown excellent performance for stereo matching. Many efforts focus on the feature extraction and similarity measurement of the matching cost computation step while less attention is paid on cost aggregation which is crucial for stereo matching. In this paper, we present a learning-based cost aggregation method for stereo matching by a novel sub-architecture in the end-to-end trainable pipeline. We reformulate the cost aggregation as a learning process of the generation and selection of cost aggregation proposals which indicate the possible cost aggregation results. The cost aggregation sub-architecture is realized by a two-stream network: one for the generation of cost aggregation proposals, the other for the selection of the proposals. The criterion for the selection is determined by the low-level structure information obtained from a light convolutional network. The two-stream network offers a global view guidance for the cost aggregation to rectify the mismatching value stemming from the limited view of the matching cost computation. The comprehensive experiments on challenge datasets such as KITTI and Scene Flow show that our method outperforms the state-of-the-art methods.


title: Symbolic Timed Observational Equivalence

abstract: Intruders can infer properties of a system by measuring the time it takes for the system to respond to some request of a given protocol, that is, by exploiting time side channels. These properties may help intruders distinguish whether a system is a honeypot or concrete system helping him avoid defense mechanisms, or track a user among others violating his privacy. Observational equivalence is the technical machinery used for verifying whether two systems are distinguishable. Moreover, efficient symbolic methods have been developed for automating the check of observational equivalence of systems. This paper introduces a novel definition of timed observational equivalence which also distinguishes systems according to their time side channels. Moreover, as our definition uses symbolic time constraints, it can be automated by using SMT-solvers.


title: Content Based Status Updates

abstract: Consider a stream of status updates generated by a source, where each update is of one of two types: high priority or ordinary (low priority). These updates are to be transmitted through a network to a monitor. However, the transmission policy of each packet depends on the type of stream it belongs to. For the low priority stream, we analyze and compare the performances of two transmission schemes: (i) Ordinary updates are served in a First-Come-First-Served (FCFS) fashion, whereas, in (ii), the ordinary updates are transmitted according to an M/G/1/1 with preemption policy. In both schemes, high priority updates are transmitted according to an M/G/1/1 with preemption policy and receive preferential treatment. An arriving priority update discards and replaces any currently-in-service high priority update, and preempts (with eventual resume for scheme (i)) any ordinary update. We model the arrival processes of the two kinds of updates, in both schemes, as independent Poisson processes. For scheme (i), we find the arrival and service rates under which the system is stable and give closed-form expressions for average peak age and a lower bound on the average age of the ordinary stream. For scheme (ii), we derive closed-form expressions for the average age and average peak age of the high priority and low priority streams. We finally show that, if the service time is exponentially distributed, the M/M/1/1 with preemption policy leads to an average age of the low priority stream higher than the one achieved using the FCFS scheme. Therefore, the M/M//1/1 with preemption policy, when applied on the low priority stream of updates and in the presence of a higher priority scheme, is not anymore the optimal transmission policy from an age point of view.


title: Status Updates in a multi-stream M/G/1/1 preemptive queue

abstract: We consider a source that collects a multiplicity of streams of updates and sends them through a network to a monitor. However, only a single update can be in the system at a time. Therefore, the transmitter always preempts the packet being served when a new update is generated. We consider Poisson arrivals for each stream and a common general service time, and refer to this system as the multi-stream M/G/1/1 queue with preemption. Using the detour flow graph method, we compute a closed form expression for the average age and the average peak age of each stream. Moreover, we deduce that although all streams are treated equally from a transmission point of view (they all preempt each other), one can still prioritize a stream from an age point of view by simply increasing its generation rate. However, this will increase the sum of the ages which is minimized when all streams have the same update rate.


title: Predicting Smartphone Battery Life based on Comprehensive and Real-time Usage Data

abstract: Smartphones and smartphone apps have undergone an explosive growth in the past decade. However, smartphone battery technology hasn't been able to keep pace with the rapid growth of the capacity and the functionality of smartphones and apps. As a result, battery has always been a bottleneck of a user's daily experience of smartphones. An accurate estimation of the remaining battery life could tremendously help the user to schedule their activities and use their smartphones more efficiently. Existing studies on battery life prediction have been primitive due to the lack of real-world smartphone usage data at scale. This paper presents a novel method that uses the state-of-the-art machine learning models for battery life prediction, based on comprehensive and real-time usage traces collected from smartphones. The proposed method is the first that identifies and addresses the severe data missing problem in this context, using a principled statistical metric called the concordance index. The method is evaluated using a dataset collected from 51 users for 21 months, which covers comprehensive and fine-grained smartphone usage traces including system status, sensor indicators, system events, and app status. We find that the remaining battery life of a smartphone can be accurately predicted based on how the user uses the device at the real-time, in the current session, and in history. The machine learning models successfully identify predictive features for battery life and their applicable scenarios.


title: A Fast Algorithm with Error Bounds for Quadrature by Expansion

abstract: Quadrature by Expansion (QBX) is a quadrature method for approximating the value of the singular integrals encountered in the evaluation of layer potentials. It exploits the smoothness of the layer potential by forming locally-valid expansion which are then evaluated to compute the near or on-surface value of the integral. Recent work towards coupling of a Fast Multipole Method (FMM) to QBX yielded a first step towards the rapid evaluation of such integrals (and the solution of related integral equations), albeit with only empirically understood error behavior. In this paper, we improve upon this approach with a modified algorithm for which we give a comprehensive analysis of error and cost in the case of the Laplace equation in two dimensions. For the same levels of (user-specified) accuracy, the new algorithm empirically has cost-per-accuracy comparable to prior approaches. We provide experimental results to demonstrate scalability and numerical accuracy.


title: Spectroscopic Observations of the Outflowing Wind in the Lensed Quasar SDSS J1001+5027

abstract: We performed spectroscopic observations of the small-separation lensed quasar SDSS J1001+5027, whose images have an angular separation $\theta \sim 2.^{\!\!\prime\prime}86$, and placed constraints on the physical properties of gas clouds in the vicinity of the quasar (i.e., in the outflowing wind launched from the accretion disk). The two cylinders of sight to the two lensed images go through the same region of the outflowing wind and they become fully separated with no overlap at a very large distance from the source ($\sim 330$ pc). We discovered a clear difference in the profile of the CIV broad absorption line (BAL) detected in the two lensed images in two observing epochs. Because the kinematic components in the BAL profile do not vary in concert, the observed variations cannot be reproduced by a simple change of ionization state. If the variability is due to gas motion around the background source (i.e., the continuum source), the corresponding rotational velocity is $v_{rot}\geq 18,000$ km/s, and their distance from the source is $r\leq 0.06$ pc assuming Keplerian motion. Among three MgII and three CIV NAL systems that we detected in the spectra, only the MgII system at $z_{abs} = 0.8716$ shows a hint of variability in its MgI profile on a rest-frame time scale of $\Delta t_{rest}$ $\leq 191$ days and an obvious velocity shear between the sightlines whose physical separation is $\sim 7$ kpc. We interpret this as the result of motion of a cosmologically intervening absorber, perhaps located in a foreground galaxy.


title: Controller Synthesis for Safety of Physically-Viable Data-Driven Models

abstract: We consider the problem of designing finite-horizon safe controllers for a dynamical system for which no explicit analytical model exists and limited data only along a single trajectory of the system are available. Given samples of the states and inputs of the system, and additional side information in terms of regularity of the evolution of the states and conservation laws, we synthesize a controller such that the evolution of the states avoid some pre-specified unsafe set over a given finite horizon. Motivated by recent results on Whitney's extension theorem, we use piecewise-polynomial approximations of the trajectories based on the data along with the regularity side information to formulate a data-driven differential inclusion model that can predict the evolution of the trajectories. For these classes of data- driven differential inclusions, we propose a safety analysis theorem based on barrier certificates. As a corollary of this theorem, we demonstrate that we can design controllers ensuring safety of the solutions to the data-driven differential inclusion over a finite horizon. From a computational standpoint, our results are cast into a set of sum-of-squares programs whenever the certificates are parametrized by polynomials of fixed degree and the sets are semi-algebraic. This computational method allows incorporating side information in terms of conservation laws and integral constraints using Positivstellensatz-type arguments.


title: Quantum Many-Body Effects in X-Ray Spectra Efficiently Computed using a Basic Graph Algorithm

abstract: The growing interest in using x-ray spectroscopy for refined materials characterization calls for accurate electronic-structure theory to interpret x-ray near-edge fine structure. In this work, we propose an efficient and unified framework to describe all the many-electron processes in a Fermi liquid after a sudden perturbation (such as a core hole). This problem has been visited by the Mahan-Nozi\'eres-De Dominicis (MND) theory, but it is intractable to implement various Feynman diagrams within first-principles calculations. Here, we adopt a non-diagrammatic approach and treat all the many-electron processes in the MND theory on an equal footing. Starting from a recently introduced determinant formalism [Phys. Rev. Lett. 118, 096402 (2017)], we exploit the linear-dependence of determinants describing different final states involved in the spectral calculations. An elementary graph algorithm, breadth-first search, can be used to quickly identify the important determinants for shaping the spectrum, which avoids the need to evaluate a great number of vanishingly small terms. This search algorithm is performed over the tree-structure of the many-body expansion, which mimics a path-finding process. We demonstrate that the determinantal approach is computationally inexpensive even for obtaining x-ray spectra of extended systems. Using Kohn-Sham orbitals from two self-consistent fields (ground and core-excited state) as input for constructing the determinants, the calculated x-ray spectra for a number of transition metal oxides are in good agreement with experiments. Many-electron aspects beyond the Bethe-Salpeter equation, as captured by this approach, are also discussed, such as shakeup excitations and many-body wave function overlap considered in Anderson's orthogonality catastrophe.


title: Singlet fermionic dark matter with Veltman conditions

abstract: We reexamine a renormalizable model of a fermionic dark matter with a gauge singlet Dirac fermion and a real singlet scalar which can ameliorate the scalar mass hierarchy problem of the Standard Model (SM). Our model setup is the minimal extension of the SM for which a realistic dark matter (DM) candidate is provided and the cancellation of one-loop quadratic divergence to the scalar masses can be achieved by the Veltman condition (VC) simultaneously. This model extension, although renormalizable, can be considered as an effective low-energy theory valid up to cut-off energies about 10 TeV. We calculate the one-loop quadratic divergence contributions of the new scalar and fermionic DM singlets, and constrain the model parameters using the VC and the perturbative unitarity conditions. Taking into account the invisible Higgs decay measurement, we show the allowed region of new physics parameters satisfying the recent measurement of relic abundance. With the obtained parameter set, we predict the elastic scattering cross section of the new singlet fermion into target nuclei for a direct detection of the dark matter. We also perform the full analysis with arbitrary set of parameters without the VC as a comparison, and discuss the implication of the constraints by the VC in detail.


title: Laplace, Residue, and Euler integral representations of GKZ hypergeometric functions

abstract: We consider four types of representations of solutions of GKZ system: series representations, Laplace integral representations, Euler integral representations, and Residue integral representations which will be introduced in this paper. In the former half of this paper, we provide a method for constructing integration cycles for Laplace, Residue, and Euler integral representations and relate them to series representations. In the latter half, we reformulate our integral representations in terms of direct images of $D$-modules and show their equivalence.


title: How to augment a small learning set for improving the performances of a CNN-based steganalyzer?

abstract: Deep learning and convolutional neural networks (CNN) have been intensively used in many image processing topics during last years. As far as steganalysis is concerned, the use of CNN allows reaching the state-of-the-art results. The performances of such networks often rely on the size of their learning database. An obvious preliminary assumption could be considering that "the bigger a database is, the better the results are". However, it appears that cautions have to be taken when increasing the database size if one desire to improve the classification accuracy i.e. enhance the steganalysis efficiency. To our knowledge, no study has been performed on the enrichment impact of a learning database on the steganalysis performance. What kind of images can be added to the initial learning set? What are the sensitive criteria: the camera models used for acquiring the images, the treatments applied to the images, the cameras proportions in the database, etc? This article continues the work carried out in a previous paper, and explores the ways to improve the performances of CNN. It aims at studying the effects of "base augmentation" on the performance of steganalysis using a CNN. We present the results of this study using various experimental protocols and various databases to define the good practices in base augmentation for steganalysis.


title: Optimal control of an evolution equation with non-smooth dissipation

abstract: In the present work we study the optimal control of an evolution equation with non-smooth dissipation. The solution mapping of this system is non-smooth and hence the analysis is quite challenging. Our approach is to regularize the dissipation via approximation by a smooth function. We derive optimality conditions for the corresponding smooth optimal control problem. Then we drive the regularization parameter to zero and obtain necessary optimality conditions for the original non-smooth problem. However, in this process we lose regularity of the adjoint variables.


title: First-order symmetry-adapted perturbation theory for multiplet splittings

abstract: We present a symmetry-adapted perturbation theory (SAPT) for the interaction of two high-spin open-shell molecules (described by their restricted open-shell Hartree-Fock determinants) resulting in low-spin states of the complex. The previously available SAPT formalisms, except for some system-specific studies for few-electron complexes, were restricted to the high-spin state of the interacting system. Thus, the new approach provides, for the first time, a SAPT-based estimate of the splittings between different spin states of the complex. We have derived and implemented the lowest-order SAPT term responsible for these splittings, that is, the first-order exchange energy. We show that within the so-called S2 approximation commonly used in SAPT (neglecting effects that vanish as fourth or higher powers of intermolecular overlap integrals), the first-order exchange energies for all multiplets are linear combinations of two matrix elements: a diagonal exchange term that determines the spin-averaged effect and a {\em spin-flip term} responsible for the splittings between the states. The numerical factors in this linear combination are determined solely by the Clebsch-Gordan coefficients: accordingly, the S2 approximation implies a Heisenberg Hamiltonian picture with a single coupling strength parameter determining all the splittings. The new approach is cast into both molecular-orbital and atomic-orbital expressions: the latter enable an efficient density-fitted implementation. We test the newly developed formalism on several open-shell complexes ranging from diatomic systems (Li+H, Mn+Mn,...) to the phenalenyl dimer.


title: On asymptotic stability of ground states of some systems of nonlinear Schr\"odinger equations

abstract: We extend to a specific class of systems of nonlinear Schr\"odinger equations (NLS) the theory of asymptotic stability of ground states already proved for the scalar NLS. Here the key point is the choice of an adequate system of modulation coordinates and the novelty, compared to the scalar NLS, is the fact that the group of symmetries of the system is non-commutative.


title: Optimal contracts under competition when uncertainty from adverse selection and moral hazard are present

abstract: In a continuous-time setting where a risk-averse agent controls the drift of an output process driven by a Brownian motion, optimal contracts are linear in the terminal output; this result is well-known in a setting with moral hazard and -under stronger assumptions - adverse selection. We show that this result continues to hold when in addition reservation utilities are type-dependent. This type of problem occurs in the study of optimal compensation problems involving competing principals.


title: Separation of Instrument Sounds using Non-negative Matrix Factorization with Spectral Envelope Constraints

abstract: Spectral envelope is one of the most important features that characterize the timbre of an instrument sound. However, it is difficult to use spectral information in the framework of conventional spectrogram decomposition methods. We overcome this problem by suggesting a simple way to provide a constraint on the spectral envelope calculated by linear prediction. In the first part of this study, we use a pre-trained spectral envelope of known instruments as the constraint. Then we apply the same idea to a blind scenario in which the instruments are unknown. The experimental results reveal that the proposed method outperforms the conventional methods.


title: 21cm Signal Recovery via the Robust Principle Component Analysis

abstract: The redshifted 21~cm signal from neutral hydrogen (HI) is potentially a very powerful probe for cosmology, but a difficulty in its observation is that it is much weaker than foreground radiation from the Milky Way as well as extragalactic radio sources. The foreground radiation at different frequencies are however coherent along one line of sight, and various methods of foreground subtraction based on this property have been proposed. In this paper, we present a new method based on the Robust Principal Component Analysis (RPCA) to subtract foreground and extract 21~cm signal, which explicitly uses both the low-rank property of the frequency covariance matrix (i.e. frequency coherence) of the foreground and the sparsity of the frequency covariance matrix of the 21~cm signal. The low-rank property of the foregrounds frequency covariance has been exploited in many previous works on foreground subtraction, but to our knowledge the sparsity of the frequency covariance of the 21~cm signal is first explored here. By exploiting both properties in the RPCA method, in principle, the foreground and signal may be separated without the signal loss problem. Our method is applicable to both small patch of sky with the flat-sky approximation, and to large area of sky where the sphericity has to be considered. It is also easy to be extended to deal with more complex conditions such as sky map with defects.


title: UV self-completion of a theory of Superfluid Dark Matter

abstract: We show that the model of superfluid dark matter developed in Refs.~\cite{Khoury:2014tka,Berezhiani:2015bqa,Berezhiani:2015pia}, which modifies the Newtonian potential and explains the galactic rotational curves, can be unitarized by the formation of classical configurations in the scattering amplitudes. The classicalization mechanism may also trigger the formation of the superfluid state from the early to the late Universe.


title: Speculose: Analyzing the Security Implications of Speculative Execution in CPUs

abstract: Whenever modern CPUs encounter a conditional branch for which the condition cannot be evaluated yet, they predict the likely branch target and speculatively execute code. Such pipelining is key to optimizing runtime performance and is incorporated in CPUs for more than 15 years. In this paper, to the best of our knowledge, we are the first to study the inner workings and the security implications of such speculative execution. We revisit the assumption that speculatively executed code leaves no traces in case it is not committed. We reveal several measurable side effects that allow adversaries to enumerate mapped memory pages and to read arbitrary memory---all using only speculated code that was never fully executed. To demonstrate the practicality of such attacks, we show how a user-space adversary can probe for kernel pages to reliably break kernel-level ASLR in Linux in under three seconds and reduce the Windows 10 KASLR entropy by 18~bits in less than a second.


title: How should a fixed budget of dwell time be spent in scanning electron microscopy to optimize image quality?

abstract: In scanning electron microscopy, the achievable image quality is often limited by a maximum feasible acquisition time per dataset. Particularly with regard to three-dimensional or large field-of-view imaging, a compromise must be found between a high amount of shot noise, which leads to a low signal-to-noise ratio, and excessive acquisition times. Assuming a fixed acquisition time per frame, we compared three different strategies for algorithm-assisted image acquisition in scanning electron microscopy. We evaluated (1) raster scanning with a reduced dwell time per pixel followed by a state-of-the-art Denoising algorithm, (2) raster scanning with a decreased resolution in conjunction with a state-of-the-art Super Resolution algorithm, and (3) a sparse scanning approach where a fixed percentage of pixels is visited by the beam in combination with state-of-the-art inpainting algorithms. Additionally, we considered increased beam currents for each of the strategies. The experiments showed that sparse scanning using an appropriate reconstruction technique was superior to the other strategies.


title: Maximal and typical nonnegative ranks of nonnegative tensors

abstract: Let $N_1, \ldots, N_d$ be positive integers with $N_1\leq\cdots\leq N_d$. Set $N=N_1\cdots N_{d-1}$. We show in this paper that an integer $r$ is a typical nonnegative rank of nonnegative tensors of format $N_1\times\cdots\times N_d$ if and only if $r\leq N$ and $r$ is greater than or equals to the generic rank of tensors over $\mathbb{C}$ of format $N_1\times\cdots\times N_d$. We also show that the maximal nonnegative rank of nonnegative tensors of format $N_1\times\cdots\times N_d$ is $N$.


title: Gene regulatory network inference: an introductory survey

abstract: Gene regulatory networks are powerful abstractions of biological systems. Since the advent of high-throughput measurement technologies in biology in the late 90s, reconstructing the structure of such networks has been a central computational problem in systems biology. While the problem is certainly not solved in its entirety, considerable progress has been made in the last two decades, with mature tools now available. This chapter aims to provide an introduction to the basic concepts underpinning network inference tools, attempting a categorisation which highlights commonalities and relative strengths. While the chapter is meant to be self-contained, the material presented should provide a useful background to the later, more specialised chapters of this book.


title: Non self-adjoint laplacians on a directed graph

abstract: We consider a non self-adjoint Laplacian on a directed graph with non symmetric edge weights. We analyse spectral properties of this Laplacian under a Kirchhoff assumption. Moreover we establish isoperimet-ric inequalities in terms of the numerical range to show the absence of the essential spectrum of the Laplacian on heavy end directed graphs.


title: Inflationary dynamics and preheating of the non-minimally coupled inflaton field in the metric and Palatini formalisms

abstract: The inflationary dynamics and preheating in a model with a non-minimally coupled inflaton field in the metric and Palatini formalisms are studied in this paper. We find that in both formalisms, irrespective of the initial conditions, our Universe will evolve into a slow-roll inflationary era and then the scalar field rolls into an oscillating phase. The value of the scalar field at the end of the inflation in the Palatini formalism is always larger than that in the metric one, which becomes more and more obvious with the increase of the absolute value of the coupling parameter $|\xi|$. During the preheating, we find that the inflaton quanta are produced explosively due to the parameter resonance and the growth of inflaton quanta will be terminated by the backreaction. With the increase of $|\xi|$, the resonance bands gradually close to the zero momentum ($k=0$), and the structure of resonance changes and becomes broader and broader in the metric formalism, while, it remains to be narrow in the Palatini formalism. The energy transfer from the inflaton field to the fluctuation becomes more and more efficient with the increase of $|\xi|$, and in the metric formalism the growth of the efficiency of energy transfer is much faster than that in the Palatini formalism. Therefore, the inflation and preheating show different characteristics in different formalisms.


title: On notions of Q-independence and Q-identical distributiveness

abstract: In a recent article A.M. Kagan and G.J.Sz\'ekely introduced a notion of Q-independent and Q-identical distributed random variables. We give a complete description of polynomials which appear in these definitions.


title: Multivariate stochastic delay differential equations and CAR representations of CARMA processes

abstract: In this study we show how to represent a continuous time autoregressive moving average (CARMA) as a higher order stochastic delay differential equation, which may be thought of as a continuous-time equivalent of the AR($\infty$) representation. Furthermore, we show how this representation gives rise to a prediction formula for CARMA processes. To be used in the above mentioned results we develop a general theory for multivariate stochastic delay differential equations, which will be of independent interest, and which will have particular focus on existence, uniqueness and representations.


title: Summary and report Working Group 6: theory and simulation, for the third edition of the European Advanced Accelerator Concept Conference

abstract: The theory and simulation working group (working group 6) of the third edition of the European Advanced Accelerator Workshop has been characterized by a strong numerical connotation. Particle in cell codes have proven to be a necessary tool to finely investigate the underlying physics both for laser and plasma wakefield acceleration processes. This year, the section has been characterised by an interest in the limitation of the numerical Cherenkov effect, the mitigation of the hose-instability, start-to-end simulations and more generally on new numerical schemes and diagnostics.


title: Laser wakefield acceleration with high-power, few-cycle mid-IR lasers

abstract: The study of laser wakefield electron acceleration (LWFA) using mid-IR laser drivers is a promising path for future laser driven electronaccelerators, when compared to traditional near-IR laser drivers uperating at 0.8-1 {\mu}m central wavelength ({\lambda}laser), as the necessary vector potential a_0 for electron injection can be achieved with smaller laser powers due to the linear dependence on {\lambda}laser. In this work, we perform 2D PIC simulations on LWFA using few-cycle high power (5-15 TW) laser systems with {\lambda}laser ranging from 0.88-10 {\mu}m. Such few-cycle systems are currently under development, aiming at Gas High Harmonics Generation applications where the favourable {\lambda}laser^2 scaling extends the range of XUV photon energies. We keep a_0 and n_e/n_cr (n_e being the plasma density and n_cr being the critical density for each {\lambda}laser) as common denominators in outr simulations, allowing for comparisons between drivers of different {\lambda}laser, with respect to the accelerated electron beam energy, charge, and conversion efficiency. While the electron energies are mainly dominated by the plasma dynamics, the laser to electron beam energy conversion efficiency shows significant enhancement with longer wavelength laser drivers.


title: Infinite Families of Equivariantly Formal Toric Orbifolds

abstract: The simplicial wedge construction on simplicial complexes and simple polytopes has been used by a variety of authors to study toric and related spaces, including non-singular toric varieties, toric manifolds, intersections of quadrics and more generally, polyhedral products. In this paper we extend the analysis to include toric orbifolds. Our main results yield infinite families of toric orbifolds, derived from a given one, whose integral cohomology is free of torsion and is concentrated in even degrees, a property which might be termed \emph{integrally equivariantly formal}. In all cases, it is possible to give a description of the cohomology ring and to relate it to the cohomology of the original orbifold.


title: Sensitivity indices for independent groups of variables

abstract: In this paper, we study sensitivity indices for independent groups of variables and we look at the particular case of block-additive models. We show in this case that most of the Sobol indices are equal to zero and that Shapley effects can be estimated more efficiently. We then apply this study to Gaussian linear models, and we provide an efficient algorithm to compute the theoretical sensitivity indices. In numerical experiments, we show that this algorithm compares favourably to other existing methods. We also use the theoretical results to improve the estimation of the Shapley effects for general models, when the inputs form independent groups of variables.


title: Hierarchical Motion Consistency Constraint for Efficient Geometrical Verification in UAV Image Matching

abstract: This paper proposes a strategy for efficient geometrical verification in unmanned aerial vehicle (UAV) image matching. First, considering the complex transformation model between correspondence set in the image-space, feature points of initial candidate matches are projected onto an elevation plane in the object-space, with assistant of UAV flight control data and camera mounting angles. Spatial relationships are simplified as a 2D-translation in which a motion establishes the relation of two correspondence points. Second, a hierarchical motion consistency constraint, termed HMCC, is designed to eliminate outliers from initial candidate matches, which includes three major steps, namely the global direction consistency constraint, the local direction-change consistency constraint and the global length consistency constraint. To cope with scenarios with high outlier ratios, the HMCC is achieved by using a voting scheme. Finally, an efficient geometrical verification strategy is proposed by using the HMCC as a pre-processing step to increase inlier ratios before the consequent application of the basic RANSAC algorithm. The performance of the proposed strategy is verified through comprehensive comparison and analysis by using real UAV datasets captured with different photogrammetric systems. Experimental results demonstrate that the generated motions have noticeable separation ability, and the HMCC-RANSAC algorithm can efficiently eliminate outliers based on the motion consistency constraint, with a speedup ratio reaching to 6 for oblique UAV images. Even though the completeness sacrifice of approximately 7 percent of points is observed from image orientation tests, competitive orientation accuracy is achieved from all used datasets. For geometrical verification of both nadir and oblique UAV images, the proposed method can be a more efficient solution.


title: Exceptional and modern intervals of the Tamari lattice

abstract: In this article we use the theory of interval-posets recently introduced by Ch{\^a}tel and Pons in order to describe some interesting families of intervals in the Tamari lattices. These families are defined as interval-posets avoiding specific configurations. At first, we consider what we call exceptional interval-posets and show that they correspond to the intervals which are obtained as images of noncrossing trees in the Dendriform operad. We also show that the exceptional intervals are exactly the intervals of the Tamari lattice induced by intervals in the poset of noncrossing partitions. In the second part we introduce the notion of modern and infinitely modern interval-posets. We show that the modern intervals are in bijection with the new intervals of the Tamari lattice in the sense of Chapoton. We deduce an intrinsic characterization of the new intervals in the Tamari lattice. Finally, we consider the family of what we call infinitely modern intervals and we we prove that there are as many infinitely modern interval-posets of size n as there are ternary trees with n inner vertices.


title: Combinatorics of compactified universal Jacobians

abstract: We use orientations on stable graphs to express the combinatorial structure of the compactified universal Jacobians in degrees g-1 and g over the moduli space of stable curves, \Mgb, and construct for them graded stratifications compatible with the one of \Mgb. In particular, for a stable curve we exhibit graded stratifications of the compactified Jacobians in terms of totally cyclic, respectively rooted, orientations on subgraphs of its dual graph.


title: Trust-Aware Decision Making for Human-Robot Collaboration: Model Learning and Planning

abstract: Trust in autonomy is essential for effective human-robot collaboration and user adoption of autonomous systems such as robot assistants. This paper introduces a computational model which integrates trust into robot decision-making. Specifically, we learn from data a partially observable Markov decision process (POMDP) with human trust as a latent variable. The trust-POMDP model provides a principled approach for the robot to (i) infer the trust of a human teammate through interaction, (ii) reason about the effect of its own actions on human trust, and (iii) choose actions that maximize team performance over the long term. We validated the model through human subject experiments on a table-clearing task in simulation (201 participants) and with a real robot (20 participants). In our studies, the robot builds human trust by manipulating low-risk objects first. Interestingly, the robot sometimes fails intentionally in order to modulate human trust and achieve the best team performance. These results show that the trust-POMDP calibrates trust to improve human-robot team performance over the long term. Further, they highlight that maximizing trust alone does not always lead to the best performance.


title: Aspects of (2,2) and (0,2) hybrid models

abstract: In this work we study the topological rings of two dimensional (2,2) and (0,2) hybrid models. In particular, we use localization to derive a formula for the correlators in both cases, focusing on the B- and B/2-twists. Although our methods apply to a vast range of hybrid CFTs, we focus on hybrid models suitable for compactifications of the heterotic string. In this case, our formula provides unnormalized Yukawa couplings of the spacetime superpotential. We apply our techniques to hybrid phases of linear models, and we find complete agreement with known results in other phases. We also obtain a prediction for a certain class of correlators involving twisted operators in (2,2) Landau-Ginzburg orbifolds. For (0,2) theories, our argument does not rely on the existence of a (2,2) locus. Finally, we derive vanishing conditions concerning worldsheet instanton corrections in (0,2) B/2-twisted hybrid models.


title: Spatio-Temporal Linkage over Location Enhanced Services

abstract: We are witnessing an enormous growth in the volume of data generated by various online services. An important portion of this data contains geographic references, since many of these services are \emph{location-enhanced} and thus produce spatio-temporal records of their usage. We postulate that the spatio-temporal usage records belonging to the same real-world entity can be matched across records from different location-enhanced services. Linking spatio-temporal records enables data analysts and service providers to obtain information that they cannot derive by analyzing only one set of usage records. In this paper, we develop a new \emph{linkage model} that can be used to match entities from two sets of spatio-temporal usage records belonging to two different location-enhanced services. This linkage model is based on the concept of $k$-$l$ \emph{diversity} --- that we developed to capture both spatial and temporal aspects of the linkage. To realize this linkage model in practice, we develop a scalable linking algorithm called \emph{ST-Link}, which makes use of effective spatial and temporal filtering mechanisms that significantly reduce the search space for matching users. Furthermore, \emph{ST-Link} utilizes sequential scan procedures to avoid random disk access and thus scales to large datasets. We evaluated our work with respect to accuracy and performance using several datasets. Experiments show that \emph{ST-Link} is effective in practice for performing spatio-temporal linkage and can scale to large datasets.


title: Generative Single Image Reflection Separation

abstract: Single image reflection separation is an ill-posed problem since two scenes, a transmitted scene and a reflected scene, need to be inferred from a single observation. To make the problem tractable, in this work we assume that categories of two scenes are known. It allows us to address the problem by generating both scenes that belong to the categories while their contents are constrained to match with the observed image. A novel network architecture is proposed to render realistic images of both scenes based on adversarial learning. The network can be trained in a weakly supervised manner, i.e., it learns to separate an observed image without corresponding ground truth images of transmission and reflection scenes which are difficult to collect in practice. Experimental results on real and synthetic datasets demonstrate that the proposed algorithm performs favorably against existing methods.


title: Self-Predicting Boolean Functions

abstract: A Boolean function $g$ is said to be an optimal predictor for another Boolean function $f$, if it minimizes the probability that $f(X^{n})\neq g(Y^{n})$ among all functions, where $X^{n}$ is uniform over the Hamming cube and $Y^{n}$ is obtained from $X^{n}$ by independently flipping each coordinate with probability $\delta$. This paper is about self-predicting functions, which are those that coincide with their optimal predictor.


title: Multiple Antennas Secure Transmission under Pilot Spoofing and Jamming Attack

abstract: Transmitter-side channel state information (CSI) of the legitimate destination plays a critical role in physical layer secure transmissions. However, channel training procedure is vulnerable to the pilot spoofing attack (PSA) or pilot jamming attack (PJA) by an active eavesdropper (Eve), which inevitably results in severe private information leakage. In this paper, we propose a random channel training (RCT) based secure downlink transmission framework for a time division duplex (TDD) multiple antennas base station (BS). In the proposed RCT scheme, multiple orthogonal pilot sequences (PSs) are simultaneously allocated to the legitimate user (LU), and the LU randomly selects one PS from the assigned PS set to transmit. Under either the PSA or PJA, we provide the detailed steps for the BS to identify the PS transmitted by the LU, and to simultaneously estimate channels of the LU and Eve. The probability that the BS makes an incorrect decision on the PS of the LU is analytically investigated. Finally, closed-form secure beamforming (SB) vectors are designed and optimized to enhance the secrecy rates during the downlink transmissions. Numerical results show that the secrecy performance is greatly improved compared to the conventional channel training scheme wherein only one PS is assigned to the LU.


title: Active repositioning of storage units in Robotic Mobile Fulfillment Systems

abstract: In our work we focus on Robotic Mobile Fulfillment Systems in e-commerce distribution centers. These systems were designed to increase pick rates by employing mobile robots bringing movable storage units (so-called pods) to pick and replenishment stations as needed, and back to the storage area afterwards. One advantage of this approach is that repositioning of inventory can be done continuously, even during pick and replenishment operations. This is primarily accomplished by bringing a pod to a storage location different than the one it was fetched from, a process we call passive pod repositioning. Additionally, this can be done by explicitly bringing a pod from one storage location to another, a process we call active pod repositioning. In this work we introduce first mechanisms for the latter technique and conduct a simulation-based experiment to give first insights of their effect.


title: Perfect codes in generalized Fibonacci cubes

abstract: The {\em Fibonacci cube} of dimension $n$, denoted as $\Gamma\_n$, is the subgraph of the $n$-cube $Q\_n$ induced by vertices with no consecutive 1's. In an article of 2016 Ashrafi and his co-authors proved the non-existence of perfect codes in $\Gamma\_n$ for $n\geq 4$. As an open problem the authors suggest to consider the existence of perfect codes in generalization of Fibonacci cubes. The most direct generalization is the family $\Gamma\_n(1^s)$ of subgraphs induced by strings without $1^s$ as a substring where $s\geq 2$ is a given integer. We prove the existence of a perfect code in $\Gamma\_n(1^s)$ for $n=2^p-1$ and $s \geq 3.2^{p-2}$ for any integer $p\geq 2$.


title: Benchmarking cross-project defect prediction approaches with costs metrics

abstract: Defect prediction can be a powerful tool to guide the use of quality assurance resources. In recent years, many researchers focused on the problem of Cross-Project Defect Prediction (CPDP), i.e., the creation of prediction models based on training data from other projects. However, only few of the published papers evaluate the cost efficiency of predictions, i.e., if they save costs if they are used to guide quality assurance efforts. Within this paper, we provide a benchmark of 26 CPDP approaches based on cost metrics. Our benchmark shows that trivially assuming everything as defective is on average better than CPDP under cost considerations. Moreover, we show that our ranking of approaches using cost metrics is uncorrelated to a ranking based on metrics that do not directly consider costs. These findings show that we must put more effort into evaluating the actual benefits of CPDP, as the current state of the art of CPDP can actually be beaten by a trivial approach in cost-oriented evaluations.


title: Gamma-ray imaging system for real-time measurements in nuclear waste characterisation

abstract: A compact, portable and large field-of-view gamma camera to identify, locate and quantify gamma-ray emitting radioisotopes in real-time has been developed. The device delivers spectroscopic and imaging capabilities, which allow one to use it in a variety of nuclear waste characterisation scenarios, such as radioactivity monitoring in nuclear power plants and more specifically for the decommissioning of nuclear facilities. The technical development of this apparatus and some examples of its application in field measurements are reported in this article. The performance of the presented gamma-camera is also benchmarked against other conventional techniques.


title: Couplings in $L^p$ distance of two Brownian motions and their L{\'e}vy area

abstract: We study co-adapted couplings of (canonical hypoelliptic) diffu-sions on the (subRiemannian) Heisenberg group, that we call (Heisenberg) Brow-nian motions and are the joint laws of a planar Brownian motion with its L{\'e}vy area. We show that contrary to the situation observed on Riemannian manifolds of non-negative Ricci curvature, for any co-adapted coupling, two Heisenberg Brownian motions starting at two given points can not stay at bounded distance for all time t $\ge$ 0. Actually, we prove the stronger result that they can not stay bounded in L p for p $\ge$ 2. We also study the coupling by reflection, and show that it stays bounded in L p for 0 $\le$ p < 1. Finally, we explain how the results generalise to the Heisenberg groups of higher dimension


title: Clinical and Non-clinical Effects on Surgery Duration: Statistical Modeling and Analysis

abstract: Surgery duration is usually used as an input to the operation room (OR) allocation and surgery scheduling problems. A good estimation of surgery duration benefits the operation planning in ORs. In contrast, we would like to investigate whether the allocation decisions in turn influence surgery duration. Using almost two years of data from a large hospital in China, we find evidence in support of our conjecture. Surgery duration decreases with the number of surgeries a surgeon performs in a day. Numerically, surgery duration will decrease by 10 minutes on average if a surgeon performs one more surgery. Furthermore, we find a non-linear relationship between surgery duration and the number of surgeries allocated to an OR. Also, a surgery's duration is affected by its position in a sequence of surgeries performed by one surgeon. In addition, surgeons exhibit different patterns on the effects of surgery type and position. Since the findings are obtained from a particular data set, We do not claim the generalizability. Instead, the analysis in this paper provides insights into surgery duration study in ORs.


title: On the goodness-of-fit of generalized linear geostatistical models

abstract: We propose a generalization of Zhang's coefficient of determination to generalized linear geostatistical models and illustrate its application to river-blindness mapping. The generalized coefficient of determination has a more intuitive interpretation than other measures of predictive performance and allows to assess the individual contribution of each explanatory variable and the random effects to spatial prediction. The developed methodology is also more widely applicable to any generalized linear mixed model.


title: Regression Based Expected Shortfall Backtesting

abstract: This paper introduces novel backtests for the risk measure Expected Shortfall (ES) following the testing idea of Mincer and Zarnowitz (1969). Estimating a regression framework for the ES stand-alone is infeasible, and thus, our tests are based on a joint regression for the Value at Risk and the ES, which allows for different test specifications. These ES backtests are the first which solely backtest the ES in the sense that they only require ES forecasts as input parameters. As the tests are potentially subject to model misspecification, we provide asymptotic theory under misspecification for the underlying joint regression. We find that employing a misspecification robust covariance estimator substantially improves the tests' performance. We compare our backtests to existing approaches and find that our tests outperform the competitors throughout all considered simulations. In an empirical illustration, we apply our backtests to ES forecasts for 200 stocks of the S&P 500 index.


title: Extended phase diagram of RNiC2 family: Linear scaling of the Peierls temperature

abstract: Physical properties for the late-lanthanide-based RNiC2 (R = Dy, Ho, Er, and Tm) ternary compounds are reported. All the compounds show antiferromagnetic ground state with the N\'eel temperature ranging from 3.4 K for HoNiC2 to 8.5 K for ErNiC2. The results of the transport and galvanomagnetic properties confirm a charge density wave state at and above room temperature with transition temperatures TCDW = 284, 335, 366, and 394 K for DyNiC2, HoNiC2, ErNiC2, and TmNiC2, respectively. The Peierls temperature TCDW scales linearly with the unit cell volume. A similar linear dependence has been observed for the temperature of the lock-in transition T1 as well. Beyond the intersection point of the trend lines, the lock-in transition is no longer observed. In this Rapid Communication we demonstrate an extended phase diagram for the RNiC2 family.


title: Tunable elastic Parity-Time symmetric structure based on the shunted piezoelectric materials

abstract: We theoretically and numerically report on tunable elastic Parity-Time (PT) symmetric structure based on shunted piezoelectric units. We show that the elastic loss and gain can be archived in piezoelectric materials when they are shunted by external circuits containing positive and negative resistances. We present and discuss, as an example, the strongly dependent relationship between the exceptional points of a three-layered system and the impedance of their external shunted circuit. The achieved results evidence the PT symmetric structures based on this proposed concept can actively be tuned without any change of their geometric configurations.


title: A Game Theoretic Approach to Hyperbolic Consensus Problems

abstract: We introduce the use of conservation laws to develop strategies in multi-player consensus games. First, basic well posedness results provide a reliable analytic setting. Then, a general non anticipative strategy is proposed through its rigorous analytic definitions and then tested by means of numerical integrations.


title: Advective accretion flow properties around rotating black holes - application to GRO J1655-40

abstract: We examine the properties of the viscous dissipative accretion flow around rotating black holes in presence of mass loss. Considering thin disc approximation, we self-consistently calculate the inflow-outflow solutions and observe that the mass outflow rates decreases with the increase of viscosity parameter ($\alpha$). Further, we carry out the model calculation of Quasi-periodic Oscillation frequency ($\nu_{\rm QPO}$) that is frequently observed in black hole sources and observe that $\nu^{\rm max}_{\rm QPO}$ increases with the increase of black hole spin ($a_k$). Then, we employ our model in order to explain the High Frequency Quasi-periodic Oscillations (HFQPOs) observed in black hole source GRO J1655-40. While doing this, we attempt to constrain the range of $a_k$ based on observed HFQPOs ($\sim$ 300 Hz and $\sim$ 450 Hz) for the black hole source GRO J1655-40.


title: Will This Video Go Viral? Explaining and Predicting the Popularity of Youtube Videos

abstract: What makes content go viral? Which videos become popular and why others don't? Such questions have elicited significant attention from both researchers and industry, particularly in the context of online media. A range of models have been recently proposed to explain and predict popularity; however, there is a short supply of practical tools, accessible for regular users, that leverage these theoretical results. HIPie -- an interactive visualization system -- is created to fill this gap, by enabling users to reason about the virality and the popularity of online videos. It retrieves the metadata and the past popularity series of Youtube videos, it employs Hawkes Intensity Process, a state-of-the-art online popularity model for explaining and predicting video popularity, and it presents videos comparatively in a series of interactive plots. This system will help both content consumers and content producers in a range of data-driven inquiries, such as to comparatively analyze videos and channels, to explain and predict future popularity, to identify viral videos, and to estimate response to online promotion.


title: Characterisation of beam driven ionisation injection in the blowout regime of Plasma Acceleration

abstract: Beam driven ionisation injection is characterised for a variety of high-Z dopant. We discuss the region of extraction and why the position where electrons are captured influences the final quality of the internally-injected bunch. The beam driven ionisation injection relies on the capability to produce a high gradient fields at the bubble closure, with magnitudes high enough to ionise by tunnelling effect the still bounded electrons (of a high-Z dopant). The ionised electrons are captured by the nonlinear plasma wave at the accelerating and focusing wake phase leading to high-brightness trailing bunches. The high transformer ratio guarantees that the ionisation only occurs at the bubble closure. The quality of the ionisation-injected trailing bunches strongly and non-linearly depends on the properties of the dopant gas (density and initial ionisation state). We use the full 3D PIC code ${\tt ALaDyn}$ to consider the highly three-dimensional nature of the effect. By means of a systematic approach we have investigated the emittance and energy spread formation and the evolution for different dopant gases and configurations.


title: Quantum phase transition in few-layer NbSe$_2$ probed through quantized conductance fluctuations

abstract: We present the first observation of dynamically modulated quantum phase transition (QPT) between two distinct charge density wave (CDW) phases in 2-dimensional 2H-NbSe$_2$. There is recent spectroscopic evidence for the presence of these two quantum phases, but its evidence in bulk measurements remained elusive. We studied suspended, ultra-thin \nbse devices fabricated on piezoelectric substrates - with tunable flakes thickness, disorder level and strain. We find a surprising evolution of the conductance fluctuation spectra across the CDW temperature: the conductance fluctuates between two precise values, separated by a quantum of conductance. These quantized fluctuations disappear for disordered and on-substrate devices. With the help of mean-field calculations, these observations can be explained as to arise from dynamical phase transition between the two CDW states. To affirm this idea, we vary the lateral strain across the device via piezoelectric medium and map out the phase diagram near the quantum critical point (QCP). The results resolve a long-standing mystery of the anomalously large spectroscopic gap in NbSe$_2$.


title: On melt solutions for the growth of CaTiO$_3$ crystals

abstract: When calcium titanate crystals are grown from stoichiometric melts, they crystallize in the cubic perovskite structure. Upon cooling to room temperature they undergo subsequent phase transitions to tetragonal and orthorhombic modifications. These phase transitions are disruptive and result in severely damaged crystals. This paper presents differential thermal analysis data for several prospective solvents, with the aim to identify a system offering the possibility to perform crystal growth of undistorted CaTiO$_3$ crystals by crystallizing them significantly below the melting point directly in the low temperature modification. From mixtures CaF$_2$:TiO$_2$:CaTiO$_3$ = 3:1:1 (molar ratio) the growth of undistorted, at least millimeter-sized CaTiO$_3$ crystals is possible.


title: Supercaloric functions for the porous medium equation

abstract: In the slow diffusion case unbounded supersolutions of the porous medium equation are of two totally different types, depending on whether the pressure is locally integrable or not. This criterion and its consequences are discussed.


title: Robust error estimation for lowest-order approximation of nearly incompressible elasticity

abstract: We consider so-called Herrmann and Hydrostatic mixed formulations of classical linear elasticity and analyse the error associated with locally stabilised $P_1-P_0$ finite element approximation. First, we prove a stability estimate for the discrete problem and establish an a priori estimate for the associated energy error. Second, we consider a residual-based a posteriori error estimator as well as a local Poisson problem estimator. We establish bounds for the energy error that are independent of the Lam\'{e} coefficients and prove that the estimators are robust in the incompressible limit. A key issue to be addressed is the requirement for pressure stabilisation. Numerical results are presented that validate the theory. The software used is available online.


title: A possible testbed for warped extra dimension from the angle of Buchdahl's limit

abstract: We consider a five dimensional AdS warped spacetime in presence of a massive scalar field in the bulk. The scalar field potential fulfills the requirement of modulus stabilization even when the effect of backreaction of the stabilizing field is taken into account. In such a scenario, we investigate the possible role of modulus field on a compact stellar structure from the perspective of four dimensional effective theory. Our result reveals that in the presence of the modulus field, the upper bound of mass-radius ratio (generally known as Buchdahl's limit) of a star can go beyond the general relativity prediction. Interestingly this provides a natural testbed for the existence of such higher dimensional modulus field.


title: Optical Follow-up of Planck Cluster Candidates with Small Instruments

abstract: We report on the search for optical counterparts of Planck Sunyaev-Zel'dovich (SZ) cluster candidates using a 0.6 meter non-professional telescope. Among the observed sources, an unconfirmed candidate, PSZ2 G156.24+22.32, is found to be associated with a region of more than 100 galaxies within a 3 arcminutes radius around the Sunyaev-Zel'dovich maximum signal coordinates. Using 14 hours of cumulated exposure over the Sloan color filters g', r', i', and, z', we estimate the photometric redshift of these galaxies at zphot=0.29 +- 0.08. Using the red-sequence galaxy method gives a photometric redshift of 0.30 +0.03 -0.05. Combined with the Planck SZ proxy mass function, this would favor a cluster of 4.4 x 10^{14} solar masses. This result suggests that a dedicated pool of observatories equipped with such instruments could collectively contribute to optical follow-up programs of massive cluster candidates at moderate redshifts.


title: Newton-Okounkov bodies and normal toric degenerations

abstract: Anderson proved that the finite generation of the value semigroup $\Gamma_{Y_\bullet}(D)$ in the construction of the Newton-Okounkov body $\Delta_{Y_\bullet}(D)$ induces a toric degeneration of the corresponding variety $X$ to some toric variety $X_0$. In this case the normalization of $X_0$ is the normal toric variety corresponding to the rational polytope $\Delta_{Y_\bullet}(D)$. Since $X_0$ is not normal in general this correspondence is rather implicit. In this article we investigate in conditions to assure that $X_0$ is normal, by comparing the Hilbert polynomial with the Ehrhart polynomial. In the case of del Pezzo surfaces this will result in an algorithm which outputs for a given divisor $D$ a flag $Y_\bullet$ such that the value semigroup in question is indeed normal. Furthermore, we will find flags on del Pezzo surfaces and on some particular weak del Pezzo surfaces which induce normal toric degenerations for all possible divisors at once. We will prove that in this case the global value semigroup $\Gamma_{Y_\bullet}(X)$ is finitely generated and normal.


title: Extending Whitney's extension theorem: nonlinear function spaces

abstract: We consider a global, nonlinear version of the Whitney extension problem for manifold-valued smooth functions on closed domains $C$, with non-smooth boundary, in possibly non-compact manifolds. Assuming $C$ is a submanifold with corners, or is compact and locally convex with rough boundary, we prove that the restriction map from everywhere-defined functions is a submersion of locally convex manifolds and so admits local linear splittings on charts. This is achieved by considering the corresponding restriction map for locally convex spaces of compactly-supported sections of vector bundles, allowing the even more general case where $C$ only has mild restrictions on inward and outward cusps, and proving the existence of an extension operator.


title: Optical left-handed metamaterials made of arrays of upright split-ring pairs

abstract: Electromagnetic metamaterials are man-made structures that have novel properties such as a negative refraction index, not attainable in naturally occurring materials. Although negative index materials (NIMs) in microwave frequencies were demonstrated in 2001, it has remained challenging to design NIMs for optical frequencies especially those with both negative permittivity and negative permeability [known as left-handed metamaterials (LHMs)]. Here, by going beyond the traditional concept of the combination of artificial electronic and magnetic meta-atoms to design NIMs, we propose a novel LHM composed of an array of simple upright split-ring pairs working in the near infrared region. Our electromagnetic simulations reveal the underlying mechanism that the coupling of the two rings can stimulate simultaneously both the electric and magnetic resonances. The proposed structure has a highest refractive index of -2, a highest figure of merit of 21, good air-matched impedance and 180 nm double negative bandwidth, which excel the performances of many previous proposals. We also numerically demonstrate the negative refraction of this metamaterial in both the single-layer form and wedge-shaped lens.


title: Improved bounds on the multicolor Ramsey numbers of paths and even cycles

abstract: We study the multicolor Ramsey numbers for paths and even cycles, $R_k(P_n)$ and $R_k(C_n)$, which are the smallest integers $N$ such that every coloring of the complete graph $K_N$ has a monochromatic copy of $P_n$ or $C_n$ respectively. For a long time, $R_k(P_n)$ has only been known to lie between $(k-1+o(1))n$ and $(k + o(1))n$. A recent breakthrough by S\'ark\"ozy and later improvement by Davies, Jenssen and Roberts give an upper bound of $(k - \frac{1}{4} + o(1))n$. We improve the upper bound to $(k - \frac{1}{2}+ o(1))n$. Our approach uses structural insights in connected graphs without a large matching. These insights may be of independent interest.


title: Gravity-capillary waves on the free surface of a liquid dielectric in a tangential electric field

abstract: Processes of propagation and interaction of nonlinear gravity-capillary waves on the free surface of a deep non-conducting liquid with high dielectric constant under the action of a tangential electric field are numerically simulated. The computational method is based on the time-dependent conformal transformation of the region occupied by the fluid into a half-plane. In the limit of a strong electric field, when the gravitational and capillary forces are negligibly small, there exists an exact analytical solution of the electro-hydrodynamic equations describing propagation without distortions of nonlinear surface waves along (or against) the electric field direction. In the situation where gravity and capillarity are taken into account, the results of numerical simulations indeed show that, for large external field, the waves traveling in a given direction tend to preserve their shape. In the limit of a strong electric field, the interaction of counter-propagating waves leads to the formation of regions, where the electrostatic and dynamic pressures undergo a discontinuity, and the curvature of the surface increases infinitely. The Fourier spectrum of the surface perturbations tends to the power-law distribution (k^-2). In the case of a finite electric field, the wave interaction results in a radiation of massive cascade of small-scale capillary waves that causes the chaotic behavior of the system. The investigated mechanism of interaction between oppositely-traveling waves can enhance development of the capillary turbulence of the fluid surface.


title: Interpretation of the vibrational spectra of glassy polymers using coarse-grained simulations

abstract: The structure and vibrational density of states (VDOS) of polymer glasses are investigated using numerical simulations based on the classical Kremer-Grest bead-spring model. We focus on the roles of chain length and bending stiffness, the latter being set by imposing three-body angular potentials along chain backbones. Upon increasing the chain length and bending stiffness, structural reorganisation leads to volumetric expansion of the material and build-up of internal stresses. The VDOS has two dominant bands: a low frequency one corresponding to inter- and intra-chain non-bonding interactions and a high frequency one corresponding principally to vibrations of bonded beads that constitute skeletal chain backbones. Upon increasing the steepness of the angular potential, vibrational modes associated with chain bending gradually move from the low-frequency to the high-frequency band. This redistribution of modes is reflected in a reduction of the so-called Boson peak upon increasing chain stiffness. Remarkably, the finer structure and the peaks of the high-frequency band, and their variations with stiffness, can, for short chains, be explained using an analytical solution derived for a model triatomic molecule. For longer chains, the qualitative evolution of the VDOS with chain stiffness is similar, although the distinct peaks observed for short chains become increasingly smoothed-out. Our findings can be used to guide a systematic approach to interpretation of Brillouin and Raman scattering spectra of glassy polymers in future work, with applications in polymer processing diagnostics.

categories: cond-mat.soft cond-mat.dis-nn cond-mat.mtrl-sci physics.chem-ph


title: On Partly Overloaded Spreading Sequences with Variable Spreading Factor

abstract: Future wireless communications systems are expected to support multi-service operation, i.e. especially multi-rate as well as multi-level quality of service (QoS) requirements. This evolution is mainly driven by the success of the Internet of Things (IoT) and the growing presence of machine type communication (MTC). Whereas in the last years information in wireless communication systems was mainly generated or at least requested by humans and was also processed by humans, we can now see a paradigm shift since so-called machine type communication is gaining growing importance. Along with these changes we also encounter changes regarding the quality of service requirements, data rate requirements, latency constraints, different duty cycles et cetera. The challenge for new communication systems will therefore be to enable different user types and their different requirements efficiently. In this paper, we present partly overloaded spreading sequences, i.e. sequences which are globally orthogonal and sequences which interfere with a subset of sequences while being orthogonal to the globally orthogonal sequences. Additionally, we are able to vary the spreading factor of these sequences, which allows us to flexibly assign appropriate sequences to different service types or user types respectively. We propose the use of these sequences for a CDMA channel access method which is able to flexibly support different traffic types.


title: Metrics for two electron random potential systems

abstract: Metrics have been used to investigate the relationship between wavefunction distances and density distances for families of specific systems. We extend this research to look at random potentials for time-dependent single electron systems, and for ground-state two electron systems. We find that Fourier series are a good basis for generating random potentials. These random potentials also yield quasi-linear relationships between the distances of ground-state densities and wavefunctions, providing a framework in which Density Functional Theory can be explored.


title: Eigenvalues of Dirichlet Laplacian within the class of open sets with constant diameter

abstract: This paper is about a shape optimization problem related to the Dirichlet Laplacian eingevalues in the Euclidean plane. More precisely we study the shape of the minimizer in the class of open sets of constant width. We prove that the disk is not a local minimizer except for a limited number of eigenvalues.


title: Deep Episodic Memory: Encoding, Recalling, and Predicting Episodic Experiences for Robot Action Execution

abstract: We present a novel deep neural network architecture for representing robot experiences in an episodic-like memory which facilitates encoding, recalling, and predicting action experiences. Our proposed unsupervised deep episodic memory model 1) encodes observed actions in a latent vector space and, based on this latent encoding, 2) infers most similar episodes previously experienced, 3) reconstructs original episodes, and 4) predicts future frames in an end-to-end fashion. Results show that conceptually similar actions are mapped into the same region of the latent vector space. Based on these results, we introduce an action matching and retrieval mechanism, benchmark its performance on two large-scale action datasets, 20BN-something-something and ActivityNet and evaluate its generalization capability in a real-world scenario on a humanoid robot.


title: Near Field Speckles The Optical Theorem Revisited

abstract: The zero angle scattered wave plays an essential role in the celebrated Optical Theorem that relates the total extinction cross section (including both scattering and absorption) to the complex amplitude of the scattered wave. In an authoritative paper where the possibilities of an experimental verification of Optical Theorem were discussed, it was pointed out that a proof could hardly come from an optics experiment. One must agree to this statement if one receives the scattering intensity from an assembly of particles as a Bragg reflection from a three dimensional random grating. Accordingly, scattered waves at any angle then behave as a circular Gaussian process, and phases are equally distributed in the interval (-{\pi}, +{\pi}), therefore giving no cue on the behaviour of the wave scattered at exactly zero angle. So a proof of the Optical Theorem would require of telling apart an incoming wave from its weak replica, and gauge the phase in between. In my PhD thesis, near field scattering method is exploited for the determination of this elusive phase of the scattered wave. Having determined the phase of the scattered wave, the experimental verification of Optical Theorem can also be attempted. In addition, the method can generate dynamic scattering data as well; it provides an elegant and efficient way to filter out stray multiple scattering artefacts and thus the method can be used even when multiple scattering is rampant.


title: First-passage times over moving boundaries for asymptotically stable walks

abstract: Let $\{S_n, n\geq1\}$ be a random walk wih independent and identically distributed increments and let $\{g_n,n\geq1\}$ be a sequence of real numbers. Let $T_g$ denote the first time when $S_n$ leaves $(g_n,\infty)$. Assume that the random walk is oscillating and asymptotically stable, that is, there exists a sequence $\{c_n,n\geq1\}$ such that $S_n/c_n$ converges to a stable law. In this paper we determine the tail behaviour of $T_g$ for all oscillating asymptotically stable walks and all boundary sequences satisfying $g_n=o(c_n)$. Furthermore, we prove that the rescaled random walk conditioned to stay above the boundary up to time $n$ converges, as $n\to\infty$, towards the stable meander.


title: Multisensor Online Transfer Learning for 3D LiDAR-based Human Detection with a Mobile Robot

abstract: Human detection and tracking is an essential task for service robots, where the combined use of multiple sensors has potential advantages that are yet to be exploited. In this paper, we introduce a framework allowing a robot to learn a new 3D LiDAR-based human classifier from other sensors over time, taking advantage of a multisensor tracking system. The main innovation is the use of different detectors for existing sensors (i.e. RGB-D camera, 2D LiDAR) to train, online, a new 3D LiDAR-based human classifier, exploiting a so-called trajectory probability. Our framework uses this probability to check whether new detections belongs to a human trajectory, estimated by different sensors and/or detectors, and to learn a human classifier in a semi-supervised fashion. The framework has been implemented and tested on a real-world dataset collected by a mobile robot. We present experiments illustrating that our system is able to effectively learn from different sensors and from the environment, and that the performance of the 3D LiDAR-based human classification improves with the number of sensors/detectors used.


title: The detection of the blazar S4 0954+65 at very-high-energy with the MAGIC telescopes during an exceptionally high optical state

abstract: The very-high-energy (VHE, $\gtrsim 100$ GeV) $\gamma$-ray MAGIC observations of the blazar S4 0954+65, were triggered by an exceptionally high flux state of emission in the optical. This blazar has a disputed redshift of z=0.368 or z$\geqslant$0.45 and an uncertain classification among blazar subclasses. The exceptional source state described here makes for an excellent opportunity to understand physical processes in the jet of S4 0954+65 and thus contribute to its classification. We investigate the multiwavelength (MWL) light curve and spectral energy distribution (SED) of the S4 0954+65 blazar during an enhanced state in February 2015 and put it in context with possible emission scenarios. We collect photometric data in radio, optical, X-ray, and $\gamma$ ray. We study both the optical polarization and the inner parsec-scale jet behavior with 43 GHz data. Observations with the MAGIC telescopes led to the first detection of S4 0954+65 at VHE. Simultaneous data with Fermi-LAT at high energy $\gamma$ ray\ (HE, 100 MeV < E < 100 GeV) also show a period of increased activity. Imaging at 43 GHz reveals the emergence of a new feature in the radio jet in coincidence with the VHE flare. Simultaneous monitoring of the optical polarization angle reveals a rotation of approximately 100$^\circ$. (...) The broadband spectrum can be modeled with an emission mechanism commonly invoked for flat spectrum radio quasars, i.e. inverse Compton scattering on an external soft photon field from the dust torus, also known as external Compton. The light curve and SED phenomenology is consistent with an interpretation of a blob propagating through a helical structured magnetic field and eventually crossing a standing shock in the jet, a scenario typically applied to flat spectrum radio quasars (FSRQs) and low-frequency peaked BL Lac objects (LBL).


title: Secure heterodyne-based quantum random number generator at 17 Gbps

abstract: Random numbers are commonly used in many different fields, ranging from simulations in fundamental science to security applications. In some critical cases, as Bell's tests and cryptography, the random numbers are required to be both secure (i.e. known only by the legitimate user) and to be provided at an ultra-fast rate (i.e. larger than Gbit/s). However, practical generators are usually considered trusted, but their security can be compromised in case of imperfections or malicious external actions. In this work we introduce an efficient protocol which guarantees security and speed in the generation. We propose a novel source-device-independent protocol based on generic Positive Operator Valued Measurements and then we specialize the result to heterodyne measurements. The security of the generated numbers is proven without any assumption on the source, which can be even fully controlled by an adversary. Furthermore, we experimentally implemented the protocol by exploiting heterodyne measurements, reaching an unprecedented secure generation rate of 17.42 Gbit/s, without the need to take into account finite-size effects. Our device combines simplicity, ultrafast-rates and high security with low cost components, paving the way to new practical solutions for random number generation.


title: Cosmic String Detection with Tree-Based Machine Learning

abstract: We explore the use of random forest and gradient boosting, two powerful tree-based machine learning algorithms, for the detection of cosmic strings in maps of the cosmic microwave background (CMB), through their unique Gott-Kaiser-Stebbins effect on the temperature anisotropies.The information in the maps is compressed into feature vectors before being passed to the learning units. The feature vectors contain various statistical measures of processed CMB maps that boost the cosmic string detectability. Our proposed classifiers, after training, give results improved over or similar to the claimed detectability levels of the existing methods for string tension, $G\mu$. They can make $3\sigma$ detection of strings with $G\mu \gtrsim 2.1\times 10^{-10}$ for noise-free, $0.9'$-resolution CMB observations. The minimum detectable tension increases to $G\mu \gtrsim 3.0\times 10^{-8}$ for a more realistic, CMB S4-like (II) strategy, still a significant improvement over the previous results.

categories: astro-ph.CO astro-ph.IM physics.data-an stat.ML


title: Holomorphic Hartree-Fock Theory: The Nature of Two-Electron Problems

abstract: We explore the existence and behaviour of holomorphic restricted Hartree-Fock (h-RHF) solutions for two-electron problems. Through algebraic geometry, the exact number of solutions with $n$ basis functions is rigorously identified as $\frac{1}{2}(3^n - 1)$, proving that states must exist for all molecular geometries. A detailed study on the h-RHF states of HZ (STO-3G) then demonstrates both the conservation of holomorphic solutions as geometry or atomic charges are varied and the emergence of complex h-RHF solutions at coalescence points. Using catastrophe theory, the nature of these coalescence points is described, highlighting the influence of molecular symmetry. The h-RHF states of HHeH$^{2+}$ and HHeH (STO-3G) are then compared, illustrating the isomorphism between systems with two electrons and two electron holes. Finally, we explore the h-RHF states of ethene (STO-3G) by considering the $\pi$ electrons as a two-electron problem, and employ NOCI to identify a crossing of the lowest energy singlet and triplet states at the perpendicular geometry.


title: Meshfree and efficient modelling of swimming cells

abstract: Locomotion in Stokes flow is an intensively-studied problem because it describes important biological phenomena such as the motility of many species' sperm, bacteria, algae and protozoa. Numerical computations can be challenging, particularly in three dimensions, due to the presence of moving boundaries and complex geometries; methods which combine ease-of-implementation and computational efficiency are therefore needed. A recently-proposed method to discretise the regularised Stokeslet boundary integral equation without the need for a connected 'mesh' is applied to the inertialess locomotion problem in Stokes flow. The mathematical formulation and key aspects of the computational implementation in MATLAB/GNU Octave are described, followed by numerical experiments with biflagellate algae and multiple uniflagellate sperm swimming between no-slip surfaces, for which both swimming trajectories and flow fields are calculated. These computational experiments required minutes of time on modest hardware; an extensible implementation is provided in a github repository. The nearest neighbour discretisation dramatically improves convergence and robustness, a key challenge in extending the regularised Stokeslet method to complicated, three dimensional, biological fluid problems.


title: Local asymptotic self-similarity for heavy tailed harmonizable fractional L\'evy motions

abstract: In this work we characterize the local asymptotic self-similarity of harmonizable fractional L\'evy motions in the heavy tailed case. The corresponding tangent process is shown to be the harmonizable fractional stable motion. In addition, we provide sufficient conditions for existence of harmonizable fractional L\'evy motions.