title: Stacked star formation rate profiles of bursty galaxies exhibit 'coherent' star formation

abstract: In a recent work based on 3200 stacked H$\alpha$ maps of galaxies at $z \sim 1$, Nelson et al.~find evidence for `coherent star formation': the stacked SFR profiles of galaxies above (below) the 'star formation main sequence' (MS) are above (below) that of galaxies on the MS at all radii. One might interpret this result as inconsistent with highly bursty star formation and evidence that galaxies evolve smoothly along the MS rather than crossing it many times. We analyze six simulated galaxies at $z\sim1$ from the Feedback in Realistic Environments (FIRE) project in a manner analogous to the observations to test whether the above interpretations are correct. The trends in stacked SFR profiles are qualitatively consistent with those observed. However, SFR profiles of individual galaxies are much more complex than the stacked profiles: the former can be flat or even peak at large radii because of the highly clustered nature of star formation in the simulations. Moreover, the SFR profiles of individual galaxies above (below) the MS are not systematically above (below) those of MS galaxies at all radii. We conclude that the time-averaged coherent star formation evident stacks of observed galaxies is consistent with highly bursty, clumpy star formation of individual galaxies and is not evidence that galaxies evolve smoothly along the MS.


title: Domain walls in the extensions of the Standard Model

abstract: Our main interest is the evolution of domain walls of the Higgs field in the early Universe. The aim of this paper is to understand how dynamics of Higgs domain walls could be influenced by yet unknown interactions from beyond the Standard Model. We assume that the Standard Model is valid up to certain, high, energy scale $\Lambda$ and use the framework of the effective field theory to describe physics below that scale. Performing numerical simulations with different values of the scale $\Lambda$ we are able to extend our previous analysis and determine its range of validity. We study domain walls interpolating between the physical electroweak vacuum and the vacuum appearing at very high field strengths. These domain walls could be formed from non-homogeneous configurations of the Higgs field produced by quantum fluctuations during inflation or thermal fluctuations during reheating. Our numerical simulations show that evolution of Higgs domain walls is rather insensitive to interactions beyond the Standard Model as long as masses of new particles are grater than $10^{12}\ \textrm{GeV}$. For lower values of $\Lambda$ the RG improved effective potential is strongly modified at field strengths crucial to the evolution of domain walls. For instance its minima become degenerate for $\Lambda$ around $10^{11}\ \textrm{GeV}$. We find that even in the case when the minima of the potential are nearly degenerate Higgs domain walls decayed shortly after their formation for generic initial conditions. On the other hand, in simulations with specifically chosen initial conditions Higgs domain walls can live longer and enter the scaling regime. We also determine the energy spectrum of gravitational waves produced by decaying domain walls of the Higgs field. For generic initial field configurations the amplitude of the signal is too small to be observed in present and planned detectors.


title: From Euclidean Sources to Lorentzian Spacetimes in Holographic Conformal Field Theories

abstract: We consider states of holographic conformal field theories constructed by adding sources for local operators in the Euclidean path integral, with the aim of investigating the extent to which arbitrary bulk coherent states can be represented by such Euclidean path-integrals in the CFT. We construct the associated dual Lorentzian spacetimes perturbatively in the sources. Extending earlier work, we provide explicit formulae for the Lorentzian fields to first order in the sources for general scalar field and metric perturbations in arbitrary dimensions. We check the results by holographically computing the Lorentzian one-point functions for the sourced operators and comparing with a direct CFT calculation. We present evidence that at the linearized level, arbitrary bulk initial data profiles can be generated by an appropriate choice of Euclidean sources. However, in order to produce initial data that is very localized, the amplitude must be taken small at the same time otherwise the required sources diverge, invalidating the perturbative approach.


title: Axisymmetric spheroidal modes of neutron stars magnetized with poloidal magnetic fields

abstract: We calculate axisymmetric spheroidal modes of neutron stars magnetized with an axisymmetric poloidal magnetic field. We use polytropes of the indices $n\sim1$ as background equilibrium models of neutron stars where we ignore the deformation due to the magnetic fields. For a poloidal magnetic field, axisymmetric normal modes of non-rotating stars are decoupled into spheroidal modes and toroidal modes, and we can treat spheroidal modes separately from toroidal modes. For the surface field strength $B_S$ ranging from $10^{14}$G to $10^{16}$G, we calculate axisymmetric spheroidal magnetic modes whose oscillation frequency is proportional to $B_S$. The typical oscillation frequency of the magnetic modes is $\sim 0.01\times\sqrt{GM/R^3}$ for $B_S\sim 10^{15}$G, where $M$ and $R$ are respectively the mass and radius of the star and $G$ is the gravitational constant. For $M=1.4M_\odot$ and $R=10^6$cm, this frequency is $\sim 20$Hz, which may explain low frequency QPOs found for SGR 1806-204 and SGR 1900+14. We also find modes of frequency $> \sqrt{GM/R^3}$ corresponding to the radial fundamental and first harmonic modes. No unstable magnetic modes are found for axisymmetric spheroidal oscillations of magnetized stars.


title: Infinitesimal symmetries for implicit planar webs

abstract: In the complex setting, let $F(x,y,y')=0$ be an analytic or algebraic differential equation with $y'$-degree $d$. We deal with the qualitative study of such equations through the geometry of the planar $d$-web generated by the generic family of integral curves. Infinitesimal symmetries of these configurations are discussed by using meromorphic connections methods "\`a la Cartan-Spencer".


title: Buses of Cuernavaca - an agent-based model for universal random matrix behavior minimizing mutual information

abstract: The public transportation system of Cuernavaca, Mexico, exhibits random matrix theory statistics [1]. In particular, the fluctuation of times between the arrival of buses on a given bus stop, follows the Wigner surmise for the Gaussian Unitary Ensemble. To model this, we propose an agent-based approach in which each bus driver tries to optimize his arrival time to the next stop with respect to an estimated arrival time of his predecessor. We choose a particular form of the associated utility function and recover the appropriate distribution in numerical experiments for a certain value of the only parameter of the model. We then investigate whether this value of the parameter is otherwise distinguished within an information theoretic approach and give numerical evidence that indeed it is associated with a minimum of averaged pairwise mutual information.


title: Decay of a Thermofield-Double State in Chaotic Quantum Systems

abstract: Scrambling in interacting quantum systems out of equilibrium is particularly effective in the chaotic regime. Under time evolution, initially localized information is said to be scrambled as it spreads throughout the entire system. This spreading can be analyzed with the spectral form factor, which is defined in terms of the analytic continuation of the partition function. The latter is equivalent to the survival probability of a thermofield double state under unitary dynamics. Using random matrices from the Gaussian unitary ensemble (GUE) as Hamiltonians for the time evolution, we obtain exact analytical expressions at finite $N$ for the survival probability. Numerical simulations of the survival probability with matrices taken from the Gaussian orthogonal ensemble (GOE) are also provided. The GOE is more suitable for our comparison with numerical results obtained with a disordered spin chain with local interactions. Common features between the random matrix and the realistic disordered model in the chaotic regime are identified. The differences that emerge as the spin model approaches a many-body localized phase are also discussed.


title: What is the Machine Learning?

abstract: Applications of machine learning tools to problems of physical interest are often criticized for producing sensitivity at the expense of transparency. To address this concern, we explore a data planing procedure for identifying combinations of variables -- aided by physical intuition -- that can discriminate signal from background. Weights are introduced to smooth away the features in a given variable(s). New networks are then trained on this modified data. Observed decreases in sensitivity diagnose the variable's discriminating power. Planing also allows the investigation of the linear versus non-linear nature of the boundaries between signal and background. We demonstrate the efficacy of this approach using a toy example, followed by an application to an idealized heavy resonance scenario at the Large Hadron Collider. By unpacking the information being utilized by these algorithms, this method puts in context what it means for a machine to learn.


title: No Sun-like dynamo on the active star $\zeta$ Andromedae from starspot asymmetry

abstract: Sunspots are cool areas caused by strong surface magnetic fields inhibiting convection. Moreover, strong magnetic fields can alter the average atmospheric structure, degrading our ability to measure stellar masses and ages. Stars more active than the Sun have more and stronger dark spots than in the solar case, including on the rotational pole itself. Doppler imaging, which has so far produced the most detailed images of surface structures on other stars than the Sun, cannot always distinguish the hemisphere in which the starspots are located, especially in the equatorial region and if the data quality is not optimal. This leads to problems in investigating the north-south distribution of starspot active latitudes (those latitudes with more spot activity), which are crucial constraints of dynamo theory. Polar spots, inferred only from Doppler tomography, could plausibly be observational artifacts, casting some doubt on their very existence. Here we report imaging of the old, magnetically-active star $\zeta$ Andromedae using long-baseline infrared interferometry. In our data, a dark polar spot is seen in each of two epochs, while lower-latitude spot structures in both hemispheres do not persist between observations revealing global starspot asymmetries. The north-south symmetry of active latitudes observed on the Sun is absent on $\zeta$ And, which hosts global spot patterns that cannot be produced by solar-type dynamos.


title: Galaxy Cluster Mass Reconstruction Project: III. The impact of dynamical substructure on cluster mass estimates

abstract: With the advent of wide-field cosmological surveys, we are approaching samples of hundreds of thousands of galaxy clusters. While such large numbers will help reduce statistical uncertainties, the control of systematics in cluster masses becomes ever more crucial. Here we examine the effects of an important source of systematic uncertainty in galaxy-based cluster mass estimation techniques: the presence of significant dynamical substructure. Dynamical substructure manifests as dynamically distinct subgroups in phase-space, indicating an 'unrelaxed' state. This issue affects around a quarter of clusters in a generally selected sample. We employ a set of mock clusters whose masses have been measured homogeneously with commonly-used galaxy-based mass estimation techniques (kinematic, richness, caustic, radial methods). We use these to study how the relation between observationally estimated and true cluster mass depends on the presence of substructure, as identified by various popular diagnostics. We find that the scatter for an ensemble of clusters does not increase dramatically for clusters with dynamical substructure. However, we find a systematic bias for all methods, such that clusters with significant substructure have higher measured masses than their relaxed counterparts. This bias depends on cluster mass: the most massive clusters are largely unaffected by the presence of significant substructure, but masses are significantly overestimated for lower mass clusters, by $\sim10\%$ at $10^{14}$ and $\geq20\%$ for $\leq10^{13.5}$. The use of cluster samples with different levels of substructure can, therefore, bias certain cosmological parameters up to a level comparable to the typical uncertainties in current cosmological studies.


title: Contemporaneous Imaging Comparisons of the Spotted Giant $\sigma$ Geminorum Using Interferometric, Spectroscopic, and Photometric Data

abstract: Nearby, active stars with relatively rapid rotation and large starspot structures offer the opportunity to compare interferometric, spectroscopic, and photometric imaging techniques. In this paper, we image a spotted star with three different methods for the first time. The giant primary star of the RS Canum Venaticorum binary $\sigma$ Geminorum ($\sigma$ Gem) was imaged for two epochs of interferometric, high-resolution spectroscopic, and photometric observations. The light curves from the reconstructions show good agreement with the observed light curves, supported by the longitudinally-consistent spot features on the different maps. However, there is strong disagreement in the spot latitudes across the methods.


title: The extragalactic ultra-high energy cosmic-ray dipole

abstract: We explore the possibility that the recently detected dipole anisotropy in the arrival directions of~$>8$~EeV ultra-high energy cosmic-rays (UHECRs) arises due to the large-scale structure (LSS). We assume that the cosmic ray sources follow the matter distribution and calculate the flux-weighted UHECRs' RMS dipole amplitude taking into account the diffusive transport in the intergalactic magnetic field (IGMF). We find that the flux-weighted RMS dipole amplitude is $\sim8$% before entering the Galaxy. The amplitude in the [4-8] EeV is only slightly lower $\sim 5$%. The required IGMF is of the order of {5-30 nG}, and the UHECR sources must be relatively nearby, within $\sim$300 Mpc. The absence of statistically significant signal in the lower energy bin can be explained if the same nuclei specie dominates the composition in both energy bins and diffusion in the Galactic magnetic field (GMF) reduces the dipole of these lower rigidity particles. Photodisintegration of higher energy UHECRs could also reduce somewhat the lower energy dipole.


title: Probing the Metal Enrichment of the Intergalactic Medium at $z=5-6$ Using the Hubble Space Telescope

abstract: We test the galactic outflow model by probing associated galaxies of four strong intergalactic CIV absorbers at $z=5$--6 using the Hubble Space Telescope (HST) ACS ramp narrowband filters. The four strong CIV absorbers reside at $z=5.74$, $5.52$, $4.95$, and $4.87$, with column densities ranging from $N_{\rm{CIV}}=10^{13.8}$ cm$^{-2}$ to $10^{14.8}$ cm$^{-2}$. At $z=5.74$, we detect an i-dropout Ly$\alpha$ emitter (LAE) candidate with a projected impact parameter of 42 physical kpc from the CIV absorber. This LAE candidate has a Ly$\alpha$-based star formation rate (SFR$_{\rm{Ly\alpha}}$) of 2 $M_\odot$ yr$^{-1}$ and a UV-based SFR of 4 $M_\odot$ yr$^{-1}$. Although we cannot completely rule out that this $i$-dropout emitter may be an [OII] interloper, its measured properties are consistent with the CIV powering galaxy at $z=5.74$. For CIV absorbers at $z=4.95$ and $z=4.87$, although we detect two LAE candidates with impact parameters of 160 kpc and 200 kpc, such distances are larger than that predicted from the simulations. Therefore we treat them as non-detections. For the system at $z=5.52$, we do not detect LAE candidates, placing a 3-$\sigma$ upper limit of SFR$_{\rm{Ly\alpha}}\approx 1.5\ M_\odot$ yr$^{-1}$. In summary, in these four cases, we only detect one plausible CIV source at $z=5.74$. Combining the modest SFR of the one detection and the three non-detections, our HST observations strongly support that smaller galaxies (SFR$_{\rm{Ly\alpha}} \lesssim 2\ M_\odot$ yr$^{-1}$) are main sources of intergalactic CIV absorbers, and such small galaxies play a major role in the metal enrichment of the intergalactic medium at $z\gtrsim5$.


title: The quantum spin quadrumer

abstract: A fundamental motif in frustrated magnetism is the fully mutually coupled cluster of $N$ spins, with each spin coupled to every other spin. Clusters with $N=2$ and $3$ have been extensively studied as building blocks of square and triangular lattice antiferromagnets. In both cases, large-$S$ semiclassical descriptions have been fruitfully constructed, providing insights into the physics of macroscopic magnetic systems. Here, we develop a semiclassical theory for the $N=4$ cluster. This problem has rich mathematical structure with a ground state space that has non-trivial topology. We show that the ground states are appropriately parametrized by a unit vector order parameter and a rotation matrix. Remarkably, in the low energy description, the physics of the cluster reduces to that of an emergent free spin-$S$ spin and a rigid rotor. This successfully explains the spectrum of the quadrumer and its associated degeneracies. However, this mapping does not hold in the vicinity of collinear ground states due to a subtle effect that arises from the non-manifold nature of the ground state space. We demonstrate this by an analysis of soft fluctuations, showing that collinear states have a larger number of soft modes. Nevertheless, as these singularities only occur on a subset of measure zero, the mapping to a spin and a rotor provides a good description of the quadrumer. We interpret thermodynamic properties of the quadrumer that are accessible in molecular magnets, in terms of the rotor and spin degrees of freedom. Our study paves the way for field theoretic descriptions of systems such as pyrochlore magnets.


title: Angular Momentum Transport in Thin Magnetically Arrested Disks

abstract: In accretion disks with large-scale ordered magnetic fields, the magnetorotational instability (MRI) is marginally suppressed, so other processes may drive angular momentum transport leading to accretion. Accretion could then be driven by large-scale magnetic fields via magnetic braking, but large-scale magnetic flux can build-up onto the black hole and within the disk leading to a magnetically-arrested disk (MAD). Such a MAD state is unstable to the magnetic Rayleigh-Taylor (RT) instability, which itself leads to vigorous turbulence and the emergence of low-density highly-magnetized bubbles. This instability was studied in a thin (ratio of half-height H to radius R, $H/R \approx 0.1$) MAD simulation, where it has a more dramatic effect on the dynamics of the disk than for thicker disks. We find that the low-density bubbles created by the magnetic RT instability decrease the stress (leading to angular momentum transport) in the disk rather than increasing magnetic torques. Indeed, we find that the dominant component of the stress is due to turbulent magnetic fields, despite the suppression of the axisymmetric MRI and the dominant presence of large-scale magnetic fields. This suggests that the magnetic RT instability plays a significant role in driving angular momentum transport in MADs.


title: The multiwavelength spectrum of NGC 3115: Hot accretion flow properties

abstract: NGC 3115 is the nearest galaxy hosting a billion solar mass black hole and is also a low-luminosity active galactic nucleus (LLAGN). X-ray observations of this LLAGN are able to spatially resolve the hot gas within the sphere of gravitational influence of the supermassive black hole. These observations make NGC 3115 an important testbed for black hole accretion theory in galactic nuclei since they constrain the outer boundary conditions of the hot accretion flow. We present a compilation of the multiwavelength spectral energy distribution (SED) of the nucleus of NGC 3115 from radio to X-rays. We report the results from modeling the observed SED with radiatively inefficient accretion flow (RIAF) models. The radio emission can be well-explained by synchrotron emission from the RIAF without the need for contribution from a relativistic jet. We obtain a tight constraint on the RIAF density profile, $\rho(r) \propto r^{-0.73_{-0.02} ^{+0.01}}$, implying that mass-loss through subrelativistic outflows from the RIAF is significant. The lower frequency radio observation requires the synchrotron emission from a nonthermal electron population in the RIAF, similarly to Sgr A*.


title: How Glassy Relaxation Slows Down by Increasing Mobility

abstract: We investigate how structural relaxation in mixtures with strong dynamical asymmetry is affected by the microscopic dynamics. Brownian and Newtonian dynamics simulations of dense mixtures of fast and slow hard spheres reveal a striking trend reversal. Below a critical density, increasing the mobility of the fast particles fluidizes the system, yet, above that critical density, the same increase in mobility strongly hinders the relaxation of the slow particles. The critical density itself does not depend on the dynamical asymmetry and can be identified with the glass-transition density of the mode-coupling theory. The asymptotic dynamics close to the critical density is universal, but strong pre-asymptotic effects prevail in mixtures with additional size asymmetry. This observation reconciles earlier findings of a strong dependence on kinetic parameters of glassy dynamics in colloid--polymer mixtures with the paradigm that the glass transition is determined by the properties of configuration space alone.


title: Flow-Sensitive Composition of Thread-Modular Abstract Interpretation

abstract: We propose a constraint-based flow-sensitive static analysis for concurrent programs by iteratively composing thread-modular abstract interpreters via the use of a system of lightweight constraints. Our method is compositional in that it first applies sequential abstract interpreters to individual threads and then composes their results. It is flow-sensitive in that the causality ordering of interferences (flow of data from global writes to reads) is modeled by a system of constraints. These interference constraints are lightweight since they only refer to the execution order of program statements as opposed to their numerical properties: they can be decided efficiently using an off-the-shelf Datalog engine. Our new method has the advantage of being more accurate than existing, flow-insensitive, static analyzers while remaining scalable and providing the expected soundness and termination guarantees even for programs with unbounded data. We implemented our method and evaluated it on a large number of benchmarks, demonstrating its effectiveness at increasing the accuracy of thread-modular abstract interpretation.


title: Pinching parameters for open (super) strings

abstract: We present an approach to the parametrization of (super) Schottky space obtained by sewing together three-punctured discs with strips. Different cubic ribbon graphs classify distinct sets of pinching parameters; we show how they are mapped onto each other. The parametrization is particularly well-suited to describing the region within (super) moduli space where open bosonic or Neveu-Schwarz string propagators become very long and thin, which dominates the IR behaviour of string theories. We show how worldsheet objects such as the Green's function converge to graph theoretic objects such as the Symanzik polynomials in the $\alpha ' \to 0$ limit, allowing us to see how string theory reproduces the sum over Feynman graphs. The (super) string measure takes on a simple and elegant form when expressed in terms of these parameters.


title: Methyl isocyanate CH3NCO: An important missing organic in current astrochemical networks

abstract: Methyl isocyanate (CH3NCO) is one of the important complex organic molecules detected on the comet 67P/Churyumov-Gerasimenko by Rosetta's Philae lander. It was also detected in hot cores around high-mass protostars along with a recent detection in the solar-type protostar IRAS 16293-2422. We propose here a gas-grain chemical model to form CH3NCO after reviewing various formation pathways with quantum chemical computations. We have used NAUTILUS 3-phase gas-grain chemical model to compare observed abundances in the IRAS 16293-2422. Our chemical model clearly indicates the ice phase origin of CH3NCO.


title: Distributed Join-the-Idle-Queue for Low Latency Cloud Services

abstract: Low latency is highly desirable for cloud services. To achieve low response time, stringent timing requirements are needed for task scheduling in a large-scale server farm spanning thousands of servers. In this paper, we conduct an in-depth analysis for distributed Join-the-Idle-Queue (JIQ), a promising new approximation of an idealized task-scheduling algorithm. In particular, we derive semi-closed form expressions for the delay performance of distributed JIQ, and we propose a new variant of distributed JIQ that offers clear advantages over alternative algorithms for large systems.


title: VENu: The Virtual Environment for Neutrinos

abstract: The Virtual Environment for Neutrinos (VENu) is a virtual reality-based visualisation of the MicroBooNE detector. MicroBooNE is a liquid-argon-based neutrino experiment, which is currently operating in Fermilab's Booster neutrino beam. The new VENu smartphone app provides informative explanations about neutrinos and uses real MicroBooNE neutrino data that can be visualised inside a virtual representation of the MicroBooNE detector. Available for both iOS and Android, the VENu app can be downloaded for free from the Apple and Google marketplaces. The app enables users to immerse themselves inside the MicroBooNE particle detector and to see particle tracks inside. This can be done in Virtual Reality mode, where the users can pair their smartphone with any consumer virtual reality headset and see the detector in 3D. To encourage learning in a fun environment, a game is also available, guiding users to learn about neutrinos and how to detect them. They can also try to "catch"' neutrinos themselves in 3D mode. The app is currently being pursued for a QuarkNet neutrino master class and outreach events at several universities and labs worldwide.

categories: physics.pop-ph cs.HC hep-ex physics.ins-det


title: Electromagnetic energy stored in inhomogeneous scattering systems

abstract: We analytically study the time-averaged electromagnetic energy stored inside scatterers containing inclusions of arbitrary shapes. Assuming the low density of inclusions, we derive the expression for the energy-transport velocity through disordered media without relying on the radiative transfer equation. Moreover, this expression is independent of the shape of scatterers. In addition, we obtain a relation between the dwell and absorption times associated with inclusions by considering the relationship between internal energy and absorption cross-section. An approximation for the electromagnetic energy stored inside a disordered medium in terms of the transport mean free path and the packing fraction is also derived. This expression suggests that the enhanced electromagnetic energy within the host medium is achieved for inclusions exhibiting negative scattering asymmetry parameters. As a result, disordered media with enhanced backscattering is expected to exhibit large quality factors.


title: Mechanism Design for Demand Response Programs with financial and non-monetary (social) Incentives

abstract: Most demand management approaches with non-mandatory policies assume full users' cooperation, which may not be the case given users' beliefs, needs and preferences. In this paper we propose a mechanism for demand management including incentives both with and without money. The mechanism is validated by means of simulation, modeling the consumers as a finite multiagent system which evolves until a stable state, and social incentives diffusion using opinion dynamics.


title: Dynamical boundary conditions in a non-cylindrical domain for the Laplace equation

abstract: In this paper, we study existence, uniqueness and asymptotic behavior of the Laplace equation with dynamical boundary conditions on regular non-cylindrical domains. We write the problem as a non-autonomous Dirichlet-to-Neumann operator and use form methods in a more general framework to accomplish our goal. A class of non-autonomous elliptic problems with dynamical boundary conditions on Lipschitz domains is also considered in this same context.


title: Monogamy of Quantum Privacy

abstract: Quantum mechanics ensures that the information stored in a quantum state is secure and the ability to send private information through a quantum channel is at least as great as the coherent information. We derive trade-off relations between quantum privacy, information gain by Eve and the disturbance caused by Eve to the quantum state that is being sent through a noisy channel. For tripartite quantum states, we show that monogamy of privacy exists in the case of a single sender and multiple receivers. When Alice prepares a tripartite entangled state and shares it with Bob and Charlie through two different noisy quantum channels, we prove that if the minimally guaranteed quantum privacy between Alice and Bob is positive, then the privacy of information between Alice and Charlie has to be negative. Thus, quantum privacy for more than two parties respects mutual exclusiveness. Then, we prove a monogamy relation for the minimally guaranteed quantum privacy for tripartite systems. We also prove a trade-off relation between the entanglement of formation across one partition and the quantum privacy along another partition. Our results show that quantum privacy cannot be freely shared among multiple parties and can have implication in future quantum networks.


title: Post-processing of the HST STIS coronagraphic observations

abstract: In the past 20 years, the Hubble Space Telescope (HST) STIS coronagraphic instrument has observed more than 100 stars, obtaining more than 4,000 readouts since its installment on HST in 1997 and the numbers are still increasing. We reduce the whole STIS coronagraphic archive at the most commonly observed positions (Wedge A0.6 and A1.0) with new post-processing methods, and present our results here. We are able to recover all of the 32 previously reported circumstellar disks, and obtain better contrast close to the star. For some of the disks, our results are limited by the over subtraction of the methods, and therefore the major regions of the disks can be recovered except the faintest regions. We also explain our efforts in the calibration of its new BAR5 occulting position, enabling STIS to explore inner regions as close as 0.2".


title: Bifurcations of a Van der Pol oscillator in a double well

abstract: We study the changes in the phase portrait of a Van der Pol oscillator in a double well with the change in the shape of the well. Both the symmetric and asymmetric double wells are considered. We find that the shape of the well plays a very important role in the dynamics of the oscillator such as with the change in shape of the well, the stability of the fixed points of the system changes as well as limit cycles appear and are being destroyed.


title: Topological Phase Transitions from Harper to Fibonacci Crystals

abstract: Topological properties of Harper and generalized Fibonacci chains are studied in crystalline cases, i.e., for rational values of the modulation frequency. The Harper and Fibonacci crystals at fixed frequency are connected by an interpolating one-parameter Hamiltonian. As the parameter is varied, one observes topological phase transitions, i.e., changes in the Chern integers of two bands due to the degeneracy of these bands at some parameter value. For small frequency, corresponding to a semiclassical regime, the degeneracies are shown to occur when the average energy of the two bands is approximately equal to the energy of the classical separatrix. Spectral and topological features of the Fibonacci crystal for small frequency leave a clear imprint on the corresponding Hofstadter butterfly for arbitrary frequency.


title: Online Learning with Randomized Feedback Graphs for Optimal PUE Attacks in Cognitive Radio Networks

abstract: In a cognitive radio network, a secondary user learns the spectrum environment and dynamically accesses the channel where the primary user is inactive. At the same time, a primary user emulation (PUE) attacker can send falsified primary user signals and prevent the secondary user from utilizing the available channel. The best attacking strategies that an attacker can apply have not been well studied. In this paper, for the first time, we study optimal PUE attack strategies by formulating an online learning problem where the attacker needs to dynamically decide the attacking channel in each time slot based on its attacking experience. The challenge in our problem is that since the PUE attack happens in the spectrum sensing phase, the attacker cannot observe the reward on the attacked channel. To address this challenge, we utilize the attacker's observation capability. We propose online learning-based attacking strategies based on the attacker's observation capabilities. Through our analysis, we show that with no observation within the attacking slot, the attacker loses on the regret order, and with the observation of at least one channel, there is a significant improvement on the attacking performance. Observation of multiple channels does not give additional benefit to the attacker (only a constant scaling) though it gives insight on the number of observations required to achieve the minimum constant factor. Our proposed algorithms are optimal in the sense that their regret upper bounds match their corresponding regret lower-bounds. We show consistency between simulation and analytical results under various system parameters.


title: ALMA 1.3 Millimeter Map of the HD 95086 System

abstract: Planets and minor bodies such as asteroids, Kuiper-belt objects and comets are integral components of a planetary system. Interactions among them leave clues about the formation process of a planetary system. The signature of such interactions is most prominent through observations of its debris disk at millimeter wavelengths where emission is dominated by the population of large grains that stay close to their parent bodies. Here we present ALMA 1.3 mm observations of HD 95086, a young early-type star that hosts a directly imaged giant planet b and a massive debris disk with both asteroid- and Kuiper-belt analogs. The location of the Kuiper-belt analog is resolved for the first time. The system can be depicted as a broad ($\Delta R/R \sim$0.84), inclined (30\arcdeg$\pm$3\arcdeg) ring with millimeter emission peaked at 200$\pm$6 au from the star. The 1.3 mm disk emission is consistent with a broad disk with sharp boundaries from 106$\pm$6 to 320$\pm$20 au with a surface density distribution described by a power law with an index of --0.5$\pm$0.2. Our deep ALMA map also reveals a bright source located near the edge of the ring, whose brightness at 1.3 mm and potential spectral energy distribution are consistent with it being a luminous star-forming galaxy at high redshift. We set constraints on the orbital properties of planet b assuming co-planarity with the observed disk.


title: Inside-Out Planet Formation. IV. Pebble Evolution and Planet Formation Timescales

abstract: Systems with tightly-packed inner planets (STIPs) are very common. Chatterjee & Tan proposed Inside-Out Planet Formation (IOPF), an in situ formation theory, to explain these planets. IOPF involves sequential planet formation from pebble-rich rings that are fed from the outer disk and trapped at the pressure maximum associated with the dead zone inner boundary (DZIB). Planet masses are set by their ability to open a gap and cause the DZIB to retreat outwards. We present models for the disk density and temperature structures that are relevant to the conditions of IOPF. For a wide range of DZIB conditions, we evaluate the gap opening masses of planets in these disks that are expected to lead to truncation of pebble accretion onto the forming planet. We then consider the evolution of dust and pebbles in the disk, estimating that pebbles typically grow to sizes of a few cm during their radial drift from several tens of AU to the inner, $\lesssim1\:$AU-scale disk. A large fraction of the accretion flux of solids is expected to be in such pebbles. This allows us to estimate the timescales for individual planet formation and entire planetary system formation in the IOPF scenario. We find that to produce realistic STIPs within reasonable timescales similar to disk lifetimes requires disk accretion rates of $\sim10^{-9}\:M_\odot\:{\rm yr}^{-1}$ and relatively low viscosity conditions in the DZIB region, i.e., Shakura-Sunyaev parameter of $\alpha\sim10^{-4}$.


title: Laboratory Investigations of Titan Haze Formation: In Situ Measurement of Gas and Particle Composition

abstract: Prior to the arrival of Cassini-Huygens, aerosol production in Titan's atmosphere was believed to begin in the stratosphere where chemical processes are predominantly initiated by FUV radiation. However, measurements taken by Cassini UVIS and CAPS indicate that haze formation initiates in the thermosphere where there is a greater flux of EUV photons and energetic particles available to initiate chemical reactions, including the destruction of N2. The discovery of previously unpredicted nitrogen species in Titan's atmosphere by Cassini INMS indicates that nitrogen participates in the chemistry to a much greater extent than was appreciated before Cassini. The degree of nitrogen in the haze is important for understanding the diversity of molecules present in Titan's atmosphere and on its surface. We have conducted a series of simulation experiments using either spark discharge or FUV photons to initiate chemistry in CH4/N2 gas mixtures (0.01% CH4/99.99% N2 to 10% CH4/90% N2). We obtained in situ real-time measurements using an HR-ToF-AMS to measure the particle composition as a function of particle size and a PIT-MS to measure the composition of gas phase products. These two techniques allow us to investigate the effect of energy source and initial CH4 concentration on the degree of nitrogen incorporation in both the gas and solid phase products. The results presented here confirm that FUV photons produce not only solid phase nitrogen bearing products but also gas phase nitrogen species. We find that in both the gas and solid phase, nitrogen is found in nitriles rather than amines and both the gas phase and solid phase products are composed primarily of molecules with a low degree of aromaticity. The UV experiments reproduce the absolute abundances measured in Titan's stratosphere for a number of gas phase species including C4H2, C6H6, HCN, CH3CN, HC3N, and C2H5CN. [Edited for length]


title: A mixed-integer branching approach for very small formulations of disjunctive constraints

abstract: An important problem in optimization is the construction of mixed-integer programming (MIP) formulations of disjunctive constraints that are both strong and small. Motivated by lower bounds on the number of integer variables that are required by traditional MIP formulations, we present a more general mixed-integer branching formulation framework. Our approach maintains favorable algorithmic properties of traditional MIP formulations: in particular, amenability to branch-and-bound and branch-and-cut algorithms. Our main technical result gives an explicit linear inequality description for both traditional MIP and mixed-integer branching formulations for a wide range of disjunctive constraints. The formulations obtained from this description have linear programming relaxations that are as strong as possible and generalize some of the most computationally effective formulations for piecewise linear functions and other disjunctive constraints. We use this result to produce a strong mixed-integer branching formulation for any disjunctive constraint that uses only two integer variables and a linear number of extra constraints. We sharpen this result for univariate piecewise linear functions and annulus constraints arising in power systems and robotics, producing strong mixed-integer branching formulations that use only two integer variables and a constant ($\leq 6$) number of general inequality constraints. Along the way, we produce two strong logarithmic-sized traditional MIP formulations for the annulus constraint using our main technical result, illustrating its broader utility in the traditional MIP setting.


title: Mirror Anomaly in Dirac Semimetals

abstract: We demonstrate that, apart from the chiral anomaly, Dirac semimetals possess another quantum anomaly, which we call the mirror anomaly, and which manifests in a singular response of the Dirac semimetal to an applied magnetic field. Namely, the anomalous Hall conductivity exhibits step-function singularities when the field is rotated. We show that this phenomenon is closely analogous to the parity anomaly of (2+1)-dimensional Dirac fermions, but with mirror symmetry, which we demonstrate emerges near any Dirac point at a time reversal invariant momentum, replacing the parity symmetry.


title: Closing in on the large-scale CMB power asymmetry

abstract: Measurements of the cosmic microwave background (CMB) temperature anisotropies have revealed a dipolar asymmetry in power at the largest scales, in apparent contradiction with the statistical isotropy of standard cosmological models. The significance of the effect is not very high, and is dependent on a posteriori choices. Nevertheless, a number of models have been proposed that produce a scale-dependent asymmetry. We confront several such models for a physical, position-space modulation with CMB temperature observations. We find that, while some models that maintain the standard isotropic power spectrum are allowed, others, such as those with modulated tensor or uncorrelated isocurvature modes, can be ruled out on the basis of the overproduction of isotropic power. This remains the case even when an extra isocurvature mode fully anti-correlated with the adiabatic perturbations is added to suppress power on large scales.


title: Quantum Mechanics and Hidden Superconformal Symmetry

abstract: Solvability of the ubiquitous quantum harmonic oscillator relies on a spectrum generating osp(1|2) superconformal symmetry. We study the problem of constructing all quantum mechanical models with a hidden osp(1|2) symmetry on a given space of states. This problem stems from interacting higher spin models coupled to gravity. In one dimension, we show that the solution to this problem is the Plyushchay family of quantum mechanical models with hidden superconformal symmetry obtained by viewing the harmonic oscillator as a one dimensional Dirac system, so that Grassmann parity equals wavefunction parity. These models--both oscillator and particle-like--realize all possible unitary irreducible representations of osp(1|2).


title: Collective sensing of $\beta$-cells generates the metabolic code at optimal islet size

abstract: Major part of a pancreatic islet is composed of beta cells that secrete insulin, a key hormone regulating influx of nutrients into all cells in a vertebrate organism to support nutrition, housekeeping or energy storage. Beta cells constantly communicate with each other using both direct, short-range interactions through gap junctions, and paracrine long-range signaling. However, how these cell interactions shape collective sensing and cell behavior in islets that leads to insulin release is unknown. When stimulated by specific ligands, primarily glucose, beta cells collectively respond with expression of a series of transient Ca$^{2+}$ changes on several temporal scales. Here we analyze a set of Ca$^{2+}$ spike trains recorded in acute rodent pancreatic tissue slice under physiological conditions. We found strongly correlated states of co-spiking cells coexisting with mostly weak pairwise correlations widespread across the islet. Furthermore, the collective Ca$^{2+}$ spiking activity in islet shows on-off intermittency with scaling of spiking amplitudes, and stimulus dependent autoassociative memory features. We use a simple spin glass-like model for the functional network of a beta cell collective to describe these findings and argue that Ca$^{2+}$ spike trains produced by collective sensing of beta cells constitute part of the islet metabolic code that regulates insulin release and limits the islet size.

categories: physics.bio-ph cond-mat.stat-mech q-bio.CB q-bio.QM


title: Bounds for the $\ell$-torsion in class groups

abstract: We prove for each integer $\ell\geq 1$ an unconditional upper bound for the size of the $\ell$-torsion subgroup $Cl_K[\ell]$ of the class group of $K$, which holds for all but a zero density set of number fields $K$ of degree $d\in\{4,5\}$ (with the additional restriction in the case $d = 4$ that the field be non-$D_4$). For sufficiently large $\ell$ this improves recent results of Ellenberg, Matchett Wood, and Pierce, and is also stronger than the best currently known pointwise bounds under GRH. Conditional on GRH and on a weak conjecture on the distribution of number fields our bounds also hold for arbitrary degrees $d$.


title: Inference of topology and the nature of synapses, and the flow of information in neuronal networks

abstract: The characterisation of neuronal connectivity is one of the most important matters in neuroscience. In this work, we show that a recently proposed informational quantity, the causal mutual information, employed with an appropriate methodology, can be used not only to correctly infer the direction of the underlying physical synapses, but also to identify their excitatory or inhibitory nature, considering easy to handle and measure bivariate time-series. The success of our approach relies on a surprising property found in neuronal networks by which non-adjacent neurons do "understand" each other (positive mutual information), however this exchange of information is not capable of causing effect (zero transfer entropy). Remarkably, inhibitory connections, responsible for enhancing synchronisation, transfer more information than excitatory connections, known to enhance entropy in the network. We also demonstrate that our methodology can be used to correctly infer directionality of synapses even in the presence of dynamic and observational Gaussian noise, and is also successful in providing the effective directionality of inter modular connectivity, when only mean fields can be measured.


title: Bridge functional for the molecular density functional theory with consistent pressure and surface tension and its importance for solvation in water

abstract: We address the problem of predicting the solvation free energy and equilibrium solvent density profile in fews minutes from the molecular density functional theory beyond the usual hypernetted-chain approximation. We introduce a bridge functional of a coarse-grained, weighted solvent density. In few minutes at most, for solutes of sizes ranging from small compounds to large proteins, we produce (i) an estimation of the free energy of solvation within 1 kcal/mol of the experimental data for the hydrophobic solutes presented here, and (ii) the solvent distribution around the solute. Contrary to previous propositions, this bridge functional is thermodynamically consistent in that it produces the correct liquid-vapor coexistence and the experimental surface tension. We show this consistency to be of crucial importance for water at room temperature and pressure. This bridge functional is designed to be simple, local, and thus numerically efficient. Finally, we illustrate this new level of molecular theory of solutions with the study of the hydration shell of a protein.

categories: physics.chem-ph cond-mat.mtrl-sci cond-mat.soft cond-mat.stat-mech physics.comp-ph


title: Performance Evaluation of Container-based Virtualization for High Performance Computing Environments

abstract: Virtualization technologies have evolved along with the development of computational environments since virtualization offered needed features at that time such as isolation, accountability, resource allocation, resource fair sharing and so on. Novel processor technologies bring to commodity computers the possibility to emulate diverse environments where a wide range of computational scenarios can be run. Along with processors evolution, system developers have created different virtualization mechanisms where each new development enhanced the performance of previous virtualized environments. Recently, operating system-based virtualization technologies captured the attention of communities abroad (from industry to academy and research) because their important improvements on performance area. In this paper, the features of three container-based operating systems virtualization tools (LXC, Docker and Singularity) are presented. LXC, Docker, Singularity and bare metal are put under test through a customized single node HPL-Benchmark and a MPI-based application for the multi node testbed. Also the disk I/O performance, Memory (RAM) performance, Network bandwidth and GPU performance are tested for the COS technologies vs bare metal. Preliminary results and conclusions around them are presented and discussed.


title: Executive stock option exercise with full and partial information on a drift change point

abstract: We analyse the optimal exercise of an executive stock option (ESO) written on a stock whose drift parameter falls to a lower value at a change point, an exponentially distributed random time independent of the Brownian motion driving the stock. Two agents, who do not trade the stock, have differing information on the change point, and seek to optimally exercise the option by maximising its discounted payoff under the physical measure. The first agent has full information, and observes the change point. The second agent has partial information and filters the change point from price observations. This scenario is designed to mimic the positions of two employees of varying seniority, a fully informed executive and a partially informed less senior employee, each of whom receives an ESO. The partial information scenario yields a model under the observation filtration $\widehat{\mathbb{F}}$ in which the stock drift becomes a diffusion driven by the innovations process, an $\widehat{\mathbb{F}}$-Brownian motion also driving the stock under $\widehat{\mathbb{F}}$, and the partial information optimal stopping value function has two spatial dimensions. We rigorously characterise the free boundary PDEs for both agents, establish shape and regularity properties of the associated optimal exercise boundaries, and prove the smooth pasting property in both information scenarios, exploiting some stochastic flow ideas to do so in the partial information case. We develop finite difference algorithms to numerically solve both agents' exercise and valuation problems and illustrate that the additional information of the fully informed agent can result in exercise patterns which exploit the information on the change point, lending credence to empirical studies which suggest that privileged information of bad news is a factor leading to early exercise of ESOs prior to poor stock price performance.


title: Resilient Learning-Based Control for Synchronization of Passive Multi-Agent Systems under Attack

abstract: In this paper, we show synchronization for a group of output passive agents that communicate with each other according to an underlying communication graph to achieve a common goal. We propose a distributed event-triggered control framework that will guarantee synchronization and considerably decrease the required communication load on the band-limited network. We define a general Byzantine attack on the event-triggered multi-agent network system and characterize its negative effects on synchronization. The Byzantine agents are capable of intelligently falsifying their data and manipulating the underlying communication graph by altering their respective control feedback weights. We introduce a decentralized detection framework and analyze its steady-state and transient performances. We propose a way of identifying individual Byzantine neighbors and a learning-based method of estimating the attack parameters. Lastly, we propose learning-based control approaches to mitigate the negative effects of the adversarial attack.


abstract: The aim of this note is to study the measure-valued Ricci tensor on smooth metric measure space with boundary, which is a generalization of Bakry-Emery's modified Ricci tensor on weighted Riemannian manifold. As an application, we offer a new approach to study curvature-dimension condition of smooth metric measure space with boundary.


title: Riemann surfaces of complex classical trajectories and tunnelling splitting in one-dimensional systems

abstract: The topology of complex classical paths is investigated to discuss quantum tunnelling splittings in one-dimensional systems. Here the Hamiltonian is assumed to be given as polynomial functions, so the fundamental group for the Riemann surface provides complete information on the topology of complex paths, which allows us to enumerate all the possible candidates contributing to the semiclassical sum formula for tunnelling splittings. This naturally leads to action relations among classically disjoined regions, revealing entirely non-local nature in the quantization condition. The importance of the proper treatment of Stokes phenomena is also discussed in Hamiltonians in the normal form.


title: Inflows, Outflows, and a Giant Donor in the Remarkable Recurrent Nova M31N 2008-12a? - Hubble Space Telescope Photometry of the 2015 Eruption

abstract: The recurrent nova M31N 2008-12a experiences annual eruptions, contains a near-Chandrasekhar mass white dwarf, and has the largest mass accretion rate in any nova system. In this paper, we present Hubble Space Telescope (HST) WFC3/UVIS photometry of the late decline of the 2015 eruption. We couple these new data with archival HST observations of the quiescent system and Keck spectroscopy of the 2014 eruption. The late-time photometry reveals a rapid decline to a minimum luminosity state, before a possible recovery / re-brightening in the run-up to the next eruption. Comparison with accretion disk models supports the survival of the accretion disk during the eruptions, and uncovers a quiescent disk mass accretion rate of the order of $10^{-6}\,M_\odot\,\mathrm{yr}^{-1}$, which may rise beyond $10^{-5}\,M_\odot\,\mathrm{yr}^{-1}$ during the super-soft source phase - both of which could be problematic for a number of well-established nova eruption models. Such large accretion rates, close to the Eddington limit, might be expected to be accompanied by additional mass loss from the disk through a wind and even collimated outflows. The archival HST observations, combined with the disk modeling, provide the first constraints on the mass donor; $L_\mathrm{donor}=103^{+12}_{-11}\,L_\odot$, $R_\mathrm{donor}=14.14^{+0.46}_{-0.47}\,R_\odot$, and $T_\mathrm{eff, donor}=4890\pm110$ K, which may be consistent with an irradiated M31 red-clump star. Such a donor would require a system orbital period $\gtrsim5$ days. Our updated analysis predicts that the M31N 2008-12a WD could reach the Chandrasekhar mass in < 20 kyr.


title: Energy Constrained Depth First Search

abstract: Depth first search is a natural algorithmic technique for constructing a closed route that visits all vertices of a graph. The length of such route equals, in an edge-weighted tree, twice the total weight of all edges of the tree and this is asymptotically optimal over all exploration strategies. This paper considers a variant of such search strategies where the length of each route is bounded by a positive integer $B$ (e.g. due to limited energy resources of the searcher). The objective is to cover all the edges of a tree $T$ using the minimum number of routes, each starting and ending at the root and each being of length at most $B$. To this end, we analyze the following natural greedy tree traversal process that is based on decomposing a depth first search traversal into a sequence of limited length routes. Given any arbitrary depth first search traversal $R$ of the tree $T$, we cover $R$ with routes $R_1,\ldots,R_l$, each of length at most $B$ such that: $R_i$ starts at the root, reaches directly the farthest point of $R$ visited by $R_{i-1}$, then $R_i$ continues along the path $R$ as far as possible, and finally $R_i$ returns to the root. We call the above algorithm \emph{piecemeal-DFS} and we prove that it achieves the asymptotically minimal number of routes $l$, regardless of the choice of $R$. Our analysis also shows that the total length of the traversal (and thus the traversal time) of piecemeal-DFS is asymptotically minimum over all energy-constrained exploration strategies. The fact that $R$ can be chosen arbitrarily means that the exploration strategy can be constructed in an online fashion when the input tree $T$ is not known in advance. Surprisingly, our results show that depth first search is efficient for energy constrained exploration of trees, even though it is known that the same does not hold for energy constrained exploration of arbitrary graphs.


title: Maat: A Platform Service for Measurement and Attestation

abstract: Software integrity measurement and attestation (M&A) are critical technologies for evaluating the trustworthiness of software platforms. To best support these technologies, next generation systems must provide a centralized service for securely selecting, collecting, and evaluating integrity measurements. Centralization of M&A avoids duplication, minimizes security risks to the system, and ensures correct ad- ministration of integrity policies and systems. This paper details the desirable features and properties of such a system, and introduces Maat, a prototype implementation of an M&A service that meets these properties. Maat is a platform service that provides a centralized policy-driven framework for determining which measurement tools and protocols to use to meet the needs of a given integrity evaluation. Maat simplifies the task of integrating integrity measurements into a range of larger trust decisions such as authentication, network access control, or delegated computations.


title: Atmospheric thermal tides and planetary spin I. The complex interplay between stratification and rotation

abstract: Thermal atmospheric tides can torque telluric planets away from spin-orbit synchronous rotation, as observed in the case of Venus. They thus participate to determine the possible climates and general circulations of the atmospheres of these planets. In this work, we write the equations governing the dynamics of thermal tides in a local vertically-stratified section of a rotating planetary atmosphere by taking into account the effects of the complete Coriolis acceleration on tidal waves. This allows us to derive analytically the tidal torque and the tidally dissipated energy, which we use to discuss the possible regimes of tidal dissipation and examine the key role played by stratification. In agreement with early studies, we find that the frequency dependence of the thermal atmospheric tidal torque in the vicinity of synchronization can be approximated by a Maxwell model. This behaviour corresponds to weakly stably stratified or convective fluid layers, as observed in ADLM2016a. A strong stable stratification allows gravity waves to propagate, which makes the tidal torque become negligible. The transition is continuous between these two regimes. The traditional approximation appears to be valid in thin atmospheres and in regimes where the rotation frequency is dominated by the forcing or the buoyancy frequencies. Depending on the stability of their atmospheres with respect to convection, observed exoplanets can be tidally driven toward synchronous or asynchronous final rotation rates. The domain of applicability of the traditional approximation is rigorously constrained by calculations.


title: Signal Yields of keV Electronic Recoils and Their Discrimination from Nuclear Recoils in Liquid Xenon

abstract: We report on the response of liquid xenon to low energy electronic recoils below 15 keV from beta decays of tritium at drift fields of 92 V/cm, 154 V/cm and 366 V/cm using the XENON100 detector. A data-to-simulation fitting method based on Markov Chain Monte Carlo is used to extract the photon yields and recombination fluctuations from the experimental data. The photon yields measured at the two lower fields are in agreement with those from literature; additional measurements at a higher field of 366 V/cm are presented. The electronic and nuclear recoil discrimination as well as its dependence on the drift field and photon detection efficiency are investigated at these low energies. The results provide new measurements in the energy region of interest for dark matter searches using liquid xenon.


title: A Practical Computational Method for the Anisotropic Redshift-Space 3-Point Correlation Function

abstract: We present an algorithm enabling computation of the anisotropic redshift-space galaxy 3-point correlation function (3PCF) scaling as $N^2$, with $N$ the number of galaxies. Our previous work showed how to compute the isotropic 3PCF with this scaling by expanding the radially-binned density field around each galaxy in the survey into spherical harmonics and combining these coefficients to form multipole moments. The $N^2$ scaling occurred because this approach never explicitly required the relative angle between a galaxy pair about the primary galaxy. Here we generalize this work, demonstrating that in the presence of azimuthally-symmetric anisotropy produced by redshift-space distortions (RSD) the 3PCF can be described by two triangle side lengths, two independent total angular momenta, and a spin. This basis for the anisotropic 3PCF allows its computation with negligible additional work over the isotropic 3PCF. We also present the covariance matrix of the anisotropic 3PCF measured in this basis. Our algorithm tracks the full 5-D redshift-space 3PCF, uses an accurate line of sight to each triplet, is exact in angle, and easily handles edge correction. It will enable use of the anisotropic large-scale 3PCF as a probe of RSD in current and upcoming large-scale redshift surveys.


title: Stable rationality of Brauer-Severi surface bundles

abstract: For sufficiently ample linear systems on rational surfaces we show that a very general associated Brauer-Severi surface bundle is not stably rational.


title: A Simple and Fast Algorithm for L1-norm Kernel PCA

abstract: We present an algorithm for L1-norm kernel PCA and provide a convergence analysis for it. While an optimal solution of L2-norm kernel PCA can be obtained through matrix decomposition, finding that of L1-norm kernel PCA is not trivial due to its non-convexity and non-smoothness. We provide a novel reformulation through which an equivalent, geometrically interpretable problem is obtained. Based on the geometric interpretation of the reformulated problem, we present a fixed-point type algorithm that iteratively computes a binary weight for each observation. As the algorithm requires only inner products of data vectors, it is computationally efficient and the kernel trick is applicable. In the convergence analysis, we show that the algorithm converges to a local optimal solution in a finite number of steps. Moreover, we provide a rate of convergence analysis, which has been never done for any L1-norm PCA algorithm, proving that the sequence of objective values converges at a linear rate. In numerical experiments, we show that the algorithm is robust in the presence of entry-wise perturbations and computationally scalable, especially in a large-scale setting. Lastly, we introduce an application to outlier detection where the model based on the proposed algorithm outperforms the benchmark algorithms.


title: Monoparametric family of metrics derived from classical Jensen-Shannon divergence

abstract: Jensen-Shannon divergence is a well known multi-purpose measure of dissimilarity between probability distributions. It has been proven that the square root of this quantity is a true metric in the sense that, in addition to the basic properties of a distance, it also satisfies the triangle inequality. In this work we extend this last result to prove that in fact it is possible to derive a monoparametric family of metrics from the classical Jensen-Shannon divergence. Motivated by our results, an application into the field of symbolic sequences segmentation is explored. Additionally, we analyze the possibility to extend this result into the quantum realm.


title: Finite-Time Distributed Linear Equation Solver for Minimum $l_1$ Norm Solutions

abstract: This paper proposes distributed algorithms for multi-agent networks to achieve a solution in finite time to a linear equation $Ax=b$ where $A$ has full row rank, and with the minimum $l_1$-norm in the underdetermined case (where $A$ has more columns than rows). The underlying network is assumed to be undirected and fixed, and an analytical proof is provided for the proposed algorithm to drive all agents' individual states to converge to a common value, viz a solution of $Ax=b$, which is the minimum $l_1$-norm solution in the underdetermined case. Numerical simulations are also provided as validation of the proposed algorithms.


title: Linearized Stability of Extreme Black Holes

abstract: Extreme black holes have been argued to be unstable, in the sense that under linearized gravitational perturbations of the extreme Kerr spacetime the Weyl scalar $\psi_4$ blows up along their event horizons at very late advanced times. We show numerically, by solving the Teukolsky equation in 2+1D, that all algebraically-independent curvature scalar polynomials approach limits that exist when advanced time along the event horizon approaches infinity. Therefore, the horizons of extreme black holes are stable against linearized gravitational perturbations. We argue that the divergence of $\psi_4$ is a consequence of the choice of a fixed tetrad, and that in a suitable dynamical tetrad all Weyl scalars, including $\psi_4$, approach their background extreme Kerr values. We make similar conclusions also for the case of scalar field perturbations of extreme Kerr.


title: A Hybrid Algorithm for Period Analysis from Multi-band Data with Sparse and Irregular Sampling for Arbitrary Light Curve Shapes

abstract: Ongoing and future surveys with repeat imaging in multiple bands are producing (or will produce) time-spaced measurements of brightness, resulting in the identification of large numbers of variable sources in the sky. A large fraction of these are periodic variables: compilations of these are of scientific interest for a variety of purposes. Unavoidably, the data-sets from many such surveys not only have sparse sampling, but also have embedded frequencies in the observing cadence that beat against the natural periodicities of any object under investigation. Such limitations can make period determination ambiguous and uncertain. For multi-band data sets with asynchronous measurements in multiple pass-bands, we want to maximally utilize the information on periodicity in a manner that is agnostic of differences in the light curve shapes across the different channels. Given large volumes of data, computational efficiency is also at a premium. This paper develops and presents a computationally economic method for determining periodicity which combines the results from two different classes of period determination algorithms. The underlying principles are illustrated through examples. The effectiveness of this approach for combining asynchronously sampled measurements in multiple observables that share an underlying fundamental frequency is also demonstrated.


title: A Distributed Algorithm for Least Square Solutions of Linear Equations

abstract: A distributed discrete-time algorithm is proposed for multi-agent networks to achieve a common least squares solution of a group of linear equations, in which each agent only knows some of the equations and is only able to receive information from its nearby neighbors. For fixed, connected, and undirected networks, the proposed discrete-time algorithm results in each agents solution estimate to converging exponentially fast to the same least squares solution. Moreover, the convergence does not require careful choices of time-varying small step sizes.


title: Nearby and vanishing cycles for perverse sheaves and D-modules

abstract: We survey nearby and vanishing cycles for both perverse sheaves and D-modules under analytic setting. Following ideas of A. Beilinson, M. Kashiwara and M. Saito, we explain in detail the proof of the comparison theorem between them in the sense of Riemann-Hilbert correspondence.


title: A Web of Hate: Tackling Hateful Speech in Online Social Spaces

abstract: Online social platforms are beset with hateful speech - content that expresses hatred for a person or group of people. Such content can frighten, intimidate, or silence platform users, and some of it can inspire other users to commit violence. Despite widespread recognition of the problems posed by such content, reliable solutions even for detecting hateful speech are lacking. In the present work, we establish why keyword-based methods are insufficient for detection. We then propose an approach to detecting hateful speech that uses content produced by self-identifying hateful communities as training data. Our approach bypasses the expensive annotation process often required to train keyword systems and performs well across several established platforms, making substantial improvements over current state-of-the-art approaches.


title: $\rm{FP}_{n}$-injective and $\rm{FP}_{n}$-flat covers and preenvelopes, and Gorenstein AC-flat covers

abstract: We prove that, for any $n \geq 2$, the classes of $\rm{FP}_{n}$-injective modules and of $\rm{FP}_n$-flat modules are both covering and preenveloping over any ring $R$. This includes the case of $\rm{FP}_{\infty}$-injective and $\rm{FP}_{\infty}$-flat modules (i.e. absolutely clean and, respectively, level modules). Then we consider a generalization of the class of (strongly) Gorenstein flat modules - the (strongly) Gorenstein AC-flat modules (cycles of exact complexes of flat modules that remain exact when tensored with any absolutely clean module). We prove that some of the properties of Gorenstein flat modules extend to the class of Gorenstein AC-flat modules; for example we show that this class is precovering over any ring $R$. We also show that (as in the case of Gorenstein flat modules) every Gorenstein AC-flat module is a direct summand of a strongly Gorenstein AC-flat module. When $R$ is such that the class of Gorenstein AC-flat modules is closed under extensions, the converse is also true. We also prove that if the class of Gorenstein AC-flat modules is closed under extensions, then this class of modules is covering.


title: First Characterization of the Neutral ISM in Two Local Volume Dwarf Galaxies

abstract: We present the first HI spectral line images of the nearby, star-forming dwarf galaxies UGC11411 and UGC 8245, acquired as part of the "Observing for University Classes" program with the Karl G. Jansky Very Large Array (VLA). These low-resolution images localize the HI gas and reveal the bulk kinematics of each system. Comparing with HST broadband and ground-based H{\alpha} imaging, we find that the ongoing star formation in each galaxy is associated with the highest HI mass surface density regions. UGC 8245 has a much lower current star formation rate than UGC 11411, which harbors very high surface brightness H{\alpha} emission in the inner disk and diffuse, lower surface brightness nebular gas that extends well beyond the stellar disk as traced by HST. We measure the dynamical masses of each galaxy and find that the halo of UGC 11411 is more than an order of magnitude more massive than the halo of UGC 8245, even though the HI and stellar masses of the sources are similar. We show that UGC8245 shares similar physical properties with other well-studied low-mass galaxies, while UGC11411 is more highly dark matter dominated. Both systems have negative peculiar velocities that are associated with a coherent flow of nearby galaxies at high supergalactic latitude.


title: On global Okounkov bodies of spherical varieties

abstract: We define and study the global Okounkov moment cone of a projective spherical variety X, generalizing both the global Okounkov body and the moment body of X defined by Kaveh and Khovanskii. Under mild assumptions on X we show that the global Okounkov moment cone of X is rational polyhedral. As a consequence, also the global Okounkov body of X, with respect to a particular valuation, is rational polyhedral.


title: Deep TAMER: Interactive Agent Shaping in High-Dimensional State Spaces

abstract: While recent advances in deep reinforcement learning have allowed autonomous learning agents to succeed at a variety of complex tasks, existing algorithms generally require a lot of training data. One way to increase the speed at which agents are able to learn to perform tasks is by leveraging the input of human trainers. Although such input can take many forms, real-time, scalar-valued feedback is especially useful in situations where it proves difficult or impossible for humans to provide expert demonstrations. Previous approaches have shown the usefulness of human input provided in this fashion (e.g., the TAMER framework), but they have thus far not considered high-dimensional state spaces or employed the use of deep learning. In this paper, we do both: we propose Deep TAMER, an extension of the TAMER framework that leverages the representational power of deep neural networks in order to learn complex tasks in just a short amount of time with a human trainer. We demonstrate Deep TAMER's success by using it and just 15 minutes of human-provided feedback to train an agent that performs better than humans on the Atari game of Bowling - a task that has proven difficult for even state-of-the-art reinforcement learning methods.


title: Universal and accessible entropy estimation using a compression algorithm

abstract: Entropy and free-energy estimation are key in thermodynamic characterization of simulated systems ranging from spin models through polymers, colloids, protein structure, and drug-design. Current techniques suffer from being model specific, requiring abundant computation resources and simulation at conditions far from the studied realization. Here, we present a universal scheme to calculate entropy using lossless compression algorithms and validate it on simulated systems of increasing complexity. Our results show accurate entropy values compared to benchmark calculations while being computationally effective. In molecular-dynamics simulations of protein folding, we exhibit unmatched detection capability of the folded states by measuring previously undetectable entropy fluctuations along the simulation timeline. Such entropy evaluation opens a new window onto the dynamics of complex systems and allows efficient free-energy calculations.

categories: cond-mat.stat-mech cond-mat.soft physics.bio-ph physics.comp-ph


title: Orthosymplectic Jordan superalgebras and the Wedderburn principal theorem (WPT)

abstract: An analogue of the Wedderbur principal theorem (WPT) is considered for finite dimensional Jordan superalgebras A with solvable radical N, such that N^2=0 and A/N is isomorphic to Josp_n|2m(F), where F is an algebraicallly closed field of characteristic zero. Let's we prove that the WPT is valid under some restrictions over the irreducible Josp_n|2m(F)-bimodules contained in N, and it is shown with counter-examples that these restrictions can not be weakened.


title: Emergent failures and cascades in power grids: a statistical physics perspective

abstract: We model power grids transporting electricity generated by intermittent renewable sources as complex networks, where line failures can emerge indirectly by noisy power input at the nodes. By combining concepts from statistical physics and the physics of power flows, and taking weather correlations into account, we rank line failures according to their likelihood and establish the most likely way such failures occur and propagate. Our insights are mathematically rigorous in a small-noise limit and are validated with data from the German transmission grid.

categories: physics.soc-ph cond-mat.dis-nn cs.SY math.PR


title: On Regularization of Second Kind Integrals

abstract: We obtain expressions for second kind integrals on non-hyperelliptic $(n,s)$-curves. Such a curve possesses a Weierstrass point at infinity which is a branch point where all sheets of the curve come together. The infinity serves as the basepoint for Abel's map, and the basepoint in the definition of the second kind integrals. We define second kind differentials as having a pole at the infinity, therefore the second kind integrals need to be regularized. We propose the regularization consistent with the structure of the field of Abelian functions on Jacobian of the curve. In this connection we introduce the notion of regularization constant, a uniquely defined free term in the expansion of the second kind integral over a local parameter in the vicinity of the infinity. This is a vector with components depending on parameters of the curve, the number of components is equal to genus of the curve. Presence of the term guarantees consistency of all relations between Abelian functions constructed with the help of the second kind integrals. We propose two methods of calculating the regularization constant, and obtain these constants for $(3,4)$, $(3,5)$, $(3,7)$, and $(4,5)$-curves. By the example of $(3,4)$-curve, we extend the proposed regularization to the case of second kind integrals with the pole at an arbitrary fixed point. Finally, we propose a scheme of obtaining addition formulas, where the second kind integrals, including the proper regularization constants, are used.


title: Looking for infrared counterparts of Fermi/LAT blazar candidates

abstract: The Fermi/LAT telescope is an efficient blazar-detector in the MeV/GeV range. More than 1100 (900) blazars detected above 100 MeV (10 GeV) are clearly associated to BL Lacertae or Flat Spectrum Radio Quasar objects in the Fermi/LAT 3FGL catalogue. This number could significantly increase if multi-wavelength counterparts could be identified for the 573 3FGL blazars with unknown type, or even for the 1010 3FGL unassociated sources which are thought to be dominated by blazars, at least at high galactic latitude. Unfortunately, the size of the Fermi/LAT error box makes multi-wavelength follow-ups difficult. We propose a method to associate "blazar-like" infrared counterparts, having coordinates with a precision of a few arcseconds, to Fermi/LAT blazars and unassociated sources. To reach this goal, we built machine-learning classifiers based on the statistical differences of magnitude measurements obtained by the WISE satellite, between a sample of well-identified infrared blazars and samples of other types of infrared sources located in regions of the sky where no known blazar is present. We provide a list of potential infrared counterparts for 3FGL blazar candidates, along with the associated number of expected false positives. This study contributes to increase the number of well-identified extragalactic blazars and also provides promising blazar targets for the Cherenkov Telescope Array.


title: Gammapy: high level data analysis for extragalactic science cases with the Cherenkov Telescope Array

abstract: The Cherenkov Telescope Array (CTA) observatory will probe the non-thermal universe above 20 GeV up to several hundreds of TeV with a significant improvement in sensitivity and angular resolution compared to current experiments. Its outstanding capabilities will allow to increase the number of extragalactic cosmic accelerators detected at very high energy (VHE) and therefore to better constrain the population of VHE accelerators and the gamma-ray absorption processes in the intergalactic medium. For the first time in the history of imaging atmospheric Cherenkov telescopes (IACTs), CTA will be an open observatory and high-level data will be made available to the astronomical community. Gammapy is an open-source Python package developed by the Cherenkov telescope community that provides tools to simulate the gamma-ray sky and analyse IACT data. The versatile architecture of, and steady user contributions to Gammapy enable a large variety of high-level data analyses. Examples of Gaammapy applications are presented, particularly in the context of extragalactic science with CTA.


title: Controlling intermediate dynamics in a family of quadratic maps

abstract: The intermediate dynamics of composed one-dimensional maps is used to multiply attractors in phase space and create multiple independent bifurcation diagrams which can split apart. Results are shown for the composition of k-paradigmatic quadratic maps with distinct values of parameters generating k-independent bifurcation diagrams with corresponding k orbital points. For specific conditions, the basic mechanism for creating the shifted diagrams is the prohibition of period doubling bifurcations transformed in saddle-node bifurcations.


title: Diagonal stability of a class of discrete-time positive switched systems with delay

abstract: A class of discrete-time nonlinear positive time-delay switched systems with sector-type nonlinearities is studied. Sufficient conditions for the existence of common and switched diagonal Lyapunov--Krasovskii functionals for this system class are derived; these are expressed as feasibility conditions for systems of linear algebraic inequalities. Corresponding spectral conditions for the existence of common L--K functionals are also described. Furthermore, it is shown that the proposed approaches can be applied to discrete-time models of digital filters and neural networks. Finally, a numerical example is given to illustrate the effectiveness of theoretical results.


title: Light stops at exceptional points

abstract: Almost twenty years ago the light was slowed down to less than $10^{-7}$ of its vacuum speed in a cloud of ultracold atoms of sodium. Upon a sudden turn-off of the coupling laser, a slow light pulse can be imprinted on cold atoms such that it can be read out and converted into photon again. In this process, the light is stopped by absorbing it and storing its shape within the atomic ensemble. Alternatively, the light can be stopped at the band edge in photonic-crystal waveguides, where the group speed vanishes. Here we extend the phenomenon of stopped light to the new field of parity-time (PT) symmetric systems. We show that zero group speed in PT symmetric optical waveguides can be achieved if the system is prepared at an exceptional point, where two optical modes coalesce. This effect can be tuned for optical pulses in a wide range of frequencies and bandwidths, as we demonstrate in a system of coupled waveguides with gain and loss.


title: X-ray luminescence computed tomography using a focused X-ray beam

abstract: Due to the low X-ray photon utilization efficiency and low measurement sensitivity of the electron multiplying charge coupled device (EMCCD) camera setup, the collimator based narrow beam X-ray luminescence computed tomography (XLCT) usually requires a long measurement time. In this paper, we, for the first time, report a focused X-ray beam based XLCT imaging system with measurements by a single optical fiber bundle and a photomultiplier tube (PMT). An X-ray tube with a polycapillary lens was used to generate a focused X-ray beam whose X-ray photon density is 1200 times larger than a collimated X-ray beam. An optical fiber bundle was employed to collect and deliver the emitted photons on the phantom surface to the PMT. The total measurement time was reduced to 12.5 minutes. For numerical simulations of both single and six fiber bundle cases, we were able to reconstruct six targets successfully. For the phantom experiment, two targets with an edge-to-edge distance of 0.4 mm and a center-to-center distance of 0.8 mm were successfully reconstructed by the measurement setup with a single fiber bundle and a PMT.


title: Stability of a Giant Connected Component in a Complex Network

abstract: We analyze the stability of the network's giant connected component under impact of adverse events, which we model through the link percolation. Specifically, we quantify the extent to which the largest connected component of a network consists of the same nodes, regardless of the specific set of deactivated links. Our results are intuitive in the case of single-layered systems: the presence of large degree nodes in a single-layered network ensures both its robustness and stability. In contrast, we find that interdependent networks that are robust to adverse events have unstable connected components. Our results bring novel insights to the design of resilient network topologies and the reinforcement of existing networked systems.


title: Imaging objects through scattering layers and around corners by retrieval of the scattered point spread function

abstract: We demonstrate a high-speed method to image objects through a thin scattering medium and around a corner. The method employs a reference object of known shape to retrieve the speckle-like point spread function of the scatterer. We extract the point spread function of the scatterer from a dynamic scene that includes a static reference object, and use this to image the dynamic objects. Sharp images are reconstructed from the transmission through a diffuser and from reflection off a rough surface. The sharp and clean reconstructed images from single shot data exemplify the robustness of the method.


title: A unified framework for polycrystal plasticity with grain boundary evolution

abstract: Plastic deformation in polycrystals is governed by the interplay between intra-granular slip and grain boundary-mediated plasticity. However, while the role played by bulk dislocations is relatively well-understood, the contribution of grain boundaries (GBs) has only recently begun to be studied. GB plasticity is known to play a key role along with bulk plasticity under a wide range of conditions, such as dynamic recovery, superplasticity, severe plastic deformation , etc., and developing models capable of simultaneously capturing GB and bulk plasticity has become a topic of high relevance. In this paper we develop a thermodynamically-consistent polycrystal plasticity model capable of simulating a variety of grain boundary-mediated plastic processes in conjunction with bulk dislocation slip. The model starts from the description of a single crystal and creates lattice strain-free polycrystalline configurations by using a specially-designed multiplicative decomposition developed by the authors. This leads to the introduction of a particular class of geometrically necessary dislocations (GND) that define fundamental GB features such as misorientation and inclination. The evolution of the system is based on an energy functional that uses a non-standard function of the GND tensor to account for the grain boundary energy, as well as for the standard elastic energy. Our implementation builds on smooth descriptions of GBs inspired on diffuse-interface models of grain evolution for numerical convenience. We demonstrate the generality and potential of the methodology by simulating a wide variety of phenomena such as shear-induced GB sliding, coupled GB motion, curvature-induced grain rotation and shrinkage, and polygonization via dislocation sub-grain formation.


title: Recognition of feature curves on 3D shapes using an algebraic approach to Hough transforms

abstract: Feature curves are largely adopted to highlight shape features, such as sharp lines, or to divide surfaces into meaningful segments, like convex or concave regions. Extracting these curves is not sufficient to convey prominent and meaningful information about a shape. We have first to separate the curves belonging to features from those caused by noise and then to select the lines, which describe non-trivial portions of a surface. The automatic detection of such features is crucial for the identification and/or annotation of relevant parts of a given shape. To do this, the Hough transform (HT) is a feature extraction technique widely used in image analysis, computer vision and digital image processing, while, for 3D shapes, the extraction of salient feature curves is still an open problem. Thanks to algebraic geometry concepts, the HT technique has been recently extended to include a vast class of algebraic curves, thus proving to be a competitive tool for yielding an explicit representation of the diverse feature lines equations. In the paper, for the first time we apply this novel extension of the HT technique to the realm of 3D shapes in order to identify and localize semantic features like patterns, decorations or anatomical details on 3D objects (both complete and fragments), even in the case of features partially damaged or incomplete. The method recognizes various features, possibly compound, and it selects the most suitable feature profiles among families of algebraic curves.


title: Distance ideals of graphs

abstract: We introduce the concept of distance ideals of graphs, which can be regarded as a generalization of the Smith normal form and the spectra of the distance matrix of a graph. We obtain a classification of the graphs with at most one trivial distance ideal.


title: Complex action suggests future-included theory

abstract: In quantum theory its action is usually taken to be real, but we can consider another theory whose action is complex. In addition, in the Feynman path integral, the time integration is usually performed over the period between the initial time $T_A$ and some specific time, say, the present time $t$. Besides such a future-not-included theory, we can consider the future-included theory, in which not only the past state $| A(T_A) \rangle$ at the initial time $T_A$ but also the future state $| B(T_B) \rangle$ at the final time $T_B$ is given at first, and the time integration is performed over the whole period from the past to the future. Thus quantum theory can be classified into four types, according to whether its action is real or not, and whether the future is included or not. We argue that, if a theory is described with a complex action, then such a theory is suggested to be the future-included theory, rather than the future-not-included theory. Otherwise persons living at different times would see different histories of the universe.


title: Possibilistic Fuzzy Local Information C-Means for Sonar Image Segmentation

abstract: Side-look synthetic aperture sonar (SAS) can produce very high quality images of the sea-floor. When viewing this imagery, a human observer can often easily identify various sea-floor textures such as sand ripple, hard-packed sand, sea grass and rock. In this paper, we present the Possibilistic Fuzzy Local Information C-Means (PFLICM) approach to segment SAS imagery into sea-floor regions that exhibit these various natural textures. The proposed PFLICM method incorporates fuzzy and possibilistic clustering methods and leverages (local) spatial information to perform soft segmentation. Results are shown on several SAS scenes and compared to alternative segmentation approaches.


title: Viscosity of a Multi-channel One-Dimensional Fermi Gas

abstract: Many one-dimensional systems of experimental interest possess multiple bands arising from shallow confining potentials. In this work, we study a gas of weakly interacting fermions and show that the bulk viscosity is dramatically altered by the occupation of more than one band. The reasons for this are two-fold: a multichannel system is more easily displaced from equilibrium and the associated relaxation processes lead to more rapid equilibration than in the single channel case. We estimate the bulk viscosity in terms of the underlying microscopic interactions. The experimental relevance of this physics is discussed in the context of quantum wires and trapped cold atomic gases.


title: Difference in proton radii of mirror nuclei as a possible surrogate for the neutron skin

abstract: It has been recently suggested that differences in the charge radii of mirror nuclei are proportional to the neutron-skin thickness of neutron-rich nuclei and to the slope of the symmetry energy $L$ [B.A. Brown Phys. Rev. Lett. 119, 122502 (2017)]. The determination of the neutron skin has important implications for nuclear physics and astrophysics. Although the use of electroweak probes provides a largely model-independent determination of the neutron skin, the experimental challenges are enormous. Thus, the possibility that differences in the charge radii of mirror nuclei may be used as a surrogate for the neutron skin is a welcome alternative. To test the validity of this assumption we perform calculations based on a set of relativistic energy density functionals that span a wide region of values of $L$. Our results confirm that the difference in charge radii of various neutron-deficient nickel isotopes and their corresponding mirror nuclei is indeed strongly correlated to both the neutron-skin thickness and $L$. Moreover, given that various neutron-star properties are also sensitive to $L$, a data-to-data relation emerges between the difference in charge radii of mirror nuclei and the radius of low-mass neutron stars.


title: On some geometrical methods leading to martingales useful in measure theory

abstract: In this paper I show a geometrical method to obtain construction of the Nikodym type set.


title: On the Time Dependence of Holographic Complexity

abstract: We evaluate the full time dependence of holographic complexity in various eternal black hole backgrounds using both the complexity=action (CA) and the complexity=volume (CV) conjectures. We conclude using the CV conjecture that the rate of change of complexity is a monotonically increasing function of time, which saturates from below to a positive constant in the late time limit. Using the CA conjecture for uncharged black holes, the holographic complexity remains constant for an initial period, then briefly decreases but quickly begins to increase. As observed previously, at late times, the rate of growth of the complexity approaches a constant, which may be associated with Lloyd's bound on the rate of computation. However, we find that this late time limit is approached from above, thus violating the bound. Adding a charge to the eternal black holes washes out the early time behaviour, i.e., complexity immediately begins increasing with sufficient charge, but the late time behaviour is essentially the same as in the neutral case. We also evaluate the complexity of formation for charged black holes and find that it is divergent for extremal black holes, implying that the states at finite chemical potential and zero temperature are infinitely more complex than their finite temperature counterparts.


title: Stochastic thermodynamics for a periodically driven single-particle pump

abstract: We present the stochastic thermodynamic analysis of a time-periodic single particle pump, including explicit results for flux, thermodynamic force, entropy production, work, heat and efficiency. These results are valid far from equilibrium. The deviations from the linear (Onsager) regime are discussed.


title: Development of a focused-X-ray luminescence tomography (FXLT) system

abstract: Biophotonics is an active research area in molecular imaging, genetic diagnosis and prognosis, with direct applicability in precision medicine. However, long-standing challenges of biophotonics are well known due to low signal-to-noise ratio and poor image quality, mainly due to strong optical scattering especially in deep tissues. Recently, X-ray luminescence computed tomography (XLCT) has emerged as a hybrid molecular imaging modality and shown great promises in overcoming the optical scattering in deep tissues. However, its high spatial resolution capacity has not been fully implemented yet. In this paper, with a superfine focused X-ray beam we design a focused-X-ray luminescence tomography (FXLT) system for spatial resolution better than 150 micrometers and molecular sensitivity of 2 micromolar. First, we describe our system design. Then, we demonstrate that the molecular sensitivity of FXLT is about 5 micromolar considering the emitted visible photons from background. Furthermore, we analyze the measurement time per scan from measured photons numbers with a fiber bundle-PMT setup, report numerical and experimental results. Finally, we specify the imaging system performance based on numerical simulations and physical measurements.