title: Towards characterization of all $3\times3$ extremal quasiconvex quadratic forms

abstract: Given a $d\times d$ quasiconvex quadratic form, $d\geq 3,$ we prove that if the determinant of its acoustic tensor is an irreducible extremal polynomial that is not identically zero, then the form itself is an extremal quasiconvex quadratic form, i.e. it loses its quasiconvexity whenever a convex quadratic form is subtracted from it. In the special case $d=3,$ we slightly weaken the condition, namely we prove, that if the determinant of the acoustic tensor of the quadratic form is an extremal polynomial that is not a perfect square, then the form itself is an extremal quadratic form. In the case $d=3$ we also prove, that if the determinant of the acoustic tensor of the form is identically zero, then the form is either an extremal or polyconvex. Also, if the determinant of the acoustic tensor of the form is a perfect square, then the form is either extremal, polycovex, or is a sum of a rank-one form and an extremal, whose acoustic tensor determinant is identically zero. Here we use the notion of extremality introduced by Milton in [\ref{bib:Mil3}]


title: Auxiliary Problem Principle of augmented Lagrangian with Varying Core Functions for Large-Scale Structured Convex Problems

abstract: The auxiliary problem principle of augmented Lagrangian (APP-AL), proposed by Cohen and Zhu (1984), aims to find the solution of a constrained optimization problem through a sequence of auxiliary problems involving augmented Lagrangian. The merits of this approach are two folds. First, the core function is usually separable, which makes the subproblems at each step decomposable and particularly attractive for parallel computing. Second, the choice of the core function is quite flexible. Consequently, by carefully specifying this function, APP-AL may reduce to some standard optimization algorithms. In this paper, we pursue enhancing such flexibility by allowing the core function to be non-identical at each step of the algorithm, and name it varying auxiliary problem principle (VAPP-AL). Depending on the problem structure, the varying core functions in VAPP-AL can be adapted to design new flexible and suitable algorithm for parallel and distributed computing. The convergence and O(1/t) convergence rate of VAPP-AL for convex problem with coupling objective and constraints is proved. Moreover, if this function is specialized to be quadratic, an o(1/t) convergence rate can be established. Interestingly, the new VAPP framework can cover several variants of Jacobian type augmented Lagrangian decomposition methods as special cases. Furthermore, our technique works for the convex problem with nonseparable objective and multi-blocks coupled linear constraints, which usually can not be handled by ADMM.


title: Molecular cooling via Raman transitions and a harmonic spectral modulation of an optical frequency comb: The role of the parity of the chirp in quantum control

abstract: A method for creation of ultracold molecules by stepwise adiabatic passage from the Feshbach state to the fundamentally ground state using an optical frequency comb is presented within a semiclassical multilevel model. The sinusoidal modulation of the spectral phase of the comb is implemented that leads to a creation of a quasi-dark dressed state having an insignificant population of the excited state manifold and, thus, efficiently mitigating decoherence in the system. In contrast, the cosine modulation does not lead to the quasi-dark state formation. The results demonstrate the importance of the parity of the spectral chirp in quantum control.


title: Conflict and Computation on Wikipedia: a Finite-State Machine Analysis of Editor Interactions

abstract: What is the boundary between a vigorous argument and a breakdown of relations? What drives a group of individuals across it? Taking Wikipedia as a test case, we use a hidden Markov model to approximate the computational structure and social grammar of more than a decade of cooperation and conflict among its editors. Across a wide range of pages, we discover a bursty war/peace structure where the systems can become trapped, sometimes for months, in a computational subspace associated with significantly higher levels of conflict-tracking "revert" actions. Distinct patterns of behavior characterize the lower-conflict subspace, including tit-for-tat reversion. While a fraction of the transitions between these subspaces are associated with top-down actions taken by administrators, the effects are weak. Surprisingly, we find no statistical signal that transitions are associated with the appearance of particularly anti-social users, and only weak association with significant news events outside the system. These findings are consistent with transitions being driven by decentralized processes with no clear locus of control. Models of belief revision in the presence of a common resource for information-sharing predict the existence of two distinct phases: a disordered high-conflict phase, and a frozen phase with spontaneously-broken symmetry. The bistability we observe empirically may be a consequence of editor turn-over, which drives the system to a critical point between them.


title: The injective Leavitt complex

abstract: For a finite quiver $Q$ without sinks, we consider the corresponding finite dimensional algebra $A$ with radical square zero. We construct an explicit compact generator for the homotopy category of acyclic complexes of injective $A$-modules. We call such a generator the injective Leavitt complex of $Q$. This terminology is justified by the following result: the differential graded endomorphism algebra of the injective Leavitt complex of $Q$ is quasi-isomorphic to the Leavitt path algebra of $Q$. Here, the Leavitt path algebra is naturally Z-graded and viewed as a differential graded algebra with trivial differential.


title: Chiral surface twists and skyrmion stability in nanolayers of cubic helimagnets

abstract: Lorentz transmission electron microscopy (LTEM) investigations of modulated states in a FeGe wedge and detailed calculations demonstrate that chiral twists arising near the surfaces of noncentrosymmetric ferromagnets (Meynell et al. Phys. Rev. B, 90, 014406 (2014)) provide a stabilization mechanism for skyrmion lattices and helicoids in cubic helimagnet nanolayers. The calculated magnetic phase diagram for free standing cubic helimagnet nanolayers shows that magnetization processes in these compounds fundamentally differ from those in bulk cubic helimagnets and are characterized by the first-order transitions between modulated phases and the formation of specific multidomain states. The paper reports LTEM observations of multidomain patterns in FeGe free-standing nanolayers.


title: Maximizing Influence in Social Networks: A Two-Stage Stochastic Programming Approach That Exploits Submodularity

abstract: We consider stochastic influence maximization problems arising in social networks. In contrast to existing studies that involve greedy approximation algorithms with a 63% performance guarantee, our work focuses on solving the problem optimally. To this end, we introduce a new class of problems that we refer to as two-stage stochastic submodular optimization models. We propose a delayed constraint generation algorithm to find the optimal solution to this class of problems with a finite number of samples. The influence maximization problems of interest are special cases of this general problem class. We show that the submodularity of the influence function can be exploited to develop strong optimality cuts that are more effective than the standard optimality cuts available in the literature. Finally, we report our computational experiments with large-scale real-world datasets for two fundamental influence maximization problems, independent cascade and linear threshold, and show that our proposed algorithm outperforms the greedy algorithm.


title: Thermal Stability Of Charged Rotating Quantum Black Holes

abstract: Criteria for thermal stability of charged rotating black holes of any dimension are derived, for horizon areas that are large relative to the Planck area (in these dimensions). The derivation is based on generic assumptions of quantum geometry based on some results of loop quantum gravity, and equilibrium statistical mechanics of the Grand Canonical ensemble. There is no explicit use of classical spacetime geometry in this analysis. The only assumption is that the mass of the black hole is a function of its horizon area, charge and angular momentum. Our stability criteria are then tested in detail against specific classical black holes in spacetime dimensions 4 and 5, whose metrics provide us with explicit relations for the dependence of the mass on the charge and angular momentum of the black holes. This enables us to predict which of these black holes are expected to be thermally unstable under Hawking radiation.


title: Spontaneous modulation of superconducting phase in Kitaev ladder

abstract: We study theoretically the two Kitaev chains put in parallel, i.e., Kitaev ladder, coupled by Josephson junction. The $\pi$-junction between the Majorana bound states at the ends of the chains competes with the usual Josephson coupling along the chain, and this frustration leads to the modulation of the phase difference of the superconducting order parameter between the two chains. We show that this modulation gives the double degeneracy of the ground states, which can be manipulated by external electric and magnetic fields.


title: Optomechanical tests of a Schr\"odinger-Newton equation for gravitational quantum mechanics

abstract: We show that optomechanical systems can test the Schr\"{o}dinger-Newton equation of gravitational quantum mechanics due to Yang et al. This equation is motivated by semiclassical gravity, a widely used theory of interacting gravitational and quantum fields. From the many-body Schr\"{o}dinger-Newton equation follows an approximate equation for the center-of-mass dynamics of macroscopic objects. This predicts a distinctive double-peaked signature in the output optical quadrature power spectral density of certain optomechanical systems. Since the Schr\"{o}dinger-Newton equation lacks free parameters, these will allow its experimental confirmation or refutation.


title: Optimal occlusion uniformly partitions red blood cells fluxes within a microvascular network

abstract: In animals, gas exchange between blood and tissues occurs in narrow vessels, whose diameter is comparable to that of a red blood cell. Red blood cells must deform to squeeze through these narrow vessels, transiently blocking or occluding the vessels they pass through. Although the dynamics of vessel occlusion have been studied extensively, it remains an open question why microvessels need to be so narrow. We study occlusive dynamics within a model microvascular network: the embryonic zebrafish trunk. We show that pressure feedbacks created when red blood cells enter the finest vessels of the trunk act together to uniformly partition red blood cells through the microvasculature. Using mathematical models as well as direct observation, we show that these occlusive feedbacks are tuned throughout the trunk network to prevent the vessels closest to the heart from short-circuiting the network. Thus occlusion is linked with another open question of microvascular function: how are red blood cells delivered at the same rate to each micro-vessel? Our analysis shows that tuning of occlusive feedbacks increase the total dissipation within the network by a factor of 11, showing that uniformity of flows rather than minimization of transport costs may be prioritized by the microvascular network.


title: Ballistic Guided Electron States in Graphene

abstract: Guiding electronic waves in a manner similar to photon transmission in optical fibers is key for developing the electron-optics toolbox. Here we outline a `weak guiding' approach, in which efficient diffraction around disorder results in low-loss, high-finesse electron guiding. We describe an implementation of this scheme for guiding along a narrow-width line gate in gapless and gapped graphene. A simple model for weak guiding, which relies on the Jackiw-Rebbi midgap states, is introduced and solved. The weak-guiding modes are shown to exist for confining potential of either sign and no matter how strong or weak. Modes evolve in a cyclic manner upon varying gate potential, repeatedly sweeping the Dirac gap and becoming dispersionless (flat band) at certain magic values of gate potential. Large mode widths facilitate diffraction around disorder in the core region, enabling exceptionally large mean free paths and long-range ballistic propagation.


title: Size-dependent concentration of N0 paramagnetic centres in HPHT nanodiamonds

abstract: Size-calibrated commercial nanodiamonds synthesized by high-pressure high-temperature (HPHT) technique were studied by high-frequency W and conventional X band electron paramagnetic resonance (EPR) spectroscopy. The numbers of spins in the studied samples were estimated. The core-shell model of the HPHT nanodiamonds was proposed to explain the observed dependence of the concentration of the N0 paramagnetic centers. Two other observed paramagnetic centers are attributed to the two types of structures in the nanodiamond shell.


title: Observation of Dirac-like band dispersion in LaAgSb$_2$

abstract: We present a combined angle-resolved photoemission spectroscopy (ARPES) and first-principles calculations study of the electronic structure of LaAgSb$_2$ in the entire first Brillouin zone. We observe a Dirac-cone-like structure in the vicinity of the Fermi level formed by the crossing of two linear energy bands, as well as the nested segments of Fermi surface pocket emerging from the cone. Our ARPES results show the close relationship of the Dirac cone to the charge-density-wave ordering, providing consistent explanations for exotic behaviors in this material.


title: Hypergraph Two-Coloring in the Streaming Model

abstract: We consider space-efficient algorithms for two-coloring $n$-uniform hypergraphs $H=(V,E)$ in the streaming model, when the hyperedges arrive one at a time. It is known that any such hypergraph with at most $0.7 \sqrt{\frac{n}{\ln n}} 2^n$ hyperedges has a two-coloring [Radhakrishnan & Srinivasan, RSA, 2000], which can be found deterministically in polynomial time, if allowed full access to the input. 1. Let $s^D(v, q, n)$ be the minimum space used by a deterministic one-pass streaming algorithm that on receiving an $n$-uniform hypergraph $H$ on $v$ vertices and $q$ hyperedges produces a proper two-coloring of $H$. We show that $s^D(n^2, q, n) = \Omega(q/n)$ when $q \leq 0.7 \sqrt{\frac{n}{\ln n}} 2^n$, and $s^D(n^2, q, n) = \Omega(\sqrt{\frac{1}{n\ln n}} 2^n)$ otherwise. 2. Let $s^R(v, q,n)$ be the minimum space used by a randomized one-pass streaming algorithm that on receiving an $n$-uniform hypergraph $H$ on $v$ vertices and $q$ hyperedges with high probability produces a proper two-coloring of $H$ (or declares failure). We show that $s^R(v, \frac{1}{10}\sqrt{\frac{n}{\ln n}} 2^n, n) = O(v \log v)$ by giving an efficient randomized streaming algorithm. The above results are inspired by the study of the number $q(n)$, the minimum possible number of hyperedges in a $n$-uniform hypergraph that is not two-colorable. It is known that $q(n) = \Omega(\sqrt{\frac{n}{\ln n}})$ [Radhakrishnan-Srinivasan] and $ q(n)= O(n^2 2^n)$ [Erdos, 1963]. Our first result shows that no space-efficient deterministic streaming algorithm can match the performance of the offline algorithm of Radhakrishnan and Srinivasan; the second result shows that there is, however, a space-efficient randomized streaming algorithm for the task.


title: Intensity Mapping Cross-Correlations: Connecting the Largest Scales to Galaxy Evolution

abstract: Intensity mapping of the neutral hydrogen (HI) is a new observational tool that can be used to efficiently map the large-scale structure of the Universe over wide redshift ranges. The power spectrum of the intensity maps contains cosmological information on the matter distribution and probes galaxy evolution by tracing the HI content of galaxies at different redshifts and the scale-dependence of HI clustering. The cross-correlation of intensity maps with galaxy surveys is a robust measure of the power spectrum which diminishes systematics caused by instrumental effects and foreground removal. We examine the cross-correlation signature at redshift z=0.9 using a variant of the semi-analytical galaxy formation model SAGE (Croton et al. 2016) applied to the Millennium simulation in order to model the HI gas of galaxies as well as their optical magnitudes based on their star-formation history. We determine the clustering of the cross-correlation power for different types of galaxies determined by their colours, acting as a proxy for their star-formation activity. We find that the cross-correlation coefficient for red quiescent galaxies falls off more quickly on smaller scales k>0.2h/Mpc than for blue star-forming galaxies. Additionally, we create a mock catalogue of highly star-forming galaxies using a selection function to mimic the WiggleZ survey, and use this to predict existing and future cross-correlation measurements of the GBT and Parkes telescope. We find that the cross-power of highly star-forming galaxies shows a higher clustering on small scales than any other galaxy type and that this significantly alters the power spectrum shape on scales k>0.2h/Mpc. We show that the cross-correlation coefficient is not negligible when interpreting the cosmological cross-power spectrum. On the other hand, it contains information about the HI content of the optically selected galaxies.


title: Effects of Disorder on a 1-D Floquet Symmetry Protected Topological Phase

abstract: Periodically driven systems, also known as Floquet systems, can realize symmetry protected topological (SPT) phases that cannot be found in equilibrium. Here, we seek to understand the effects of strong disorder on such SPT phases, working with a one-dimensional Floquet system belonging to Altland-Zirnbauer class BDI. Using the transfer matrix method, we find that the disordered system hosts an array of trivial, topological, and Floquet topological phases. We then explore the phase diagram in the case of uniformly distributed bonds. Our analytic results are confirmed by a numeric computation of the bulk topological invariants. Although all states are generically localized, near phase transitions the localization length of topological edge states diverges. We derive the critical exponents governing these divergences, finding that they take on equilibrium values, despite the fact that we are considering phases unique to non-equilibrium. We also discuss how our work sets the stage for studying interacting Floquet SPT phases.


title: Weak antilocalization and electron-electron interaction in coupled multiple-channel transport in a Bi$_2$Se$_3$ thin film

abstract: Electron transport properties of a topological insulator Bi$_2$Se$_3$ thin film are studied in Hall-bar geometry. The film with a thickness of 10 nm is grown by van der Waals epitaxy on fluorophlogopite mica and Hall-bar devices are fabricated from the as-grown film directly on the mica substrate. Weak antilocalization and electron-electron interaction effects are observed and analyzed at low temperatures. The phase-coherence length extracted from the measured weak antilocalization characteristics shows a strong power-law increase with decreasing temperature and the transport in the film is shown to occur via coupled multiple (topological surface and bulk states) channels. The conductivity of the film shows a logarithmically decrease with decreasing temperature and thus the electron-electron interaction plays a dominant role in quantum corrections to the conductivity of the film at low temperatures.


title: An $r-$process macronova/kilonova in GRB 060614: evidence for the merger of a neutron star-black hole binary

abstract: After the jet break at $t\sim 1.4$ days, the optical afterglow emission of the long-short burst GRB 060614 can be divided into two components. One is the power-law decaying forward shock afterglow emission. The other is an excess of flux in several multi-band photometric observations, which emerges at $\sim$4 days after the burst, significantly earlier than that observed for a supernova associated with a long-duration GRB. At $t>13.6$ days, the F814W-band flux drops faster than $t^{-3.2}$. Moreover, the spectrum of the excess component is very soft and the luminosity is extremely low. These observed signals are incompatible with those from weak supernovae, but the ejection of $\sim 0.1~M_\odot$ of $r-$process material from a black hole-neutron star merger, as recently found in some numerical simulations, can produce it. If this interpretation is correct, it represents the first time that a multi-epoch/band lightcurve of a Li-Paczynski macronova (also known as kilonova) has been obtained and black hole-neutron star mergers are sites of significant production of $r-$process elements.


title: Gauge color codes in two dimensions

abstract: We present a family of quantum error-correcting codes that support a universal set of transversal logic gates using only local operations on a two-dimensional array of physical qubits. The construction is a subsystem version of color codes where gauge fixing through local measurements dynamically determines which gates are transversal. Although the operations are local, the underlying code is not topological in structure, which is how the construction circumvents no-go constraints imposed by the Bravyi-K\"onig and Pastawski-Yoshida theorems. We provide strong evidence that the encoding has no error threshold in the conventional sense, though it is still possible to have logical gates with error probability much lower than that of physical gates.


title: Symplectic methods based on Pad$\acute{e}$ approximation for some stochastic Hamiltonian systems

abstract: In this article, we introduce a kind of numerical schemes, based on Pad$\acute{e}$ approximation, for two stochastic Hamiltonian systems which are treated separately. For the linear stochastic Hamiltonian systems, it is shown that the applied Pad$\acute e$ approximations $P_{(k,k)}$ give numerical solutions that inherit the symplecticity and the proposed numerical schemes based on $P_{(r,s)}$ are of mean-square order $\frac{r+s}{2}$ under appropriate conditions. In case of the special stochastic Hamiltonian systems with additive noises, the numerical method using two kinds of Pad$\acute e$ approximation $P_{(\hat r,\hat s)}$ and $P_{(\check r,\check s)}$ has mean-square order $\check r+\check s+1$ when $\hat r+\hat s=\check r+\check s+2$. Moreover, the numerical solution is symplectic if $\hat r=\hat s$.


title: Brown's lemma in second-order arithmetic

abstract: We show that Brown's lemma is equivalent to Sigma02-induction over RCA0* and that the finite version of Brown's lemma is provable in RCA0 but not in RCA0*.


title: A conditionally exactly solvable generalization of the potential step

abstract: Motivated by the interest in non-relativistic quantum mechanics for determining exact solutions to the Schrodinger equation we give two potentials that are conditionally exactly solvable. The two potentials are partner potentials and we obtain that each linearly independent solution of the Schrodinger equation includes two hypergeometric functions. Furthermore we calculate their reflection and transmission amplitudes. Finally we discuss some additional properties of these potentials.


title: Collisional cooling of light ions by co-trapped heavy atoms

abstract: We experimentally demonstrate cooling of trapped ions by collisions with co-trapped, higher mass neutral atoms. It is shown that the lighter $^{39}$K$^{+}$ ions, created by ionizing $^{39}$K atoms in a magneto-optical trap (MOT), when trapped in an ion trap and subsequently allowed to cool by collisions with ultracold, heavier $^{85}$Rb atoms in a MOT, exhibit a longer trap lifetime than without the localized $^{85}$Rb MOT atoms. A similar cooling of trapped $^{85}$Rb$^{+}$ ions by ultracold $^{133}$Cs atoms in a MOT is also demonstrated in a different experimental configuration to validate this mechanism of ion cooling by localized and centered ultracold neutral atoms. Our results suggest that cooling of ions by localized cold atoms holds for any mass ratio, thereby enabling studies on a wider class of atom-ion systems irrespective of their masses.

categories: physics.atom-ph cond-mat.quant-gas physics.chem-ph physics.optics


title: Coronal Magnetic Fields Derived from Simultaneous Microwave and EUV Observations and Comparison with the Potential Field Model

abstract: We estimated the accuracy of coronal magnetic fields derived from radio observations by comparing them to potential field calculations and the DEM measurements using EUV observations. We derived line of sight component of the coronal magnetic field from polarization observations of the thermal bremsstrahlung in the NOAA active region 11150, observed around 3:00 UT on February 3, 2011 using the Nobeyama Radioheliograph at 17 GHz. Because the thermal bremsstrahlung intensity at 17 GHz includes both chromospheric and coronal components, we extracted only the coronal component by measuring the coronal emission measure in EUV observations. In addition, we derived only the radio polarization component of the corona by selecting the region of coronal loops and weak magnetic field strength in the chromosphere along the line of sight. The upper limit of the coronal longitudinal magnetic fields were determined as 100 - 210 G. We also calculated the coronal longitudinal magnetic fields from the potential field extrapolation using the photospheric magnetic field obtained from the Helioseismic and Magnetic Imager (HMI). However, the calculated potential fields were certainly smaller than the observed coronal longitudinal magnetic field. This discrepancy between the potential and the observed magnetic field strengths can be explained consistently by two reasons; (a) the underestimation of the coronal emission measure resulting from the limitation of the temperature range of the EUV observations, (b) the underestimation of the coronal magnetic field resulting from the potential field assumption.


title: Towards Lattice QCD Baryon Forces at the Physical Point: First Results

abstract: Lattice QCD calculations of baryon forces are performed for the first time with (almost) physical quark masses. $N_f = 2+1$ dynamical clover fermion gauge configurations are generated at the lattice spacing of $a \simeq 0.085$ fm on a $(96 a)^4 \simeq (8.2 {\rm fm})^4$ lattice with quark masses corresponding to $(m_\pi, m_K) \simeq (146, 525)$ MeV. Baryon forces are calculated using the time-dependent HAL QCD method. In this report, we study $\Xi\Xi$ and $NN$ systems both in $^1S_0$ and $^3S_1$-$^3D_1$ channels, and the results for the central and tensor forces as well as phase shifts in the $\Xi\Xi$ $(^1S_0)$ channel are presented.


title: Parameterized Algorithms on Perfect Graphs for deletion to $(r,\ell)$-graphs

abstract: For fixed integers $r,\ell \geq 0$, a graph $G$ is called an {\em $(r,\ell)$-graph} if the vertex set $V(G)$ can be partitioned into $r$ independent sets and $\ell$ cliques. The class of $(r, \ell)$ graphs generalizes $r$-colourable graphs (when $\ell =0)$ and hence not surprisingly, determining whether a given graph is an $(r, \ell)$-graph is \NP-hard even when $r \geq 3$ or $\ell \geq 3$ in general graphs. When $r$ and $\ell$ are part of the input, then the recognition problem is NP-hard even if the input graph is a perfect graph (where the {\sc Chromatic Number} problem is solvable in polynomial time). It is also known to be fixed-parameter tractable (FPT) on perfect graphs when parameterized by $r$ and $\ell$. I.e. there is an $f(r+\ell) \cdot n^{\Oh(1)}$ algorithm on perfect graphs on $n$ vertices where $f$ is some (exponential) function of $r$ and $\ell$. In this paper, we consider the parameterized complexity of the following problem, which we call {\sc Vertex Partization}. Given a perfect graph $G$ and positive integers $r,\ell,k$ decide whether there exists a set $S\subseteq V(G)$ of size at most $k$ such that the deletion of $S$ from $G$ results in an $(r,\ell)$-graph. We obtain the following results: \begin{enumerate} \item {\sc Vertex Partization} on perfect graphs is FPT when parameterized by $k+r+\ell$. \item The problem does not admit any polynomial sized kernel when parameterized by $k+r+\ell$. In other words, in polynomial time, the input graph can not be compressed to an equivalent instance of size polynomial in $k+r+\ell$. In fact, our result holds even when $k=0$. \item When $r,\ell$ are universal constants, then {\sc Vertex Partization} on perfect graphs, parameterized by $k$, has a polynomial sized kernel. \end{enumerate}


title: Ab initio derivation of multipolar expansion of optical force

abstract: Like many other physical quantities, the optical force can be expanded using multipole expansion, which has been done in [Nat. Photon. 5, 531], up to electric octupole order. However, in that study, the existence of radiation multipoles were pre-assumed, and the role of the fundamental building units, charges, are not evident. Here, we derive the same multipolar expression of optical force by treating the particles as a collection of point charges or point dipoles, which results in more transparent physics and mathematics.


title: Preconditioned Stochastic Gradient Descent

abstract: Stochastic gradient descent (SGD) still is the workhorse for many practical problems. However, it converges slow, and can be difficult to tune. It is possible to precondition SGD to accelerate its convergence remarkably. But many attempts in this direction either aim at solving specialized problems, or result in significantly more complicated methods than SGD. This paper proposes a new method to estimate a preconditioner such that the amplitudes of perturbations of preconditioned stochastic gradient match that of the perturbations of parameters to be optimized in a way comparable to Newton method for deterministic optimization. Unlike the preconditioners based on secant equation fitting as done in deterministic quasi-Newton methods, which assume positive definite Hessian and approximate its inverse, the new preconditioner works equally well for both convex and non-convex optimizations with exact or noisy gradients. When stochastic gradient is used, it can naturally damp the gradient noise to stabilize SGD. Efficient preconditioner estimation methods are developed, and with reasonable simplifications, they are applicable to large scaled problems. Experimental results demonstrate that equipped with the new preconditioner, without any tuning effort, preconditioned SGD can efficiently solve many challenging problems like the training of a deep neural network or a recurrent neural network requiring extremely long term memories.


title: A rephasing invariant study of neutrino mixing

abstract: We derive a set of renormalization group equations (RGE) for Dirac neutrinos using a rephasing invariant parametrization. The symmetric properties of these equations under flavor permutation facilitate the derivation of some exact and approximate RGE invariants. Even though the complete analytical solutions for the RGE are unavailable, we provide a numerical example that illustrate the evolution of the neutrino mixing parameters.


title: On the Cohen-Macaulayness of tangent cones of monomial curves in $\mathbb{A}^{4}(K)$

abstract: In this paper we give necessary and sufficient conditions for the Cohen-Macaulayness of the tangent cone of a monomial curve in the 4-dimensional affine space. We study particularly the case where $C$ is a Gorenstein non-complete intersection monomial curve.


title: Compressed Dynamic Mode Decomposition for Background Modeling

abstract: We introduce the method of compressed dynamic mode decomposition (cDMD) for background modeling. The dynamic mode decomposition (DMD) is a regression technique that integrates two of the leading data analysis methods in use today: Fourier transforms and singular value decomposition. Borrowing ideas from compressed sensing and matrix sketching, cDMD eases the computational workload of high resolution video processing. The key principal of cDMD is to obtain the decomposition on a (small) compressed matrix representation of the video feed. Hence, the cDMD algorithm scales with the intrinsic rank of the matrix, rather then the size of the actual video (data) matrix. Selection of the optimal modes characterizing the background is formulated as a sparsity-constrained sparse coding problem. Our results show, that the quality of the resulting background model is competitive, quantified by the F-measure, Recall and Precision. A GPU (graphics processing unit) accelerated implementation is also presented which further boosts the computational efficiency of the algorithm.


title: Electronic structures of [111]-oriented free-standing InAs and InP nanowires

abstract: We report on a theoretical study of the electronic structures of the [111]-oriented, free-standing, zincblende InAs and InP nanowires with hexagonal cross sections by means of an atomistic $sp^{3}s^{*} $, spin-orbit interaction included, nearest-neighbor, tight-binding method. The band structures and the band state wave functions of these nanowires are calculated and the symmetry properties of the bands and band states are analyzed based on the $C_{3v}$ double point group. It is shown that all bands of these nanowires are doubly degenerate at the $\Gamma$-point and some of these bands will split into non-degenerate bands when the wave vector $k$ moves away from the $\Gamma$-point as a manifestation of spin-splitting due to spin-orbit interaction. It is also shown that the lower conduction bands of these nanowires all show simple parabolic dispersion relations, while the top valence bands show complex dispersion relations and band crossings. The band state wave functions are presented by the spatial probability distributions and it is found that all the band states show $2\pi/3$-rotation symmetric probability distributions. The effects of quantum confinement on the band structures of the [111]-oriented InAs and InP nanowires are also examined and an empirical formula for the description of quantization energies of the lowest conduction band and the highest valence band is presented. The formula can simply be used to estimate the enhancement of the band gaps of the nanowires at different sizes as a result of quantum confinement.


title: A review of mu-tau flavor symmetry in neutrino physics

abstract: Behind the observed pattern of lepton flavor mixing is a partial or approximate mu-tau flavor symmetry --- a milestone on our road to the true origin of neutrino masses and flavor structures. In this review article we first describe the features of mu-tau permutation and reflection symmetries, and then explore their various consequences on model building and neutrino phenomenology. We pay particular attention to soft mu-tau symmetry breaking, which is crucial for our deeper understanding of the fine effects of flavor mixing and CP violation.


title: Watch-Bot: Unsupervised Learning for Reminding Humans of Forgotten Actions

abstract: We present a robotic system that watches a human using a Kinect v2 RGB-D sensor, detects what he forgot to do while performing an activity, and if necessary reminds the person using a laser pointer to point out the related object. Our simple setup can be easily deployed on any assistive robot. Our approach is based on a learning algorithm trained in a purely unsupervised setting, which does not require any human annotations. This makes our approach scalable and applicable to variant scenarios. Our model learns the action/object co-occurrence and action temporal relations in the activity, and uses the learned rich relationships to infer the forgotten action and the related object. We show that our approach not only improves the unsupervised action segmentation and action cluster assignment performance, but also effectively detects the forgotten actions on a challenging human activity RGB-D video dataset. In robotic experiments, we show that our robot is able to remind people of forgotten actions successfully.


title: Higher Groupoid Actions, Bibundles, and Differentiation

abstract: In this thesis, we employ simplicial methods to study actions, principal bundles, and bibundles of higher groupoids. Roughly, we use Kan fibrations to model actions of higher groupoids, we use pairs of a Kan fibration and a special acyclic fibration to model principal bundles of higher groupoids, we use inner Kan fibrations over the interval to model bibundles of higher groupoids. In particular, we show that our definitions given by the simplicial method agree with those given by the categorification approach to actions, principal bundles, and bibundles of 2-groupoids. In addition, we use the simplicial technique to prove a theorem on differentiation of higher Lie groupoids, which shows that the differentiation functor sends a higher Lie groupoid to a higher Lie algebroid.


title: Refinement rings, localization and diagonal reduction of matrices

abstract: In this paper we prove that if $R$ is a commutative refinement ring and $M$, $N$ are two $R$-modules then, $M\cong N$ if and only if for every maximal ideal $m$ of $R$, $M_m\cong N_m$. We prove if $R$ is a refinement ring, then every regular matrix over $\frac{R}{J(R)}$ admits a diagonal reduction iff every regular matrix over $R$ admits a diagonal reduction.


title: Anderson Localization in Degenerate Spin-Orbit Coupled Fermi Gases with Disorder

abstract: Competition between superconductivity and disorder plays an essential role in understanding the metal-insulator transition. Based on the Bogoliubov-de Gennes equation, we studied an s-wave superconductor with both spin-orbit coupling and disorder are presented. With increasing the strength of disorder, the mean superconducting order parameter will vanish while the energy gap will persist, which indicates that the system undergoes a transition from a superconducting state to an insulating state. This can be confirmed by calculating the inverse participation ratio. We also find that, if the strength of disorder is small, the superconducting order parameter and the energy gap will decrease if we increase the strength of spin-orbit coupling and Zeeman field. In the large disorder limits, increasing the strength of spin-orbit coupling will increase the mean superconducting order parameter. This phenomenon shows that the system is more insensitive to disorder if the spin-orbit coupling is presented. Numerical computing also shows that the whole system breaks up into several superconducting islands instead of being superconductivity.


title: Ingham type inequalities in lattices

abstract: A classical theorem of Ingham extended Parseval's formula of the trigonometrical system to arbitrary families of exponentials satisfying a uniform gap condition. Later his result was extended to several dimensions, but the optimal integration domains have only been determined in very few cases. The purpose of this paper is to determine the optimal connected integration domains for all regular two-dimensional lattices.


title: Quenching effect in below-threshold high harmonic generation

abstract: We theoretically demonstrate the quenching effect in below-threshold high harmonic generation (HHG) by using the time-dependent density-functional theory (TDDFT) and solving the time-dependent Schr\"{o}dinger equation (TDSE). It is shown that the HHG is substantially suppressed in particular harmonic orders in the below-threshold region when multi-electron interaction comes into play. The position of the suppression is determined by the energy gap between the highest occupied orbital and the higher-lying orbital of the target. We show that the quenching effect is due to a new class of multi-electron dynamics involving electron-electron energy transfer, which is analog to the fluorescence quenching owing to the energy transfer between molecules in fluorescent material. This work reveals the important role of the multi-electron interaction on HHG especially in the below-threshold region.


title: The Globalization of Academic Entrepreneurship? The Recent Growth (2009-2014) in University Patenting Decomposed

abstract: The contribution of academia to US patents has become increasingly global. Following a pause, with a relatively flat rate, from 1998 to 2008, the long-term trend of university patenting rising as a share of all patenting has resumed, driven by the internationalization of academic entrepreneurship and the persistence of US university technology transfer. We disaggregate this recent growth in university patenting at the US Patent and Trademark Organization (USPTO) in terms of nations and patent classes. Foreign patenting in the US has almost doubled during the period 2009-2014, mainly due to patenting by universities in Taiwan, Korea, China, and Japan. These nations compete with the US in terms of patent portfolios, whereas most European countries--with the exception of the UK--have more specific portfolios, mainly in the bio-medical fields. In the case of China, Tsinghua University holds 63% of the university patents in USPTO, followed by King Fahd University with 55.2% of the national portfolio.


title: Entire slice regular functions

abstract: Entire functions in one complex variable are extremely relevant in several areas ranging from the study of convolution equations to special functions. An analog of entire functions in the quaternionic setting can be defined in the slice regular setting, a framework which includes polynomials and power series of the quaternionic variable. In the first chapters of this work we introduce and discuss the algebra and the analysis of slice regular functions. In addition to offering a self-contained introduction to the theory of slice-regular functions, these chapters also contain a few new results (for example we complete the discussion on lower bounds for slice regular functions initiated with the Ehrenpreis-Malgrange, by adding a brand new Cartan-type theorem). The core of the work is Chapter 5, where we study the growth of entire slice regular functions, and we show how such growth is related to the coefficients of the power series expansions that these functions have. It should be noted that the proofs we offer are not simple reconstructions of the holomorphic case. Indeed, the non-commutative setting creates a series of non-trivial problems. Also the counting of the zeros is not trivial because of the presence of spherical zeros which have infinite cardinality. We prove the analog of Jensen and Carath\'eodory theorems in this setting.


title: Strain and Vector-Magnetic-Field Tuning of the Anomalous Phase in Sr3Ru2O7

abstract: Symmetry is a basic paradigm for analysis of problems in solid state physics, and to controllably break point-group symmetries is a powerful probe of novel forms of order. In this work, we apply two methods, controlled lattice strain and vector magnetic field, to study a well-known correlated electron phase that appears at low temperatures in the layered material Sr3Ru2O7. This phase responds strongly to C4-symmetry-breaking fields, suggesting spontaneous C4 symmetry breaking. However, in this work we find that the data are more easily explained with microscopic coexistence of two orthogonal order parameter components than spontaneous symmetry breaking. The techniques presented in this work will be relevant for other systems where spontaneous C4 symmetry breaking is suspected.


title: Physical Limits on Bacterial Navigation in Dynamic Environments

abstract: Many chemotactic bacteria inhabit environments in which chemicals appear as localized pulses and evolve by processes such as diffusion and mixing. We show that, in such environments, physical limits on the accuracy of temporal gradient sensing govern when and where bacteria can accurately measure the cues they use to navigate. Chemical pulses are surrounded by a predictable dynamic region, outside which bacterial cells cannot resolve gradients above noise. The outer boundary of this region initially expands in proportion to $\sqrt{t}$, before rapidly contracting. Our analysis also reveals how chemokinesis - the increase in swimming speed many bacteria exhibit when absolute chemical concentration exceeds a threshold - may serve to enhance chemotactic accuracy and sensitivity when the chemical landscape is dynamic. More generally, our framework provides a rigorous method for partitioning bacteria into populations that are "near" and "far" from chemical hotspots in complex, rapidly evolving environments such as those that dominate aquatic ecosystems.


title: Crossings states and sets of states in P\'olya random walks

abstract: We consider the P\'olya random walk in $\mathbb{Z}^2$. The paper establishes a number of results for the distributions and expectations of the number of usual (undirected) and specifically defined in the paper up- and down-directed state-crossings and different sets of states crossings. One of the most important results of this paper is that the expected number of undirected state-crossings $\mathbf{n}$ is equal to 1 for any state $\mathbf{n}\in\mathbb{Z}^2\setminus\{\mathbf{0}\}$. As well, the results of the paper are extended to $d$-dimensional random walks, $d\geq2$, in bounded areas.


title: On the Relation between two Rotation Metrics

abstract: In their work "Global Optimization through Rotation Space Search", Richard Hartley and Fredrik Kahl introduce a global optimization strategy for problems in geometric computer vision, based on rotation space search using a branch-and-bound algorithm. In its core, Lemma 2 of their publication is the important foundation for a class of global optimization algorithms, which is adopted over a wide range of problems in subsequent publications. This lemma relates a metric on rotations represented by rotation matrices with a metric on rotations in axis-angle representation. This work focuses on a proof for this relationship, which is based on Rodrigues' Rotation Theorem for the composition of rotations in axis-angle representation.


title: Jeans Analysis of Bok globules in $f(R)$ gravity

abstract: We examine the effects of $f(R)$ gravity on Jeans analysis of collapsing dust clouds. We provide a method for testing modified gravity models by their effects on star formation as the presence of $f(R)$ gravity is found to modify the limit for collapse. In this analysis we add perturbations to a de Sitter background. Depending on the characteristics of a chosen $f(R)$ model, the appearance of new limits is possible. The physicality of these limits is further examined. We find the asymptotic Jeans masses for $f(R)$ theories compared to standard Jeans mass. Through this ratio, the effects of the $f(R)$ modified Jeans mass for viable theories are examined in molecular clouds. Bok globules have a mass range comparable to Jeans masses in question and are therefore used for comparing different $f(R)$ models. Viable theories are found to assist in star formation.


title: A non-boost-invariant solution of relativistic hydrodynamics in 1+3 dimensions

abstract: We present a new solution of relativistic hydrodynamics in 1+3 dimensions which depends on both the transverse coordinate and rapidity. At early times the flow expands dominantly longitudinally in a non-boost-invariant manner, and at late times it expands nearly spherically. These two regimes are shown to be related by symmetry. The effect of viscosity is also discussed.


title: Amortized Averaging Algorithms for Approximate Consensus

abstract: We introduce a new class of distributed algorithms for the approximate consensus problem in dynamic rooted networks, which we call amortized averaging algorithms. They are deduced from ordinary averaging algorithms by adding a value-gathering phase before each value update. This allows their decision time to drop from being exponential in the number $n$ of processes to being linear under the assumption that each process knows $n$. In particular, the amortized midpoint algorithm, which achieves a linear decision time, works in completely anonymous dynamic rooted networks where processes can exchange and store continuous values, and under the assumption that the number of processes is known to all processes. We then study the way amortized averaging algorithms degrade when communication graphs are from time to time non rooted, or with a wrong estimate of the number of processes. Finally, we analyze the amortized midpoint algorithm under the additional constraint that processes can only store and send quantized values, and get as a corollary that the 2-set consensus problem is solvable in linear time in any rooted dynamic network model when allowing all decision values to be in the range of initial values.


title: Shear mixing in stellar radiative zones - II. Robustness of numerical simulations

abstract: Recent numerical simulations suggest that the model by Zahn (1992, A&A, 265, 115) for the turbulent mixing of chemical elements due to differential rotation in stellar radiative zones is valid. We investigate the robustness of this result with respect to the numerical configuration and Reynolds number of the flow. We compare results from simulations performed with two different numerical codes, including one that uses the shearing-box formalism. We also extensively study the dependence of the turbulent diffusion coefficient on the turbulent Reynolds number. The two numerical codes used in this study give consistent results. The turbulent diffusion coefficient is independent of the size of the numerical domain if at least three large turbulent structures fit in the box. Generally, the turbulent diffusion coefficient depends on the turbulent Reynolds number. However, our simulations suggest that an asymptotic regime is obtained when the turbulent Reynolds number is larger than $10^3$. Shear mixing in the regime of small P\'eclet numbers can be investigated numerically both with shearing-box simulations and simulations using explicit forcing. Our results suggest that Zahn's model is valid at large turbulent Reynolds numbers.


title: Vacuum Selection on Axionic Landscapes

abstract: We compute the distribution of minima that are reached dynamically on multi-field axionic landscapes, both numerically and analytically. Such landscapes are well suited for inflationary model building due to the presence of shift symmetries and possible alignment effects (the KNP mechanism). The resulting distribution of dynamically reached minima differs considerably from the naive expectation based on counting all vacua. These differences are more pronounced in the presence of many fields due to dynamical selection effects: while low lying minima are preferred as fields roll down the potential, trajectories are also more likely to get trapped by one of the many nearby minima. We show that common analytic arguments based on random matrix theory in the large $D$-limit to estimate the distribution of minima are insufficient for quantitative arguments pertaining to the dynamically reached ones. This discrepancy is not restricted to axionic potentials. We provide an empirical expression for the expectation value of such dynamically reached minimas' height and argue that the cosmological constant problem is not alleviated in the absence of anthropic arguments. We further comment on the likelihood of inflation on axionic landscapes in the large D-limit.


title: Conversions of Bound Muons: Lepton Flavour Violation from Doubly Charged Scalars

abstract: We present the first detailed computation of the conversion of a bound muon into an electron mediated by a doubly charged $SU(2)$ singlet scalar. Although such particles are not too exotic, up to now their contribution to $\mu$-$e$ conversion is unknown. We close this gap by presenting a detailed calculation, which will allow the reader not only to fully comprehend the discussion but also to generalise our results to similar cases if needed. We furthermore compare the predictions, for both the general case and for an example model featuring a neutrino mass at 2-loop level, to current experimental data and future sensitivities. We show that, depending on the explicit values of the couplings as well as on the actual future limits on the branching ratio, $\mu$-$e$ conversion may potentially yield a lower limit on the doubly charged singlet scalar mass which is stronger than what could be obtained by colliders. Our results considerably strengthen the case for low-energy lepton flavour violation searches being a very valuable addition to collider experiments.


title: Random Sampling with Removal

abstract: We study randomized algorithms for constrained optimization, in abstract frameworks that include, in strictly increasing generality: convex programming; LP-type problems; violator spaces; and a setting we introduce, consistent spaces. Such algorithms typically involve a step of finding the optimal solution for a random sample of the constraints. They exploit the condition that, in finite dimension $\delta$, this sample optimum violates a provably small expected fraction of the non-sampled constraints, with the fraction decreasing in the sample size~$r$. We extend such algorithms by considering the technique of removal, where a fixed number $k$ of constraints are removed from the sample according to a fixed rule, with the goal of improving the solution quality. This may have the effect of increasing the number of violated non-sampled constraints. We study this increase, and bound it in a variety of general settings. This work is motivated by, and extends, results on removal as proposed for chance-constrained optimization. For many relevant values of $r$, $\delta$, and $k$, we prove matching upper and lower bounds for the expected number of constraints violated by a random sample, after the removal of $k$ constraints. For a large range of values of $k$, the new upper bounds improve the previously best bounds for LP-type problems, which moreover had only been known in special cases, and not in the generality we consider. Moreover, we show that our results extend from finite to infinite spaces, for chance-constrained optimization.


title: Velocity distribution in active particles systems

abstract: We derive an analytic expression for the distribution of velocities of multiple interacting active particles which we test by numerical simulations. In clear contrast with equilibrium we find that the velocities are coupled to positions. Our model shows that, even for two particles only, the individual velocities display a variance depending on the interparticle separation and the emergence of correlations between the velocities of the particles. When considering systems composed of many particles we find an analytic expression connecting the overall velocity variance to density, at the mean-field level, and to the pair distribution function valid in the limit of small noise correlation times. Finally we discuss the intriguing analogies and main differences between our effective free energy functional and the theoretical scenario proposed so far for phase-separating active particles.


title: An Asymptotic-Preserving Method for a Relaxation of the Navier-Stokes-Korteweg Equations

abstract: The Navier-Stokes-Korteweg (NSK) equations are a classical diffuse-interface model for compressible two-phase flow. As direct numerical simulations based on the NSK system are quite expensive and in some cases even impossible, we consider a relaxation of the NSK system, for which robust numerical methods can be designed. However, time steps for explicit numerical schemes depend on the relaxation parameter and therefore numerical simulations in the relaxation limit are very inefficient. To overcome this restriction, we propose an implicit-explicit asymptotic-preserving finite volume method. We prove that the new scheme provides a consistent discretization of the NSK system in the relaxation limit and demonstrate that it is capable of accurately and efficiently computing numerical solutions of problems with realistic density ratios and small interfacial widths.


title: Quantum oscillations and the Fermi-surface topology of the Weyl semimetal NbP

abstract: The Weyl semimetal NbP was found to exhibit topological Fermi arcs and exotic magneto-transport properties. Here, we report on magnetic quantum-oscillation measurements on NbP and construct the 3D Fermi surface with the help of band-structure calculations. We reveal a pair of spin-orbit-split electron pockets at the Fermi energy and a similar pair of hole pockets, all of which are strongly anisotropic. The Fermi surface well explains the linear magnetoresistance observed in high magnetic fields by the quantum-limit scenario. The Weyl points that are located in the $k_z \approx \pi/c$ plane are found to exist 5 meV above the Fermi energy. Therefore, we predict that the chiral anomaly effect can be realized in NbP by electron doping to drive the Fermi energy to the Weyl points.


title: Direct photon measurement in Pb-Pb collisions at $\sqrt{s_{\rm NN}}$ = 2.76 TeV with ALICE

abstract: The ALICE experiment has measured the direct photon spectra in Pb-Pb collisions at $\sqrt{s_{NN}} = 2.76$ TeV for three different centrality selections. The measurement was performed emplying a method utilizing conversion of photons into $e^+e^-$ pairs in the detector material, and a method using the PHOS calorimeter. The two measurements were combined in order to measure direct photons over a broad transverse momentum range of $0.9 < p_{\rm T} < 14$ GeV/$c$. A direct photon signal was observed in the thermal region $0.9 < p_{\rm T} < 2.1$ GeV/$c$, with a significance of $2.6\sigma$ in the 20\% most central collisions. The corresponding direct photon spectrum can be described by an exponential with an inverse slope parameter of $(304 \pm 11^{\rm stat} \pm 40^{\rm syst})$ MeV, without subtracting the hard scattering component, which is somewhat larger than the inverse slope parameter of the direct photon spectrum at RHIC energies. Within uncertainties, the data can be described by different models for photon production in heavy-ion collisions. These proceedings provide a summary of the results published in [1] and [2].


title: Semiclassical treatment for molecular rotation spectra in high electric fields

abstract: Molecular rotation spectra, generated by the coupling of the molecular electric-dipole moments to an external time-dependent electric field, are discussed in a few particular conditions which can be of some experimental interest. First, the spherical-pendulum molecular model is reviewed, with the aim of introducing an approximate method which consists in the separation of the azimuthal and zenithal motions. Second, rotation spectra are considered in the presence of a static electric field. Two particular cases are analyzed, corresponding to strong and weak fields. In both cases the dipoles perform rotations and vibrations about equilibrium positions, which may exhibit parametric resonances. For strong fields a large macroscopic electric polarization may appear. This situation may be relevant for polar matter (like pyroelectrics, ferroelectrics), or for heavy impurities embedded in a polar solid. The dipolar interaction is analyzed in polar condensed matter, where it is shown that new polarization modes appear for a spontaneous macroscopic electric polarization (these modes are tentatively called "dipolons"); one of the polarization modes is related to parametric resonances. The extension of these considerations to magnetic dipoles is briefly discussed. The treatment is extended to strong electric fields which oscillate with a high frequency, as those provided by high-power lasers. It is shown that the effect of such fields on molecular dynamics is governed by a much weaker, effective, renormalized, static electric field.


title: Collective and Single-particle Motion in Beyond Mean Field Approaches

abstract: We present a novel nuclear energy density functional method to calculate spectroscopic properties of atomic nuclei. Intrinsic nuclear quadrupole deformations and rotational frequencies are considered simultaneously as the degrees of freedom within a symmetry conserving configuration mixing framework. The present method allows the study of nuclear states with collective and single-particle character. We calculate the fascinating structure of the semi-magic 44S nucleus as a first application of the method, obtaining an excellent quantitative agreement both with the available experimental data and with state-of-the-art shell model calculations.


title: Apsidal motion and a light curve solution for thirteen LMC eccentric eclipsing binaries

abstract: New CCD observations for thirteen eccentric eclipsing binaries from the Large Magellanic Cloud were carried out using the Danish 1.54-meter telescope located at the La Silla Observatory in Chile. These systems were observed for their times of minima and 56 new minima were obtained. These are needed for accurate determination of the apsidal motion. Besides that, in total 436 times of minima were derived from the photometric databases OGLE and MACHO. The (O-C) diagrams of minima timings for these B-type binaries were analysed and the parameters of the apsidal motion were computed. The light curves of these systems were fitted using the program PHOEBE, giving the light curve parameters. We derived for the first time the relatively short periods of the apsidal motion ranging from 21 to 107 years. The system OGLE-LMC-ECL-07902 was also analysed using the spectra and radial velocities, resulting in masses 6.8, and 4.4 M0 for the eclipsing components. For one system (OGLE-LMC-ECL-20112), the third-body hypothesis was also used for describing the residuals after subtraction of the apsidal motion, resulting in period of about 22 years. For several systems an additional third light was also detected, which makes these systems suspicious of triplicity. Based on data collected with the Danish 1.54-m telescope at the ESO La Silla Observatory.


title: Two applications of the spectrum of numbers

abstract: Let the base $\beta$ be a complex number, $|\beta|>1$, and let $A \subset \C$ be a finite alphabet of digits. The \emph{$A$-spectrum} of $\beta$ is the set $S_{A}(\beta) = \{\sum_{k=0}^n a_k\beta^k \mid n \in \mathbb{N}, \ a_k \in {A}\}$. We show that the spectrum $S_{{A}}(\beta)$ has an accumulation point if and only if $0$ has a particular $(\beta, A)$-representation, said to be \emph{rigid}. The first application is restricted to the case that $\beta >1 $ and the alphabet is $A=\{-M, \ldots, M\}$, $M \ge 1$ integer. We show that the set $Z_{\beta,M}$ of infinite $(\beta, A)$-representations of $0$ is recognizable by a finite B\"uchi automaton if and only if the spectrum $S_A(\beta)$ has no accumulation point. Using a result of Akiyama-Komornik and Feng, this implies that $Z_{\beta, M}$ is recognizable by a finite B\"uchi automaton for any positive integer $M \ge \lceil \beta \rceil -1$ if and only if $\beta$ is a Pisot number. This improves the previous bound $M \ge \lceil \beta \rceil $. For the second application the base and the digits are complex. We consider the on-line algorithm for division of Trivedi and Ercegovac generalized to a complex numeration system. In on-line arithmetic the operands and results are processed in a digit serial manner, starting with the most significant digit. The divisor must be far from $0$, which means that no prefix of the $(\beta,A)$-representation of the divisor can be small. The numeration system $(\beta,A)$ is said to \emph{allow preprocessing} if there exists a finite list of transformations on the divisor which achieve this task. We show that $(\beta,A )$ allows preprocessing if and only if the spectrum $S_{{A}}(\beta)$ has no accumulation point.


title: Network theory for inhomogeneous thermoelectrics

abstract: The Onsager-de Groot-Callen transport theory, implemented as a network model, is used to simulate the transient Harman method, which is widely used experimentally to determine all thermoelectric transport coefficients in a single measurement setup. It is shown that this method systematically overestimates the Seebeck coefficient for samples composed of two different materials. As a consequence, the figure of merit is also overestimated, if the thermal coupling of the measurement setup to the environment is weak. For a mixture of metal and semiconductor particles near metal percolation the figure of merit obtained by the Harman method is more than 100 \% too large. For a correct interpretation of the experimental data, information on composition and microstructure of the sample are indispensable.


title: Tessellating the moduli space of strictly convex projective structures on the once-punctured torus

abstract: We show that associating the Euclidean cell decomposition due to Cooper and Long to each point of the moduli space of framed strictly convex real projective structures of finite volume on the once-punctured torus gives this moduli space a natural cell decomposition. The proof makes use of coordinates due to Fock and Goncharov, the action of the mapping class group as well as algorithmic real algebraic geometry. We also show that the decorated moduli space of framed strictly convex real projective structures of finite volume on the thrice-punctured sphere has a natural cell decomposition.


title: Growth and cogrowth of normal subgroups of a free group

abstract: We give a sufficient condition for a sequence of normal subgroups of a free group to have the property that both, their growths tend to the upper bound and their cogrowths tend to the lower bound. The condition is represented by planarity of the quotient graphs of the tree.


title: Finite-temperature effects on interacting bosonic 1D systems in disordered lattices

abstract: We analyze the finite-temperature effects on the phase diagram describing the insulating properties of interacting 1D bosons in a quasi-periodic lattice. We examine thermal effects by comparing experimental results to exact diagonalization for small-sized systems and to density-matrix renormalization group (DMRG) computations. At weak interactions, we find short thermal correlation lengths, indicating a substantial impact of temperature on the system coherence. Conversely, at strong interactions, the obtained thermal correlation lengths are significantly larger than the localization length, and the quantum nature of the T=0 Bose glass phase is preserved up to a crossover temperature that depends on the disorder strength. Furthermore, in the absence of disorder, we show how quasi-exact finite-T DMRG computations, compared to experimental results, can be employed to estimate the temperature, which is not directly accessible in the experiment.


title: Status of the ANAIS Dark Matter Project at the Canfranc Underground Laboratory

abstract: The ANAIS (Annual modulation with NaI(Tl) Scintillators) experiment aims at the confirmation of the DAMA/LIBRA signal using the same target and technique at the Canfranc Underground Laboratory (LSC). Along 2016, 112.5 kg of ultra pure NaI(Tl) crystals will be installed at LSC in a 3x3 modules matrix configuration. The ANAIS-25 and ANAIS-37 set-ups have been taking data at the LSC testing the detector performance, the DAQ and analysis systems, and assessing the background. Main results coming from both set-ups will be summarized in this paper, focusing on the excellent detector performance and background understanding. Prospects for the experiment will be also briefly revised.


title: Perturbed Newtonian description of the Lema\^itre model with non-negligible pressure

abstract: We study the validity of the Newtonian description of cosmological perturbations using the Lemaitre model, an exact spherically symmetric solution of Einstein's equation. This problem has been investigated in the past for the case of a dust fluid. Here, we extend the previous analysis to the more general case of a fluid with non-negligible pressure, and, for the numerical examples, we consider the case of radiation (P=\rho/3). We find that, even when the density contrast has a nonlinear amplitude, the Newtonian description of the cosmological perturbations using the gravitational potential \psi and the curvature potential \phi is valid as long as we consider sub-horizon inhomogeneities. However, the relation \psi+\phi={\cal O}(\phi^2), which holds for the case of a dust fluid, is not valid for a relativistic fluid and effective anisotropic stress is generated. This demonstrates the usefulness of the Lemaitre model which allows us to study in an exact nonlinear fashion the onset of anisotropic stress in fluids with non-negligible pressure. We show that this happens when the characteristic scale of the inhomogeneity is smaller than the sound horizon and that the deviation is caused by the nonlinear effect of the fluid's fast motion. We also find that \psi+\phi= \max[{\cal O}(\phi^2),{\cal O}(c_s^2\phi \, \delta)] for an inhomogeneity with density contrast \delta whose characteristic scale is smaller than the sound horizon, unless w is close to -1, where w and c_s are the equation of state parameter and the sound speed of the fluid, respectively. On the other hand, we expect \psi+\phi={\cal O}(\phi^2) to hold for an inhomogeneity whose characteristic scale is larger than the sound horizon, unless the amplitude of the inhomogeneity is large and w is close to -1.


title: Gamma-Ray Bursts at high and very high energies

abstract: Gamma-Ray Bursts (GRBs) are extra-galactic and extremely energetic transient emissions of gamma rays, which are thought to be associated with the death of massive stars or the merger of compact objects in binary systems. Their huge luminosities involve the presence a newborn stellar-mass black hole emitting a relativistic collimated outflow, which accelerates particles and produces non-thermal emissions from the radio domain to the highest energies. In this article, I review recent progresses in the understanding of GRB jet physics above 100 MeV, based on Fermi observations of bright GRBs. I discuss the physical implications of these observations and their impact on GRB modeling, and I present some prospects for GRB observation at very high energies in the near future.


title: Non-homogeneous random walks on a half strip with generalized Lamperti drifts

abstract: We study a Markov chain on $\mathbb{R}_+ \times S$, where $\mathbb{R}_+$ is the non-negative real numbers and $S$ is a finite set, in which when the $\mathbb{R}_+$-coordinate is large, the $S$-coordinate of the process is approximately Markov with stationary distribution $\pi_i$ on $S$. If $\mu_i(x)$ is the mean drift of the $\mathbb{R}_+$-coordinate of the process at $(x,i) \in \mathbb{R}_+ \times S$, we study the case where $\sum_{i} \pi_i \mu_i (x) \to 0$, which is the critical regime for the recurrence-transience phase transition. If $\mu_i(x) \to 0$ for all $i$, it is natural to study the Lamperti case where $\mu_i(x) = O(1/x)$; in that case the recurrence classification is known, but we prove new results on existence and non-existence of moments of return times. If $\mu_i (x) \to d_i$ for $d_i \neq 0$ for at least some $i$, then it is natural to study the generalized Lamperti case where $\mu_i (x) = d_i + O (1/x)$. By exploiting a transformation which maps the generalized Lamperti case to the Lamperti case, we obtain a recurrence classification and existence of moments results for the former. The generalized Lamperti case is seen to be more subtle, as the recurrence classification depends on correlation terms between the two coordinates of the process.


title: Trigonometric dictionary based codec for music compression with high quality recovery

abstract: A codec for compression of music signals is proposed. The method belongs to the class of transform lossy compression. It is conceived to be applied in the high quality recovery range though. The transformation, endowing the codec with its distinctive feature, relies on the ability to construct high quality sparse approximation of music signals. This is achieved by a redundant trigonometric dictionary and a dedicated pursuit strategy. The potential of the approach is illustrated by comparison with the OGG Vorbis format, on a sample consisting of clips of melodic music. The comparison evidences remarkable improvements in compression performance for the identical quality of the decompressed signal.


title: Dynamical polaron ansatz: a theoretical tool for the ultra-strong coupling regime of circuit QED

abstract: In this work we develop a semi-analytical variational ansatz to study the properties of few photon excitations interacting with a collection of quantum emitters in regimes that go beyond the rotating wave approximation. This method can be used to approximate both the static and dynamical properties of a superconducting qubit in an open transmission line, including the spontaneous emission spectrum and the resonances in scattering experiments. The approximations are quantitatively accurate for rather strong couplings, as shown by a direct comparison to Matrix-Product-State numerical methods, and provide also a good qualitative description for stronger couplings well beyond the Markovian regime.


title: Holomorphic functions on the symmetrized bidisk II -- Interpolating sequences and the Toeplitz corona theorem

abstract: This note contains two new theorems about bounded holomorphic functions on the symmetrized bidisk -- a characterization of interpolating sequences and a Toeplitz corona theorem.


title: Numerical stability of iterative refinement with a relaxation for linear systems

abstract: Stability analysis of Wilkinson's iterative refinement with a relaxation IR(omega) for solving linear systems is given. It extends existing results for omega=1, i.e., for Wilkinson's iterative refinement. We assume that all computations are performed in fixed (working) precision arithmetic. Numerical tests were done in MATLAB to illustrate our theoretical results. A particular emphasis is given on convergence of iterative refinement with a relaxation. Our tests confirm that the choice omega=1 is the best choice from the point of numerical stability.


title: Task-related edge density (TED) - a new method for revealing large-scale network formation in fMRI data of the human brain

abstract: The formation of transient networks in response to external stimuli or as a reflection of internal cognitive processes is a hallmark of human brain function. However, its identification in fMRI data of the human brain is notoriously difficult. Here we propose a new method of fMRI data analysis that tackles this problem by considering large-scale, task-related synchronisation networks. Networks consist of nodes and edges connecting them, where nodes correspond to voxels in fMRI data, and the weight of an edge is determined via task-related changes in dynamic synchronisation between their respective times series. Based on these definitions, we developed a new data analysis algorithm that identifies edges in a brain network that differentially respond in unison to a task onset and that occur in dense packs with similar characteristics. Hence, we call this approach "Task-related Edge Density" (TED). TED proved to be a very strong marker for dynamic network formation that easily lends itself to statistical analysis using large scale statistical inference. A major advantage of TED compared to other methods is that it does not depend on any specific hemodynamic response model, and it also does not require a presegmentation of the data for dimensionality reduction as it can handle large networks consisting of tens of thousands of voxels. We applied TED to fMRI data of a fingertapping task provided by the Human Connectome Project. TED revealed network-based involvement of a large number of brain areas that evaded detection using traditional GLM-based analysis. We show that our proposed method provides an entirely new window into the immense complexity of human brain function.


title: Non-compact Newton boundary and Whitney equisingularity for non-isolated singularities

abstract: In an unpublished lecture note, J. Brian\c{c}on observed that if $\{f_t\}$ is a family of isolated complex hypersurface singularities such that the Newton boundary of $f_t$ is independent of $t$ and $f_t$ is non-degenerate, then the corresponding family of hypersurfaces $\{f_t^{-1}(0)\}$ is Whitney equisingular (and hence topologically equisingular). A first generalization of this assertion to families with non-isolated singularities was given by the second author under a rather technical condition. In the present paper, we give a new generalization under a simpler condition.


title: The $2\pi$ Subsystem in Diffractively Produced $\pi^-\pi^+\pi^-$ at COMPASS

abstract: The COMPASS experiment at CERN has collected a large dataset of $50$ million $\pi^-\pi^+\pi^-$ events produced diffractively from a proton target using a $190\,\mathrm{GeV}/c$ pion beam. The partial-wave analysis (PWA) of these high-precision data reveals previously unseen details but is limited in parts by systematic effects. The PWA is based on the isobar model, in which multi-particle decays are described as a chain of subsequent two-body decays. Here, fixed mass distributions for the appearing intermediate resonances, the so-called isobars, are assumed. These shapes, which e.g. may be parametrized by Breit-Wigner amplitudes, represent prior knowledge that has to be put into the analysis model and may therefore introduce a model dependence, thus increasing systematic uncertainties. We present a novel method, which allows to extract isobar amplitudes directly from the data in a more model-independent way. As a first application, diffractively produced $\pi^-\pi^+\pi^-$ events are analyzed. Here, the focus lies in particular on the scalar $\pi^+\pi^-$ subsystem, where in a previous analysis a signal for a new axial-vector state $a_1(1420)$ was found in the $f_0(980)\pi$ decay mode.


title: A Highly Literate Approach to Ontology Building

abstract: Ontologies present an attractive technology for describing bio-medicine, because they can be shared, and have rich computational properties. However, they lack the rich expressivity of English and fit poorly with the current scientific "publish or perish" model. While, there have been attempts to combine free text and ontologies, most of these perform \textit{post-hoc} annotation of text. In this paper, we introduce our new environment which borrows from literate programming, to allow an author to co-develop both text and ontological description. We are currently using this environment to document the Karyotype Ontology which allows rich descriptions of the chromosomal complement in humans. We explore some of the advantages and difficulties of this form of ontology development.


title: In-Beam Background Suppression Shield

abstract: The long (3ms) proton pulse of the European Spallation Source (ESS) gives rise to unique and potentially high backgrounds for the instrument suite. In such a source an instrument capabilities will be limited by it's Signal to Noise (S/N) ratio. The instruments with a direct view of the moderator, which do not use a bender to help mitigate the fast neutron background, are the most challenging. For these beam lines we propose the innovative shielding of placing blocks of material directly into the guide system, which allow a minimum attenuation of the cold and thermal fluxes relative to the background suppression. This shielding configuration has been worked into a beam line model using Geant4. We study particularly the advantages of single crystal sapphire and silicon blocks .


title: OFDM Channel Estimation via Phase Retrieval

abstract: Pilot-aided channel estimation is nowadays a standard component in each wireless receiver enabling coherent transmission of complex-valued constellations, only affected by noise and interference. Whenever these disturbances are sufficiently small and long data frames are used, high data rates can be achieved and the resource overhead due to the pilots vanishes asymptotically. On the other, it is expected that for the next generation of mobile networks not only data rate is in the main focus but also low latency, short and sporadic messages, massive connectivity, distributed & adhoc processing and robustness with respect to asynchronism. Therefore a review of several well-established principles in communication has been started already. A particular implication when using complex-valued pilots is that these values have to be known at the receiver and therefore these resources can not be used simultaneously for user data. For an OFDM-like multicarrier scheme this means that pilot tones (usually placed equidistantly according to the Nyquist theorem) are allocated with globally known amplitudes and phases to reconstruct the channel impulse response. Phases are designed and allocated globally which is in contrast to a distributed infrastructure. In this work we present therefore a new phaseless pilot scheme where only pilot amplitudes need to be known at the receiver, i.e., phases are available again and can be used for various other purposes. The idea is based on a phase retrieval result for symmetrized and zero-padded magnitude Fourier measurements obtained by two of the authors. The phases on the pilot tones can now be used to carry additional user-specific data or compensate for other signal characteristics, like the PAPR.


title: A manifestly Lorentz covariant, interacting and non-commutative Dirac equation

abstract: We propose a manifestly Lorentz covariant, non-commutative Dirac equation for charged particles interacting with an electromagnetic field. The equation is formulated on the operator level, but operators are not composed through the normal operator product, but a modified product that restores the Lorentz covariance. This equation is solved for the free particle and a particle moving in a constant magnetic field. An abstract action, constructed on the operator level, that yields this equation as equation of motion is also derived. To relate this formalism to current formulations of non-commutative quantum field theories, this action is written in a coherent state basis, leading to an action in 4-dimensional Minkowski space-time. The resulting action differs from existing non-commutative actions, but still exhibits non-commutative effects through non-locality.


title: Involutions on the affine Grassmannian and moduli spaces of principal bundles

abstract: Let $G$ be a simply connected semisimple group over $\mathbb{C}$. We show that a certain involution of an open subset of the affine Grassmannian of $G$, defined previously by Achar and the author, corresponds to the action of the nontrivial Weyl group element of $\mathrm{SL}(2)$ on the framed moduli space of $\mathbb{G}_m$-equivariant principal $G$-bundles on $\mathbb{P}^2$. As a result, the fixed-point set of the involution can be partitioned into strata indexed by conjugacy classes of homomorphisms $N\to G$ where $N$ is the normalizer of $\mathbb{G}_m$ in $\mathrm{SL}(2)$. When $G=\mathrm{SL}(r)$, the strata are Nakajima quiver varieties $\mathfrak{M}_0^{\mathrm{reg}}(\mathbf{v},\mathbf{w})$ of type D.


title: The Fixed Initial Credit Problem for Partial-Observation Energy Games is Ack-complete

abstract: In this paper we study two-player games with asymmetric partial observation and an energy objective. Such games are played on a weighted automaton by Eve, choosing actions, and Adam, choosing a transition labelled with the given action. Eve attempts to maintain the sum of the weights (of the transitions taken) non-negative while Adam tries to do the opposite. Eve does not know the exact state of the game, she is only given an equivalence class of states which contains it. In contrast, Adam has full observation. We show the fixed initial credit problem for these games is Ack-complete.


title: Gas expulsion in massive star clusters? Constraints from observations of young and gas-free objects

abstract: Gas expulsion is a central concept in some of the models for multiple populations and the light-element anticorrelations in globular clusters. If the star formation efficiency was around 30 per cent and the gas expulsion happened on the crossing timescale, this process could expel preferentially stars born with the chemical composition of the proto-cluster gas, while stars with special composition born in the centre would remain bound. Recently, a sample of extragalactic, gas-free, young massive clusters has been identified that has the potential to test the conditions for gas expulsion. We compute a large number of thin shell models, and calculate if the Rayleigh-Taylor instability is able to disrupt the shell before it reaches the escape speed. We show that the success of gas expulsion depends on the compactness index of a star cluster C5, proportionate to stellar mass over half-mass radius. For given C5, a certain critical, local star formation efficiency is required to remove the rest of the gas. Common stellar feedback processes may not lead to gas expulsion with significant loss of stars above C5 = 1. Considering pulsar winds and hypernovae, the limit increases to C5 = 30. If successful, gas expulsion generally takes place on the crossing timescale. Some observed young massive clusters have 1 < C5 < 10 and are gas-free at 10 Myr. This suggests that gas expulsion does not affect their stellar mass significantly, unless powerful pulsar winds and hypernovae are common in such objects. By comparison to observations, we show that C5 is a better predictor for the expression of multiple populations than stellar mass. The best separation between star clusters with and without multiple populations is achieved by a stellar winds-based gas expulsion model, where gas expulsion would occur exclusively in star clusters without multiple populations.


title: Particle acceleration, magnetization and radiation in relativistic shocks

abstract: What are the mechanisms of particle acceleration and radiation, as well as magnetic field build up and decay in relativistic shocks are open questions with important implications to various phenomena in high energy astrophysics. While the Weibel instability is possibly responsible for magnetic field build up and diffusive shock acceleration is a model for acceleration, both have problems and current PIC simulation show that particles are accelerated only under special conditions and the magnetic field decays on a short length scale. We present here a novel model for the structure and the emission of highly relativistic collisionless shocks. The model takes into account (and is based on) non-local energy and momentum transport across the shock front via emission and absorption of high-energy photons. This leads to a pre-acceleration of the fluid and pre-amplificaiton of the magnetic fields in the upstream region. Both have drastic implications on the shock structure. The model explains the persistence of the shock generated magnetic field at large distances from the shock front. The dissipation of this magnetic field results in a continuous particle acceleration within the downstream region. The model suggests two non-uniform emission zones (the downstream and the upstream), that give rise to three emission components with different spectral and temporal properties. A unique feature of the model is the existence of an "attractor", toward which any shock will evolve. This enables us to estimate from first principles the synchrotron and inverse Compton spectrum of the downstream emission. The model is applicable to any relativistic shock, but its distinctive features show up only for large compactness. We demonstrate that prompt and afterglow Gamma-Ray Bursts' shocks satisfy the relevant conditions and compare their observations with the predictions.


title: Measurement of photoelectron yield of the prototype for the CDEX-10 liquid argon detector

abstract: The China Dark Matter Experiment (CDEX) is a low background experiment at China Jinping Underground Laboratory (CJPL) designed to directly detect dark matter with a high-purity Germanium (HPGe) detector. In the second phase CDEX-10 with a 10 kg Germanium array detector system, the liquid argon (LAr) anti-compton active shielding and cooling system is proposed. For purpose of studying the properties of LAr detector, a prototype with an active volume of 7 liters of liquid argon was built and operated. The photoelectron yields, as a critically important parameter for the prototype detector, have been measured to be 0.051-0.079 p.e./keV for 662 keV Gamma lines at different positions. The good agreement between the experimental and simulation results has provided a quite reasonable understanding and determination of the important parameters such as the Surviving Fraction of the Ar2 excimers, the absorption length for 128 nm photons in liquid argon, the reflectivity of Teflon and so on.


title: High-scale SUSY from an R-invariant New Inflation in the Landscape

abstract: We provide an anthropic reason that the supersymmetry breaking scale is much higher than the electroweak scale as indicated by the null result of collider experiments and observed 125 GeV Higgs boson. We focus on a new inflation model as a typical low-scale inflation model that may be expected in the string landscape. In this model, the R-symmetry is broken at the minimum of the inflaton potential and its breaking scale is related to the reheating temperature. Once we admit that the anthropic principle requires thermal leptogenesis, we obtain a lower bound on gravitino mass, which is related to R-symmetry breaking scale. This scenario and resulting gravitino mass predict the consistent amplitude of density perturbations. We also find that string axions and saxions are consistently implemented in this scenario.


title: Fredholm operators on $C^*$-algebras

abstract: The aim of this note is to generalize the notion of Fredholm operator to an arbitrary $C^*$-algebra. Namely, we define "finite type" elements in an axiomatic way, and also we define Fredholm type element $a$ as such element of a given $C^*$-algebra for which there are finite type elements $p$ and $q$ such that $(1-q)a(1-p)$ is "invertible". We derive index theorem for such operators. In applications we show that classical Fredholm operators on a Hilbert space, Fredholm operators in the sense of Breuer, Atiyah and Singer on a properly infinite von Neumann algebra, and Fredholm operators on Hilbert $C^*$-modules over an unital $C^*$-algebra in the sense of Mishchenko and Fomenko are special cases of our theory.


title: Analysis of electrolyte transport through charged nanopores

abstract: We revisit the classical problem of flow of electrolyte solutions through charged capillary nanopores or nanotubes as described by the capillary pore model (also called "space charge" theory). This theory assumes very long and thin pores and uses a one-dimensional flux-force formalism which relates fluxes (electrical current, salt flux, fluid velocity) and driving forces (difference in electric potential, salt concentration, pressure). We analyze the general case with overlapping electric double layers in the pore and a nonzero axial salt concentration gradient. The 3$\times$3 matrix relating these quantities exhibits Onsager symmetry and we report a significant new simplification for the diagonal element relating axial salt flux to the gradient in chemical potential. We prove that Onsager symmetry is preserved under changes of variables, which we illustrate by transformation to a different flux-force matrix given by Gross and Osterle (1968). The capillary pore model is well-suited to describe the nonlinear response of charged membranes or nanofluidic devices for electrokinetic energy conversion and water desalination, as long as the transverse ion profiles remain in local quasi-equilibrium. As an example, we evaluate electrical power production from a salt concentration difference by reverse electrodialysis, using an efficiency vs. power diagram. We show that since the capillary pore model allows for axial gradients in salt concentration, partial loops in current, salt flux or fluid flow can develop in the pore. Predictions for macroscopic transport properties using a reduced model where the potential and concentration are assumed to be invariant with radial coordinate ("uniform potential" or "fine capillary pore" model), are close to results of the full model.


title: Pseudo-exponential maps, variants, and quasiminimality

abstract: We give a construction of quasiminimal fields equipped with pseudo-analytic maps, generalising Zilber's pseudo-exponential function. In particular we construct pseudo-exponential maps of simple abelian varieties, including pseudo-$\wp$-functions for elliptic curves. We show that the complex field with the corresponding analytic function is isomorphic to the pseudo-analytic version if and only the appropriate version of Schanuel's conjecture is true and the corresponding version of the strong exponential-algebraic closedness property holds. Moreover, we relativize the construction to build a model over a fairly arbitrary countable subfield and deduce that the complex exponential field is quasiminimal if it is exponentially-algebraically closed. This property asks only that the graph of exponentiation have non-trivial intersection with certain algebraic varieties but does not require genericity of these points. Furthermore Schanuel's conjecture is not required as a condition for quasiminimality.


title: Velocity profile of granular flows down a heap described by dimensional analysis

abstract: Gravity-driven thick granular flows are relevant to many industrial and geophysical processes. In particular, it is important to know and understand the particle velocity distributions as we get deeper into the flow from the free surface. In this paper, we use dimensional analysis as a tool to reproduce the velocity profile experimentally reported for granular flows down a confined heap for the so-called flowing layer and the creep layer underneath: the grains velocity first decrease linearly from the free surface, and then exponentially.


title: The XXL Survey. II. The bright cluster sample: catalogue and luminosity function

abstract: Context. The XXL Survey is the largest survey carried out by the XMM-Newton satellite and covers a total area of 50 square degrees distributed over two fields. It primarily aims at investigating the large-scale structures of the Universe using the distribution of galaxy clusters and active galactic nuclei as tracers of the matter distribution. Aims. This article presents the XXL bright cluster sample, a subsample of 100 galaxy clusters selected from the full XXL catalogue by setting a lower limit of $3\times 10^{-14}\,\mathrm{erg \,s^{-1}cm^{-2}}$ on the source flux within a 1$^{\prime}$ aperture. Methods. The selection function was estimated using a mixture of Monte Carlo simulations and analytical recipes that closely reproduce the source selection process. An extensive spectroscopic follow-up provided redshifts for 97 of the 100 clusters. We derived accurate X-ray parameters for all the sources. Scaling relations were self-consistently derived from the same sample in other publications of the series. On this basis, we study the number density, luminosity function, and spatial distribution of the sample. Results. The bright cluster sample consists of systems with masses between $M_{500}=7\times 10^{13}$ and $3\times 10^{14} M_\odot$, mostly located between $z=0.1$ and 0.5. The observed sky density of clusters is slightly below the predictions from the WMAP9 model, and significantly below the predictions from the Planck 2015 cosmology. In general, within the current uncertainties of the cluster mass calibration, models with higher values of $\sigma_8$ and/or $\Omega_m$ appear more difficult to accommodate. We provide tight constraints on the cluster differential luminosity function and find no hint of evolution out to $z\sim1$. We also find strong evidence for the presence of large-scale structures in the XXL bright cluster sample and identify five new superclusters.


title: Opto-mechanical system for the controlled growth of granular piles

abstract: The formation of sandpiles is a subject if interest in the field of granular matter, but most experiments are unable to control the kinetic energy of the grains added to the apex of the growing pile. In this short article, we describe an opto-mechanical system to grow granular piles by adding sand from a fixed distance above the apex of the pile. To achieve that, we use a laser beam that detects the position of the apex and, as the pile grows, sends a feedback signal to a mechanical device able to control the height of the sand delivey system above the pile.


title: Multistability of Phase-Locking and Topological Winding Numbers in Locally Coupled Kuramoto Models on Single-Loop Networks

abstract: Determining the number of stable phase-locked solutions for locally coupled Kuramoto models is a long-standing mathematical problem with important implications in biology, condensed matter physics and electrical engineering among others. We investigate Kuramoto models on networks with various topologies and show that different phase-locked solutions are related to one another by loop currents. The latter take only discrete values, as they are characterized by topological winding numbers. This result is generically valid for any network, and also applies beyond the Kuramoto model, as long as the coupling between oscillators is antisymmetric in the oscillators' coordinates. Motivated by these results we further investigate loop currents in Kuramoto-like models. We consider loop currents in nonoriented $n$-node cycle networks with nearest-neighbor coupling. Amplifying on earlier works, we give an algebraic upper bound $\mathcal{N} \le 2 \, {\rm Int}[n/4]+1$ for the number $\cal N$ of different, linearly stable phase-locked solutions. We show that the number of different stable solutions monotonically decreases as the coupling strength is decreased. Furthermore stable solutions with a single angle difference exceeding $\pi/2$ emerge as the coupling constant $K$ is reduced, as smooth continuations of solutions with all angle differences smaller than $\pi/2$ at higher $K$. In a cycle network with nearest-neighbor coupling we further show that phase-locked solutions with two or more angle differences larger than $\pi/2$ are all linearly unstable. We point out similarities between loop currents and vortices in superfluids and superconductors as well as persistent currents in superconducting rings and two-dimensional Josephson junction arrays.