title: Unveiling the structural arrangements responsible for the atomic dynamics in metallic glasses during physical aging

abstract: Understanding and controlling physical aging, i.e. the spontaneous temporal evolution of out-of-equilibrium systems, represents one of the greatest tasks in material science. Recent studies have revealed the existence of a complex atomic motion in metallic glasses, with different aging regimes in contrast with the typical continuous aging observed in macroscopic quantities. By combining dynamical and structural synchrotron techniques, for the first time we directly connect previously identified microscopic structural mechanisms with the peculiar atomic motion, providing a broader unique view of their complexity. We show that the atomic scale is dominated by the interplay between two processes: rearrangements releasing residual stresses related to a cascade mechanism of relaxation, and medium range ordering processes, which do not affect the local density, likely due to localized relaxations of liquid-like regions. As temperature increases, a surprising additional secondary relaxation process sets in, together with a faster medium range ordering, likely precursors of crystallization.


title: Solvation of the Morpholinium Cation in Acetonitrile. Effect of an Anion

abstract: Ionic liquids constitute a fast growing class of compounds finding multiple applications in science and technology. Morpholinium-based ionic liquids (MBILs) and their mixtures with polar molecular co-solvents are interesting as sustainable electrolyte systems for electrochemistry. We investigate local structures of protic and propton-free morpholinium cations in acetonitrile (ACN) using semi-empirical molecular dynamics (MD) simulations. An impact of an anion (acetate) on the cation solvation regularities is discussed. Unlike oxygen, nitrogen of the morpholine ring is a strong electrophilic binding center. This site is responsible for the interactions of the cation with the solvent and with the anion. In protic MBILs, the role of nitrogen is delegated to the proton, which is linked to nitrogen. The acetate anion weakens solvation of the cation due to occupation of space near nitrogen or proton. The analysis reveals a favorable solvation of MBILs in ACN, which is a prerequisite for a new high-performance electrolyte system. The reported structural data are validated through a point-to-point comparison with the MP2 post-Hartree-Fock theory.


title: An introduction to the analysis and implementation of sparse grid finite element methods

abstract: Our goal is to present an elementary approach to the analysis and programming of sparse grid finite element methods. This family of schemes can compute accurate solutions to partial differential equations, but using far fewer degrees of freedom than their classical counterparts. After a brief discussion of the classical Galerkin finite element method with bilinear elements, we give a short analysis of what is probably the simplest sparse grid method: the two-scale technique of Lin et al. (2001). We then demonstrate how to extend this to a multiscale sparse grid method which, up to choice of basis, is equivalent to the hierarchical approach, as described by, e.g., Bungartz and Griebel (2004). However, by presenting it as an extension of the two-scale method, we can give an elementary treatment of its analysis and implementation. For each method considered, we provide MATLAB code, and a comparison of accuracy and computational costs.


title: Double-Trace Deformations and Entanglement Entropy in AdS

abstract: We compute the bulk entanglement entropy of a massive scalar field in a Poincare AdS with the Dirichlet and Neumann boundary condition when we trace out a half space. Moreover, by taking into account the quantum back reaction to the minimal surface area, we calculate how much the entanglement entropy changes under a double-trace deformation of a holographic CFT. In the AdS3/CFT2 setup, our result agrees with the known result in 2d CFTs. In higher dimensions, our results offer holographic predictions.


title: A Comparison of Higher-Order Weak Numerical Schemes for Stopped Stochastic Differential Equations

abstract: We review, implement, and compare numerical integration schemes for spatially bounded diffusions stopped at the boundary which possess a convergence rate of the discretization error with respect to the timestep $h$ higher than ${\cal O}(\sqrt{h})$. We address specific implementation issues of the most general-purpose of such schemes. They have been coded into a single Matlab program and compared, according to their accuracy and computational cost, on a wide range of problems in up to ${\mathbb R}^{48}$. The paper is self-contained and the code will be made freely downloadable.


title: Sensitivity of the JEM-EUSO telescope to gravity effects in neutrino-induced air showers

abstract: We examine the JEM-EUSO sensitivity to gravity effects in the context of Randall-Sundrum (RS) model with a single extra dimension and small curvature of the metric. Exchanges of reggeized Kaluza-Klein gravitons in the $t$-channel contribute to the inelastic cross-section for scattering of ultra-high-energy neutrinos off nucleons. Such effects can be detected in deeply penetrating quasi-horizontal air showers induced by interactions of cosmic neutrinos with atmospheric nucleons. For this reason, we calculate the expected number of quasi-horizontal air showers at the JEM-EUSO observatory as a function of two parameters of the RS model.


title: The Sagnac Effect in Optical Lattices with Laser-Assisted Tunneling

abstract: We propose a scheme to realize the rotation sensing using optical lattices with laser-assisted tunneling. We demonstrate that the competition between the rotation and the spin-orbit coupling governs the spin-dependent response of the cyclotron dynamics of the spin-orbit coupled bosons. The Sagnac-type cumulative phase can be read out from the envelop of a beat-frequency particle current in the spin-balanced system and enhanced by the cyclotron motion. We also theoretically show that the sensitivity limit of the spin-orbit-coupled system to the rotational motion can reach 7*10^-8 rad s^-1 Hz^1/2.


title: Hydrodynamical noise and Gubser flow

abstract: Hydrodynamical noise is introduced on top of Gubser's analytical solution to viscous hydrodynamics. With respect to the ultra-central collision events of Pb-Pb, p-Pb and p-p at the LHC energies, we solve the evolution of noisy fluid systems and calculate the radial flow velocity correlations. We show that the absolute amplitude of the hydrodynamical noise is determined by the multiplicity of the collision event. The evolution of azimuthal anisotropies, which is related to the generation of harmonic flow, receives finite enhancements from hydrodynamical noise. Although it is strongest in the p-p systems, the effect of hydrodynamical noise on flow harmonics is found to be negligible, especially in the ultra-central Pb-Pb collisions. For the short-range correlations, hydrodynamical noise contributes to the formation of a near-side peak on top of the correlation structure originated from initial state fluctuations. The shape of the peak is affected by the strength of hydrodynamical noise, whose height and width grow from the Pb-Pb system to the p-Pb and p-p systems.


title: On the Grain-Modified Magnetic Diffusivities in Protoplanetary Disks

abstract: Weakly ionized protoplanetary disks (PPDs) are subject to non-ideal-magnetohydrodynamic (MHD) effects including Ohmic resistivity, the Hall effect and ambipolar diffusion (AD), and the resulting magnetic diffusivities ($\eta_O, \eta_H$ and $\eta_A$) largely control the disk gas dynamics. The presence of grains not only strongly reduces disk ionization fraction, but also modify the scalings of $\eta_H$ and $\eta_A$ with magnetic field strength. We derive analytically asymptotic expressions of $\eta_H$ and $\eta_A$ in both strong and weak field limits and show that towards strong field, $\eta_H$ can change sign (at a threshold field strength $B_{\rm th}$), mimicking a flip of field polarity, and AD is substantially reduced. Applying to PPDs, we find that when small $\sim0.1$ ($0.01$)$\mu$m grains are sufficiently abundant [mass ratio $\sim0.01$ ($10^{-4}$)], $\eta_H$ can change sign up to $\sim2-3$ scale heights above midplane at modest field strength (plasma $\beta\sim100$) over a wide range of disk radii. Reduction of AD is also substantial towards the AD dominated outer disk and may activate the magneto-rotational instability. We further perform local non-ideal MHD simulations of the inner disk (within 10 AU) and show that with sufficiently abundant small grains, magnetic field amplification due to the Hall-shear instability saturates at very low level near the threshold field strength $B_{\rm th}$. Together with previous studies, we conclude by discussing the grain-abundance-dependent phenomenology of PPD gas dynamics.


title: Limit cycles bifurcating from a period annulus in continuous piecewise differential systems with three zones

abstract: We study a class of planar continuous piecewise linear vector fields with three zones. Using the Poincar\'e map and some techniques for proving the existence of limit cycles for smooth differential systems, we prove that this class admits at least two limit cycles that appear by perturbations of a period annulus. Moreover, we describe the bifurcation of the limit cycles for this class through two examples of two-parameter families of piecewise linear vector fields with three zones.


title: Aligned natural inflation with modulations

abstract: The weak gravity conjecture applied for the aligned natural inflation indicates that generically there can be a modulation of the inflaton potential, with a period determined by sub-Planckian axion scale. We study the oscillations in the primordial power spectrum induced by such modulation, and discuss the resulting observational constraints on the model.


title: Complete positivity, finite temperature effects, and additivity of noise for time-local qubit dynamics

abstract: We present a general model of qubit dynamics which entails pure dephasing and dissipative time-local master equations. This allows us to describe the combined effect of thermalisation and dephasing beyond the usual Markovian approximation. We investigate the complete positivity conditions and introduce a heuristic model that is always physical and provides the correct Markovian limit. We study the effects of temperature on the non-Markovian behaviour of the system and show that the noise additivity property discussed by Yu and Eberly in Ref. [1] holds beyond the Markovian limit.


title: Some Dynamic Market Models

abstract: In this text, we study the temporal behavior of markets using models expressible as ordinary differential equations. The markets studied are those where each customer buys only one copy of the good, for example, subscription of smartphone service, journals and newspapers, and goods such as books, music and games. The underlying model is the diffusion model of Frank Bass. Evolution of markets with no competitors and markets with several competitors are analyzed where, in particulat, the effects of churning upon the market evolution is investigated. Analytic solutions are given for the temporal evolution of several types of interactive games.


title: H\"older-type inequalities and their applications to concentration and correlation bounds

abstract: Let $Y_v, v\in V,$ be $[0,1]$-valued random variables having a dependency graph $G=(V,E)$. We show that \[ \mathbb{E}\left[\prod_{v\in V} Y_{v} \right] \leq \prod_{v\in V} \left\{ \mathbb{E}\left[Y_v^{\frac{\chi_b}{b}}\right] \right\}^{\frac{b}{\chi_b}}, \] where $\chi_b$ is the $b$-fold chromatic number of $G$. This inequality may be seen as a dependency-graph analogue of a generalised H\"older inequality, due to Helmut Finner. Additionally, we provide applications of H\"older-type inequalities to concentration and correlation bounds for sums of weakly dependent random variables.


title: Parameter estimation in atomic spectroscopy using exceptional points

abstract: The dynamics of spontaneous emission of an atomic system is studied in the framework of an open quantum system. The resulting quantum master equation for the atomic system is non hermitian. The generator $\mathcal{L}$ can possess non-hermitian degeneracies, i.e. exceptional points ($EP$), for specific values of the external laser driving amplitude and detuning from the atomic lines. We suggest to employ the special properties of these $EP$s for accurate parameter estimation. The method is demonstrated for the atomic spectrum of $S \rightarrow P$ transitions of $~^{85}$Rb and $~^{40}$Ca$~^+$.


title: Excitation of an Atomic Transition with a Vortex Laser Beam

abstract: Photons carry one unit of angular momentum associated with their spin~\cite{Beth1936}. Structured vortex beams carry additional orbital angular momentum which can also be transferred to matter~\cite{Allen1992}. This extra twist has been used for example to drive motion of microscopic particles in optical tweezers as well as to create vortices in degenerate quantum gases~\cite{He1995,Andersen2006}. Here we demonstrate the transfer of optical orbital angular momentum from the transverse spatial structure of the beam to the internal (electronic) degrees of freedom of an atom. Probing a quadrupole transition of a single trapped $^{40}$Ca$^+$ ion localized at the center of the vortex, we observe strongly modified selection rules, accounting for both the photon spin and the vorticity of the field. In particular, we show that an atom can absorb two quanta of angular momentum from a single photon even when rotational symmetry is conserved. In contrast to previous findings~\cite{Araoka2005,Loeffler2011a,Mathevet2013}, our experiment allows for conditions where the vorticity of the laser beam determines the optical excitation, contributing to the long-standing discussion on whether the orbital angular momentum of photons can be transferred to atomic internal degrees of freedom~\cite{VanEnk1994,Babiker2002,Jauregui2004, Schmiegelow2012, Mondal2014, Scholz-Marggraf2014} and paves the way for its use to tailor light-matter interactions.


title: Improving the performance of the linear systems solvers using CUDA

abstract: Parallel computing can offer an enormous advantage regarding the performance for very large applications in almost any field: scientific computing, computer vision, databases, data mining, and economics. GPUs are high performance many-core processors that can obtain very high FLOP rates. Since the first idea of using GPU for general purpose computing, things have evolved and now there are several approaches to GPU programming: CUDA from NVIDIA and Stream from AMD. CUDA is now a popular programming model for general purpose computations on GPU for C/C++ programmers. A great number of applications were ported to CUDA programming model and they obtain speedups of orders of magnitude comparing to optimized CPU implementations. In this paper we present an implementation of a library for solving linear systems using the CCUDA framework. We present the results of performance tests and show that using GPU one can obtain speedups of about of approximately 80 times comparing with a CPU implementation.


title: A Generalization of a Bin-Based Modification to the Stochastic Simulation Algorithm

abstract: In a 1996 paper, See$\beta$elberg, Trautmann and Thorn modified Gillespie's (1975) Monte Carlo algorithm which is used to stochastically simulate the collision and coalescence process. Their modification reduces the storage requirements of the simulation by several orders of magnitude. However, their modification creates unphysical and potentially fatal conditions when used with common initial distributions. We identify those conditions and propose a solution to maintain physically real conditions for all state variables during the evolution of any initial distribution.


title: Multi-Agent Continuous Transportation with Online Balanced Partitioning

abstract: We introduce the concept of continuous transportation task to the context of multi-agent systems. A continuous transportation task is one in which a multi-agent team visits a number of fixed locations, picks up objects, and delivers them to a final destination. The goal is to maximize the rate of transportation while the objects are replenished over time. Examples of problems that need continuous transportation are foraging, area sweeping, and first/last mile problem. Previous approaches typically neglect the interference and are highly dependent on communications among agents. Some also incorporate an additional reconnaissance agent to gather information. In this paper, we present a hybrid of centralized and distributed approaches that minimize the interference and communications in the multi-agent team without the need for a reconnaissance agent. We contribute two partitioning-transportation algorithms inspired by existing algorithms, and contribute one novel online partitioning-transportation algorithm with information gathering in the multi-agent team. Our algorithms have been implemented and tested extensively in the simulation. The results presented in this paper demonstrate the effectiveness of our algorithms that outperform the existing algorithms, even without any communications between the agents and without the presence of a reconnaissance agent.


title: Nearest Neighbor search in Complex Network for Community Detection

abstract: Nearest neighbor search is a basic computational tool used extensively in almost research domains of computer science specially when dealing with large amount of data. However, the use of nearest neighbor search is restricted for the purpose of algorithmic development by the existence of the notion of nearness among the data points. The recent trend of research is on large, complex networks and their structural analysis, where nodes represent entities and edges represent any kind of relation between entities. Community detection in complex network is an important problem of much interest. In general, a community detection algorithm represents an objective function and captures the communities by optimizing it to extract the interesting communities for the user. In this article, we have studied the nearest neighbor search problem in complex network via the development of a suitable notion of nearness. Initially, we have studied and analyzed the exact nearest neighbor search using metric tree on proposed metric space constructed from complex network. After, the approximate nearest neighbor search problem is studied using locality sensitive hashing. For evaluation of the proposed nearest neighbor search on complex network we applied it in community detection problem. The results obtained using our methods are very competitive with most of the well known algorithms exists in the literature and this is verified on collection of real networks. On the other-hand, it can be observed that time taken by our algorithm is quite less compared to popular methods.


title: Noisy Submodular Maximization via Adaptive Sampling with Applications to Crowdsourced Image Collection Summarization

abstract: We address the problem of maximizing an unknown submodular function that can only be accessed via noisy evaluations. Our work is motivated by the task of summarizing content, e.g., image collections, by leveraging users' feedback in form of clicks or ratings. For summarization tasks with the goal of maximizing coverage and diversity, submodular set functions are a natural choice. When the underlying submodular function is unknown, users' feedback can provide noisy evaluations of the function that we seek to maximize. We provide a generic algorithm -- \submM{} -- for maximizing an unknown submodular function under cardinality constraints. This algorithm makes use of a novel exploration module -- \blbox{} -- that proposes good elements based on adaptively sampling noisy function evaluations. \blbox{} is able to accommodate different kinds of observation models such as value queries and pairwise comparisons. We provide PAC-style guarantees on the quality and sampling cost of the solution obtained by \submM{}. We demonstrate the effectiveness of our approach in an interactive, crowdsourced image collection summarization application.


title: Face Alignment Across Large Poses: A 3D Solution

abstract: Face alignment, which fits a face model to an image and extracts the semantic meanings of facial pixels, has been an important topic in CV community. However, most algorithms are designed for faces in small to medium poses (below 45 degree), lacking the ability to align faces in large poses up to 90 degree. The challenges are three-fold: Firstly, the commonly used landmark-based face model assumes that all the landmarks are visible and is therefore not suitable for profile views. Secondly, the face appearance varies more dramatically across large poses, ranging from frontal view to profile view. Thirdly, labelling landmarks in large poses is extremely challenging since the invisible landmarks have to be guessed. In this paper, we propose a solution to the three problems in an new alignment framework, called 3D Dense Face Alignment (3DDFA), in which a dense 3D face model is fitted to the image via convolutional neutral network (CNN). We also propose a method to synthesize large-scale training samples in profile views to solve the third problem of data labelling. Experiments on the challenging AFLW database show that our approach achieves significant improvements over state-of-the-art methods.


title: Center-to-limb variation of intensity and polarization in continuum spectra of FGK stars for spherical atmospheres

abstract: One of the necessary parameters needed for the interpretation of the light curves of transiting exoplanets or eclipsing binaries, as well as interferometric measurements of a star or microlensing events is how the intensity and polarization of a light change from the center to the limb. Scattering and absorption processes in stellar atmosphere affect both the center-to limb variation of intensity (CLVI) and polarization (CLVP). In this paper, we present a study of the CLVI and CLVP in continuum spectra considering different contributions of scattering and absorption opacity for different spectral type stars with spherical atmospheres. We solve the polarized radiative transfer equation in the presence of continuum scattering, considering spherical stellar model atmospheres. We developed two independent codes based on Feautrier and short characteristics methods to cross-check our results. We calculate the CLVI and CLVP in continuum for the Phoenix grid of spherical stellar model atmospheres for a range of $T_{eff} = 4000 - 7000 \rm K$, $\log g = 1.0 - 5.5$ and $\lambda = 4000 - 7000 \rm \AA$, which are tabulated and available at the CDS. For sub-giant and dwarf stars ($\log g = 3.0 - 4.5$), lower $\log g$ and lower $T_{eff}$ of a star lead to higher limb polarization of the star. For giant and supergiant stars ($\log g = 1.0 - 2.5$), the highest effective temperature yields the largest polarization. By decreasing of the $T_{eff}$ of a star down to $4500 - 5500 \rm K$ (depending on $\log g$) the limb polarization decreases and reaches a local minimum. It increases again down to $T_{eff}$ of $4000 \rm K$. For the most compact dwarf stars ($\log g = 5.0 - 5.5$) the limb polarization degree shows a maximum for models with $T_{eff}$ in the range $4200 - 4600 \rm K$ (depending on $\log g$) and decreases toward higher and lower temperatures.


title: The ambient obstruction tensor and conformal holonomy

abstract: For a conformal manifold, we describe a new relation between the ambient obstruction tensor of Fefferman and Graham and the holonomy of the normal conformal Cartan connection. This relation allows us to prove several results on the vanishing and the rank of the obstruction tensor, for example for conformal structures admitting twistor spinors or normal conformal Killing forms. As our main tool we introduce the notion of a conformal holonomy distribution and show that its integrability is closely related to the exceptional conformal structures in dimensions five and six that were found by Nurowski and Bryant.


title: Atomtronics

abstract: The progress achieved in micro-fabricating potential for cold atoms has defined a new field in quantum technology - Atomtronics - where a variety of 'atom circuits' of very different spatial shapes and depth have been devised for atom manipulation, with a precision that nowadays is approaching that of lithographic techniques. Atomtronic setups are characterized by enhanced flexibility and control of the fundamental mechanisms underlying their functionalities and by the reduced decoherence rate that is typical of cold-atom systems. Such an approach is expected to be instrumental for the realization of quantum devices of a radically new type and, at the same time, to enlarge the scope of cold atom quantum simulators. In this article we give a short overview of the field and draw a roadmap for potential future directions.


title: Testing Hadronic Interaction Models using a Highly Granular Silicon-Tungsten Calorimeter

abstract: The CALICE collaboration has published a detailed study of hadronic interactions using data recorded with the highly granular CALICE silicon-tungsten electromagnetic calorimeter (Si-W ECAL). Approximately 350,000 selected negative pion events at energies between 2 and 10 GeV have been studied. The predictions of several physics models available within the Geant4 simulation tool kit are compared to this data. A reasonable overall description of the data is observed; the Monte Carlo predictions are within 20 % of the data, and for many observables much closer. The largest quantitative discrepancies are found in the longitudinal and transverse distributions of the reconstructed energy. Based on the good control of the data set and general observables, the next step is to achieve a deeper understanding of hadronic interactions by studying the interaction zone and by reconstructing secondaries that emerge from the hadronic interaction in the Si-W ECAL.


title: Positive Discrete Spectrum of the Evolutionary Operator of Supercritical Branching Walks with Heavy Tails

abstract: We consider a continuous-time symmetric supercritical branching random walk on a multidimensional lattice with a finite set of the particle generation centres, i.e. branching sources. The main object of study is the evolutionary operator for the mean number of particles both at an arbitrary point and on the entire lattice. The existence of positive eigenvalues in the spectrum of an evolutionary operator results in an exponential growth of the number of particles in branching random walks, called supercritical in the such case. For supercritical branching random walks, it is shown that the amount of positive eigenvalues of the evolutionary operator, counting their multiplicity, does not exceed the amount of branching sources on the lattice, while the maximal of these eigenvalues is always simple. We demonstrate that the appearance of multiple lower eigenvalues in the spectrum of the evolutionary operator can be caused by a kind of `symmetry' in the spatial configuration of branching sources. The presented results are based on Green's function representation of transition probabilities of an underlying random walk and cover not only the case of the finite variance of jumps but also a less studied case of infinite variance of jumps.


title: Convex Optimization Based State Estimation against Sparse Integrity Attacks

abstract: We consider the problem of robust estimation in the presence of integrity attacks. There are m sensors monitoring the state and p of them are under attack. The malicious measurements collected by the compromised sensors can be manipulated arbitrarily by the attacker. The classical estimators such as the least squares estimator may not provide a reliable estimate under the so-called (p,m)-sparse attack. In this work, we are not restricting our efforts in studying whether any specific estimator is resilient to the attack or not, but instead we aim to present some generic sufficient and necessary conditions for robustness by considering a general class of convex optimization based estimators. The sufficient and necessary conditions are shown to be tight, with a trivial gap.


title: Excited Baryon Structure Using Exclusive Reactions with CLAS12

abstract: Studying excited nucleon structure through exclusive electroproduction reactions is an important avenue for exploring the nature of the non-perturbative strong interaction. Electrocouplings for $N^*$ states in the mass range below 1.8~GeV have been determined from analyses of CLAS $\pi N$, $\eta N$, and $\pi \pi N$ data. This work made it clear that consistency of independent analyses of exclusive channels with different couplings and non-resonant backgrounds but the same $N^*$ electro-excitation amplitudes, is essential to have confidence in the extracted results. In terms of hadronic coupling, many high-lying $N^*$ states preferentially decay through the $\pi \pi N$ channel instead of $\pi N$. Data from the $KY$ channels will therefore be critical to provide an independent analysis to compare the extracted electrocouplings for the high-lying $N^*$ states against those determined from the $\pi N$ and $\pi \pi N$ channels. A program to study excited $N^*$ decays to non-strange and strange exclusive final states using CLAS12 will measure differential cross sections to be used as input to extract the $\gamma_vNN^*$ transition form factors for the most prominent $N^*$ states in the range of invariant energy $W$ up 3~GeV in the virtually unexplored domain of momentum transfers $Q^2$ up to 12~GeV$^2$.


title: Different kinds of chimera death states in nonlocally coupled oscillators

abstract: We investigate the significance of nonisochronicity parameter in a network of nonlocally coupled Stuart-Landau oscillators with symmetry breaking form. We observe that the presence of nonisochronicity parameter leads to structural changes in the chimera death region while varying the strength of the interaction. This gives rise to the existence of different types of chimera death states such as multi-chimera death state, type-I periodic chimera death (PCD) state and, type-II periodic chimera death (PCD) state. We also find that the number of periodic domains in both the types of PCD states exponentially decreases with an increase of coupling range and obeys a power law under nonlocal coupling. Additionally, we also analyze the structural changes of chimera death states by reducing the system of dynamical equations to a phase model through the phase approximation. We also briefly study the role of nonisochronicity parameter on chimera states, where the existence of multi-chimera state with respect to the coupling range is pointed out. Moreover, we also analyze the robustness of the chimera death state to perturbations in the natural frequencies of the oscillators.


title: Supercongruences on some binomial sums involving Lucas sequences

abstract: In this paper, we confirm several conjectured congruences of Sun concerning the divisibility of binomial sums. For example, with help of a quadratic hypergeometric transformation, we prove that $$ \sum_{k=0}^{p-1}\binom{p-1}k\binom{2k}k^2\frac{P_k}{8^k}\equiv0\pmod{p^2} $$ for any prime $p\equiv 7\mod{8}$, where $P_k$ is the $k$-th Pell number. Further, we also propose three new congruences of the same type.


title: A predictive structural model of the primate connectome

abstract: Anatomical connectivity imposes strong constraints on brain function, but there is no general agreement about principles that govern its organization. Based on extensive quantitative data we tested the power of three models to predict connections of the primate cerebral cortex: architectonic similarity (structural model), spatial proximity (distance model) and thickness similarity (thickness model). Architectonic similarity showed the strongest and most consistent influence on connection features. This parameter was strongly associated with the presence or absence of inter-areal connections and when integrated with spatial distance, the model allowed predicting the existence of projections with very high accuracy. Moreover, architectonic similarity was strongly related to the laminar pattern of projections origins, and the absolute number of cortical connections of an area. By contrast, cortical thickness similarity and distance were not systematically related to connection features. These findings suggest that cortical architecture provides a general organizing principle for connections in the primate brain.


title: Towards analytical model optimization in atmospheric tomography

abstract: Modern ground-based telescopes rely on a technology called adaptive optics (AO) in order to compensate for the loss of image quality caused by atmospheric turbulence. Next-generation AO systems designed for a wide field of view require a stable and high-resolution reconstruction of the refractive index fluctuations in the atmosphere. By introducing a novel Bayesian method, we address the problem of estimating an atmospheric turbulence strength profile and reconstructing the refractive index fluctuations simultaneously, where we only use wavefront measurements of incoming light from guide stars. Most importantly, we demonstrate how this method can be used for model optimization as well. We propose two different algorithms for solving the maximum a posteriori estimate: the first approach is based on alternating minimization and has the advantage of integrability into existing atmospheric tomography methods. In the second approach, we formulate a convex non-differentiable optimization problem, which is solved by an iterative thresholding method. This approach clearly illustrates the underlying sparsity-enforcing mechanism for the strength profile. By introducing a tuning/regularization parameter, an automated model reduction of the layer structure of the atmosphere is achieved. Using numerical simulations, we demonstrate the performance of our method in practice.


title: Tensor-optimized antisymmetrized molecular dynamics in nuclear physics

abstract: We develop a new formalism to treat nuclear many-body systems using bare nucleon-nucleon interaction. It has become evident that the tensor interaction plays important role in nuclear many-body systems due to the role of the pion in strongly interacting system. We take the antisymmetrized molecular dynamics (AMD) as a basic framework and add a tensor correlation operator acting on the AMD wave function using the concept of the tensor-optimized shell model (TOSM). We demonstrate a systematical and straightforward formulation utilizing the Gaussian integration and differentiation method and the antisymmetrization technique to calculate all the matrix elements of the many-body Hamiltonian. We can include the three-body interaction naturally and calculate the matrix elements systematically in the progressive order of the tensor correlation operator. We call the new formalism "tensor-optimized antisymmetrized molecular dynamics".


title: Critical current density and trapped field in HTS with asymmetric magnetization loops

abstract: Applications of the extended critical state model are considered. The trapped magnetic field, the penetration field and the field dependence of the critical current density are analysed. The critical current density and the trapped field in superconducting grains depend on the grain size. Asymmetry of the hysteresis curves relative to the M = 0 axis is related to the scale of the current circulation.


title: Pipelined, Flexible Krylov Subspace Methods

abstract: We present variants of the Conjugate Gradient (CG), Conjugate Residual (CR), and Generalized Minimal Residual (GMRES) methods which are both pipelined and flexible. These allow computation of inner products and norms to be overlapped with operator and nonlinear or nondeterministic preconditioner application.The methods are hence aimed at hiding network latencies and synchronizations which can become computational bottlenecks in Krylov methods on extreme-scale systems or in the strong-scaling limit. The new variants are not arithmetically equivalent to their base flexible Krylov methods, but are chosen to be similarly performant in a realistic use case, the application of strong nonlinear preconditioners to large problems which require many Krylov iterations. We provide scalable implementations of our methods as contributions to the PETSc package and demonstrate their effectiveness with practical examples derived from models of mantle convection and lithospheric dynamics with heterogeneous viscosity structure. These represent challenging problems where multiscale nonlinear preconditioners are required for the current state-of-the-art algorithms, and are hence amenable to acceleration with our new techniques. Large-scale tests are performed in the strong-scaling regime on a contemporary leadership supercomputer, where speedups approaching, and even exceeding $2\times$ can be observed. We conclude by analyzing our new methods with a performance model targeted at future exascale machines.


title: Molecular Gas and Star Formation in the Cartwheel

abstract: Atacama Large Millimeter/submillimeter Array (ALMA) 12CO(J=1-0) observations are used to study the cold molecular ISM of the Cartwheel ring galaxy and its relation to HI and massive star formation (SF). CO moment maps find $(2.69\pm0.05)\times10^{9}$ M$_{\odot}$ of H$_2$ associated with the inner ring (72%) and nucleus (28%) for a Galactic I(CO)-to-N(H2) conversion factor ($\alpha_{\rm CO}$). The spokes and disk are not detected. Analysis of the inner ring's CO kinematics show it to be expanding ($V_{\rm exp}=68.9\pm4.9$ km s$^{-1}$) implying an $\approx70$ Myr age. Stack averaging reveals CO emission in the starburst outer ring for the first time, but only where HI surface density ($\Sigma_{\rm HI}$) is high, representing $M_{\rm H_2}=(7.5\pm0.8)\times10^{8}$ M$_{\odot}$ for a metallicity appropriate $\alpha_{\rm CO}$, giving small $\Sigma_{\rm H_2}$ ($3.7$ M$_{\odot}$ pc$^{-2}$), molecular fraction ($f_{\rm mol}=0.10$), and H$_2$ depletion timescales ($\tau_{\rm mol} \approx50-600$ Myr). Elsewhere in the outer ring $\Sigma_{\rm H_2}\lesssim 2$ M$_{\odot}$ pc$^{-2}$, $f_{\rm mol}\lesssim 0.1$ and $\tau_{\rm mol}\lesssim 140-540$ Myr (all $3\sigma$). The inner ring and nucleus are H$_2$-dominated and are consistent with local spiral SF laws. $\Sigma_{\rm SFR}$ in the outer ring appears independent of $\Sigma_{\rm H_2}$, $\Sigma_{\rm HI}$ or $\Sigma_{\rm HI+H_2}$. The ISM's long confinement in the robustly star forming rings of the Cartwheel and AM0644-741 may result in either a large diffuse H$_2$ component or an abundance of CO-faint low column density molecular clouds. The H$_2$ content of evolved starburst rings may therefore be substantially larger. Due to its lower $\Sigma_{\rm SFR}$ and age the Cartwheel's inner ring has yet to reach this state. Alternately, the outer ring may trigger efficient SF in an HI-dominated ISM.


title: Sensitivity Study of Searching for $\tau^- \to \gamma \mu^-$ at HIEPA

abstract: The charged lepton flavor violation process is a clean and sensitive probe of new physics beyond the Standard Model. A sensitivity study is performed to the process $\tau^- \to \gamma \mu^-$ based on a 3~fb$^{-1}$ inclusive Monte Carlo sample of $e^+e^-$ collisions at a center-of-mass energy of 4.26 or 4.6~GeV, in the framework of the BESIII software system. The 90\% confidence level upper limits on $\BR(\tau^- \to \gamma \mu^-)$ are estimated assuming no signal is produced. We also obtain the sensitivity on $\BR(\tau^- \to \gamma \mu^-)$ as a function of the integrated luminosity, to serve as a reference for the HIEPA being proposed in China. It is found that 6.34~ab$^{-1}$ are needed to reach the current best upper limit of $4.4\times 10^{-8}$ and about 2510~ab$^{-1}$ are needed to reach a sensitivity of $10^{-9}$ if the detector design is similar to that of BESIII.


title: Leading and following sunspots: their magnetic properties and ultra-violet emission above them

abstract: Using SDO/HMI and SDO/AIA data for sunspot groups of the 24th solar cycle, we analyzed magnetic properties and He II 304 emission in leading and following sunspots separately. Simultaneous examination of umbral magnetic properties and atmospheric characteristics above the umbrae draws on average differences in He II 304 contrast over the umbrae of leading and following spots we discovered earlier for solar cycle 23 sunspot groups based on SOHO data as well as on the hypothetical relationship between contrast asymmetry and magnetic field asymmetry in umbrae. We use a more accurate and faster algorithm for solving the pi-uncertainty problem of the transverse magnetic field direction in this research producing new results on differences in magnetic field properties between magneto-conjugated leaders and followers. We found that, in ~78% of the cases, the minimum (over the umbra area) angle between the magnetic field line and the normal to the solar surface, a_min, is smaller in the leading spots, so the magnetic field there is more vertical than that in the counterpart following spot. It was also found that umbral area-averaged angle <a> in ~83% of the spot groups examined is smaller in the leader compared to the follower and the maximum and mean magnetic flux densities inside the umbra depend on the umbral area.


title: Some Results on Skorokhod Embedding and Robust Hedging with Local Time

abstract: In this paper, we provide some results on Skorokhod embedding with local time and its applications to the robust hedging problem in finance. First we investigate the robust hedging of options depending on the local time by using the recently introduced stochastic control approach, in order to identify the optimal hedging strategies, as well as the market models that realize the extremal no-arbitrage prices. As a by-product, the optimality of Vallois' Skorokhod embeddings is recovered. In addition, under appropriate conditions, we derive a new solution to the two-marginal Skorokhod embedding as a generalization of the Vallois solution. It turns out from our analysis that one needs to relax the monotonicity assumption on the embedding functions in order to embed a larger class of marginal distributions. Finally, in a full-marginal setting where the stopping times given by Vallois are well-ordered, we construct a remarkable Markov martingale which provides a new example of fake Brownian motion.


title: Remarks on the nonexistence of biharmonic maps

abstract: In this short note we study nonexistence result of biharmonic maps from a complete Riemannian manifold into a Riemannian manifold with nonpositive sectional curvature. Assume that $\phi:(M,g)\to (N, h)$ is a biharmonic map, where $(M, g)$ is a complete Riemannian manifold and $(N,h)$ a Riemannian manifold with nonpositive sectional curvature, we will prove that $\phi$ is a harmonic map if one of the following conditions holds: (i) $|d\phi|$ is bounded in $L^q(M)$ and $ \int_M|\tau(\phi)|^pdv_g<\infty, $ for some $1\leq q\leq\infty$, $1< p<\infty$; or (ii) $Vol(M)=\infty$ and $ \int_M|\tau(\phi)|^pdv_g<\infty, $ for some $1< p<\infty$. In addition if $N$ has negative sectional curvature, we assume that $rank\phi(q)\geq2$ for some $q\in M$ and $\int_M|\tau(\phi)|^pdv_g<\infty, $ for some $1< p<\infty$. These results improve the related theorems due to Baird et al.(cf. \cite{BFO}), Nakauchi et al.(cf. \cite{NUG}), Maeta(cf. \cite{Ma}) and Luo(cf. \cite{Luo}).


title: Chern mosaic - topology of chiral superconductivity on ferromagnetic adatom lattices

abstract: In this work we will explore the properties of superconducting surfaces decorated by two-dimensional ferromagnetic adatom lattices. As discovered recently [R\"ontynen and Ojanen, Phys. Rev. Lett. \textbf{114}, 236803 (2015)], in the presence of a Rashba spin-orbit coupling these systems may support topological superconductivity with complex phase diagrams and high Chern numbers. We show how the long-range hopping nature of the effective low-energy theory generically gives rise to a phase diagram covered by a \emph{Chern mosaic} -- a rich pattern of topological phases with large Chern numbers. We study different lattice geometries and the dependence of energy gaps and abundance of different phases as a function of system parameters. Our findings establish the studied system as one of the richest platforms for topological matter known to date.


title: Construction of Unit-Memory MDS Convolutional Codes

abstract: Maximum-distance separable (MDS) convolutional codes form an optimal family of convolutional codes, the study of which is of great importance. There are very few general algebraic constructions of MDS convolutional codes. In this paper, we construct a large family of unit-memory MDS convolutional codes over $\F$ with flexible parameters. Compared with previous works, the field size $q$ required to define these codes is much smaller. The construction also leads to many new strongly-MDS convolutional codes, an important subclass of MDS convolutional codes proposed and studied in \cite{GL2}. Many examples are presented at the end of the paper.


title: The semiclassical propagator in Fock space: dynamical echo and many-body interference

abstract: We present a semiclassical approach to many-body quantum propagation in terms of coherent sums over quantum amplitudes associated with the solutions of corresponding classical nonlinear wave equations. This approach adequately describes interference effects in the many-body space of interacting bosonic systems. The main quantity of interest, the transition amplitude between Fock states when the dynamics is driven by both single-particlecontributions and many-body interactions of similar magnitude, is non-perturbatively constructed in the spirit of Gutzwiller's derivation of the van Vleck propagator from the path integral representation of the time evolution operator, but lifted to the space of symmetrized many-body states. Effects beyond mean-field, here representing the classical limit of the theory, are semiclassically described by means of interfering amplitudes where the action and stability of the classical solutions enter. In this way, a genuinely many-body echo phenomenon, coherent backscattering in Fock space, is presented arising due to coherent quantum interference between classical solutions related by time reversal.


title: On the well-posedness of the Holm-Staley b-family of equations

abstract: In this paper we consider the Holm-Staley $b$-family of equations in the Sobolev spaces $H^s(\mathbb R)$ for $s > 3/2$. Using a geometric approach we show that, for any value of the parameter $b$, the corresponding solution map,$u(0) \mapsto u(T)$, is nowhere locally uniformly continuous.


title: Does Gaussian Approximation Work Well for The Long-Length Polar Code Construction?

abstract: Gaussian approximation (GA) is widely used to construct polar codes. However when the code length is long, the subchannel selection inaccuracy due to the calculation error of conventional approximate GA (AGA), which uses a two-segment approximation function, results in a catastrophic performance loss. In this paper, new principles to design the GA approximation functions for polar codes are proposed. First, we introduce the concepts of polarization violation set (PVS) and polarization reversal set (PRS) to explain the essential reasons that the conventional AGA scheme cannot work well for the long-length polar code construction. In fact, these two sets will lead to the rank error of subsequent subchannels, which means the orders of subchannels are misaligned, which is a severe problem for polar code construction. Second, we propose a new metric, named cumulative-logarithmic error (CLE), to quantitatively evaluate the remainder approximation error of AGA in logarithm. We derive the upper bound of CLE to simplify its calculation. Finally, guided by PVS, PRS and CLE bound analysis, we propose new construction rules based on a multi-segment approximation function, which obviously improve the calculation accuracy of AGA so as to ensure the excellent performance of polar codes especially for the long code lengths. Numerical and simulation results indicate that the proposed AGA schemes are critical to construct the high-performance polar codes.


title: Predicting Relevance based on Assessor Disagreement: Analysis and Practical Applications for Search Evaluation

abstract: Evaluation of search engines relies on assessments of search results for selected test queries, from which we would ideally like to draw conclusions in terms of relevance of the results for general (e.g., future, unknown) users. In practice however, most evaluation scenarios only allow us to conclusively determine the relevance towards the particular assessor that provided the judgments. A factor that cannot be ignored when extending conclusions made from assessors towards users, is the possible disagreement on relevance, assuming that a single gold truth label does not exist. This paper presents and analyzes the Predicted Relevance Model (PRM), which allows predicting a particular result's relevance for a random user, based on an observed assessment and knowledge on the average disagreement between assessors. With the PRM, existing evaluation metrics designed to measure binary assessor relevance, can be transformed into more robust and effectively graded measures that evaluate relevance towards a random user. It also leads to a principled way of quantifying multiple graded or categorical relevance levels for use as gains in established graded relevance measures, such as normalized discounted cumulative gain (nDCG), which nowadays often use heuristic and data-independent gain values. Given a set of test topics with graded relevance judgments, the PRM allows evaluating systems on different scenarios, such as their capability of retrieving top results, or how well they are able to filter out non-relevant ones. Its use in actual evaluation scenarios is illustrated on several information retrieval test collections.


title: Multiple Changepoint Detection with Partial Information on Changepoint Times

abstract: This paper proposes a new minimum description length procedure to detect multiple changepoints in time series data when some times are a priori thought more likely to be changepoints. This scenario arises with temperature time series homogenization pursuits, our focus here. Our Bayesian procedure constructs a natural prior distribution for the situation, and is shown to estimate the changepoint locations consistently, with an optimal convergence rate. Our methods substantially improve changepoint detection power when prior information is available. The methods are also tailored to bivariate data, allowing changes to occur in one or both component series.


title: Extension of structure integration to magnetic system

abstract: We extend our previously developed method, "structure integration", to evaluate free energy directly for magnetic large systems on given lattice. The present method express the density of states(DOS) as the parameters independent of system but depend on lattice, and replace the DOS unknown a priori as that known a priori. Through two-dimensional square lattice Ising model, we find that the present method can evaluate magnetic free energy efficiently and accurately above critical temperature without iterative method like Monte Carlo simulation.


title: Quasicircles of dimension 1+k^2 do not exist

abstract: A well-known theorem of S. Smirnov states that the Hausdorff dimension of a $k$-quasicircle is at most $1+k^2$. Here, we show that the precise upper bound $D(k) = 1+\Sigma^2 k^2 + \mathcal O(k^{8/3-\varepsilon})$ where $\Sigma^2$ is the maximal asymptotic variance of the Beurling transform, taken over the unit ball of $L^\infty$. The quantity $\Sigma^2$ was introduced in a joint work with K. Astala, A. Per\"al\"a and I. Prause where it was proved that $0.879 < \Sigma^2 \le 1$, while recently, H. Hedenmalm discovered that surprisingly $\Sigma^2 <1$. We deduce the asymptotic expansion of $D(k)$ from a more general statement relating the universal bounds for the integral means spectrum and the asymptotic variance of conformal maps. Our proof combines fractal approximation techniques with the classical argument of J. Becker and Ch. Pommerenke for estimating integral means.


title: Cleaning graphene : a first quantum/classical molecular dynamics approach

abstract: Graphene outstanding properties created a huge interest in the condensed matter community and unprecedented fundings at the international scale in the hope of application developments. Recently, there have been several reports of incomplete removal of the polymer resists used to transfer as-grown graphene from one substrate to another, resulting in altered graphene transport properties. Finding a large-scale solution to clean graphene from adsorbed residues is highly desirable and one promising possibility would be to use hydrogen plasmas. In this spirit, we couple here quantum and classical molecular dynamics simulations to explore the kinetic energy ranges required by atomic hydrogen to selectively etch a simple residue, a CH3 group, without irreversibly damaging the graphene. For incident energies in the 2-15 eV range, the CH3 radical can be etched by forming a volatile CH4 compound which leaves the surface, either in the CH4 form or breaking into CH3+H fragments, without further defect formation. At this energy, adsorption of H atoms on graphene is possible and further annealing will be required to recover pristine graphene.


title: Statistical Properties of Gray Color Categorization: Asymptotics of Psychometric Function

abstract: The results of our experiments on categorical perception of different shades of gray are reported. A special color generator was created for conducting the experiments on categorizing a random sequence of colors into two classes, light-gray and dark-gray. The collected data are analyzed based on constructing (i) the asymptotics of the corresponding psychometric functions and (ii) the mean decision time in categorizing a given shade of gray depending on the shade brightness (shade number). Conclusions about plausible mechanisms governing categorical perception, at least for the analyzed system, are drawn.


title: Theta characteristics of hyperelliptic graphs

abstract: We study theta characteristics of hyperelliptic metric graphs of genus $g$ with no bridge edges. These graphs have a harmonic morphism of degree two to a metric tree that can be lifted to morphism of degree two of a hyperelliptic curve $X$ over $K$ to the projective line, with $K$ an algebraically closed field of char$(K) \not =2$, complete with respect to a non-Archimedean valuation, with residue field $k$ of char$(k)\not=2$. The hyperelliptic curve has $2^{2g}$ theta characteristics. We show that for each effective theta characteristics on the graph, $2^{g-1}$ even and $2^{g-1}$ odd theta characteristics on the curve specialize to it; and $2^g$ even theta characteristics on the curve specialize to the unique not effective theta characteristics on the graph.


title: Nonlocal and multipoint boundary value problems for linear evolution equations

abstract: We derive the solution representation for a large class of nonlocal boundary value problems for linear evolution PDEs with constant coefficients in one space variable. The prototypical such PDE is the heat equation, for which problems of this form model physical phenomena in chemistry and for which we formulate and prove a full result. We also consider the third order case, which is much less studied and has been shown by the authors to have very different structural properties in general. The nonlocal conditions we consider can be reformulated as \emph{multipoint conditions}, and then an explicit representation for the solution of the problem is obtained by an application of the Fokas transform method. The analysis is carried out under the assumption that the problem being solved is well posed, i.e.\ that it admits a unique solution. For the second order case, we also give criteria that guarantee well-posedness.


title: Characterizing Model-Theoretic Dividing Lines via Collapse of Generalized Indiscernibles

abstract: We use the notion of collapse of generalized indiscernible sequences to classify various model theoretic dividing lines. In particular, we use collapse of n-multi-order indiscernibles to characterize op-dimension n; collapse of function-space indiscernibles (i.e. parameterized equivalence relations) to characterize rosy theories; and finally, convex equivalence relation indiscernibles to characterize NTP2 theories.


title: Debris Engine: A Potential Thruster for Space Debris Removal

abstract: We present a design concept for a space engine that can continuously remove the orbit debris by using the debris as a propellant. Space robotic cleaner is adopted to capture the targeting debris and to transfer them into the engine. Debris with larger size is first disintegrated into small pieces by using a mechanical method. The planetary ball mill is then adopted to grind the pieces into micrometer or smaller powder. The energy needed in this process is get from the nuclear and solar power. By the effect of gamma-ray photoelectric or the behavior of tangently rub of tungsten needles, the debris powered is charged. This behavior can be used to speed up the movement of powder in a tandem electrostatic particle accelerator. By ejecting the high-temperture and high-pressure charged powered from the nozzle of the engine,the continuously thrust is obtained. This thrust can be used to perform orbital maneuver and debris rendezvous for the spacecraft and robotic cleaner. The ejected charged particle will be blown away from the circumterrestrial orbit by the solar wind. By digesting the space debris, we obtain not only the previous thrust but also the clean space. In the near future, start trek will not just a dream, human exploration will extend to deep universe. The analysis shown, the magnitude of the specific impulse for debris engine is determined by the accelerating electrostatic potential and the charge-to-mass ratio of the powder.


title: NetVLAD: CNN architecture for weakly supervised place recognition

abstract: We tackle the problem of large scale visual place recognition, where the task is to quickly and accurately recognize the location of a given query photograph. We present the following three principal contributions. First, we develop a convolutional neural network (CNN) architecture that is trainable in an end-to-end manner directly for the place recognition task. The main component of this architecture, NetVLAD, is a new generalized VLAD layer, inspired by the "Vector of Locally Aggregated Descriptors" image representation commonly used in image retrieval. The layer is readily pluggable into any CNN architecture and amenable to training via backpropagation. Second, we develop a training procedure, based on a new weakly supervised ranking loss, to learn parameters of the architecture in an end-to-end manner from images depicting the same places over time downloaded from Google Street View Time Machine. Finally, we show that the proposed architecture significantly outperforms non-learnt image representations and off-the-shelf CNN descriptors on two challenging place recognition benchmarks, and improves over current state-of-the-art compact image representations on standard image retrieval benchmarks.


title: Multifractal analysis of weighted networks by a modified sandbox algorithm

abstract: Complex networks have attracted growing attention in many fields. As a generalization of fractal analysis, multifractal analysis (MFA) is a useful way to systematically describe the spatial heterogeneity of both theoretical and experimental fractal patterns. Some algorithms for MFA of unweighted complex networks have been proposed in the past a few years, including the sandbox (SB) algorithm recently employed by our group. In this paper, a modified SB algorithm (we call it SBw algorithm) is proposed for MFA of weighted networks.First, we use the SBw algorithm to study the multifractal property of two families of weighted fractal networks (WFNs): "Sierpinski" WFNs and "Cantor dust" WFNs. We also discuss how the fractal dimension and generalized fractal dimensions change with the edge-weights of the WFN. From the comparison between the theoretical and numerical fractal dimensions of these networks, we can find that the proposed SBw algorithm is efficient and feasible for MFA of weighted networks. Then, we apply the SBw algorithm to study multifractal properties of some real weighted networks ---collaboration networks. It is found that the multifractality exists in these weighted networks, and is affected by their edge-weights.


title: On the total $(k,r)$-domination number of random graphs

abstract: A subset $S$ of a vertex set of a graph $G$ is a total $(k,r)$-dominating set if every vertex $u \in V(G)$ is within distance $k$ of at least $r$ vertices in $S$. The minimum cardinality among all total $(k,r)$-dominating sets of $G$ is called the total $(k,r)$-domination number of $G$, denoted by $\gamma^{t}_{(k,r)}(G)$. We previously gave an upper bound on $\gamma^{t}_{(2,r)}(G(n,p))$ in random graphs with non-fixed $p \in (0,1)$. In this paper we generalize this result to give an upper bound on $\gamma^{t}_{(k,r)}(G(n,p))$ in random graphs with non-fixed $p \in (0,1)$ for $k\geq 3$ as well as present an upper bound on $\gamma^{t}_{(k,r)}(G)$ in graphs with large girth.


title: Two weight inequalities for bilinear forms

abstract: Let $1\le p_0<p,q <q_0\le \infty$. Given a pair of weights $(w,\sigma)$ and a sparse family $\mathcal S$, we study the two weight inequality for the following bi-sublinear form \[ B(f, g)= \sum_{Q\in\mathcal S}\langle |f|^{p_0}\rangle_Q^{\frac 1{p_0}} \langle|g|^{q_0'}\rangle_Q^{\frac 1{q_0'}}\lambda_Q\le \mathcal N\|f\|_{L^{p}(w)}\|g\|_{L^{q'}(\sigma)}. \] When $\lambda_Q=|Q|$ and $p=q$, Bernicot, Frey and Petermichl showed that $B(f,g)$ dominates $\langle Tf, g\rangle$ for a large class of singular non-kernel operators. We give a characterization for the above inequality and then obtain the mixed $A_p$-$A_\infty$ estimates and the corresponding entropy bounds when $\lambda_Q=|Q|$ and $p=q$. We also proposed a new conjecture which implies both the one supremum conjecture and the separated bump conjecture.


title: Full escape of mass for the diagonal group

abstract: Building on the work of Cassels we prove the existence of infinite families of compact orbits of the diagonal group in the space of lattices which accumulate only on the divergent orbit of the standard lattice. As a consequence, we prove the existence of full escape of mass for sequences of such orbits, sharpening known results about possible partial escape of mass. The topological complexity of the visits of these orbits to a given compact set in the space of lattices is analyzed. Along the way we develop new tools to compute regulators of orders in a number fields.


title: Skew-morphisms of cyclic $p$-groups

abstract: Let $G$ be a finite group having a factorisation $G=AB$ into subgroups $A$ and $B$ with $B$ cyclic and $A\cap B=1,$ and let $b$ be a generator of $B$. The associated skew-morphism is the bijective mapping $f:A \to A$ well defined by the equality $baB=f(a)B$ where $a\in A$. In this paper, we shall classify all skew-morphisms of cyclic $p$-groups where $p$ is an odd prime.


title: New results on maximal partial line spreads in PG(5,q)

abstract: In this work, we prove the existence of maximal partial line spreads in PG(5,q) of size q^3+q^2+kq+1, with 1 \leq k \leq (q^3-q^2)/(q+1), k an integer. Moreover, by a computer search, we do this for larger values of k, for q \leq 7. Again by a computer search, we find the sizes for the largest maximal partial line spreads and many new results for q \leq 5.


title: Direct numerical simulation of single mode three-dimensional Rayleigh-Taylor experiments

abstract: The single-mode Rayleigh-Taylor instability (smRTI) is well defined, poorly understood, and applicable to many fluid flows directly and through its relationship to multi-mode Rayleigh-Taylor models. This study reproduces three low-Atwood smRTI experimental runs (Wilkinson and Jacobs, 2007) in a specialized version of the Nek5000 spectral element code. The simulations use the initial amplitude, wavelength, acceleration, Atwood number, and viscosity from the three specific experiments and impose no-slip and no-flux boundaries on the velocity and scalar, respectively. The simulations are shown to reproduce the linear, saturation, stagnation, and re-acceleration phases of the smRTI seen in the experiments. Additionally, access to the full velocity and scalar fields demonstrates three different finite size effects: wall drag, wall lift, and a long wavelength mode along the diagonal. One of the simulations is extended by a factor of two in the vertical direction and the resulting late-time dynamics reach Froude numbers around 1.8, higher than previously reported. Finally, inspection of the span-wise flow reveals secondary flows of the first kind that transport the scalar from the bubble-spike interfaces into the bubble and spike centers. The agreement between simulations and experiments inspires confidence in the spectral element method for studying the Rayleigh-Taylor instability.


title: Theoretical Results on Neutrinos

abstract: In this talk, I first summarize our current knowledge about the fundamental properties of neutrinos and emphasize the remaining unsolved problems in neutrino physics. Then, recent theoretical results on neutrino mass models are introduced. Different approaches to understanding tiny neutrino masses, lepton flavor mixing and CP violation are presented. Finally, I report briefly some new progress in the studies of astrophysical neutrinos, including keV sterile neutrinos, supernova neutrinos and ultrahigh-energy cosmic neutrinos.


title: Notes on the SL(2,R) CFT

abstract: This is a set of notes which reviews and addresses issues in the SL(2,R) conformal field theory, while working primarily in a basis of vertex operators of definite weight under the affine algebra. Following a review of the H3 coset model and the spectrum of the SL(2,R) CFT, derivations are given for the selection rules and the reduction of descendant to primary correlators involving arbitrary spectral flowed modules. Together with a description of the equivalence of correlation functions of fixed total spectral flow, this leads to an enumeration of the minimal set of three-point amplitudes required to determine the fusion rules. The corresponding fusion rules for elements of the spectrum are compared to those which arise through analytic continuation of the OPE of primary fields in the H3 model, and the apparent failure of the closure of the associated SL(2,R) OPE is described in detail. A discussion of the consistency of the respective fusion rules with the conjectured structure of equivalent factorizations of SL(2,R) four-point correlators is presented.


title: Probing the Local Bubble with Diffuse Interstellar Bands (DIBs)

abstract: The Sun lies in the middle of an enormous cavity of a million degree gas, known as the Local Bubble. The Local Bubble is surrounded by a wall of denser neutral and ionized gas. The Local Bubble extends around 100 pc in the plane of Galaxy and hundreds of parsecs vertically, but absorption-line surveys of neutral sodium and singly-ionized calcium have revealed a highly irregular structure and the presence of neutral clouds within an otherwise tenuous and hot gas. We have undertaken an all-sky, European-Iranian survey of the Local Bubble in the absorption of a number of diffuse interstellar bands (DIBs) to offer a novel view of our neighbourhood. Our dedicated campaigns with ESO's New Technology Telescope and the ING's Isaac Newton Telescope comprise high signal-to-noise, medium-resolution spectra, concentrating on the 5780 and 5797 \AA\ bands which trace ionized/irradiated and neutral/shielded environments, respectively; their carriers are unknown but likely to be large carbonaceous molecules. With about 660 sightlines towards early-type stars distributed over distances up to about 200 pc, our data allow us to reconstruct the first ever 3D DIB map of the Local Bubble, which we present here. While we confirm our expectations that the 5780 \AA\ DIB is relatively strong compared to the 5797 \AA\ DIB in hot/irradiated regions such as which prevail within the Local Bubble and its walls, and the opposite is true for cooler/shielded regions beyond the confines of the Local Bubble, we unexpectedly also detect DIB cloudlets inside of the Local Bubble. These results reveal new insight into the structure of the Local Bubble, as well as helping constrain our understanding of the carriers of the DIBs.


title: Direct Detection of the L-Dwarf Donor in WZ Sagittae

abstract: Analysis of a large set of phase-resolved $K$-band spectra of the cataclysmic variable WZ Sge shows that the secondary star of this system appears to be an L-dwarf. Previous $K$-band spectra of WZ Sge found that the CO overtone bandheads were in emission. We show that absorption from the $^{\rm 12}$CO$_{\rm (2,0)}$ bandhead of the donor star creates a dip in the $^{\rm 12}$CO$_{\rm (2,0)}$ emission feature. Measuring the motion of this feature over the orbital period, we construct a radial velocity curve that gives a velocity amplitude of K$_{\rm abs}$ = 520 $\pm$ 35 km s$^{\rm -1}$, consistent with the previously published values for this parameter.


title: Application of planar Randers geodesics with river-type perturbation in search models

abstract: We consider remodeling the planar search patterns, in the presence of the river-type perturbation represented by the weak vector field, basing on the time-optimal paths as Finslerian solutions to the Zermelo navigation problem via Randers metric.


title: IPOPv2 online service for the generation of opacity tables

abstract: In the framework of the present phase -- IPOPv2 -- of the international Opacity Project (OP), a new web service has been implemented based on the latest release of the OP opacities. The user may construct online opacity tables to be conveniently included in stellar evolution codes in the format most commonly adopted by stellar physicists, namely the OPAL format. This facility encourages the use and comparison of both the OPAL and OP data sets in applications. The present service allows for the calculation of multi-element mixtures containing the 17 species (H, He, C, N, O, Ne, Na, Mg, Al, Si, S, Ar, Ca, Cr, Mn, Fe and Ni) considered by the OP, and underpins the latest release of OP opacities. This new service provides tables of Rosseland mean opacites using OP atomic data. We provide an alternative to the OPAL opacity services allowing direct comparison as well as study of the effect of uncertainties in stellar modeling due to mean opacities.


title: A Python Extension for the Massively Parallel Multiphysics Simulation Framework waLBerla

abstract: We present a Python extension to the massively parallel HPC simulation toolkit waLBerla. waLBerla is a framework for stencil based algorithms operating on block-structured grids, with the main application field being fluid simulations in complex geometries using the lattice Boltzmann method. Careful performance engineering results in excellent node performance and good scalability to over 400,000 cores. To increase the usability and flexibility of the framework, a Python interface was developed. Python extensions are used at all stages of the simulation pipeline: They simplify and automate scenario setup, evaluation, and plotting. We show how our Python interface outperforms the existing text-file-based configuration mechanism, providing features like automatic nondimensionalization of physical quantities and handling of complex parameter dependencies. Furthermore, Python is used to process and evaluate results while the simulation is running, leading to smaller output files and the possibility to adjust parameters dependent on the current simulation state. C++ data structures are exported such that a seamless interfacing to other numerical Python libraries is possible. The expressive power of Python and the performance of C++ make development of efficient code with low time effort possible.


title: R-Matrix incorporating Time (RMT) method for H2+ in short and intense laser fields

abstract: In this work we develop an approach for a molecular hydrogen ion (H2+) in the Born-Oppenheimer approximation while exposed to intense short-pulse radiation. Our starting point is the R-Matrix incorporating Time (RMT) formulation for atomic hydrogen [L. A. A. Nikolopoulos et al, Phys. Rev. A 78, 063420 (2008)] which has proven to be successful at treating multi-electron atomic systems efficiently and to high accuracy [L. R. Moore et al J. Mod. Opt. 58,1132, (2011)]. The present study on H2+ has been performed with a similar objective of developing an ab initio method for solving the Time-dependent Schr\"odinger Equation (TDSE) for multi-electron diatomic molecules exposed to an external time-dependent potential field. The theoretical formulation is developed in detail for the molecular hydrogen ion where all the multi-electron and inter-nuclei complications are absent. As in the atomic case, the configuration space of the electron's coordinates are separated artificially over two regions; the inner (I) and outer (II) regions. In the region I the time-dependent wavefunction is expanded on the eigenstate basis corresponding to the molecule's Hamiltonian augmented by Bloch operators, while in region II a grid representation is used. We demonstrate the independence of our results on the introduced artificial boundary-surface by calculating observables that are directly accessed experimentally and also by showing gauge-dependent quantities are also invariant with the region I box size. We also compare with other theoretical works and emphasize cases where basis-set approaches are currently very computationally expensive or intractable in terms of computational resources.


title: Input Sparsity Time Low-Rank Approximation via Ridge Leverage Score Sampling

abstract: We present a new algorithm for finding a near optimal low-rank approximation of a matrix $A$ in $O(nnz(A))$ time. Our method is based on a recursive sampling scheme for computing a representative subset of $A$'s columns, which is then used to find a low-rank approximation. This approach differs substantially from prior $O(nnz(A))$ time algorithms, which are all based on fast Johnson-Lindenstrauss random projections. It matches the guarantees of these methods while offering a number of advantages. Not only are sampling algorithms faster for sparse and structured data, but they can also be applied in settings where random projections cannot. For example, we give new single-pass streaming algorithms for the column subset selection and projection-cost preserving sample problems. Our method has also been used to give the fastest algorithms for provably approximating kernel matrices [MM16].


title: Nonlocal spin torques in Rashba quantum wires with steep magnetic textures

abstract: We provide a general procedure to calculate the current-induced spin-transfer torque which acts on a general steep magnetic texture due to the exchange interaction with an applied spin-polarized current. As an example, we consider a one-dimensional ferromagnetic quantum wire and also include a Rashba spin-orbit interaction. The spin-transfer torque becomes generally spatially non-local. Likewise, the Rashba spin-orbit interaction induces a spatially nonlocal field-like nonequilibrium spin-transfer torque. We also find a spatially varying nonadiabaticity parameter and markedly different domain wall dynamics for very steep textures as compared to wide domain walls.


title: Extending the Li&Ma method to include PSF information

abstract: The so called Li&Ma formula is still the most frequently used method for estimating the significance of observations carried out by Imaging Atmospheric Cherenkov Telescopes. In this work a straightforward extension of the method for point sources that profits from the good imaging capabilities of current instruments is proposed. It is based on a likelihood ratio under the assumption of a well-known PSF and a smooth background. Its performance is tested with Monte Carlo simulations based on real observations and its sensitivity is compared to standard methods which do not incorporate PSF information. The gain of significance that can be attributed to the inclusion of the PSF is around of 10% and can be boosted if a background model is assumed or a finer binning is used.


title: Noise Induced Phenomena in the Dynamics of Two Competing Species

abstract: Noise through its interaction with the nonlinearity of the living systems can give rise to counter-intuitive phenomena. In this paper we shortly review noise induced effects in different ecosystems, in which two populations compete for the same resources. We also present new results on spatial patterns of two populations, while modeling real distributions of anchovies and sardines. The transient dynamics of these ecosystems are analyzed through generalized Lotka-Volterra equations in the presence of multiplicative noise, which models the interaction between the species and the environment. We find noise induced phenomena such as quasi-deterministic oscillations, stochastic resonance, noise delayed extinction, and noise induced pattern formation. In addition, our theoretical results are validated with experimental findings. Specifically the results, obtained by a coupled map lattice model, well reproduce the spatial distributions of anchovies and sardines, observed in a marine ecosystem. Moreover, the experimental dynamical behavior of two competing bacterial populations in a meat product and the probability distribution at long times of one of them are well reproduced by a stochastic microbial predictive model.


title: Automatic Reasoning on Recursive Data-Structures with Sharing

abstract: We consider the problem of automatically verifying programs which manipulate arbitrary data structures. Our specification language is expressive, contains a notion of \emph{separation}, and thus enables a precise specification of \emph{frames}. The main contribution then is a program verification method which combines strongest postcondition reasoning in the form symbolic execution, unfolding recursive definitions of the data structure in question, and a new frame rule to achieve \emph{local reasoning} so that proofs can be compositional. Finally, we present an implementation of our verifier, and demonstrate automation on a number of representative programs. In particular, we present the first automatic proof of a classic graph marking algorithm, paving the way for dealing with a class of programs which traverse a complex data structure.


title: Automorphism groups of Cayley graphs generated by block transpositions and regular Cayley maps

abstract: This paper deals with the Cayley graph $\mathrm{Cay}(\mathrm{Sym}_n,T_n),$ where the generating set consists of all block transpositions. A motivation for the study of these particular Cayley graphs comes from current research in Bioinformatics. As the main result, we prove that Aut$(\mathrm{Cay}(\mathrm{Sym}_n,T_n))$ is the product of the left translation group by a dihedral group $\mathsf{D}_{n+1}$ of order $2(n+1)$. The proof uses several properties of the subgraph $\Gamma$ of $\mathrm{Cay}(\mathrm{Sym}_n,T_n)$ induced by the set $T_n$. In particular, $\Gamma$ is a $2(n-2)$-regular graph whose automorphism group is $\mathsf{D}_{n+1},$ $\Gamma$ has as many as $n+1$ maximal cliques of size $2,$ and its subgraph $\Gamma(V)$ whose vertices are those in these cliques is a $3$-regular, Hamiltonian, and vertex-transitive graph. A relation of the unique cyclic subgroup of $\mathsf{D}_{n+1}$ of order $n+1$ with regular Cayley maps on $\mathrm{Sym}_n$ is also discussed. It is shown that the product of the left translation group by the latter group can be obtained as the automorphism group of a non-$t$-balanced regular Cayley map on $\mathrm{Sym}_n$.


title: Degree of commutativity of infinite groups

abstract: We prove that, in a finitely generated residually finite group of subexponential growth, the proportion of commuting pairs is positive if and only if the group is virtually abelian. In particular, this covers the case where the group has polynomial growth (i.e., virtually nilpotent groups, where the hypothesis of residual finiteness is always satisfied). We also show that, for non-elementary hyperbolic groups, the proportion of commuting pairs is always zero.


title: Spectroscopic survey of emission-line stars. I. B[e] stars

abstract: Emission-line stars are typically surrounded by dense circumstellar material, often in form of rings or disc-like structures. Line emission from forbidden transitions trace a diversity of density and temperature regimes. Of particular interest are the forbidden lines of [O I] {\lambda}{\lambda}6300, 6364 and [Ca II] {\lambda}{\lambda}7291, 7324. They arise in complementary, high-density environments, such as the inner-disc regions around B[e] supergiants. To study physical conditions traced by these lines and to investigate how common they are, we initiated a survey of emission-line stars. Here, we focus on a sample of nine B[e] stars in different evolutionary phases. Emission of the [O I] lines is one of the characteristics of B[e] stars. We find that four of the objects display [Ca II] line emission: for the B[e] supergiants V1478 Cyg and 3 Pup the kinematics obtained from the [O I] and [Ca II] line profiles agrees with a Keplerian rotating disc scenario; the forbidden lines of the compact planetary nebula OY Gem display no kinematical broadening beyond spectral resolution; the LBV candidate V1429 Aql shows no [O I] lines, but the profile of its [Ca II] lines suggests that the emission originates in its hot, ionized circumbinary disc. As none of the B[e] stars of lower mass displays [Ca II] line emission, we conclude that these lines are more likely observable in massive stars with dense discs, supporting and strengthening the suggestion that their appearance requires high-density environments.


title: Synthesizing Omnidirectional Antenna Patterns, Received Power and Path Loss from Directional Antennas for 5G Millimeter-Wave Communications

abstract: Omnidirectional path loss models are vital for radiosystem design in wireless communication systems, as they allow engineers to perform network simulations for systems with arbitrary antenna patterns. At millimeter-wave frequencies, channel measurements are frequently conducted using steerable highgain directional antennas at both the transmitter and receiver to make up for the significant increase in free space path loss at these frequencies compared to traditional cellular systems that operate at lower frequencies. The omnidirectional antenna pattern, and resulting omnidirectional received power must therefore be synthesized from many unique pointing angles, where the transmit and receive antennas are rotated over many different azimuth and elevation planes. In this paper, the equivalent omnidirectional antenna pattern and omnidirectional received power are synthesized by summing the received powers from all measured unique pointing angles obtained at antenna halfpower beamwidth step increments in the azimuth and elevation planes, and this method is validated by demonstrating that the synthesized omnidirectional received power and path loss are independent of antenna beamwidth, through theoretical analyses and millimeter-wave propagation measurements using antennas with different beamwidths. The method in this paper is shown to provide accurate results while enhancing the measurement range substantially through the use of directional antennas.


title: Estimating long-term behavior of periodically driven flows without trajectory integration

abstract: Periodically driven flows are fundamental models of chaotic behavior and the study of their transport properties is an active area of research. A well-known analytic construction is the augmentation of phase space with an additional time dimension; in this augmented space, the flow becomes autonomous or time-independent. We prove several results concerning the connections between the original time-periodic representation and the time-extended representation, focusing on transport properties. In the deterministic setting, these include single-period outflows and time-asymptotic escape rates from time-parameterized families of sets. We also consider stochastic differential equations with time-periodic advection term. In this stochastic setting one has a time-periodic generator (the differential operator given by the right-hand-side of the corresponding time-periodic Fokker-Planck equation). We define in a natural way an autonomous generator corresponding to the flow on time-extended phase space. We prove relationships between these two generator representations and use these to quantify decay rates of observables and to determine time-periodic families of sets with slow escape rate. Finally, we use the generator on the time-extended phase space to create efficient numerical schemes to implement the various theoretical constructions. These ideas build on the work of Froyland et al (SINUM, 2013), and no expensive time integration is required. We introduce an efficient new hybrid approach, which treats the space and time dimensions separately.


title: Modal "self-coupling" as a sensitive probe for nanomechanical detection

abstract: We present a high-sensitivity measurement technique for mechanical nanoresonators. Due to intrinsic nonlinear effects, different flexural modes of a nanobeam can be coupled while driving each of them on resonance. This mode-coupling scheme is dispersive and one mode resonance shifts with respect to the motional amplitude of the other. The same idea can be implemented on a {\it single} mode, exciting it with two slightly detuned signals. This two-tone scheme is used here to measure the resonance lineshape of one mode through a frequency shift in the response of the device. The method acts as an amplitude-to-frequency transduction which ultimately suffers only from phase noise of the local oscillator used and of the nanomechanical device itself. We also present a theory which reproduces the data without free parameters.


title: The number of trees in a graph

abstract: Let $T$ be a tree with $t$ edges. We show that the number of isomorphic (labeled) copies of $T$ in a graph $G = (V,E)$ of minimum degree at least $t$ is at least \[2|E| \prod_{v \in V} (d(v) - t + 1)^{\frac{(t-1)d(v)}{2|E|}}.\] Consequently, any $n$-vertex graph of average degree $d$ and minimum degree at least $t$ contains at least $$nd(d-t+1)^{t-1}$$ isomorphic (labeled) copies of $T$. This answers a question of Dellamonica et. al. (where the above statement was proved when $T$ is the path with three edges) while extending an old result of Erd\H os and Simonovits.


title: Learning Simple Algorithms from Examples

abstract: We present an approach for learning simple algorithms such as copying, multi-digit addition and single digit multiplication directly from examples. Our framework consists of a set of interfaces, accessed by a controller. Typical interfaces are 1-D tapes or 2-D grids that hold the input and output data. For the controller, we explore a range of neural network-based models which vary in their ability to abstract the underlying algorithm from training instances and generalize to test examples with many thousands of digits. The controller is trained using $Q$-learning with several enhancements and we show that the bottleneck is in the capabilities of the controller rather than in the search incurred by $Q$-learning.


title: A phase-field model of relaxor ferroelectrics based on random field theory

abstract: A mechanically coupled phase-field model is proposed for the first time to simulate the peculiar behavior of relaxor ferroelectrics. Based on the random field theory for relaxors, local random fields are introduced to characterize the effect of chemical disorder. This generic model is developed from a thermodynamic framework and the microforce theory and is implemented by a nonlinear finite element method. Simulation results show that the model can reproduce relaxor features, such as domain miniaturization, small remnant polarization and large piezoelectric response. In particular, the influence of random field strength on these features are revealed. Simulation results on domain structure and hysteresis behavior are discussed and compared with related experimental results.


title: Atomic Quadrupole Moment Measurement Using Dynamic Decoupling

abstract: We present a method that uses dynamic decoupling of a multi-level quantum probe to distinguish small frequency shifts that depend on $m^2_{j}$, where $m^2_{j}$ is the angular momentum of level $\left|j\right\rangle$ along the quantization axis, from large noisy shifts that are linear in $m_{j}$, such as those due to magnetic field noise. Using this method we measured the electric quadrupole moment of the $4D_{\frac{5}{2}}$ level in $^{88}Sr^{+}$ to be $2.973^{+0.026}_{-0.033}\, ea_{0}^{2}$. Our measurement improves the uncertainty of this value by an order of magnitude and thus helps mitigate an important systematic uncertainty in $^{88}Sr^{+}$ based optical atomic clocks and verifies complicated many-body quantum calculations.


title: The difference between two random mixed quantum states: exact and asymptotic spectral analysis

abstract: We investigate the spectral statistics of the difference of two density matrices, each of which is independently obtained by partially tracing a random bipartite pure quantum state. We first show how a closed-form expression for the exact joint eigenvalue probability density function for arbitrary dimensions can be obtained from the joint probability density function of the diagonal elements of the difference matrix, which is straightforward to compute. Subsequently, we use standard results from free probability theory to derive a relatively simple analytic expression for the asymptotic eigenvalue density (AED) of the difference matrix ensemble, and using Carlson's theorem, we obtain an expression for its absolute moments. These results allow us to quantify the typical asymptotic distance between the two random mixed states using various distance measures; in particular, we obtain the almost sure asymptotic behavior of the operator norm distance and the trace distance.


title: Parabolic models for chemotaxis on weighted networks

abstract: In this work we consider the Keller-Segel model for chemotaxis on networks, both in the doubly parabolic case and in the parabolic-elliptic one. Introducing appropriate transition conditions at vertices, we prove the existence of a time global and spatially continuous solution for each of the two systems. The main tool is the use of the explicit formula for the fundamental solution of the heat equation on a weighted graph and of the corresponding sharp estimates.


title: Nonlinear parametric amplification in a tri-port nanoelectromechanical device

abstract: We report on measurements performed at low temperatures on a nanoelectromechanical system (NEMS) under (capacitive) parametric pumping. The excitations and detection schemes are purely electrical, and enable in the present experiment the straightforward measurement of forces down to about a femtonewton, for displacements of an Angstr\"om, using standard room temperature electronics. We demonstrate that a small (linear) force applied on the device can be amplified up to more than a 100 times, while the system is {\it truly moving}. We explore the dynamics up to about 50$~$nm deflections for cantilevers about 200$~$nm thick by 3$~$$\mu$m long oscillating at a frequency of 7$~$MHz. We present a generic modeling of nonlinear parametric amplification, and give analytic theoretical solutions enabling the fit of experimental results. We finally discuss the practical limits of the technique, with a particular application: the measurement of {\it anelastic damping} in the metallic coating of the device with an exceptional resolution of about 0.5$~$\%.


title: Sparse Linear Models applied to Power Quality Disturbance Classification

abstract: Power quality (PQ) analysis describes the non-pure electric signals that are usually present in electric power systems. The automatic recognition of PQ disturbances can be seen as a pattern recognition problem, in which different types of waveform distortion are differentiated based on their features. Similar to other quasi-stationary signals, PQ disturbances can be decomposed into time-frequency dependent components by using time-frequency or time-scale transforms, also known as dictionaries. These dictionaries are used in the feature extraction step in pattern recognition systems. Short-time Fourier, Wavelets and Stockwell transforms are some of the most common dictionaries used in the PQ community, aiming to achieve a better signal representation. To the best of our knowledge, previous works about PQ disturbance classification have been restricted to the use of one among several available dictionaries. Taking advantage of the theory behind sparse linear models (SLM), we introduce a sparse method for PQ representation, starting from overcomplete dictionaries. In particular, we apply Group Lasso. We employ different types of time-frequency (or time-scale) dictionaries to characterize the PQ disturbances, and evaluate their performance under different pattern recognition algorithms. We show that the SLM reduce the PQ classification complexity promoting sparse basis selection, and improving the classification accuracy.


title: Magnetic neutron scattering by magnetic vortices in thin submicron-sized soft ferromagnetic cylinders

abstract: Using analytical expressions for the magnetization textures of thin submicron-sized magnetic cylinders in vortex state, we derive closed-form algebraic expressions for the ensuing small-angle neutron scattering (SANS) cross sections. Specifically, for the perpendicular and parallel scattering geometries, we have computed the magnetic SANS cross sections for the case of small vortex-center displacements without formation of magnetic charges on the side faces of the cylinder. The results represent a significant qualitative and quantitative step forward in SANS-data analysis on isolated magnetic nanoparticle systems, which are commonly assumed to be homogeneously or stepwise-homogeneously magnetized. We suggest a way to extract the fine details of the magnetic vortex structure during the magnetization process from the SANS measurements in order to help resolving the long standing question of the magnetic vortex displacement mode.


title: Exploring the Inner Edge of the Habitable Zone with Fully Coupled Oceans

abstract: Rotation in planetary atmospheres plays an important role in regulating atmospheric and oceanic heat flow, cloud formation and precipitation. Using the Goddard Institute for Space Studies (GISS) three dimension General Circulation Model (3D-GCM) we investigate how the effects of varying rotation rate and increasing the incident stellar flux on a planet set bounds on a planet's habitable zone with its parent star. From ensemble climate simulations we identify which factors are the primary controllers of uncertainty in setting these bounds. This is shown in particular for fully coupled ocean (FCO) runs -- some of the first that have been utilized in this context. Results with a Slab Ocean (SO) of 100m mixed layer depth are compared with a similar study by Yang et al. 2014, which demonstrates consistency across models. However, there are clear differences for rotations rates of 1-16x present Earth sidereal day lengths between the 100m SO and FCO models, which points to the necessity of using FCOs whenever possible. The latter was recently demonstrated quite clearly by Hu & Yang 2014 in their aquaworld study with a FCO when compared with similar mixed layer ocean studies and by Cullum et al. 2014. We also show how these results have implications for Venus in the early history of our Solar System since even at this time Venus received more solar flux than Earth does today while it may still have had a slow retrograde rotation. The Venus runs utilize a 2.9Gya solar spectrum generated with the code of Claire et al. 2012, a modern Venus topography with an ocean filling the lowlands (giving an equivalent depth of 310 meters if spread across the entire surface), atmosphere of 1 bar N2, CO2=0.4mb, CH4=0.001mb and present day orbital parameters, radius, & gravity. We demonstrate that ancient Venus could have had quite moderate surface temperatures given these assumptions.


title: Towards spin turbulence of light: Spontaneous disorder and chaos in cavity-polariton systems

abstract: Recent advances in nanophotonics have brought about coherent light sources with chaotic circular polarization; a low-dimensional chaotic evolution of optical spin was evidenced in laser diodes. Here we propose a mechanism that gives rise to light with a spatiotemporal spin chaos resembling turbulent states in hydrodynamics. The spin-chaotic radiation is emitted by exciton polaritons under resonant optical pumping in arbitrarily sized planar microcavities, including, as a limiting case, pointlike systems with only three degrees of freedom. The underlying mechanism originates in the interplay between spin symmetry breakdown and scattering into pairs of Bogolyubov excitations. As a practical matter, it opens up the way for spin modulation of light on the scale of picoseconds and micrometers.


title: Cyclic groups are CI-groups for balanced configurations

abstract: In this paper it is shown that every finite cyclic group satisfies the CI-property for the class of balanced configurations.


title: Polarized neutron channeling as a tool for the investigations of weakly magnetic thin films

abstract: We present and apply a new method to measure directly weak magnetization in thin films. The polarization of a neutron beam channeling through a thin film structure is measured after exiting the structure edge as a microbeam. We have applied the method to a tri-layer thin film structure acting as a planar waveguide for polarized neutrons. The middle guiding layer is a rare earth based ferrimagnetic material TbCo5 with a low magnetization of about 20 mT. We demonstrate that the channeling method is more sensitive than the specular neutron reflection method.