title: A Subject-Specific Four-Degree-of-Freedom Foot Interface to Control a Robot Arm

abstract: In robotic surgery, the surgeon controls robotic instruments using dedicated interfaces. One critical limitation of current interfaces is that they are designed to be operated by only the hands. This means that the surgeon can only control at most two robotic instruments at one time while many interventions require three instruments. This paper introduces a novel four-degree-of-freedom foot-machine interface which allows the surgeon to control a third robotic instrument using the foot, giving the surgeon a "third hand". This interface is essentially a parallel-serial hybrid mechanism with springs and force sensors. Unlike existing switch-based interfaces that can only un-intuitively generate motion in discrete directions, this interface allows intuitive control of a slave robotic arm in continuous directions and speeds, naturally matching the foot movements with dynamic force & position feedbacks. An experiment with ten naive subjects was conducted to test the system. In view of the significant variance of motion patterns between subjects, a subject-specific mapping from foot movements to command outputs was developed using Independent Component Analysis (ICA). Results showed that the ICA method could accurately identify subjects' foot motion patterns and significantly improve the prediction accuracy of motion directions from 68% to 88% as compared with the forward kinematics-based approach. This foot-machine interface can be applied for the teleoperation of industrial/surgical robots independently or in coordination with hands in the future.


title: Defect removal by solvent vapor annealing in thin films of lamellar diblock copolymers

abstract: Solvent vapor annealing (SVA) is known to be a simple, low-cost and highly efficient technique to produce defect-free diblock copolymer (BCP) thin films. Not only can the solvent weaken the BCP segmental interactions, but it can vary the characteristic spacing of the BCP microstructures. We carry out systematic theoretical studies on the effect of adding solvent into lamellar BCP thin films on the defect removal close to the BCP order-disorder transition. We find that the increase of the lamellar spacing, as is induced by addition of solvent, facilitates more efficient removal of defects. The stability of a particular defect in a lamellar BCP thin film is given in terms of two key controllable parameters: the amount of BCP swelling and solvent evaporation rate. Our results highlight the SVA mechanism for obtaining defect-free BCP thin films, as is highly desired in nanolithography and other industrial applications.


title: Rank-3 antisymmetric matter on 5-brane webs

abstract: We discuss Type IIB 5-brane configurations for 5d $\mathcal{N}=1$ gauge theories with hypermultiplets in the rank-3 antisymmetric representation and with various other hypermultiplets, which flow to a UV fixed point at the infinite coupling. We propose 5-brane web diagrams for the theories of $SU(6)$ and $Sp(3)$ gauge groups with rank-3 antisymmetric matter and check our proposed 5-brane webs against several consistency conditions implied from the one-loop corrected prepotential. Using the obtained 5-brane webs for rank-3 antisymmetric matter, we apply the topological vertex method to compute the partition function for one of these $SU(6)$ gauge theories.


title: Multi-Prototype Networks for Unconstrained Set-based Face Recognition

abstract: In this paper, we study the challenging unconstrained set-based face recognition problem where each subject face is instantiated by a set of media (images and videos) instead of a single image. Naively aggregating information from all the media within a set would suffer from the large intra-set variance caused by heterogeneous factors (e.g., varying media modalities, poses and illuminations) and fail to learn discriminative face representations. A novel Multi-Prototype Network (MPNet) model is thus proposed to learn multiple prototype face representations adaptively from the media sets. Each learned prototype is representative for the subject face under certain condition in terms of pose, illumination and media modality. Instead of handcrafting the set partition for prototype learning, MPNet introduces a Dense SubGraph (DSG) learning sub-net that implicitly untangles inconsistent media and learns a number of representative prototypes. Qualitative and quantitative experiments clearly demonstrate superiority of the proposed model over state-of-the-arts.


title: Anomalous $H\to ZZ \to 4\ell$ decay and its interference effects at the LHC

abstract: We calculate the spinor helicity amplitudes of anomalous $H\to ZZ \to 4\ell$ decay. After embedding these analytic formulas into the $\texttt{MCFM}$ package, we study the interference effects between the anomalous $gg\to H\to ZZ \to 4\ell$ process and the SM processes, which are indispensable in the Higgs off-shell region. Subsequently, the constraints on the anomalous couplings are estimated using LHC experimental data.


title: Giant Negative Thermal Expansion Induced by the Synergistic Effects of Ferroelectrostriction and Spin-Crossover in PbTiO3-Based Perovskites

abstract: The discovery of unusual negative thermal expansion (NTE) provides the opportunity to control the common but much desired property of thermal expansion, which is valuable not only in scientific interests but also in practical applications. However, most of the available NTE materials are limited to a narrow temperature range, and the NTE effect is generally weakened by means of various modifications. Here, we report an enhanced NTE effect that occurs over a wide temperature range (\alpha V = -5.24 * 10-5 {\deg}C^-1, 25-575 {\deg} C), and this NTE effect is accompanied by an abnormal enhanced tetragonality, a large spontaneous polarization, and a G-type antiferromagnetic ordering in the present perovskite-type ferroelectric of (1-x)PbTiO3-xBiCoO3. Specifically, for the composition of 0.5PbTiO3-0.5BiCoO3, an extensive volumetric contraction of ~4.8 % has been observed near the Curie temperature of 700 {\deg}C, which represents the highest level in PbTiO3-based ferroelectrics. According to our experimental and theoretical results, the giant NTE originates from a synergistic effect of the ferroelectrostriction and spin-crossover of cobalt on the crystal lattice. The actual NTE mechanism is contrasted with previous functional NTE materials, in which the NTE is simply coupled with one ordering such as electronic, magnetic, or ferroelectric ordering. The present study sheds light on the understanding of NTE mechanisms and it attests that NTE could be simultaneouslycoupled with different orderings, which will pave a new way toward the design of large NTE materials.


title: Steady streaming viscometry of Newtonian liquids in microfluidic devices

abstract: We report a novel technique capable of measuring the kinematic shear viscosity of Newtonian liquids with steady streaming in microfluidic devices. This probe-free microrheological method utilizes sub-kilohertz liquid oscillation frequencies around a cylindrical obstacle, ensuring that the inner streaming layer is comparable in size to the cylinder radius. To calibrate the viscometer, the evolution of the inner streaming layer as a function of oscillation frequency for a liquid of known viscosity is characterized using standard particle tracking techniques. Once calibrated, we show how the steady streaming viscometer can be used to measure low-viscosity liquids and volatile liquids.


title: Scaling laws for frictional granular materials confined by constant pressure under oscillatory shear

abstract: Herein, we numerically study the rheology of a two-dimensional frictional granular system confined by constant pressure under oscillatory shear. Several scaling laws for the storage and loss moduli against the scaled strain amplitude have been found. The scaling laws in plastic regime for large strain amplitude can be understood by the angular distributions of the contact force. The scaling exponents are estimated by considering the physical mechanism.


title: Scaling Limits of Wide Neural Networks with Weight Sharing: Gaussian Process Behavior, Gradient Independence, and Neural Tangent Kernel Derivation

abstract: Several recent trends in machine learning theory and practice, from the design of state-of-the-art Gaussian Process to the convergence analysis of deep neural nets (DNNs) under stochastic gradient descent (SGD), have found it fruitful to study wide random neural networks. Central to these approaches are certain scaling limits of such networks. We unify these results by introducing a notion of a straightline \emph{tensor program} that can express most neural network computations, and we characterize its scaling limit when its tensors are large and randomized. From our framework follows (1) the convergence of random neural networks to Gaussian processes for architectures such as recurrent neural networks, convolutional neural networks, residual networks, attention, and any combination thereof, with or without batch normalization; (2) conditions under which the \emph{gradient independence assumption} -- that weights in backpropagation can be assumed to be independent from weights in the forward pass -- leads to correct computation of gradient dynamics, and corrections when it does not; (3) the convergence of the Neural Tangent Kernel, a recently proposed kernel used to predict training dynamics of neural networks under gradient descent, at initialization for all architectures in (1) without batch normalization. Mathematically, our framework is general enough to rederive classical random matrix results such as the semicircle and the Marchenko-Pastur laws, as well as recent results in neural network Jacobian singular values. We hope our work opens a way toward design of even stronger Gaussian Processes, initialization schemes to avoid gradient explosion/vanishing, and deeper understanding of SGD dynamics in modern architectures.

categories: cs.NE cond-mat.dis-nn cs.LG math-ph math.MP stat.ML


title: Structure and dynamics of a self-propelled semiflexible filament

abstract: We investigate structural and dynamical properties of a self-propelled filament using coarse-grained Brownian dynamics simulations. A self-propulsion force is applied along the bond vectors, i.e., tangent to the filament and their locations are considered in two different manners. In case one, force is applied to all beads of the filament, which is termed as homogeneous self-propulsion. Here, we obtain a monotonic decrease in the flexibility of the filament with P\'eclet number. Hence, radius of gyration also displays the same trend. Moreover, the radius of gyration of the filament shows universal dependence for various bending rigidities with flexure number. The effective diffusivity of the filament shows enhancement with the active force and it increases linearly with force and bending rigidity. In case two, self-propulsion force is applied only to few bond vectors. The location of active forces is chosen in a periodic manner starting from the tail of the filament and leaving the front end without force. In this case, filament acquires various structures such as rod-like, helical, circular, and folded states. The transition from several states is understood in terms of tangent-tangent correlation, bending energy and torsional order parameter. The helical state is identified through a crossover from exponential to oscillatory behavior of the tangent-tangent correlation. A sudden increase in the bending energy separates a helical to a folded states of the filament.


title: Trajectory Optimization in UAV-Assisted Cellular Networks under Mission Duration Constraint

abstract: In this paper, we address the problem of finding the optimal trajectory for an unmanned aerial vehicle (UAV) for improving the wireless coverage of a terrestrial cellular network. In particular, we consider a UAV that is tasked to travel from one point to another within a given time constraint, and it can also simultaneously assist the cellular network by providing wireless coverage during its mission. Considering an interference limited downlink of a cellular network, we formulate an optimization problem for maximizing the proportional-fair (PF) data rate of the cellular network and explore dynamic programming (DP) technique for finding the optimum UAV trajectory. We also explore the optimal UAV trajectories associated with maximum sum-rate and 5th percentile spectral efficiency (5pSE) rate and compare the capacity and coverage performance of the three approaches. Numerical simulations show that the maximum sum-rate trajectory provides the best per user capacity whereas, the optimal PF trajectory yields higher coverage probability than the other two trajectories. The optimal trajectories are generally infeasible to follow exactly as the UAVs can not take sharp turns due to kinematic constraints. Hence, we generate smooth trajectories using Bezier curve.


title: Wireless Traffic Prediction with Scalable Gaussian Process: Framework, Algorithms, and Verification

abstract: The cloud radio access network (C-RAN) is a promising paradigm to meet the stringent requirements of the fifth generation (5G) wireless systems. Meanwhile, wireless traffic prediction is a key enabler for C-RANs to improve both the spectrum efficiency and energy efficiency through load-aware network managements. This paper proposes a scalable Gaussian process (GP) framework as a promising solution to achieve large-scale wireless traffic prediction in a cost-efficient manner. Our contribution is three-fold. First, to the best of our knowledge, this paper is the first to empower GP regression with the alternating direction method of multipliers (ADMM) for parallel hyper-parameter optimization in the training phase, where such a scalable training framework well balances the local estimation in baseband units (BBUs) and information consensus among BBUs in a principled way for large-scale executions. Second, in the prediction phase, we fuse local predictions obtained from the BBUs via a cross-validation based optimal strategy, which demonstrates itself to be reliable and robust for general regression tasks. Moreover, such a cross-validation based optimal fusion strategy is built upon a well acknowledged probabilistic model to retain the valuable closed-form GP inference properties. Third, we propose a C-RAN based scalable wireless prediction architecture, where the prediction accuracy and the time consumption can be balanced by tuning the number of the BBUs according to the real-time system demands. Experimental results show that our proposed scalable GP model can outperform the state-of-the-art approaches considerably, in terms of wireless traffic prediction performance.


title: Explicit lower bounds on strong simulation of quantum circuits in terms of $T$-gate count

abstract: We investigate Clifford+$T$ quantum circuits with a small number of $T$-gates. Using the sparsification lemma, we identify time complexity lower bounds in terms of $T$-gate count below which a strong simulator would improve on the state-of-the-art $3$-SAT solving.


title: An efficient methodology to estimate the parameters of a two-dimensional chirp signal model

abstract: In various capacities of statistical signal processing two-dimensional (2-D) chirp models have been considered significantly, particularly in image processing$-$ to model gray-scale and texture images, magnetic resonance imaging, optical imaging etc. In this paper we address the problem of estimation of the unknown parameters of a 2-D chirp model under the assumption that the errors are independently and identically distributed (i.i.d.). The key attribute of the proposed estimation procedure is that it is computationally more efficient than the least squares estimation method. Moreover, the proposed estimators are observed to have the same asymptotic properties as the least squares estimators, thus providing computational effectiveness without any compromise on the efficiency of the estimators. We extend the propounded estimation method to provide a sequential procedure to estimate the unknown parameters of a 2-D chirp model with multiple components and under the assumption of i.i.d. errors we study the large sample properties of these sequential estimators. Simulation studies and a synthetic data analysis show that the proposed estimators perform satisfactorily.


title: Compositional differences between the component stars of eclipsing close binary systems showing chemical peculiarities

abstract: A spectroscopic study was carried out on the surface chemical abundances of CNO and several heavier elements in the primary and secondary components of 5 eclipsing close binaries around A-type (AR Aur, beta Aur, YZ Cas, WW Aur, and RR Lyn), in order to investigate the nature of chemical differences between both components (being comparatively slow rotators alike due to tidal synchronization). Regarding the systems comprising similar components, beta Aur and WW Aur were confirmed to exhibit no compositional difference between the primary and secondary both showing almost the same Am anomaly, though the chemical peculiarities seen in the component stars of AR Aur show distinct differences (HgMn star and Am star). In contrast, as to the systems (YZ Cas and RR Lyn) consisting of considerably different (A and early-F) components, the surface abundances are markedly different between the primary (Am) and secondary (normal). These observational results may indicate Teff-dependent characteristics regarding the chemical anomalies of non-magnetic stars on the upper main sequence: (1) In the effective temperature range of 10000K > Teff > 7000K, rotational velocity is the most important factor for determining the extent of Am peculiarity. (2) However, the emergence of Am phenomenon seems to have a lower Teff limit at ~7000K, below which no abundance anomaly is observed regardless of stellar rotation. (3) The transition from Am anomaly (mild deficiency in CNO) to HgMn anomaly (unusually large N depletion) is likely to take place as Teff increases from ~10000K to ~11000K.


title: Evolutionary Algorithms for the Chance-Constrained Knapsack Problem

abstract: Evolutionary algorithms have been applied to a wide range of stochastic problems. Motivated by real-world problems where constraint violations have disruptive effects, this paper considers the chance-constrained knapsack problem (CCKP) which is a variance of the binary knapsack problem. The problem aims to maximize the profit of selected items under a constraint that the knapsack capacity bound is violated with a small probability. To tackle the chance constraint, we introduce how to construct surrogate functions by applying well-known deviation inequalities such as Chebyshev's inequality and Chernoff bounds. Furthermore, we investigate the performance of several deterministic approaches and introduce a single- and multi-objective evolutionary algorithm to solve the CCKP. In the experiment section, we evaluate and compare the deterministic approaches and evolutionary algorithms on a wide range of instances. Our experimental results show that a multi-objective evolutionary algorithm outperforms its single-objective formulation for all instances and performance better than deterministic approaches according to the computation time. Furthermore, our investigation points out in which circumstances to favour Chebyshev's inequality or the Chernoff bound when dealing with the CCKP.


title: Efficient Cross-Validation for Semi-Supervised Learning

abstract: Manifold regularization, such as laplacian regularized least squares (LapRLS) and laplacian support vector machine (LapSVM), has been widely used in semi-supervised learning, and its performance greatly depends on the choice of some hyper-parameters. Cross-validation (CV) is the most popular approach for selecting the optimal hyper-parameters, but it has high complexity due to multiple times of learner training. In this paper, we provide a method to approximate the CV for manifold regularization based on a notion of robust statistics, called Bouligand influence function (BIF). We first provide a strategy for approximating the CV via the Taylor expansion of BIF. Then, we show how to calculate the BIF for general loss function,and further give the approximate CV criteria for model selection in manifold regularization. The proposed approximate CV for manifold regularization requires training only once, hence can significantly improve the efficiency of traditional CV. Experimental results show that our approximate CV has no statistical discrepancy with the original one, but much smaller time cost.


title: Interacting non-Hermitian ultracold atoms in a harmonic trap: Two-body exact solution and high-order exceptional point

abstract: We study interacting ultracold atoms in a three-dimensional (3D) harmonic trap with spin-selective dissipations, which can be effectively described by non-Hermitian parity-time ($\mathcal{PT}$) symmetric Hamiltonians. By solving the non-Hermitian two-body problem of spin-1/2 (spin-1) bosons in a 3D harmonic trap exactly, we find that the system can exhibit third-order (fifth-order) exceptional point (EP) with ultra-sensitive cube-root (fifth-root) spectral response due to interaction anisotropies in spin channels. We also present the general principle for the creation of high-order EPs and their spectral sensitivities with arbitrary particle number $N$ and arbitrary spin $s$. Generally, with spin-independent interactions, the EP order of bosons can be as high as $2Ns+1$, and the spectral response around EP can be as sensitive as $\sim \epsilon^{1/(2ks+1)}$ under a $k$-body interaction anisotropy. Moreover, we propose to detect the ultra-sensitive spectral response through the probability dynamics of certain state. These results suggest a convenient route towards more powerful sensor devices in spinor cold atomic systems.


title: Gas-phase functionalization of macroscopic carbon nanotube fiber fabrics: reaction control, electrochemical properties and use for flexible supercapacitors

abstract: The assembly of aligned carbon nanotubes (CNT) into fibers (CNTF) is a convenient approach to exploit and apply the unique physico-chemical properties of CNTs in many fields. CNT functionalization has been extensively used for their implementation into composites and devices. However, CNTF functionalization is still in its infancy, due to the challenges associated with preservation of CNTF morphology. Here, we report a thorough study of the gas-phase functionalization of CNTF fabrics using ozone that was generated in situ from a UV-source. By contrast with liquid-based oxidation methods, this gas-phase approach preserves CNTF morphology, whilst notably increasing its hydrophilicity. The functionalized material is thoroughly characterized by Raman, XPS, TEM and SEM. Its newly acquired hydrophilicity enables CNTF electrochemical characterization in aqueous media, which was not possible for the pristine material. Through comparison of electrochemical measurements in aqueous electrolyte and ionic liquid we decouple the effects of functionalization on pseudocapacitive reactions and quantum capacitance. The functionalized CNTF fabric is successfully used as active material and current collector in all-solid supercapacitor flexible devices with ionic liquid-based polymer electrolyte.


title: Controlling injection using a magnetic field to produce sub-femtosecond bunches in the laser wakefield accelerator

abstract: It is shown that electron injection into a laser-driven plasma bubble can be manipulated by applying an external magnetic field in the presence of a plasma density gradient. The down-ramp of the density-tailored plasma locally reduces the plasma wave phase velocity, which triggers injection. The longitudinal magnetic field dynamically induces an expanding hole in the electron density distribution at the rear of the wake bubble, which reduces the peak electron velocity in its vicinity. Electron injection is suppressed when the electron velocity drops below the phase velocity, which depends on the size of the density hole. This enables the start and end of electron injection to be independently controlled, which allows generation of sub-femtosecond electron bunches with peak currents of a few kilo-Ampere, for an applied magnetic field of $\sim 10$ Tesla.


title: On Trivial Extensions and Higher Preprojective Algebras

abstract: In this paper, we show that for a Koszul $n$-homogeneous algebra $\Lambda$, the quadratic dual of certain twisted trivial extension is the $(n+1)$-preprojective algebra of its quadratic dual, that is, $ (\Delta_{\nu}\Lambda)^{!,op} \simeq\Pi( \Lambda^{ !, op })$. This is applied to the $\tau$-slice algebras of stable $n$-translation algebras and gives a noncommutative version of Bernstein-Gelfand-Gelfand correspondence for such algebras.


title: Lepton flavor violation induced by neutral and doubly-charged scalars at future lepton colliders

abstract: New physics scenarios beyond the Standard Model (SM) for neutrino mass mechanism often necessitate the existence of a neutral scalar $H$ and/or doubly-charged scalar $H^{\pm\pm}$, which couple to the SM charged leptons in a flavor violating way, while evading all existing constraints. Such scalars could be effectively produced at future lepton colliders like CEPC, ILC, FCC-ee and CLIC, either on-shell or off-shell, and induce striking charged lepton flavor violating (LFV) signals. We find that a large parameter space of the scalar masses and the LFV couplings can be probed at lepton colliders, well beyond the current low-energy constraints in the lepton sector. The neutral scalar explanation of the muon $g-2$ anomaly could also be directly tested.


title: Distributed Online Linear Regression

abstract: We study online linear regression problems in a distributed setting, where the data is spread over a network. In each round, each network node proposes a linear predictor, with the objective of fitting the \emph{network-wide} data. It then updates its predictor for the next round according to the received local feedback and information received from neighboring nodes. The predictions made at a given node are assessed through the notion of regret, defined as the difference between their cumulative network-wide square errors and those of the best off-line network-wide linear predictor. Various scenarios are investigated, depending on the nature of the local feedback (full information or bandit feedback), on the set of available predictors (the decision set), and the way data is generated (by an oblivious or adaptive adversary). We propose simple and natural distributed regression algorithms, involving, at each node and in each round, a local gradient descent step and a communication and averaging step where nodes aim at aligning their predictors to those of their neighbors. We establish regret upper bounds typically in ${\cal O}(T^{3/4})$ when the decision set is unbounded and in ${\cal O}(\sqrt{T})$ in case of bounded decision set.


title: Indirect Microwave Holography with Resolution Enhancement in Metallic Imaging

abstract: The development of compact indirect microwave holographic set-up by the implementation of low cost, specifically designed directive antennae as transmitter and receiver is proposed. Microwave holograms are recorded by 2D scanning over a plane using motorized translation stage. The recorded interference pattern i.e. holograms are then processed numerically to reconstruct the amplitude and phase information employing the angular spectrum diffraction method. The quality of the reconstructed amplitude image is further enhanced through the deep neural network, in order to combat with the low resolution of reconstructed images. The qualitative experimental results exploit the possibility of developing the miniaturized, and low cost indirect microwave holographic system for near field applications.


title: Observational constraints on the oscillating dark energy cosmologies

abstract: In this study we combine the background and the growth rate data in order to study the ability of the two oscillating dark energy parameterizations, to fit the observational data. Using the likelihood and MCMC method we try to explore the posterior space and put constraints on the free parameters of the models. Based on the values of the well known Akaike and Bayesian information criteria we find that both of oscillating dark energy models considered in this work are disfavored by the combined (background+growth rate) data. Although using the expansion data we can not reject oscillating dark energy models, the combined analysis provides strong evidences against these models.


title: On some properties of relative capacity and thinness in weighted variable exponent Sobolev spaces

abstract: In this paper, we define weighted relative $p(.)$-capacity and discuss properties of capacity in the space $W_{\vartheta }^{1,p(.)}(\mathbb{R}^{n}).$ Also, we investigate some properties of weighted variable Sobolev capacity. It is shown that there is a relation between these two capacities. Moreover, we introduce a thinness in sense to this new defined relative capacity and prove an equivalence statement for this thinness.


title: Cavities, shocks and a cold front around 3C 320

abstract: We present results obtained from the analysis of a total of 110 ks Chandra observations of 3C 320 FR II radio galaxy, located at the centre of a cluster of galaxies at a redshift $z=0.342$. A pair of X-ray cavities have been detected at an average distance of $\sim$38 kpc along the East and West directions with the cavity energy, age and total power equal to $\sim$7.7$\times$10$^{59}$ erg, $\sim$7$\times$10$^7$ yr and $\sim$3.5$\times$10$^{44}$ erg s$^{-1}$, respectively. The cooling luminosity within the cooling radius of $\sim$100 kpc was found to be $L_{cool} \sim8.5\times10^{43}$ erg s$^{-1}$. Comparison of these two estimates implies that the cavity power is sufficiently high to balance the radiative loss. A pair of weak shocks have also been evidenced at distances of $\sim$47 kpc and $\sim$76 kpc surrounding the radio bubbles. Using the observed density jumps of $\sim$1.8 and $\sim$2.1 at shock locations along the East and West directions, we estimate the Mach numbers ($\mathcal{M}$) to be $\sim$1.6 and $\sim$1.8, respectively. A sharp surface brightness edge was also detected at relatively larger radius ($\sim$80 kpc) along the South direction. Density jump at this surface brightness edge was estimated to be $\sim$1.6 and is probably due to the presence of a cold front in this cluster. The far-infrared luminosity yielded the star formation rate of 51 M$_{\odot}$ yr$^{-1}$ and is 1/4$^{th}$ of the cooling rate ($\dot{M}$ $\sim$ 192 M$_{\odot}$ yr$^{-1}$).


title: Sample-Optimal Parametric Q-Learning Using Linearly Additive Features

abstract: Consider a Markov decision process (MDP) that admits a set of state-action features, which can linearly express the process's probabilistic transition model. We propose a parametric Q-learning algorithm that finds an approximate-optimal policy using a sample size proportional to the feature dimension $K$ and invariant with respect to the size of the state space. To further improve its sample efficiency, we exploit the monotonicity property and intrinsic noise structure of the Bellman operator, provided the existence of anchor state-actions that imply implicit non-negativity in the feature space. We augment the algorithm using techniques of variance reduction, monotonicity preservation, and confidence bounds. It is proved to find a policy which is $\epsilon$-optimal from any initial state with high probability using $\widetilde{O}(K/\epsilon^2(1-\gamma)^3)$ sample transitions for arbitrarily large-scale MDP with a discount factor $\gamma\in(0,1)$. A matching information-theoretical lower bound is proved, confirming the sample optimality of the proposed method with respect to all parameters (up to polylog factors).


title: ProtoDUNE and a Dual-phase LArTPC

abstract: The four 10 kt Liquid Argon Time Projection Chambers (LArTPCs) of the future DUNE experiment will enable precise measurements of the oscillation parameters and the discovery of CP violation for leptons, thanks to their excellent 3D imaging capabilities and calorimetric capabilities. One or more modules of the DUNE detector may exploit a Dual Phase (DP) LArTPC that, relying on the extraction of the charge produced in the liquid volume and its subsequent multiplication in argon gas, may offer a robust and competitive signal-to-noise ratio and a fully active volume. In 2018 and 2019, the ProtoDUNE experiment at CERN will validate the designs of the DUNE far detector, showing the feasibility of large scale Dual-Phase LArTPC and providing precious insight on the DUNE physics potential.


title: New coordinates for a simpler canonical derivation of the Hawking effect

abstract: In order to achieve a Hamiltonian-based canonical derivation of the Hawking effect, one usually faces multiple hurdles. Firstly, the spacetime foliation using Schwarzschild time does not lead to hyper-surfaces which are always spacelike. Secondly, the null coordinates which are frequently used in covariant approach, do not lead to a true matter Hamiltonian. Recently, an exact canonical derivation was presented using the so-called near-null coordinates. However, there too one faces the difficulty of having to deal with non-vanishing matter diffeomorphism generator as the spatial decomposition involves a non-zero shift vector. Here we introduce a new set of coordinates which allows one to perform an exact canonical derivation of Hawking effect without having to deal with matter diffeomorphism generator.


title: On the Expressive Power of Kernel Methods and the Efficiency of Kernel Learning by Association Schemes

abstract: We study the expressive power of kernel methods and the algorithmic feasibility of multiple kernel learning for a special rich class of kernels. Specifically, we define \emph{Euclidean kernels}, a diverse class that includes most, if not all, families of kernels studied in literature such as polynomial kernels and radial basis functions. We then describe the geometric and spectral structure of this family of kernels over the hypercube (and to some extent for any compact domain). Our structural results allow us to prove meaningful limitations on the expressive power of the class as well as derive several efficient algorithms for learning kernels over different domains.


title: Mathematical Notions vs. Human Perception of Fairness: A Descriptive Approach to Fairness for Machine Learning

abstract: Fairness for Machine Learning has received considerable attention, recently. Various mathematical formulations of fairness have been proposed, and it has been shown that it is impossible to satisfy all of them simultaneously. The literature so far has dealt with these impossibility results by quantifying the tradeoffs between different formulations of fairness. Our work takes a different perspective on this issue. Rather than requiring all notions of fairness to (partially) hold at the same time, we ask which one of them is the most appropriate given the societal domain in which the decision-making model is to be deployed. We take a descriptive approach and set out to identify the notion of fairness that best captures \emph{lay people's perception of fairness}. We run adaptive experiments designed to pinpoint the most compatible notion of fairness with each participant's choices through a small number of tests. Perhaps surprisingly, we find that the most simplistic mathematical definition of fairness---namely, demographic parity---most closely matches people's idea of fairness in two distinct application scenarios. This conclusion remains intact even when we explicitly tell the participants about the alternative, more complicated definitions of fairness, and we reduce the cognitive burden of evaluating those notions for them. Our findings have important implications for the Fair ML literature and the discourse on formalizing algorithmic fairness.


title: \'{E}tale coverings in codimension 1 with applications to Mori Dream Spaces

abstract: The present paper is devoted to developing relations between Galois \'etale coverings in codimension 1 and \'etale fundamental groups in codimension 1 of algebraic varieties, aimed to studying the topology of Mori dream spaces. In particular, the universal \'etale covering in codimension 1 of a non-degenerate toric variety and a canonical Galois \'etale covering in codimension 1 of a Mori dream space (MDS) are exhibited. Sufficient conditions for the latter being either still a MDS or the universal \'etale covering in codimension 1 are given. As an application, a canonical toric embedding of K3 universal coverings, of Enriques surfaces which are Mori dream, is described.


title: Constructing Antidictionaries in Output-Sensitive Space

abstract: A word $x$ that is absent from a word $y$ is called minimal if all its proper factors occur in $y$. Given a collection of $k$ words $y_1,y_2,\ldots,y_k$ over an alphabet $\Sigma$, we are asked to compute the set $\mathrm{M}^{\ell}_{y_{1}\#\ldots\#y_{k}}$ of minimal absent words of length at most $\ell$ of word $y=y_1\#y_2\#\ldots\#y_k$, $\#\notin\Sigma$. In data compression, this corresponds to computing the antidictionary of $k$ documents. In bioinformatics, it corresponds to computing words that are absent from a genome of $k$ chromosomes. This computation generally requires $\Omega(n)$ space for $n=|y|$ using any of the plenty available $\mathcal{O}(n)$-time algorithms. This is because an $\Omega(n)$-sized text index is constructed over $y$ which can be impractical for large $n$. We do the identical computation incrementally using output-sensitive space. This goal is reasonable when $||\mathrm{M}^{\ell}_{y_{1}\#\ldots\#y_{N}}||=o(n)$, for all $N\in[1,k]$. For instance, in the human genome, $n \approx 3\times 10^9$ but $||\mathrm{M}^{12}_{y_{1}\#\ldots\#y_{k}}|| \approx 10^6$. We consider a constant-sized alphabet for stating our results. We show that all $\mathrm{M}^{\ell}_{y_{1}},\ldots,\mathrm{M}^{\ell}_{y_{1}\#\ldots\#y_{k}}$ can be computed in $\mathcal{O}(kn+\sum^{k}_{N=1}||\mathrm{M}^{\ell}_{y_{1}\#\ldots\#y_{N}}||)$ total time using $\mathcal{O}(\mathrm{MaxIn}+\mathrm{MaxOut})$ space, where $\mathrm{MaxIn}$ is the length of the longest word in $\{y_1,\ldots,y_{k}\}$ and $\mathrm{MaxOut}=\max\{||\mathrm{M}^{\ell}_{y_{1}\#\ldots\#y_{N}}||:N\in[1,k]\}$. Proof-of-concept experimental results are also provided confirming our theoretical findings and justifying our contribution.


title: The Kolmogorov-Riesz Theorem and Some Compactness Criterions of Bounded Subsets in Weighted Variable Exponent Amalgam and Sobolev Spaces

abstract: We study totally bounded subsets in weighted variable exponent amalgam and Sobolev spaces. Moreover, this paper includes several detailed generalized results of some compactness criterions in these spaces.


title: Detailed study of detection method for ultra-low-frequency gravitational waves with pulsar spin-down rate statistics

abstract: A new detection method for gravitational waves (GWs) with ultra-low frequencies ($f_{\rm GW} \lesssim 10^{-10}~{\rm Hz}$), which is much lower than the range of pulsar timing arrays (PTAs), was proposed in Yonemaru et al. (2016). This method utilizes the statistical properties of spin-down rates of milli-second pulsars (MSPs) and the sensitivity was evaluated in Yonemaru et al. (2018). There, some simplifying assumptions, such as neglect of the "pulsar term" and spatially uniform distribution of MSPs, were adopted and the sensitivity on the time derivative of GW amplitude was estimated to be $10^{-19}~{\rm s}^{-1}$ independent of the direction, polarization and frequency of GWs. In this paper, extending the previous analysis, realistic simulations are performed to evaluate the sensitivity more reasonably. We adopt a model of 3-dimensional pulsar distribution in our Galaxy and take the pulsar term into account. As a result, we obtain expected sensitivity as a function of the direction, polarization and frequency of GWs. The dependence on GW frequency is particularly significant and the sensitivity becomes worse by a few orders for $< 10^{-12}~{\rm Hz}$ compared to the previous estimates.


title: Strain Controlling Catalytic Efficiency of Water Oxidation for Ni1-xFexOOH alloy

abstract: A catalyst surface may be exposed to strain due to application of load or interfacing with a substrate with large lattice mismatch. In order to test the effect of strain on catalytic efficiency, we use Density Functional Theory +U (DFT+U) to model water oxidation on expanded and contracted surfaces of the Ni1-xFexOOH alloy, one of the best known catalysts for water splitting. We find that a low amount of iron content has a similar effect as of applying compressive strain. Due to high oxidation state of Fe4+ at the active site, the Fe-O bond is shorter than in pure FeOOH, which is beneficial for extracting electrons from states delocalized on Fe and Ni atoms. At 33% of Fe content the efficiency is even better since Fe3+ is at the active site and can easily change oxidation state during the reaction. However, the efficiency drops at higher iron percentages since the surface is unstable and may form FeOOH aggregates. We find that the best performance is obtained at 33% iron content and 5% expansion. Therefore, in addition to changing iron content, strain can also be used as a control handle for improving water splitting catalysis.


title: Concurrent Computing with Shared Replicated Memory

abstract: The behavioural theory of concurrent systems states that any concurrent system can be captured by a behaviourally equivalent concurrent Abstract State Machine (cASM). While the theory in general assumes shared locations, it remains valid, if different agents can only interact via messages, i.e. sharing is restricted to mailboxes. There may even be a strict separation between memory managing agents and other agents that can only access the shared memory by sending query and update requests to the memory agents. This article is dedicated to an investigation of replicated data that is maintained by a memory management subsystem, whereas the replication neither appears in the requests nor in the corresponding answers. We show how the behaviour of a concurrent system with such a memory management can be specified using concurrent communicating ASMs. We provide several refinements of a high-level ground model addressing different replication policies and internal messaging between data centres. For all these refinements we analyse their effects on the runs such that decisions concerning the degree of consistency can be consciously made.


title: SaGe: Web Preemption for Public SPARQL Query Services

abstract: To provide stable and responsive public SPARQL query services, data providers enforce quotas on server usage. Queries which exceed these quotas are interrupted and deliver partial results. Such interruption is not an issue if it is possible to resume queries execution afterward. Unfortunately, there is no preemption model for the Web that allows for suspending and resuming SPARQL queries. In this paper, we propose SaGe: a SPARQL query engine based on Web preemption. SaGe allows SPARQL queries to be suspended by the Web server after a fixed time quantum and resumed upon client request. Web preemption is tractable only if its cost in time is negligible compared to the time quantum. The challenge is to support the full SPARQL query language while keeping the cost of preemption negligible. Experimental results demonstrate that SaGe outperforms existing SPARQL query processing approaches by several orders of magnitude in term of the average total query execution time and the time for first results.


title: Self-similar expansion model of cylindrical flux ropes combined with Alfven wave current system

abstract: Magnetic flux ropes in space are generally connected to some regions electromagnetically. We consider the whole closed current system of the expanding flux ropes including the electric current associated with them. By combining the theories regarding the self-similar expansion of cylindrical flux ropes and the Alfven wave current system, we examine conditions under which the electric current matches. These matching conditions are satisfied when the time dependence of the current flowing in the closed circuit agrees with that which maintains the expanding flux rope. In consequence, we encountered three possible forms of expansion. The two-step eruption of solar filaments may be interpreted as a transition from one form of expansion to another. If this process works, increasing the diffusion outside of the flux rope is necessary to trigger the transition.


title: An overlapping decomposition framework for wave propagation in heterogeneous and unbounded media: Formulation, analysis, algorithm, and simulation

abstract: A natural medium for wave propagation comprises a coupled bounded heterogeneous region and an unbounded homogeneous free-space. Frequency-domain wave propagation models in the medium, such as the variable coefficient Helmholtz equation, include a faraway decay radiation condition (RC). It is desirable to develop algorithms that incorporate the full physics of the heterogeneous and unbounded medium wave propagation model, and avoid an approximation of the RC. In this work we first present and analyze an overlapping decomposition framework that is equivalent to the full-space heterogeneous-homogenous continuous model, governed by the Helmholtz equation with a spatially dependent refractive index and the RC. Our novel overlapping framework allows the user to choose two free boundaries, and gain the advantage of applying established high-order finite and boundary element methods (FEM and BEM) to simulate an equivalent coupled model. The coupled model comprises auxiliary interior bounded heterogeneous and exterior unbounded homogeneous Helmholtz problems. A smooth boundary can be chosen for simulating the exterior problem using a spectrally accurate BEM, and a simple boundary can be used to setup a high-order FEM for the interior problem. Thanks to the spectral accuracy of the exterior computational model, the resulting coupled system in the overlapping region is relatively very small. Using the decomposed equivalent framework, we develop a novel overlapping FEM-BEM algorithm for simulating the acoustic or electromagnetic wave propagation in two dimensions. Our FEM-BEM algorithm for the full-space model incorporates the RC exactly. Numerical experiments demonstrate the efficiency of the FEM-BEM approach for simulating smooth and non-smooth wave fields, with the latter induced by a complex heterogeneous medium and a discontinuous refractive index.


title: SECTOR: A Neural Model for Coherent Topic Segmentation and Classification

abstract: When searching for information, a human reader first glances over a document, spots relevant sections and then focuses on a few sentences for resolving her intention. However, the high variance of document structure complicates to identify the salient topic of a given section at a glance. To tackle this challenge, we present SECTOR, a model to support machine reading systems by segmenting documents into coherent sections and assigning topic labels to each section. Our deep neural network architecture learns a latent topic embedding over the course of a document. This can be leveraged to classify local topics from plain text and segment a document at topic shifts. In addition, we contribute WikiSection, a publicly available dataset with 242k labeled sections in English and German from two distinct domains: diseases and cities. From our extensive evaluation of 20 architectures, we report a highest score of 71.6% F1 for the segmentation and classification of 30 topics from the English city domain, scored by our SECTOR LSTM model with bloom filter embeddings and bidirectional segmentation. This is a significant improvement of 29.5 points F1 compared to state-of-the-art CNN classifiers with baseline segmentation.


title: Analysis of the Block Coordinate Descent Method for Linear Ill-Posed Problems

abstract: Block coordinate descent (BCD) methods approach optimization problems by performing gradient steps along alternating subgroups of coordinates. This is in contrast to full gradient descent, where a gradient step updates all coordinates simultaneously. BCD has been demonstrated to accelerate the gradient method in many practical large-scale applications. Despite its success no convergence analysis for inverse problems is known so far. In this paper, we investigate the BCD method for solving linear inverse problems. As main theoretical result, we show that for operators having a particular tensor product form, the BCD method combined with an appropriate stopping criterion yields a convergent regularization method. To illustrate the theory, we perform numerical experiments comparing the BCD and the full gradient descent method for a system of integral equations. We also present numerical tests for a non-linear inverse problem not covered by our theory, namely one-step inversion in multi-spectral X-ray tomography.


title: Fibonacci Sequences And Real Quadratic p-Rational Fields

abstract: We study the p-rationality of real quadratic fields in terms of generalized Fibonacci numbers and their periods modulo positive integers.


title: Sparsity pattern of the self-energy for classical and quantum impurity problems

abstract: We prove that for various impurity models, in both classical and quantum settings, the self-energy matrix is a sparse matrix with a sparsity pattern determined by the impurity sites. In the quantum setting, such a sparsity pattern has been known since Feynman. Indeed, it underlies several numerical methods for solving impurity problems, as well as many approaches to more general quantum many-body problems, such as the dynamical mean field theory. The sparsity pattern is easily motivated by a formal perturbative expansion using Feynman diagrams. However, to the extent of our knowledge, a rigorous proof has not appeared in the literature. In the classical setting, analogous considerations lead to a perhaps less-known result, i.e., that the precision matrix of a Gibbs measure of a certain kind differs only by a sparse matrix from the precision matrix of a corresponding Gaussian measure. Our argument for this result mainly involves elementary algebraic manipulations and is in particular non-perturbative. Nonetheless, the proof can be robustly adapted to various settings of interest in physics, including quantum systems (both fermionic and bosonic) at zero and finite temperature, non-equilibrium systems, and superconducting systems.


title: Interplay of the spin density wave and a possible Fulde-Ferrell-Larkin-Ovchinnikov state in $\mathrm{CeCoIn_5}$ in rotating magnetic field

abstract: The $d$-wave superconductor $\mathrm{CeCoIn_5}$ has been proposed as a strong candidate for supporting the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) state near the low-temperature boundary of its upper critical field. Neutron diffraction, however, finds spin-density-wave (SDW) order in this part of the phase diagram for field in the $a$-$b$ plane, and evidence for the SDW disappears as the applied field is rotated toward the tetragonal $c$ axis. It is important to understand the interplay between the SDW and a possible FFLO state in $\mathrm{CeCoIn_5}$, as the mere existence of an SDW does not necessary exclude an FFLO state. Here, based on a model constructed on the basis of available experiments, we show that an FFLO state competes with an SDW phase. The SDW state in $\mathrm{CeCoIn_5}$ is stabilized when the field is directed close to the $a$-$b$ plane. When the field is rotated toward the $c$ axis, the FFLO state emerges, and the SDW phase disappears. In the FFLO state, the nodal planes with extra quasiparticles (where the superconducting order parameter is zero) are perpendicular to the field, and in the SDW phase, the quasiparticle density of states is reduced. We test this model prediction by measuring heat transported by normal quasiparticles in the superconducting state. As a function of field, we observe a reduction of thermal conductivity for field close to the $a$-$b$ plane and an enhancement of thermal conductivity when field is close to the $c$ axis, consistent with theoretical expectations. Our modeling and experiments, therefore, indicate the existence of the FFLO state when field is parallel to the $c$ axis.


title: Quantum computing, Seifert surfaces and singular fibers

abstract: The fundamental group $\pi_1(L)$ of a knot or link $L$ may be used to generate magic states appropriate for performing universal quantum computation and simultaneously for retrieving complete information about the processed quantum states. In this paper, one defines braids whose closure is the $L$ of such a quantum computer model and computes their Seifert surfaces and the corresponding Alexander polynomial. In particular, some $d$-fold coverings of the trefoil knot, with $d=3$, $4$, $6$ or $12$, define appropriate links $L$ and the latter two cases connect to the Dynkin diagrams of $E_6$ and $D_4$, respectively. In this new context, one finds that this correspondence continues with the Kodaira's classification of elliptic singular fibers. The Seifert fibered toroidal manifold $\Sigma'$, at the boundary of the singular fiber $\tilde {E_8}$, allows possible models of quantum computing.


title: Characterizations of umbilic hypersurfaces in warped product manifolds

abstract: We consider closed orientable hypersurfaces in a wide class of warped product manifolds, which include space forms, deSitter-Schwarzschild and Reissner-Nordstr\"{o}m manifolds. By using a new integral formula or Brendle's Heintze-Karcher type inequality, we present some new characterizations of umbilic hypersurfaces. These results can be viewed as generalizations of the classical Jellet-Liebmann theorem and the Alexandrov theorem in Euclidean space.


title: Emergence and expansion of cosmic space in an accelerating BIon

abstract: We generalize the Padmanabhan [arXiv:hep-th/1206.4916] mechanism to an accelerating BIon and show that the difference between the number of degrees of freedom on the boundary surface and the number of degrees of freedom in a bulk region causes the accelerated expansion of a BIon. We also consider the evolution of a universe which emerges on this BIon, and obtain its Hubble parameter and energy density.


title: Particle motion driven by non-uniform thermodynamic forces

abstract: We present a complete reciprocal description of particle motion inside multi-component fluids that extends the conventional Onsager formulation of non-equilibrium transport to systems where the thermodynamic forces are non-uniform on the colloidal scale. Based on the dynamic length and time scale separation in suspensions, the particle flux is shown to be related to the volume-averaged coupling between the Stokes flow tensor and the thermodynamic force density acting on the fluid. The flux is then expressed in terms of thermodynamic quantities that can be computed from the interfacial properties and equation of state of the colloids. Our results correctly describe diffusion and sedimentation, and suggest that force-free phoretic motion can occur even in the absence of interfacial interactions, provided that the thermodynamic gradients are non-uniform at the colloidal surface. In particular, we derive an explicit hydrodynamic form for the phoretic force resulting from these non-uniform gradients. The form is validated by the recovery of the Henry function for electrophoresis and the Ruckenstein term for thermophoresis.


title: Modeling interactions between transportation networks and territories: a co-evolution approach

abstract: Interactions between transportation networks and territories are the subject of open scientific debates, in particular regarding the possible existence of structuring effects of networks, and linked to crucial practical issues of territorial development. We propose an entry on these through co-evolution, and more particularly by the modeling of co-evolution processes between transportation networks and territories. We construct a multi-disciplinary definition of co-evolution which is proper to territorial systems and which can be tested empirically. We then develop the lessons learnt from the development of two types of models, macroscopic interaction models in systems of cities and mesoscopic morphogenesis models through co-evolution. This research opens the perspective of multi-scale models that could be applied to territorial prospective.


title: Optimal measurement strategies for fast entanglement detection

abstract: With the advance of quantum information technology, the question of how to most efficiently test quantum circuits is becoming of increasing relevance. Here we introduce the statistics of lengths of measurement sequences that allows one to certify entanglement across a given bi-partition of a multi-qubit system over the possible sequence of measurements of random unknown states, and identify the best measurement strategies in the sense of the (on average) shortest measurement sequence of (multi-qubit) Pauli measurements. The approach is based on the algorithm of truncated moment sequences that allows one to deal naturally with incomplete information, i.e. information that does not fully specify the quantum state. We find that the set of measurements corresponding to diagonal matrix elements of the moment matrix of the state are particularly efficient. For symmetric states their number grows only like the third power of the number $N$ of qubits. Their efficiency grows rapidly with $N$, leaving already for $N=4$ less than a fraction $10^{-6}$ of randomly chosen entangled states undetected.


title: Low-energy limit of the O(4) quark-meson model

abstract: We study the generation of low-energy couplings induced by quantum fluctuations within the O(4)-symmetric quark-meson model. To this end, we compute the functional renormalization group flow of the linearly realized quark-meson model including higher-derivative interactions and subsequently transform the resulting effective action into a nonlinear effective pion action. The latter is referred to as the low-energy limit of the O(4) quark-meson model. The present study may be considered as a preparatory work for the dynamical generation of low-energy couplings from functional QCD fluctuations in order to determine meaningful renormalization scales for purely pionic models.


title: Task-based Augmented Contour Trees with Fibonacci Heaps

abstract: This paper presents a new algorithm for the fast, shared memory, multi-core computation of augmented contour trees on triangulations. In contrast to most existing parallel algorithms our technique computes augmented trees, enabling the full extent of contour tree based applications including data segmentation. Our approach completely revisits the traditional, sequential contour tree algorithm to re-formulate all the steps of the computation as a set of independent local tasks. This includes a new computation procedure based on Fibonacci heaps for the join and split trees, two intermediate data structures used to compute the contour tree, whose constructions are efficiently carried out concurrently thanks to the dynamic scheduling of task parallelism. We also introduce a new parallel algorithm for the combination of these two trees into the output global contour tree. Overall, this results in superior time performance in practice, both in sequential and in parallel thanks to the OpenMP task runtime. We report performance numbers that compare our approach to reference sequential and multi-threaded implementations for the computation of augmented merge and contour trees. These experiments demonstrate the run-time efficiency of our approach and its scalability on common workstations. We demonstrate the utility of our approach in data segmentation applications.


title: Weak solutions of semilinear elliptic equations with Leray-Hardy potential and measure data

abstract: We study existence and stability of solutions of (E 1) --$\Delta$u + $\mu$ |x| 2 u + g(u) = $\nu$ in $\Omega$, u = 0 on $\partial$$\Omega$, where $\Omega$ is a bounded, smooth domain of R N , N $\ge$ 2, containing the origin, $\mu$ $\ge$ -- (N --2) 2 4 is a constant, g is a nondecreasing function satisfying some integral growth assumption and $\nu$ is a Radon measure on $\Omega$. We show that the situation differs according $\nu$ is diffuse or concentrated at the origin. When g is a power we introduce a capacity framework to find necessary and sufficient condition for solvability.


title: Tutorial on the analytical calculation of optical forces on spherical particles in optical tweezers

abstract: Arthur Ashkin was awarded the 2018 Nobel prize in physics for the invention of optical tweezers. Since the first publication in 1986 Optical Tweezers have been used as a tool to measure forces and rheological properties of microscopic systems. For the calibration of these measurements, the knowledge of the forces is fundamental. However, it is still common to deduce the optical forces from assumptions based on the particle size with respect to the trapping laser wavelength. This show the necessity to develop a complete and accurate electromagnetic model that does not depend on early approximations of the force model. Furthermore, the model we developed has several advantages, such as morphology--dependent resonances, size dependence for large spheres, and multipole effects from smaller particles, just to name a few. In this tutorial, we review and discuss the physical modeling of optical forces in optical tweezers, which are the resultant forces exerted by a trapping beam on a sphere of any size and composition.


title: Scanning Tunneling Microscopy of an Air Sensitive Dichalcogenide Through an Encapsulating Layer

abstract: Many atomically thin exfoliated 2D materials degrade when exposed to ambient conditions. They can be protected and investigated by means of transport and optical measurements if they are encapsulated between chemically inert single layers in the controlled atmosphere of a glove box. Here, we demonstrate that the same encapsulation procedure is also compatible with scanning tunneling microscopy (STM) and spectroscopy (STS). To this end, we report a systematic STM/STS investigation of a model system consisting of an exfoliated 2H-NbSe2 crystal capped with a protective 2H-MoS2 monolayer. We observe different electronic coupling between MoS2 and NbSe2, from a strong coupling when their lattices are aligned within a few degrees to 2 essentially no coupling for 30{\deg} misaligned layers. We show that STM always probes intrinsic NbSe2 properties such as the superconducting gap and charge density wave at low temperature when setting the tunneling bias inside the MoS2 band gap, irrespective of the relative angle between the NbSe2 and MoS2 lattices. This study demonstrates that encapsulation is fully compatible with STM/STS investigations of 2D materials.


title: Rewriting Abstract Structures: Materialization Explained Categorically

abstract: The paper develops an abstract (over-approximating) semantics for double-pushout rewriting of graphs and graph-like objects. The focus is on the so-called materialization of left-hand sides from abstract graphs, a central concept in previous work. The first contribution is an accessible, general explanation of how materializations arise from universal properties and categorical constructions, in particular partial map classifiers, in a topos. Second, we introduce an extension by enriching objects with annotations and give a precise characterization of strongest post-conditions, which are effectively computable under certain assumptions.


title: Self-Supervised Surgical Tool Segmentation using Kinematic Information

abstract: Surgical tool segmentation in endoscopic images is the first step towards pose estimation and (sub-)task automation in challenging minimally invasive surgical operations. While many approaches in the literature have shown great results using modern machine learning methods such as convolutional neural networks, the main bottleneck lies in the acquisition of a large number of manually-annotated images for efficient learning. This is especially true in surgical context, where patient-to-patient differences impede the overall generalizability. In order to cope with this lack of annotated data, we propose a self-supervised approach in a robot-assisted context. To our knowledge, the proposed approach is the first to make use of the kinematic model of the robot in order to generate training labels. The core contribution of the paper is to propose an optimization method to obtain good labels for training despite an unknown hand-eye calibration and an imprecise kinematic model. The labels can subsequently be used for fine-tuning a fully-convolutional neural network for pixel-wise classification. As a result, the tool can be segmented in the endoscopic images without needing a single manually-annotated image. Experimental results on phantom and in vivo datasets obtained using a flexible robotized endoscopy system are very promising.


title: On Nonconvex Optimization for Machine Learning: Gradients, Stochasticity, and Saddle Points

abstract: Gradient descent (GD) and stochastic gradient descent (SGD) are the workhorses of large-scale machine learning. While classical theory focused on analyzing the performance of these methods in convex optimization problems, the most notable successes in machine learning have involved nonconvex optimization, and a gap has arisen between theory and practice. Indeed, traditional analyses of GD and SGD show that both algorithms converge to stationary points efficiently. But these analyses do not take into account the possibility of converging to saddle points. More recent theory has shown that GD and SGD can avoid saddle points, but the dependence on dimension in these analyses is polynomial. For modern machine learning, where the dimension can be in the millions, such dependence would be catastrophic. We analyze perturbed versions of GD and SGD and show that they are truly efficient---their dimension dependence is only polylogarithmic. Indeed, these algorithms converge to second-order stationary points in essentially the same time as they take to converge to classical first-order stationary points.


title: Group level MEG/EEG source imaging via optimal transport: minimum Wasserstein estimates

abstract: Magnetoencephalography (MEG) and electroencephalogra-phy (EEG) are non-invasive modalities that measure the weak electromagnetic fields generated by neural activity. Inferring the location of the current sources that generated these magnetic fields is an ill-posed inverse problem known as source imaging. When considering a group study, a baseline approach consists in carrying out the estimation of these sources independently for each subject. The ill-posedness of each problem is typically addressed using sparsity promoting regularizations. A straightforward way to define a common pattern for these sources is then to average them. A more advanced alternative relies on a joint localization of sources for all subjects taken together, by enforcing some similarity across all estimated sources. An important advantage of this approach is that it consists in a single estimation in which all measurements are pooled together, making the inverse problem better posed. Such a joint estimation poses however a few challenges, notably the selection of a valid regularizer that can quantify such spatial similarities. We propose in this work a new procedure that can do so while taking into account the geometrical structure of the cortex. We call this procedure Minimum Wasserstein Estimates (MWE). The benefits of this model are twofold. First, joint inference allows to pool together the data of different brain geometries, accumulating more spatial information. Second, MWE are defined through Optimal Transport (OT) metrics which provide a tool to model spatial proximity between cortical sources of different subjects, hence not enforcing identical source location in the group. These benefits allow MWE to be more accurate than standard MEG source localization techniques. To support these claims, we perform source localization on realistic MEG simulations based on forward operators derived from MRI scans. On a visual task dataset, we demonstrate how MWE infer neural patterns similar to functional Magnetic Resonance Imaging (fMRI) maps.


title: Lower Bound Convex Programs for Exact Sparse Optimization

abstract: In exact sparse optimization problems on Rd (also known as sparsity constrained problems), one looks for solution that have few nonzero components. In this paper, we consider problems where sparsity is exactly measured either by the nonconvex l0 pseudonorm (and not by substitute penalty terms) or by the belonging of the solution to a finite union of subsets. Due to the combinatorial nature of the sparsity constraint, such problems do not generally display convexity properties, even if the criterion to minimize is convex. In the most common approach to tackle them, one replaces the sparsity constraint by a convex penalty term, supposed to induce sparsity. Thus doing, one loses the original exact sparse optimization problem, but gains convexity. However, by doing so, it is not clear that one obtains a lower bound of the original exact sparse optimization problem. In this paper, we propose another approach, where we lose convexity but where we gain at keeping the original exact sparse optimization formulation, by displaying lower bound convex minimization programs. For this purpose , we introduce suitable conjugacies, induced by a novel class of one-sided linear couplings. Thus equipped, we present a systematic way to design norms and lower bound convex minimization programs over their unit ball. The family of norms that we display encompasses most of the sparsity inducing norms used in machine learning. Therefore, our approach provides foundation and interpretation for their use.


title: Weighted Stochastic Field Exponent Sobolev Spaces and Nonlinear Degenerated Elliptic Problem

abstract: In this study, we consider weighted stochastic field exponent function spaces $L_{\vartheta }^{p(.,.)}\left( D\times \Omega \right) $ and $W_{\vartheta }^{k,p(.,.)}\left( D\times \Omega \right) $. Also, we investigate some basic properties and embeddings of these spaces. Finally, we present an application of these spaces to the stochastic partial differential equations with stochastic field growth.


title: Inexact Convex Relaxations for AC Optimal Power Flow: Towards AC Feasibility

abstract: Convex relaxations of AC optimal power flow (AC-OPF) problems have attracted significant interest as in several instances they provably yield the global optimum to the original non-convex problem. If, however, the relaxation is inexact, the obtained solution is not AC-feasible. The quality of the obtained solution is essential for several practical applications of AC-OPF, but detailed analyses are lacking in existing literature. This paper aims to cover this gap. We provide an in-depth investigation of the solution characteristics when convex relaxations are inexact, we assess the most promising AC feasibility recovery methods for large-scale systems, and we propose two new metrics that lead to a better understanding of the quality of the identified solutions. We perform a comprehensive assessment on 96 different test cases, ranging from 14 to 3120 buses, and we show the following: (i) Despite an optimality gap of less than 1%, several test cases still exhibit substantial distances to both AC feasibility and local optimality and the newly proposed metrics characterize these deviations. (ii) Penalization methods fail to recover an AC-feasible solution in 15 out of 45 cases, and using the proposed metrics, we show that most failed test instances exhibit substantial distances to both AC-feasibility and local optimality. For failed test instances with small distances, we show how our proposed metrics inform a fine-tuning of penalty weights to obtain AC-feasible solutions. (iii) The computational benefits of warm-starting non-convex solvers have significant variation, but a computational speedup exists in over 75% of the cases.


title: A Suitable Conjugacy for the l0 Pseudonorm

abstract: The so-called l0 pseudonorm on R d counts the number of nonzero components of a vector. It is well-known that the l0 pseudonorm is not convex, as its Fenchel biconjugate is zero. In this paper, we introduce a suitable conjugacy, induced by a novel coupling, Caprac, having the property of being constant along primal rays, like the l0 pseudonorm. The Caprac coupling belongs to the class of one-sided linear couplings, that we introduce. We show that they induce conjugacies that share nice properties with the classic Fenchel conjugacy. For the Caprac conjugacy, induced by the coupling Caprac, we prove that the l0 pseudonorm is equal to its biconjugate: hence, the l0 pseudonorm is Caprac-convex in the sense of generalized convexity. As a corollary, we show that the l0 pseudonorm coincides, on the sphere, with a convex lsc function. We also provide expressions for conjugates in terms of two families of dual norms, the 2-k-symmetric gauge norms and the k-support norms.


title: Truth-preservation under fuzzy pp-formulas

abstract: How can non-classical logic contribute to the analysis of complexity in computer science? In this paper, we give a step towards this question, taking a logical model-theoretic approach to the analysis of complexity in fuzzy constraint satisfaction. We study fuzzy positive-primitive sentences, and we present an algebraic characterization of classes axiomatized by this kind of sentences in terms of homomorphisms and direct products. The ultimate goal is to study the expressiveness and reasoning mechanisms of non-classical languages, with respect to constraint satisfaction problems and, in general, in modelling decision scenarios. fuzzy constraint satisfaction, fuzzy logics, model theory.


title: The Odds are Odd: A Statistical Test for Detecting Adversarial Examples

abstract: We investigate conditions under which test statistics exist that can reliably detect examples, which have been adversarially manipulated in a white-box attack. These statistics can be easily computed and calibrated by randomly corrupting inputs. They exploit certain anomalies that adversarial attacks introduce, in particular if they follow the paradigm of choosing perturbations optimally under p-norm constraints. Access to the log-odds is the only requirement to defend models. We justify our approach empirically, but also provide conditions under which detectability via the suggested test statistics is guaranteed to be effective. In our experiments, we show that it is even possible to correct test time predictions for adversarial attacks with high accuracy.


title: $L^1$-Poincar\'e inequalities for differential forms on Euclidean spaces and Heisenberg groups

abstract: In this paper, we prove interior Poincar{\'e} and Sobolev inequalities in Euclidean spaces and in Heisenberg groups, in the limiting case where the exterior (resp. Rumin) differential of a differential form is measured in L 1 norm. Unlike for L p , p > 1, the estimates are doomed to fail in top degree. The singular integral estimates are replaced with inequalities which go back to Bourgain-Brezis in Euclidean spaces, and to Chanillo-van Schaftingen in Heisenberg groups.


title: Petascale Cloud Supercomputing for Terapixel Visualization of a Digital Twin

abstract: Background: Photo-realistic terapixel visualization is computationally intensive and to date there have been no such visualizations of urban digital twins, the few terapixel visualizations that exist have looked towards space rather than earth. Objective: our aims are: creating a scalable cloud supercomputer software architecture for visualization; a photo-realistic terapixel 3D visualization of urban IoT data supporting daily updates; a rigorous evaluation of cloud supercomputing for our application. Method: we migrated the Blender Cycles path tracer to the public cloud within a new software framework designed to scale to petaFLOP performance. Results: we demonstrate we can compute a terapixel visualization in under one hour, the system scaling at 98% efficiency to use 1024 public cloud GPU nodes delivering 14 petaFLOPS. The resulting terapixel image supports interactive browsing of the city and its data at a wide range of sensing scales. Conclusion: The GPU compute resource available in the cloud is greater than anything available on our national supercomputers providing access to globally competitive resources. The direct financial cost of access, compared to procuring and running these systems, was low. The indirect cost, in overcoming teething issues with cloud software development, should reduce significantly over time.


title: From delayed minimization to the harmonic map heat equation

abstract: In the context of cell motility modelling and more particularly related to the Filament Based Lamelipodium Model [Manhart et al 2015 & 2017], this work deals with a rigorous mathematical proof of convergence between solutions of two problems : we start from a microscopic description of adhesions using a delayed and constrained vector valued equation with spacial diffusion and show the convergence towards the corresponding friction limit. The convergence is performed with respect to the bond characteristic lifetime $\varepsilon$ whose inverse is also proportional to the stifness of the bonds. The originality of this work is the extension of gradient flow techniques to our setting. Namely, the discrete finite difference term in the gradient flow energy is here replaced by a delay term which complicates greatly the mathematical analysis. Contrarily to the standard approach [Oelz SeMa 2011], compactness in time is not provided by the energy minimization process~: a series of past times are taken into account in our discrete energy. A supplementary equation on the time derivative is obtained requiring uniform estimate with respect to $\varepsilon$ of the Lagrange multiplier and provides compactness. Due to the non-linearity induced by the constraint, a specific stability estimate useful in our previous works, is not at hand here. Numerical simulations even showed that this estimate does not hold. Nevertheless, transposing our delay operator, we succeed in proving convergence under slightly weaker hypotheses. The result relies on a careful initial layer analysis, extending [Milisic Esaim Proc 2018] to the space dependent setting.


title: Existence of Weak Solutions for $p(.)$-Laplacian Equation via Compact Embeddings of the Double Weighted Variable Exponent Sobolev Spaces

abstract: In this study, we define double weighted variable exponent Sobolev spaces $W^{1,q(.),p(.)}\left( \Omega ,\vartheta _{0},\vartheta \right) $ with respect to two different weight functions. Also, we investigate the basic properties of this spaces. Moreover, we discuss the existence of weak solutions for weighted Dirichlet problem of $p(.)$-Laplacian equation \begin{equation*} \left\{ \begin{array}{cc} -\text{div}\left( \vartheta (x)\left\vert \nabla f\right\vert ^{p(x)-2}\nabla f\right) =\vartheta _{0}(x)\left\vert f\right\vert ^{q(x)-2}f & x\in \Omega \\ f=0 & x\in \partial \Omega \end{array} \right. \end{equation*} under some conditions of compact embedding involving the double weighted variable exponent Sobolev spaces.


title: X-ray pumping of the Th-229 nuclear clock isomer

abstract: Thorium-229 is a unique case in nuclear physics: it presents a metastable first excited state Th-229m, just a few electronvolts above the nuclear ground state. This so-called isomer is accessible by VUV lasers, which allows transferring the amazing precision of atomic laser spectroscopy to nuclear physics. Being able to manipulate the Th-229 nuclear states at will opens up a multitude of prospects, from studies of the fundamental interactions in physics to applications as a compact and robust nuclear clock. However, direct optical excitation of the isomer or its radiative decay back to the ground state has not yet been observed, and a series of key nuclear structure parameters such as the exact energies and half-lives of the low-lying nuclear levels of Th-229 are yet unknown. Here we present the first active optical pumping into Th-229m. Our scheme employs narrow-band 29 keV synchrotron radiation to resonantly excite the second excited state, which then predominantly decays into the isomer. We determine the resonance energy with 0.07 eV accuracy, measure a half-life of 82.2 ps, an excitation linewidth of 1.70 neV, and extract the branching ratio of the second excited state into the ground and isomeric state respectively. These measurements allow us to re-evaluate gamma spectroscopy data that have been collected over 40~years.


title: Diagnosing the Clumpy Protoplanetary Disk of the UXor Type Young Star GM Cephei

abstract: UX Orionis stars (UXors) are Herbig Ae/Be or T Tauri stars exhibiting sporadic occultation of stellar light by circumstellar dust. GM\,Cephei is such a UXor in the young ($\sim4$~Myr) open cluster Trumpler\,37, showing prominent infrared excess, emission-line spectra, and flare activity. Our photometric monitoring (2008--2018) detects (1)~an $\sim$3.43~day period, likely arising from rotational modulation by surface starspots, (2)~sporadic brightening on time scales of days due to accretion, (3)~irregular minor flux drops due to circumstellar dust extinction, and (4)~major flux drops, each lasting for a couple of months with a recurrence time, though not exactly periodic, of about two years. The star experiences normal reddening by large grains, i.e., redder when dimmer, but exhibits an unusual "blueing" phenomenon in that the star turns blue near brightness minima. The maximum extinction during relatively short (lasting $\leq 50$~days) events, is proportional to the duration, a consequence of varying clump sizes. For longer events, the extinction is independent of duration, suggestive of a transverse string distribution of clumps. Polarization monitoring indicates an optical polarization varying $\sim3\%$--8$\%$, with the level anticorrelated with the slow brightness change. Temporal variation of the unpolarized and polarized light sets constraints on the size and orbital distance of the circumstellar clumps in the interplay with the young star and scattering envelope. These transiting clumps are edge-on manifestations of the ring- or spiral-like structures found recently in young stars with imaging in infrared of scattered light, or in submillimeter of thermalized dust emission.


title: Build-up of Vibron-Mediated Electron Correlations in Molecular Junctions

abstract: We investigate on the same footing the time-dependent electronic transport properties and vibrational dynamics of a molecular junction. We show that fluctuations of both the molecular vibron displacement and the electronic current across the junction undergo damped oscillations towards the steady-state. We assign the former to the onset of electron tunneling events assisted by vibron-emission. The time-dependent build-up of electron-hole correlations is revealed as a departure of the charge-transfer statistics from the generalized-binomial one after a critical time tc. The phonon-back action on the tunneling electrons is shown to amplify and accelerate this build-up mechanism.


title: Reassessment of the basis of cell size control based on analysis of cell-to-cell variability

abstract: Fundamental mechanisms governing cell size control and homeostasis are still poorly understood. The relationship between sizes at division and birth in single cells is used as a metric to categorize the basis of size homeostasis [1-3]. Cells dividing at a fixed size regardless of birth size (sizer) are expected to show a division-birth slope of 0, whereas cells dividing after growing for a fixed size increment (adder) have an expected slope of +1 [4]. These two theoretical values are, however, rarely experimentally observed. For example, rod-shaped fission yeast $\it{Schizosaccharomyces}$ $\it{pombe}$ cells, which divide at a fixed surface area [5, 6], exhibit a division-birth slope for cell lengths of 0.25$\pm$0.02, significantly different from the expected sizer value of zero. Here we investigate possible reasons for this discrepancy by developing a mathematical model of sizer control including the relevant sources of variation. Our results support $\it{pure}$ sizer control and show that deviation from zero slope is exaggerated by measurement of an inappropriate geometrical quantity (e.g., length instead of area), combined with cell-to-cell radius variability. The model predicts that mutants with greater errors in size sensing or septum positioning paradoxically appear to behave as better sizers. Furthermore, accounting for cell width variability, we show that pure sizer control can in some circumstances reproduce the apparent adder behaviour observed in $\it{E. coli}$. These findings demonstrate that analysis of geometric variation can lead to new insights into cell size control.


title: Bayesian inference using synthetic likelihood: asymptotics and adjustments

abstract: Implementing Bayesian inference is often computationally challenging in applications involving complex models, and sometimes calculating the likelihood itself is difficult. Synthetic likelihood is one approach for carrying out inference when the likelihood is intractable, but it is straightforward to simulate from the model. The method constructs an approximate likelihood by taking a vector summary statistic as being multivariate normal, with the unknown mean and covariance matrix estimated by simulation for any given parameter value. Our article makes three contributions. The first shows that if the summary statistic satisfies a central limit theorem, then the synthetic likelihood posterior is asymptotically normal and yields credible sets with the correct level of frequentist coverage. This result is similar to that obtained by approximate Bayesian computation. The second contribution compares the computational efficiency of Bayesian synthetic likelihood and approximate Bayesian computation using the acceptance probability for rejection and importance sampling algorithms with a "good" proposal distribution. We show that Bayesian synthetic likelihood is computationally more efficient than approximate Bayesian computation, and behaves similarly to regression-adjusted approximate Bayesian computation. Based on the asymptotic results, the third contribution proposes using adjusted inference methods when a possibly misspecified form is assumed for the covariance matrix of the synthetic likelihood, such as diagonal or a factor model, to speed up the computation. The methodology is illustrated with some simulated and real examples.


title: On the achromatic number of signed graphs

abstract: In this paper, we generalize the concept of complete coloring and achromatic number to 2-edge-colored graphs and signed graphs. We give some useful relationships between different possible definitions of such achromatic numbers and prove that computing any of them is NP-complete.


title: Mass-conserving self-similar solutions to coagulation-fragmentation equations

abstract: Existence of mass-conserving self-similar solutions with a sufficiently small total mass is proved for a specific class of homogeneous coagulation and fragmentation coefficients. The proof combines a dynamical approach to construct such solutions for a regularised coagulation-fragmentation equation in scaling variables and a compactness method.


title: Volume forms on moduli spaces of d-differentials

abstract: Given $d\in \mathbb{N}$, $g\in \mathbb{N} \cup\{0\}$, and an integral vector $\kappa=(k_1,\dots,k_n)$ such that $k_i>-d$ and $k_1+\dots+k_n=d(2g-2)$, let $\Omega^d\mathcal{M}_{g,n}(\kappa)$ denote the moduli space of meromorphic $d$-differentials on Riemann surfaces of genus $g$ whose zeros and poles have orders prescribed by $\kappa$. We show that $\Omega^d\mathcal{M}_{g,n}(\kappa)$ carries a canonical volume form that is parallel with respect to its affine complex manifold structure, and that the total volume of $\mathbb{P}\Omega^d\mathcal{M}_{g,n}(\kappa)=\Omega^d\mathcal{M}_{g,n}/\mathbb{C}^*$ with respect to the measure induced by this volume form is finite.


title: On the generalisation of cohomology with compact support to non-finite type schemes

abstract: In this article we extend Deligne's construction of Grothendieck's six operations on the derived category of torsion sheaves over the \'etale site of a scheme for morphisms of finite type to a larger class of morphisms. This class includes profinite \'etale coverings as well as separated morphisms perfectly of finite type


title: Relative rationality: Is machine rationality subjective?

abstract: Rational decision making in its linguistic description means making logical decisions. In essence, a rational agent optimally processes all relevant information to achieve its goal. Rationality has two elements and these are the use of relevant information and the efficient processing of such information. In reality, relevant information is incomplete, imperfect and the processing engine, which is a brain for humans, is suboptimal. Humans are risk averse rather than utility maximizers. In the real world, problems are predominantly non-convex and this makes the idea of rational decision-making fundamentally unachievable and Herbert Simon called this bounded rationality. There is a trade-off between the amount of information used for decision-making and the complexity of the decision model used. This explores whether machine rationality is subjective and concludes that indeed it is.


title: Delayed Outburst of H 1743--322 in 2003 and relation with its other outbursts

abstract: The Galactic transient black hole candidate H 1743--322 exhibited a long duration outburst in 2003 after more than two and a half decades of inactivity. The 2003 event was extensively studied in multi-wavelength bands by many groups. The striking feature is that the total energy released is extremely high as compared to that in tens of outbursts which followed. In this paper, we look at this event and study both the spectral and temporal properties of the source using two component advective flow (TCAF) paradigm. We extract accretion flow parameters for each observation from spectral properties of the decay phase and determine the mass of the black hole. We computed the energy released during all the known outbursts since 2003 and showed that on an average, the energy release in an outburst is proportional to the duration of the quiescent state just prior to it, with the exception of the 2004 outburst. A constant rate of supply of matter from the companion cannot explain the energy release in 2004 outburst. However, if the energy release of 2003 is incomplete and the leftover is released in 2004, then the companion's rate of matter supply can be constant since 1977 till date. We believe that erratic behaviour of viscosity at the accumulation radius $X_p$ of matter as well as location the $X_p$ itself, rather than the random variation of mass transfer rate from the companion, could be responsible for non-uniformity in outburst pattern. We discuss several factors on which the waiting time and duration of the next outburst could depend.


title: Dynamic Non-Diagonal Regularization in Interior Point Methods for Linear and Convex Quadratic Programming

abstract: In this paper, we present a dynamic non-diagonal regularization for interior point methods. The non-diagonal aspect of this regularization is implicit, since all the off-diagonal elements of the regularization matrices are cancelled out by those elements present in the Newton system, which do not contribute important information in the computation of the Newton direction. Such a regularization has multiple goals. The obvious one is to improve the spectral properties of the Newton system solved at each iteration of the interior point method. On the other hand, the regularization matrices introduce sparsity to the aforementioned linear system, allowing for more efficient factorizations. We also propose a rule for tuning the regularization dynamically based on the properties of the problem, such that sufficiently large eigenvalues of the non-regularized system are perturbed insignificantly. This alleviates the need of finding specific regularization values through experimentation, which is the most common approach in literature. We provide perturbation bounds for the eigenvalues of the non-regularized system matrix and then discuss the spectral properties of the regularized matrix. Finally, we demonstrate the efficiency of the method applied to solve standard small and medium-scale linear and convex quadratic programming test problems.


title: Two-dimensional pattern formation in ionic liquids confined between graphene walls

abstract: We perform molecular dynamics simulations of ionic liquids confined between graphene walls under a large variety of conditions (pure ionic liquids, mixtures with water and alcohols, mixtures with lithium salts and defective graphene walls). Our results show that the formation of striped and hexagonal patterns in the Stern layer can be considered as a general feature of ionic liquids at electrochemical interfaces, the transition between patterns being controlled by the net balance of charge in the innermost layer of adsorbed molecules. This explains previously reported experimental and computational results and, for the first time, why these pattern changes are triggered by any perturbation of the charge density at the innermost layer of the electric double layer (voltage and composition changes, and vacancies at the electrode walls, among others), which may help tuning electrode-ionic liquid interfaces. Using Monte Carlo simulations we show that such structures can be reproduced by a simple two-dimensional lattice model with only nearest-neighbour interactions, governed by highly screened ionic interactions and short-range and excluded volume interactions. We also show that the results of our simulations are consistent with those inferred from the Landau-Brazovskii theory of pattern formation in self-assembling systems. The presence of these patterns at the ionic liquid graphene-electrode interfaces may have a strong impact on the process of ionic transfer from the bulk mixtures to the electrodes, on the differential capacitance of the electrode-electrolyte double layer or on the rates of redox reactions at the electrodes, among other physicochemical properties, and is therefore an effect of great technological interest.


title: Differentials and distances in probabilistic coherence spaces

abstract: In probabilistic coherence spaces, a denotational model of probabilistic functional languages, mor-phisms are analytic and therefore smooth. We explore two related applications of the corresponding derivatives. First we show how derivatives allow to compute the expectation of execution time in the weak head reduction of probabilistic PCF (pPCF). Next we apply a general notion of "local" differential of morphisms to the proof of a Lipschitz property of these morphisms allowing in turn to relate the observational distance on pPCF terms to a distance the model is naturally equipped with. This suggests that extending probabilistic programming languages with derivatives, in the spirit of the differential lambda-calculus, could be quite meaningful.


title: Waves Interacting With A Partially Immersed Obstacle In The Boussinesq Regime

abstract: This paper is devoted to the derivation and mathematical analysis of a wave-structure interaction problem which can be reduced to a transmission problem for a Boussinesq system. Initial boundary value problems and transmission problems in dimension d= 1 for $2\times2$ hyperbolic systems are well understood. However, for many applications, and especially for the description of surface water waves, dispersive perturbations of hyperbolic systems must be considered. We consider here a configuration where the motion of the waves is governed by a Boussinesq system (a dispersive perturbation of the hyperbolic nonlinear shallow water equations), and in the presence of a fixed partially immersed obstacle. We shall insist on the differences and similarities with respect to the standard hyperbolic case, and focus our attention on a new phenomenon, namely, the apparition of a dispersive boundary layer. In order to obtain existence and uniform bounds on the solutions over the relevant time scale, a control of this dispersive boundary layer and of the oscillations in time it generates is necessary. This analysis leads to a new notion of compatibility condition that is shown to coincide with the standard hyperbolic compatibility conditions when the dispersive parameter is set to zero. To the authors' knowledge, this is the first time that these phenomena (likely to play a central role in the analysis of initial boundary value problems for dispersive perturbations of hyperbolic systems) are exhibited.


title: A Partial Data Problem in Linear Elasticity

abstract: We discuss the determination of the Lam\'e parameters of an elastic material by the means of boundary measurements. We will combine previous results of Eskin-Ralston and Isakov to prove inverse results in the case of bounded domains with partial data. Moreover, we generalise these results to infinite cylinders.


title: Mv-algebras And Partially Cyclically Ordered Groups

abstract: We prove that there exists a functorial correspondence between MV-algebras and partially cyclically ordered groups which are wound round of lattice-ordered groups. It follows that some results about cyclically ordered groups can be stated in terms of MV-algebras. For example, the study of groups together with a cyclic order allows to get a first-order characterization of groups of unimodular complex numbers and of finite cyclic groups. We deduce a characterization of pseudofinite MV-chains and of pseudo-simple MV-chains (i.e. which share the same first-order properties as some simple ones). We can generalize these results to some non-lineraly ordered MV-algebras, for example hyper-archimedean MV-algebras.