title: 2-D model of the global ionospheric conductor connected with the magnetospheric conductors

abstract: A model of the ionospheric global conductor is designed. The ionospheric conductor is considered in the framework of a two-dimensional approximation based on high conductivity in the direction of the magnetic field. Under this assumption the magnetic field lines are equipotential, and the charge transfer between them is determined only by integral Pedersen and Hall conductivities. The model is constructed as the first approximation in the small parameter expansion of the solution of the three-dimensional problems of electrical conductivity. The small parameter is the ratio of Pedersen and field-aligned conductivities. The space distributions of the Pedersen and Hall conductivities are calculated using the empirical models IRI, MSISE, IGRF and applied to construct the maps of the integral conductivities. The parts of the magnetosphere with high conductivity across the magnetic field lines, namely, the cusps and the plasma layer are analyzed. It is shown that the connection of these magnetospheric conductors to the ionosphere in parallel makes the auroral zones equipotential regions. As a consequence, for the ionospheric electric fields, which generators are located in the ionosphere or in the atmosphere, the global problem of electrical conductivity is separated into three independent boundary value problems in three regions: two polar caps and the main part of the ionosphere which includes the mid- and low-latitude parts of the ionosphere. The model can be used for the analysis of the ionospheric part of the Global Electric Circuit, for calculation of the ionospheric dynamo electric field and as a fragment in more complex ionospheric and magnetospheric models.


title: Stereo obstacle detection for unmanned surface vehicles by IMU-assisted semantic segmentation

abstract: A new obstacle detection algorithm for unmanned surface vehicles (USVs) is presented. A state-of-the-art graphical model for semantic segmentation is extended to incorporate boat pitch and roll measurements from the on-board inertial measurement unit (IMU), and a stereo verification algorithm that consolidates tentative detections obtained from the segmentation is proposed. The IMU readings are used to estimate the location of horizon line in the image, which automatically adjusts the priors in the probabilistic semantic segmentation model. We derive the equations for projecting the horizon into images, propose an efficient optimization algorithm for the extended graphical model, and offer a practical IMU-camera-USV calibration procedure. Using an USV equipped with multiple synchronized sensors, we captured a new challenging multi-modal dataset, and annotated its images with water edge and obstacles. Experimental results show that the proposed algorithm significantly outperforms the state of the art, with nearly 30% improvement in water-edge detection accuracy, an over 21% reduction of false positive rate, an almost 60% reduction of false negative rate, and an over 65% increase of true positive rate, while its Matlab implementation runs in real-time.


title: Non-rigid Object Tracking via Deep Multi-scale Spatial-temporal Discriminative Saliency Maps

abstract: In this paper, we propose a novel effective non-rigid object tracking framework based on the spatial-temporal consistent saliency detection. In contrast to most existing trackers that utilize a bounding box to specify the tracked target, the proposed framework can extract accurate regions of the target as tracking outputs. It achieves a better description of the non-rigid objects and reduces the background pollution for the tracking model. Furthermore, our model has several unique features. First, a tailored fully convolutional neural network (TFCN) is developed to model the local saliency prior for a given image region, which not only provides the pixel-wise outputs but also integrates the semantic information. Second, a novel multi-scale multi-region mechanism is proposed to generate local saliency maps that effectively consider visual perceptions with different spatial layouts and scale variations. Subsequently, local saliency maps are fused via a weighted entropy method, resulting in a final discriminative saliency map. Finally, we present a non-rigid object tracking algorithm based on the predicted saliency maps. By utilizing a spatial-temporal consistent saliency map (STCSM), we conduct target-background classification and use a simple fine-tuning scheme for online updating. Extensive experiments demonstrate that the proposed algorithm achieves competitive performance in both saliency detection and visual tracking, especially outperforming other related trackers on the non-rigid object tracking datasets.


title: A New Approach for Measuring Power Spectra and Reconstructing Time Series in Active Galactic Nuclei

abstract: We provide a new approach to measure power spectra and reconstruct time series in active galactic nuclei (AGNs) based on the fact that the Fourier transform of AGN stochastic variations is a series of complex Gaussian random variables. The approach parameterizes a stochastic series in frequency domain and transforms it back to time domain to fit the observed data. The parameters and their uncertainties are derived in a Bayesian framework, which also allows us to compare the relative merits of different power spectral density models. The well-developed fast Fourier transform algorithm together with parallel computation enable an acceptable time complexity for the approach.


title: Power law analysis for temperature dependence of magnetocrystalline anisotropy constants of Nd$_2$Fe$_{14}$B magnets

abstract: Phenomenological analysis for the temperature dependence of the magnetocrystalline anisotropy (MA) in rare earth magnets is presented. We define phenomenological power laws applicable to compound magnets using the Zener theory, apply these laws to the magnetocrystalline anisotropy constants (MACs) of Nd$_2$Fe$_{14}$B magnets. The results indicate that the MACs obey the power law well, and a general understanding for the temperature-dependent MA in rare earth magnets is obtained through the analysis. Furthermore, to examine the validity of the power law, we discuss the temperature dependence of the MACs in Dy$_2$Fe$_{14}$B and Y$_2$Fe$_{14}$B magnets as examples wherein it is difficult to interpret the MA using the power law.


title: Electric quadrupole and magnetic dipole coupling in plasmonic nanoparticle arrays

abstract: Collective resonances in plasmonic nanoparticle arrays with electric dipole moment oriented along the lattice wave propagation are theoretically investigated. The role of electric quadrupole (EQ) and magnetic dipole (MD) moments of gold nanoparticles in the resonant features of the arrays is analyzed. We perform both semi-analytical calculations of coupled multipole equations and rigorous numerical simulations varying contributions of the electric and magnetic multipoles by changing particle size and shape (spheres and disks). The arrays in homogeneous and non-homogeneous environments are considered. We find that even very weak non-resonant EQ and MD moments of a single particle are significantly enhanced in the periodic lattice at the wavelength of collective (lattice) resonance excitation. Importantly, we show that in the infinite arrays, the EQ and MD moments of nanoparticles are coupled and affect each other resonant contributions. We also demonstrate that at the lattice-resonance wavelength, the enhanced EQ and MD moments have contributions to reflection comparable to the dipole one resulting in a significant decrease of reflection and providing the satisfaction of the generalized Kerker condition for reflection suppression.


title: Well-posedness to the continuous coagulation processes with collision-induced multiple fragmentation

abstract: An existence result on weak solutions to the continuous coagulation equation with collision-induced multiple fragmentation is established for certain classes of unbounded coagulation, collision and breakup kernels. In this model, a pair of particles can coagulate into a larger one if their confrontation is a completely inelastic collision; otherwise, one of them will split into many smaller particles due to a destructive collision. In the present work, both coagulation and fragmentation processes are considered to be intrinsically nonlinear. The breakup kernel may have a possibility to attain a singularity at the origin. The proof is based on the classical weak L^1 compactness method applied to suitably chosen approximating equations. In addition, we study the uniqueness of weak solutions under additional growth conditions on collision and breakup kernels which mainly relies on the integrability of higher moments. Finally, it is obtained that the unique weak solution is mass-conserving.


title: Entangled systems are unbounded sources of nonlocal correlations and of certified random numbers

abstract: The outcomes of local measurements made on entangled systems can be certified to be random provided that the generated statistics violate a Bell inequality. This way of producing randomness relies only on a minimal set of assumptions because it is independent of the internal functioning of the devices generating the random outcomes. In this context it is crucial to understand both qualitatively and quantitatively how the three fundamental quantities -- entanglement, non-locality and randomness -- relate to each other. To explore these relationships, we consider the case where repeated (non projective) measurements are made on the physical systems, each measurement being made on the post-measurement state of the previous measurement. In this work, we focus on the following questions: For systems in a given entangled state, how many nonlocal correlations in a sequence can we obtain by measuring them repeatedly? And from this generated sequence of non-local correlations, how many random numbers is it possible to certify? In the standard scenario with a single measurement in the sequence, it is possible to generate non-local correlations between two distant observers only and the amount of random numbers is very limited. Here we show that we can overcome these limitations and obtain any amount of certified random numbers from an entangled pair of qubit in a pure state by making sequences of measurements on it. Moreover, the state can be arbitrarily weakly entangled. In addition, this certification is achieved by near-maximal violation of a particular Bell inequality for each measurement in the sequence. We also present numerical results giving insight on the resistance to imperfections and on the importance of the strength of the measurements in our scheme.


title: Semi-explicit solutions to the water-wave dispersion relation and their role in the nonlinear Hamiltonian coupled-mode theory

abstract: The Hamiltonian coupled-mode theory (HCMT), recently derived by Athanassoulis and Papoutsellis [1], provides an efficient new approach for solving fully nonlinear water-wave problems over arbitrary bathymetry. In HCMT, heavy use is made of the roots of a local, water-wave dispersion relation with varying parameter, which have to be calculated at every horizontal position and every time instant. Thus, fast and accurate calculation of these roots, valid for all possible values of the varying parameter, are of fundamental importance. In this paper, new, semi-explicit and highly accurate root-finding formulae are derived, especially for the roots corresponding to evanescent modes. The derivation is based on the successive application of a Picard-type iteration and the Householder's root finding method. Explicit approximate formulae of very good accuracy are obtained, which are adequate to support HCMT for many types of applications. In most demanding cases, e.g. very steep, deep-water waves, machine-accurate determination of the required roots is achieved by no more than three iterations, using the explicit forms as initial values. Exploiting this root-finding procedure in the HCMT, results in an efficient, numerical solver able to treat fully nonlinear water waves over arbitrary bathymetry. Applications to demanding nonlinear problems demonstrate the efficiency and the robustness of the present approach.


title: Topological phases of non-Hermitian systems

abstract: Recent experimental advances in controlling dissipation have brought about unprecedented flexibility in engineering non-Hermitian Hamiltonians in open classical and quantum systems. A particular interest centers on the topological properties of non-Hermitian systems, which exhibit unique phases with no Hermitian counterparts. However, no systematic understanding in analogy with the periodic table of topological insulators and superconductors has been achieved. In this paper, we develop a coherent framework of topological phases of non-Hermitian systems. After elucidating the physical meaning and the mathematical definition of non-Hermitian topological phases, we start with one-dimensional lattices, which exhibit topological phases with no Hermitian counterparts and are found to be characterized by an integer topological winding number even with no symmetry constraint, reminiscent of the quantum Hall insulator in Hermitian systems. A system with a nonzero winding number, which is experimentally measurable from the wave-packet dynamics, is shown to be robust against disorder, a phenomenon observed in the Hatano-Nelson model with asymmetric hopping amplitudes. We also unveil a novel bulk-edge correspondence that features an infinite number of (quasi-)edge modes. We then apply the K-theory to systematically classify all the non-Hermitian topological phases in the Altland-Zirnbauer classes in all dimensions. The obtained periodic table unifies time-reversal and particle-hole symmetries, leading to highly nontrivial predictions such as the absence of non-Hermitian topological phases in two dimensions. We provide concrete examples for all the nontrivial non-Hermitian AZ classes in zero and one dimensions. In particular, we identify a Z2 topological index for arbitrary quantum channels. Our work lays the cornerstone for a unified understanding of the role of topology in non-Hermitian systems.

categories: cond-mat.mes-hall cond-mat.dis-nn cond-mat.quant-gas cond-mat.stat-mech quant-ph


title: An Observationally-Constrained Model of a Flux Rope that Formed in the Solar Corona

abstract: Coronal mass ejections (CMEs) are large-scale eruptions of plasma from the coronae of stars. Understanding the plasma processes involved in CME initiation has applications to space weather forecasting and laboratory plasma experiments. James et al. (Sol. Phys. 292, 71, 2017) used EUV observations to conclude that a magnetic flux rope formed in the solar corona above NOAA Active Region 11504 before it erupted on 14 June 2012 (SOL2012-06-14). In this work, we use data from the Solar Dynamics Observatory to model the coronal magnetic field of the active region one hour prior to eruption using a nonlinear force-free field extrapolation, and find a flux rope reaching a maximum height of 150 Mm above the photosphere. Estimations of the average twist of the strongly asymmetric extrapolated flux rope are between 1.35 and 1.88 turns, depending on the choice of axis, although the erupting structure was not observed to kink. The decay index near the apex of the axis of the extrapolated flux rope is comparable to typical critical values required for the onset of the torus instability, so we suggest that the torus instability drove the eruption.


title: Incremental and Iterative Learning of Answer Set Programs from Mutually Distinct Examples

abstract: Over the years the Artificial Intelligence (AI) community has produced several datasets which have given the machine learning algorithms the opportunity to learn various skills across various domains. However, a subclass of these machine learning algorithms that aimed at learning logic programs, namely the Inductive Logic Programming algorithms, have often failed at the task due to the vastness of these datasets. This has impacted the usability of knowledge representation and reasoning techniques in the development of AI systems. In this research, we try to address this scalability issue for the algorithms that learn answer set programs. We present a sound and complete algorithm which takes the input in a slightly different manner and performs an efficient and more user controlled search for a solution. We show via experiments that our algorithm can learn from two popular datasets from machine learning community, namely bAbl (a question answering dataset) and MNIST (a dataset for handwritten digit recognition), which to the best of our knowledge was not previously possible. The system is publicly available at https://goo.gl/KdWAcV. This paper is under consideration for acceptance in TPLP.


title: Near Isometric Terminal Embeddings for Doubling Metrics

abstract: Given a metric space $(X,d)$, a set of terminals $K\subseteq X$, and a parameter $t\ge 1$, we consider metric structures (e.g., spanners, distance oracles, embedding into normed spaces) that preserve distances for all pairs in $K\times X$ up to a factor of $t$, and have small size (e.g. number of edges for spanners, dimension for embeddings). While such terminal (aka source-wise) metric structures are known to exist in several settings, no terminal spanner or embedding with distortion close to 1, i.e., $t=1+\epsilon$ for some small $0<\epsilon<1$, is currently known. Here we devise such terminal metric structures for {\em doubling} metrics, and show that essentially any metric structure with distortion $1+\epsilon$ and size $s(|X|)$ has its terminal counterpart, with distortion $1+O(\epsilon)$ and size $s(|K|)+1$. In particular, for any doubling metric on $n$ points, a set of $k=o(n)$ terminals, and constant $0<\epsilon<1$, there exists: (1) A spanner with stretch $1+\epsilon$ for pairs in $K\times X$, with $n+o(n)$ edges. (2) A labeling scheme with stretch $1+\epsilon$ for pairs in $K\times X$, with label size $\approx \log k$. (3) An embedding into $\ell_\infty^d$ with distortion $1+\epsilon$ for pairs in $K\times X$, where $d=O(\log k)$. Moreover, surprisingly, the last two results apply if only $K$ is a doubling metric, while $X$ can be arbitrary.


title: High-pressure phase diagram of hydrogen and deuterium sulfides from first principles: structural and vibrational properties including quantum and anharmonic effects

abstract: We study the structural and vibrational properties of the high-temperature superconducting sulfur trihydride and trideuteride in the high-pressure $Im\bar{3}m$ and $R3m$ phases by first-principles density-functional-theory calculations. On lowering pressure, the rhombohedral transition $Im\bar{3}m \rightarrow R3m$ is expected, with hydrogen bond desymmetrization and occurrence of trigonal lattice distortion. In hydrostatic conditions we find that, contrary to what suggested in some recent experiments, if the rhombohedral distortion exists it affects mainly the hydrogen-bonds, whereas the resulting cell distortion is minimal. We estimate that the occurrence of a stress anisotropy of approximately $10\%$ could explain this discrepancy. Assuming hydrostatic conditions, we calculate the critical pressure at which the rhombohedral transition occurs. Quantum and anharmonic effects, which are relevant in this system, are included at nonperturbative level with the stochastic self-consistent harmonic approximation (SSCHA). Within this approach, we determine the transition pressure by calculating the free energy Hessian. We find that quantum anharmonic effects are responsible for a strong reduction of the critical pressure with respect to the one obtained with the classical harmonic approach. Moreover, we observe a prominent isotope effect, as we estimate higher pressure transition for D${}_3$S than for H${}_3$S. Finally, within SSCHA we calculate the anharmonic phonon spectral functions in the $Im\bar{3}m$ phase. The strong anharmonicity of the system is confirmed by the occurrence of very large anharmonic broadenings leading to complex non-Lorentzian line shapes. However, for the vibrational spectra at zone center, accessible e.g. by infrared spectroscopy, the broadenings are very small (linewidth at most around 2~meV) and anharmonic phonon quasiparticles are well defined.


title: The continuous coagulation and nonlinear multiple fragmentation equation

abstract: The present paper deals with the existence and uniqueness of global classical solutions to the continuous coagulation and nonlinear multiple fragmentation equations for large classes of unbounded coagulation, collision and breakup kernels. In addition, it is shown that solutions are mass conserving. The coagulation and breakup kernels may have singularities on both the co-ordinate axes whereas the collision kernel grows up to bilinearity.


title: The exterior derivative of the Lee form of almost Hermitian manifolds

abstract: The exterior derivative $d \theta$ of the Lee form $\theta$ of almost Hermitian manifolds is studied. If $\omega$ is the K\"ahler two-form, it is proved that the $\mathbb{R}\omega$-component of $d\theta$ is always zero. expressions for the other components, in $[\lambda_0^{1,1}]$ and in $[[ \lambda^{2,0} ]]$, of $d\theta$ are also obtained. They are given in terms of the intrinsic torsion. Likewise, it is described some interrelations between the Lee form and $U(n)$-components of the Riemannian curvature tensor.


title: Robustness of classifiers to uniform $\ell\_p$ and Gaussian noise

abstract: We study the robustness of classifiers to various kinds of random noise models. In particular, we consider noise drawn uniformly from the $\ell\_p$ ball for $p \in [1, \infty]$ and Gaussian noise with an arbitrary covariance matrix. We characterize this robustness to random noise in terms of the distance to the decision boundary of the classifier. This analysis applies to linear classifiers as well as classifiers with locally approximately flat decision boundaries, a condition which is satisfied by state-of-the-art deep neural networks. The predicted robustness is verified experimentally.


title: Thermal Model Description of p--Pb Collisions at $\sqrt{s_{NN}} $ = 5.02 TeV

abstract: The ALICE data on light flavor hadron production obtained in $p-Pb$ collisions at $\sqrt{s_{NN}} $ = 5.02 TeV are studied in the thermal model using the canonical approach with exact strangeness conservation. The chemical freeze-out temperature is independent of centrality except for the lowest multiplicity bin, with values close to 160 MeV but consistent with those obtained in $Pb-Pb$ collisions at $\sqrt{s_{NN}}$ = 2.76 TeV. The value of the strangeness non-equilibrium factor $\gamma_s$ is slowly increasing with multiplicity from 0.9 to 0.96, i.e. it is always very close to full chemical equilibrium.


title: On higher dimensional singularities for the fractional Yamabe problem: a non-local Mazzeo-Pacard program

abstract: We consider the problem of constructing solutions to the fractional Yamabe problem that are singular at a given smooth sub-manifold, and we establish the classical gluing method of Mazzeo and Pacard for the scalar curvature in the fractional setting. This proof is based on the analysis of the model linearized operator, which amounts to the study of an ODE, and thus our main contribution here is the development of new methods coming from conformal geometry and scattering theory for the study of non-local ODEs. No traditional phase-plane analysis is available here. Instead, first, we provide a rigorous construction of radial fast-decaying solutions by a blow-up argument and a bifurcation method. Second, we use conformal geometry to rewrite this non-local ODE, giving a hint of what a non-local phase-plane analysis should be. Third, for the linear theory, we examine a fractional Schr\"{o}dinger equation with a Hardy type critical potential. We construct its Green's function, deduce Fredholm properties, and analyze its asymptotics at the singular points in the spirit of Frobenius method. Surprisingly enough, a fractional linear ODE may still have a two-dimensional kernel as in the second order case.


title: Evolution in complex objects

abstract: This paper describes work carried out on a model for the evolution of graph classes in complex objects. By defining evolution rules and propagation strategies on graph classes, we aim to define a user-definable means to manage data evolution model which tackles the complex nature of the classes managed, using the concepts defined in object systems. So, depending on their needs and on those of the targeted application, designers can choose the evolution mechanism they consider to suit them best. They can either create new evolutions or reuse predefined ones to respond to a given need.


title: Options for encoding names for data linking at the Australian Bureau of Statistics

abstract: Publicly, ABS has said it would use a cryptographic hash function to convert names collected in the 2016 Census of Population and Housing into an unrecognisable value in a way that is not reversible. In 2016, the ABS engaged the University of Melbourne to provide expert advice on cryptographic hash functions to meet this objective. For complex unit-record level data, including Census data, auxiliary data can be often be used to link individual records, even without names. This is the basis of ABS's existing bronze linking. This means that records can probably be re-identified without the encoded name anyway. Protection against re-identification depends on good processes within ABS. The undertaking on the encoding of names should therefore be considered in the full context of auxiliary data and ABS processes. There are several reasonable interpretations: 1. That the encoding cannot be reversed except with a secret key held by ABS. This is the property achieved by encryption (Option 1), if properly implemented; 2. That the encoding, taken alone without auxiliary data, cannot be reversed to a single value. This is the property achieved by lossy encoding (Option 2), if properly implemented; 3. That the encoding doesn't make re-identification easier, or increase the number of records that can be re-identified, except with a secret key held by ABS. This is the property achieved by HMAC-based linkage key derivation using subsets of attributes (Option 3), if properly implemented. We explain and compare the privacy and accuracy guarantees of five possible approaches. Options 4 and 5 investigate more sophisticated options for future data linking. We also explain how some commonly-advocated techniques can be reversed, and hence should not be used.


title: Algorithmic Boundedness-From-Below Conditions for Generic Scalar Potentials

abstract: Checking that a scalar potential is bounded from below (BFB) is an ubiquitous and notoriously difficult task in many models with extended scalar sectors. Exact analytic BFB conditions are known only in simple cases. In this work, we present a novel approach to algorithmically establish the BFB conditions for any polynomial scalar potential. The method relies on elements of multivariate algebra, in particular, on resultants and on the spectral theory of tensors, which is being developed by the mathematical community. We give first a pedagogical introduction to this approach, illustrate it with elementary examples, and then present the working Mathematica implementation publicly available at GitHub. Due to the rapidly increasing complexity of the problem, we have not yet produced ready-to-use analytical BFB conditions for new multi-scalar cases. But we are confident that the present implementation can be dramatically improved and may eventually lead to such results.


title: AFLOW-SYM: Platform for the complete, automatic and self-consistent symmetry analysis of crystals

abstract: Determination of the symmetry profile of structures is a persistent challenge in materials science. Results often vary amongst standard packages, hindering autonomous materials development by requiring continuous user attention and educated guesses. Here, we present a robust procedure for evaluating the complete suite of symmetry properties, featuring various representations for the point-, factor-, space groups, site symmetries, and Wyckoff positions. The protocol determines a system-specific mapping tolerance that yields symmetry operations entirely commensurate with fundamental crystallographic principles. The self consistent tolerance characterizes the effective spatial resolution of the reported atomic positions. The approach is compared with the most used programs and is successfully validated against the space group information provided for over 54,000 entries in the Inorganic Crystal Structure Database. Subsequently, a complete symmetry analysis is applied to all 1.7$+$ million entries of the AFLOW data repository. The AFLOW-SYM package has been implemented in, and made available for, public use through the automated, $\textit{ab-initio}$ framework AFLOW.


title: Cambrian acyclic domains: counting $c$-singletons

abstract: We study the size of certain acyclic domains that arise from geometric and combinatorial constructions. These acyclic domains consist of all permutations visited by commuting equivalence classes of maximal reduced decompositions if we consider the symmetric group and, more generally, of all c-singletons of a Cambrian lattice associated to the weak order of a finite Coxeter group. For this reason, we call these sets Cambrian acyclic domains. Extending a closed formula of Galambos--Reiner for a particular acyclic domain called Fishburn's alternating scheme, we provide explicit formulae for the size of any Cambrian acyclic domain and characterize the Cambrian acyclic domains of minimum or maximum size.


title: Nautilus multi-grain model: Importance of cosmic-ray-induced desorption in determining the chemical abundances in the ISM

abstract: Species abundances in the interstellar medium (ISM) strongly depend on the chemistry occurring at the surfaces of the dust grains. To describe the complexity of the chemistry, various numerical models have been constructed. In most of these models, the grains are described by a single size of 0.1$\mu$m. We study the impact on the abundances of many species observed in the cold cores by considering several grain sizes in the Nautilus multi-grain model. We used grain sizes with radii in the range of $0.005\mu$m to $0.25\mu$m. We sampled this range in many bins. We used the previously published, MRN and WD grain size distributions to calculate the number density of grains in each bin. Other parameters such as the grain surface temperature or the cosmic-ray-induced desorption rates also vary with grain sizes. We present the abundances of various molecules in the gas phase and also on the dust surface at different time intervals during the simulation. We present a comparative study of results obtained using the single grain and the multi-grain models. We also compare our results with the observed abundances in TMC-1 and L134N clouds. We show that the grain size, the grain size dependent surface temperature and the peak surface temperature induced by cosmic ray collisions, play key roles in determining the ice and the gas phase abundances of various molecules. We also show that the differences between the MRN and the WD models are crucial for better fitting the observed abundances in different regions in the ISM. We show that the small grains play a very important role in the enrichment of the gas phase with the species which are mainly formed on the grain surface, as non-thermal desorption induced by collisions of cosmic ray particles is very efficient on the small grains.


title: Learning to Route with Sparse Trajectory Sets---Extended Version

abstract: Motivated by the increasing availability of vehicle trajectory data, we propose learn-to-route, a comprehensive trajectory-based routing solution. Specifically, we first construct a graph-like structure from trajectories as the routing infrastructure. Second, we enable trajectory-based routing given an arbitrary (source, destination) pair. In the first step, given a road network and a collection of trajectories, we propose a trajectory-based clustering method that identifies regions in a road network. If a pair of regions are connected by trajectories, we maintain the paths used by these trajectories and learn a routing preference for travel between the regions. As trajectories are skewed and sparse, many region pairs are not connected by trajectories. We thus transfer routing preferences from region pairs with sufficient trajectories to such region pairs and then use the transferred preferences to identify paths between the regions. In the second step, we exploit the above graph-like structure to achieve a comprehensive trajectory-based routing solution. Empirical studies with two substantial trajectory data sets offer insight into the proposed solution, indicating that it is practical. A comparison with a leading routing service offers evidence that the paper's proposal is able to enhance routing quality. This is an extended version of "Learning to Route with Sparse Trajectory Sets" [1], to appear in IEEE ICDE 2018.


title: Unified Dark Energy and Dark Matter from Dynamical Space Time

abstract: A unification of dark matter and dark energy based on a dynamical space time theory is suggested. By introducing a dynamical space time vector field $\chi_\mu$ as a Lagrange multiplier, a conservation of an energy momentum tensor $T^{\mu\nu}_{(\chi)}$ is implemented. This Lagrangian generalizes the "Unified dark energy and dark matter from a scalar field different from quintessence" [Phys.RevD 81, 043520 (2010)] which did not consider a Lagrangian formulation. This generalization allows the solutions which were found previously, but in addition to that also non singular bouncing solutions that rapidly approach to the $\Lambda$CDM model. The dynamical time vector field exactly coincides with the cosmic time for the a $\Lambda$CDM solution and suffers a slight shift (advances slower) with respect to the cosmic time in the region close to the bounce for the bouncing non singular solutions. In addition we introduced some exponential potential which could enter into the $T^{\mu\nu}_{(\chi)}$ stress energy tensor or coupled directly to the measure $\sqrt{-g}$, gives a possible interaction between DE and DM and could explain the coincidence problem.


title: Shared Services Center for E-Government Policy

abstract: It is a general opinion that applicative cooperation represents a useful vehicle for the development of e-government. At the architectural level, solutions for applicative cooperation are quite stable, but organizational and methodological problems prevent the expected and needed development of cooperation among different administrations. Moreover, the introduction of the digital government requires a considerable involvement of resources that can be unsustainable for small public administrations. This work shows how the above mentioned problems can be (partially) solved with the introduction of a Shared Services Center (SSC).


title: Tapir: Automation Support of Exploratory Testing Using Model Reconstruction of the System Under Test

abstract: For a considerable number of software projects, the creation of effective test cases is hindered by design documentation that is either lacking, incomplete or obsolete. The exploratory testing approach can serve as a sound method in such situations. However, the efficiency of this testing approach strongly depends on the method, the documentation of explored parts of a system, the organization and distribution of work among individual testers on a team, and the minimization of potential (very probable) duplicities in performed tests. In this paper, we present a framework for replacing and automating a portion of these tasks. A screen-flow-based model of the tested system is incrementally reconstructed during the exploratory testing process by tracking testers' activities. With additional metadata, the model serves for an automated navigation process for a tester. Compared with the exploratory testing approach, which is manually performed in two case studies, the proposed framework allows the testers to explore a greater extent of the tested system and enables greater detection of the defects present in the system. The results show that the time efficiency of the testing process improved with framework support. This efficiency can be increased by team-based navigational strategies that are implemented within the proposed framework, which is documented by another case study presented in this paper.


title: A new optimization problem in FSO communication system

abstract: According to the physical phenomena of atmospheric channels and wave propagation, performance of wireless communication systems can be optimized by simply adjusting its parameters. This way is more economically favorable than consuming power or using processing techniques. In this paper for the first time an optimization problem is developed on the performance of free-space optical multi-input multi-output (FSO-MIMO) communication system. Also it is the first time that the optimization of FSO is developed under saturated atmospheric turbulences. In order to get closer to the actual results, the effect of pointing error is taken into considerations. Assuming MPSK, DPSK modulation schemes, new closed-form expressions are derived for Bit Error Rate (BER) of the proposed structure. Furthermore, an optimization is developed taking into account the beam width as the variable parameter, and BER as the objective function, there is no constraint in this system. The obtained results can be a useful outcome for FSO-MIMO system designers in order to limit effects of pointing error as well as atmospheric turbulences and thus achieves optimum performance.


title: Community Detection with Metadata in a Network of Biographies of Western Art Painters

abstract: In this work we look at the structure of the influences between Western art painters as revealed by their biographies on Wikipedia. We use a modified version of modularity maximisation with metadata to detect a partition of artists into communities based on their artistic genre and school in which they belong. We then use this community structure to discuss how influential artists reached beyond their own communities and had a lasting impact on others, by proposing modifications on standard centrality measures.


title: A Version of $\kappa$-Miller Forcing

abstract: Let $\kappa$ be an uncountable cardinal such that $2^{<\kappa} = \kappa$ or just ${\rm cf}(\kappa) > \omega$, $2^{2^{<\kappa}}= 2^\kappa$, and $([\kappa]^\kappa, \supseteq)$ collapses $2^\kappa$ to $\omega$. We show under these assumptions the $\kappa$-Miller forcing with club many splitting nodes collapses $2^\kappa$ to $\omega$ and adds a $\kappa$-Cohen real.


title: Invariant surfaces in Euclidean space with a log-linear density

abstract: A $\lambda$-translating soliton with density vector $\vec{v}$ is a surface in Euclidean space whose mean curvature $H$ satisfies $2H=2\lambda+\langle N,\vec{v}\rangle$, where $N$ is the Gauss map. We classify all $\lambda$-translating solitons that are invariant by a one-parameter group of translations and a one-parameter group of rotations.


title: Bosonic quantum Hall states in single-layer two-dimensional optical lattices

abstract: Quantum Hall (QH) states of 2D single layer optical lattices are examined using Bose-Hubbard model (BHM) in presence of artificial gauge field. We study the QH states of both the homogeneous and inhomogeneous systems. For the homogeneous case we use cluster Gutzwiller mean field (CGMF) theory with cluster sizes ranging from $2\times 2$ to $5\times 5$. We, then, consider the inhomogeneous case, which is relevant to experimental realization. In this case, we use CGMF and exact diagonalization (ED). The ED studies are using lattice sizes ranging from $3\times 3$ to $4\times 12$. Our results show that the geometry of the QH states are sensitive to the magnetic flux $\alpha$ and cluster sizes. For homogeneous system, among various combinations of $1/5\leqslant \alpha\leqslant 1/2$ and filling factor $\nu$, only the QH state of $\alpha=1/4$ with $\nu=1/2$, $1$, $3/2$ and $2$ occur as ground states. For other combinations, the competing superfluid (SF) state is the ground state and QH state is metastable. For BHM with envelope potential all the QH states observed in homogeneous system exist for box potentials, but none for the harmonic potential. The QH states also persist for very shallow Gaussian envelope potential. As a possible experimental signature we study the two point correlations of the QH and SF states.


title: Electronic properties of phosphorene and graphene nanoribbons with edge vacancies in magnetic field

abstract: The graphene and phosphorene nanostructures have a big potential application in a large area of actuals research in physics. However, their methods of synthesis still do not allow the production of perfect materials with an intact molecular structure. In this paper, the occurrence of atomic vacancies was considered in the edge structure of the zigzag phosphorene and graphene nanoribbons. For different concentrations of these edge vacancies, their influence on the metallic properties was investigated. The calculations were performed for different sizes of the unit cell. Furthermore, for a smaller size, the influence of a uniform magnetic field was added.


title: Joint Antenna Selection and Phase-Only Beamforming Using Mixed-Integer Nonlinear Programming

abstract: In this paper, we consider the problem of joint antenna selection and analog beamformer design in downlink single-group multicast networks. Our objective is to reduce the hardware costs by minimizing the number of required phase shifters at the transmitter while fulfilling given distortion limits at the receivers. We formulate the problem as an L0 minimization problem and devise a novel branch-and-cut based algorithm to solve the resulting mixed-integer nonlinear program to optimality. We also propose a suboptimal heuristic algorithm to solve the above problem approximately with a low computational complexity. Computational results illustrate that the solutions produced by the proposed heuristic algorithm are optimal in most cases. The results also indicate that the performance of the optimal methods can be significantly improved by initializing with the result of the suboptimal method.


title: Decomposition of a graph into two disjoint odd subgraphs

abstract: An odd (resp. even) subgraph in a multigraph is its subgraph in which every vertex has odd (resp. even) degree. We say that a multigraph can be decomposed into two odd subgraphs if its edge set can be partitioned into two sets so that both form odd subgraphs. In this paper we give a necessary and sufficient condition for the decomposability of a multigraph into two odd subgraphs. We also present a polynomial time algorithm for finding such a decomposition or showing its non-existence. We also deal with the case of the decomposability into an even subgraph and an odd subgraph.


title: Rodin's formula in arbitrary codimension

abstract: We extend the Rodin's formula for $p$--modulus of the family of curves in $\mathbb{R}^n$ to arbitrary codimension. The proof relies on the formula for the $p$--modulus of family of level sets of a submersion and an algebraic lemma relating Jacobi matrices of considered maps. We state appropriate examples.


title: Tip-induced Superconductivity Coexisting with Preserved Topological Properties in Line-nodal Semimetal ZrSiS

abstract: ZrSiS was recently shown to be a new material with topologically non-trivial band structure which exhibits multiple Dirac nodes and a robust linear band dispersion up to an unusually high energy of 2\,eV. Such a robust linear dispersion makes the topological properties of ZrSiS insensitive to perturbations like carrier doping or lattice distortion. Here we show that a novel superconducting phase with a remarkably high $T_c$ of 7.5\,K can be induced in single crystals of ZrSiS by a non-superconducting metallic tip of Ag. From first-principles calculations we show that the observed superconducting phase might originate from dramatic enhancement of density of states due to the presence of a metallic tip on ZrSiS. Our calculations also show that the emerging tip-induced superconducting phase co-exists with the well preserved topological properties of ZrSiS.


title: Compact $\lambda$-translating solitons with boundary

abstract: A $\lambda$-translating soliton with density vector $\vec{v}$ is a surface $\Sigma$ in Euclidean space ${\mathbb R}^3$ whose mean curvature $H$ satisfies $2H=2\lambda+\langle N,\vec{v}\rangle$, where $N$ is the Gauss map of $\Sigma$. In this article we study the shape of a compact $\lambda$-translating soliton in terms of its boundary. If $\Gamma$ is a given closed curve, we deduce under what conditions on $\lambda$ there exists a compact $\lambda$-translating soliton $\Sigma$ with boundary $\Gamma$ and we provide estimates of the surface area in relation with the height of $\Sigma$. Finally we study the shape of $\Sigma$ related with the one of $\Gamma$, in particular, we give conditions that assert that $\Sigma$ inherits the symmetries of its boundary $\Gamma$.


title: Multidimensional multiscale scanning in Exponential Families: Limit theory and statistical consequences

abstract: We consider the problem of finding anomalies in a $d$-dimensional field of independent random variables $\{Y_i\}_{i \in \left\{1,...,n\right\}^d}$, each distributed according to a one-dimensional natural exponential family $\mathcal F = \left\{F_\theta\right\}_{\theta \in\Theta}$. Given some baseline parameter $\theta_0 \in\Theta$, the field is scanned using local likelihood ratio tests to detect from a (large) given system of regions $\mathcal{R}$ those regions $R \subset \left\{1,...,n\right\}^d$ with $\theta_i \neq \theta_0$ for some $i \in R$. We provide a unified methodology which controls the overall family wise error (FWER) to make a wrong detection at a given error rate. Fundamental to our method is a Gaussian approximation of the distribution of the underlying multiscale test statistic with explicit rate of convergence. From this, we obtain a weak limit theorem which can be seen as a generalized weak invariance principle to non identically distributed data and is of independent interest. Furthermore, we give an asymptotic expansion of the procedures power, which yields minimax optimality in case of Gaussian observations.


title: Multicolor photometric behavior of the young stellar object V1704 Cygni

abstract: Results from BVRI photometric observations of the pre-main sequence star V1704 Cyg collected during the time period from August 2010 to December 2017 are presented. The star is located in the star-forming HII region IC 5070 and it exhibits photometric variability in all-optical passbands. After analyzing the obtained data, V1704 Cyg is classified as a classical T Tauri star.


title: Generating High-Quality Query Suggestion Candidates for Task-Based Search

abstract: We address the task of generating query suggestions for task-based search. The current state of the art relies heavily on suggestions provided by a major search engine. In this paper, we solve the task without reliance on search engines. Specifically, we focus on the first step of a two-stage pipeline approach, which is dedicated to the generation of query suggestion candidates. We present three methods for generating candidate suggestions and apply them on multiple information sources. Using a purpose-built test collection, we find that these methods are able to generate high-quality suggestion candidates.


title: Robust estimators in a generalized partly linear regression model under monotony constraints

abstract: In this paper, we consider the situation in which the observations follow an isotonic generalized partly linear model. Under this model, the mean of the responses is modelled, through a link function, linearly on some covariates and nonparametrically on an univariate regressor in such a way that the nonparametric component is assumed to be a monotone function. A class of robust estimates for the monotone nonparametric component and for the regression parameter, related to the linear one, is defined. The robust estimators are based on a spline approach combined with a score function which bounds large values of the deviance. As an application, we consider the isotonic partly linear log--Gamma regression model. Through a Monte Carlo study, we investigate the performance of the proposed estimators under a partly linear log-Gamma regression model with increasing nonparametric component.


title: Differential technique for the covariant orbital angular momentum operators

abstract: The orbital angular momentum operator expansion turns to be a powerful tool to construct the fully covariant partial wave amplitudes of hadron decay reactions and hadron photo- and electroproduction processes. In this paper we consider a useful development of the orbital angular momentum operator expansion method. We present the differential technique allowing the direct calculation of convolutions of two orbital angular momentum operators with an arbitrary number of open Lorentz indices. This differential technique greatly simplifies calculations when the reaction subject to the partial wave analysis involves high spin particles in the initial and/or final states. We also present a useful generalization of the orbital angular momentum operators.


title: Gravitational decoupled anisotropies in compact stars

abstract: Simple generic extensions of isotropic Durgapal--Fuloria stars to the anisotropic domain are presented. These anisotropic solutions are obtained by guided minimal deformations over a self gravitating isotropic system. When the isotropic and the anisotropic sector interacts in a purely gravitational manner, the conditions to decouple both sectors by means of the minimal geometric deformation approach are satisfied. Hence the anisotropic field equations are isolated resulting a more treatable set. The simplicity of the equations allows one to manipulate the anisotropies that can be implemented in a systematic way to obtain different realistic models for anisotropic configurations. Later on, observational effects of such anisotropies when measuring the redshift are discussed. To conclude, the application of the method over anisotropic solutions is generalized. In this manner, different anisotropic sectors can be isolated of each other and modeled in a simple and systematic way. Besides, a generic property of the minimal geometric deformation approach, its noncommutativity, is discussed. This property duplicates the solutions obtained through this approach; anisotropies applied in the reversed order give different physically acceptable configurations.


title: On the permanent of Sylvester-Hadamard matrices

abstract: We prove a conjecture due to Wanless about the permanent of Hadamard matrices in the particular case of Sylvester-Hadamard matrices. Namely we show that for all n greater or equal to 2, the dyadic valuation of the permanent of the Sylvester-Hadamard matrix of order n is equal to the dyadic valuation of n!. As a consequence, the permanent of the Sylvester-Hadamard matrix of order n doesn't vanish for n greater or equal to 2.


title: Correlation of tunnel magnetoresistance with the magnetic properties in perpendicular CoFeB-based junctions with exchange bias

abstract: We investigate the dependence of magnetic properties on the post-annealing temperature/time, the thickness of soft ferromagnetic electrode and Ta dusting layer in the pinned electrode as well as their correlation with the tunnel magnetoresistance ratio, in a series of perpendicular magnetic tunnel junctions of materials sequence Ta/Pd/IrMn/CoFe/Ta$(\textit{x})$/CoFeB/MgO$(\textit{y})$/CoFeB$(\textit{z})$/Ta/Pd. We obtain a large perpendicular exchange bias of 79.6$\,$kA/m for $x=0.3\,$nm. For stacks with $z=1.05\,$nm, the magnetic properties of the soft electrode resemble the characteristics of superparamagnetism. For stacks with $x=0.4\,$nm, $y=2\,$nm, and $z=1.20\,$nm, the exchange bias presents a significant decrease at post annealing temperature $T_\textrm{ann}=330\,^{\circ}$C for 60 min, while the interlayer exchange coupling and the saturation magnetization per unit area sharply decay at $T_\textrm{ann}=340\,^{\circ}$C for 60 min. Simultaneously, the tunnel magnetoresistance ratio shows a peak of $65.5\%$ after being annealed at $T_\textrm{ann}=300\,^{\circ}$C for 60 min, with a significant reduction down to $10\%$ for higher annealing temperatures ($T_\textrm{ann}\geq330\,^{\circ}$C) and down to $14\%$ for longer annealing times ($T_\textrm{ann}=300\,^{\circ}$C for 90 min). We attribute the large decrease of tunnel magnetoresistance ratio to the loss of exchange bias in the pinned electrode.


title: Strong Valley Zeeman Effect of Dark Excitons in Monolayer Transition Metal Dichalcogenides in a Tilted Magnetic Field

abstract: The dependence of the excitonic photoluminescence (PL) spectrum of monolayer transition metal dichalcogenides (TMDs) on the tilt angle of an applied magnetic field is studied. Starting from a four-band Hamiltonian we construct a theory which quantitatively reproduces the available experimental PL spectra for perpendicular and in-plane magnetic fields. In the presence of a tilted magnetic field, we demonstrate that the dark exciton PL peaks brighten due to the in-plane component of the magnetic field and split for light with different circular polarization as a consequence of the perpendicular component of the magnetic field. This splitting is more than twice as large as the splitting of the bright exciton peaks in tungsten-based TMDs. We propose an experimental setup that will allow to access the predicted splitting of the dark exciton peaks in the PL spectrum.


title: The use of sampling weights in the M-quantile random-effects regression: an application to PISA mathematics scores

abstract: M-quantile random-effects regression represents an interesting approach for modelling multilevel data when the interest of researchers is focused on the conditional quantiles. When data are based on complex survey designs, sampling weights have to be incorporate in the analysis. A pseudo-likelihood approach for accommodating sampling weights in the M-quantile random-effects regression is presented. The proposed methodology is applied to the Italian sample of the "Program for International Student Assessment 2015" survey in order to study the gender gap in mathematics at various quantiles of the conditional distribution. Findings offer a possible explanation of the low share of females in "Science, Technology, Engineering and Mathematics" sectors.


title: Employment of Multiple Algorithms for Optimal Path-based Test Selection Strategy

abstract: Executing various sequences of system functions in a system under test represents one of the primary techniques in software testing. The natural way to create effective, consistent and efficient test sequences is to model the system under test and employ an algorithm to generate the tests that satisfy a defined test coverage criterion. Several criteria of test set optimality can be defined. In addition, to optimize the test set from an economic viewpoint, the priorities of the various parts of the system model under test must be defined. Using this prioritization, the test cases exercise the high priority parts of the system under test more intensely than those with low priority. Evidence from the literature and our observations confirm that finding a universal algorithm that produces an optimal test set for all test coverage and test set optimality criteria is a challenging task. Moreover, for different individual problem instances, different algorithms provide optimal results. In this paper, we present a path-based strategy to perform optimal test selection. The strategy first employs a set of current algorithms to generate test sets; then, it assesses the optimality of each test set by the selected criteria, and finally, chooses the optimal test set. The experimental results confirm the validity and usefulness of this strategy. For individual instances of 50 system under test models, different algorithms provided optimal results; these results varied by the required test coverage level, the size of the priority parts of the model, and the selected test set optimality criteria.


title: Heat transport in insulators from ab initio Green-Kubo theory

abstract: The Green-Kubo theory of thermal transport has long be considered incompatible with modern simulation methods based on electronic-structure theory, because it is based on such concepts as energy density and current, which are ill-defined at the quantum-mechanical level. Besides, experience with classical simulations indicates that the estimate of heat-transport coefficients requires analysing molecular trajectories that are more than one order of magnitude longer than deemed feasible using ab initio molecular dynamics. In this paper we report on recent theoretical advances that are allowing one to overcome these two obstacles. First, a general gauge invariance principle has been established, stating that thermal conductivity is insensitive to many details of the microscopic expression for the energy density and current from which it is derived, thus permitting to establish a rigorous expression for the energy flux from Density-Functional Theory, from which the conductivity can be computed in practice. Second, a novel data analysis method based on the statistical theory of time series has been proposed, which allows one to considerably reduce the simulation time required to achieve a target accuracy on the computed conductivity. These concepts are illustrated in detail, starting from a pedagogical introduction to the Green-Kubo theory of linear response and transport, and demonstrated with a few applications done with both classical and quantum-mechanical simulation methods.


title: Dependence of $T_c$ on the $q-\omega$ structure of the spin-fluctuation spectrum

abstract: A phenomenological spin-fluctuation analysis (Ref. 1), based upon inelastic neutron scattering (INS) and angular resolved photoemission spectroscopy (ARPES) data for ${\rm YBCO}_{6.6}(T_c=61K)$, is used to calculate the functional derivative of the d-wave eigenvalue $\lambda_d$ of the linearized gap equation with respect to the imaginary part of the spin susceptibility $\chi''(q,\omega)$ at 70K. For temperatures near $T_c$, the variation of $T_c$ with respect to $\chi''(q,\omega)$ is proportional to this functional derivative. Based on this, we discuss how different parts of the $q$ and $\omega$ dependent spin-fluctuation spectrum of YBCO$_{6.6}$ contribute to $T_c$.


title: Sounderfeit: Cloning a Physical Model with Conditional Adversarial Autoencoders

abstract: An adversarial autoencoder conditioned on known parameters of a physical modeling bowed string synthesizer is evaluated for use in parameter estimation and resynthesis tasks. Latent dimensions are provided to capture variance not explained by the conditional parameters. Results are compared with and without the adversarial training, and a system capable of "copying" a given parameter-signal bidirectional relationship is examined. A real-time synthesis system built on a generative, conditioned and regularized neural network is presented, allowing to construct engaging sound synthesizers based purely on recorded data.


title: Iterate averaging as regularization for stochastic gradient descent

abstract: We propose and analyze a variant of the classic Polyak-Ruppert averaging scheme, broadly used in stochastic gradient methods. Rather than a uniform average of the iterates, we consider a weighted average, with weights decaying in a geometric fashion. In the context of linear least squares regression, we show that this averaging scheme has a the same regularizing effect, and indeed is asymptotically equivalent, to ridge regression. In particular, we derive finite-sample bounds for the proposed approach that match the best known results for regularized stochastic gradient methods.


title: Towards an Understanding of Entity-Oriented Search Intents

abstract: Entity-oriented search deals with a wide variety of information needs, from displaying direct answers to interacting with services. In this work, we aim to understand what are prominent entity-oriented search intents and how they can be fulfilled. We develop a scheme of entity intent categories, and use them to annotate a sample of queries. Specifically, we annotate unique query refiners on the level of entity types. We observe that, on average, over half of those refiners seek to interact with a service, while over a quarter of the refiners search for information that may be looked up in a knowledge base.


title: Ruijsenaars-Schneider three-body models with N=2 supersymmetry

abstract: The Ruijsenaars-Schneider models are conventionally regarded as relativistic generalizations of the Calogero integrable systems. Surprisingly enough, their supersymmetric generalizations escaped attention. In this work, N=2 supersymmetric extensions of the rational and hyperbolic Ruijsenaars-Schneider three-body models are constructed within the framework of the Hamiltonian formalism. It is also known that the rational model can be described by the geodesic equations associated with a metric connection. We demonstrate that the hyperbolic systems are linked to non-metric connections.


title: Learning Topic Models by Neighborhood Aggregation

abstract: Topic models are frequently used in machine learning owing to their high interpretability and modular structure. However, extending a topic model to include a supervisory signal, to incorporate pre-trained word embedding vectors and to include a nonlinear output function is not an easy task because one has to resort to a highly intricate approximate inference procedure. The present paper shows that topic modeling with pre-trained word embedding vectors can be viewed as implementing a neighborhood aggregation algorithm where messages are passed through a network defined over words. From the network view of topic models, nodes correspond to words in a document and edges correspond to either a relationship describing co-occurring words in a document or a relationship describing the same word in the corpus. The network view allows us to extend the model to include supervisory signals, incorporate pre-trained word embedding vectors and include a nonlinear output function in a simple manner. In experiments, we show that our approach outperforms the state-of-the-art supervised Latent Dirichlet Allocation implementation in terms of held-out document classification tasks.


title: Intrinsic Motivation and Mental Replay enable Efficient Online Adaptation in Stochastic Recurrent Networks

abstract: Autonomous robots need to interact with unknown, unstructured and changing environments, constantly facing novel challenges. Therefore, continuous online adaptation for lifelong-learning and the need of sample-efficient mechanisms to adapt to changes in the environment, the constraints, the tasks, or the robot itself are crucial. In this work, we propose a novel framework for probabilistic online motion planning with online adaptation based on a bio-inspired stochastic recurrent neural network. By using learning signals which mimic the intrinsic motivation signalcognitive dissonance in addition with a mental replay strategy to intensify experiences, the stochastic recurrent network can learn from few physical interactions and adapts to novel environments in seconds. We evaluate our online planning and adaptation framework on an anthropomorphic KUKA LWR arm. The rapid online adaptation is shown by learning unknown workspace constraints sample-efficiently from few physical interactions while following given way points.


title: Finding Top-k Optimal Sequenced Routes -- Full Version

abstract: Motivated by many practical applications in logistics and mobility-as-a-service, we study the top-k optimal sequenced routes (KOSR) querying on large, general graphs where the edge weights may not satisfy the triangle inequality, e.g., road network graphs with travel times as edge weights. The KOSR querying strives to find the top-k optimal routes (i.e., with the top-k minimal total costs) from a given source to a given destination, which must visit a number of vertices with specific vertex categories (e.g., gas stations, restaurants, and shopping malls) in a particular order (e.g., visiting gas stations before restaurants and then shopping malls). To efficiently find the top-k optimal sequenced routes, we propose two algorithms PruningKOSR and StarKOSR. In PruningKOSR, we define a dominance relationship between two partially-explored routes. The partially-explored routes that can be dominated by other partially-explored routes are postponed being extended, which leads to a smaller searching space and thus improves efficiency. In StarKOSR, we further improve the efficiency by extending routes in an A* manner. With the help of a judiciously designed heuristic estimation that works for general graphs, the cost of partially explored routes to the destination can be estimated such that the qualified complete routes can be found early. In addition, we demonstrate the high extensibility of the proposed algorithms by incorporating Hop Labeling, an effective label indexing technique for shortest path queries, to further improve efficiency. Extensive experiments on multiple real-world graphs demonstrate that the proposed methods significantly outperform the baseline method. Furthermore, when k=1, StarKOSR also outperforms the state-of-the-art method for the optimal sequenced route queries.


title: Spanned lines and Langer's inequality

abstract: We collect some results in combinatorial geometry that follow from an inequality of Langer in algebraic geometry. Langer's inequality gives a lower bound on the number of incidences between a point set and its spanned lines, and was recently used by Han to improve the constant in the weak Dirac conjecture. Here we observe that this inequality also leads to improved constants in Beck's theorem, which states that a finite point set in the real or complex plane has many points on a line or spans many lines. Most of the proofs that we use are not original, and the goal of this note is mainly to carefully record the quantitative results in one place. We also include some discussion of possible further improvements to these statements.


title: Microscopic characterisation of laser-written phenomena for component-wise testing of photonic integrated circuits

abstract: Photonic integrated circuits (PICs) directly written with a femtosecond laser have shown great potential in many areas such as quantum information processing (QIP). Many applications, like photon-based quantum computing, demand the up-scaling of PICs and ever-higher optical performance, such as controllable polarisation dependence and lower loss. In order to overcome current limitations in fabrication precision, repeatability and material uniformity, a solution for non-destructive testing of large-scale PICs in a component-wise manner is desired to meet those ever-stricter demands. Here we demonstrate a solution for non-destructive component-wise testing by predicting the performance of a PIC component based on imaging with an adaptive optical third-harmonic-generation (THG) three-dimensional (3D) microscope. The 3D THG imaging can be performed on any component or part of it inside multi-component PIC. Moreover, through discovering new phenomena we also demonstrated that 3D THG microscopy provides a new pathway towards studying the fundamentals of light-matter interaction in transparent materials.


title: On the classification of almost contact metric manifolds

abstract: On connected manifolds of dimension higher than three, the non-existence of $132$ Chinea and Gonz\'alez-D\'avila types of almost contact metric structures is proved. This is a consequence of some interrelations among components of the intrinsic torsion of an almost contact metric structure. Such interrelations allow to describe the exterior derivatives of some relevant forms in the context of almost contact metric geometry.


title: Structure and Supersaturation for Intersecting Families

abstract: The extremal problems regarding the maximum possible size of intersecting families of various combinatorial objects have been extensively studied. In this paper, we investigate supersaturation extensions, which in this context ask for the minimum number of disjoint pairs that must appear in families larger than the extremal threshold. We study the minimum number of disjoint pairs in families of permutations and in $k$-uniform set families, and determine the structure of the optimal families. Our main tool is a removal lemma for disjoint pairs. We also determine the typical structure of $k$-uniform set families without matchings of size $s$ when $n \ge 2sk + 38s^4$, and show that almost all $k$-uniform intersecting families on vertex set $[n]$ are trivial when $n\ge (2+o(1))k$.


title: Neutrinoless double-beta decay with massive scalar emission

abstract: Searches for neutrino-less double-beta decay ($0\nu2\beta$) place an important constraint on models where light fields beyond the Standard Model participate in the neutrino mass mechanism. While $0\nu2\beta$ experimental collaborations often consider various massless majoron models, including various forms of majoron couplings and multi-majoron final-state processes, none of these searches considered the scenario where the "majoron" $\phi$ is not massless, $m_\phi\sim$~MeV, of the same order as the $Q$-value of the $0\nu2\beta$ reaction. We consider this parameter region and estimate $0\nu2\beta\phi$ constraints for $m_\phi$ of order MeV. The constraints are affected not only by kinematical phase space suppression but also by a change in the signal to background ratio characterizing the search. As a result, $0\nu2\beta\phi$ constraints for $m_\phi>0$ diminish significantly below the reaction threshold. This has phenomenological implications, which we illustrate focusing on high-energy neutrino telescopes. Our results motivate a dedicated analysis by $0\nu2\beta$ collaborations, analogous to the dedicated analyses targeting massless majoron models.


title: On Rational Delegations in Liquid Democracy

abstract: Liquid democracy is a proxy voting method where proxies are delegable. We propose and study a game-theoretic model of liquid democracy to address the following question: when is it rational for a voter to delegate her vote? We study the existence of pure-strategy Nash equilibria in this model, and how group accuracy is affected by them. We complement these theoretical results by means of agent-based simulations to study the effects of delegations on group's accuracy on variously structured social networks.


title: SparCML: High-Performance Sparse Communication for Machine Learning

abstract: Applying machine learning techniques to the quickly growing data in science and industry requires highly-scalable algorithms. Large datasets are most commonly processed "data parallel" distributed across many nodes. Each node's contribution to the overall gradient is summed using a global allreduce. This allreduce is the single communication and thus scalability bottleneck for most machine learning workloads. We observe that frequently, many gradient values are (close to) zero, leading to sparse of sparsifyable communications. To exploit this insight, we analyze, design, and implement a set of communication-efficient protocols for sparse input data, in conjunction with efficient machine learning algorithms which can leverage these primitives. Our communication protocols generalize standard collective operations, by allowing processes to contribute arbitrary sparse input data vectors. Our generic communication library, SparCML, extends MPI to support additional features, such as non-blocking (asynchronous) operations and low-precision data representations. As such, SparCML and its techniques will form the basis of future highly-scalable machine learning frameworks.


title: Equalizer 2.0 - Convergence of a Parallel Rendering Framework

abstract: Developing complex, real world graphics applications which leverage multiple GPUs and computers for interactive 3D rendering tasks is a complex task. It requires expertise in distributed systems and parallel rendering in addition to the application domain itself. We present a mature parallel rendering framework which provides a large set of features, algorithms and system integration for a wide range of real-world research and industry applications. Using the Equalizer parallel rendering framework, we show how a wide set of generic algorithms can be integrated in the framework to help application scalability and development in many different domains, highlighting how concrete applications benefit from the diverse aspects and use cases of Equalizer. We present novel parallel rendering algorithms, powerful abstractions for large visualization setups and virtual reality, as well as new experimental results for parallel rendering and data distribution.


title: The Best of Both Worlds: Asymptotically Efficient Mechanisms with a Guarantee on the Expected Gains-From-Trade

abstract: The seminal impossibility result of Myerson and Satterthwaite (1983) states that for bilateral trade, there is no mechanism that is individually rational (IR), incentive compatible (IC), weakly budget balanced, and efficient. This has led follow-up work on two-sided trade settings to weaken the efficiency requirement and consider approximately efficient simple mechanisms, while still demanding the other properties. The current state-of-the-art of such mechanisms for two-sided markets can be categorized as giving one (but not both) of the following two types of approximation guarantees on the gains from trade: a constant ex-ante guarantee, measured with respect to the second-best efficiency benchmark, or an asymptotically optimal ex-post guarantee, measured with respect to the first-best efficiency benchmark. Here the second-best efficiency benchmark refers to the highest gains from trade attainable by any IR, IC and weakly budget balanced mechanism, while the first-best efficiency benchmark refers to the maximum gains from trade (attainable by the VCG mechanism, which is not weakly budget balanced). In this paper, we construct simple mechanisms for double-auction and matching markets that simultaneously achieve both types of guarantees: these are ex-post IR, Bayesian IC, and ex-post weakly budget balanced mechanisms that 1) ex-ante guarantee a constant fraction of the gains from trade of the second-best, and 2) ex-post guarantee a realization-dependent fraction of the gains from trade of the first-best, such that this realization-dependent fraction converges to 1 (full efficiency) as the market grows large.


title: Chemical Heredity as Group Selection at the Molecular Level

abstract: Many examples of cooperation exist in biology. In chemical systems however, which can sometimes be quite complex, we do not appear to observe intricate cooperative interactions. A key question for the origin of life, is then how can molecular cooperation first arise in an abiotic system prior to the emergence of biological replication. We postulate that selection at the molecular level is a driving force behind the complexification of chemical systems, particularly during the origins of life. In the theory of multilevel selection the two selective forces are: within-group and between-group, where the former tends to favor "selfish" replication of individuals and the latter favor cooperation between individuals enhancing the replication of the group as a whole. These forces can be quantified using the Price equation, which is a standard tool used in evolutionary biology to quantify evolutionary change. Our central claim is that replication and heredity in chemical systems are subject to selection, and quantifiable using the multilevel Price equation. We demonstrate this using the Graded Autocatalysis Replication Domain computer model, describing simple protocell composed out of molecules and its replication, which respectively analogue to the group and the individuals. In contrast to previous treatments of this model, we treat the lipid molecules themselves as replicating individuals and the protocells they form as groups of individuals. Our goal is to demonstrate how evolutionary biology tools and concepts can be applied in chemistry and we suggest that molecular cooperation may arise as a result of group selection. Further, the biological relation of parent-progeny is proposed to be analogue to the reactant-product relation in chemistry, thus allowing for tools from evolutionary biology to be applied to chemistry and would deepen the connection between chemistry and biology.


title: Dynamics of driven translocation of semiflexible polymers

abstract: We study translocation of semiflexible polymers driven by force $f_d$ inside a nanometer-scale pore using our three-dimensional Langevin dynamics model. We show that the translocation time $\tau$ increases with increasing bending rigidity $\kappa$. Similarly, the exponent $\beta$ for the scaling of $\tau$ with polymer length $N$, $\tau \sim N^\beta$, increases with increasing $\kappa$ as well as with increasing $f_d$. By comparing waiting times between semiflexible and fully flexible polymers we show that for realistic $f_d$ translocation dynamics is to a large extent, but not completely, determined by the polymer's elastic length measured in number of Kuhn segments $N_{\rm Kuhn}$. Unlike in driven translocation of flexible polymers, friction related to the polymer segment on the trans side has a considerable effect on the resulting dynamics. This friction is intermittently reduced by buckling of the polymer segment in the vicinity of the pore opening on the trans side. We show that in the experimentally relevant regime, where viscosity is higher than in computer simulation models, the probability for this buckling increases with increasing $f_d$, giving rise to larger contribution to trans side friction at small $f_d$. Similarly to flexible polymers, we find significant center of mass diffusion of the cis side polymer segment. This speeds up translocation, which effect is larger for smaller $f_d$. However, this speed-up is smaller than the slowing down due to the trans side friction. At large enough $N_{\rm Kuhn}$, the roles can be seen to be reversed and the dynamics of flexible polymers be reached. However, for example, polymers used in translocation experiments of DNA, are elastically so short, that the finite-length dynamics outlined here applies.


title: Complex-valued Neural Networks with Non-parametric Activation Functions

abstract: Complex-valued neural networks (CVNNs) are a powerful modeling tool for domains where data can be naturally interpreted in terms of complex numbers. However, several analytical properties of the complex domain (e.g., holomorphicity) make the design of CVNNs a more challenging task than their real counterpart. In this paper, we consider the problem of flexible activation functions (AFs) in the complex domain, i.e., AFs endowed with sufficient degrees of freedom to adapt their shape given the training data. While this problem has received considerable attention in the real case, a very limited literature exists for CVNNs, where most activation functions are generally developed in a split fashion (i.e., by considering the real and imaginary parts of the activation separately) or with simple phase-amplitude techniques. Leveraging over the recently proposed kernel activation functions (KAFs), and related advances in the design of complex-valued kernels, we propose the first fully complex, non-parametric activation function for CVNNs, which is based on a kernel expansion with a fixed dictionary that can be implemented efficiently on vectorized hardware. Several experiments on common use cases, including prediction and channel equalization, validate our proposal when compared to real-valued neural networks and CVNNs with fixed activation functions.


title: MEC-assisted End-to-End Latency Evaluations for C-V2X Communications

abstract: The efficient design of fifth generation (5G) mobile networks is driven by the need to support the dynamic proliferation of several vertical market segments. Considering the automotive sector, different Cellular Vehicle-to-Everything (C-V2X) use cases have been identified by the industrial and research world, referring to infotainment, automated driving and road safety. A common characteristic of these use cases is the need to exploit collective awareness of the road environment towards satisfying performance requirements. One of these requirements is the End-to-End (E2E) latency when, for instance, Vulnerable Road Users (VRUs) inform vehicles about their status (e.g., location) and activity, assisted by the cellular network. In this paper, focusing on a freeway-based VRU scenario, we argue that, in contrast to conventional, remote cloud-based cellular architecture, the deployment of Multi-access Edge Computing (MEC) infrastructure can substantially prune the E2E communication latency. Our argument is supported by an extensive simulation-based performance comparison between the conventional and the MEC-assisted network architecture.


title: Approximate controllabilty from the exterior of space-time fractional diffusive equations

abstract: Let $\Om\subset\RR^N$ a bounded domain with a Lipschitz continuous boundary. We study the controllability of the space-time fractional diffusion equation \begin{equation*} \begin{cases} \mathbb D_t^\alpha u+(-\Delta)^su=0\;\;&\mbox{ in }\;(0,T)\times\Omega\\ u=g &\mbox{ in }\;(0,T)\times(\RR^N\setminus\Omega)\\ u(0,\cdot)=u_0&\mbox{ in }\;\Omega, \end{cases} \end{equation*} where $u=u(t,x)$ is the state to be controlled and $g=g(t,x)$ is the control function which is localized in a subset $\mathcal O$ of $\Omc$. Here, $0<\alpha\le 1$, $0<s<1$ and $T>0$ be real numbers. After giving an explicit representation of solutions, we show that the system is always approximately controllable for every $T>0$, $u_0\in L^2(\Omega)$ and $g\in \mathcal D((0,T)\times\mathcal O)$ where $\mathcal O\subset(\RR^N\setminus\bOm)$ is any open set. The results obtained are sharp in the sense that such a system is never null controllable if $0<\alpha<1$. The proof of our result is based on a new unique continuation principle for the eigenvalues problem associated with the fractional Laplace operator subject to the zero exterior boundary condition that we have established.


title: Optical Chirality of Time-Harmonic Wavefields for Classification of Scatterers

abstract: We derive expressions for the scattering, extinction and conversion of the chirality of monochromatic light scattered by bodies which are characterized by a T-matrix. In analogy to the conditions obtained from the conservation of energy, these quantities enable the classification of arbitrary scattering objects due to their full, i.e. either chiral or achiral, electromagnetic response. To this end, we put forward and determine the concepts of duality and breaking of duality symmetry, anti-duality, helicity variation, helicity annhiliation and the breaking of helicity annihilation. Different classes, such as chiral and dual scatterers, are illustrated in this analysis with model examples of spherical and non-spherical shape. As for spheres, these concepts are analysed by considering non-Rayleigh dipolar dielectric particles of high refractive index, which, having a strong magnetic response to the incident wavefield, offer an excellent laboratory to test and interpret such changes in the chirality of the illumination. In addition, comparisons with existing experimental data are made.


title: Study of the filamentary infrared dark cloud G192.76+00.10 in the S254-S258 OB complex

abstract: We present results of a high resolution study of the filamentary infrared dark cloud G192.76+00.10 in the S254-S258 OB complex in several molecular species tracing different physical conditions. These include three isotopologues of carbon monoxide (CO), ammonia (NH$_3$), carbon monosulfide (CS). The aim of this work is to study the general structure and kinematics of the filamentary cloud, its fragmentation and physical parameters. The gas temperature is derived from the NH$_3 $ $(J,K) = (1,1), (2,2)$ and $^{12}$CO(2--1) lines and the $^{13}$CO(1--0), $^{13}$CO(2--1) emission is used to investigate the overall gas distribution and kinematics. Several dense clumps are identified from the CS(2--1) data. Values of the gas temperature lie in the ranges $10-35$ K, column density $N(\mathrm{H}_2)$ reaches the value 5.1 10$^{22}$ cm$^{-2}$. The width of the filament is of order 1 pc. The masses of the dense clumps range from $ \sim 30 $ M$_\odot$ to $ \sim 160 $ M$_\odot$. They appear to be gravitationally unstable. The molecular emission shows a gas dynamical coherence along the filament. The velocity pattern may indicate longitudinal collapse.


title: The tangent bundle of a model category

abstract: This paper studies the homotopy theory of parameterized spectrum objects in a model category from a global point of view. More precisely, for a model category $\mathcal{M}$ satisfying suitable conditions, we construct a relative model category $\mathcal{TM} \to \mathcal{M}$, called the tangent bundle, whose fibers are models for spectra in the various over-categories of $\mathcal{M}$, and which presents the $\infty$-categorical tangent bundle. Moreover, the tangent bundle $\mathcal{TM}$ inherits an enriched model structure when such a structure exists on $\mathcal{M}$. This additional structure is used in subsequent work to identify the tangent bundles of algebras over an operad and of enriched categories.


title: QuEST and High Performance Simulation of Quantum Computers

abstract: We introduce QuEST, the Quantum Exact Simulation Toolkit, and compare it to ProjectQ, qHipster and a recent distributed implementation of Quantum++. QuEST is the first open source, OpenMP and MPI hybridised, GPU accelerated simulator of universal quantum circuits. Embodied as a C library, it is designed so that a user's code can be deployed seamlessly to any platform from a laptop to a supercomputer. QuEST is capable of simulating generic quantum circuits of general single-qubit gates and multi-qubit controlled gates, on pure and mixed states, represented as state-vectors and density matrices, and under the presence of decoherence. Using the ARCUS Phase-B and ARCHER supercomputers, we benchmark QuEST's simulation of random circuits of up to 38 qubits, distributed over up to 2048 compute nodes, each with up to 24 cores. We directly compare QuEST's performance to ProjectQ's on single machines, and discuss the differences in distribution strategies of QuEST, qHipster and Quantum++. QuEST shows excellent scaling, both strong and weak, on multicore and distributed architectures.


title: Approximating the nearest stable discrete-time system

abstract: In this paper, we consider the problem of stabilizing discrete-time linear systems by computing a nearby stable matrix to an unstable one. To do so, we provide a new characterization for the set of stable matrices. We show that a matrix $A$ is stable if and only if it can be written as $A=S^{-1}UBS$, where $S$ is positive definite, $U$ is orthogonal, and $B$ is a positive semidefinite contraction (that is, the singular values of $B$ are less or equal to 1). This characterization results in an equivalent non-convex optimization problem with a feasible set on which it is easy to project. We propose a very efficient fast projected gradient method to tackle the problem in variables $(S,U,B)$ and generate locally optimal solutions. We show the effectiveness of the proposed method compared to other approaches.


title: Pseudo-Harmonic Maps From Complete Noncompact Pseudo-Hermitian Manifolds To Regular Balls

abstract: In this paper, we give an estimate of sub-Laplacian of Riemannian distance functions in pseudo-Hermitian geometry which plays a similar role as Laplacian comparison theorem in Riemannian geometry, and deduce a prior horizontal gradient estimate of pseudo-harmonic maps from pseudo-Hermitian manifolds to regular balls of Riemannian manifolds. As an application, Liouville theorem is established under the conditions of nonnegative pseudo-Hermitian Ricci curvature and vanishing pseudo-Hermitian torsion. Moreover, we obtain the existence of pseudo-harmonic maps from complete noncompact pseudo-Hermitian manifolds to regular balls of Riemannian manifolds.


title: Horizontal Gradient Estimate of Positive Pseudo-Harmonic Functions on Complete Noncompact Pseudo-Hermitian Manifolds

abstract: In this paper, we will give a horizontal gradient estimate of positive solutions of $\Delta_b u = - \lambda u$ on complete noncompact pseudo-Hermitian manifolds. As a consequence, we recapture the Liouville theorem of positive pseudo-harmonic functions on Sasakian manifolds with nonnegative pseudo-Hermitian Ricci curvature.


title: Synchronizing the Smallest Possible System

abstract: We investigate the minimal Hilbert-space dimension for a system to be synchronized. We first show that qubits cannot be synchronized due to the lack of a limit cycle. Moving to larger spin values, we demonstrate that a single spin 1 can be phase-locked to a weak external signal of similar frequency and exhibits all the standard features of the theory of synchronization. Our findings rely on the Husimi Q representation based on spin coherent states which we propose as a tool to obtain a phase portrait.


title: Are Two (Samples) Really Better Than One? On the Non-Asymptotic Performance of Empirical Revenue Maximization

abstract: The literature on "mechanism design from samples," which has flourished in recent years at the interface of economics and computer science, offers a bridge between the classic computer-science approach of worst-case analysis (corresponding to "no samples") and the classic economic approach of average-case analysis for a given Bayesian prior (conceptually corresponding to the number of samples tending to infinity). Nonetheless, the two directions studied so far are two extreme and almost diametrically opposed directions: that of asymptotic results where the number of samples grows large, and that where only a single sample is available. In this paper, we take a first step toward understanding the middle ground that bridges these two approaches: that of a fixed number of samples greater than one. In a variety of contexts, we ask what is possibly the most fundamental question in this direction: "are two samples really better than one sample?". We present a few surprising negative results, and complement them with our main result: showing that the worst-case, over all regular distributions, expected-revenue guarantee of the Empirical Revenue Maximization algorithm given two samples is greater than that of this algorithm given one sample. The proof is technically challenging, and provides the first result that shows that some deterministic mechanism constructed using two samples can guarantee more than one half of the optimal revenue.


title: A note on mass-conserving solutions to the coagulation-fragmentation equation by using non-conservative approximation

abstract: In general, the non-conservative approximation of coagulation-fragmentation equations (CFEs) may lead to the occurrence of gelation phenomenon. In this article, it is shown that the non-conservative approximation of CFEs can also provide the existence of mass conserving solutions to CFEs for large classes of unbounded coagulation and fragmentation kernels.


title: Comparative analysis of SVIT and Skype features in e-Learning process

abstract: The article discusses capabilities of the SVIT and Skype programs in terms of use in e-Learning process. Also various types of connections, image and sound quality were investigated.


title: Inverse analysis of traction-separation relationship based on sequentially linear approach

abstract: Traction-separation relationship is an important material characteristic describing the fracture behaviour of quasi-brittle solids. A new numerical scheme for identification of the traction-separation relation by inverse analysis of data obtained from various types of fracture tests is proposed. Due to employing the concept of sequentially linear analysis, the method exhibits a superior numerical stability and versatility. The applicability and effectiveness of the proposed method is demonstrated on examples involving identification of the traction-separation relationship using experimental data from various test configurations.


title: Charged pion masses under strong magnetic fields in the NJL model

abstract: The behavior of charged pion masses in the presence of a static uniform magnetic field is studied in the framework of the two-flavor NJL model, using a magnetic field-independent regularization scheme. Analytical calculations are carried out employing the Ritus eigenfunction method, which allows us to properly take into account the presence of Schwinger phases in the quark propagators. Numerical results are obtained for definite model parameters, comparing the predictions of the model with present lattice QCD results.


title: 2VRP: a benchmark problem for small but rich VRPs

abstract: We consider a 2-vehicle routing problem (2VRP) which can be viewed as a building block for the variety of vehicle routing problems (VRP). As a simplified version of the 2VRP, we consider a 2-period balanced travelling salesman problem (2TSP) and describe a polynomially solvable case of this NP-hard problem. For the 2VRP with general settings, we suggest a framework based on the Held and Karp dynamic programming algorithm. Our algorithms based on this framework show an exceptionally good performance on the published test data. Our approach can be easily extended to a variety of constraints/attributes in the VRP, hence the wording "small but rich" in the title of our paper. We also introduce a new methodological approach: we use easy solvable special cases for generating test instances and then use these instances in computational experiments.


title: Unrecognized Astrometric Confusion in the Galactic Centre

abstract: The Galactic Centre is a highly crowded stellar field and frequent unrecognized events of source confusion, which involve undetected faint stars, are expected to introduce astrometric noise on a sub-mas level. This confusion noise is the main non-instrumental effect limiting the astrometric accuracy and precision of current near-infrared imaging observations and the long-term monitoring of individual stellar orbits in the vicinity of the central supermassive black hole. We self-consistently simulate the motions of the known and the yet unidentified stars to characterize this noise component and show that a likely consequence of source confusion is a bias in estimates of the stellar orbital elements, as well as the inferred mass and distance of the black hole, in particular if stars are being observed at small projected separations from it, such as the star S2 during pericentre passage. Furthermore, we investigate modeling the effect of source confusion as an additional noise component that is time-correlated, demonstrating a need for improved noise models to obtain trustworthy estimates of the parameters of interest (and their uncertainties) in future astrometric studies.


title: Random Forest Classification of Stars in the Galactic Centre

abstract: Near-infrared high-angular resolution imaging observations of the Milky Way's nuclear star cluster have revealed all luminous members of the existing stellar population within the central parsec. Generally, these stars are either evolved late-type giants or massive young, early-type stars. We revisit the problem of stellar classification based on intermediate-band photometry in the K-band, with the primary aim of identifying faint early-type candidate stars in the extended vicinity of the central massive black hole. A random forest classifier, trained on a subsample of spectroscopically identified stars, performs similarly well as competitive methods (F1=0.85), without involving any model of stellar spectral energy distributions. Advantages of using such a machine-trained classifier are a minimum of required calibration effort, a predictive accuracy expected to improve as more training data becomes available, and the ease of application to future, larger data sets. By applying this classifier to archive data, we are also able to reproduce the results of previous studies of the spatial distribution and the K-band luminosity function of both the early- and late-type stars.


title: Spatially Resolved Thermodynamic Integration: An Efficient Method to Compute Chemical Potentials of Dense Fluids

abstract: Many popular methods for the calculation of chemical potentials rely on the insertion of test particles into the target system. In the case of liquids and liquid mixtures, this procedure increases in difficulty upon increasing density or concentration, and the use of sophisticated enhanced sampling techniques becomes inevitable. In this work we propose an alternative strategy, spatially resolved thermodynamic integration, or SPARTIAN for short. Here, molecules are described with atomistic resolution in a simulation subregion, and as ideal gas particles in a larger reservoir. All molecules are free to diffuse between subdomains adapting their resolution on the fly. To enforce a uniform density profile across the simulation box, a single-molecule external potential is computed, applied, and identified with the difference in chemical potential between the two resolutions. Since the reservoir is represented as an ideal gas bath, this difference exactly amounts to the excess chemical potential of the target system. The present approach surpasses the high density/concentration limitation of particle insertion methods because the ideal gas molecules entering the target system region spontaneously adapt to the local environment. The ideal gas representation contributes negligibly to the computational cost of the simulation, thus allowing one to make use of large reservoirs at minimal expenses. The method has been validated by computing excess chemical potentials for pure Lennard-Jones liquids and mixtures, SPC and SPC/E liquid water, and aqueous solutions of sodium chloride. The reported results well reproduce literature data for these systems.