title: Highly responsive ground state of PbTaSe$_2$: structural phase transition and evolution of superconductivity under pressure

abstract: Transport and magnetic studies of PbTaSe$_2$ under pressure suggest existence of two superconducting phases with the low temperature phase boundary at $\sim 0.25$ GPa that is defined by a very sharp, first order, phase transition. The first order phase transition line can be followed via pressure dependent resistivity measurements, and is found to be near 0.12 GPa near room temperature. Transmission electron microscopy and x-ray diffraction at elevated temperatures confirm that this first order phase transition is structural and occurs at ambient pressure near $\sim 425$ K. The new, high temperature / high pressure phase has a similar crystal structure and slightly lower unit cell volume relative to the ambient pressure, room temperature structure. Based on first-principles calculations this structure is suggested to be obtained by shifting the Pb atoms from the $1a$ to $1e$ Wyckoff position without changing the positions of Ta and Se atoms. PbTaSe$_2$ has an exceptionally pressure sensitive, structural phase transition with $\Delta T_s/\Delta P \approx - 1700$ K/GPa near 4 K, this first order transition causes an $\sim 1$ K ($\sim 25 \%$) step - like decrease in $T_c$ as pressure is increased through 0.25 GPa.


title: Fault diagnosability of data center networks

abstract: The data center networks $D_{n,k}$, proposed in 2008, has many desirable features such as high network capacity. A kind of generalization of diagnosability for network $G$ is $g$-good-neighbor diagnosability which is denoted by $t_g(G)$. Let $\kappa^g(G)$ be the $R^g$-connectivity. Lin et. al. in [IEEE Trans. on Reliability, 65 (3) (2016) 1248--1262] and Xu et. al in [Theor. Comput. Sci. 659 (2017) 53--63] gave the same problem independently that: the relationship between the $R^g$-connectivity $\kappa^g(G)$ and $t_g(G)$ of a general graph $G$ need to be studied in the future. In this paper, this open problem is solved for general regular graphs. We firstly establish the relationship of $\kappa^g(G)$ and $t_g(G)$, and obtain that $t_g(G)=\kappa^g(G)+g$ under some conditions. Secondly, we obtain the $g$-good-neighbor diagnosability of $D_{k,n}$ which are $t_g(D_{k,n})=(g+1)(k-1)+n+g$ for $1\leq g\leq n-1$ under the PMC model and the MM model, respectively. Further more, we show that $D_{k,n}$ is tightly super $(n+k-1)$-connected for $n\geq 2$ and $k\geq 2$ and we also prove that the largest connected component of the survival graph contains almost all of the remaining vertices in $D_{k,n}$ when $2k+n-2$ vertices removed.


title: Machine learning based compact photonic structure design for strong light confinement

abstract: We present a novel approach based on machine learning for designing photonic structures. In particular, we focus on strong light confinement that allows the design of an efficient free-space-to-waveguide coupler which is made of Si- slab overlying on the top of silica substrate. The learning algorithm is implemented using bitwise square Si- cells and the whole optimized device has a footprint of $\boldsymbol{2 \, \mu m \times 1\, \mu m}$, which is the smallest size ever achieved numerically. To find the effect of Si- slab thickness on the sub-wavelength focusing and strong coupling characteristics of optimized photonic structure, we carried out three-dimensional time-domain numerical calculations. Corresponding optimum values of full width at half maximum and coupling efficiency were calculated as $\boldsymbol{0.158 \lambda}$ and $\boldsymbol{-1.87\,dB}$ with slab thickness of $\boldsymbol{280nm}$. Compared to the conventional counterparts, the optimized lens and coupler designs are easy-to-fabricate via optical lithography techniques, quite compact, and can operate at telecommunication wavelengths. The outcomes of the presented study show that machine learning can be beneficial for efficient photonic designs in various potential applications such as polarization-division, beam manipulation and optical interconnects.


title: Phenomenological forecasting of disease incidence using heteroskedastic Gaussian processes: a dengue case study

abstract: In 2015 the US federal government sponsored a dengue forecasting competition using historical case data from Iquitos, Peru and San Juan, Puerto Rico. Competitors were evaluated on several aspects of out-of-sample forecasts including the targets of peak week, peak incidence during that week and total season incidence across each of several seasons. Our team was one of the top performers of that competition, outperforming all other teams in multiple targets/locals. In this paper we report on our methodology, a large component of which, surprisingly, ignores the known biology of epidemics at large---in particular relationships between dengue transmission and environmental factors---and instead relies on flexible nonparametric nonlinear Gaussian process (GP) regression fits that "memorize" the trajectories of past seasons, and then "match" the dynamics of the unfolding season to past ones in real-time. Our phenomenological approach has advantages in situations where disease dynamics are less well understood, e.g., at sites with shorter histories of disease (such as Iquitos), or where measurements and forecasts of ancillary covariates like precipitation are unavailable and/or where the strength of association with cases are as yet unknown. In particular, we show that the GP approach generally outperforms a more classical generalized linear (autoregressive) model (GLM) that we developed to utilize abundant covariate information. We illustrate variations of our method(s) on the two benchmark locales alongside a full summary of results submitted by other contest competitors.


title: Comparison of different bonding techniques for efficient strain transfer using piezoelectric actuators

abstract: In this paper strain transfer efficiencies from single crystalline piezoelectric lead magnesium niobate-lead titanate (PMN-PT) substrate to a GaAs semiconductor membrane bonded on top are investigated using state-of-the-art x-ray diffraction (XRD) techniques and finite-element-method (FEM) simulations. Two different bonding techniques are studied, namely gold-thermo-compression and polymer-based SU8 bonding. Our results show a much higher strain-transfer for the "soft" SU8 bonding in comparison to the "hard" bonding via gold-thermo-compression. A comparison between the XRD results and FEM simulations allows to explain this unexpected result with the presence of complex interface structures between the different layers.


title: Symmetry breaking for orthogonal groups and a conjecture by B. Gross and D. Prasad

abstract: We consider irreducible unitary representations $A_i$ of G=SO(n+1,1) with the same infinitesimal character as the trivial representation and representations $B_j$ of H=SO(n,1) with the same properties and discuss H-equivariant homomorphisms Hom_H($A_i,B_j$). For tempered representations our results confirm the predictions of conjectures by B. Gross and D. Prasad.


title: Un sous-syst\`eme d'un syst\`eme plat de dimension diff\'erentielle 2 est plat. A subsystem of a flat system of differential dimension 2 is flat

abstract: Un sous-syst\`eme de dimension diff\'erentielle au plus 2 d'une extension plate est plate. Si un tel syst\`eme plat est stationnaire, il admet des sorties plates ind\'ependantes du temps. A subsystem of a flat system of differential dimension at most 2 is flat. Furthermore, if such a flat system is stationary, we show that there exist flat utputs not depending on the time.


title: A new look at Bernoulli's inequality

abstract: In this work, a generalization of the well known Bernoulli inequality is obtained by using the theory of discrete fractional calculus. As far as we know our approach is novel.


title: Ultra-Flat Galaxies Selected from RFGC Catalog. III. Star Formation Rate

abstract: We examine the star formation properties of galaxies with very thin disks selected from the Revised Flat Galaxy Catalog (RFGC). The sample contains 333 ultra-flat galaxies (UFG) at high Galactic latitudes, $|b|>10^{\circ}$, with a blue major angular diameter of $a\geq 1.2^{\prime}$, blue and red apparent axial ratios of $(a/b)_b > 10$, $(a/b)_r > 8.5$ and radial velocities within 10000~km s$^{-1}$. As a control sample for them we use a population of 722 more thick RFGC galaxies with $(a/b)_b > 7$, situated in the same volume. The UFG distribution over the sky indicates them as a population of quite isolated galaxies. We found that the specific star formation rate, sSFR FUV, determined via the FUV GALEX flux, increases steadily from the early type to late type disks for both the UFG and RFGC-UFG samples, showing no significant mutual difference within each morphological type T. The population of UFG disks has the average H,I-mass-to-stellar-mass ratio by $(0.25\pm0.03)$ dex higher than that of RFGC--UFG galaxies. Being compared with arbitrary orientated disks of the same type, the ultra-flat edge-on galaxies reveal that their total H,I mass is hidden by self-absorption on the average by approximately 0.20 dex. We demonstrate that using the robust stellar mass estimate via $\langle B-K \rangle$-color and galaxy type T for the thin disks, together with a nowaday accounting for internal extinction, yields their sSFR quantities definitely lying below the limit of $-9.4$ dex,(yr$^{-1}$). The collected observational data on UFG disks imply that their average star formation rate in the past has been approximately three times the current SFR. The UFG galaxies have also sufficient amount of gas to support their observed SFR over the following nearly 9 Gyrs.


title: Multiple soliton compression stages from competing plasma nonlinearities in mid-IR gas-filled hollow-core fibers

abstract: We investigate numerically soliton-plasma interaction in a noble-gas-filled silica hollow-core anti-resonant fiber pumped in the mid-IR at 3.0 {\mu}m. We observe multiple soliton self-compression stages due to distinct stages where either the self-focusing or the self-defocusing nonlinearity dominates. Specifically, the parameters may be tuned so the competing plasma self-defocusing nonlinearity only dominates over the Kerr self-focusing nonlinearity around the soliton self-compression stage, where the increasing peak intensity on the leading pulse edge initiates a competing self-defocusing plasma nonlinearity acting nonlocally on the trailing edge, effectively preventing soliton-formation there. As the plasma switches off after the self-compression stage, self-focusing dominates again, initiating another soliton self-compression stage in the trailing edge. This process is accompanied by supercontinuum generation spanning 1-4 {\mu}m. The technique could be exploited to generate an ultrafast sequence of several few-cycle pulses.


title: Proof Diagrams for Multiplicative Linear Logic: Syntax and Semantics

abstract: Proof nets are a syntax for linear logic proofs which gives a coarser notion of proof equivalence with respect to syntactic equality together with an intuitive geometrical representation of proofs. In this paper we give an alternative $2$-dimensional syntax for multiplicative linear logic derivations. The syntax of string diagrams authorizes the definition of a framework where the sequentializability of a term, i.e. deciding whether the term corresponds to a correct derivation, can be verified in linear time. Furthermore, we can use this syntax to define a denotational semantics for multiplicative linear logic with units by means of equivalence classes of proof diagrams modulo a terminating rewriting.


title: Regularity of anisotropic refinable functions

abstract: This paper presents a detailed regularity analysis of anisotropic wavelet frames and subdivision. In the univariate setting, the smoothness of wavelet frames and subdivision is well understood by means of the matrix approach. In the multivariate setting, this approach has been extended only to the special case of isotropic refinement with the dilation matrix all of whose eigenvalues are equal in the absolute value. The general anisotropic case has resisted to be fully understood: the matrix approach can determine whether a refinable function belongs to $C(\mathbb{R}^s)$ or $L_p(\mathbb{R}^s)$, $1 \le p < \infty$, but its H\"older regularity remained mysteriously unattainable. It this paper we show how to compute the H\"older regularity in $C(\mathbb{R}^s)$ or $L_p(\mathbb{R}^s)$, $1 \le p < \infty$. In the anisotropic case, our expression for the exact H\"older exponent of a refinable function reflects the impact of the variable moduli of the eigenvalues of the corresponding dilation matrix. In the isotropic case, our results reduce to the well-known facts from the literature. We provide an efficient algorithm for determining the finite set of the restricted transition matrices whose spectral properties characterize the H\"older exponent of the corresponding refinable function. We also analyze the higher regularity, the local regularity, the corresponding moduli of continuity, and the rate of convergence of the corresponding subdivision schemes. We illustrate our results with several examples.


title: Optical selection rules for excitonic Rydberg series in the massive Dirac cones of hexagonal 2D materials

abstract: We investigate the optical transition selection rules for excitonic Rydberg series formed in massive Dirac cones. The entanglement of the exciton envelop function with the pseudospin texture leads to anomalous selection rules for one-photon generation of excitons, where $d$-orbitals can be excited with the opposite helicity selection rule from the $s$-orbitals in a given valley. The trigonal warping effects in realistic hexagonal lattices further renders more excited states bright, where $p$-orbitals can also be accessed by one-photon excitation with the opposite valley selection rules to the $s$-orbitals. The one-photon generation of exciton in the various states and the intra-excitonic transition between these states are both dictated by the discrete in-plane rotational symmetry of the lattices, and our results show that in hexagonal 2D materials the symmetry allowed transitions are enabled when trigonal warping effects are included in the massive Dirac fermion model. In monolayer transition metal dichalcogenides where excitons can be generated by visible light and intra-excitonic transitions can be induced by infrared light, we give the strength of these optical transitions, estimated using modified hydrogen-like envelope functions combined with the optical transition matrix elements between the Bloch states calculated at various $k$ points.


title: Project of refurbishment of the GILDA beamline at the ESRF

abstract: This document presents the refurbishment project for GILDA, the Italian CRG beamline at the ESRF. After a short introduction to the beamline and its achievements in the last years, the base ideas and concepts driving the design are exposed whereas a more detailed technical description of the new instrument is given in a further section. An analysis of the costs and a schedule for the implementation of the program is also presented in the last part of the document.


title: Generic Dynamical Phase Transition in One-Dimensional Bulk-Driven Lattice Gases with Exclusion

abstract: Dynamical phase transitions are crucial features of the fluctuations of statistical systems, corresponding to boundaries between qualitatively different mechanisms of maintaining unlikely values of dynamical observables over long periods of time. They manifest themselves in the form of non-analyticities in the large deviation function of those observables. In this paper, we look at bulk-driven exclusion processes with open boundaries. It is known that the standard asymmetric simple exclusion process exhibits a dynamical phase transition in the large deviations of the current of particles flowing through it. That phase transition has been described thanks to specific calculation methods relying on the model being exactly solvable, but more general methods have also been used to describe the extreme large deviations of that current, far from the phase transition. We extend those methods to a large class of models based on the ASEP, where we add arbitrary spatial inhomogeneities in the rates and short-range potentials between the particles. We show that, as for the regular ASEP, the large deviation function of the current scales differently with the size of the system if one considers very high or very low currents, pointing to the existence of a dynamical phase transition between those two regimes: high current large deviations are extensive in the system size, and the typical states associated to them are Coulomb gases, which are correlated ; low current large deviations do not depend on the system size, and the typical states associated to them are anti-shocks, consistently with a hydrodynamic behaviour. Finally, we illustrate our results numerically on a simple example, and we interpret the transition in terms of the current pushing beyond its maximal hydrodynamic value, as well as relate it to the appearance of Tracy-Widom distributions in the relaxation statistics of such models.


title: A new magnetic phase in the nickelate perovskite TlNiO$_3$

abstract: The RNiO$_3$ perovskites are known to order antiferromagnetically below a material-dependent N\'eel temperature $T_\text{N}$. We report experimental evidence indicating the existence of a second magnetically-ordered phase in TlNiO$_3$ above $T_\text{N} = 104$ K, obtained using nuclear magnetic resonance and muon spin rotation spectroscopy. The new phase, which persists up to a temperature $T_\text{N}^* = 202$ K, is suppressed by the application of an external magnetic field of approximately 1 T. It is not yet known if such a phase also exists in other perovskite nickelates.


title: An ultraproduct method via left reversible semigroups to study Bruck's generalized conjecture

abstract: We use a method similar to ultraproducts to study the common fixed point of a left reversible semitopological semigroup acting on a Banach space. As an application, we prove a Bruck-type theorem for nearly uniformly convex Banach spaces to the effect that such spaces have weak fixed point property for left reversible semigroups.


title: The Average Dimension of the Hermitian Hull of Constayclic Codes over Finite Fields

abstract: The hulls of linear and cyclic codes have been extensively studied due to their wide applications. In this paper, the average dimension of the Hermitian hull of constacyclic codes of length $n$ over a finite field $\mathbb{F}_{q^2}$ is determined together with some upper and lower bounds. It turns out that either the average dimension of the Hermitian hull of constacyclic codes of length $n$ over $\mathbb{F}_{q^2}$ is zero or it grows the same rate as $n$. Comparison to the average dimension of the Euclidean hull of cyclic codes is discussed as well.


title: New Beam Tracking Technique for Millimeter Wave-band Communications

abstract: In this paper, we propose an efficient beam tracking method for mobility scenario in mmWave-band communications. When the position of the mobile changes in mobility scenario, the base-station needs to perform beam training frequently to track the time-varying channel, thereby spending significant resources for training beams. In order to reduce the training overhead, we propose a new beam training approach called "beam tracking" which exploits the continuous nature of time varying angle of departure (AoD) for beam selection. We show that transmission of only two training beams is enough to track the time-varying AoD at good accuracy. We derive the optimal selection of beam pair which minimizes Cramer-Rao Lower Bound (CRLB) for AoD estimation averaged over statistical distribution of the AoD. Our numerical results demonstrate that the proposed beam tracking scheme produces better AoD estimation than the conventional beam training protocol with less training overhead.


title: On self-affine sets

abstract: We survey the dimension theory of self-affine sets for general mathematical audience. The article is in Finnish.


title: A level training set with both a computer-based control and a compact controller

abstract: In engineering education, the combination with theoretical education and practical education is an essential problem.


title: Gravitational Focusing and the Star Cluster Initial Mass Function

abstract: We discuss the possibility that gravitational focusing, is responsible for the power-law mass function of star clusters $N(\log M) \propto M^{-1}$. This power law can be produced asymptotically when the mass accretion rate of an object depends upon the mass of the accreting body as $\dot{M} \propto M^2$. While Bondi-Hoyle-Littleton accretion formally produces this dependence on mass in a uniform medium, realistic environments are much more complicated. However, numerical simulations in SPH allowing for sink formation yield such an asymptotic power-law mass function. We perform pure N-body simulations to isolate the effects of gravity from those of gas physics and to show that clusters naturally result with the power-law mass distribution. We also consider the physical conditions necessary to produce clusters on appropriate timescales. Our results help support the idea that gravitationally-dominated accretion is the most likely mechanism for producing the cluster mass function.


title: Two classes of quadratic vector fields for which the Kahan discretization is integrable

abstract: Applying Kahan's discretization to the reduced Nahm equations, we obtain two classes of integrable mappings.


title: The stifness of the supranuclear equation of state (once again)

abstract: We revisit the present status of the stiffness of the supranuclear equations of state, particularly the solutions that increase the stiffness in the presence of hyperons, the putative transition to a quark matter phase and the robustness of massive compact star observations.


title: Correspondence on "Single-shot simulations of dynamic quantum many-body systems"

abstract: Comment on a letter to Nature Physics, where Sakmann and Kasevich claim to solve the many-body time dependent Schr\"odinger equation to simulate single experimental runs of interacting quantum systems.


title: Initial simulation study on high-precision radio measurements of the depth of shower maximum with SKA1-low

abstract: As LOFAR has shown, using a dense array of radio antennas for detecting extensive air showers initiated by cosmic rays in the Earth's atmosphere makes it possible to measure the depth of shower maximum for individual showers with a statistical uncertainty less than $20\,g/cm^2$. This allows detailed studies of the mass composition in the energy region around $10^{17}\,eV$ where the transition from a Galactic to an Extragalactic origin could occur. Since SKA1-low will provide a much denser and very homogeneous antenna array with a large bandwidth of $50-350\,MHz$ it is expected to reach an uncertainty on the $X_{\max}$ reconstruction of less than $10\,g/cm^2$. We present first results of a simulation study with focus on the potential to reconstruct the depth of shower maximum for individual showers to be measured with SKA1-low. In addition, possible influences of various parameters such as the numbers of antennas included in the analysis or the considered frequency bandwidth will be discussed.


title: Musical intervals under 12-note equal temperament: a geometrical interpretation

abstract: Musical intervals in multiple of semitones under 12-note equal temperament, or more specifically pitch-class subsets of assigned cardinality ($n$-chords) are conceived as positive integer points within an Euclidean $n$-space. The number of distinct $n$-chords is inferred from combinatorics with the extension to $n=0$, involving an Euclidean 0-space. The number of repeating $n$-chords, or points which are turned into themselves during a circular permutation, $T_n$, of their coordinates, is inferred from algebraic considerations. Finally, the total number of $n$-chords and the number of $T_n$ set classes are determined. Palindrome and pseudo palindrome $n$-chords are defined and included among repeating $n$-chords, with regard to an equivalence relation, $T_n/T_nI$, where reflection is added to circular permutation. To this respect, the number of $T_n$ set classes is inferred concerning palindrome and pseudo palindrome $n$-chords and the remaining $n$-chords. The above results are reproduced within the framework of a geometrical interpretation, where positive integer points related to $n$-chords of cardinality, $n$, belong to a regular inclined $n$-hedron, $\Psi_{12}^n$, the vertexes lying on the coordinate axes of a Cartesian orthogonal reference frame at a distance, $x_i=12$, $1\le i\le n$, from the origin. Considering $\Psi_{12}^n$ as special cases of lattice polytopes, the number of related nonnegative integer points is also determined for completeness. A comparison is performed with the results inferred from group theory.


title: Paley and the Paley graphs

abstract: This paper discusses some aspects of the history of the Paley graphs and their automorphism groups.


title: DARE Mission Design: Low RFI Observations from a Low-Altitude Frozen Lunar Orbit

abstract: The Dark Ages Radio Experiment (DARE) seeks to study the cosmic Dark Ages approximately 80 to 420 million years after the Big Bang. Observations require truly quiet radio conditions, shielded from Sun and Earth electromagnetic (EM) emissions, on the far side of the Moon. DARE's science orbit is a frozen orbit with respect to lunar gravitational perturbations. The altitude and orientation of the orbit remain nearly fixed indefinitely, maximizing science time without the need for maintenance. DARE's observation targets avoid the galactic center and enable investigation of the universe's first stars and galaxies.


title: Targeting pure quantum states by strong noncommutative dissipation

abstract: We propose a solution to the problem of realizing a predefined and arbitrary pure quantum state, based on the simultaneous presence of coherent and dissipative dynamics, noncommuting on the target state and in the limit of strong dissipation. More precisely, we obtain a necessary and sufficient criterion whereby the nonequilibrium steady state (NESS) of an open quantum system described by a Lindblad master equation approaches a target pure state in the Zeno regime, i.e., for infinitely large dissipative coupling. We also provide an explicit formula for the characteristic dissipative strength beyond which the purity of the NESS becomes effective, thus paving the way to an experimental implementation of our criterion. For an illustration, we deal with targeting a Bell state, an arbitrary pure state of $N$ qubits, and a spin-helix state of $N$ qubits.

categories: quant-ph cond-mat.stat-mech math-ph math.MP


title: Low-Dose CT with a Residual Encoder-Decoder Convolutional Neural Network (RED-CNN)

abstract: Given the potential X-ray radiation risk to the patient, low-dose CT has attracted a considerable interest in the medical imaging field. The current main stream low-dose CT methods include vendor-specific sinogram domain filtration and iterative reconstruction, but they need to access original raw data whose formats are not transparent to most users. Due to the difficulty of modeling the statistical characteristics in the image domain, the existing methods for directly processing reconstructed images cannot eliminate image noise very well while keeping structural details. Inspired by the idea of deep learning, here we combine the autoencoder, the deconvolution network, and shortcut connections into the residual encoder-decoder convolutional neural network (RED-CNN) for low-dose CT imaging. After patch-based training, the proposed RED-CNN achieves a competitive performance relative to the-state-of-art methods in both simulated and clinical cases. Especially, our method has been favorably evaluated in terms of noise suppression, structural preservation and lesion detection.


title: Local oscillations in moderately dense sequences of primes

abstract: The distribution of differences of consecutive members of sequences of primes is investigated. A quantitative measure for oscillations among these differences is the curvature of the sequence. If the sequence is not too sparse, then sharp estimates for its curvature are provided.


title: Attacking the V: On the Resiliency of Adaptive-Horizon MPC

abstract: We introduce the concept of a V-formation game between a controller and an attacker, where controller's goal is to maneuver the plant (a simple model of flocking dynamics) into a V-formation, and the goal of the attacker is to prevent the controller from doing so. Controllers in V-formation games utilize a new formulation of model-predictive control we call Adaptive-Horizon MPC (AMPC), giving them extraordinary power: we prove that under certain controllability assumptions, an AMPC controller is able to attain V-formation with probability 1. We define several classes of attackers, including those that in one move can remove R birds from the flock, or introduce random displacement into flock dynamics. We consider both naive attackers, whose strategies are purely probabilistic, and AMPC-enabled attackers, putting them on par strategically with the controllers. While an AMPC-enabled controller is expected to win every game with probability 1, in practice, it is resource-constrained: its maximum prediction horizon and the maximum number of game execution steps are fixed. Under these conditions, an attacker has a much better chance of winning a V-formation game. Our extensive performance evaluation of V-formation games uses statistical model checking to estimate the probability an attacker can thwart the controller. Our results show that for the bird-removal game with R = 1, the controller almost always wins (restores the flock to a V-formation). For R = 2, the game outcome critically depends on which two birds are removed. For the displacement game, our results again demonstrate that an intelligent attacker, i.e. one that uses AMPC in this case, significantly outperforms its naive counterpart that randomly executes its attack.


title: (G, \mu)-displays and Rapoport-Zink spaces

abstract: Let (G, \mu) be a pair of a reductive group G over the p-adic integers and a minuscule cocharacter {\mu} of G defined over an unramified extension. We introduce and study "(G, \mu)-displays" which generalize Zink's Witt vector displays. We use these to define certain Rapoport-Zink formal schemes purely group theoretically, i.e. without p-divisible groups.


title: Modeling of radio emission from a particle cascade in a magnetic field and its experimental validation

abstract: The SLAC T-510 experiment was designed to compare controlled laboratory measurements of radio emission of particle showers to predictions using particle-level simulations, which are relied upon in ultra-high-energy cosmic-ray air shower detection. Established formalisms for the simulation of radio emission physics, the "endpoint" formalism and the "ZHS" formalism, lead to results which can be explained by a superposition of magnetically induced transverse current radiation and charge-excess radiation due to the Askaryan effect. Here, we present the results of Geant4 simulations for the SLAC T-510 experiment, taking into account the details of the experimental setup (beam energy, target geometry and material, magnetic field configuration, and refraction effects) and their comparison to measured data with respect to e.g. signal polarisation, linearity with magnetic field, and angular distribution. We find that the microscopic calculations reproduce the measurements within uncertainties and describe the data well.


title: Note on "Average resistance of toroidal graphs" by Rossi, Frasca and Fagnani

abstract: In our recent paper W.S. Rossi, P. Frasca and F. Fagnani, "Average resistance of toroidal graphs", SIAM Journal on Control and Optimization, 53(4):2541--2557, 2015, we studied how the average resistances of $d$-dimensional toroidal grids depend on the graph topology and on the dimension of the graph. Our results were based on the connection between resistance and Laplacian eigenvalues. In this note, we contextualize our work in the body of literature about random walks on graphs. Indeed, the average effective resistance of the $d$-dimensional toroidal grid is proportional to the mean hitting time of the simple random walk on that grid. If $d\geq3 $, then the average resistance can be bounded uniformly in the number of nodes and its value is of order $1/d$ for large $d$.


title: Poincar\'e-Einstein metrics and Yamabe invariants

abstract: In this note we prove the existence of infinitely many positive conformal classes on $S^7$ which cannot be the conformal infinity of a Poincar\'e-Einstein metric on the ball $B^8$. We also prove a sharp inequality between the Yamabe invariant of the conformal infinity and the Yamabe invariant of the interior (after a suitable compactification).


title: Spin-polarized current in non-collinear antiferromagnets

abstract: Noncollinear antiferromagnets, such as Mn$_3$Sn and Mn$_3$Ir, were recently shown to be analogous to ferromagnets in that they have a large anomalous Hall effect. Here we show that these materials are similar to ferromagnets in another aspect: the charge current in these materials is spin-polarized. In addition, we show that the same mechanism that leads to the spin-polarized current also leads to a transverse spin current, which has a distinct symmetry and origin from the conventional spin Hall effect. We illustrate the existence of the spin-polarized current and the transverse spin current by performing \emph{ab initio} microscopic calculations and by analyzing the symmetry. We discuss possible applications of these novel spin currents, such as an antiferromagnetic metallic or tunneling junction.


title: Computational Efficiency and Amdahl's law for the Adaptive Resolution Simulation Technique

abstract: We discuss the computational performance of the adaptive resolution technique in molecular simulation when it is compared with equivalent full coarse-grained and full atomistic simulations. We show that an estimate of its efficiency, within $10-15\%$ accuracy, is given by the Amdahl's Law adapted to the specific quantities involved in the problem. The derivation of the predictive formula is general enough that it may be applied to the general case of molecular dynamics approaches where a reduction of degrees of freedom in a multiscale fashion occurs.


title: Mean square in the prime geodesic theorem

abstract: We prove upper bounds for the mean square of the remainder in the prime geodesic theorem, for every cofinite Fuchsian group, which improve on average on the best known pointwise bounds. The proof relies on the Selberg trace formula. For the modular group we prove a refined upper bound by using the Kuznetsov trace formula.


title: Cascading Failures in Interdependent Systems: Impact of Degree Variability and Dependence

abstract: We study cascading failures in a system comprising interdependent networks/systems, in which nodes rely on other nodes both in the same system and in other systems to perform their function. The (inter-)dependence among nodes is modeled using a dependence graph, where the degree vector of a node determines the number of other nodes it can potentially cause to fail in each system through aforementioned dependency. In particular, we examine the impact of the variability and dependence properties of node degrees on the probability of cascading failures. We show that larger variability in node degrees hampers widespread failures in the system, starting with random failures. Similarly, positive correlations in node degrees make it harder to set off an epidemic of failures, thereby rendering the system more robust against random failures.


title: A charged anisotropic well-behaved Adler-Finch-Skea solution Satisfying Karmarkar Condition

abstract: In the present article, we discover a new well-behaved charged anisotropic solution of Einstein-Maxwell's field equations. We ansatz the metric potential $g_{00}$ of the form given by Maurya el al. (arXiv:1607.05582v1) with $n=2$. In their article it is mentioned that for $n=2$ the solution is not well-behaved for neutral configuration as the speed of sound is non-decreasing radially outward. However, the solution can represent a physically possible configuration with the inclusion of some net electric charged i.e. the solution can become a well-behaved solution with decreasing sound speed radially outward for a charged configuration. Due to the inclusion of electric charged the solution leads to a very stiff equation of state (EoS) with the velocity of sound at the center $v_{r0}^2=0.819, ~v_{t0}^2=0.923$ and the compactness parameter $u=0.823$ is closed to the Buchdahl limit 0.889. This stiff EoS support a compact star configuration of mass $5.418M_\odot$ and radius of $10.1 km$.


title: Phenomenological Universalities: Coherence, Supersymmetry and Growth

abstract: The phenomenological universalities (PU) are extended to include time-depended quantum oscillatory phenomena, coherence and supersymmetry. It will be proved that this approach generates minimum uncertainty coherent states of time-dependent oscillators, which in the dissociation (classical) limit reduce to the functions describing growth (regression) of the systems evolving over time. The results obtained reveal existence of a new class of macroscopic quantum (or quasi-quantum) phenomena, which may play a vital role in coherent formation of the specific growth patterns in complex systems.


title: The Final Theory of Physics - a Tautology?

abstract: We acuminate the idea of a final theory of physics in order to analyze its logical implications and consequences. It is argued that the rationale of a final theory is the principle of sufficient reason. This implies that a final theory of physics, presumed such a theory is possible, does not allow to incorporate substantial (non-trivial) propositions unless they are logically or mathematically deduced. Differences between physics and mathematics are discussed with emphasis on the role of physical constants. It is shown that it is logically impossible to introduce constants on the fundamental level of a final theory. The most fundamental constants emerging within a final theory are constants of motion. It is argued that the only possibility to formulate a final theory is necessarily a tautology: A final theory of physics can only be derived from those presumptions about reality that are inherent in the idea and practice of physics itself. It is argued that a final theory is based on the notion of objectivity, but it is logically impossible that an ideal final theory supports realism.


title: Plasmonics of Topological Insulators at Optical Frequencies

abstract: The development of nanoplasmonic devices, such as plasmonic circuits and metamaterial superlenses in the visible to ultraviolet frequency range, is hampered by the lack of low-loss plasmonic media. Recently, strong plasmonic response was reported in a certain class of topological insulators. Here, we present a first-principles density functional theory analysis of the dielectric functions of topologically insulating quaternary (Bi,Sb)2(Te,Se)3 trichalcogenide compounds. Bulk plasmonic properties, dominated by interband transitions, are observed from 2-3 eV and extend to higher frequencies. Moreover, trichalcogenide compounds are better plasmonic media than gold and silver at blue and UV wavelengths. By analysing thin slabs, we also show that these materials exhibit topologically protected surface states, which are capable of supporting propagating plasmon polariton modes over an extremely broad spectral range, from the visible to the mid-infrared and beyond, owing to a combination of inter- and intra-surface band transitions.


title: A New Paradigm for Turbulent Transport Across a Steep Gradient in Toroidal Plasmas

abstract: First principle gyrokinetic simulation of the edge turbulent transport in toroidal plasmas finds a reverse trend in the turbulent transport coefficients under strong gradients. It is found that there exist both linear and nonlinear critical gradients for the nonmonotonicity of transport characteristics. The discontinuity of transport flux slope around the turning gradient shows features of second order phase transition. Under strong gradient the most unstable modes are in non-ground eigenstates with unconventional mode structures, which significantly reduces the effective correlation length and thus reverse the transport trend. Our results suggest a completely new mechanism for the low to high confinement mode transition without invoking shear flow or zonal flow.


title: Pushing for weighted tree automata

abstract: A weight normalization procedure, commonly called pushing, is introduced for weighted tree automata (wta) over commutative semifields. The normalization preserves the recognized weighted tree language even for nondeterministic wta, but it is most useful for bottom-up deterministic wta, where it can be used for minimization and equivalence testing. In both applications a careful selection of the weights to be redistributed followed by normalization allows a reduction of the general problem to the corresponding problem for bottom-up deterministic unweighted tree automata. This approach was already successfully used by Mohri and Eisner for the minimization of deterministic weighted string automata. Moreover, the new equivalence test for two wta $M$ and $M'$ runs in time $\mathcal O((\lvert M \rvert + \lvert M'\rvert) \cdot \log {(\lvert Q\rvert + \lvert Q'\rvert)})$, where $Q$ and $Q'$ are the states of $M$ and $M'$, respectively, which improves the previously best run-time $\mathcal O(\lvert M \rvert \cdot \lvert M'\rvert)$.


title: Associated production of a dilepton and a $\Upsilon(J/\psi)$ at the LHC as a probe of gluon transverse momentum dependent distributions

abstract: We discuss the impact on the study of gluon transverse momentum dependent distributions (TMDs) of the associated production of a lepton pair and a $\Upsilon$ (or a $J/\psi$) in unpolarised proton-proton collisions, $pp\to {\cal Q} \, \ell \bar{\ell} X$, at LHC energies, where one can assume that such final states are dominantly induced by gluon fusion. If the transverse momentum of the quarkonium-dilepton system - namely, the transverse momentum imbalance of the quarkonium state and the lepton pair - is small, the corresponding cross sections can be calculated within the framework of TMD factorisation. Using the helicity formalism, we show in detail how these cross sections are connected to the moments of two independent TMDs: the distribution of unpolarised gluons, $f_1^g$, and the distribution of linearly polarised gluons, $h_1^{\perp g}$. We complete our exhaustive derivation of these general relations with a phenomenological analysis of the feasibility of the TMD extraction, as well as some outlooks.


title: Inflationary $\alpha$-attractor cosmology: A global dynamical systems perspective

abstract: We study flat FLRW $\alpha$-attractor $\mathrm{E}$- and $\mathrm{T}$-models by introducing a dynamical systems framework that yields regularized unconstrained field equations on two-dimensional compact state spaces. This results in both illustrative figures and a complete description of the entire solution spaces of these models, including asymptotics. In particular, it is shown that observational viability, which requires a sufficient number of $e$-folds, is associated with a solution given by a one-dimensional center manifold of a past asymptotic de Sitter state, where the center manifold structure also explains why nearby solutions are attracted to this `inflationary attractor solution.' A center manifold expansion yields a description of the inflationary regime with arbitrary analytic accuracy, where the slow-roll approximation asymptotically describes the tangency condition of the center manifold at the asymptotic de Sitter state.

categories: gr-qc astro-ph.CO hep-th math-ph math.MP


title: Pixel-wise Ear Detection with Convolutional Encoder-Decoder Networks

abstract: Object detection and segmentation represents the basis for many tasks in computer and machine vision. In biometric recognition systems the detection of the region-of-interest (ROI) is one of the most crucial steps in the overall processing pipeline, significantly impacting the performance of the entire recognition system. Existing approaches to ear detection, for example, are commonly susceptible to the presence of severe occlusions, ear accessories or variable illumination conditions and often deteriorate in their performance if applied on ear images captured in unconstrained settings. To address these shortcomings, we present in this paper a novel ear detection technique based on convolutional encoder-decoder networks (CEDs). For our technique, we formulate the problem of ear detection as a two-class segmentation problem and train a convolutional encoder-decoder network based on the SegNet architecture to distinguish between image-pixels belonging to either the ear or the non-ear class. The output of the network is then post-processed to further refine the segmentation result and return the final locations of the ears in the input image. Different from competing techniques from the literature, our approach does not simply return a bounding box around the detected ear, but provides detailed, pixel-wise information about the location of the ears in the image. Our experiments on a dataset gathered from the web (a.k.a. in the wild) show that the proposed technique ensures good detection results in the presence of various covariate factors and significantly outperforms the existing state-of-the-art.


title: Approximate Variational Estimation for a Model of Network Formation

abstract: We develop approximate estimation methods for exponential random graph models (ERGMs), whose likelihood is proportional to an intractable normalizing constant. The usual approach approximates this constant with Monte Carlo simulations, however convergence may be exponentially slow. We propose a deterministic method, based on a variational mean-field approximation of the ERGM's normalizing constant. We compute lower and upper bounds for the approximation error for any network size, adapting nonlinear large deviations results. This translates into bounds on the distance between true likelihood and mean-field likelihood. Monte Carlo simulations suggest that in practice our deterministic method performs better than our conservative theoretical approximation bounds imply, for a large class of models.


title: Up and down quark masses and corrections to Dashen's theorem from lattice QCD and quenched QED

abstract: We present a determination of the corrections to Dashen's theorem and of the individual up and down quark masses from a lattice calculation based on quenched QED and $N_f=2+1$ QCD simulations with 5 lattice spacings down to 0.054 fm. The simulations feature lattice sizes up to 6 fm and average up-down quark masses all the way down to their physical value. For the parameter which quantifies violations to Dashens's theorem we obtain $\epsilon=0.73(2)(5)(17)$, where the first error is statistical, the second is systematic, and the third is an estimate of the QED quenching error. For the light quark masses we obtain, $m_u=2.27(6)(5)(4) \, MeV$ and $m_d=4.67(6)(5)(4) \, MeV$ in the $\overline{MS}$ scheme at $2 \, GeV$ and the isospin breaking ratios $m_u/m_d=0.485(11)(8)(14)$, $R=38.2(1.1)(0.8)(1.4)$ and $Q=23.4(0.4)(0.3)(0.4)$. Our results exclude the $m_u=0$ solution to the strong CP problem by more than 24 standard deviations.


title: Charm hadrons and lattice QCD

abstract: We review two lattice calculations involving charm hadrons: a determination of the vector and scalar form factors of the semileptonic $D \to \pi \ell \nu$ decays, which are relevant for the extraction of the CKM matrix element $\lvert V_{cd} \rvert$ from experimental data, and a calculation of the matrix elements of four-fermion operators relevant to the description of the neutral D mixing in the Standard Model and its extensions. Both analyses are based on the gauge configurations produced by the European Twisted Mass Collaboration with $N_f = 2 + 1 + 1$ flavors of dynamical quarks. We simulated at three different values of the lattice spacing and with pion masses as small as 210 MeV.


title: Transaction Support over Redis: An Overview

abstract: This document outlines the approach to supporting cross-node transactions over a Redis cluster.


title: Dynamic load balancing for large-scale adaptive finite element computation

abstract: For the parallel computation of partial differential equations, one key is the grid partitioning. It requires that each process owns the same amount of computations, and also, the partitioning quality should be proper to reduce the communications among processes. When calculating the partial differential equations using adaptive finite element methods, the grid and the basis functions adjust in each iteration, which introduce load balancing issues. The grid should be redistributed dynamically. This paper studies dynamic load balancing algorithms and the implementation on the adaptive finite element platform PHG. The numerical experiments show that algorithms studied in this paper have good partitioning quality, and they are efficient.


title: Soft Extragalactic X-Ray Binaries at the Eddington Threshold

abstract: The luminosity range at and just below the 10^39 erg/s cut-off for defining ultraluminous X-ray sources (ULXs) is a little-explored regime. It none-the-less hosts a large number of X-ray sources, and has great potential for improving our understanding of sources with ~Eddington accretion rates. We select a sample of four sources in this Eddington Threshold regime with good data for further study; these objects possess a variety of soft spectral shapes. We perform X-ray spectral and timing analysis on the XMM-Newton and Chandra data for these objects to gain insight into their accretion mechanisms, and also examine their optical counterparts using HST images. NGC 300 X-1 is a highly luminous and well-known example of the canonical steep power-law accretion state. M51 ULS exhibits a cool blackbody-like spectrum and is consistent with being an ultraluminous supersoft source (ULS), possibly a super-Eddington accreting object viewed at a high inclination through an optically thick outflowing wind. NGC 4395 ULX-1 and NGC 6946 ULX-1 have unusually steep power-law tails, for which we discuss a variety of possible physical mechanisms and links to similar features in Galactic microquasars, and we conclude that these sources are likely intermediate objects between the soft ultraluminous regime of ULXs and classic ULSs.


title: Geant4 simulations of radio signals from particle showers for the SLAC T-510 experiment

abstract: The SLAC T-510 experiment was designed to reproduce the physics of radio emission from air showers caused by ultra-high energy cosmic rays in a controlled lab experiment with the goal to test established formalisms for simulation of radio emission physics: the "end-point" formalism and the "ZHS" formalism. Simulation results derived with these formalisms can be explained by a superposition of magnetically induced transverse current radiation and the Askaryan (charge-excess) effect. Here, we present results of Geant4 simulations of the experiment with both formalisms, taking into account the details of the experimental setup (beam energy, target geometry and material, magnetic field configuration, and refraction effects) to test this hypothesis


title: Lorentzian surfaces and the curvature of the Schmidt metric

abstract: The b-boundary is a mathematical tool used to attach a topological boundary to incomplete Lorentzian manifolds using a Riemaniann metric called the Schmidt metric on the frame bundle. In this paper, we give the general form of the Schmidt metric in the case of Lorentzian surfaces. Furthermore, we write the Ricci scalar of the Schmidt metric in terms of the Ricci scalar of the Lorentzian manifold and give some examples. Finally, we discuss some applications to general relativity.


title: Exclusive neutrino-production of a charmed meson

abstract: We calculate the leading order in $\alpha_s$ QCD amplitude for exclusive neutrino and antineutrino production of a $D$ pseudoscalar charmed meson on an unpolarized nucleon. We work in the framework of the collinear QCD approach where generalized parton distributions (GPDs) factorize from perturbatively calculable coefficient functions. We include both $O(m_c)$ terms in the coefficient functions and $O(M_D)$ mass term contributions in the heavy meson distribution amplitudes. We emphasize the sensitivity of specific observables on the transversity quark GPDs.


title: On SGD's Failure in Practice: Characterizing and Overcoming Stalling

abstract: Stochastic Gradient Descent (SGD) is widely used in machine learning problems to efficiently perform empirical risk minimization, yet, in practice, SGD is known to stall before reaching the actual minimizer of the empirical risk. SGD stalling has often been attributed to its sensitivity to the conditioning of the problem; however, as we demonstrate, SGD will stall even when applied to a simple linear regression problem with unity condition number for standard learning rates. Thus, in this work, we numerically demonstrate and mathematically argue that stalling is a crippling and generic limitation of SGD and its variants in practice. Once we have established the problem of stalling, we generalize an existing framework for hedging against its effects, which (1) deters SGD and its variants from stalling, (2) still provides convergence guarantees, and (3) makes SGD and its variants more practical methods for minimization.


title: A Hybrid Evolutionary Algorithm Based on Solution Merging for the Longest Arc-Preserving Common Subsequence Problem

abstract: The longest arc-preserving common subsequence problem is an NP-hard combinatorial optimization problem from the field of computational biology. This problem finds applications, in particular, in the comparison of arc-annotated Ribonucleic acid (RNA) sequences. In this work we propose a simple, hybrid evolutionary algorithm to tackle this problem. The most important feature of this algorithm concerns a crossover operator based on solution merging. In solution merging, two or more solutions to the problem are merged, and an exact technique is used to find the best solution within this union. It is experimentally shown that the proposed algorithm outperforms a heuristic from the literature.


title: EPFL Lectures on General Relativity as a Quantum Field Theory

abstract: These notes are an introduction to General Relativity as a Quantum Effective Field Theory, following the material given in a short course on the subject at EPFL. The intent is to develop General Relativity starting from a quantum field theoretic viewpoint, and to introduce some of the techniques needed to understand the subject.


title: Finite-state Independence and Normal Sequences

abstract: We consider the previously defined notion of finite-state independence and we focus specifically on normal words. We characterize finite-state independence of normal words in three different ways, using three different kinds of asynchronous deterministic finite automata with two input tapes containing infinite words. Based on one of the characterizations we give an algorithm to construct a pair of finite-state independent normal words.


title: Optimality of integrability estimates for advection-diffusion equations

abstract: We discuss $L^p$ integrability estimates for the solution $u$ of the advection-diffusion equation $\partial_t u + \mathrm{div} (bu) = \Delta u$, where the velocity field $b \in L^r_t L^q_x$. We first summarize some classical results proving such estimates for certain ranges of the exponents $r$ and $q$. Afterwards we prove the optimality of such ranges by means of new original examples.


title: Counterintuitive Reconstruction of the Polar O-Terminated ZnO Surface With Zinc Vacancies and Hydrogen

abstract: Understanding the structure of ZnO surface reconstructions and their resultant properties is crucial to the rational design of ZnO-containing devices ranging from optoelectronics to catalysts. Here, we are motivated by recent experimental work which showed a new surface reconstruction containing Zn vacancies ordered in a Zn(3x3) pattern in the subsurface of (0001)-O terminated ZnO. A reconstruction with Zn vacancies on (0001)-O is surprising and counterintuitive because Zn vacancies enhance the surface dipole rather than reduce it. In this work, we show using Density Functional Theory (DFT) that subsurface Zn vacancies can form on (0001)-O when coupled with adsorption of surface H and are in fact stable under a wide range of common conditions. We also show these vacancies have a significant ordering tendency and that Sb-doping created subsurface inversion domain boundaries (IDBs) enhances the driving force of Zn vacancy alignment into large domains of the Zn(3x3) reconstruction.


title: Excitonic linewidth approaching the homogeneous limit in MoS2 based van der Waals heterostructures : accessing spin-valley dynamics

abstract: The strong light matter interaction and the valley selective optical selection rules make monolayer (ML) MoS2 an exciting 2D material for fundamental physics and optoelectronics applications. But so far optical transition linewidths even at low temperature are typically as large as a few tens of meV and contain homogenous and inhomogeneous contributions. This prevented in-depth studies, in contrast to the better-characterized ML materials MoSe2 and WSe2. In this work we show that encapsulation of ML MoS2 in hexagonal boron nitride can efficiently suppress the inhomogeneous contribution to the exciton linewidth, as we measure in photoluminescence and reflectivity a FWHM down to 2 meV at T = 4K. This indicates that surface protection and substrate flatness are key ingredients for obtaining stable, high quality samples. Among the new possibilities offered by the well-defined optical transitions we measure the homogeneous broadening induced by the interaction with phonons in temperature dependent experiments. We uncover new information on spin and valley physics and present the rotation of valley coherence in applied magnetic fields perpendicular to the ML.


title: The problem of characterising stationary data for the vacuum Einstein equations

abstract: We take a step towards characterising stationary data for the vacuum Einstein equations, by finding a necessary condition on initial data for which the evolution is a solution of the vacuum equations admitting a Killing vector, which is time-like at least in some region of the Cauchy development.


title: Hybrid Fuel Cells Power for Long Duration Robot Missions in Field Environments

abstract: Mobile robots are often needed for long duration missions. These include search and rescue, sentry, repair, surveillance and entertainment. Current power supply technology limit walking and climbing robots from many such missions. Internal combustion engines have high noise and emit toxic exhaust while rechargeable batteries have low energy densities and high rates of self-discharge. In theory, fuel cells do not have such limitations. In particular Proton Exchange Membrane (PEMs) can provide very high energy densities, are clean and quiet. However, PEM fuel cells are found to be unreliable due to performance degradation. This can be mitigated by protecting the fuel cell in a fuel-cell battery hybrid configuration using filtering electronics that ensure the fuel cell is isolated from electrical noise and a battery to isolate it from power surges. Simulation results are presented for a HOAP 2 humanoid robot that suggests a fuel cell powered hybrid power supply superior to conventional batteries.


title: Resolvability in c.c.c. generic extensions

abstract: Every crowded space $X$ is ${\omega}$-resolvable in the c.c.c generic extension $V^{Fn(|X|,2})$ of the ground model. We investigate what we can say about ${\lambda}$-resolvability in c.c.c-generic extensions for ${\lambda}>{\omega}$? A topological space is "monotonically $\omega_1$-resolvable" if there is a function $f:X\to {\omega_1}$ such that $$\{x\in X: f(x)\ge {\alpha} \}\subset^{dense}X $$ for each ${\alpha}<{\omega_1}$. We show that given a $T_1$ space $X$ the following statements are equivalent: (1) $X$ is ${\omega}_1$-resolvable in some c.c.c-generic extension, (2) $X$ is monotonically $\omega_1$-resolvable. (3) $X$ is ${\omega}_1$-resolvable in the Cohen-generic extension $V^{Fn({\omega_1},2)}$. We investigate which spaces are monotonically $\omega_1$-resolvable. We show that if a topological space $X$ is c.c.c, and ${\omega}_1\le \Delta(X)\le |X|<{\omega}_{\omega}$, then $X$ is monotonically $\omega_1$-resolvable. On the other hand, it is also consistent, modulo the existence of a measurable cardinal, that there is a space $Y$ with $|Y|=\Delta(Y)=\aleph_\omega$ which is not monotonically $\omega_1$-resolvable. The characterization of ${\omega_1}$-resolvability in c.c.c generic extension raises the following question: is it true that crowded spaces from the ground model are ${\omega}$-resolvable in $V^{Fn({\omega},2)}$? We show that (i) if $V=L$ then every crowded c.c.c. space $X$ is ${\omega}$-resolvable in $V^{Fn({\omega},2)}$, (ii) if there is no weakly inaccesssible cardinals, then every crowded space $X$ is ${\omega}$-resolvable in $V^{Fn({\omega}_1,2)}$. On the other hand, it is also consistent that there is a crowded space $X$ with $|X|=\Delta(X)={\omega_1}$ such that $X$ remains irresolvable after adding a single Cohen real.


title: Variable dispersion beta regressions with parametric link functions

abstract: This paper presents a new class of regression models for continuous data restricted to the interval $(0,1)$, such as rates and proportions. The proposed class of models assumes a beta distribution for the variable of interest with regression structures for the mean and dispersion parameters. These structures consider covariates, unknown regression parameters, and parametric link functions. Link functions depend on parameters that model the relationship between the random component and the linear predictors. The symmetric and assymetric Aranda-Ordaz link functions are considered in details. Depending on the parameter values, these link functions refer to particular cases of fixed links such as logit and complementary log-log functions. Joint estimation of the regression and link function parameters is performed by maximum likelihood. Closed-form expressions for the score function and Fisher's information matrix are presented. Aspects of large sample inferences are discussed, and some diagnostic measures are proposed. A Monte Carlo simulation study is used to evaluate the finite sample performance of point estimators. Finally, a practical application that employs real data is presented and discussed.


title: Robust iterative schemes for non-linear poromechanics

abstract: We consider a non-linear extension of Biot's model for poromechanics, wherein both the fluid flow and mechanical deformation are allowed to be non-linear. We perform an implicit discretization in time (backward Euler) and propose two iterative schemes for solving the non-linear problems appearing within each time step: a splitting algorithm extending the undrained split and fixed stress methods to non-linear problems, and a monolithic L-scheme. The convergence of both schemes is shown rigorously. Illustrative numerical examples are presented to confirm the applicability of the schemes and validate the theoretical results.


title: Epigenome-wide association study and integrative analysis with the transcriptome based on GWAS summary statistics

abstract: The past decade has seen a rapid growth in omics technologies. Genome-wide association studies (GWAS) have uncovered susceptibility variants for a variety of complex traits. However, the functional significance of most discovered variants are still not fully understood. On the other hand, there is increasing interest in exploring the role of epigenetic variations such as DNA methylation in disease pathogenesis. In this work, we present a general framework for epigenome-wide association study and integrative analysis with the transcriptome based on GWAS summary statistics and data from methylation and expression quantitative trait loci (QTL) studies. The framework is based on Mendelian randomization, which is much less vulnerable to confounding and reverse causation compared to conventional studies. The framework was applied to five complex diseases. We first identified loci that are differentially methylated due to genetic variations, and then developed several approaches for joint testing with the GWAS-imputed transcriptome. We discovered a number of novel candidate genes that are not implicated in the original GWAS studies. We also observed strong evidence (lowest p = 2.01e-184) for differential expression among the top genes mapped to methylation loci. The framework proposed here opens a new way of analyzing GWAS summary data and will be useful for gaining deeper insight into disease mechanisms.


title: Continuum limit and universality of the Columbia plot

abstract: Results on the thermal transition of QCD with 3 degenerate flavors, in the lower-left corner of the Columbia plot, are puzzling. The transition is expected to be first-order for massless quarks, and to remain so for a range of quark masses until it turns second-order at a critical quark mass. But this critical quark mass and resulting "pion" mass disagree violently between Wilson and staggered fermions at finite lattice spacing, and decrease sharply with the lattice spacing, for staggered fermions at least. To clarify this puzzle and eliminate potential systematic effects from rooting, we study the 4-flavor theory with staggered fermions, on lattices with 4 to 10 time-slices. Our results are qualitatively similar to the 3-flavor case, so that rooting is not an issue. However, dramatic cutoff effects are visible, even on our finest lattices. Universality implies that cutoff effects for Wilson fermions are even more dramatic. In order to obtain a first-order thermal transition in the continuum theory, extremely light quarks are needed.


title: Room-temperature 1.54 $\mu$m photoluminescence of Er:O$_x$ centers at extremely low concentration in silicon

abstract: The demand for single photon sources at $\lambda~=~1.54~\mu$m, which follows from the consistent development of quantum networks based on commercial optical fibers, makes Er:O$_x$ centers in Si still a viable resource thanks to the optical transition of $Er^{3+}~:~^4I_{13/2}~\rightarrow~^4I_{15/2}$. Yet, to date, the implementation of such system remains hindered by its extremely low emission rate. In this Letter, we explore the room-temperature photoluminescence (PL) at the telecomm wavelength of very low implantation doses of $Er:O_x$ in $Si$. The emitted photons, excited by a $\lambda~=~792~nm$ laser in both large areas and confined dots of diameter down to $5~\mu$m, are collected by an inverted confocal microscope. The lower-bound number of detectable emission centers within our diffraction-limited illumination spot is estimated to be down to about 10$^4$, corresponding to an emission rate per individual ion of about $4~\times~10^{3}$ photons/s.


title: Intrusion and extrusion of a liquid on nanostructured surfaces

abstract: Superhydrophobicity is connected to the presence of gas pockets within surface asperities. Upon increasing the pressure this "suspended" state may collapse, causing the complete wetting of the rough surface. In order to quantitatively characterize this process on nanostructured surfaces, we perform rare-event atomistic simulations at different pressures and for several texture geometries. Such approach allows us to identify for each pressure the stable and metastable states and the free energy barriers separating them. Results show that, by starting from the superhydrophobic state and increasing the pressure, the suspended state abruptly collapses at a critical intrusion pressure. If the pressure is subsequently decreased, the system remains trapped in the metastable state corresponding to the wet surface. The liquid can be extruded from the nanostructures only at very negative pressures, by reaching the critical extrusion pressure (spinodal for the confined liquid). The intrusion and extrusion curves form a hysteresis cycle determined by the large free energy barriers separating the suspended and wet state. These barriers, which grow very quickly for pressures departing from the intrusion/extrusion pressure, are shown to strongly depend on the texture geometry.


title: All entangled pure quantum states violate the bilocality inequality

abstract: The nature of quantum correlations in networks featuring independent sources of entanglement remains poorly understood. Here, focusing on the simplest network of entanglement swapping, we start a systematic characterization of the set of quantum states leading to violation of the so-called "bilocality" inequality. First, we show that all possible pairs of entangled pure states can violate the inequality. Next, we derive a general criterion for violation for arbitrary pairs of mixed two-qubit states. Notably, this reveals a strong connection between the CHSH Bell inequality and the bilocality inequality, namely that any entangled state violating CHSH also violates the bilocality inequality. We conclude with a list of open questions.


title: Anisotropic hydrodynamic turbulence in accretion disks

abstract: Recently, the vertical shear instability (VSI) has become an attractive purely hydrodynamic candidate for the anomalous angular momentum transport required for weakly ionized accretion disks. In direct three-dimensional numerical simulations of VSI turbulence in disks, a meridional circulation pattern was observed that is opposite to the usual viscous flow behavior. Here, we investigate whether this feature can possibly be explained by an anisotropy of the VSI turbulence. Using three-dimensional hydrodynamical simulations, we calculate the turbulent Reynolds stresses relevant for angular momentum transport for a representative section of a disk. We find that the vertical stress is significantly stronger than the radial stress. Using our results in viscous disk simulations with different viscosity coefficients for the radial and vertical direction, we find good agreement with the VSI turbulence for the stresses and meridional flow; this provides additional evidence for the anisotropy. The results are important with respect to the transport of small embedded particles in disks.


title: Control of a Bucket-Wheel for Surface Mining of Asteroids and Small-Bodies

abstract: Near Earth Asteroids (NEAs) are thought to contain a wealth of resources, including water, iron, titanium, nickel, platinum and silicates. Future space missions that can exploit these resources by performing In-Situ Resource Utilization (ISRU) gain substantial benefit in terms of range, payload capacity and mission flexibility. Compared to the Moon or Mars, the milligravity on some asteroids demands a fraction of the energy for digging and accessing hydrated regolith just below the surface. However, asteroids and small-bodies, because of their low gravity present a major challenge in landing, surface excavation and resource capture. These challenges have resulted in adoption of a "touch and go techniques", like the upcoming Osiris-rex sample-return mission. Previous asteroid excavation efforts have focused on discrete capture events (an extension of sampling technology) or whole-asteroid capture and processing. This paper analyzes the control of a bucket-wheel design for asteroid or small-body excavation. Our study focuses on system design of two counter rotating bucket-wheels that are attached to a hovering spacecraft. Regolith is excavated and heated to 1000 C to extract water. The water in turn is electrolyzed to produce hydrogen and oxygen for rocket fuel. We analyze control techniques to maximize traction of the bucket-wheels on the asteroid surface and minimize lift-off the surface, together with methods to dig deeper into the asteroid surface. Our studies combine analytical models, with simulation and hardware testing. For initial evaluation of material-spacecraft dynamics and mechanics, we assume lunar-like regolith for bulk density, particle size and cohesion. Our early studies point towards a promising pathway towards refinement of this technology for demonstration aboard a future space mission.


title: Aspects of Cubical Higher Category Theory

abstract: In this article we show how to build main aspects of our paper on globular weak $(\infty,n)$-categories, but now for the cubical geometry. Thus we define a monad on the category $\mathbb{C}\mathbb{S}ets$ of cubical sets which algebras are models of cubical weak $\infty$-categories. Also for each $n\in\mathbb{N}$ we define a monad on $\mathbb{C}\mathbb{S}ets$ which algebras are models of cubical weak $(\infty,n)$-categories. And finally we define a monad on the category $\mathbb{C}\mathbb{S}ets^2$ which algebras are models of cubical weak $\infty$-functors, and a monad on the category $\mathbb{C}\mathbb{S}ets^4$ which algebras are models of cubical weak natural $\infty$-transformations.


title: Effective equilibrium states in the colored-noise model for active matter II. A unified framework for phase equilibria, structure and mechanical properties

abstract: Active particles driven by colored noise can be approximately mapped onto a system that obeys detailed balance. The effective interactions which can be derived for such a system allow to describe the structure and phase behavior of the active fluid by means of an effective free energy. In this paper we explain why the related thermodynamic results for pressure and interfacial tension do not represent the results one would measure mechanically. We derive a dynamical density functional theory, which in the steady state simultaneously validates the use of effective interactions and provides access to mechanical quantities. Our calculations suggest that in the colored-noise model the mechanical pressure in coexisting phases might be unequal and the interfacial tension can become negative.


title: Siamese Network of Deep Fisher-Vector Descriptors for Image Retrieval

abstract: This paper addresses the problem of large scale image retrieval, with the aim of accurately ranking the similarity of a large number of images to a given query image. To achieve this, we propose a novel Siamese network. This network consists of two computational strands, each comprising of a CNN component followed by a Fisher vector component. The CNN component produces dense, deep convolutional descriptors that are then aggregated by the Fisher Vector method. Crucially, we propose to simultaneously learn both the CNN filter weights and Fisher Vector model parameters. This allows us to account for the evolving distribution of deep descriptors over the course of the learning process. We show that the proposed approach gives significant improvements over the state-of-the-art methods on the Oxford and Paris image retrieval datasets. Additionally, we provide a baseline performance measure for both these datasets with the inclusion of 1 million distractors.


title: Block circulant and Toeplitz structures in the linearized Hartree-Fock equation on finite lattices: tensor approach

abstract: This paper introduces and analyses the new grid-based tensor approach to approximate solution of the elliptic eigenvalue problem for the 3D lattice-structured systems. We consider the linearized Hartree-Fock equation over a spatial $L_1\times L_2\times L_3$ lattice for both periodic and non-periodic problem setting, discretized in the localized Gaussian-type orbitals basis. In the periodic case, the Galerkin system matrix obeys a three-level block-circulant structure that allows the FFT-based diagonalization, while for the finite extended systems in a box (Dirichlet boundary conditions) we arrive at the perturbed block-Toeplitz representation providing fast matrix-vector multiplication and low storage size. The proposed grid-based tensor techniques manifest the twofold benefits: (a) the entries of the Fock matrix are computed by 1D operations using low-rank tensors represented on a 3D grid, (b) in the periodic case the low-rank tensor structure in the diagonal blocks of the Fock matrix in the Fourier space reduces the conventional 3D FFT to the product of 1D FFTs. Lattice type systems in a box with Dirichlet boundary conditions are treated numerically by our previous tensor solver for single molecules, which makes possible calculations on rather large $L_1\times L_2\times L_3$ lattices due to reduced numerical cost for 3D problems. The numerical simulations for both box-type and periodic $L\times 1\times 1$ lattice chain in a 3D rectangular "tube" with $L$ up to several hundred confirm the theoretical complexity bounds for the block-structured eigenvalue solvers in the limit of large $L$.


title: BFR: a Bloom Filter-based Routing Approach for Information-Centric Networks

abstract: Locating the demanded content is one of the major challenges in Information-Centric Networking (ICN). This process is known as content discovery. To facilitate content discovery, in this paper we focus on Named Data Networking (NDN) and propose a novel routing scheme for content discovery, called Bloom Filter-based Routing (BFR), which is fully distributed, content oriented, and topology agnostic at the intra-domain level. In BFR, origin servers advertise their content objects using Bloom filters. We compare the performance of the proposed BFR with flooding and shortest path content discovery approaches. BFR outperforms its counterparts in terms of the average round-trip delay, while it is shown to be very robust to false positive reports from Bloom filters. Also, BFR is much more robust than shortest path routing to topology changes. BFR strongly outperforms flooding and performs almost equal with shortest path routing with respect to the normalized communication costs for data retrieval and total communication overhead for forwarding Interests. All the three approaches achieve similar mean hit distance. The signalling overhead for content advertisement in BFR is much lower than the signalling overhead for calculating shortest paths in the shortest path approach. Finally, BFR requires small storage overhead for maintaining content advertisements.


title: Anisotropic blockade using pendular Rydberg butterfly molecules

abstract: We propose an experiment to demonstrate a novel blockade mechanism caused by long-range anisotropic interactions in an ultracold dipolar gas composed of the recently observed "butterfly" Rydberg molecules. At the blockade radius, the strong intermolecular interaction between two adjacent molecules shifts their molecular states out of resonance with the photoassociation laser, preventing their simultaneous excitation. When the molecules are prepared in a quasi-one-dimensional (Q1D) trap, the interaction's strength can be tuned via a weak external field. The molecular density thus depends strongly on the angle between the trap axis and the field. The available Rydberg and molecular states provide a wide range of tunability.


title: A note on local rigidity

abstract: The aim of this note is to give a geometric proof for classical local rigidity of lattices in semisimple Lie groups. We are reproving well known results in a more geometric (and hopefully clearer) way.


title: The inversion of motion of bristle bots: analytical and experimental analysis

abstract: Bristle bots are vibration-driven robots actuated by the motion of an internal oscillating mass. Vibrations are translated into directed locomotion due to the alternating friction resistance between robots' bristles and the substrate during oscillations. Bristle bots are, in general, unidirectional locomotion systems. In this paper we demonstrate that motion direction of vertically vibrated bristle systems can be controlled by tuning the frequency of their oscillatory actuation. We report theoretical and experimental results obtained by studying an equivalent system, consisting of an inactive robot placed on a vertically vibrating substrate.


title: Parametric representations and boundary fixed points of univalent self-maps of the unit disk

abstract: A classical result in the theory of Loewner's parametric representation states that the semigroup $\mathfrak U_*$ of all conformal self-maps $\phi$ of the unit disk $\mathbb{D}$ normalized by $\phi(0) = 0$ and $\phi'(0) > 0$ can be obtained as the reachable set of the Loewner - Kufarev control system $$ \frac{\mathrm{d} w_t}{\mathrm{d} t}=G_t\circ w_t,\quad t\geqslant0,\qquad w_0=\mathsf{id}_{\mathbb{D}}, $$ where the control functions $t\mapsto G_t\in\mathsf{Hol}(\mathbb{D},\mathbb{C})$ form a certain convex cone. Here we extend this result to the semigroup $\mathfrak U[F]$ consisting of all conformal $\phi:\mathbb{D}\to\mathbb{D}$ whose set of boundary regular fixed points contains a given finite set $F\subset\partial\mathbb{D}$ and to its subsemigroup $\mathfrak U_\tau[F]$ formed by $\mathsf{id}_{\mathbb{D}}$ and all $\phi\in\mathfrak U[F]\setminus\{\mathsf{id}_{\mathbb{D}}\}$ with the prescribed boundary Denjoy - Wolff point $\tau\in\partial\mathbb{D}\setminus F$. This completes the study launched in [P. Gumenyuk, Preprint 2016, ArXiv:1603.04043], where the case of interior Denjoy - Wolff point $\tau\in\mathbb{D}$ was considered.


title: Ratios of helicity amplitudes for exclusive $\rho^0$ electroproduction on transversely polarized protons

abstract: Exclusive $\rho^0$-meson electroproduction is studied by the HERMES experiment, using the 27.6 GeV longitudinally polarized electron/positron beam of HERA and a transversely polarized hydrogen target, in the kinematic region 1.0 GeV$^2$<Q$^2$<7.0 GeV$^2$, 3.0 GeV<W<6.3 GeV, and -t'<0.4 GeV$^2$. Using an unbinned maximum-likelihood method, 25 parameters are extracted. These determine the real and imaginary parts of the ratios of several helicity amplitudes describing $\rho^0$-meson production by a virtual photon. The denominator of those ratios is the dominant amplitude, the nucleon-helicity-non-flip amplitude $F_{0\frac{1}{2}0\frac{1}{2}}$, which describes the production of a longitudinal $\rho^{0}$-meson by a longitudinal virtual photon. The ratios of nucleon-helicity-non-flip amplitudes are found to be in good agreement with those from the previous HERMES analysis. The transverse target polarization allows for the first time the extraction of ratios of a number of nucleon-helicity-flip amplitudes to $F_{0\frac{1}{2}0\frac{1}{2}}$. Results obtained in a handbag approach based on generalized parton distributions taking into account the contribution from pion exchange are found to be in good agreement with these ratios. Within the model, the data favor a positive sign for the $\pi-\rho$ transition form factor. By also exploiting the longitudinal beam polarization, a total of 71 $\rho^0$ spin-density matrix elements is determined from the extracted 25 parameters, in contrast to only 53 elements as directly determined in earlier analyses.


title: Quasi-Spectral Sparse Bi-Global Stability Analysis of Compressible Channel Flow over Complex Impedance

abstract: We have developed a fully sparse, compact-scheme based biglobal stability analysis numerical solver applied, for the scope of the current paper, to the investigation of the effects of impedance boundary conditions (IBCs) on the structure of a fully developed compressible turbulent channel flow. A sixth-order compact finite difference scheme is used to discretize the linearized Navier-Stokes equations leading to a Generalized Eigenvalue Problem (GEVP). Sparsity is retained by explicitly introducing derivatives of the perturbation as additional unknowns, increasing the overall problem size (number of columns $\times$ number of rows) while significantly reducing the number of non-zeros and the computational cost with respect to traditional implementations yielding otherwise dense matrix blocks. The resulting GEVP is coded in Python and solved employing an Message Passing Interface (MPI) parallelized PETSc-based sparse eigenvalue solver adopting a modified Arnoldi algorithm. Base flow is taken from impermeable isothermal-wall turbulent channel flow simulations at bulk Reynolds number, $Re_b= 6900$ and Mach number, $M_b = 0.85$. The eigenvalue spectrum calculated in the biglobal stability analysis shows distinct groups of modes associated with a discrete set of streamwise wave numbers accommodated by computational domain. An iterative strategy for the imposition of the complex IBCs, which are a nonlinear function of the real-valued (Fourier) frequency in the GEVP, has been devised. The adopted IBC specifically represents an array of sub-surface-mounted Helmholtz cavities with resonant frequency, $f_{res}$, covered by a porous sheet with permeability inversely proportional to the impedance resistance $R$. The tunable resonant frequency has been shown to be an attractor for the instability, yielding a single unstable mode at that frequency.


title: Optimal weak value measurements: Pure states

abstract: We apply the notion of \emph{optimality} of measurements for state determination(tomography) as originally given by Wootters and Fields to \emph{weak value tomography} of \emph{pure states}. They defined measurements to be optimal if they 'minimised' the effects of statistical errors. For technical reasons they actually maximised the state averaged information, precisely quantified as the negative logarithm of 'error volume'. In this paper we optimise both the state averaged information as well as error volumes. We prove, for Hilbert spaces of arbitrary (finite) dimensionality, that varieties of weak value measurements are optimal when the post-selected bases are \emph{mutually unbiased} with respect to the eigenvectors of the observable being measured. We prove a number of important results about the geometry of state spaces when expressed through the weak values as coordinates. We derive an expression for the Ka\"ehler potential for the N-dimensional case with the help of which we give an exact treatment of the arbitrary-spin case.


title: Two-body decays of gluino at full one-loop level in the quark-flavour violating MSSM

abstract: We study the two-body decays of the gluino at full one-loop level in the Minimal Supersymmetric Standard Model with quark-flavour violation (QFV) in the squark sector. The renormalization is done in the $\overline{\rm DR}$ scheme. The gluon and photon radiations are included by adding the corresponding three-body decay widths. We discuss the dependence of the gluino decay widths on the QFV parameters. The main dependence stems from the $\tilde c_R - \tilde t_R$ mixing in the decays to up-type squarks, and from the $\tilde s_R - \tilde b_R$ mixing in the decays to down-type squarks due to the strong constraints from B-physics on the other quark-flavour mixing parameters. The full one-loop corrections to the gluino decay widths are mostly negative and of the order of about -10\%. The QFV part stays small in the total width but can vary up to -8\% for the decay width into the lightest squark. For the corresponding branching ratio the effect is somehow washed out by at least a factor of two. The electroweak corrections can be as large as 35\% of the SUSY QCD corrections.


title: Entangling two atoms of different isotopes via Rydberg blockade

abstract: Quantum entanglement is crucial for simulating and understanding exotic physics of strongly correlated many-body systems, such as high--temperature superconductors, or fractional quantum Hall states. The entanglement of non-identical particles exhibits richer physics of strong many-body correlations and offers more opportunities for quantum computation, especially with neutral atoms where in contrast to ions the interparticle interaction is widely tunable by Feshbach resonances. Moreover, the inter-species entanglement forms a basis for the properties of various compound systems, ranging from Bose-Bose mixtures to photosynthetic light-harvesting complexes. So far, the inter-species entanglement has only been obtained for trapped ions. Here we report on the experimental realization of entanglement of two neutral atoms of different isotopes. A ${}^{87}\mathrm{Rb}$ atom and a ${}^{85}\mathrm{Rb}$ atom are confined in two single--atom optical traps separated by 3.8 $\mu$m. Creating a strong Rydberg blockade, we demonstrate a heteronuclear controlled--NOT (C--NOT) quantum gate and generate a heteronuclear entangled state, with raw fidelities $0.73 \pm 0.01$ and $0.59 \pm 0.03$, respectively. Our work, together with the technologies of single--qubit gate and C--NOT gate developed for identical atoms, can be used for simulating any many--body system with multi-species interactions. It also has applications in quantum computing and quantum metrology, since heteronuclear systems exhibit advantages in low crosstalk and in memory protection.


title: Effect of iron oxide loading on magnetoferritin structure in solution as revealed by SAXS and SANS

abstract: Synthetic biological macromolecule of magnetoferritin containing an iron oxide core inside a protein shell (apoferritin) is prepared with different content of iron. Its structure in aqueous solution is analyzed by small-angle synchrotron X-ray (SAXS) and neutron (SANS) scattering. The loading factor (LF) defined as the average number of iron atoms per protein is varied up to LF=800. With an increase of the LF, the scattering curves exhibit a relative increase in the total scattered intensity, a partial smearing and a shift of the match point in the SANS contrast variation data. The analysis shows an increase in the polydispersity of the proteins and a corresponding effective increase in the relative content of magnetic material against the protein moiety of the shell with the LF growth. At LFs above ~150, the apoferritin shell undergoes structural changes, which is strongly indicative of the fact that the shell stability is affected by iron oxide presence.


title: Distribution of phenotype sizes in sequence-to-structure genotype-phenotype maps

abstract: An essential quantity to ensure evolvability of populations is the navigability of the genotype space. Navigability relies on the existence of sufficiently large genotype networks, that is ensembles of sequences with the same phenotype that guarantee an efficient random drift through sequence space. The number of sequences compatible with a given structure (e.g. the number of RNA sequences folding into a particular secondary structure, or the number of DNA sequences coding for the same protein structure) is astronomically large in all functional molecules investigated. However, an exhaustive experimental or computational study of all RNA folds or all protein structures becomes impossible even for moderately long sequences. Here, we analytically derive the distribution of phenotype sizes for a hierarchy of models which successively incorporate features of increasingly realistic sequence-to-structure genotype-phenotype maps. The main feature of these models relies on the characterization of each phenotype through a prototypical sequence whose sites admit a variable fraction of letters of the alphabet. Our models interpolate between two limit distributions: a powerlaw distribution, when the ordering of sites in the prototypical sequence is strongly constrained, and a lognormal distribution, as suggested for RNA, when different orderings of the same set of sites yield different phenotypes. Our main result is the qualitative and quantitative identification of those features of the sequence-to-structure map that lead to different distributions of phenotype sizes.


title: Towards Radiative Transitions in Charmonium

abstract: We present preliminary calculations towards radiative transitions in charmonium using anisotropic $N_f = 2 + 1$ dynamical ensembles generated by the Hadron Spectrum Collaboration. With the use of newer technologies we aim to investigate transitions between states, including potential exotic charmonium states, lying higher in the spectrum than in previous studies. A crucial ingredient in this work is the use of variationally optimised interpolating operators which allow for a reliable determination of the three-point correlation functions needed. Using these operators, we perform first calculations of relevant three-point correlation functions before discussing future directions.


title: The non-cooperative tile assembly model is not intrinsically universal or capable of bounded Turing machine simulation

abstract: The field of algorithmic self-assembly is concerned with the computational and expressive power of nanoscale self-assembling molecular systems. In the well-studied cooperative, or temperature 2, abstract tile assembly model it is known that there is a tile set to simulate any Turing machine and an intrinsically universal tile set that simulates the shapes and dynamics of any instance of the model, up to spatial rescaling. It has been an open question as to whether the seemingly simpler noncooperative, or temperature 1, model is capable of such behaviour. Here we show that this is not the case, by showing that there is no tile set in the noncooperative model that is intrinsically universal, nor one capable of time-bounded Turing machine simulation within a bounded region of the plane. Although the noncooperative model intuitively seems to lack the complexity and power of the cooperative model it was not obvious how to prove this. One reason is that there have been few tools to analyse the structure of complicated paths in the plane. This paper provides a number of such tools. A second reason is that almost every obvious and small generalisation to the model (e.g. allowing error, 3D, non-square tiles, signals/wires on tiles, tiles that repel each other, parallel synchronous growth) endows it with great computational, and sometimes simulation, power. Our main results show that all of these generalisations provably increase computational and/or simulation power. Our results hold for both deterministic and nondeterministic noncooperative systems. Our first main result stands in stark contrast with the fact that for both the cooperative tile assembly model, and for 3D noncooperative tile assembly, there are respective intrinsically universal tilesets. Our second main result gives a new technique (reduction to simulation) for proving negative results about computation in tile assembly.


title: Topological Principles of Control in Dynamical Network Systems

abstract: Networked systems display complex patterns of interactions between a large number of components. In physical networks, these interactions often occur along structural connections that link components in a hard-wired connection topology, supporting a variety of system-wide dynamical behaviors such as synchronization. While descriptions of these behaviors are important, they are only a first step towards understanding the relationship between network topology and system behavior, and harnessing that relationship to optimally control the system's function. Here, we use linear network control theory to analytically relate the topology of a subset of structural connections (those linking driver nodes to non-driver nodes) to the minimum energy required to control networked systems. As opposed to the numerical computations of control energy, our accurate closed-form expressions yield general structural features in networks that require significantly more or less energy to control, providing topological principles for the design and modification of network behavior. To illustrate the utility of the mathematics, we apply this approach to high-resolution connectomes recently reconstructed from drosophila, mouse, and human brains. We use these principles to show that connectomes of increasingly complex species are wired to reduce control energy. We then use the analytical expressions we derive to perform targeted manipulation of the brain's control profile by removing single edges in the network, a manipulation that is accessible to current clinical techniques in patients with neurological disorders. Cross-species comparisons suggest an advantage of the human brain in supporting diverse network dynamics with small energetic costs, while remaining unexpectedly robust to perturbations. Our results ground the expectation of a system's dynamical behavior in its network architecture.