title: Anomalous diffusion in convergence to effective ergodicity

abstract: Power-law exponents in the convergence to effective ergodicity is quantified for Ising-Lenz model in one dimension. Modified Thirumalai-Mountain (TM) metric for magnetisation is computed for the range of temperature values under strongly correlated dynamics. In producing evolution of TM metric over time, time-averaged dynamics is generated by using Metropolis and Glauber single-spin-flip dynamics, and ensemble-averaged dynamics with an exact solution. Superdiffusive behaviour is numerically identified in the parameter regimes studied, i.e., power-law exponents, $\alpha > 1.0$.


title: Scalable image coding based on epitomes

abstract: In this paper, we propose a novel scheme for scalable image coding based on the concept of epitome. An epitome can be seen as a factorized representation of an image. Focusing on spatial scalability, the enhancement layer of the proposed scheme contains only the epitome of the input image. The pixels of the enhancement layer not contained in the epitome are then restored using two approaches inspired from local learning-based super-resolution methods. In the first method, a locally linear embedding model is learned on base layer patches and then applied to the corresponding epitome patches to reconstruct the enhancement layer. The second approach learns linear mappings between pairs of co-located base layer and epitome patches. Experiments have shown that significant improvement of the rate-distortion performances can be achieved compared to an SHVC reference.


title: On Hadamard products of linear varieties

abstract: In this paper we address the Hadamard product of linear varieties not necessarily in general position. In $\mathbb{P}^2$ we obtain a complete description of the possible outcomes. In particular, in the case of two disjoint finite sets X and X' of collinear points, we get conditions for their hadamard product to be either a collinear finite set of points or a grid of | X|| X'| points. In $\mathbb{P}^3,$ under suitable conditions (which we prove to be generic), we show that the hadamard product of X and X' consists of | X||X'| points on the two different rulings of a non-degenerate quadric and we compute its Hilbert function in the case |X| = |X'|.


title: Energetic electron precipitation into the middle atmosphere - Constructing the loss cone fluxes from MEPED POES

abstract: The impact of energetic electron precipitation (EEP) on the chemistry of the middle atmosphere (50-90 km) is still an outstanding question as accurate quantification of EEP is lacking due to instrumental challenges and insufficient pitch angle coverage of current particle detectors. The Medium Energy Proton and Electron Detectors (MEPED) instrument on board the NOAA/Polar Orbiting Environmental Satellites(POES) and MetOp spacecraft has two sets of electron and proton telescopes pointing close to zenith ($0 ^{\circ}$) and in the horizontal plane ($90 ^{\circ}$). Using measurements from either the $0 ^{\circ}$ or $90 ^{\circ}$ telescope will underestimate or overestimate the bounce loss cone flux, respectively, as the energetic electron fluxes are often strongly anisotropic with decreasing fluxes toward the center of the loss cone. By combining the measurements from both telescopes with electron pitch angle distributions from theory of wave-particle interactions in the magnetosphere, a complete bounce loss cone flux is constructed for each of the electron energy channels >50 keV, >100 keV, and >300 keV. We apply a correction method to remove proton contamination in the electron counts. We also account for the relativistic (>1000 keV) electrons contaminating the proton detector at subauroral latitudes. This gives us full range coverage of electron energies that will be deposited in the middle atmosphere. Finally, we demonstrate the method's applicability on strongly anisotropic pitch angle distributions during a weak geomagnetic storm in February 2008. We compare the electron fluxes and subsequent energy deposition estimates to OH observations from the Microwave Limb Sounder on the Aura satellite substantiating that the estimated fluxes are representative for the true precipitating fluxes impacting the atmosphere.


title: Spatially and spectrally resolved quantum path interference with chirped driving pulses

abstract: We measure spectrally and spatially resolved high-order harmonics generated in argon using chirped multi-cycle laser pulses. Using a stable, high-repetition rate laser we observe detailed interference structures in the far-field. The structures are of two kinds; off-axis interference from the long trajectory only and on-axis interference including the short and long trajectories. The former is readily visible in the far-field spectrum, modulating both the spectral and spatial profile. To access the latter, we vary the chirp of the fundamental, imparting different phases on the different trajectories, thereby changing their relative phase. Using this method together with an analytical model, we are able to explain the on-axis behaviour and access the dipole phase parameters for the short (\(\alpha_s\)) and long (\(\alpha_l\)) trajectories. The extracted results compare very well with phase parameters calculated by solving the time-dependent Schr\"odinger equation. Going beyond the analytical model, we are also able to successfully reproduce the off-axis interference structure.


title: Reviving Threshold-Moving: a Simple Plug-in Bagging Ensemble for Binary and Multiclass Imbalanced Data

abstract: Class imbalance presents a major hurdle in the application of data mining methods. A common practice to deal with it is to create ensembles of classifiers that learn from resampled balanced data. For example, bagged decision trees combined with random undersampling (RUS) or the synthetic minority oversampling technique (SMOTE). However, most of the resampling methods entail asymmetric changes to the examples of different classes, which in turn can introduce its own biases in the model. Furthermore, those methods require a performance measure to be specified a priori before learning. An alternative is to use a so-called threshold-moving method that a posteriori changes the decision threshold of a model to counteract the imbalance, thus has a potential to adapt to the performance measure of interest. Surprisingly, little attention has been paid to the potential of combining bagging ensemble with threshold-moving. In this paper, we present probability thresholding bagging (PT-bagging), a versatile plug-in method that fills this gap. Contrary to usual rebalancing practice, our method preserves the natural class distribution of the data resulting in well calibrated posterior probabilities. We also extend the proposed method to handle multiclass data. The method is validated on binary and multiclass benchmark data sets. We perform analyses that provide insights into the proposed method.


title: How to Compute Halting

abstract: A consistently specified halting function may be computed.


title: First detection at 5.5 and 9 GHz of the radio relics in bullet cluster with ATCA

abstract: We present here results from observations at 5.5 and 9 GHz of the Bullet cluster 1E 0657-55.8 with the Australia Telescope Compact Array (ATCA). Our results show detection of diffuse emission in the cluster. Our findings are consistent with the previous observations by Shimwell et al. (2014) and Shimwell et al. (2015) at 1.1-3.1 GHz. Morphology of diffuse structures (relic regions A and B and the radio halo) are consistent with those reported by the previous study. Our results indicate steepening in the spectral index at higher frequencies (at and greater than 5.0 GHz) for region A. The spectrum can be fit well by a broken power law. We discuss the possibility of a few recent theoretical models explaining this break in the power law spectrum, and find that a modified Diffusive Shock Acceleration (DSA) model or a turbulent reacceleration model may be relevant. Deep radio observations at high frequencies (at and greater than 5 GHz) are required for a detailed comparison with this model.


title: The Hydrogen Epoch of Reionization Array (HERA). Improvement of the antenna response with a matching network and scientific impacts

abstract: The Hydrogen Epoch of Reionization Array (HERA) is a new powerful radio-telescope, dedicated to the study of the early universe. Its main goal is to characterise the period of the universe where the first galaxies and stars started to form, by studying the evolution of the 21-cm emission signal from neutral hydrogen, during the Epoch of Reionization. In this article, we present an electromagnetic and electrical co-simulation of the antenna performed with CST. We focus our analysis on the characterisation of the chromatic effects caused by the antenna, in particular multiple reflections of the received signal between the feed and the dish. These reflections can have an important impact on the scientific results, and it is crucial to keep them as low as possible. Therefore, we are currently developing a matching circuit which aims to improve the impedance matching between the feed and the front end.


title: Nonlinear Luttinger liquid: Exact result for the Green function in terms of the fourth Painlev\'e transcendent

abstract: We show that exact time dependent single particle Green function in the Imambekov-Glazman theory of nonlinear Luttinger liquids can be written, for any value of the Luttinger parameter, in terms of a particular solution of the Painlev\'e IV equation. Our expression for the Green function has a form analogous to the celebrated Tracy-Widom result connecting the Airy kernel with Painlev\'e II. The asymptotic power law of the exact solution as a function of a single scaling variable $x/\sqrt{t}$ agrees with the mobile impurity results. The full shape of the Green function in the thermodynamic limit is recovered with arbitrary precision via a simple numerical integration of a nonlinear ODE.


title: Observations on the Halting Problem

abstract: The Halting Problem is ill-conceived and ill-defined.


title: The nanophysics age and its new perspectives

abstract: The nanophysics is halfway between the size scales of quantum mechanics and macroscopic physics governed by the laws of Newton and Einstein. The correct definition of nanophysics is the physics of structures and artefacts with dimensions in the nanometer range or of phenomena occurring in nanoseconds. Modern physical methods whose fundamental are developed in physics laboratories have become critically important in nanoscience. Nanophysics brings together multiple disciplines, using theoretical and experimental methods to determine the physical properties of materials in the nanoscale size range. Interesting properties include the structural, electronic, optical, and thermal behavior of nanomaterials, electrical and thermal conductivity, the forces between nanoscale objects, and the transition between classical and quantum behavior. Nanophysics has now become an independent branch of physics, simultaneously expanding into many new areas and playing a vital role in fields that were once the domain of engineering, chemical, or life sciences. Nanoscience and nanotechnology are all about relating and exploiting phenomena for materials having one, two or three dimensions reduced to the nanoscale. Breakthroughs in nanotechnology require a firm grounding in the principles of nanophysics. It is intended to fulfill a crucial purpose. Nanophysics aims to connect scientists with disparate interests to begin interdisciplinary projects and incorporate the theory and methodology of other fields into their work.


title: Asymmetric hysteresis for probing Dzyaloshinskii-Moriya interaction

abstract: The interfacial Dzyaloshinskii-Moriya interaction (DMI) is intimately related to the prospect of superior domain-wall dynamics and the formation of magnetic skyrmions. Although some experimental efforts have been recently proposed to quantify these interactions and the underlying physics, it is still far from trivial to address the interfacial DMI. Inspired by the reported tilt of the magnetization of the side edge of a thin film structure, we here present a quasi-static, straightforward measurement tool. By using laterally asymmetric triangular-shaped microstructures, it is demonstrated that interfacial DMI combined with an in-plane magnetic field yields a unique and significant shift in magnetic hysteresis. By systematic variation of the shape of the triangular objects combined with a droplet model for domain nucleation, a robust value for the strength and sign of interfacial DMI is obtained. This method gives immediate and quantitative access to DMI, enabling a much faster exploration of new DMI systems for future nanotechnology.


title: New Insight into the Ground State of FePc: A Diffusion Monte Carlo Study

abstract: We have applied DMC to evaluate relative stability of the possible electronic configurations of an isolated FePc under $D_{4h}$ symmetry, considering some fixed nodes generated from different methods. They predict $A_{2g}$ ground state consistently, supporting preceding DFT studies, with confidence overcoming the ambiguity about exchange-correlation (XC) functionals. By comparing DMC with several XC, we clarified the importance of the short range exchange to describe the relative stability. We examined why the predicted $A_{2g}$ is excluded from possible ground states in the recent ligand field based model. Simplified assumptions made in the superposition model are identified to give unreasonably less energy gain for $A_{2g}$ when compared with the reality. The state is found to have possible reasons for the stabilization, reducing the occupations from an unstable anti-bonding orbital, avoiding double occupation of a spatially localized orbital, and gaining exchange energy by putting a triplet spin pair in degenerate orbitals.


title: On parametrized verification of asynchronous, shared-memory pushdown systems

abstract: We consider the model of parametrized asynchronous shared-memory pushdown systems, as introduced in [Hague'11]. In a series of recent papers it has been shown that reachability in this model is PSPACE-complete [Esparza, Ganty, Majumdar'13] and that liveness is decidable in NEXPTIME [Durand-Gasselin, Esparza, Ganty, Majumdar'15]. We show here that the liveness problem is PSPACE-complete. We also introduce the universal reachability problem. We show that it is decidable, and coNEXPTIME-complete. Finally, using these results, we prove that the verifying regular properties of traces of executions, satisfying some stuttering condition, is also decidable in NEXPTIME for this model.


title: Domain Specific Language for Modular Knitting Pattern Definitions: Purl

abstract: Purl is a language to be used for modular definition and verification of knitting patterns. The syntax is similar to the standard knitting pattern notation provided by the Craft Yarn Council. Purl provides constructs not available in the standard notation to allow reuse of segments of patterns. This report describes the basics of knitting and hand-knitting patterns. A knitting pattern language more terse than the standard notation is presented with the implementation of a compiler to this standard.


title: A note on the upper critical field of Sr$_2$RuO$_4$ under strain

abstract: In the light of the recently discussed mechanism for suppression of superconductivity in multi-orbital systems called inter-orbital effect, here we extend our analysis of the upper critical field in Sr$_2$RuO$_4$ now under strain. We show that the presence of the standard orbital effect and the new mechanism can consistently account for the qualitative changes observed in the upper critical field for the strained system, in particular, combining the overall enhancement of the critical field with the reduction in the anisotropy. The proposed picture holds for a triplet superconducting state, showing that a singlet state is not the only possibility to account for the observations. We suggest further experiments in order to clarify our understanding about the superconducting phase of Sr$_2$RuO$_4$.


title: The principle of the mutual energy

abstract: Advanced potential solution of Maxwell equations isn't often accepted. We have proven if without advanced potential, it is not possible to satisfy the Maxwell equations. We also shown that it is not the Poynting vector related energy current transferring energy in the space and it is the mutual energy really did that. A important result of the mutual energy theorem is that the advanced potential can suck energy from the transmitter. This energy is equal to the energy received at the receiver. Hence a transmitter can not send any energy out without the receiver. For two remote objects, the energy is transferred only can by the mutual energy of a retarded potential from the source together with an advanced potential from the sink. If the sucked energy is discrete, the summation of mutual energy current of the infinite background atoms or currents, which can be seen as receivers, is a random process. This means that the photon energy sent by the transmitter is actually grabbed by the receiver. Hence the photon from very beginning knows their destination. This receiver send advanced potential to the transmitter. This explanation also avoided the wave function collapse. The retarded potential first reached the receiver, cause the current in the receiver, the current of receiver send a advanced potential to the transmitter with a reversed time, in the same time, a photon minus-time-instantly runs from receiver to transmitter. In our normal feeling, the photon is still runs from the transmitter to the receiver with a positive time. How to transfer superluminal signal using advanced potential is also discussed.Maxwell equations as principle of the theory of the electromagnetic fields is replaced by the mutual energy principle.


title: Glueball Spectra from a Matrix Model of Pure Yang-Mills Theory

abstract: We present variational estimates for the low-lying energies of a simple matrix model that approximates $SU(3)$ Yang-Mills theory on a three-sphere of radius $R$. By fixing the ground state energy, we obtain the (integrated) renormalization group (RG) equation for the Yang-Mills coupling $g$ as a function of $R$. This RG equation allows to estimate the masses of other glueball states, which we find to be in excellent agreement with lattice simulations.


title: Spatio-temporal interference of photo electron wave packets and time scale of non-adiabatic transition in high-frequency regime

abstract: The method of the envelope Hamiltonian [K. Toyota, U. Saalmann, and J. M. Rost, New J. Phys. {\bf 17}, 073005~(2015)] is applied to further study a detachment dynamics of a model negative ion in one-dimension in high-frequency regime. This method is based on the Floquet approach, but the time-dependency of an envelope function is explicitly kept for arbitrary pulse durations. Therefore, it is capable of describing not only a photo absorption/emission but also a non-adiabatic transition which is induced by the time-varying envelope of the pulse. It was shown that the envelope Hamiltonian accurately retrieves the results obtained by the time-dependent Schr\"odinger equation, and underlying physics were well understood by the adiabatic approximation based on the envelope Hamiltonian. In this paper, we further explore two more aspects of the detachment dynamics, which were not done in our previous work. First, we find out features of both a {\it spatial} and {\it temporal} interference of photo electron wave packets in a photo absorption process. We conclude that both the interference mechanisms are universal in ionization dynamics in high-frequency regime. To our knowledge, it is first time that both the interference mechanisms in high-frequency regime are extracted from the first principle. Second, we extract a pulse duration which maximize a yield of the non-adiabatic transition as a function of a pulse duration. It is shown that it becomes maximum when the pulse duration is comparable to a time-scale of an electron.


title: X-Events and Their Interpretation

abstract: We point out that when doubly-charged bileptons are pair produced at the LHC, kinematics dictate that they are both almost at rest in the lab frame and therefore their decays lead to final state muons in a characteristic X-shape with only very tiny track curvature because of the high muon energies. Such X-events have essentially no standard model background and provide a smoking gun for the 331 model.


title: Comment on "Scientific Regress" (First Things May 2016)

abstract: Response to William A. Wilson on the limits and fallibility of science.


title: The Story of the Neutrino

abstract: This is an elementary review of the history and physics of neutrinos. The story of the discovery of neutrino mass through neutrino oscillations is described in some detail. Experiments on solar neutrinos and atmospheric neutrinos played an important part. Recent advances are summarized and future developments are indicated.


title: Extraction of harmonics from trigonometric polynomials by amplitude and phase operators

abstract: Extraction of harmonics of a given order from real trigonometric polynomials (signals) is one of the main problems in harmonic analysis. It has many applications in physics, radio and electrical engineering, in particular, in filtration of harmonic signals of different nature. There exist many methods (mainly, approximative) for solution of this problem. The most common ones are spectral methods based on Fourier transform and other resonance principles. In this paper we propose a new method for extracting harmonics by amplitude and phase transformation of trigonometric polynomials. These transformations use the two simplest operations -- multiplication by a real constant and phase shift -- to obtain polynomials similar to the initial ones. A harmonic is extracted by an amplitude and phase operator that simply overlays (sums up) a finite number of such similar polynomials. The overlay method enables us to obtain precise analytical formulas for calculating harmonics of a given order.


title: Valley-contrasting orbital angular momentum in photonic valley crystals

abstract: Valley, as a degree of freedom, has been exploited to realize valley-selective Hall transport and circular dichroism in two-dimensional layered materials. On the other hand, orbital angular momentum of light with helical phase distribution has attracted great attention for its unprecedented opportunity to optical communicagtions, atom trapping, and even nontrivial topology engineering. Here, we reveal valley-contrasting orbital angular momentum in all-dielectric photonic valley crystals. Selective excitation of valley chiral bulk states is realized by sources carrying orbital angular momentum with proper chirality. Valley dependent edge states, predictable by nonzero valley Chern number, enable to suppress the inter-valley scattering along zigzag boundary, leading to broadband robust transmission in Z-shape bend without corner morphological optimization. Our work may open up a new door towards the discovery of novel quantum states and the manipulation of spin-orbit interaction of light in nanophotonics.


title: Learning Nash Equilibrium for General-Sum Markov Games from Batch Data

abstract: This paper addresses the problem of learning a Nash equilibrium in $\gamma$-discounted multiplayer general-sum Markov Games (MG). A key component of this model is the possibility for the players to either collaborate or team apart to increase their rewards. Building an artificial player for general-sum MGs implies to learn more complex strategies which are impossible to obtain by using techniques developed for two-player zero-sum MGs. In this paper, we introduce a new definition of $\epsilon$-Nash equilibrium in MGs which grasps the strategy's quality for multiplayer games. We prove that minimizing the norm of two Bellman-like residuals implies the convergence to such an $\epsilon$-Nash equilibrium. Then, we show that minimizing an empirical estimate of the $L_p$ norm of these Bellman-like residuals allows learning for general-sum games within the batch setting. Finally, we introduce a neural network architecture named NashNetwork that successfully learns a Nash equilibrium in a generic multiplayer general-sum turn-based MG.


title: Using Sequence Ensembles for Seeding Alignments of MinION Sequencing Data

abstract: Oxford Nanopore MinION sequencer is currently the smallest sequencing device available. While being able to produce very long reads (reads of up to 100~kbp were reported), it is prone to high sequencing error rates of up to 30%. Since most of these errors are insertions or deletions, it is very difficult to adapt popular seed-based algorithms designed for aligning data sets with much lower error rates. Base calling of MinION reads is typically done using hidden Markov models. In this paper, we propose to represent each sequencing read by an ensemble of sequences sampled from such a probabilistic model. This approach can improve the sensitivity and false positive rate of seeding an alignment compared to using a single representative base call sequence for each read.


title: Irreversible Brownian heat engine

abstract: We model a Brownian heat engine as a Brownian particle that hops in a periodic ratchet potential where the ratchet potential is coupled with a linearly decreasing background temperature. It is shown that the efficiency of such Brownian heat engine is far from Carnot efficiency even at quaistatic limit. At quasistatic limit, the efficiency of the heat engine approaches the efficiency of endoreversible engine $\eta=1-\sqrt{{T_{c}/T_{h}}}$ \cite{c18}. On the other hand, the maximum power efficiency of the engine approaches $\eta^{MAX}=1-({T_{c}/T_{h}})^{1\over 4}$. Moreover, the dependence of the current as well as the efficiency on the model parameters is explored analytically by omitting the heat exchange via the kinetic energy. In this case we show that the optimized efficiency always lies between the efficiently at quaistatic limit and the efficiency at maximum power. On the other hand, the efficiency at maximum power is always less than the optimized efficiency since the fast motion of the particle comes at the expense of the energy cost. The role of time on the performance of the motor is also explored via numerical simulations. Our numerical results depict that the time $t$ as well as the external load dictate the direction of the particle velocity. Moreover the performance of the heat engine improves with time. At large $t$ (steady state),the velocity, the efficiency and the coefficient of performance of the refrigerator attain their maximum value.


title: Optimal Control of a Tuberculosis Model with State and Control Delays

abstract: We introduce delays in a tuberculosis (TB) model, representing the time delay on the diagnosis and commencement of treatment of individuals with active TB infection. The stability of the disease free and endemic equilibriums is investigated for any time delay. Corresponding optimal control problems, with time delays in both state and control variables, are formulated and studied. Although it is well-known that there is a delay between two to eight weeks between TB infection and reaction of body's immune system to tuberculin, delays for the active infected to be detected and treated, and delays on the treatment of persistent latent individuals due to clinical and patient reasons, which clearly justifies the introduction of time delays on state and control measures, our work seems to be the first to consider such time-delays for TB and apply time-delay optimal control to carry out the optimality analysis.


title: Epimenides, G\"odel, Turing: an Eternal G\"olden Tangle

abstract: The Halting Problem is a version of the Liar's Paradox.


title: Magnetization switching of FePt nanoparticle recording medium by femtosecond laser pulses

abstract: Manipulation of magnetization with ultrashort laser pulses is promising for information storage device applications. The dynamic of the magnetization response depends on the energy transfer from the photons to the spins during the initial laser excitation. A material of special interest for magnetic storage is FePt nanoparticles , on which optical writing with optical angular momentum was demonstrated recently by Lambert et al., although the mechanism remained unclear. Here we investigate experimentally and theoretically the all-optical switching of FePt nanoparticles. We show that the magnetization switching is a stochastic process. We develop a complete multiscale model which allows us to optimize the number of laser shots needed to write the magnetization of high anisotropy FePt nanoparticles in our experiments. We conclude that only angular momentum induced optically by the inverse Faraday effect will provide switching with one single femtosecond laser pulse.


title: Increasing robustness of pairwise methods for effective connectivity in Magnetic Resonance Imaging by using fractional moment series of BOLD signal distributions

abstract: Estimating causal interactions in the brain from functional magnetic resonance imaging (fMRI) data remains a challenging task. Multiple studies have demonstrated that all current approaches to determine direction of connectivity perform poorly even when applied to synthetic fMRI datasets. Recent advances in this field include methods for pairwise inference, which involve creating a sparse connectome in the first step, and then using a classifier in order to determine the directionality of connection between of every pair of nodes in the second step. In this work, we introduce an advance to the second step of this procedure, by building a classifier based on fractional moments of the BOLD distribution combined into cumulants. The classifier is trained on datasets generated under the Dynamic Causal Modeling (DCM) generative model. The directionality is inferred based upon statistical dependencies between the two node time series, e.g. assigning a causal link from time series of low variance to time series of high variance. Our approach outperforms or performs as well as other methods for effective connectivity when applied to the benchmark datasets. Crucially, it is also more resilient to confounding effects such as differential noise level across different areas of the connectome.


title: Perspective: Ultrafast magnetism and THz spintronics

abstract: This year the discovery of femtosecond demagnetization by laser pulses is 20 years old. For the first time this milestone work by Bigot and coworkers gave insight in a very direct way into the time scales of microscopic interactions that connect the spin and electron system. While intense discussions in the field were fueled by the complexity of the processes in the past, it now became evident that it is a puzzle of many different parts. Rather than giving an overview that has been presented in previous reviews on ultrafast processes in ferromagnets, this perspective will show that with our current depth of knowledge the first real applications are on their way: THz spintronics and all-optical spin manipulation are becoming more and more feasible. The aim of this perspective is to point out where we can connect the different puzzle pieces of understanding gathered over 20 years to develop novel applications. based on many observations in a large number of experiments. Differences in the theoretical models arise from the localized and delocalized nature of ferromagnetism. Transport effects are intrinsically non-local in spintronic devices and at interfaces. We review the need for multiscale modeling to address processes starting from electronic excitation of the spin system on the picometer length scale and sub-femtosecond time scale, to spin wave generation, and towards the modeling of ultrafast phase transitions that altogether determine the response time of the ferromagnetic system. Today, our current understanding gives rise to the first real applications of ultrafast spin physics for ultrafast magnetism control: THz spintronic devices. This makes the field of ultrafast spin-dynamics an emerging topic open for many researchers right now.


title: On finite generation of the section ring of the determinant of cohomology line bundle

abstract: For $C$ a stable curve of arithmetic genus $g\ge 2$, and $\mathcal{D}$ the determinant of cohomology line bundle on $\operatorname{Bun}_{\operatorname{SL}(r)}(C)$, we show the section ring for the pair $(\operatorname{Bun}_{\operatorname{SL}(r)}(C), \mathcal{D})$ is finitely generated. Three applications are given.


title: Hilbertian Interpolation

abstract: I want to prove that all classical techniques of interpolation and approximation as Lagrange, Taylor, Hermite interpolations Beziers interpolants, Quasi interpolants, Box splines and others (radial splines, simplicial splines) are derived from a \textbf{unique} simple hilbertian scheme. For sake of simplicity, we shall consider only elementary examples which could be easily generalized.


title: On Banachic Kernels and Approximation Theory

abstract: In this paper, I generalize a previous one about hilbertian kernels and approximation theory


title: Traces of Besov, Triebel-Lizorkin and Sobolev spaces on metric spaces

abstract: We establish trace theorems for function spaces defined on general Ahlfors regular metric spaces $Z$. The results cover the Triebel-Lizorkin spaces and the Besov spaces for smoothness indices $s<1,$ as well as the first order Haj{\l}asz-Sobolev space $M^{1,p}(Z)$. They generalize the classical results from the Euclidean setting, since the traces of these function spaces onto any closed Ahlfors regular subset $F \subset Z$ are Besov spaces defined intrinsically on $F$. Our method employs the definitions of the function spaces via hyperbolic fillings of the underlying metric space.


title: Maximal entropy distribution functions from generalized R\'enyi entropy

abstract: New class of reference distribution functions for numerical approximation of the solution of the Fokker-Planck equations associated to the charged particle dynamics in tokamak are studied. The reference distribution functions are obtained by maximization of the generalized Renyi entropy under scale-invariant restrictions. Explicit analytic form, with algebraic decay, that is a generalization of the previous distribution with exponential tails was derived.


title: Impulse control maximising average cost per unit time: a non-uniformly ergodic case

abstract: This paper studies maximisation of an average-cost-per-unit-time ergodic functional over impulse strategies controlling a Feller-Markov process. The uncontrolled process is assumed to be ergodic but, unlike the extant literature, the convergence to invariant measure does not have to be uniformly geometric in total variation norm; in particular, we allow for non-uniform geometric or polynomial convergence. Cost of an impulse may be unbounded, e.g., proportional to the distance the process is shifted. We show that the optimal value does not depend on the initial state and provide optimal or $\ve$-optimal strategies.


title: Spatial Variations of the Extinction Law in the Galactic Disk from Infrared Observations

abstract: Infrared photometry in the J (1.2 microns), H (1.7 microns), Ks (2.2 microns) bands from the 2MASS catalogue and in the W1 (3.4 microns), W2 (4.6 microns), W3 (12 microns), W4 (22 microns) bands from the WISE catalogue is used to reveal the spatial variations of the interstellar extinction law in the infrared near the midplane of the Galaxy by the method of extrapolation of the extinction law applied to clump giants. The variations of the coefficients E(H-W1)/E(H-Ks), E(H-W2)/E(H-Ks), E(H-W3)/E(H-Ks), and E(H-W4)/E(H-Ks) along the line of sight in 2 deg per 2 deg squares of the sky centered at b=0 and l=20, 30, ..., 330, 340 deg as well as in several 4 deg per 4 deg squares with |b|=10 are considered. The results obtained here agree with those obtained by Zasowski et al. in 2009 using 2MASS and Spitzer-IRAC photometry for the same longitudes and similar photometric bands, confirming their main result: in the inner (relative to the Sun) Galactic disk, the fraction of fine dust increases with Galactocentric distance (or the mean dust grain size decreases). However, in the outer Galactic disk that was not considered by Zasowski et al., this trend is reversed: at the disk edge, the fraction of coarse dust is larger than that in the solar neighborhood. This general Galactic trend seems to be explained by the influence of the spiral pattern: its processes sort the dust by size and fragment it so that coarse and fine dust tend to accumulate, respectively, at the outer and inner (relative to the Galactic center) edges of the spiral arms. As a result, fine dust may exist only in the part of the Galactic disk far from both the Galactic center and the edge, while coarse dust dominates at the Galactic center, at the disk edge, and outside the disk.


title: Recurrent Neural Networks for Dialogue State Tracking

abstract: This paper discusses models for dialogue state tracking using recurrent neural networks (RNN). We present experiments on the standard dialogue state tracking (DST) dataset, DSTC2. On the one hand, RNN models became the state of the art models in DST, on the other hand, most state-of-the-art models are only turn-based and require dataset-specific preprocessing (e.g. DSTC2-specific) in order to achieve such results. We implemented two architectures which can be used in incremental settings and require almost no preprocessing. We compare their performance to the benchmarks on DSTC2 and discuss their properties. With only trivial preprocessing, the performance of our models is close to the state-of- the-art results.


title: Three-dimensional simulations of turbulent convective mixing in ONe and CO classical nova explosions

abstract: Classical novae are thermonuclear explosions that take place in the envelopes of accreting white dwarfs in binary systems. The material piles up under degenerate conditions, driving a thermonuclear runaway. The energy released by the suite of nuclear processes operating at the envelope heats the material up to peak temperatures about 100 - 400 MK. During these events, about 10-3 - 10-7 Msun, enriched in CNO and, sometimes, other intermediate-mass elements (e.g., Ne, Na, Mg, Al) are ejected into the interstellar medium. To account for the gross observational properties of classical novae (in particular, the large concentrations of metals spectroscopically inferred in the ejecta), models require mixing between the (solar-like) material transferred from the secondary and the outermost layers (CO- or ONe-rich) of the underlying white dwarf. Recent multidimensional simulations have demonstrated that Kelvin-Helmholtz instabilities can naturally produce self-enrichment of the accreted envelope with material from the underlying white dwarf at levels that agree with observations. However, the feasibility of this mechanism has been explored in the framework of CO white dwarfs, while mixing with different substrates still needs to be properly addressed. Three-dimensional simulations of mixing at the core-envelope interface during nova outbursts have been performed with the multidimensional code FLASH, for two types of substrates: CO- and ONe-rich. We show that the presence of an ONe-rich substrate, as in "neon novae", yields larger metallicity enhancements in the ejecta, compared to CO,rich substrates (i.e., non-neon novae). A number of requirements and constraints for such 3-D simulations (e.g., minimum resolution, size of the computational domain) are also outlined.


title: Influence of Fermi arc states and double Weyl node on tunneling in a Dirac semimetal

abstract: Most theoretical studies of tunneling in Dirac and the closely related Weyl semimetals have modeled these materials as single Weyl nodes described by the three-dimensional Dirac equation $H = v_f \vec{p}\cdot\vec{\sigma}$. The influence of scattering between the different valleys centered around different Weyl nodes, and the Fermi arc states which connect these nodes are hence not evident from these studies. In this work we study the tunneling in a thin film system of the Dirac semimetal $\text{Na}_3\text{Bi}$ consisting of a central segment with a gate potential, sandwiched between identical semi-infinite source and drain segments. The model Hamiltonian we use for $\text{Na}_3\text{Bi}$ gives, for each spin, two Weyl nodes separated in $k$-space symmetrically about $k_z=0$. The presence of a top and bottom surface in the thin film geometry results in the appearance of Fermi arc states and energy subbands. We show that (for each spin) the presence of two Weyl nodes and the Fermi arc states result in enhanced transmission oscillations, and finite transmission even when the energy falls within the \textit{bulk} band gap in the central segment respectively. These features are not evident in single Weyl node models.


title: Hardy inequality and fractional Leibnitz rule for perturbed Hamiltonians on the line

abstract: We consider the following perturbed Hamiltonian $\mathcal{H}= -\partial_x^2 + V(x)$ on the real line. The potential $V(x)$ is a real - valued function of short range type. We study the equivalence of classical homogeneous Sobolev type spaces $\dot{H}^s_p$, $p \in (1,\infty)$ and the corresponding perturbed homogeneous Sobolev spaces associated with the perturbed Hamiltonian. It is shown that the assumption zero is not a resonance guarantees that the perturbed and unperturbed homogeneous Sobolev norms of order $s = \gamma - 1 \in [0,1/p)$ are equivalent. As a corollary, the corresponding wave operators leave classical homogeneous Sobolev spaces of order $s \in [0,1/p)$ invariant.


title: Design and Experimental Characterization of Spiral Resonator variants-TTSR , NBSR

abstract: In this paper we have carried out numerical simulation and experimental characterization of two variants of SR (Spiral resonator) viz. TTSR (Two-Turn Spiral Resonator) and NBSR (Non Bianisotropic spiral resonator). The numerical simulation of these develop magnetic inclusion structures are done with commercially available Finite Element Method (FEM) based Ansys HFSS-13 simulation software, and transmission experiment using S-parameter retrieval technique. The experimental characterization of those fabricated structures is used in the parallel plate waveguide experiment (not reported earlier) and then validate our experimental results with the Ansys HFSS-13 simulation software; which in not available in any published literature.


title: Dynamical phase transition in large deviation statistics of the Kardar-Parisi-Zhang equation

abstract: We study the short-time behavior of the probability distribution $\mathcal{P}(H,t)$ of the surface height $h(x=0,t)=H$ in the Kardar-Parisi-Zhang (KPZ) equation in $1+1$ dimension. The process starts from a stationary interface: $h(x,t=0)$ is given by a realization of two-sided Brownian motion constrained by $h(0,0)=0$. We find a singularity of the large deviation function of $H$ at a critical value $H=H_c$. The singularity has the character of a second-order phase transition. It reflects spontaneous breaking of the reflection symmetry $x \leftrightarrow -x$ of optimal paths $h(x,t)$ predicted by the weak-noise theory of the KPZ equation. At $|H|\gg |H_c|$ the corresponding tail of $\mathcal{P}(H)$ scales as $-\ln \mathcal{P} \sim |H|^{3/2}/t^{1/2}$ and agrees, at any $t>0$, with the proper tail of the Baik-Rains distribution, previously observed only at long times. The other tail of $\mathcal{P}$ scales as $-\ln \mathcal{P} \sim |H|^{5/2}/t^{1/2}$ and coincides with the corresponding tail for the sharp-wedge initial condition.


title: Stability of closed solutions to the VFE hierarchy with application to the Hirota equation

abstract: The Vortex Filament Equation (VFE) is part of an integrable hierarchy of filament equations. Several equations part of this hierarchy have been derived to describe vortex filaments in various situations. Inspired by these results, we develop a general framework for studying the existence and the linear stability of closed solutions of the VFE hierarchy. The framework is based on the correspondence between the VFE and the nonlinear Schr\"odinger (NLS) hierarchies. Our results show that it is possible to establish a connection between the AKNS Floquet spectrum and the stability properties of the solutions of the filament equations. We apply our machinery to solutions of the filament equation associated to the Hirota equation. We also discuss how our framework applies to soliton solutions.


title: Alternating Anderson-Richardson method: An efficient alternative to preconditioned Krylov methods for large, sparse linear systems

abstract: We present the Alternating Anderson-Richardson (AAR) method: an efficient and scalable alternative to preconditioned Krylov solvers for the solution of large, sparse linear systems on high performance computing platforms. Specifically, we generalize the recently proposed Alternating Anderson-Jacobi (AAJ) method (Pratapa et al., J. Comput. Phys. (2016), 306, 43--54) to include preconditioning, discuss efficient parallel implementation, and provide serial MATLAB and parallel C/C++ implementations. In serial applications to nonsymmetric systems, we find that AAR is comparably robust to GMRES, using the same preconditioning, while often outperforming it in time to solution; and find AAR to be more robust than Bi-CGSTAB for the problems considered. In parallel applications to the Helmholtz and Poisson equations, we find that AAR shows superior strong and weak scaling to GMRES, Bi-CGSTAB, and Conjugate Gradient (CG) methods, using the same preconditioning, with consistently shorter times to solution at larger processor counts. Finally, in massively parallel applications to the Poisson equation, on up to 110,592 processors, we find that AAR shows superior strong and weak scaling to CG, with shorter minimum time to solution. We thus find that AAR offers a robust and efficient alternative to current state-of-the-art solvers, with increasing advantages as the number of processors grows.


title: Dynamic Watermarking: Active Defense of Networked Cyber-Physical Systems

abstract: The coming decades may see the large scale deployment of networked cyber-physical systems to address global needs in areas such as energy, water, healthcare, and transportation. However, as recent events have shown, such systems are vulnerable to cyber attacks. Being safety critical, their disruption or misbehavior can cause economic losses or injuries and loss of life. It is therefore important to secure such networked cyber-physical systems against attacks. In the absence of credible security guarantees, there will be resistance to the proliferation of cyber-physical systems, which are much needed to meet global needs in critical infrastructures and services. This paper addresses the problem of secure control of networked cyber-physical systems. This problem is different from the problem of securing the communication network, since cyber-physical systems at their very essence need sensors and actuators that interface with the physical plant, and malicious agents may tamper with sensors or actuators, as recent attacks have shown. We consider physical plants that are being controlled by multiple actuators and sensors communicating over a network, where some sensors could be "malicious," meaning that they may not report the measurements that they observe. We address a general technique by which the actuators can detect the actions of malicious sensors in the system, and disable closed-loop control based on their information. This technique, called "watermarking," employs the technique of actuators injecting private excitation into the system which will reveal malicious tampering with signals. We show how such an active defense can be used to secure networked systems of sensors and actuators.


title: The radio evolution of supernova SN$\,$2008iz in M$\,$82

abstract: We report on multi-frequency Very Large Array (VLA) and Very Long Baseline Interferometry (VLBI) radio observations for a monitoring campaign of supernova SN$\,$2008iz in the nearby irregular galaxy M$\,$82. We fit two models to the data, a simple time power-law, S$\propto$ t$^\beta$, and a simplified Weiler model, yielding decline indices, $\beta$ of -1.22$\pm$0.07 (days 100-1500) and -1.41$\pm$0.02 (days 76-2167), respectively. The late-time radio light curve evolution shows flux-density flares at $\sim\,$970 and $\sim\,$1400 days which are a factor of $\sim\,$2 and $\sim\,$4 higher than the expected flux, respectively. The later flare, besides being brighter, does not show signs of decline at least from results examined so far (2014 January 23, day 2167). We derive the spectral index, $\alpha$, S$\propto$ $\nu^\alpha$ for frequencies 1.4 to 43$\,$GHz for SN$\,$2008iz during the period from $\sim\,$430 to 2167 days after the supernova explosion. The value of $\alpha$ shows no signs of evolution and remains steep $\approx$$-$1 throughout the period, unlike that of SN$\,$1993J which started flattening at $\sim\,$day 970. From the 4.8 and 8.4$\,$GHz VLBI images, the supernova expansion is seen to start with shell like structure that gets more and more asymmetric, then breaks up in the later epochs with bright structures dominating the southern part of the ring. This structural evolution differs significantly from SN$\,$1993J, which remains circularly symmetric over 4000 days after the explosion. The VLBI 4.8 and 8.4$\,$GHz images are used to derive a deceleration index, $m$, for SN$\,$2008iz, of 0.86$\,\pm\,$0.02, and the average expansion velocity between days 73 and 1400 as (12.1$\pm$0.2)$\times$10$^{3}\,$km$\,$s$^{-1}$.


title: Monotonicity-preserving finite element schemes based on differentiable nonlinear stabilization

abstract: In this work, we propose a nonlinear stabilization technique for scalar conservation laws with implicit time stepping. The method relies on an artificial diffusion method, based on a graph-Laplacian operator. It is nonlinear, since it depends on a shock detector. The same shock detector is used to gradually lump the mass matrix. The resulting method is LED, positivity preserving, linearity preserving, and also satisfies a global DMP. Lipschitz continuity has also been proved. However, the resulting scheme is highly nonlinear, leading to very poor nonlinear convergence rates. We propose a smooth version of the scheme, which leads to twice differentiable nonlinear stabilization schemes. It allows one to straightforwardly use Newton's method and obtain quadratic convergence. In the numerical experiments, steady and transient linear transport, and transient Burgers' equation have been considered in 2D. Using the Newton method with a smooth version of the scheme we can reduce 10 to 20 times the number of iterations of Anderson acceleration with the original non-smooth scheme. In any case, these properties are only true for the converged solution, but not for iterates. In this sense, we have also proposed the concept of projected nonlinear solvers, where a projection step is performed at the end of every nonlinear iteration onto a FE space of admissible solutions. The space of admissible solutions is the one that satisfies the desired monotonic properties (maximum principle or positivity).


title: A new method to assign galaxy cluster membership using photometric redshifts

abstract: We introduce a new effective strategy to assign group and cluster membership probabilities $P_{mem}$ to galaxies using photometric redshift information. Large dynamical ranges both in halo mass and cosmic time are considered. The method takes the magnitude distribution of both cluster and field galaxies as well as the radial distribution of galaxies in clusters into account using a non-parametric formalism and relies on Bayesian inference to take photometric redshift uncertainties into account. We successfully test the method against 1,208 galaxy clusters within redshifts $z=0.05-2.58$ and masses $10^{13.29-14.80}~M_\odot$ drawn from wide field simulated galaxy mock catalogs developed for the Euclid mission. Median purity $(55^{+17}_{-15})\%$ and completeness $(95^{+5}_{-10})\%$ are reached for galaxies brighter than 0.25$L_\ast$ within $r_{200}$ of each simulated halo and for a statistical photometric redshift accuracy $\sigma((z_s-z_p)/(1+z_s))=0.03$. The mean values $p=56\%$ and $c=93\%$ have sub-percent uncertainties. Accurate photometric redshifts ($\sigma((z_s-z_p)/(1+z_s))\lesssim0.05$) and robust estimates for the cluster redshift and the center coordinates are required. The method is applied to derive accurate richness estimates. A statistical comparison between the true ($N_{\rm true}$) vs estimated richness ($\lambda=\sum P_{mem}$) yields on average to unbiased results, $Log(\lambda/N_{\rm true})=-0.0051\pm0.15$. The scatter around the mean of the logarithmic difference between $\lambda$ and the halo mass is 0.10~dex, for massive halos $\gtrsim10^{14.5}~M_\odot$. Our estimates could be useful to calibrate independent cluster mass estimates from weak lensing, SZ, and X-ray studies. Our method can be applied to any list of galaxy clusters or groups in both present and forthcoming surveys such as SDSS, CFHTLS, DES, LSST, and Euclid.


title: Plastic response and correlations in athermally sheared amorphous solids

abstract: The onset of irreversible deformation in low-temperature amorphous solids is due to the accumulation of elementary events, consisting of spacially and temporally localized atomic rearrangements involving only a few tens of atoms. Recently, numerical and experimental work addressed the issue of spatio-temporal correlations between these plastic events. Here, we provide further insight into these correlations by investigating, via molecular dynamics (MD) simulations, the plastic response of a two-dimensional amorphous solid to artificially triggered local shear transformations. We show that while the plastic response is virtually absent in as-quenched configurations, it becomes apparent if a shear strain was previously imposed on the system. Plastic response has a four-fold symmetry which is characteristic of the shear stress redistribution following the local transformation. At high shear rate we report evidence for a fluctuation-dissipation relation, connecting plastic response and correlation, which seems to break down if lower shear rates are considered.


title: Vortex Lattices in the Bose-Fermi Superfluid Mixture

abstract: In this letter we show that the vortex lattice structure in the Bose-Fermi superfluid mixture can undergo a sequence of structure transitions when the Fermi superfluid is tuned from the BCS regime to the BEC regime. This is due to different vortex core structure of the Fermi superfluid in the BCS regime and in the BEC regime. In the former the vortex core is nearly filled, while the density at the vortex core gradually decreases until it empties out at the BEC regime. Therefore, with the density-density interaction between the Bose and the Fermi superfluids, the two sets of vortex lattices interact stronger in the BEC regime that yields the structure transition of vortex lattices. In view of recent realization of this superfluid mixture and vortices therein, our theoretical predication can be verified experimentally in near future.


title: Simulations of 3D-Si sensors for the innermost layer of the ATLAS pixel upgrade

abstract: The LHC is expected to reach luminosities up to 3000fb-1 and the innermost layer of the ATLAS upgrade plans to cope with higher occupancy and to decrease the pixel size. 3D-Si sensors are a good candidate for the innermost layer of the ATLAS pixel upgrade since they exhibit good performance under high fluences and the new designs will have smaller pixel size to fulfill the electronics expectations. This paper reports TCAD simulations of the 3D-Si sensors designed at IMB-CNM with non passing-through columns that are being fabricated for the next innermost layer of the ATLAS pixel upgrade, shows the charge collection response before and after irradiation, and the response of 3D-Si sensors located at large $\eta$ angles.


title: Modeling the shock-cloud interaction in SN 1006: unveiling the origin of nonthermal X-ray and gamma-ray emission

abstract: The supernova remnant SN 1006 is a source of high-energy particles and its southwestern limb is interacting with a dense ambient cloud, thus being a promising region for gamma-ray hadronic emission. We aim at describing the physics and the nonthermal emission associated with the shock-cloud interaction to derive the physical parameters of the cloud (poorly constrained by the data analysis), to ascertain the origin of the observed spatial variations in the spectral properties of the X-ray synchrotron emission, and to predict spectral and morphological features of the resulting gamma-ray emission. We performed 3-D magnetohydrodynamic simulations modeling the evolution of SN 1006 and its interaction with the ambient cloud, and explored different model setups. By applying the REMLIGHT code on the model results, we synthesized the synchrotron X-ray emission, and compared it with actual observations, to constrain the parameters of the model. We also synthesized the leptonic and hadronic gamma-ray emission from the models, deriving constraints on the energy content of the hadrons accelerated at the southwestern limb. We found that the impact of the SN 1006 shock front with a uniform cloud with density 0.5 cm^-3 can explain the observed morphology, the azimuthal variations of the cutoff frequency of the X-ray synchrotron emission, and the shock proper motion in the interaction region. Our results show that the current upper limit for the total hadronic energy in the southwestern limb is 2.5e49 erg.


title: Variational Geometric Approach to Generalized Differential and Fenchel Conjugate Calculi in Convex Analysis

abstract: This paper develops a geometric approach of variational analysis for the case of convex objects considered in locally convex topological spaces and also in Banach space settings. Besides deriving in this way new results of convex calculus, we present an overview of some known achievements with their unified and simplified proofs based on the developed geometric variational schemes.


title: Device-independence for two-party cryptography and position verification

abstract: Quantum communication has demonstrated its usefulness for quantum cryptography far beyond quantum key distribution. One domain is two-party cryptography, whose goal is to allow two parties who may not trust each other to solve joint tasks. Another interesting application is position-based cryptography whose goal is to use the geographical location of an entity as its only identifying credential. Unfortunately, security of these protocols is not possible against an all powerful adversary. However, if we impose some realistic physical constraints on the adversary, there exist protocols for which security can be proven, but these so far relied on the knowledge of the quantum operations performed during the protocols. In this work we give device-independent security proofs of two-party cryptography and Position Verification for memoryless devices under different physical constraints on the adversary. We assess the quality of the devices by observing a Bell violation and we show that security can be attained for any violation of the Clauser-Holt-Shimony-Horne inequality.


title: Comparison of OFDM and Single-Carrier for Large-Scale Antenna Systems

abstract: Large-scale antenna (LSA) or massive multiple-input multiple-output (MIMO) has gained a lot of attention due to its potential to significantly improve system throughput. As a natural evolution from traditional MIMO-orthogonal frequency division multiplexing (OFDM), LSA has been combined with OFDM to deal with frequency selectivity of wireless channels in most existing works. As an alternative approach, single-carrier (SC) has also been proposed for LSA systems due to its low implementation complexity. In this article, a comprehensive comparison between LSA-OFDM and LSA-SC is presented, which is of interest to the waveform design for the next generation wireless systems.


title: On quadratic polynomial mappings $f: \Bbb C^2 \to \Bbb C^2$

abstract: We show that up to linear equivalence, there is only finitely many polynomial quadratic mappings $f:\Bbb C^2\to\Bbb C^2$ and $f:\Bbb R^2\to \Bbb R^2.$ We list all possibilities.


title: Reweighting QCD matrix-element and parton-shower calculations

abstract: We present the implementation and validation of the techniques used to efficiently evaluate parametric and perturbative theoretical uncertainties in matrix-element plus parton-shower simulations within the Sherpa event-generator framework. By tracing the full $\alpha_s$ and PDF dependences, including the parton-shower component, as well as the fixed-order scale uncertainties, we compute variational event weights on-the-fly, thereby greatly reducing the computational costs to obtain theoretical-uncertainty estimates.


title: Dual and chiral objects for optical activity in general scattering directions

abstract: Optically active artificial structures have attracted tremendous research attention. Such structures must meet two requirements: Lack of spatial inversion symmetries and, a condition usually not explicitly considered, the structure shall preserve the helicity of light, which implies that there must be a vanishing coupling between the states of opposite polarization handedness among incident and scattered plane waves. Here, we put forward and demonstrate that a unit cell made from chiraly arranged electromagnetically dual scatterers serves exactly this purpose. We prove this by demonstrating optical activity of such unit cell in general scattering directions.


title: Experimental evidence of independence of nuclear de-channeling length on the particle charge sign

abstract: Under coherent interactions, particles undergo correlated collisions with the crystal lattice and their motion result in confinement in the fields of atomic planes, i.e. particle channeling. Other than coherently interacting with the lattice, particles also suffer incoherent interactions with individual nuclei and may leave their bounded motion, i.e., they de-channel. This latter is the main limiting factor for applications of coherent interactions in crystal-assisted particle steering. We experimentally investigated the nature of dechanneling of 120 GeV/c $e^{-}$ and $e^{+}$ in a bent silicon crystal at H4-SPS external line at CERN. We found out that while channeling efficiency differs significantly for $e^{-}$ ($2\pm2$ $\%$) and $e^{+}$ ($54\pm2$ $\%$), their nuclear dechanneling length is comparable, $(0.6\pm0.1)$ mm for $e^{-}$ and $(0.7\pm0.3)$ mm for $e^{+}$. The experimental proof of the equality of the nuclear dechanneling length for positrons and electrons is interpreted in terms of similar dynamics undergone by the channeled particles in the field of nuclei no matter of their charge.


title: On the (non-)uniqueness of the Levi-Civita solution in the Einstein-Hilbert-Palatini formalism

abstract: We study the most general solution for affine connections that are compatible with the variational principle in the Palatini formalism for the Einstein-Hilbert action (with possible minimally coupled matter terms). We find that there is a family of solutions generalising the Levi-Civita connection, characterised by an arbitrary, non-dynamical vector field ${\cal A}_\mu$. We discuss the mathematical properties and the physical implications of this family and argue that, although there is a clear mathematical difference between these new Palatini connections and the Levi-Civita one, both unparametrised geodesics and the Einstein equation are shared by all of them. Moreover, the Palatini connections are characterised precisely by these two properties, as well as by other properties of its parallel transport. Based on this, we conclude that physical effects associated to the choice of one or the other will not be distinguishable, at least not at the level of solutions or test particle dynamics. We propose a geometrical interpretation for the existence and unobservability of the new solutions.


title: Stochastic Models of Vesicular Sorting in Cellular Organelles

abstract: The proper sorting of membrane components by regulated exchange between cellular organelles is crucial to intra-cellular organisation. This process relies on the budding and fusion of transport vesicles, and should be strongly influenced by stochastic fluctuations considering the relatively small size of many organelles. We identify the perfect sorting of two membrane components initially mixed in a single compartment as a first passage process, and we show that the mean sorting time exhibits two distinct regimes as a function of the ratio of vesicle fusion to budding rates. Low ratio values lead to fast sorting, but results in a broad size distribution of sorted compartments dominated by small entities. High ratio values result in two well defined sorted compartments but is exponentially slow. Our results suggest an optimal balance between vesicle budding and fusion for the rapid and efficient sorting of membrane components, and highlight the importance of stochastic effects for the steady-state organisation of intra-cellular compartments.


title: Survey of Nodeless Regular Almost-Everywhere Holomorphic Solutions for Exactly Solvable Gauss-Reference Liouville Potentials on the Line I.Subsets of Nodeless Jacobi-Seed Solutions Co-Existent with Discrete Energy Spectrum

abstract: The paper collates a complete list of nodeless regular almost-everywhere holomorphic (AEH) solutions for a subset of rational canonical Sturm-Liouville equations (RCSLEs) exactly quantized on a finite interval by classical Jacobi polynomials. The subset was constrained by the requirement that the appropriate Liouville transformation results in the Schrodinger equation on the line. The common remarkable feature of the selected nodeless solutions co-existent with the discrete energy spectrum is that they can be used as seed functions for multi-step 'canonical Liouville-Darboux transformations' (CLDTs) to convert the Gauss-Reference (GRef) potential (appearing in the resultant Schrodinger equation) into its isospectral rational SUSY partners conditionally exactly quantized by the so-called 'Jacobi-Seed' Heine polynomials.


title: Bayesian analysis of immune response dynamics with sparse time series data

abstract: In vaccine development, the temporal profiles of relative abundance of subtypes of immune cells (T-cells) is key to understanding vaccine efficacy. Complex and expensive experimental studies generate very sparse time series data on this immune response. Fitting multi-parameter dynamic models of the immune response dynamics-- central to evaluating mechanisms underlying vaccine efficacy-- is challenged by data sparsity. The research reported here addresses this challenge. For HIV/SIV vaccine studies in macaques, we: (a) introduce novel dynamic models of progression of cellular populations over time with relevant, time-delayed components reflecting the vaccine response; (b) define an effective Bayesian model fitting strategy that couples Markov chain Monte Carlo (MCMC) with Approximate Bayesian Computation (ABC)-- building on the complementary strengths of the two approaches, neither of which is effective alone; (c) explore questions of information content in the sparse time series for each of the model parameters, linking into experimental design and model simplification for future experiments; and (d) develop, apply and compare the analysis with samples from a recent HIV/SIV experiment, with novel insights and conclusions about the progressive response to the vaccine, and how this varies across subjects.


title: Figure-eight choreographies of the equal mass three-body problem with Lennard-Jones-type potentials

abstract: We report on figure-eight choreographic solutions to a system of three identical particles interacting through a potential of Lennard-Jones type, $1/r^{12}-1/r^6$ where $r$ is a distance between the particles. By numerical search, we found there are a multitude of such solutions. A series of them are close to a figure-eight solutions to a homogeneous system with no $1/r^{12}$ term in the potential. The rest are very different from them and have several points with large curvatures in their figure-eight orbits, at which particles are repelled. Here figure-eight choreographies are the periodic motion whose shape is symmetric in both horizontal and vertical axis, starting with an isosceles triangle configuration and going back to an isosceles triangle configuration with opposite direction through Euler configuration. Thus the lobe of this figure-eight may be complex shape and needs not to be convex.


title: A Dissipative Systems Theory for FDTD with Application to Stability Analysis and Subgridding

abstract: This paper establishes a far-reaching connection between the Finite-Difference Time-Domain method (FDTD) and the theory of dissipative systems. The FDTD equations for a rectangular region are written as a dynamical system having the magnetic and electric fields on the boundary as inputs and outputs. Suitable expressions for the energy stored in the region and the energy absorbed from the boundaries are introduced, and used to show that the FDTD system is dissipative under a generalized Courant-Friedrichs-Lewy condition. Based on the concept of dissipation, a powerful theoretical framework to investigate the stability of FDTD methods is devised. The new method makes FDTD stability proofs simpler, more intuitive, and modular. Stability conditions can indeed be given on the individual components (e.g. boundary conditions, meshes, embedded models) instead of the whole coupled setup. As an example of application, we derive a new subgridding method with material traverse, arbitrary grid refinement, and guaranteed stability. The method is easy to implement and has a straightforward stability proof. Numerical results confirm its stability, low reflections, and ability to handle material traverse.


title: A user's guide to cloning systems

abstract: In joint work of the author with Stefan Witzel, a procedure was developed for building new examples of groups in the extended family of R. Thompson's groups, using what we termed \emph{cloning systems}. These new Thompson-like groups can be thought of as limits of families of groups, though unlike other limiting processes, e.g., direct limits, these tend to be well behaved with respect to finiteness properties. In this expository note, we distill the crucial parts of that 50-page paper into a more digestible form, for those curious to understand the construction but less curious about the gritty details. We also give one new example, of a cloning system involving signed symmetric groups.


title: Beyond Linear Fields: the Lie-Taylor Expansion

abstract: The work extends the linear fields' solution of compressible nonlinear magnetohydrodynamics~(MHD) to the case where the magnetic field depends on superlinear powers of position vector, usually but not always, expressed in Cartesian components. Implications of the resulting Lie-Taylor series expansion for physical applicability of the Dolzhansky-Kirchhoff~(D-K) equations are found to be positive. It is demonstrated how resistivity may be included in the D-K model. Arguments are put forward that the D-K equations may be regarded as illustrating properties of nonlinear MHD in the same sense that the Lorenz equations inform about the onset of convective turbulence. It is suggested that the Lie-Taylor series approach may lead to valuable insights into other fluid models.


title: Complexity Bounds of Constant-Space Quantum Computation

abstract: We realize constant-space quantum computation by measure-many two-way quantum finite automata and evaluate their language recognition power by analyzing patterns of their exotic behaviors and by exploring their structural properties. In particular, we show that, when the automata halt "in finite steps" along all computation paths, they must terminate in worst-case liner time. In the bounded-error probability case, the acceptance of the automata depends only on the computation paths that terminate within exponentially many steps even if not all computation paths may terminate. We also present a classical simulation of those automata on two-way multi-head probabilistic finite automata with cut points. Moreover, we discuss how the recognition power of the automata varies as the automata's acceptance criteria change to error free, one-sided error, bounded error, and unbounded error by comparing the complexity of their computational powers. We further note that, with the use of arbitrary complex transition amplitudes, two-way unbounded-error quantum finite automata and two-way bounded-error 2-head quantum finite automata can recognize certain non-recursive languages, whereas two-way error-free quantum finite automata recognize only recursive languages.


title: Exciting Reflectionless Unidirectional Edge Modes in a Reciprocal Photonic Topological Insulator Medium

abstract: Photonic topological insulators are an interesting class of materials whose photonic band structure can have a bandgap in the bulk while supporting topologically protected unidirectional edge modes. Recent studies [1-6] on bianisotropic metamaterials that emulate the electronic quantum spin Hall effect using its electromagnetic analog are examples of such systems with relatively simple and elegant design. In this paper, we present a rotating magnetic dipole antenna, composed of two perpendicularly oriented coils, that can efficiently excite the unidirectional topologically protected surface waves in the bianisotropic metawaveguide (BMW) structure recently realized by Ma, et al. [1], despite the fact that the BMW medium does not break time-reversal invariance. In addition to achieving high directivity, the antenna can be tuned continuously to excite reflectionless edge modes to the two opposite directions with various amplitude ratios. We demonstrate its performance through experiment and compare to simulation results.


title: Robust Massively Parallel Sorting

abstract: We investigate distributed memory parallel sorting algorithms that scale to the largest available machines and are robust with respect to input size and distribution of the input elements. The main outcome is that four sorting algorithms cover the entire range of possible input sizes. For three algorithms we devise new low overhead mechanisms to make them robust with respect to duplicate keys and skewed input distributions. One of these, designed for medium sized inputs, is a new variant of quicksort with fast high-quality pivot selection. At the same time asymptotic analysis provides performance guarantees and guides the selection and configuration of the algorithms. We validate these hypotheses using extensive experiments on 7 algorithms, 10 input distributions, up to 262144 cores, and varying input sizes over 9 orders of magnitude. For difficult input distributions, our algorithms are the only ones that work at all. For all but the largest input sizes, we are the first to perform experiments on such large machines at all and our algorithms significantly outperform the ones one would conventionally have considered.


title: $D^*$ $\Xi N$ bound state in strange three-body systems

abstract: The recent update of the strangeness $-2$ ESC08c Nijmegen potential incorporating the NAGARA and KISO events predicts a $\Xi N$ bound state, $D^*$, in the $^3S_1 (I=1)$ channel. We study if the existence of this two-body bound state could give rise to stable three-body systems. For this purpose we solve the bound state problem of three-body systems where the $\Xi N$ state is merged with $N$'s, $\Lambda$'s, $\Sigma's$ or $\Xi$'s, making use of the most recent updates of the two-body ESC08c Nijmegen potentials. We found that there appear stable states in the $\Xi NN$ and $\Xi \Xi N$ systems, the $\Xi \Lambda N$ and $\Xi \Sigma N$ systems being unbound.


title: Degree-correlation, robustness, and vulnerability in finite scale-free networks

abstract: Many naturally occurring networks have a power-law degree distribution as well as a non-zero degree correlation. Despite this, most studies analyzing the robustness to random node-deletion and vulnerability to targeted node-deletion have concentrated only on power-law degree distribution and ignored degree correlation. This study looks specifically at the effect degree-correlation has on robustness and vulnerability in scale-free networks. Our results confirm Newman's finding that positive degree-correlation increases robustness and decreases vulnerability. However, we found that networks with positive degree-correlation are more vulnerable to random node-deletion than to targeted deletion methods that utilize knowledge of initial node-degree only. Targeted deletion sufficiently alters the topology of the network to render this method less effective than uniform random methods unless changes in topology are accounted for. This result indicates the importance of degree correlation in certain network applications.


title: A note on the scaling limits of random P\'olya trees

abstract: Panagiotou and Stufler (arXiv:1502.07180v2) recently proved one important fact on their way to establish the scaling limits of random P\'{o}lya trees: a uniform random P\'{o}lya tree of size $n$ consists of a conditioned critical Galton-Watson tree $C_n$ and many small forests, where with probability tending to one as $n$ tends to infinity, any forest $F_n(v)$, that is attached to a node $v$ in $C_n$, is maximally of size $\vert F_n(v)\vert=O(\log n)$. Their proof used the framework of a Boltzmann sampler and deviation inequalities. In this paper, first, we employ a unified framework in analytic combinatorics to prove this fact with additional improvements on the bound of $\vert F_n(v)\vert$, namely $\vert F_n(v)\vert=\Theta(\log n)$. Second, we give a combinatorial interpretation of the rational weights of these forests and the defining substitution process in terms of automorphisms associated to a given P\'{o}lya tree. Finally, we derive the limit probability that for a random node $v$ the attached forest $F_n(v)$ is of a given size.


title: Time-dependent global sensitivity analysis with active subspaces for a lithium ion battery model

abstract: Renewable energy researchers use computer simulation to aid the design of lithium ion storage devices. The underlying models contain several physical input parameters that affect model predictions. Effective design and analysis must understand the sensitivity of model predictions to changes in model parameters, but global sensitivity analyses become increasingly challenging as the number of input parameters increases. Active subspaces are part of an emerging set of tools for discovering and exploiting low-dimensional structures in the map from high-dimensional inputs to model outputs. We extend linear and quadratic model-based heuristic for active sub- space discovery to time-dependent processes and apply the resulting technique to a lithium ion battery model. The results reveal low-dimensional structure and sensitivity metrics that a designer may exploit to study the relationship between parameters and predictions.


title: Characterising the physical and chemical properties of a young Class 0 protostellar core embedded in the Orion B9 filament

abstract: The present study aims to characterise the physical and chemical properties of the protostellar core Orion B9-SMM3. The APEX telescope was used to perform a follow-up molecular line survey of SMM3. The following species were identified from the frequency range 218.2-222.2 GHz: $^{13}$CO, C$^{18}$O, SO, para-H$_2$CO, and E$_1$-type CH$_3$OH. The on-the-fly mapping observations at 215.1-219.1 GHz revealed that SMM3 is associated with a dense gas core as traced by DCO$^+$ and p-H$_2$CO. Altogether three different p-H$_2$CO transitions were detected with clearly broadened linewidths (8.2-11 km s$^{-1}$ in FWHM). The derived p-H$_2$CO rotational temperature, $64\pm15$ K, indicates the presence of warm gas. We also detected a narrow p-H$_2$CO line (FWHM=0.42 km s$^{-1}$) at the systemic velocity. The p-H$_2$CO abundance for the broad component appears to be enhanced by two orders of magnitude with respect to the narrow line value ($\sim3\times10^{-9}$ versus $\sim2\times10^{-11}$). The detected methanol line shows a linewidth similar to those of the broad p-H$_2$CO lines, which indicates their coexistence. The CO isotopologue data suggest that the CO depletion factor decreases from $\sim27\pm2$ towards the core centre to a value of $\sim8\pm1$ towards the core edge. In the latter position, the N$_2$D$^+$/N$_2$H$^+$ ratio is revised down to $0.14\pm0.06$. The origin of the subfragments inside the SMM3 core we found previously can be understood in terms of the Jeans instability if non-thermal motions are taken into account. The estimated fragmentation timescale, and the derived chemical abundances suggest that SMM3 is a few times $10^5$ yr old, in good agreement with its Class 0 classification inferred from the spectral energy distribution analysis. The broad p-H$_2$CO and CH$_3$OH lines, and the associated warm gas provide the first clear evidence of a molecular outflow driven by SMM3.


title: Search for bottom squark pair production in proton--proton collisions at $\sqrt{s}=$13 TeV with the ATLAS detector

abstract: The result of a search for pair production of the supersymmetric partner of the Standard Model bottom quark ($\tilde{b}_1$) is reported. The search uses 3.2 fb$^{-1}$ of $pp$ collisions at $\sqrt{s}=$13 TeV collected by the ATLAS experiment at the Large Hadron Collider in 2015. Bottom squarks are searched for in events containing large missing transverse momentum and exactly two jets identified as originating from $b$-quarks. No excess above the expected Standard Model background yield is observed. Exclusion limits at 95% confidence level on the mass of the bottom squark are derived in phenomenological supersymmetric $R$-parity-conserving models in which the $\tilde{b}_1$ is the lightest squark and is assumed to decay exclusively via $\tilde{b}_1 \rightarrow b \tilde{\chi}_1^0$, where $\tilde{\chi}_1^0$ is the lightest neutralino. The limits significantly extend previous results; bottom squark masses up to 800 (840) GeV are excluded for the $\tilde{\chi}_1^0$ mass below 360 (100) GeV whilst differences in mass above 100 GeV between the $\tilde{b}_1$ and the $\tilde{\chi}_1^0$ are excluded up to a $\tilde{b}_1$ mass of 500 GeV.


title: Measure Algebras on Homogeneous Spaces

abstract: For a locally compact group $G$ and a compact subgroup $H$, we show that the Banach space $M(G/H)$ may be considered as a quotient space of $M(G)$. Also, we define a convolution on $M(G/H)$ which makes it into a Banach algebra. It may be identified with a closed subalgebra of the involutive Banach algebra $M(G)$, and there is no involution on $M(G/H)$ compatible with this identification unless $H$ is a normal subgroup of $G$. In other words, $M(G/H)$ is a $*$-Banach subalgebra of $M(G)$ only if $H$ is a normal subgroup of $G$. As well, it is a unital Banach algebra just when $H$ is a normal subgroup. Furthermore, when $G/H$ is attached to a strongly quasi-invariant measure, $L^1(G/H)$ is a Banach subspace of $M(G/H)$. Using the restriction of the convolution on $M(G/H)$, we obtain a Banach algebra $L^1(G/H)$, which may be considered as a Banach subalgebra of $L^1(G)$, with a right approximate identity. It has no involution and no left approximate identity except for a normal subgroup $H$. Consequently, the Banach algebra $L^1(G/H)$ is amenable if and only if $H$ is a normal subgroup and $G$ is amenable.


title: Control of Spin Helix Symmetry in Semiconductor Quantum Wells by Crystal Orientation

abstract: We investigate the possibility of spin-preserving symmetries due to the interplay of Rashba and Dresselhaus spin-orbit coupling in n-doped zinc-blende semiconductor quantum wells of general crystal orientation. It is shown that a conserved spin operator can be realized if and only if at least two growth-direction Miller indices agree in modulus. The according spin-orbit field has in general both in-plane and out-of-plane components and is always perpendicular to the shift vector of the corresponding persistent spin helix. We also analyze higher-order effects arising from the Dresselhaus term, and the impact of our results on weak (anti)localization corrections.


title: Mass density of the Earth from a Gravito-Electro-Magnetic 5D vacuum

abstract: We calculate the mass density of the Earth using a Gravito-Electro-Magnetic theory on an extended 5D Schwarzschild-de Sitter metric, in which we define the vacuum. Our results are in very good agreement with that of the Dziewonski-Anderson model.


title: Performance of dynamical decoupling in bosonic environments and under pulse-timing fluctuations

abstract: We study the suppression of qubit dephasing through Uhrig dynamical decoupling (UDD) in nontrivial environments modeled within the spin-boson formalism. In particular, we address the case of (i) a qubit coupled to a bosonic bath with power-law spectral density, and (ii) a qubit coupled to a single harmonic oscillator that dissipates energy into a bosonic bath, which embodies an example of a structured bath for the qubit. We then model the influence of random time jitter in the UDD protocol by sorting pulse-application times from Gaussian distributions centered at appropriate values dictated by the optimal protocol. In case (i) we find that, when few pulses are applied and a sharp cutoff is considered, longer coherence times and robust UDD performances (against random timing errors) are achieved for a super-Ohmic bath. On the other hand, when an exponential cutoff is considered a super-Ohmic bath is undesirable. In case (ii) the best scenario is obtained for an overdamped harmonic motion. Our study provides relevant information for the implementation of optimized schemes for the protection of quantum states from decoherence.


title: "Show me the cup": Reference with Continuous Representations

abstract: One of the most basic functions of language is to refer to objects in a shared scene. Modeling reference with continuous representations is challenging because it requires individuation, i.e., tracking and distinguishing an arbitrary number of referents. We introduce a neural network model that, given a definite description and a set of objects represented by natural images, points to the intended object if the expression has a unique referent, or indicates a failure, if it does not. The model, directly trained on reference acts, is competitive with a pipeline manually engineered to perform the same task, both when referents are purely visual, and when they are characterized by a combination of visual and linguistic properties.


title: Could Segue 1 be a destroyed star cluster? - a dynamical perspective

abstract: We attempt to find a progenitor for the ultra-faint object Segue 1 under the assumption that it formed as a dark matter free star cluster in the past. We look for orbits, using the elongation of Segue 1 on the sky as a tracer of its path. Those orbits are followed backwards in time to find the starting points of our N-body simulations. The successful orbit, with which we can reproduce Segue 1 has a proper motion of mu_alpha = -0.19 mas/yr and mu_delta = -1.9 mas/yr, placing Segue 1 near its apo-galacticon today. Our best fitting model has an initial mass of 6224 Msun and an initial scale-length of 5.75 pc.


title: Interface propagation in fiber bundles: Local, mean-field and intermediate range-dependent statistics

abstract: The fiber bundle model is essentially an array of elements that break when sufficient load is applied on them. With a local loading mechanism, this can serve as a model for a one-dimensional interface separating the broken and unbroken parts of a solid in mode-I fracture. The interface can propagate through the system depending on the loading rate and disorder present in the failure thresholds of the fibers. In the presence of a quasi-static drive, the intermittent dynamics of the interface mimic front propagation in disordered media. Such situations appear in diverse physical systems such as mode-I crack propagation, domain wall dynamics in magnets, charge density waves, contact line in wetting etc. We study the effect of the range of interaction, i.e. the neighborhood of the interface affected following a local perturbation, on the statistics of the intermittent dynamics of the front. There exists a crossover from local to global behavior as the range of interaction grows and a continuously varying `universality' in the intermediate range. This means that the interaction range is a relevant parameter of any resulting physics. This is particularly relevant in view of the fact that there is a scatter in the experimental observations of the exponents, in even idealized experiments on fracture fronts, and also a possibility in changing the interaction range in real samples.


title: Comment on "Frustration and Multicriticality in the Antiferromagnetic Spin-1 Chain"

abstract: The phase diagram of the spin-1 chain with bilinear-biquadratic and next-nearest neighbor inter- actions, recently investigated by Pixley, Shashi and Nevidomskyy [Phys. Rev. B 90, 214426 (2014)], has been revisited in the light of results we have recently obtained on a similar model. Combining extensive Density Matrix Renormalization Group (DMRG) simulations with conformal-field theory arguments, we confirm the presence of the three phases identified by Pixley et al, a Haldane phase, a next-nearest neighbor (NNN) Haldane phase, and a dimerized phase, but we come to significantly different conclusions regarding the nature of the phase transitions to the dimerized phase: i) We provide numerical evidence of a continuous Ising transition between the NNN-Haldane phase and the dimerized phase; ii) We show that the tri-critical end point, where the continuous transition between the Haldane phase and the dimerized phase turns into a first order transition, is distinct from the triple point where the three phases meet; iii) Finally, we demonstrate that the tri-critical end point is in the same Wess-Zumino-Witten (WZW) SU(2) level 2 universality class as the continuous transition line that ends at this point


title: Simulating generic spin-boson models with matrix product states

abstract: The global coupling of few-level quantum systems ("spins") to a discrete set of bosonic modes is a key ingredient for many applications in quantum science, including large-scale entanglement generation, quantum simulation of the dynamics of long-range interacting spin models, and hybrid platforms for force and spin sensing. We present a general numerical framework for treating the out-of-equilibrium dynamics of such models based on matrix product states. Our approach applies for generic spin-boson systems: it treats any spatial and operator dependence of the two-body spin-boson coupling and places no restrictions on relative energy scales. We show that the full counting statistics of collective spin measurements and infidelity of quantum simulation due to spin-boson entanglement, both of which are difficult to obtain by other techniques, are readily calculable in our approach. We benchmark our method using a recently developed exact solution for a particular spin-boson coupling relevant to trapped ion quantum simulators. Finally, we show how decoherence can be incorporated within our framework using the method of quantum trajectories, and study the dynamics of an open-system spin-boson model with spatially non-uniform spin-boson coupling relevant for trapped atomic ion crystals in the presence of molecular ion impurities.


title: The time evolution of gaps in tidal streams

abstract: We model the time evolution of gaps in tidal streams caused by the impact of a dark matter subhalo, while both orbit a spherical gravitational potential. To this end, we make use of the simple behaviour of orbits in action-angle space. A gap effectively results from the divergence of two nearby orbits whose initial phase-space separation is, for very cold thin streams, largely given by the impulse induced by the subhalo. We find that in a spherical potential the size of a gap increases linearly with time for sufficiently long timescales. We have derived an analytic expression that shows how the growth rate depends on the mass of the perturbing subhalo, its scale and its relative velocity with respect to the stream. We have verified these scalings using N-body simulations and find excellent agreement. For example, a subhalo of mass 10^8 Msun directly impacting a very cold thin stream on an inclined orbit can induce a gap that may reach a size of several tens of kpc after a few Gyr. The gap size fluctuates importantly with phase on the orbit, and it is largest close to pericentre. This indicates that it may not be fully straightforward to invert the spectrum of gaps present in a stream to recover the mass spectrum of the subhalos.


title: Assessment of Models of Galactic Thermal Dust Emission Using COBE/FIRAS and COBE/DIRBE Observations

abstract: Accurate modeling of the spectrum of thermal dust emission at millimeter wavelengths is important for improving the accuracy of foreground subtraction for CMB measurements, for improving the accuracy with which the contributions of different foreground emission components can be determined, and for improving our understanding of dust composition and dust physics. We fit four models of dust emission to high Galactic latitude COBE/FIRAS and COBE/DIRBE observations from 3 millimeters to 100 microns and compare the quality of the fits. We consider the two-level systems model because it provides a physically motivated explanation for the observed long wavelength flattening of the dust spectrum and the anticorrelation between emissivity index and dust temperature. We consider the model of Finkbeiner, Davis, and Schlegel because it has been widely used for CMB studies, and the generalized version of this model recently applied to Planck data by Meisner and Finkbeiner. For comparison we have also fit a phenomenological model consisting of the sum of two graybody components. We find that the two-graybody model gives the best fit and the FDS model gives a significantly poorer fit than the other models. The Meisner and Finkbeiner model and the two-level systems model remain viable for use in Galactic foreground subtraction, but the FIRAS data do not have sufficient signal-to-noise ratio to provide a strong test of the predicted spectrum at millimeter wavelengths.


title: Nonlocal Infrared Modifications of Gravity. A Review

abstract: We review an approach developed in the last few years by our group in which GR is modified in the infrared, at an effective level, by nonlocal terms associated to a mass scale. We begin by recalling the notion of quantum effective action and its associated nonlocalities, illustrating some of their features with the anomaly-induced effective actions in $D=2$ and $D=4$. We examine conceptual issues of nonlocal theories such as causality, degrees of freedoms and ghosts, stressing the importance of the fact that these nonlocalities only emerge at the effective level. We discuss a particular class of nonlocal theories where the nonlocal operator is associated to a mass scale, and we show that they perform very well in the comparison with cosmological observations, to the extent that they fit CMB, supernovae, BAO and structure formation data at a level fully competitive with $\Lambda$CDM, with the same number of free parameters. We explore some extensions of these `minimal' models, and we finally discuss some directions of investigation for deriving the required effective nonlocality from a fundamental local QFT.


title: Phenomenology of SU(5) low-energy realizations: the diphoton excess and Higgs flavor violation

abstract: We discuss different $\rm SU(5)$ low-energy realizations and illustrate their use with the diphoton excess and Higgs flavor violation, which require new physics at the TeV scale. In particular, we study two scenarios for a $750$ GeV resonance: in the first one the resonance belongs to the adjoint of $\rm SU(5)$, being either an $\rm SU(2)_L$ singlet or a triplet, while in the second case the signal is due to the CP-even and CP-odd states of a new $\rm SU(2)_L$ Higgs doublet belonging to a ${\bf 45}_H$ or a ${\bf 70}_H$ representations, giving rise to a two-Higgs doublet model at low energies. We study the fine-tuning needed for the desired members of the multiplets to be light enough, while having the rest at the GUT scale. In these scenarios, the production and decay into photons of the new resonance are mediated by the leptoquarks (LQ) present in these large $\rm SU(5)$ representations. We analyse the phenomenology of such scenarios, focusing on the most relevant predictions that can help to disentangle the different models, like decays into gauge bosons, Standard Model (SM) fermions and LQs pair production. In the case of the ${\bf 45}_H$ (the Georgi-Jarlskog model), we also study the possibility to have Higgs flavor violation. We find that $B_s$ mixing limits (in addition to $\tau\rightarrow \mu \gamma$) always imply that $\rm{BR}(h\rightarrow \tau\mu,\,bs)\lesssim 10^{-5}$.