title: Near UV Imager with an MCP Based Photon Counting Detector

abstract: We are developing a compact UV Imager using light weight components, that can be flown on a small CubeSat or a balloon platform. The system has a lens-based optics that can provide an aberration-free image over a wide field of view. The backend instrument is a photon counting detector with off-the-shelf MCP, CMOS sensor and electronics. We are using a Z-stack MCP with a compact high voltage power supply and a phosphor screen anode, which is read out by a CMOS sensor and the associated electronics. The instrument can be used to observe solar system objects and detect bright transients from the upper atmosphere with the help of CubeSats or high altitude balloons. We have designed the imager to be capable of working in direct frame transfer mode as well in the photon-counting mode for single photon event detection. The identification and centroiding of each photon event are done using an FPGA-based data acquisition and real-time processing system.


title: On the intersection of tame subgroups in groups acting on trees

abstract: Let $G$ be a group acting on a tree $T$ with finite edge stabilizers of bounded order. We provide, in some very interesting cases, upper bounds for the complexity of the intersection $H\cap K$ of two tame subgroups $H$ and $K$ of $G$ in terms of the complexities of $H$ and $K$. In particular, we obtain bounds for the Kurosh rank $Kr(H\cap K)$ of the intersection in terms of Kurosh ranks $Kr(H)$ and $Kr(K)$, in the case where $H$ and $K$ act freely on the edges of $T$.


title: Flat polarization-controlled cylindrical lens based on the Pancharatnam-Berry geometric phase

abstract: The working principle of ordinary refractive lenses can be explained in terms of the space-variant optical phase retardations they introduce, which reshape the optical wavefront curvature and hence affect the subsequent light propagation. These phases, in turn, are due to the varying optical path length seen by light at different transverse positions relative to the lens centre. A similar lensing behavior can however be obtained when the optical phases are introduced by an entirely different mechanism. Here, we consider the "geometric phases" that arise from the polarization transformations occurring in anisotropic optical media, named after Pancharatnam and Berry. The medium anisotropy axis is taken to be space-variant in the transverse plane and the resulting varying geometric phases give rise to the wavefront reshaping and lensing effect, which however depends also on the input polarization. We describe the realization and characterization of a cylindrical geometric-phase lens that is converging for a given input circular polarization state and diverging for the orthogonal one, which provides one of the simplest possible examples of optical element based on geometric phases. The demonstrated lens is flat and only few microns thick (not including the supporting substrates); moreover, its working wavelength can be tuned and the lensing can be switched on and off by the action of an external control electric field. Other kinds of lenses or more general phase elements inducing different wavefront distortions can be obtained by a similar approach. Besides their potential for optoelectronic technology, these devices offer good opportunities for introducing college-level students to an advanced topic of modern physics, such as the Berry phase, with the help of interesting optical demonstrations.


title: Residue-type indices and holomorphic foliations

abstract: We investigate residue-type indices for germs of holomorphic foliations in the plane and characterize second type foliations - those not containing tangent saddle-nodes in the reduction of singularities - by an expression involving the Baum-Bott, variation and polar excess indices. These local results are applied in the study of logarithmic foliations on compact complex surfaces.


title: Deep artifact learning for compressed sensing and parallel MRI

abstract: Purpose: Compressed sensing MRI (CS-MRI) from single and parallel coils is one of the powerful ways to reduce the scan time of MR imaging with performance guarantee. However, the computational costs are usually expensive. This paper aims to propose a computationally fast and accurate deep learning algorithm for the reconstruction of MR images from highly down-sampled k-space data. Theory: Based on the topological analysis, we show that the data manifold of the aliasing artifact is easier to learn from a uniform subsampling pattern with additional low-frequency k-space data. Thus, we develop deep aliasing artifact learning networks for the magnitude and phase images to estimate and remove the aliasing artifacts from highly accelerated MR acquisition. Methods: The aliasing artifacts are directly estimated from the distorted magnitude and phase images reconstructed from subsampled k-space data so that we can get an aliasing-free images by subtracting the estimated aliasing artifact from corrupted inputs. Moreover, to deal with the globally distributed aliasing artifact, we develop a multi-scale deep neural network with a large receptive field. Results: The experimental results confirm that the proposed deep artifact learning network effectively estimates and removes the aliasing artifacts. Compared to existing CS methods from single and multi-coli data, the proposed network shows minimal errors by removing the coherent aliasing artifacts. Furthermore, the computational time is by order of magnitude faster. Conclusion: As the proposed deep artifact learning network immediately generates accurate reconstruction, it has great potential for clinical applications.


title: On Parameterized Complexity of Group Activity Selection Problems on Social Networks

abstract: In Group Activity Selection Problem (GASP), players form coalitions to participate in activities and have preferences over pairs of the form (activity, group size). Recently, Igarashi et al. have initiated the study of group activity selection problems on social networks (gGASP): a group of players can engage in the same activity if the members of the group form a connected subset of the underlying communication structure. Igarashi et al. have primarily focused on Nash stable outcomes, and showed that many associated algorithmic questions are computationally hard even for very simple networks. In this paper we study the parameterized complexity of gGASP with respect to the number of activities as well as with respect to the number of players, for several solution concepts such as Nash stability, individual stability and core stability. The first parameter we consider in the number of activities. For this parameter, we propose an FPT algorithm for Nash stability for the case where the social network is acyclic and obtain a W[1]-hardness result for cliques (i.e., for classic GASP); similar results hold for individual stability. In contrast, finding a core stable outcome is hard even if the number of activities is bounded by a small constant, both for classic GASP and when the social network is a star. Another parameter we study is the number of players. While all solution concepts we consider become polynomial-time computable when this parameter is bounded by a constant, we prove W[1]-hardness results for cliques (i.e., for classic GASP).


title: First-Order Logic with Counting: At Least, Weak Hanf Normal Forms Always Exist and Can Be Computed!

abstract: We introduce the logic FOCN(P) which extends first-order logic by counting and by numerical predicates from a set P, and which can be viewed as a natural generalisation of various counting logics that have been studied in the literature. We obtain a locality result showing that every FOCN(P)-formula can be transformed into a formula in Hanf normal form that is equivalent on all finite structures of degree at most d. A formula is in Hanf normal form if it is a Boolean combination of formulas describing the neighbourhood around its tuple of free variables and arithmetic sentences with predicates from P over atomic statements describing the number of realisations of a type with a single centre. The transformation into Hanf normal form can be achieved in time elementary in $d$ and the size of the input formula. From this locality result, we infer the following applications: (*) The Hanf-locality rank of first-order formulas of bounded quantifier alternation depth only grows polynomially with the formula size. (*) The model checking problem for the fragment FOC(P) of FOCN(P) on structures of bounded degree is fixed-parameter tractable (with elementary parameter dependence). (*) The query evaluation problem for fixed queries from FOC(P) over fully dynamic databases of degree at most d can be solved efficiently: there is a dynamic algorithm that can enumerate the tuples in the query result with constant delay, and that allows to compute the size of the query result and to test if a given tuple belongs to the query result within constant time after every database update.


title: Semiclassical expansion of the Bethe scalar products in the XXZ spin chain

abstract: In this note we adress the problem of performing the semiclassical expansion of the scalar product of Bethe states in the case of the XXZ spin chain. Our approach closely follows the one developped in [1]: after expressing the scalar products in terms of a central quantity - the B functional, we reinterpret it as a grand partition function of a gas of particles with pairwise interaction potential. This work has been done as part of the author's master thesis under the supervision of Ivan Kostov and Didina Serban at Institut de Physique Th\'eorique (CEA Saclay) in September 2015 - February 2016, and present a consice summary of the results obtained at that time.


title: Converging expansions for Lipschitz self-similar perforations of a plane sector

abstract: In contrast with the well-known methods of matching asymptotics and multiscale (or compound) asymptotics, the " functional analytic approach " of Lanza de Cristoforis (Analysis 28, 2008) allows to prove convergence of expansions around interior small holes of size $\epsilon$ for solutions of elliptic boundary value problems. Using the method of layer potentials, the asymptotic behavior of the solution as $\epsilon$ tends to zero is described not only by asymptotic series in powers of $\epsilon$, but by convergent power series. Here we use this method to investigate the Dirichlet problem for the Laplace operator where holes are collapsing at a polygonal corner of opening $\omega$. Then in addition to the scale $\epsilon$ there appears the scale $\eta = \epsilon^{\pi/\omega}$. We prove that when $\pi/\omega$ is irrational, the solution of the Dirichlet problem is given by convergent series in powers of these two small parameters. Due to interference of the two scales, this convergence is obtained, in full generality, by grouping together integer powers of the two scales that are very close to each other. Nevertheless, there exists a dense subset of openings $\omega$ (characterized by Diophantine approximation properties), for which real analyticity in the two variables $\epsilon$ and $\eta$ holds and the power series converge unconditionally. When $\pi/\omega$ is rational, the series are unconditionally convergent, but contain terms in log $\epsilon$.


title: Enhanced steady-state dissolution flux in reactive convective dissolution

abstract: Chemical reactions can accelerate, slow down or even be at the very origin of the development of dissolution-driven convection in partially miscible stratifications, when they impact the density profile in the host fluid phase. We numerically analyze the dynamics of this reactive convective dissolution in the fully developed non-linear regime for a phase A dissolving into a host layer containing a dissolved reactant B. We show that for a general A+B$\rightarrow$C reaction in solution, the dynamics vary with the Rayleigh numbers of the chemical species, i.e. with the nature of the chemicals in the host phase. Depending on whether the reaction slows down, accelerates or is at the origin of the development of convection, the spatial distributions of species A, B or C, the dissolution flux and the reaction rate are different. We show that chemical reactions enhance the steady-state flux as they consume A and can induce more intense convection than in the absence of reactions. This result is important in the context of CO$_2$ geological sequestration where quantifying the storage rate of CO$_2$ dissolving into the host oil or aqueous phase is crucial to assess the efficiency and the safety of the project.


title: Finite Blaschke products with prescribed critical points, Stieltjes polynomials, and moment problems

abstract: The determination of a finite Blaschke product from its critical points is a well-known problem with interrelations to other topics. Though existence and uniqueness of solutions are established for long, we present several new aspects which have not yet been explored to their full extent. In particular, we show that the following three problems are equivalent: (i) determining a finite Blaschke product from its critical points, (ii) finding the equilibrium position of moveable point charges interacting with a special configuration of fixed charges, (iii) solving a moment problem for the canonical representation of power moments on the real axis. These equivalences are not only of theoretical interest, but also open up new perspectives for the design of algorithms. For instance, the second problem is closely linked to the determination of certain Stieltjes and Van Vleck polynomials for a second order ODE and allows the description of solutions as global minimizers of an energy functional.


title: On the Behavior of Convolutional Nets for Feature Extraction

abstract: Deep neural networks are representation learning techniques. During training, a deep net is capable of generating a descriptive language of unprecedented size and detail in machine learning. Extracting the descriptive language coded within a trained CNN model (in the case of image data), and reusing it for other purposes is a field of interest, as it provides access to the visual descriptors previously learnt by the CNN after processing millions of images, without requiring an expensive training phase. Contributions to this field (commonly known as feature representation transfer or transfer learning) have been purely empirical so far, extracting all CNN features from a single layer close to the output and testing their performance by feeding them to a classifier. This approach has provided consistent results, although its relevance is limited to classification tasks. In a completely different approach, in this paper we statistically measure the discriminative power of every single feature found within a deep CNN, when used for characterizing every class of 11 datasets. We seek to provide new insights into the behavior of CNN features, particularly the ones from convolutional layers, as this can be relevant for their application to knowledge representation and reasoning. Our results confirm that low and middle level features may behave differently to high level features, but only under certain conditions. We find that all CNN features can be used for knowledge representation purposes both by their presence or by their absence, doubling the information a single CNN feature may provide. We also study how much noise these features may include, and propose a thresholding approach to discard most of it. All these insights have a direct application to the generation of CNN embedding spaces.


title: Statistical biases in measurements with multiple candidates

abstract: Many measurements at collider experiments study physics candidates that are a subset of a collision event. The presence of multiple such candidates in a given event can cause raw biases which are large compared to typical statistical uncertainties. Selecting a single candidate is common practice but only helps if the likelihood of selecting the true candidate is very high. Otherwise the precision of the measurement can be affected, and biases can be generated, even if none are present in the data sample prior to this operation. This paper aims at describing the problem in a systematic way. It sets definitions, provides examples of potential biases using pseudoexperiments and gives recommendations.


title: Magnetic field induced polarization enhancement in monolayers of tungsten dichalcogenides: Effects of temperature

abstract: Optical orientation of localized/bound excitons is shown to be effectively enhanced by the application of magnetic fields as low as 20 mT in monolayer WS$_2$. At low temperatures, the evolution of the polarization degree of different emission lines of monolayer WS$_2$ with increasing magnetic fields is analyzed and compared to similar results obtained on a WSe$_2$ monolayer. We study the temperature dependence of this effect up to $T=60$ K for both materials, focusing on the dynamics of the valley pseudospin relaxation. A rate equation model is used to analyze our data and from the analysis of the width of the polarization deep in magnetic field we conclude that the competition between the dark exciton pseudospin relaxation and the decay of the dark exciton population into the localized states are rather different in these two materials which are representative of the two extreme cases for the ratio of relaxation rate and depolarization rate.


title: Extensions of Operators, Liftings of Monads and Distributive Laws

abstract: In a previous study, the algebraic formulation of the First Fundamental Theorem of Calculus (FFTC) is shown to allow extensions of differential and Rota-Baxter operators on the one hand, and to give rise to categorical explanations using the ideas of liftings of monads and comonads, and mixed distributive laws on the other. Generalizing the FFTC, we consider in this paper a class of constraints between a differential operator and a Rota-Baxter operator. For a given constraint, we show that the existences of extensions of differential and Rota-Baxter operators, of liftings of monads and comonads, and of mixed distributive laws are equivalent.


title: Gaia Data Release 1. Open cluster astrometry: performance, limitations, and future prospects

abstract: Context. The first Gaia Data Release contains the Tycho-Gaia Astrometric Solution (TGAS). This is a subset of about 2 million stars for which, besides the position and photometry, the proper motion and parallax are calculated using Hipparcos and Tycho-2 positions in 1991.25 as prior information. Aims. We investigate the scientific potential and limitations of the TGAS component by means of the astrometric data for open clusters. Methods. Mean cluster parallax and proper motion values are derived taking into account the error correlations within the astrometric solutions for individual stars, an estimate of the internal velocity dispersion in the cluster, and, where relevant, the effects of the depth of the cluster along the line of sight. Internal consistency of the TGAS data is assessed. Results. Values given for standard uncertainties are still inaccurate and may lead to unrealistic unit-weight standard deviations of least squares solutions for cluster parameters. Reconstructed mean cluster parallax and proper motion values are generally in very good agreement with earlier Hipparcos-based determination, although the Gaia mean parallax for the Pleiades is a significant exception. We have no current explanation for that discrepancy. Most clusters are observed to extend to nearly 15 pc from the cluster centre, and it will be up to future Gaia releases to establish whether those potential cluster-member stars are still dynamically bound to the clusters. Conclusions. The Gaia DR1 provides the means to examine open clusters far beyond their more easily visible cores, and can provide membership assessments based on proper motions and parallaxes. A combined HR diagram shows the same features as observed before using the Hipparcos data, with clearly increased luminosities for older A and F dwarfs.


title: Analysis of a fractional-step scheme for the P1 radiative diffusion model

abstract: We address in this paper a nonlinear parabolic system, which is built to retain the main mathematical difficulties of the P1 radiative diffusion physical model. We propose a finite volume fractional-step scheme for this problem enjoying the following properties. First, we show that each discrete solution satisfies a priori L -estimates, through a discrete maxi- mum principle; by a topological degree argument, this yields the existence of a solution, which is proven to be unique. Second, we establish uniform (with respect to the size of the meshes and the time step) L2 -bounds for the space and time translates; this proves, by the Kolmogorov theorem, the relative compactness of any sequence of solutions obtained through a sequence of discretizations the time and space steps of which tend to zero; the limits of converging subsequences are then shown to be a solution to the continuous problem. Estimates of time translates of the discrete solutions are obtained through the formalization of a generic argument, interesting for its own sake.


title: CM points on Shimura curves and $p$-adic binary quadratic forms

abstract: We prove that the set of CM points on the Shimura curve associated to an Eichler order inside an indefinite quaternion $\mathbb{Q}$-algebra, is in bijection with the set of certain classes of $p$-adic binary quadratic forms, where $p$ is a prime dividing the discriminant of the quaternion algebra. The classes of $p$-adic binary quadratic forms are obtain by the action of a discrete and cocompact subgroup of $\mathrm{PGL}_{2}(\mathbb{Q}_{p})$ arising from the $p$-adic uniformization of the Shimura curve. We finally compute families of $p$-adic binary quadratic forms associated to an infinite family of Shimura curves studied in a previous paper of Amor\'os-Milione. This extends results of Alsina-Bayer to the $p$-adic context.


title: Fibre-wise linear Poisson structures related to W*-algebras

abstract: In this paper we investigate fiber-wise linear complex Banach sub-Poisson structures defined canonically by the structure of a W*-algebra M. In particular we show that these structures are arranged in the short exact sequence of complex Banach sub-Poisson VB-groupoids with the groupoid of partially invertible elements of M as the side groupoid.


title: Deep Learning with Domain Adaptation for Accelerated Projection-Reconstruction MR

abstract: Purpose: The radial k-space trajectory is a well-established sampling trajectory used in conjunction with magnetic resonance imaging. However, the radial k-space trajectory requires a large number of radial lines for high-resolution reconstruction. Increasing the number of radial lines causes longer acquisition time, making it more difficult for routine clinical use. On the other hand, if we reduce the number of radial lines, streaking artifact patterns are unavoidable. To solve this problem, we propose a novel deep learning approach with domain adaptation to restore high-resolution MR images from under-sampled k-space data. Methods: The proposed deep network removes the streaking artifacts from the artifact corrupted images. To address the situation given the limited available data, we propose a domain adaptation scheme that employs a pre-trained network using a large number of x-ray computed tomography (CT) or synthesized radial MR datasets, which is then fine-tuned with only a few radial MR datasets. Results: The proposed method outperforms existing compressed sensing algorithms, such as the total variation and PR-FOCUSS methods. In addition, the calculation time is several orders of magnitude faster than the total variation and PR-FOCUSS methods.Moreover, we found that pre-training using CT or MR data from similar organ data is more important than pre-training using data from the same modality for different organ. Conclusion: We demonstrate the possibility of a domain-adaptation when only a limited amount of MR data is available. The proposed method surpasses the existing compressed sensing algorithms in terms of the image quality and computation time.


title: Band splitting and Weyl nodes in trigonal tellurium studied by angle-resolved photoemission spectroscopy and density functional theory

abstract: We have performed high-resolution angle-resolved photoemission spectroscopy (ARPES) on trigonal tellurium consisting of helical chains in the crystal. Through the band-structure mapping in the three-dimensional Brillouin zone, we found a definitive evidence for the band splitting originating from the chiral nature of crystal. A direct comparison of the band dispersion between the ARPES results and the first-principles band-structure calculations suggests the presence of Weyl nodes and tiny spin-polarized hole pockets around the H point. The present result opens a pathway toward studying the interplay among crystal symmetry, band structure, and exotic physical properties in chiral crystals.


title: Model Spaces for Risk Measures

abstract: We show how risk measures originally defined in a model free framework in terms of acceptance sets and reference assets imply a meaningful underlying probability structure. Hereafter we construct a maximal domain of definition of the risk measure respecting the underlying ambiguity profile. We particularly emphasise liquidity effects and discuss the correspondence between properties of the risk measure and the structure of this domain as well as subdifferentiability properties. Keywords: Model free risk assessment, extension of risk measures, continuity properties of risk measures, subgradients.


title: Multiplicative Weights Update with Constant Step-Size in Congestion Games: Convergence, Limit Cycles and Chaos

abstract: The Multiplicative Weights Update (MWU) method is a ubiquitous meta-algorithm that works as follows: A distribution is maintained on a certain set, and at each step the probability assigned to element $\gamma$ is multiplied by $(1 -\epsilon C(\gamma))>0$ where $C(\gamma)$ is the "cost" of element $\gamma$ and then rescaled to ensure that the new values form a distribution. We analyze MWU in congestion games where agents use \textit{arbitrary admissible constants} as learning rates $\epsilon$ and prove convergence to \textit{exact Nash equilibria}. Our proof leverages a novel connection between MWU and the Baum-Welch algorithm, the standard instantiation of the Expectation-Maximization (EM) algorithm for hidden Markov models (HMM). Interestingly, this convergence result does not carry over to the nearly homologous MWU variant where at each step the probability assigned to element $\gamma$ is multiplied by $(1 -\epsilon)^{C(\gamma)}$ even for the most innocuous case of two-agent, two-strategy load balancing games, where such dynamics can provably lead to limit cycles or even chaotic behavior.


title: Inline self-diffraction dispersion-scan of over octave-spanning pulses in the single-cycle regime

abstract: We present an implementation of dispersion-scan based on self-diffraction (SD d-scan) and apply it to the measurement of over octave-spanning sub-4-fs pulses. The results are compared with second-harmonic generation (SHG) d-scan. The efficiency of the SD process is derived theoretically and compared with the spectral response retrieved by the d-scan algorithm. The new SD d-scan has a robust inline setup and enables measuring pulses with over-octave spectra, single-cycle durations and wavelength ranges beyond those of SHG crystals, such as the ultraviolet and the deep-ultraviolet.


title: Reversed interplay of quantum interference and which-way information in multi-photon entangled states

abstract: We report experimental studies of the multi-photon quantum interference of a two-mode three-photon entangled Fock state $|2, 1\rangle$ + $|1, 2\rangle$ impinging on a two-port balanced beam splitter. When the distinguishability between the two input paths is increased, we observe a reduction followed by a resurgence of the quantum interference signal. We ascribe this unusual behavior to the competition among contributions from distinct numbers of interfering photons. Our theoretical analysis shows that this phenomenon will occur for any entangled Fock-state input of the form $|N, M\rangle$ + $|M, N\rangle$ where M, N > 0. Our results are an indication that wave-particle duality may give rise to a wide range of, largely unexplored, phenomena in multi-particle interference.


title: Dynamic State Warping

abstract: The ubiquity of sequences in many domains enhances significant recent interest in sequence learning, for which a basic problem is how to measure the distance between sequences. Dynamic time warping (DTW) aligns two sequences by nonlinear local warping and returns a distance value. DTW shows superior ability in many applications, e.g. video, image, etc. However, in DTW, two points are paired essentially based on point-to-point Euclidean distance (ED) without considering the autocorrelation of sequences. Thus, points with different semantic meanings, e.g. peaks and valleys, may be matched providing their coordinate values are similar. As a result, DTW is sensitive to noise and poorly interpretable. This paper proposes an efficient and flexible sequence alignment algorithm, dynamic state warping (DSW). DSW converts each time point into a latent state, which endows point-wise autocorrelation information. Alignment is performed by using the state sequences. Thus DSW is able to yield alignment that is semantically more interpretable than that of DTW. Using one nearest neighbor classifier, DSW shows significant improvement on classification accuracy in comparison to ED (70/85 wins) and DTW (74/85 wins). We also empirically demonstrate that DSW is more robust and scales better to long sequences than ED and DTW.


title: Symmetric Laplacians, Quantum Density Matrices and their Von-Neumann Entropy

abstract: We show that the (normalized) symmetric Laplacian of a simple graph can be obtained from the partial trace over a pure bipartite quantum state that resides in a bipartite Hilbert space (one part corresponding to the vertices, the other corresponding to the edges). This suggests an interpretation of the symmetric Laplacian's Von Neumann entropy as a measure of bipartite entanglement present between the two parts of the state. We then study extreme values for a connected graph's generalized R\'enyi-$p$ entropy. Specifically, we show that (1) the complete graph achieves maximum entropy, (2) the $2$-regular graph: a) achieves minimum R\'enyi-$2$ entropy among all $k$-regular graphs, b) is within $\log 4/3$ of the minimum R\'enyi-$2$ entropy and $\log4\sqrt{2}/3$ of the minimum Von Neumann entropy among all connected graphs, c) achieves a Von Neumann entropy less than the star graph. Point $(2)$ contrasts sharply with similar work applied to (normalized) combinatorial Laplacians, where it has been shown that the star graph almost always achieves minimum Von Neumann entropy. In this work we find that the star graph achieves maximum entropy in the limit as the number of vertices grows without bound. Keywords: Symmetric; Laplacian; Quantum; Entropy; Bounds; R\'enyi.


title: Why is it hard to beat $O(n^2)$ for Longest Common Weakly Increasing Subsequence?

abstract: The Longest Common Weakly Increasing Subsequence problem (LCWIS) is a variant of the classic Longest Common Subsequence problem (LCS). Both problems can be solved with simple quadratic time algorithms. A recent line of research led to a number of matching conditional lower bounds for LCS and other related problems. However, the status of LCWIS remained open. In this paper we show that LCWIS cannot be solved in strongly subquadratic time unless the Strong Exponential Time Hypothesis (SETH) is false. The ideas which we developed can also be used to obtain a lower bound based on a safer assumption of NC-SETH, i.e. a version of SETH which talks about NC circuits instead of less expressive CNF formulas.


title: Auxetic behavior and acoustic properties of microstructured piezoelectric strain sensors

abstract: The use of multifunctional composite materials adopting piezo-electric periodic cellular lattice structures with auxetic elastic behavior is a recent and promising solution in the design of piezoelectric sensors. In the present work, periodic anti-tetrachiral auxetic lattice structures, characterized by different geometries, are taken into account and the mechanical and piezoelectrical response are investigated. The equivalent piezoelectric properties are obtained adopting a first order computational homogenization approach, generalized to the case of electro-mechanical coupling, and various polarization directions are adopted. Two examples of in-plane and out-of-plane strain sensors are proposed as design concepts. Moreover, a piezo-elasto-dynamic dispersion analysis adopting the Floquet-Bloch decomposition is performed. The acoustic behavior of the periodic piezoelectric material with auxetic topology is studied and possible band gaps are detected.


title: Coulomb interactions and screening effects in few-layer black phosphorus: a tight-binding consideration beyond the long-wavelength limit

abstract: Coulomb interaction and its screening play an important role in many physical phenomena of materials ranging from optical properties to many-body effects including superconductivity. Here, we report on a systematic study of dielectric screening in few-layer black phosphorus (BP), a two-dimensional material with promising electronic and optical characteristics. We use a combination of a tight-binding model and rigorously determined bare Coulomb interactions, which allows us to consider relevant microscopic effects beyond the long-wavelength limit. We calculate the dielectric function of few-layer BP in the random phase approximation and show that it exhibits strongly anisotropic behavior even in the static limit. We also estimate the strength of effective local and non-local Coulomb interactions and determine their doping dependence. We find that the $p_z$ states responsible for low-energy excitations in BP provide a moderate contribution to the screening, weakening the on-site Coulomb interaction by less that a factor of two. Finally, we calculate the full plasmon spectrum of few-layer BP and discuss the effects beyond long-wavelengths.


title: Grazing Incidence Small Angle X-Ray Scattering (GISAXS) on Small Targets Using Large Beams

abstract: GISAXS is often used as a versatile tool for the contactless and destruction-free investigation of nanostructured surfaces. However, due to the shallow incidence angles, the footprint of the X-ray beam is significantly elongated, limiting GISAXS to samples with typical target lengths of several millimetres. For many potential applications, the production of large target areas is impractical, and the targets are surrounded by structured areas. Because the beam footprint is larger than the targets, the surrounding structures contribute parasitic scattering, burying the target signal. In this paper, GISAXS measurements of isolated as well as surrounded grating targets in Si substrates with line lengths from $50\,{\rm\mu m}$ down to $4\,{\rm\mu m}$ are presented. For the isolated grating targets, the changes in the scattering patterns due to the reduced target length are explained. For the surrounded grating targets, the scattering signal of a $15\,{\rm\mu m}\,\times\,15\,{\rm\mu m}$ target grating structure is separated from the scattering signal of $100\,{\rm\mu m}\,\times\,100\,{\rm\mu m}$ nanostructured surroundings by producing the target with a different orientation with respect to the predominant direction of the surrounding structures. The described technique allows to apply GISAXS, e.g. for characterization of metrology fields in the semiconductor industry, where up to now it has been considered impossible to use this method due to the large beam footprint.


title: Tetramer Bound States in Heteronuclear Systems

abstract: We calculate the universal spectrum of trimer and tetramer states in heteronuclear mixtures of ultracold atoms with different masses in the vicinity of the heavy-light dimer threshold. To extract the energies, we solve the three- and four-body problem for simple two- and three-body potentials tuned to the universal region using the Gaussian expansion method. We focus on the case of one light particle of mass $m$ and two or three heavy bosons of mass $M$ with resonant heavy-light interactions. We find that trimer and tetramer cross into the heavy-light dimer threshold at almost the same point and that as the mass ratio $M/m$ decreases, the distance between the thresholds for trimer and tetramer states becomes smaller. We also comment on the possibility of observing exotic three-body states consisting of a dimer and two atoms in this region and compare with previous work.


title: Runtime Optimization of Join Location in Parallel Data Management Systems

abstract: Applications running on parallel systems often need to join a streaming relation or a stored relation with data indexed in a parallel data storage system. Some applications also compute UDFs on the joined tuples. The join can be done at the data storage nodes, corresponding to reduce side joins, or by fetching data from the storage system to compute nodes, corresponding to map side join. Both may be suboptimal: reduce side joins may cause skew, while map side joins may lead to a lot of data being transferred and replicated. In this paper, we present techniques to make runtime decisions between the two options on a per key basis, in order to improve the throughput of the join, accounting for UDF computation if any. Our techniques are based on an extended ski-rental algorithm and provide worst-case performance guarantees with respect to the optimal point in the space considered by us. Our techniques use load balancing taking into account the CPU, network and I/O costs as well as the load on compute and storage nodes. We have implemented our techniques on Hadoop, Spark and the Muppet stream processing engine. Our experiments show that our optimization techniques provide a significant improvement in throughput over existing techniques.


title: A note on minimum linear arrangement for BC graphs

abstract: A linear arrangement is a labeling or a numbering or a linear ordering of the vertices of a graph. In this paper we solve the minimum linear arrangement problem for bijective connection graphs (for short BC graphs) which include hypercubes, M\"{o}bius cubes, crossed cubes, twisted cubes, locally twisted cube, spined cube, $Z$-cubes, etc. as the subfamilies.


title: On the intersection graph of ideals of $\mathbb{Z}_m$

abstract: Let $m>1$ be an integer, and let $I(\mathbb{Z}_m)^*$ be the set of all non-zero proper ideals of $\mathbb{Z}_m$. The intersection graph of ideals of $\mathbb{Z}_m$, denoted by $G(\mathbb{Z}_m)$, is a graph with vertices $I(\mathbb{Z}_m)^*$ and two distinct vertices $I,J\in I(\mathbb{Z}_m)^*$ are adjacent if and only if $I\cap J\neq 0$. Let $n>1$ be an integer and $\mathbb{Z}_n$ be a $\mathbb{Z}_m$-module. In this paper, we introduce and study a kind of graph structure of $\mathbb{Z}_m$, denoted by $G_n(\mathbb{Z}_m)$. It is the undirected graph with the vertex set $I(\mathbb{Z}_m)^*$, and two distinct vertices $I$ and $J$ are adjacent if and only if $I\mathbb{Z}_n\cap J\mathbb{Z}_n\neq 0$. Clearly, $G_m(\mathbb{Z}_m)=G(\mathbb{Z}_m)$. We obtain some graph theoretical properties of $G_n(\mathbb{Z}_m)$ and we compute some of its numerical invariants, namely girth, independence number, domination number, maximum degree and chromatic index. We also determine all integer numbers $n$ and $m$ for which $G_n(\mathbb{Z}_m)$ is Eulerian.


title: Charge exchange in the ultraviolet: implication for interacting clouds in the core of NGC 1275

abstract: Charge exchange emission is known to provide a key diagnostic to the interface between hot and cold matter in many astrophysical environments. Most of the recent charge exchange studies focus on its emission in the X-ray band, but few on the UV part, although the latter can also provide a powerful probe of the charge exchange process. An atomic calculation, as well as an application to observed data, are presented to explore and describe the potential use of the UV data for the study of cosmic charge exchange. Using the newest charge exchange model in the SPEX code v3.03, we re-analyze an archival Hubble STIS data of the central region of NGC 1275. The NGC 1275 spectrum shows hints for three possible weak lines at about 1223.6~{\AA}, 1242.4~{\AA}, and 1244.0~{\AA}, each with a significance of about $2-3\sigma$. The putative features are best explained by charge exchange between highly ionized hydrogen, neon, and sulfur with neutral matter. The wavelengths of the charge exchange lines are found robustly with uncertainties $\leq 0.3$~{\AA}. The possible charge exchange emission shows a line-of-sight velocity offset of about $-3400$ km s$^{-1}$ with respect to the NGC 1275 nucleus, which resembles one of the Ly$\alpha$ absorbers reported in Baum et al. (2005). This indicates that the charge exchange lines might be emitted as the same position of the absorber, which could be ascribed to outflowing gas from the nucleus.

categories: astro-ph.GA astro-ph.CO astro-ph.HE physics.atom-ph


title: Equivalence of Lattice Orbit Polytopes

abstract: Let $G$ be a finite permutation group acting on $\mathbb{R}^d$ by permuting coordinates. A core point (for $G$) is an integral vector $z\in \mathbb{Z}^d$ such that the convex hull of the orbit $Gz$ contains no other integral vectors but those in the orbit $Gz$. Herr, Rehn and Sch\"urmann considered the question for which groups there are infinitely many core points up to translation equivalence, that is, up to translation by vectors fixed by the group. In the present paper, we propose a coarser equivalence relation for core points called normalizer equivalence. These equivalence classes often contain infinitely many vectors up to translation, for example when the group admits an irrational invariant subspace or an invariant irreducible subspace occurring with multiplicity greater than $1$. We also show that the number of core points up to normalizer equivalence is finite if $G$ is a so-called QI-group. These groups include all transitive permutation groups of prime degree. We give an example to show how the concept of normalizer equivalence can be used to simplify integer convex optimization problems.


title: Nucleon structure functions from lattice operator product expansion

abstract: Deep-inelastic scattering, in the laboratory and on the lattice, is most instructive for understanding how the nucleon is built from quarks and gluons. The long-term goal is to compute the associated structure functions from first principles. So far this has been limited to model calculations. In this Letter we propose a new method to compute the structure functions directly from the virtual, all-encompassing Compton amplitude, utilizing the operator product expansion. This overcomes issues of renormalization and operator mixing, which so far have hindered lattice calculations of power corrections and higher moments.


title: New zirconium hydrides predicted by structure search method based on first principles calculations

abstract: The formation of precipitated zirconium (Zr) hydrides is closely related to the hydrogen embrittlement problem for the cladding materials of pressured water reactors (PWR). In this work, we systematically investigated the crystal structures of zirconium hydride (ZrHx) with different hydrogen concentrations (x = 0~2, atomic ratio) by combining the basin hopping algorithm with first principles calculations. We conclude that the P3m1 {\zeta}-ZrH0.5 is dynamically unstable, while a novel dynamically stable P3m1 ZrH0.5 structure was discovered in the structure search. The stability of bistable P42/nnm ZrH1.5 structures and I4/mmm ZrH2 structures are also revisited. We find that the P42/nnm (c/a > 1) ZrH1.5 is dynamically unstable, while the I4/mmm (c/a = 1.57) ZrH2 is dynamically stable.The P42/nnm (c/a < 1) ZrH1.5 might be a key intermediate phase for the transition of {\gamma}->{\delta}->{\epsilon} phases. Additionally, by using the thermal dynamic simulations, we find that {\delta}-ZrH1.5 is the most stable structure at high temperature while ZrH2 is the most stable hydride at low temperature. Slow cooling process will promote the formation of {\delta}-ZrH1.5, and fast cooling process will promote the formation of {\gamma}-ZrH. These results may help to understand the phase transitions of zirconium hydrides.


title: Improper multiferroicity and colossal dielectric constants in Bi$_{2}$CuO$_{4}$

abstract: The layered cuprate Bi$_{2}$CuO$_{4}$ is investigated using magnetic, dielectric and pyroelectric measurements. This system is observed to be an improper multiferroic, with a robust ferroelectric state being established near the magnetic transition. Magnetic and dielectric measurements indicate the presence of a region above the antiferromagnetic Neel temperature with concomitant polar and magnetic short range order. Bi$_{2}$CuO$_{4}$ is also seen to exhibit colossal dielectric constants at higher temperatures with clearly distinguishable grain and grain boundary contributions, both of which exhibit non-Debye relaxation.


title: Quantum Surface and Intertwiner Dynamics in Loop Quantum Gravity

abstract: We introduce simple generic models of surface dynamics in loop quantum gravity (LQG). A quantum surface is defined as a set of elementary patches of area glued together. We provide it with an extra structure of locality (nearest neighbors), thought of as induced by the whole spin network state defining the 3d bulk geometry around the quantum surface. Here, we focus on classical surface dynamics, using a spinorial description of surface degrees of freedom. We introduce two classes of dynamics, to be thought as templates for future investigation of LQG dynamics with in mind the dynamics of quantum black holes. The first defines global dynamics of the closure defect of the surface, with two basic toy-models, either a dissipative dynamics relaxing towards the closure constraint or a Hamiltonian dynamics precessing the closure defect. The second class of dynamics describes the isolated regime, when both area and closure defect are conserved throughout the evolution. The surface dynamics is implemented through U(N) transformations and generalizes to a Bose-Hubbard Hamiltonian with a local quadratic potential interaction. We briefly discuss the implications of modeling the quantum black hole dynamics by a surface Bose-Hubbard model.


title: A free boundary optimization problem for the $\infty$-Laplacian

abstract: We study a free boundary optimization problem in heat conduction, ruled by the infinity-Laplace operator, with lower temperature bound and a volume constraint. We obtain existence and regularity results and derive geometric properties for the solution and the free boundaries.


title: Compactly supported linearised observables in single-field inflation

abstract: We investigate the gauge-invariant observables constructed by smearing the graviton and inflaton fields by compactly supported tensors at linear order in general single-field inflation. These observables correspond to gauge-invariant quantities that can be measured locally. In particular, we show that these observables are equivalent to (smeared) local gauge-invariant observables such as the linearised Weyl tensor, which have better infrared properties than the graviton and inflaton fields. Special cases include the equivalence between the compactly supported gauge-invariant graviton observable and the smeared linearised Weyl tensor in Minkowski and de Sitter spaces. Our results indicate that the infrared divergences in the tensor and scalar perturbations in single-field inflation have the same status as in de Sitter space and are both a gauge artefact, in a certain technical sense, at tree level.


title: Manipulation of entanglement sudden death in an all-optical setup

abstract: The unavoidable and irreversible interaction between an entangled quantum system and its environment causes decoherence of the individual qubits as well as degradation of the entanglement between them. Entanglement sudden death (ESD) is the phenomenon wherein disentanglement happens in finite time even when individual qubits decohere only asymptotically in time due to noise. Prolonging the entanglement is essential for the practical realization of entanglement-based quantum information and computation protocols. For this purpose, the local NOT operation in the computational basis on one or both qubits has been proposed. Here, we formulate an all-optical experimental set-up involving such NOT operations that can hasten, delay, or completely avert ESD, all depending on when it is applied during the process of decoherence. Analytical expressions for these are derived in terms of parameters of the initial state's density matrix, whether for pure or mixed entangled states. After a discussion of the schematics of the experiment, the problem is theoretically analyzed, and simulation results of such manipulations of ESD are presented.


title: A nonlinear Monte Carlo model of super-diffusive shock acceleration with magnetic field amplification

abstract: Fast collisionless shocks in cosmic plasmas convert their kinetic energy flow into the hot downstream thermal plasma with a substantial fraction of energy going into a broad spectrum of superthermal charged particles and magnetic fluctuations. The superthermal particles can penetrate into the shock upstream region producing an extended shock precursor. The cold upstream plasma flow is decelerated by the force provided by the superthermal particle pressure gradient. In high Mach number collisionless shocks, efficient particle acceleration is likely coupled with turbulent magnetic field amplification (MFA) generated by the anisotropic distribution of accelerated particles. This anisotropy is determined by the fast particle transport making the problem strongly nonlinear and multi-scale. Here, we present a nonlinear Monte Carlo model of collisionless shock structure with super-diffusive propagation of high-energy Fermi accelerated particles coupled to particle acceleration and MFA which affords a consistent description of strong shocks. A distinctive feature of the Monte Carlo technique is that it includes the full angular anisotropy of the particle distribution at all precursor positions. The model reveals that the super-diffusive transport of energetic particles (i.e., Levy-walk propagation) generates a strong quadruple anisotropy in the precursor particle distribution. The resultant pressure anisotropy of the high-energy particles produces a non-resonant mirror-type instability which amplifies compressible wave modes with wavelengths longer than the gyroradii of the highest energy protons produced by the shock.


title: FeUdal Networks for Hierarchical Reinforcement Learning

abstract: We introduce FeUdal Networks (FuNs): a novel architecture for hierarchical reinforcement learning. Our approach is inspired by the feudal reinforcement learning proposal of Dayan and Hinton, and gains power and efficacy by decoupling end-to-end learning across multiple levels -- allowing it to utilise different resolutions of time. Our framework employs a Manager module and a Worker module. The Manager operates at a lower temporal resolution and sets abstract goals which are conveyed to and enacted by the Worker. The Worker generates primitive actions at every tick of the environment. The decoupled structure of FuN conveys several benefits -- in addition to facilitating very long timescale credit assignment it also encourages the emergence of sub-policies associated with different goals set by the Manager. These properties allow FuN to dramatically outperform a strong baseline agent on tasks that involve long-term credit assignment or memorisation. We demonstrate the performance of our proposed system on a range of tasks from the ATARI suite and also from a 3D DeepMind Lab environment.


title: Piezoelectricity in two-dimensional materials: a comparative study between lattice dynamics and ab-initio calculations

abstract: Elastic constant C_{11} and piezoelectric stress constant e_{1,11} of two-dimensional (2D) dielectric materials comprising h-BN, 2H MoS2 and other transition metal dichalcogenides (TMDCs) and -dioxides (TMDOs) are calculated using lattice dynamical theory. The results are compared with corresponding quantities obtained by ab-initio calculations. We identify the difference between clamped-ion and relaxed-ion contributions with the dependence on inner strains which are due to the relative displacements of the ions in the unit cell. Lattice dynamics allows to express the inner strains contributions in terms of microscopic quantities such as effective ionic charges and optoacoustical couplings, which allows us to clarify differences in the piezoelectric behavior between h- BN versus MoS2. Trends in the different microscopic quantities as functions of atomic composition are discussed.


title: Electronic and atomic kinetics in solids irradiated with free-electron lasers or swift-heavy ions

abstract: In this brief review we discuss the transient processes in solids under irradiation with femtosecond X-ray free-electron-laser (FEL) pulses and swift-heavy ions (SHI). Both kinds of irradiation produce highly excited electrons in a target on extremely short timescales. Transfer of the excess electronic energy into the lattice may lead to observable target modifications such as phase transitions and damage formation. Transient kinetics of material excitation and relaxation under FEL or SHI irradiation are comparatively discussed. The same origin for the electronic and atomic relaxation in both cases is demonstrated. Differences in these kinetics introduced by the geometrical effects ({\mu}m-size of a laser spot vs nm-size of an ion track) and initial irradiation (photoabsorption vs an ion impact) are analyzed. The basic mechanisms of electron transport and electron-lattice coupling are addressed. Appropriate models and their limitations are presented. Possibilities of thermal and nonthermal melting of materials under FEL and SHI irradiation are discussed.


title: Nonlinear Model Predictive Control for Multi-Micro Aerial Vehicle Robust Collision Avoidance

abstract: Multiple multirotor Micro Aerial Vehicles sharing the same airspace require a reliable and robust collision avoidance technique. In this paper we address the problem of multi-MAV reactive collision avoidance. A model-based controller is employed to achieve simultaneously reference trajectory tracking and collision avoidance. Moreover, we also account for the uncertainty of the state estimator and the other agents position and velocity uncertainties to achieve a higher degree of robustness. The proposed approach is decentralized, does not require collision-free reference trajectory and accounts for the full MAV dynamics. We validated our approach in simulation and experimentally.


title: A simple thermodynamic model for the hydrogen phase diagram

abstract: We describe a classical thermodynamic model that reproduces the main features of the solid hydrogen phase diagram. In particular, we show how the general structure types that are found by electronic structure calculations and the quantum nature of the protons can also be understood from a classical viewpoint. The model provides a picture not only of crystal structure, but also for the anomalous melting curve and insights into isotope effects, liquid metallisation and InfraRed activity. The existence of a classical picture for this most quantum of condensed matter systems provides a surprising extension of the correspondence principle of quantum mechanics, in particular the equivalent effects of classical and quantum uncertainty.


title: Give Me Some Slack: Efficient Network Measurements

abstract: Many networking applications require timely access to recent network measurements, which can be captured using a sliding window model. Maintaining such measurements is a challenging task due to the fast line speed and scarcity of fast memory in routers. In this work, we study the impact of allowing \emph{slack} in the window size on the asymptotic requirements of sliding window problems. That is, the algorithm can dynamically adjust the window size between $W$ and $W(1+\tau)$ where $\tau$ is a small positive parameter. We demonstrate this model's attractiveness by showing that it enables efficient algorithms to problems such as MAX and GENERAL-SUM that require $\Omega(W)$ bits even for constant factor approximations in the exact sliding window model. Additionally, for problems that admit sub-linear approximation algorithms such as BASIC-SUMMING and COUNT-DISTINCT, the slack model enables a further asymptotic improvement.


title: A compact Verilog-A ReRAM switching model

abstract: The translation of emerging application concepts that exploit Resistive Random Access Memory (ReRAM) into large-scale practical systems requires realistic, yet computationally efficient, empirical models that can capture all observed physical devices. Here, we present a Verilog-A ReRAM model built upon experimental routines performed on TiOx-based prototypes. This model was based on custom biasing protocols, specifically designed to reveal device switching rate dependencies on a) bias voltage and b) initial resistive state. Our model is based on the assumption that a stationary switching rate surface m(R,v) exists for sufficiently low voltage stimulation. The proposed model comes in compact form as it is expressed by a simple voltage dependent exponential function multiplied with a voltage and initial resistive state dependent second order polynomial expression, which makes it suitable for fast and/or large-scale simulations.


title: Sum-set Inequalities from Aligned Image Sets: Instruments for Robust GDoF Bounds

abstract: We present sum-set inequalities specialized to the generalized degrees of freedom (GDoF) framework. These are information theoretic lower bounds on the entropy of bounded density linear combinations of discrete, power-limited dependent random variables in terms of the joint entropies of arbitrary linear combinations of new random variables that are obtained by power level partitioning of the original random variables. These bounds generalize the aligned image sets approach, and are useful instruments to obtain GDoF characterizations for wireless networks, especially with multiple antenna nodes, subject to arbitrary channel strength and channel uncertainty levels. To demonstrate the utility of these bounds, we consider a non-trivial instance of wireless networks - a two user interference channel with different number of antennas at each node, and different levels of partial channel knowledge available to the transmitters. We obtain tight GDoF characterization for specific instance of this channel with the aid of sum-set inequalities.


title: Preduals for spaces of operators involving Hilbert spaces and trace-class operators

abstract: Continuing the study of preduals of spaces $\mathcal{L}(H,Y)$ of bounded, linear maps, we consider the situation that $H$ is a Hilbert space. We establish a natural correspondence between isometric preduals of $\mathcal{L}(H,Y)$ and isometric preduals of $Y$. The main ingredient is a Tomiyama-type result which shows that every contractive projection that complements $\mathcal{L}(H,Y)$ in its bidual is automatically a right $\mathcal{L}(H)$-module map. As an application, we show that isometric preduals of $\mathcal{L}(\mathcal{S}_1)$, the algebra of operators on the space of trace-class operators, correspond to isometric preduals of $\mathcal{S}_1$ itself (and there is an abundance of them). On the other hand, the compact operators are the unique predual of $\mathcal{S}_1$ making its multiplication separately weak* continuous.


title: A Survey on Content-Aware Video Analysis for Sports

abstract: Sports data analysis is becoming increasingly large-scale, diversified, and shared, but difficulty persists in rapidly accessing the most crucial information. Previous surveys have focused on the methodologies of sports video analysis from the spatiotemporal viewpoint instead of a content-based viewpoint, and few of these studies have considered semantics. This study develops a deeper interpretation of content-aware sports video analysis by examining the insight offered by research into the structure of content under different scenarios. On the basis of this insight, we provide an overview of the themes particularly relevant to the research on content-aware systems for broadcast sports. Specifically, we focus on the video content analysis techniques applied in sportscasts over the past decade from the perspectives of fundamentals and general review, a content hierarchical model, and trends and challenges. Content-aware analysis methods are discussed with respect to object-, event-, and context-oriented groups. In each group, the gap between sensation and content excitement must be bridged using proper strategies. In this regard, a content-aware approach is required to determine user demands. Finally, the paper summarizes the future trends and challenges for sports video analysis. We believe that our findings can advance the field of research on content-aware video analysis for broadcast sports.


title: On the equivalence between bumblebee models and electrodynamics in a non-linear gauge

abstract: Bumblebee models are effective field theories describing a vector field with a nonzero vacuum expectation value that spontaneously breaks Lorentz invariance. They provide an alternative way of exploring the similarities between theories with spontaneous Lorentz symmetry breaking and gauge theories. The equivalence between bumblebee models with suitable conditions and standard electrodynamics in a non-linear gauge $A_\mu A^\mu-b^2=0$ is taken for granted; however, this point is very subtle and has not yet been fully addressed. The main goal of this paper is to fill in this gap. More precisely, here we study the relation between a bumblebee model, with a smooth potential of the form $V(B _{\mu}) = V (B _{\mu} B ^{\mu} - b ^{2})$, and standard electrodynamics in the non-linear gauge $A _{\mu} A ^{\mu} - b ^{2} = 0$, both at the classical and quantum levels. Using the Dirac's method we show that after introducing Dirac brackets with suitable initial conditions, the classical dynamics of the bumblebee model corresponds to that of standard electrodynamics in the aforementioned non-linear gauge. In the quantum case we demonstrate that perturbative calculations of Feynman amplitudes to any physical process in each model are indistinguishable. To do this, we show that the Feynman rules and propagators of standard electrodynamics in the non-linear gauge and those describing the bumblebee model are the same.


title: Room temperature zero field skyrmions in Fe-based thin film stacks

abstract: A new paradigm is required to facilitate the demand for the huge data storage capacity and faster data processing in the future. Nano structures such as magnetic skyrmions have been predicted to address this issues as these vortex structures are the smallest particle-like magnetic features and are topologically protected from crystallographic defects or magnetic disorder. We report here stable skyrmions at room temperature and with zero applied field in ebeam evaporated Ir Fe Ir MgO thin film stacks. Micromagnetic simulations show that these skyrmions are induced by interfacial Dzyaloshinskii Moriya interactions between the ferromagnetic and heavy metal ultra thin layers values in the range between 3.1 to 3.6. We also show the field dependencies of the skyrmion width and density. The room temperature zero field width is 110 nm, approximately comparable to the value showed by Fert et al and Boulle et al for the multi stacks thin film. DFT calculations corroborate our experimental results by predicting a DMI value of 3.67.


title: Vibronic exciton theory of singlet fission. I. Linear absorption and the anatomy of the correlated triplet pair state

abstract: Recent time-resolved spectroscopic experiments have indicated that vibronic coupling plays a vital role in facilitating the process of singlet fission. In this work, which forms the first article of a series, we set out to unravel the mechanisms underlying singlet fission through a vibronic exciton theory. We formulate a model in which both electronic and vibrational degrees of freedom are treated microscopically and non-perturbatively. Using pentacene as a prototypical material for singlet fission, we subject our theory to comparison with measurements on polarization-resolved absorption of single crystals, and employ our model to characterize the excited states underlying the absorption band. Special attention is given to convergence of photophysical observables with respect to the basis size employed, through which we determine the optimal basis for more expensive calculations to be presented in subsequent work. We furthermore evaluate the energetic separation between the optically prepared singlet excited state and the correlated triplet pair state, as well as provide a real-space characterization of the latter, both of which are of key importance in the discussion of fission dynamics. We discuss our results in the context of recent experimental studies.


title: Vibronic exciton theory of singlet fission. II. Two-dimensional spectroscopic detection of the correlated triplet pair state

abstract: Singlet fission, the molecular process through which photons are effectively converted into pairs of lower energy triplet excitons, holds promise as a means of boosting photovoltaic device efficiencies. In the preceding article of this series, we formulated a vibronic theory of singlet fission, inspired by previous experimental and theoretical studies suggesting that vibronic coupling plays an important role in fission dynamics. Here, we extend our model in order to simulate two-dimensional electronic spectra, through which the theory is further validated based on a comparison to recent measurements on pentacene crystals. Moreover, by means of such spectral simulations, we provide new insights into the nature of the correlated triplet pair state, the first product intermediate in the fission process. In particular, we address a disagreement in the literature regarding the identification, energies, and transition dipole moments of its optical transitions towards higher-lying triplet states.


title: $O(N)$ Iterative and $O(NlogN)$ Fast Direct Volume Integral Equation Solvers with a Minimal-Rank ${\cal H}^2$-Representation for Large-Scale $3$-D Electrodynamic Analysis

abstract: Linear complexity iterative and log-linear complexity direct solvers are developed for the volume integral equation (VIE) based general large-scale electrodynamic analysis. The dense VIE system matrix is first represented by a new cluster-based multilevel low-rank representation. In this representation, all the admissible blocks associated with a single cluster are grouped together and represented by a single low-rank block, whose rank is minimized based on prescribed accuracy. From such an initial representation, an efficient algorithm is developed to generate a minimal-rank ${\cal H}^2$-matrix representation. This representation facilitates faster computation, and ensures the same minimal rank's growth rate with electrical size as evaluated from singular value decomposition. Taking into account the rank's growth with electrical size, we develop linear-complexity ${\cal H}^2$-matrix-based storage and matrix-vector multiplication, and thereby an $O(N)$ iterative VIE solver regardless of electrical size. Moreover, we develop an $O(NlogN)$ matrix inversion, and hence a fast $O(NlogN)$ \emph{direct} VIE solver for large-scale electrodynamic analysis. Both theoretical analysis and numerical simulations of large-scale $1$-, $2$- and $3$-D structures on a single-core CPU, resulting in millions of unknowns, have demonstrated the low complexity and superior performance of the proposed VIE electrodynamic solvers. %The algorithms developed in this work are kernel-independent, and hence applicable to other IE operators as well.


title: Vanishing Viscosity Limit for Incompressible Viscoelasticity in Two Dimensions

abstract: This paper studies the inviscid limit of the two-dimensional incompressible viscoelasticity, which is a system coupling a Navier-Stokes equation with a transport equation for the deformation tensor. The existence of global smooth solutions near the equilibrium with a fixed positive viscosity was known since the work of F. H. Lin, C. Liu, and P. Zhang in "On hydrodynamics of viscoelastic fluids". The inviscid case was solved recently by the second author Z. Lei. in "Global well-posedness of incompressible elastodynamics in two dimensions". While the latter was solely based on the techniques from the studies of hyperbolic equations, and hence the 2D problem is in general more challenge than that in higher dimensions, the former was relied crucially upon a dissipative mechanism. Indeed, after a symmetrization and a linearization around the equilibrium, the system of the incompressible viscoelasticity reduces to an incompressible system of damped wave equations for both the fluid velocity and the deformation tensor. These two approaches are not compatible. In this paper, we prove global existence of solutions, uniformly in both time $t \in [0, \infty)$ and viscosity $\mu \geq 0$. This allows us to justify in particular the vanishing viscosity limit for all time. In order to overcome difficulties coming from the incompatibility between the purely hyperbolic limiting system and the systems with additional parabolic viscous perturbations, we introduce in this paper a rather robust method which may apply to a wide class of physical systems of similar nature. Roughly speaking, the method works in two dimensional case whenever the hyperbolic system satisfies intrinsically a "Strong Null Condition". For dimensions not less than three, the usual null condition is sufficient for this method to work.


title: Geometric unification theory of the grand unification and gravitational interactions and new physics

abstract: This paper discovers geometric unification theory of the grand unification and gravitational interactions and their new physics according to the general fiber bundle theory, symmetry and so on. Consequently, the research of this paper is based on the exact scientific bases of mathematics and physics. The Lagrangians of the grand unification and gravitational interactions are unifiedly deduced from quantitative causal principle (QCP) and satisfy the gauge invariant principle of general gauge fields interacting with Fermion and/or boson fields. The geometry and physics meanings of gauge invariant property of different physical systems are revealed, and it is discovered that all the Lagrangians of the four fundamental physics interactions are composed of the invariant quantities in corresponding spacetime structures. The difficulties that fundamental physics interactions and Noether theorem are not able to be unifiedly given and investigated are overcome, the geometric unification theory and origins of the four fundamental physics interactions and Noether theorem are shown by QCP, their two-order general Euler-Lagrange Equations and corresponding Noether conservation currents are derived in general curved spacetime. This paper further deduces QCP from symmetric principle. Consequently, geometric unification theory of the grand unification and gravitation theories and Noether theorem based on symmetric principle and the new physics are given in this paper. This paper further gives the unification of QCP and symmetric principle. Thus, this paper opens a door to both study and give new developments of geometric unification theory of physics laws, and using the new geometric unification theory, a lot of research works about different branches of physics can be anew done and expressed simpler with different symmetric characters.


title: Inner structure of ZnO microspheres fabricated via laser ablation in superfluid helium

abstract: ZnO microspheres fabricated via laser ablation in superfluid helium were found to have bubble-like voids. Even a microsphere demonstrating clear whispering gallery mode resonances in the luminescence had voids. Our analysis confirmed that the voids are located away from the surface and have negligible or little effect on the whispering gallery mode resonances since the electromagnetic energy localizes near the surface of these microspheres. The existence of the voids indicates that helium gas or any evaporated target material was present within the molten microparticles during the microsphere formation.


title: Low-acceleration dwarf galaxies as tests of quantised inertia

abstract: Dwarf satellite galaxies of the Milky Way appear to be gravitationally bound, but their stars' orbital motion seems too fast to allow this given their visible mass. This is akin to the larger-scale galaxy rotation problem. In this paper, a modification of inertia called quantised inertia or MiHsC (Modified inertia due to a Hubble-scale Casimir effect) which correctly predicts larger galaxy rotations without dark matter is tested on eleven dwarf satellite galaxies of the Milky Way, for which mass and velocity data are available. Quantised inertia slightly outperforms MoND (Modified Newtonian Dynamics) in predicting the velocity dispersion of these systems, and has the fundamental advantage over MoND that it does not need an adjustable parameter.


title: The Lie-Trotter integrator in the dynamics of the symmetric free rigid body

abstract: The numerical integration plays a fundamental role in understanding the behaviour of many mechanical systems. In this paper some important aspects of the mechanical integrators on the dynamics of a mechanical system are studied. More specific, we have shown that if that the Lie-Trotter integrator is obtained, in case of Euler equations for the dynamics of symmetric free rigid body, then it is a Poisson integrator. At the end of the paper some important remarks are presented.


title: Size Dependence of Nanoscale Wear of Silicon Carbide

abstract: Nanoscale, single-asperity wear of single-crystal silicon carbide (sc-SiC) and nanocrystalline silicon carbide (nc-SiC) is investigated using single-crystal diamond nanoindenter tips and nanocrystalline diamond atomic force microscopy (AFM) tips under dry conditions, and the wear behavior is compared to that of single-crystal silicon with both thin and thick native oxide layers. We discovered a transition in the relative wear resistance of the SiC samples compared to that of Si as a function of contact size. With larger nanoindenter tips (tip radius around 370 nm), the wear resistances of both sc-SiC and nc-SiC are higher than that of Si. This result is expected from the Archard's equation because SiC is harder than Si. However, with the smaller AFM tips (tip radius around 20 nm), the wear resistances of sc-SiC and nc-SiC are lower than that of Si, despite the fact that the contact pressures are comparable to those applied with the nanoindenter tips, and the plastic zones are well-developed in both sets of wear experiments. We attribute the decrease in the relative wear resistance of SiC compared to that of Si to a transition from a wear regime dominated by the materials' resistance to plastic deformation (i.e., hardness) to a regime dominated by the materials' resistance to interfacial shear. This conclusion is supported by our AFM studies of wearless friction, which reveal that the interfacial shear strength of SiC is higher than that of Si. The contributions of surface roughness and surface chemistry to differences in interfacial shear strength are also discussed.


title: Two dimensional ice mapping of molecular cores

abstract: We present maps of the column densities of H$_2$O, CO$_2$, and CO ices towards the molecular cores B~35A, DC~274.2-00.4, BHR~59, and DC~300.7-01.0. These ice maps, probing spatial distances in molecular cores as low as 2200~AU, challenge the traditional hypothesis that the denser the region observed, the more ice is present, providing evidence that the relationships between solid molecular species are more varied than the generic picture we often adopt to model gas-grain chemical processes and explain feedback between solid phase processes and gas phase abundances. We present the first combined solid-gas maps of a single molecular species, based upon observations of both CO ice and gas phase C$^{18}$O towards B~35A, a star-forming dense core in Orion. We conclude that molecular species in the solid phase are powerful tracers of "small scale" chemical diversity, prior to the onset of star formation. With a component analysis approach, we can probe the solid phase chemistry of a region at a level of detail greater than that provided by statistical analyses or generic conclusions drawn from single pointing line-of-sight observations alone.


title: Radial metal abundance profiles in the intra-cluster medium of cool-core galaxy clusters, groups, and ellipticals

abstract: The hot intra-cluster medium (ICM) permeating galaxy clusters and groups is not pristine, as it is continuously enriched by metals synthesised in Type Ia (SNIa) and core-collapse (SNcc) supernovae since the major epoch of star formation (z ~ 2-3). The cluster/group enrichment history and the mechanisms responsible for releasing and mixing the metals can be probed via the radial distribution of SNIa and SNcc products within the ICM. In this paper, we use deep XMM-Newton/EPIC observations from a sample of 44 nearby cool-core galaxy clusters, groups, and ellipticals (CHEERS) to constrain the average radial O, Mg, Si, S, Ar, Ca, Fe, and Ni abundance profiles. The radial distributions of all these elements, averaged over a large sample for the first time, represent the best constrained profiles available currently. We find an overall decrease of the Fe abundance with radius out to ~$0.9 r_{500}$ and ~$0.6 r_{500}$ for clusters and groups, respectively, in good agreement with predictions from the most recent hydrodynamical simulations. The average radial profiles of all the other elements (X) are also centrally peaked and, when rescaled to their average central X/Fe ratios, follow well the Fe profile out to at least ~0.5$r_{500}$. Using two sets of SNIa and SNcc yield models reproducing well the X/Fe abundance pattern in the core, we find that, as predicted by recent simulations, the relative contribution of SNIa (SNcc) to the total ICM enrichment is consistent with being uniform at all radii, both for clusters and groups. In addition to implying that the central metal peak is balanced between SNIa and SNcc, our results suggest that the enriching SNIa and SNcc products must share the same origin, and that the delay between the bulk of the SNIa and SNcc explosions must be shorter than the timescale necessary to diffuse out the metals.


title: Testing isotropy in the Two Micron All-Sky redshift survey with information entropy

abstract: We use information entropy to test the isotropy in the nearby galaxy distribution mapped by the Two Micron All-Sky redshift survey (2MRS). We find that the galaxy distribution is highly anisotropic on small scales. The radial anisotropy gradually decreases with increasing length scales and the observed anisotropy is consistent with that expected for an isotropic Poisson distribution beyond a length scale of $90 \, h^{-1}\, {\rm Mpc}$. Using mock catalogues from N-body simulations, we find that the galaxy distribution in the 2MRS exhibits a degree of anisotropy compatible with that of the $\Lambda$CDM model after accounting for the clustering bias of the 2MRS galaxies. We also quantify the polar and azimuthal anisotropies and identify two directions $(l,b)=(150^{\circ}, -15^{\circ})$, $(l,b)=(310^{\circ},-15^{\circ})$ which are significantly anisotropic compared to the other directions in the sky. We suggest that their preferential orientations on the sky may indicate a possible alignment of the Local Group with two nearby large scale structures. Despite the differences in the degree of anisotropy on small scales, we find that the galaxy distributions in both the 2MRS and the $\Lambda$CDM model are isotropic on a scale of $90 \, h^{-1}\, {\rm Mpc}$.


title: Synchrotron spectrum of fast cooling electrons in GRBs

abstract: We discuss the synchrotron emission of fast cooling electrons in shocks. The fast cooling electrons behind the shocks can generate a position-dependent inhomogeneous electron distribution if they do not have enough time to mix homogeneously. This would lead to a very different synchrotron spectrum in low-frequency bands to that in the homogeneous case due to the synchrotron absorption. In this paper, we calculate the synchrotron spectrum in the inhomogeneous case in a gamma-ray burst (GRB). Both the forward shock and the reverse shock are considered. We find for the reverse shock dominated case, we would expect a "reverse shock bump" in the low-frequency spectrum. The spectral bump is due to the combination synchrotron absorption in both the forward and reverse shock regions. In the forward shock spectrum in the low frequencies has two unconventional segments with spectral slopes of $\lesssim1$ and $11/8$. The slope of $11/8$ has been found by some authors, while the slope of $\lesssim1$ is new, which is due to the approximately constant electron temperature in the optically thick region. In the future, simultaneous observations in multiple bands (especially in the low frequency bands) in the GRB early afterglow or prompt emission phases will possibly reveal these spectral characteristics and enable us to identify the reverse shock component and distinguish between the forward and reverse shock emissions. This also may be a method with which to diagnose the electron distribution status (homogeneous or inhomogeneous) after fast cooling in the relativistic shock region.


title: The frozen nucleon approximation in two-particle two-hole response functions

abstract: We present a fast and efficient method to compute the inclusive two-particle two-hole (2p-2h) electroweak responses in the neutrino and electron quasielastic inclusive cross sections. The method is based on two approximations. The first neglects the motion of the two initial nucleons below the Fermi momentum, which are considered to be at rest. This approximation, which is reasonable for high values of the momentum transfer, turns out also to be quite good for moderate values of the momentum transfer $q\gtrsim k_F$. The second approximation involves using in the "frozen" meson-exchange currents (MEC) an effective $\Delta$-propagator averaged over the Fermi sea. Within the resulting "frozen nucleon approximation", the inclusive 2p-2h responses are accurately calculated with only a one-dimensional integral over the emission angle of one of the final nucleons, thus drastically simplifying the calculation and reducing the computational time. The latter makes this method especially well-suited for implementation in Monte Carlo neutrino event generators.


title: Weighted Growing Simplicial Complexes

abstract: Simplicial complexes describe collaboration networks, protein interaction networks and brain networks and in general network structures in which the interactions can include more than two nodes. In real applications, often simplicial complexes are weighted. Here we propose a non-equilibrium model for weighted growing simplicial complexes. The proposed dynamics is able to generate weighted simplicial complexes with a rich interplay between weights and topology emerging not just at the level of nodes and links, but also at the level of faces of higher dimension.


title: Arrested relaxation in an isolated molecular ultracold plasma

abstract: Spontaneous avalanche to plasma splits the core of an ellipsoidal Rydberg gas of nitric oxide. Ambipolar expansion first quenches the electron temperature of this core plasma. Then, long-range, resonant charge transfer from ballistic ions to frozen Rydberg molecules in the wings of the ellipsoid quenches the centre-of-mass ion/Rydberg molecule velocity distribution. This sequence of steps gives rise to a remarkable mechanics of self-assembly, in which the kinetic energy of initially formed hot electrons and ions drives an observed separation of plasma volumes. These dynamics adiabatically sequester energy in a reservoir of mass transport, starting a process that anneals separating volumes to form an apparent glass of strongly coupled ions and electrons. Short-time electron spectroscopy provides experimental evidence for complete ionization. The long lifetime of this system, particularly its stability with respect to recombination and neutral dissociation, suggests that this transformation affords a robust state of arrested relaxation, far from thermal equilibrium.

categories: physics.plasm-ph cond-mat.dis-nn cond-mat.mes-hall physics.atom-ph


title: Periodic and quasi-periodic attractors for the spin-orbit evolution of Mercury with a realistic tidal torque

abstract: Mercury is entrapped in a 3:2 resonance: it rotates on its axis three times for every two revolutions it makes around the Sun. It is generally accepted that this is due to the large value of the eccentricity of its orbit. However, the mathematical model originally introduced to study its spin-orbit evolution proved not to be entirely convincing, because of the expression commonly used for the tidal torque. Only recently, in a series of papers mainly by Efroimsky and Makarov, a different model for the tidal torque has been proposed, which has the advantages of being more realistic, and of providing a higher probability of capture in the 3:2 resonance with respect to the previous models. On the other hand, a drawback of the model is that the function describing the tidal torque is not smooth and consists of a superposition of kinks, so that both analytical and numerical computations turn out to be rather delicate: indeed, standard perturbation theory based on power series expansion cannot be applied and the implementation of a fast algorithm to integrate the equations of motion numerically requires a high degree of care. In this paper, we make a detailed study of the spin-orbit dynamics of Mercury, as predicted by the realistic model: In particular, we present numerical and analytical results about the nature of the librations of Mercury's spin in the 3:2 resonance. The results provide evidence that the librations are quasi-periodic in time.

categories: math.DS astro-ph.EP math-ph math.CA math.MP


title: Multicast Transmissions in Directional mmWave Communications

abstract: Multicast transmissions have been widely analyzed in traditional networks as a way to improve spectrum efficiency when multiple users are interested in the same data. However, their application to mmWave communications has been studied only marginally so far. The goal of this paper is to partially fill this gap by investigating optimal and suboptimal multicast schemes for mmWave communications with directional beams. In particular, we propose a Markov setup to model the retransmission status of the unsuccessfully transmitted packets and, because of the computational complexity of the optimal solution, we introduce a suboptimal hierarchical optimization procedure, which is much easier to derive. Finally, we numerically show that restricting the link to unicast beams is strongly suboptimal, especially when many packets have to be transmitted.


title: Lifting of curves with automorphisms

abstract: The lifting problem for curves with automorphisms asks whether we can lift a smooth projective characteristic p curve with a group G of automorphisms to characteristic zero. This was solved by Grothendieck when G acts with prime-to-p stabilizers, and there has been much progress over the last few decades in the wild case. We survey the techniques and obstructions for this lifting problem, aiming at a reader whose background is limited to scheme theory at the level of Hartshorne's book. Throughout, we include numerous examples and clarifying remarks. We also provide a list of open questions.


title: Tree Notation: an antifragile program notation

abstract: This paper presents Tree Notation, a new simple, universal syntax. Language designers can invent new programming languages, called Tree Languages, on top of Tree Notation. Tree Languages have a number of advantages over traditional programming languages. We include a Visual Abstract to succinctly display the problem and discovery. Then we describe the problem--the BNF to abstract syntax tree (AST) parse step--and introduce the novel solution we discovered: a new family of 2D programming languages that are written directly as geometric trees.


title: Evidence that particle acceleration in hotspots of FR II galaxies is not constrained by synchrotron cooling

abstract: We study the hotspots of powerful radiogalaxies, where electrons accelerated at the jet termination shock emit synchrotron radiation. The turnover of the synchrotron spectrum is typically observed between infrared and optical frequencies, indicating that the maximum energy of non-thermal electrons accelerated at the shock is ~TeV for a canonical magnetic field of ~100 micro Gauss. We show that this maximum energy cannot be constrained by synchrotron losses as usually assumed, unless the jet density is unreasonably large and most of the jet upstream energy goes to non-thermal particles. We test this result by considering a sample of hotspots observed at radio, infrared and optical wavelengths.


title: Highly excited electronic image states of metallic nanorings

abstract: We study electronic image states around a metallic nanoring and show that the interplay between the attractive polarization force and a repulsive centrifugal force gives rise to Rydberg-like image states trapped several nanometers away from the surface. The nanoring is modeled as a perfectly conducting isolated torus whose electrostatic image potential is derived analytically. The image states are computed via a two-dimensional finite-difference scheme as solutions of the effective Schr\"odinger equation describing the outer electron subject to this image potential. These findings demonstrate not only the existence of detached image states around nanorings but allow us also to provide general criteria on the ring geometry, i.e. the aspect ratio of the torus, that need to be fulfilled in order to support such states.


title: Why Smart Appliances May Result in a Stupid Energy Grid?

abstract: This article discusses unexpected consequences of idealistic conceptions about the modernization of power grids. We will focus our analysis on demand-response policies based on automatic decisions by the so-called smart home appliances. Following the usual design approach, each individual appliance has access to a universal signal (e.g. grid frequency or electricity price) that is believed to indicate the system state. Such information is then used as the basis of the appliances' individual decisions. While each single device has a negligible impact in the system, the aggregate effect of the distributed appliances' reactions is expect to bring improvements in the system efficiency; this effect is the demand-response policy goal. The smartness of such an ideal system, composed by isolated appliances with their individual decisions, but connected in the same physical grid, may worsen the system stability. This first-sight undesirable outcome comes as a consequence of synchronization among agents that are subject to the same signal. We argue that this effect is in fact byproduct of methodological choices, which are many times implicit. To support this claim, we employ a different approach that understands the electricity system as constituted by physical, informational and regulatory (networked and structured) layers that cannot be reduced to only one or two of them, but have to be viewed as an organic whole. By classifying its structure under this lens, more appropriate management tools can be designed by looking at the system totality in action. Two examples are provided to illustrate the strength of this modeling.


title: Learning Identifiable Gaussian Bayesian Networks in Polynomial Time and Sample Complexity

abstract: Learning the directed acyclic graph (DAG) structure of a Bayesian network from observational data is a notoriously difficult problem for which many hardness results are known. In this paper we propose a provably polynomial-time algorithm for learning sparse Gaussian Bayesian networks with equal noise variance --- a class of Bayesian networks for which the DAG structure can be uniquely identified from observational data --- under high-dimensional settings. We show that $O(k^4 \log p)$ number of samples suffices for our method to recover the true DAG structure with high probability, where $p$ is the number of variables and $k$ is the maximum Markov blanket size. We obtain our theoretical guarantees under a condition called Restricted Strong Adjacency Faithfulness, which is strictly weaker than strong faithfulness --- a condition that other methods based on conditional independence testing need for their success. The sample complexity of our method matches the information-theoretic limits in terms of the dependence on $p$. We show that our method out-performs existing state-of-the-art methods for learning Gaussian Bayesian networks in terms of recovering the true DAG structure while being comparable in speed to heuristic methods.


title: Local Photoemission Spectra and Effects of Spatial Inhomogeneity in the BCS-BEC Crossover Regime of a Trapped Ultracold Fermi Gas

abstract: We theoretically investigate single particle excitations in the BCS (Bardeen-Cooper-Schrieffer)-BEC (Bose-Einstein-Condensation) crossover regime of an ultracold Fermi gas. Including strong pairing fluctuations within a $T$-matrix approximation, as well as effects of a harmonic trap potential in the local density approximation, we calculate the local photoemission spectrum in the normal state. Recently, JILA group has measured this quantity in a $^{40}$K Fermi gas, in order to examine homogeneous single-particle properties of this system. Comparing our results with this experiment, we show that, this attempt indeed succeeds under the JILA's experimental condition. However, we also find that the current local photoemission spectroscopy still has room for improvement, in order to examine the pseudogap phenomenon predicted in the BCS-BEC crossover region. Since ultracold Fermi gases are always in a trap, our results would be useful in applying this system to various homogeneous Fermi systems, as a quantum simulator.


title: Structural studies of the bond-orientational order and hexatic-smectic transition in liquid crystals of various compositions

abstract: We report on X-ray studies of freely suspended hexatic films of three different liquid crystal compounds. By applying angular X-ray cross-correlation analysis (XCCA) to the measured diffraction patterns the parameters of the bond-orientational (BO) order in the hexatic phase were directly determined. The temperature evolution of the BO order parameters was analyzed on the basis of the multicritical scaling theory (MCST). Our results confirmed the validity of the MCST in the whole temperature range of existence of the hexatic phase for all three compounds. The temperature dependence of the BO order parameters in the vicinity of the hexatic-smectic transition was fitted by a conventional power law with a critical exponent $\beta\approx0.1$ of extremely small value. We found that the temperature dependence of higher order harmonics of the BO order scales as the powers of the first harmonic, with exponent equal to harmonic number. This indicates a nonlinear coupling of the BO order parameters of different order. It is shown that compounds of various composition, possessing different phase sequences, display the same thermodynamic behavior in the hexatic phase and in the vicinity of the smectic-hexatic phase transition.


title: The affine approach to homogeneous geodesics in homogeneous Finsler spaces

abstract: In a recent paper, it was claimed that any homogeneous Finsler space of odd dimension admits a homogeneous geodesic through any point. For the proof, the algebraic method dealing with the reductive decomposition of the Lie algebra of the isometry group was used. However, the proof contains a serious gap. In the present paper, homogeneous geodesics in Finsler homogeneous spaces are studied using the affine method, which was developed in earlier papers by the author. The mentioned statement is proved correctly and it is further proved that any homogeneous Berwald space or homogeneous reversible Finsler space admits a homogeneous geodesic through any point.


title: Everware toolkit. Supporting reproducible science and challenge-driven education

abstract: Modern science clearly demands for a higher level of reproducibility and collaboration. To make research fully reproducible one has to take care of several aspects: research protocol description, data access, environment preservation, workflow pipeline, and analysis script preservation. Version control systems like git help with the workflow and analysis scripts part. Virtualization techniques like Docker or Vagrant can help deal with environments. Jupyter notebooks are a powerful platform for conducting research in a collaborative manner. We present project Everware that seamlessly integrates git repository management systems such as Github or Gitlab, Docker and Jupyter helping with a) sharing results of real research and b) boosts education activities. With the help of Everware one can not only share the final artifacts of research but all the depth of the research process. This been shown to be extremely helpful during organization of several data analysis hackathons and machine learning schools. Using Everware participants could start from an existing solution instead of starting from scratch. They could start contributing immediately. Everware allows its users to make use of their own computational resources to run the workflows they are interested in, which leads to higher scalability of the toolkit.