title: Mechanically Verified Calculational Abstract Interpretation

abstract: Calculational abstract interpretation, long advocated by Cousot, is a technique for deriving correct-by-construction abstract interpreters from the formal semantics of programming languages. This paper addresses the problem of deriving correct-by-verified-construction abstract interpreters with the use of a proof assistant. We identify several technical challenges to overcome with the aim of supporting verified calculational abstract interpretation that is faithful to existing pencil-and-paper proofs, supports calculation with Galois connections generally, and enables the extraction of verified static analyzers from these proofs. To meet these challenges, we develop a theory of Galois connections in monadic style that include a specification effect. Effectful calculations may reason classically, while pure calculations have extractable computational content. Moving between the worlds of specification and implementation is enabled by our metatheory. To validate our approach, we give the first mechanically verified proof of correctness for Cousot's "Calculational design of a generic abstract interpreter." Our proof "by calculus" closely follows the original paper-and-pencil proof and supports the extraction of a verified static analyzer.


title: The Dirac Operator on Regular Metric Trees

abstract: A metric tree is a tree whose edges are viewed as line segments of positive length. The Dirac operator on such tree is the operator which operates on each edge, complemented by the matching conditions at the vertices which were given by Bolte and Harrison \cite{BolteHarrison2003}. The spectrum of Dirac operator can be quite different, reflecting geometry of the tree. We discuss a special case of trees, namely the so-called regular trees. They possess a rich group of symmetries. This allows one to construct an orthogonal decomposition of the space $L^2(\Gamma)$ which reduces the Dirac. Based upon this decomposition, a detailed spectral analysis of Dirac operator on the regular metric trees is possible.


title: Thermodynamics of the classical spin-ice model with nearest neighbour interactions using the Wang-Landau algorithm

abstract: In this article we study the classical nearest-neighbour spin-ice model (nnSI) by means of Monte Carlo simulations, using the Wang-Landau algorithm. The nnSI describes several of the salient features of the spin-ice materials. Despite its simplicity it exhibits a remarkably rich behaviour. The model has been studied using a variety of techniques, thus it serves as an ideal benchmark to test the capabilities of the Wang Landau algorithm in magnetically frustrated systems. We study in detail the residual entropy of the nnSI and, by introducing an applied magnetic field in two different crystallographic directions ([111] and [100],) we explore the physics of the kagome-ice phase, the transition to full polarisation, and the three dimensional Kasteleyn transition. In the latter case, we discuss how additional constraints can be added to the Hamiltonian, by taking into account a selective choice of states in the partition function and, then, show how this choice leads to the realization of the ideal Kasteleyn transition in the system.


title: Predicting Scheduling Failures in the Cloud

abstract: Cloud Computing has emerged as a key technology to deliver and manage computing, platform, and software services over the Internet. Task scheduling algorithms play an important role in the efficiency of cloud computing services as they aim to reduce the turnaround time of tasks and improve resource utilization. Several task scheduling algorithms have been proposed in the literature for cloud computing systems, the majority relying on the computational complexity of tasks and the distribution of resources. However, several tasks scheduled following these algorithms still fail because of unforeseen changes in the cloud environments. In this paper, using tasks execution and resource utilization data extracted from the execution traces of real world applications at Google, we explore the possibility of predicting the scheduling outcome of a task using statistical models. If we can successfully predict tasks failures, we may be able to reduce the execution time of jobs by rescheduling failed tasks earlier (i.e., before their actual failing time). Our results show that statistical models can predict task failures with a precision up to 97.4%, and a recall up to 96.2%. We simulate the potential benefits of such predictions using the tool kit GloudSim and found that they can improve the number of finished tasks by up to 40%. We also perform a case study using the Hadoop framework of Amazon Elastic MapReduce (EMR) and the jobs of a gene expression correlations analysis study from breast cancer research. We find that when extending the scheduler of Hadoop with our predictive models, the percentage of failed jobs can be reduced by up to 45%, with an overhead of less than 5 minutes.


title: Systematically characterizing regions of the First Fermi-LAT SNR Catalog

abstract: While supernova remnants (SNRs) are widely thought to be powerful cosmic-ray accelerators, indirect evidence comes from a small number of well-studied cases. Here we systematically determine the gamma-ray emission detected by the Fermi Large Area Telescope (LAT) from all known Galactic SNRs, disentangling them from the sea of cosmic-ray generated photons in the Galactic plane. Using LAT data we have characterized the 1-100 GeV emission in 279 regions containing SNRs, accounting for systematic uncertainties caused by source confusion and instrumental response. We have also developed a method to explore some systematic effects on SNR properties caused by the modeling of the interstellar emission (IEM). The IEM contributes substantially to gamma-ray emission in the regions where SNRs are located. To explore the systematics we consider different model construction methods, different model input parameters, and independently fit the model components to the gamma-ray data. We will describe this analysis method in detail. In the First Fermi LAT SNR Catalog there are 30 sources classified as SNRs, using spatial overlap with the radio position. For all the remaining regions we evaluated upper limits on SNRs' emission. In this work we will present a study of the aggregate characteristics of SNRs, such as comparisons between GeV and radio sizes as well as fluxes and spectral indexes and with TeV.


title: Extensions of the universal theta divisor

abstract: The Jacobian varieties of smooth curves fit together to form a family, the universal Jacobian, over the moduli space of smooth marked curves, and the theta divisors of these curves form a divisor in the universal Jacobian. In this paper we describe how to extend these families over the moduli space of stable marked curves (or rather an open subset thereof) using a stability parameter. We then prove a wall-crossing formula describing how the theta divisor varies with the stability parameter. We use that result to analyze a divisor on the moduli space of smooth marked curves that has recently been studied by Grushevsky-Zakharov, Hain and M\"uller. In particular, we compute the pullback of the theta divisor studied in Alexeev's work on stable abelic varieties and in Caporaso's work on theta divisors of compactified Jacobians.


title: Test for exotic isoscalar resonance dominating $D^0 \to \pi^+ \pi^- \pi^0$ decays

abstract: The decay $D^0 \to \pi^+ \pi^- \pi^0$ appears to be dominated by $\rho \pi$ states in a configuration of zero total isotopic spin. The spin $J$, parity $P$, and charge-conjugation eigenvalue $C$ of this final state are therefore $J^{PC} = 0^{--}$, which cannot be formed of a quark $q$ and antiquark $\bar q$. If a resonance near $M(D^0)$ dominates the final state, it must be a {\it hybrid} composed of a quark-antiquark pair and a constituent gluon, or a {\it tetraquark} $q q \bar q \bar q$. A test for this resonance in electroproduction is proposed.


title: Low-rank Solutions of Linear Matrix Equations via Procrustes Flow

abstract: In this paper we study the problem of recovering a low-rank matrix from linear measurements. Our algorithm, which we call Procrustes Flow, starts from an initial estimate obtained by a thresholding scheme followed by gradient descent on a non-convex objective. We show that as long as the measurements obey a standard restricted isometry property, our algorithm converges to the unknown matrix at a geometric rate. In the case of Gaussian measurements, such convergence occurs for a $n_1 \times n_2$ matrix of rank $r$ when the number of measurements exceeds a constant times $(n_1+n_2)r$.


title: Direct Method on Stochastic Maximum Principle for Optimization with Recursive Utilities

abstract: We obtain the variational equations for backward stochastic differential equations in recursive stochastic optimal control problems, and then get the maximum principle which is novel. The control domain need not be convex, and the generator of the backward stochastic differential equation can contain z.


title: Higgs decay to two photons at one-loop in the SMEFT

abstract: We present the calculation of the $\rm CP$ conserving contributions to $\Gamma(h \rightarrow \gamma \gamma)$, from dimension six operators at one-loop order, in the linear Standard Model Effective Field Theory. We discuss the impact of these corrections on interpreting current and future experimental bounds on this decay.


title: The Segal--Bargmann transform for odd-dimensional hyperbolic spaces

abstract: We develop isometry and inversion formulas for the Segal--Bargmann transform on odd-dimensional hyperbolic spaces that are as parallel as possible to the dual case of odd-dimensional spheres.


title: Extreme violation of local realism in quantum hypergraph states

abstract: Hypergraph states form a family of multiparticle quantum states that generalizes the well-known concept of Greenberger-Horne-Zeilinger states, cluster states, and more broadly graph states. We study the nonlocal properties of quantum hypergraph states. We demonstrate that the correlations in hypergraph states can be used to derive various types of nonlocality proofs, including Hardy-type arguments and Bell inequalities for genuine multiparticle nonlocality. Moreover, we show that hypergraph states allow for an exponentially increasing violation of local realism which is robust against loss of particles. Our results suggest that certain classes of hypergraph states are novel resources for quantum metrology and measurement-based quantum computation.


title: Curvature of the chiral pseudo-critical line in QCD: continuum extrapolated results

abstract: We determine the curvature of the pseudo-critical line of strong interactions by means of numerical simulations at imaginary chemical potentials. We consider $N_f=2+1$ stout improved staggered fermions with physical quark masses and the tree level Symanzik gauge action, and explore four different sets of lattice spacings, corresponding to $N_t = 6,8,10,12$, in order to extrapolate results to the continuum limit. Our final estimate is $\kappa = 0.0135(20)$.


title: BioNetGen 2.2: Advances in Rule-Based Modeling

abstract: BioNetGen is an open-source software package for rule-based modeling of complex biochemical systems. Version 2.2 of the software introduces numerous new features for both model specification and simulation. Here, we report on these additions, discussing how they facilitate the construction, simulation, and analysis of larger and more complex models than previously possible.


title: Coronal energy input and dissipation in a solar active region 3D MHD model

abstract: Context. We have conducted a 3D MHD simulation of the solar corona above an active region (AR) in full scale and high resolution, which shows coronal loops, and plasma flows within them, similar to observations. Aims. We want to find the connection between the photospheric energy input by field-line braiding with the coronal energy conversion by Ohmic dissipation of induced currents. Methods. To this end we compare the coronal energy input and dissipation within our simulation domain above different fields of view, e.g. for a small loops system in the AR core. We also choose an ensemble of field lines to compare, e.g., the magnetic energy input to the heating per particle along these field lines. Results. We find an enhanced Ohmic dissipation of currents in the corona above areas that also have enhanced upwards-directed Poynting flux. These regions coincide with the regions where hot coronal loops within the AR core are observed. The coronal density plays a role in estimating the coronal temperature due to the generated heat input. A minimum flux density of about 200 Gauss is needed in the photosphere to heat a field line to coronal temperatures of about 1 MK. Conclusions. This suggests that the field-line braiding mechanism provides the coronal energy input and that the Ohmic dissipation of induced currents dominates the coronal heating mechanism.


title: Scaling symmetry and scalar hairy Lifshitz black holes

abstract: By utilizing the scaling symmetry of the reduced action for planar black holes, we obtain the corresponding conserved charge. We use the conserved charge to find the generalized Smarr relation of static hairy planar black holes in various dimensions. Our results not only reproduce the relation in the various known cases but also give the new relation in the Lifshitz planar black holes with the scalar hair.


title: Design and fabrication of a data logger for atmospheric pressure, temperature and relative humidity for gas-filled detector development

abstract: A novel instrument has been developed to monitor and record the ambient pa- rameters such as temperature, atmospheric pressure and relative humidity. These parameters are very essential for understanding the characteristics such as gain of gas filled detectors like Gas Electron Multiplier (GEM) and Multi Wire Propor- tional Counter (MWPC). In this article the details of the design, fabrication and operation processes of the device has been presented.


title: On the Possibility of Improving the Orbits of Satellites Based on Observations of Isolated X-ray Pulsars

abstract: At present, there is a great worldwide interest in the development of technologies that allow information about the X-ray emission from pulsating cosmic sources to be used to obtain navigation solutions for deep-space spacecraft. In this paper, we illustrate the technique for determining the spatial position of a spacecraft based on the already existing data from the RXTE X-ray space observatory. We show that the spacecraft position toward the Crab pulsar can be determined using an X-ray detector with an effective area of about 0.6 sq.m in the energy range 3-15 keV with an accuracy up to 730 m in a signal integration time of 1000 s. Extending the energy range to 1 keV (the efficiency of the RXTE/PCA spectrometer decreases dramatically at energies below 3 keV) will allow a spacecraft position accuracy of 400-450 m to be achieved at the same effective area and up to 300-350 m when using detectors with an effective area of ~1 sq.m in the energy range 1-10 keV.


title: JSKETCH: Sketching for Java

abstract: Sketch-based synthesis, epitomized by the SKETCH tool, lets developers synthesize software starting from a partial program, also called a sketch or template. This paper presents JSKETCH, a tool that brings sketch-based synthesis to Java. JSKETCH's input is a partial Java program that may include holes, which are unknown constants, expression generators, which range over sets of expressions, and class generators, which are partial classes. JSKETCH then translates the synthesis problem into a SKETCH problem; this translation is complex because SKETCH is not object-oriented. Finally, JSKETCH synthesizes an executable Java program by interpreting the output of SKETCH.


title: Micro moon versus macro moon: Brightness and size

abstract: The moon, moonlight, phases of the moon and its relatively simple recurring cycle has been of interest since time immemorial to the human beings, navigators, astronomers and astrologers. The fact that its orbit is elliptical as well its plane is inclined with the plane of rotation of the earth gives rise to new moon to full moon and solar and lunar eclipses. During the phase of the full moon, the luminous flux and its apparent size will depend on its distance from the earth. In case it is at farthest point known as lunar apogee causes smallest full moon or micro full moon and if it is closest to us termed as lunar perigee will result in macro full moon, also known as super moon, a term coined by astrologer Richard Nolle in 1979. The theoretical expressions for the lunar luminous fluxes on the earth representing the power of lunar light the earth intercepts in the direction normal to the incidence over an area of one square meter are derived for two extreme positions lunar apogee and lunar perigee. The expressions for the apparent sizes of full moons corresponding to said positions are also mentioned. It is found that full perigee moon is about 29 percent brighter and 14 percent bigger than the full apogee moon consistent with the reported values.


title: Coupling quantum circuits to magnetic molecular qubits

abstract: This thesis explores the coupling of magnetic systems to quantum circuits in the context of quantum computing applications. In particular we study the coupling of superconducting coplanar waveguide resonators to Single Molecule Magnets (SMMs) . The combination of approaches from the fields of Cavity Quantum electrodynamics (QED) and Circuit QED with those from the field of molecular magnetism con provide unique opportunities for quantum computing. We investigate the necessary conditions for coupling single spins and spin ensembles to resonators and what characteristics SMMs should have in order to provide interesting alternatives as quantum bits. We present test measurements of several magnetic samples using both broadband spectroscopy with open waveguides and EPR spectroscopy using coplanar waveguide resonators. We also design, fabricate, and test nanometric constrictions in superconducting resonators with the objective of improving their coupling to single spins. We evaluate the performance of these constricted resonators in comparison to unmodified resonators.


title: Tentative detection of ethylene glycol toward W51/e2 and G34.3+0.2

abstract: How complex organic - and potentially prebiotic - molecules are formed in regions of low- and high-mass star-formation remains a central question in astrochemistry. In particular, with just a few sources studied in detail, it is unclear what role environment plays in complex molecule formation. In this light, a comparison of relative abundances of related species between sources might be useful to explain observed differences. We seek to measure the relative abundance between three important complex organic molecules, ethylene glycol ((CH$_2$OH)$_2$), glycolaldehyde (CH$_2$OHCHO) and methyl formate (HCOOCH$_3$), toward high-mass protostars and thereby provide additional constraints on their formation pathways. We use IRAM 30-m single dish observations of the three species toward two high-mass star-forming regions - W51/e2 and G34.3+0.2 - and report a tentative detection of (CH2OH)2 toward both sources. Assuming that (CH$_2$OH)$_2$, CH$_2$OHCHO and HCOOCH$_3$ spatially coexist, relative abundance ratios, HCOOCH$_3$/(CH$_2$OH)$_2$, of 31 and 35 are derived for G34.3+0.2 and W51/e2, respectively. CH$_2$OHCHO is not detected, but the data provide lower limits to the HCOOCH$_3$/CH$_2$OHCHO abundance ratios of $\ge$193 for G34.3+0.2 and $\ge$550 for W51/e2. A comparison of these results to measurements from various sources in the literature indicates that the source luminosities may be correlated with the HCOOCH$_3$/(CH$_2$OH)$_2$ and HCOOCH$_3$/CH$_2$OHCHO ratios. This apparent correlation may be a consequence of the relative timescales each source spend at different temperatures-ranges in their evolution. Furthermore, we obtain lower limits to the ratio of (CH$_2$OH)$_2$/CH2OHCHO for G34.3+0.2 ($\ge$6) and W51/e2 ($\ge$16). This result confirms that a high (CH$_2$OH)$_2$/CH$_2$OHCHO abundance ratio is not a specific property of comets, as previously speculated.


title: Quantum digital signature scheme

abstract: Digital signatures are the building blocks of modern communication to prevent masquerading by any party other than recipients, repudiation by signatory and forgery by any individual recipient. Digital signature scheme is said to be standard if the signature (a) is a pattern depending upon the message to be signed, (b) is built upon some information publically known and unique to the signatory, (c) can be stored by all the recipients. While classical methods provide computational security only, quantum mechanics guarantees information-theoretically secure and standard digital signature schemes. However, standard quantum digital signature schemes are based on quantum one-way functions and hence require long term quantum memory for storing quantum signatures, which is not practically feasible yet. We demonstrate here a standard quantum digital signature scheme by replacing quantum one-way functions with multiparty controlled EPR channels. It allows signatory to generate non-locally correlated quantum signatures, instead of multiple copies of a unique quantum state, and assures security against any individual since others have non-locally correlated information.


title: Swift J1112.2-8238: A Candidate Relativistic Tidal Disruption Flare

abstract: We present observations of Swift J1112.2-8238, and identify it as a candidate relativistic tidal disruption flare (rTDF). The outburst was first detected by Swift/BAT in June 2011 as an unknown, long-lived (order of days) $\gamma$-ray transient source. We show that its position is consistent with the nucleus of a faint galaxy for which we establish a likely redshift of $z=0.89$ based on a single emission line that we interpret as the blended [OII]$\lambda3727$ doublet. At this redshift, the peak X/$\gamma$-ray luminosity exceeded $10^{47}$ ergs s$^{-1}$, while a spatially coincident optical transient source had $i^{\prime} \sim 22$ (M$_g \sim -21.4$ at $z=0.89$) during early observations, $\sim 20$ days after the Swift trigger. These properties place Swift J1112.2-8238 in a very similar region of parameter space to the two previously identified members of this class, Swift J1644+57 and Swift J2058+0516. As with those events the high-energy emission shows evidence for variability over the first few days, while late time observations, almost 3 years post-outburst, demonstrate that it has now switched off. Swift J1112.2-8238 brings the total number of such events observed by Swift to three, interestingly all detected by Swift over a $\sim$3 month period ($<3\%$ of its total lifetime as of March 2015). While this suggests the possibility that further examples may be uncovered by detailed searches of the BAT archives, the lack of any prime candidates in the years since 2011 means these events are undoubtedly rare.


title: Enhancement of superconductivity at the onset of charge-density-wave order in a metal

abstract: We analyze superconductivity in the cuprates near the onset of an incommensurate charge density wave (CDW) order with momentum ${\bf Q} = (Q,0)/(0,Q)$, as observed in the experiments. We first consider a semi-phenomenological charge-fermion model in which hot fermions, separated by ${\bf Q}$, attract each other by exchanging soft CDW fluctuations. We find that in a quantum-critical region near CDW transition, $T_c = A {\bar g}_c$, where ${\bar g}_c$ is charge-fermion coupling and $A$ is the prefactor which we explicitly compute. We then consider the particular microscopic scenario in which CDW order parameter emerges as a composite field made out of primary spin-density-wave fields. We show that charge-fermion coupling ${\bar g}_c$ is of order of spin-fermion coupling ${\bar g}_s$. As the consequence, superconducting $T_c$ is substantially enhanced near the onset of CDW order. Finally we analyze the effect of an external magnetic field $H$. We show that, as $H$ increases, optimal $T_c$ decreases and the superconducting dome becomes progressively more confined to the CDW quantum-critical point. These results are consistent with the experiments.


title: Improving reconstruction of the baryon acoustic peak : the effect of local environment

abstract: Precise measurements of the baryon acoustic oscillation (BAO) scale as a standard ruler in the clustering pattern of large-scale structure is a central goal of current and future galaxy surveys. The BAO peak may be sharpened using the technique of density-field reconstruction, in which the bulk displacements of galaxies are estimated using a Zel'dovitch approximation. We use numerical simulations to demonstrate how the accuracy of this approximation depends strongly on local environment, and how this information may be used to construct an improved BAO measurement through environmental re-weighting and using higher-order perturbation theory. We outline further applications of the displacement field for testing cosmological models.


title: The Star-Formation Main Sequence: The Dependence of Specific Star Formation Rate and Its Dispersion on Galaxy Stellar Mass

abstract: The dispersion of the star-formation main sequence (SFMS) reflects the diversity of star formation histories and variation in star formation rates (SFRs) in star-forming galaxies (SFGs) with similar stellar masses ($M^\ast$). We examine the dispersion of local SFMS using a complete sample of Sloan Digital Sky Survey galaxies at 0.01$<z<$0.03 with $\log(M^\ast/M_\odot)>$8.8. The SFRs are estimated from H$\alpha$ in combination with 22$\mu m$ observation from WISE. The catalog of bulge+disk decomposition from Simard et al. (2011) is available for the sample galaxies. We measure the dispersion of specific SFR (SSFR) as a function of $M^*$. We confirm that the dispersion increases with $M^*$ from 0.37$\pm0.01$dex at $\log(M^\ast/M_\odot)<$9.6 to 0.51$\pm0.02$dex at $\log(M^\ast/M_\odot)>$10.2. Despite star formation is mostly associated with disks, the dispersion of disk SSFR still increases with $M^*$. We conclude that the presence of bulges/bars is likely responsible for the large dispersion of SSFR in massive SFGs while low-mass SFGs are mostly disk-dominated and thus with small dispersion. Our results suggest that star formation on galactic scales is dramatically affected by central dense structures through both enhancing and/or quenching processes; while lower-mass SFGs tend to have less bursty star formation histories. However, the dispersion of SSFR becomes significantly smaller and remains constant when only disk-dominated SFGs are counted. This finding implies that the impact of stochastic stellar feedback on star formation is likely to follow the same pattern in all disk galaxies, showing no correlation with halo potential.


title: Patterned Dynamics of Delay-Coupled Swarms with Random Communication Graphs

abstract: Swarm and modular robotics are an emerging area in control of autonomous systems. However, coordinating a large group of interacting autonomous agents requires careful consideration of the logistical issues involved. In particular, inter-agent communication generally involves time delay, and bandwidth restrictions limit the number of neighbors with which each agent in the swarm can communicate. In this paper, we analyze coherent pattern dynamics of groups of delay-coupled agents, where the communication network is an Erdos-Renyi graph. We show that overall motion patterns for a globally-coupled swarm persist under decreasing network connectivity, and derive the bifurcation structure scaling re- lations for the emergence of different swarming behaviors as a function of the average network degree. We show excellent agreement between the theoretical scaling results and numerical simulations.


title: Optical imaging polarimetry of the LkCa 15 protoplanetary disk with SPHERE ZIMPOL

abstract: We present the first optical (590--890 nm) imaging polarimetry observations of the pre-transitional protoplanetary disk around the young solar analog LkCa 15, addressing a number of open questions raised by previous studies. We detect the previously unseen far side of the disk gap, confirm the highly eccentric scattered-light gap shape that was postulated from near-infrared imaging, at odds with the symmetric gap inferred from millimeter interferometry. Furthermore, we resolve the inner disk for the first time and trace it out to 30 AU. This new source of scattered light may contribute to the near-infrared interferometric signal attributed to the protoplanet candidate LkCa 15 b, which lies embedded in the outer regions of the inner disk. Finally, we present a new model for the system architecture of LkCa 15 that ties these new findings together. These observations were taken during science verification of SPHERE ZIMPOL and demonstrate this facility's performance for faint guide stars under adverse observing conditions.


title: A Ground-Based Albedo Upper Limit for HD 189733b from Polarimetry

abstract: We present 50 nights of polarimetric observations of HD 189733 in $B$ band using the POLISH2 aperture-integrated polarimeter at the Lick Observatory Shane 3-m telescope. This instrument, commissioned in 2011, is designed to search for Rayleigh scattering from short-period exoplanets due to the polarized nature of scattered light. Since these planets are spatially unresolvable from their host stars, the relative contribution of the planet-to-total system polarization is expected to vary with an amplitude of order 10 parts per million (ppm) over the course of the orbit. Non-zero and also variable at the 10 ppm level, the inherent polarization of the Lick 3-m telescope limits the accuracy of our measurements and currently inhibits conclusive detection of scattered light from this exoplanet. However, the amplitude of observed variability conservatively sets a $3 \sigma$ upper limit to the planet-induced polarization of the system of 58 ppm in $B$ band, which is consistent with a previous upper limit from the POLISH instrument at the Palomar Observatory 5-m telescope (Wiktorowicz 2009). A physically-motivated Rayleigh scattering model, which includes the depolarizing effects of multiple scattering, is used to conservatively set a $3 \sigma$ upper limit to the geometric albedo of HD 189733b of $A_g < 0.37$. This value is consistent with the value $A_g = 0.226 \pm 0.091$ derived from occultation observations with HST STIS (Evans et al. 2013), but it is inconsistent with the large $A_g = 0.61 \pm 0.12$ albedo reported by (Berdyugina et al. 2011).


title: Green Bank Telescope Constraints on Dark Matter Annihilation in Segue I

abstract: We use a non-detection in $\nu = 1.4\,$GHz Green Bank Telescope observations of the ultra-faint dwarf spheroidal galaxy Segue I, which could be immersed in a non-negligible halo magnetic field of the Milky Way, to place bounds on particle dark matter properties. We model the galaxy using an Einasto dark matter profile, and compute the expected synchrotron flux from dark matter annihilation as a function of the magnetic field strength $B$, diffusion coefficient $D_0$, and particle mass $m_\chi$ for different annihilation channels. The data strongly disfavor annihilations to $e^+e^-$ for $m_\chi \lesssim 50\,$GeV, but are not sensitive to the $b \bar b$ channel. Adopting a fiducial $B \sim 2\,\mu$G inferred from Segue I's proximity to the Milky Way, our models of annihilation to $\tau^+\tau^-$ with $m_\chi = 30\,$GeV require an intermediate value of $D_0$ for consistency with the data. The most compelling limits are obtained for WIMP annihilation to $\mu^+\mu^-$: we exclude $m_\chi \lesssim 30\,$GeV$\,\rightarrow\mu^+\mu^-$ at 95% confidence, unless $D_0$ exceeds the Milky Way value or $B$ is significantly smaller than we have assumed.


title: NIHAO IV: Core creation and destruction in dark matter density profiles across cosmic time

abstract: We use the NIHAO simulations to investigate the effects of baryonic physics on the time evolution of Dark Matter central density profiles. The sample is made of $\approx 70$ independent high resolution hydrodynamical simulations of galaxy formation and covers a wide mass range: 1e10< Mhalo <1e12, i.e., from dwarfs to L* . We confirm previous results on the dependence of the inner dark matter density slope, $\alpha$, on the ratio between stellar-to-halo mass. We show that this relation holds approximately at all redshifts (with an intrinsic scatter of ~0.18 in $\alpha$). This implies that in practically all haloes the shape of their inner density profile changes quite substantially over cosmic time, as they grow in stellar and total mass. Thus, depending on their final stellar-to-halo mass ratio, haloes can either form and keep a substantial density core (size~1 kpc), or form and then destroy the core and re-contract the halo, going back to a cuspy profile, which is even steeper than CDM predictions for massive galaxies (~1e12 Msun). We show that results from the NIHAO suite are in good agreement with recent observational measurements of $\alpha$ in dwarf galaxies. Overall our results suggest that the notion of a universal density profile for dark matter haloes is no longer valid in the presence of galaxy formation.


title: B fields in OB stars (BOB): Detection of a strong magnetic field in the O9.7 V star HD54879

abstract: The number of magnetic stars detected among massive stars is small; nevertheless, the role played by the magnetic field in stellar evolution cannot be disregarded. Links between line profile variability, enhancements/depletions of surface chemical abundances, and magnetic fields have been identified for low-mass B-stars, but for the O-type domain this is almost unexplored. Based on FORS2 and HARPS spectropolarimetric data, we present the first detection of a magnetic field in HD54879, a single slowly rotating O9.7 V star. Using two independent and different techniques we obtained the firm detection of a surface average longitudinal magnetic field with a maximum amplitude of about 600 G, in modulus. A quantitative spectroscopic analysis of the star with the stellar atmosphere code FASTWIND results in an effective temperature and a surface gravity of 33000$\pm1000$ K and 4.0$\pm0.1$ dex. The abundances of carbon, nitrogen, oxygen, silicon, and magnesium are found to be slightly lower than solar, but compatible within the errors. We investigate line-profile variability in HD54879 by complementing our spectra with spectroscopic data from other recent OB-star surveys. The photospheric lines remain constant in shape between 2009 and 2014, although H$\alpha$ shows a variable emission. The H$\alpha$ emission is too strong for a standard O9.7 V and is probably linked to the magnetic field and the presence of circumstellar material. Its normal chemical composition and the absence of photospheric line profile variations make HD54879 the most strongly magnetic, non-variable single O-star detected to date.


title: How to Recover a Qubit That Has Fallen Into a Black Hole

abstract: We demonstrate an algorithm for the retrieval of a qubit, encoded in spin angular momentum, that has been dropped into a no-firewall black hole. Retrieval is achieved analogously to quantum teleportation by collecting Hawking radiation and performing measurements on the black hole. Importantly, these methods only require the ability to perform measurements from outside the event horizon.


title: Dynamical formation & scattering of hierarchical triples: Cross sections, Kozai-Lidov oscillations, and collisions

abstract: Dynamical scattering of binaries and triple systems of stars, planets, and compact objects may produce highly inclined triple systems subject to Kozai-Lidov (KL) oscillations, potentially leading to collisions, mergers, Type Ia supernovae, and other phenomena. We present the results of more than 400 million gravitational scattering experiments of binary-binary, triple-single, and triple-binary scattering. We compute the cross sections for all possible outcomes and explore their dependencies on incoming velocity, mass, semi-major axis, and eccentricity, including analytic fits and discussion of the velocity dependence. For the production of new triple systems by scattering we find that compact triples are preferred, with ratios of outer to inner semi-major axes of ~few--100, flat or quasi-thermal eccentricity distributions, and flat distributions in cosine of the mutual inclination. Dynamically formed triples are thus subject to strong KL oscillations, the "eccentric Kozai mechanism," and non-secular effects. For single and binary flyby encounters with triple systems, we compute the cumulative cross section for changes to the mutual inclination, eccentricity, and semi-major axis ratio. We apply these results to scattering events in the field, open clusters, and globular clusters, and explore the implications for Type Ia supernovae via collisions and mergers, stellar collisions, and the lifetime and dynamical isolation of triple systems undergoing KL oscillations. An Appendix provides an analysis of the velocity dependence of the collision cross section in binary-single scattering.


title: A timing constraint on the (total) mass of the Large Magellanic Cloud

abstract: This paper explores the effect of the LMC on the mass estimates obtained from the timing argument. We show that accounting for the presence of the LMC systematically lowers the Local Group mass ($M_{\rm LG}$) derived from the relative motion of the Milky Way--Andromeda pair. Motivated by this result we apply a Bayesian technique devised by Pe\~narrubia et al. (2014) to simultaneously fit (i) distances and velocities of galaxies within 3~Mpc and (ii) the relative motion between the Milky Way and Andromeda derived from HST observations, with the LMC mass ($M_{\rm LMC}$) as a free parameter. Our analysis returns a Local Group mass $M_{\rm LG}=2.64^{+0.42}_{-0.38}\times 10^{12}M_\odot$ at a 68\% confidence level. The masses of the Milky Way, $M_{\rm MW}=1.04_{-0.23}^{+0.26}\times 10^{12}M_\odot$, and Andromeda, $M_{\rm M31}=1.33_{-0.33}^{+0.39}\times 10^{12}M_\odot$, are consistent with previous estimates that neglect the impact of the LMC on the observed Hubble flow. We find a (total) LMC mass $M_{\rm LMC}=0.25_{-0.08}^{+0.09}\times 10^{12}M_\odot$, which is indicative of an extended dark matter halo and supports the scenario where this galaxy is just past its first pericentric approach. Consequently, these results suggest that the LMC may induce significant perturbations on the Galactic potential.


title: Channel Blockade in a Two-Path Triple-Quantum-Dot System

abstract: Electronic transport through a two-path triple-quantum-dot system with two source leads and one drain is studied. By separating the conductance of the two double dot paths, we are able to observe double dot and triple dot physics in transport and study the interaction between the paths. We observe channel blockade as a result of inter-channel Coulomb interaction. The experimental results are understood with the help of a theoretical model which calculates the parameters of the system, the stability regions of each state and the full dynamical transport in the triple dot resonances.


title: Observational Signatures of Planets in Protoplanetary Disks II: Spiral Arms Observed in Scattered Light Imaging Can be Induced by Planets

abstract: Using 3D global hydro simulations coupled with radiative transfer calculations, we study the appearance of density waves induced by giant planets in direct imaging observations at near infrared wavelengths. We find that a 6 MJ planet in a typical disk around a 1 M_sun star can produce prominent and detectable spiral arms both interior and exterior to its orbit. The inner arms have (1) two well separated arms in roughly m=2 symmetry, (2) exhibit ~10-15 degrees pitch angles, (3) ~180-270 degrees extension in the azimuthal direction, and (4) ~150% surface brightness enhancement, all broadly consistent with observed spiral arms in the SAO 206462 and MWC 758 systems. The outer arms cannot explain observations as they are too tightly wound given typical disk scale height. We confirm previous results that the outer density waves excited by a 1 MJ planet exhibit low contrast in the IR and are practically not detectable. We also find that 3D effects of the waves are important. Compared to isothermal models, density waves in adiabatic disks exhibit weaker contrast in surface density but stronger contrast in scattered light images, due to a more pronounced vertical structure in the former caused by shock heating. To drive observed pairs of arms with an external companion on a circular orbit, a massive planet, possibly a brown dwarf, is needed at around [r~0.7", PA~10 degrees] (position angle PA from north to east) in SAO 206462 and [r~0.6", PA~10 degrees] in MWC 758. Their existence may be confirmed by direct imaging planet searches.


title: The emergence of 4-cycles in polynomial maps over the extended integers

abstract: Let $f(x) \in \mathbb{Z}[x]$; for each integer $\alpha$ it is interesting to consider the number of iterates $n_{\alpha}$, if possible, needed to satisfy $f^{n_{\alpha}}(\alpha) = \alpha$. The sets $\{\alpha, f(\alpha), \ldots, f^{n_{\alpha} - 1}(\alpha), \alpha\}$ generated by the iterates of $f$ are called cycles. For $\mathbb{Z}[x]$ it is known that cycles of length 1 and 2 occur, and no others. While much is known for extensions to number fields, we concentrate on extending $\mathbb{Z}$ by adjoining reciprocals of primes. Let $\mathbb{Z}[1/p_1, \ldots, 1/p_n]$ denote $\mathbb{Z}$ extended by adding in the reciprocals of the $n$ primes $p_1, \ldots, p_n$ and all their products and powers with each other and the elements of $\mathbb{Z}$. Interestingly, cycles of length 4, called 4-cycles, emerge for polynomials in $\mathbb{Z}\left[1/p_1, \ldots, 1/p_n\right][x]$ under the appropriate conditions. The problem of finding criteria under which 4-cycles emerge is equivalent to determining how often a sum of four terms is zero, where the terms are $\pm 1$ times a product of elements from the list of $n$ primes. We investigate conditions on sets of primes under which 4-cycles emerge. We characterize when 4-cycles emerge if the set has one or two primes, and (assuming a generalization of the ABC conjecture) find conditions on sets of primes guaranteed not to cause 4-cycles to emerge.


title: Determining Optimal Test Functions for Bounding the Average Rank in Families of $L$-Functions

abstract: Given an $L$-function, one of the most important questions concerns its vanishing at the central point; for example, the Birch and Swinnerton-Dyer conjecture states that the order of vanishing there of an elliptic curve $L$-function equals the rank of the Mordell-Weil group. The Katz and Sarnak Density Conjecture states that this and other behavior is well-modeled by random matrix ensembles. This correspondence is known for many families when the test functions are suitably restricted. For appropriate choices, we obtain bounds on the average order of vanishing at the central point in families. In this note we report on progress in determining the optimal test functions for the various classical compact groups for different support restrictions, and discuss how this relates to improved rank bounds.


title: The Structure of Spiral Shocks Excited by Planetary-mass Companions

abstract: Direct imaging observations have revealed spiral structures in protoplanetary disks. Previous studies have suggested that planet-induced spiral arms cannot explain some of these spiral patterns, due to the large pitch angle and high contrast of the spiral arms in observations. We have carried out three dimensional (3-D) hydrodynamical simulations to study spiral wakes/shocks excited by young planets. We find that, in contrast with linear theory, the pitch angle of spiral arms does depend on the planet mass, which can be explained by the non-linear density wave theory. A secondary (or even a tertiary) spiral arm, especially for inner arms, is also excited by a massive planet. With a more massive planet in the disk, the excited spiral arms have larger pitch angle and the separation between the primary and secondary arms in the azimuthal direction is also larger. We also find that although the arms in the outer disk do not exhibit much vertical motion, the inner arms have significant vertical motion, which boosts the density perturbation at the disk atmosphere. Combining hydrodynamical models with Monte-Carlo radiative transfer calculations, we find that the inner spiral arms are considerably more prominent in synthetic near-IR images using full 3-D hydrodynamical models than images based on 2-D models assuming vertical hydrostatic equilibrium, indicating the need to model observations with full 3-D hydrodynamics. Overall, companion-induced spiral arms not only pinpoint the companion's position but also provide three independent ways (pitch angle, separation between two arms, and contrast of arms) to constrain the companion's mass.


title: WIMP Dark Matter and Unitarity-Conserving Inflation via a Gauge Singlet Scalar

abstract: A gauge singlet scalar with non-minimal coupling to gravity can drive inflation and later freeze out to become cold dark matter. We explore this idea by revisiting inflation in the singlet direction (S-inflation) and Higgs Portal Dark Matter in light of the Higgs discovery, limits from LUX and observations by Planck. We show that large regions of parameter space remain viable, so that successful inflation is possible and the dark matter relic abundance can be reproduced. Moreover, the scalar singlet can stabilise the electroweak vacuum and at the same time overcome the problem of unitarity-violation during inflation encountered by Higgs Inflation, provided the singlet is a real scalar. The 2-$\sigma$ Planck upper bound on $n_s$ imposes that the singlet mass is below 2 TeV, so that almost the entire allowed parameter range can be probed by XENON1T.


title: Detecting the Companions and Ellipsoidal Variations of RS CVn Primaries: II. omicron Draconis, a Candidate for Recent Low-Mass Companion Ingestion

abstract: To measure the stellar and orbital properties of the metal-poor RS CVn binary o Draconis (o Dra), we directly detect the companion using interferometric observations obtained with the Michigan InfraRed Combiner at Georgia State University's Center for High Angular Resolution Astronomy (CHARA) Array. The H-band flux ratio between the primary and secondary stars is the highest confirmed flux ratio (370 +/- 40) observed with long-baseline optical interferometry. These detections are combined with radial velocity data of both the primary and secondary stars, including new data obtained with the Tillinghast Reflector Echelle Spectrograph on the Tillinghast Reflector at the Fred Lawrence Whipple Observatory and the 2-m Tennessee State University Automated Spectroscopic Telescope at Fairborn Observatory. We determine an orbit from which we find model-independent masses and ages of the components (M_A = 1.35 +\- 0.05 M_Sun, M_B = 0.99 +\- 0.02 M_Sun, system age = 3.0 -\+ 0.5 Gyr). An average of a 23-year light curve of o Dra from the Tennessee State University Automated Photometric Telescope folded over the orbital period newly reveals eclipses and the quasi-sinusoidal signature of ellipsoidal variations. The modeled light curve for our system's stellar and orbital parameters confirm these ellipsoidal variations due to the primary star partially filling its Roche lobe potential, suggesting most of the photometric variations are not due to stellar activity (starspots). Measuring gravity darkening from the average light curve gives a best-fit of beta = 0.07 +\- 0.03, a value consistent with conventional theory for convective envelope stars. The primary star also exhibits an anomalously short rotation period, which, when taken with other system parameters, suggests the star likely engulfed a low-mass companion that had recently spun-up the star.


title: Geodesic motion in equal angular momenta Myers-Perry-AdS spacetimes

abstract: We study the geodesic motion of massive and massless test particles in the background of equally spinning Myers-Perry-anti-de Sitter (AdS) black holes in five dimensions. By adopting a coordinate system that makes manifest the cohomogeneity-1 property of these spacetimes, the equations of motion simplify considerably. This allows us to easily separate the radial motion from the angular part and to obtain solutions for angular trajectories in a compact closed form. For the radial motion we focus our attention on spherical orbits. In particular, we determine the timelike innermost stable circular orbits (ISCOs) for these asymptotically AdS spacetimes, as well as the location of null circular orbits. We find that the ISCO dives below the ergosurface for black holes rotating close to extremality and merges with the event horizon exactly at extremality, in analogy with the four-dimensional Kerr case. For sufficiently massive black holes in AdS there exists a spin parameter range in which the background spacetime is stable against superradiance and the ISCO lies inside the ergoregion. Our results for massless geodesics show that there are no stable circular null orbits outside the horizon, but there exist such orbits inside the horizon, as well as around over-extremal spacetimes, i.e., naked singularities. We also discuss how these orbits deform from the static to the rotating case.


title: Universality Classes of Scale Invariant Inflation

abstract: We investigate the inflationary implications of extensions of Poincare symmetry. The simplest constructions with local scale invariance lead to universal predictions: the spectral index is $n_s = 1-2/N$, in excellent agreement with Planck data, while the tensor-to-scalar ratio is determined by a free parameter to $r = 12 \alpha / N^2$. For the special value $\alpha=1$ one finds symmetry enhancement to the full conformal group. We show that these findings hold both for two-derivative scalar-tensor theories as well as higher-derivative gravity. Therefore scale invariance underlies a promising set of inflationary models.


title: Using chemical tagging to redefine the interface of the Galactic disk and halo

abstract: We present a chemical abundance distribution study in 14 $\alpha$, odd-Z, even-Z, light, and Fe-peak elements of approximately 3200 intermediate metallicity giant stars from the APOGEE survey. The main aim of our analysis is to explore the Galactic disk-halo transition region between -1.20 $<$ [Fe/H] $<$ -0.55 as a means to study chemical difference (and similarities) between these components. In this paper, we show that there is an $\alpha$-poor and $\alpha$-rich sequence within both the metal-poor and intermediate metallicity regions. Using the Galactic rest-frame radial velocity and spatial positions, we further separate our sample into the canonical Galactic components. We then studied the abundances ratios, of Mg, Ti, Si, Ca, O, S, Al, C+N, Na, Ni, Mn, V, and K for each of the components and found the following: (1) the $\alpha$-poor halo subgroup is chemically distinct in the $\alpha$-elements (particularly O, Mg, and S), Al, C+N, and Ni from the $\alpha$-rich halo, consistent with the literature confirming the existence of an $\alpha$-poor accreted halo population; (2) the canonical thick disk and halo are not chemically distinct in all elements indicating a smooth transition between the thick disk and halo; (3) a subsample of the $\alpha$-poor stars at metallicities as low as [Fe/H] $\sim$ -0.85 dex are chemically and dynamically consistent with the thin disk indicating that the thin disk may extend to lower metallicities than previously thought, and (4) that the location of the most metal-poor thin disk stars are consistent with a negative radial metallicity gradient. Finally, we used our analysis to suggest a new set of chemical abundance planes ([$\alpha$/Fe], [C+N/Fe], [Al/Fe], and [Mg/Mn]) that may be able to chemically label the Galactic components in a clean and efficient way independent of kinematics.


title: Connecting Galaxies with Halos Across Cosmic Time: Stellar mass assembly distribution modeling of galaxy statistics

abstract: In this work, I explore an empirically motivated model for investigating the relationship between galaxy stellar masses, star formation rates and their halo masses and mass accretion histories. The core statistical quantity in this model is the stellar mass assembly distribution, $P(dM_{*}/dt|\mathbf{X},a)$, which specifies the probability density distribution of stellar mass assembly rates given a set of halo properties $\mathbf{X}$ and epoch $a$. Predictions from this model are obtained by integrating the stellar mass assembly distribution (SMAD) over halo merger trees, easily obtained from modern, high-resolution $N$-body simulations. Further properties of the galaxies hosted by the halos can be obtained by post-processing the stellar mass assembly histories with stellar population synthesis models. In my particular example implementation of this model, I use the \citet{behroozi13a} constraint on the median stellar mass assembly rates of halos as a function of their mass and redshift to construct an example parameterization of $P(dM_{*}/dt|\mathbf{X},a)$. This SMAD is then integrated over individual halo mass accretion histories from $N$-body merger trees starting at z = 4, using simple rules to account for merging halos. I find that this a simple model can reproduce qualitatively the bimodal features of the low-redshift galaxy population, including the qualitative split in the two-point clustering as a function of specific star formation rate. These results indicate that models which directly couple halo and galaxy growth through simple efficiency functions can naturally predict the star formation rate bimodality in higher-order statistics of the galaxy field, such as its two-point correlations or galactic conformity signals.


title: Monte Carlo Neutrino Transport Through Remnant Disks from Neutron Star Mergers

abstract: We present Sedonu, a new open source, steady-state, special relativistic Monte Carlo (MC) neutrino transport code, available at bitbucket.org/srichers/sedonu. The code calculates the energy- and angle-dependent neutrino distribution function on fluid backgrounds of any number of spatial dimensions, calculates the rates of change of fluid internal energy and electron fraction, and solves for the equilibrium fluid temperature and electron fraction. We apply this method to snapshots from two-dimensional simulations of accretion disks left behind by binary neutron star mergers, varying the input physics and comparing to the results obtained with a leakage scheme for the case of a central black hole and a central hypermassive neutron star. Neutrinos are guided away from the densest regions of the disk and escape preferentially around 45 degrees from the equatorial plane. Neutrino heating is strengthened by MC transport a few scale heights above the disk midplane near the innermost stable circular orbit, potentially leading to a stronger neutrino-driven wind. Neutrino cooling in the dense midplane of the disk is stronger when using MC transport, leading to a globally higher cooling rate by a factor of a few and a larger leptonization rate by an order of magnitude. We calculate neutrino pair annihilation rates and estimate that an energy of 2.8e46 erg is deposited within 45 degrees of the symmetry axis over 300 ms when a central BH is present. Similarly, 1.9e48 erg is deposited over 3 s when an HMNS sits at the center, but neither estimate is likely to be sufficient to drive a GRB jet.


title: Multi-Dimensional Geometric Complexity in Urban Transportation Systems

abstract: Transportation networks serve as windows into the complex world of urban systems. By properly characterizing a road network, we can therefore better understand its encompassing urban system. This study offers a geometrical approach towards capturing inherent properties of urban road networks. It offers a robust and efficient methodology towards defining and extracting three relevant indicators of road networks: area, line, and point thresholds, through measures of their grid equivalents. By applying the methodology to 50 U.S. urban systems, we successfully observe differences between eastern versus western, coastal versus inland, and old versus young, cities. Moreover, we show that many socio-economic characteristics as well as travel patterns within urban systems are directly correlated with their corresponding area, line, and point thresholds.


title: High temperature magnetic stabilization of cobalt nanoparticles by an antiferromagnetic proximity effect

abstract: Thermal activation tends to destroy the magnetic stability of small magnetic nanoparticles, with crucial implications in ultra-high density recording among other applications. Here we demonstrate that low blocking temperature ferromagnetic (FM) Co nanoparticles (TB<70 K) become magnetically stable above 400 K when embedded in a high N\'eel temperature antiferromagnetic (AFM) NiO matrix. The origin of this remarkable TB enhancement is due to a magnetic proximity effect between a thin CoO shell (with low N\'eel temperature, TN; and high anisotropy, KAFM) surrounding the Co nanoparticles and the NiO matrix (with high TN but low KAFM). This proximity effect yields an effective AFM with an apparent TN beyond that of bulk CoO, and an enhanced anisotropy compared to NiO. In turn, the Co core FM moment is stabilized against thermal fluctuations via core-shell exchange-bias coupling, leading to the observed TB increase. Mean-field calculations provide a semi-quantitative understanding of this magnetic- proximity stabilization mechanism.


title: Seismic modelling of the rotating, slowly pulsating B-type star HD 21071

abstract: Interpretation of the oscillation spectrum of the slowly pulsating B-type star HD21071 is presented. We show that non-rotating models cannot account for the two highest amplitude frequencies and taking into account the effects of rotation is necessary. Rotating seismic models are constructed using various chemical compositions, opacity data, core overshooting parameters and rotational velocities. There are prospects for seismic modelling of SPB stars, even if no asymptotic pattern is observed in their oscillation spectra, provided an unambiguous mode identification is doable and the effects of rotation are properly included.


title: Motion of particles on a Four-Dimensional Asymptotically AdS Black Hole with Scalar Hair

abstract: Motivated by black hole solutions with matter fields outside their horizon, we study the effect of these matter fields in the motion of massless and massive particles. We consider as background a four-dimensional asymptotically AdS black hole with scalar hair. The geodesics are studied numerically and we discuss about the differences in the motion of particles between the four-dimensional asymptotically AdS black holes with scalar hair and their no-hair limit, that is, Schwarzschild AdS black holes. Mainly, we found that there are bounded orbits like planetary orbits in this background. However, the periods associated to circular orbits are modified by the presence of the scalar hair. Besides, we found that some classical tests such as perihelion precession, deflection of light and gravitational time delay have the standard value of general relativity plus a correction term coming from the cosmological constant and the scalar hair. Finally, we found a specific value of the parameter associated to the scalar hair, in order to explain the discrepancy between the theory and the observations, for the perihelion precession of Mercury and light deflection.


title: Sun-Like Magnetic Cycles in the Rapidly-Rotating Young Solar Analog HD 30495

abstract: A growing body of evidence suggests that multiple dynamo mechanisms can drive magnetic variability on different timescales, not only in the Sun but also in other stars. Many solar activity proxies exhibit a quasi-biennial ($\sim$2 year) variation, which is superimposed upon the dominant 11 year cycle. A well-characterized stellar sample suggests at least two different relationships between rotation period and cycle period, with some stars exhibiting long and short cycles simultaneously. Within this sample, the solar cycle periods are typical of a more rapidly rotating star, implying that the Sun might be in a transitional state or that it has an unusual evolutionary history. In this work, we present new and archival observations of dual magnetic cycles in the young solar analog HD 30495, an $\sim$1 Gyr-old G1.5V star with a rotation period near 11 days. This star falls squarely on the relationships established by the broader stellar sample, with short-period variations at $\sim$1.7 years and a long cycle of $\sim$12 years. We measure three individual long-period cycles and find durations ranging from 9.6-15.5 years. We find the short-term variability to be intermittent, but present throughout the majority of the time series, though its occurrence and amplitude are uncorrelated with the longer cycle. These essentially solar-like variations occur in a Sun-like star with more rapid rotation, though surface differential rotation measurements leave open the possibility of a solar equivalence.


title: Comparative Traffic Performance Analysis of Urban Transportation Network Structures

abstract: The network structure of an urban transportation system has a significant impact on its traffic performance. This study uses network indicators along with several traffic performance measures including speed, trip length, travel time, and traffic volume, to compare a selection of seven transportation networks with a variety of structures and under different travel demand conditions. The selected network structures are: modified linear, branch, grid, 3-directional grid, 1-ring web, 2-ring web, and radial. For the analysis, a base origin-destination matrix is chosen, to which different growth factors are applied in order to simulate various travel demand conditions. Results show that overall the 2-ring web network offers the most efficient traffic performance, followed by the grid and the 1-ring networks. A policy application of this study is that the branch, 3-directional grid, and radial networks are mostly suited for small cities with uncongested traffic conditions. In contrast, the 2-ring web, grid, and 1-ring web networks are better choices for large urban areas since they offer more connectivity, thus allowing them to perform efficiently under congested traffic conditions.


title: Continuous matrix product states solution for the mixing/demixing transition in one-dimensional quantum fields

abstract: We solve the mixing-demixing transition in repulsive one-dimensional bose-bose mixtures. This is done numerically by means of the continuous matrix product states variational ansatz. We show that the effective low-energy bosonization theory is able to detect the transition whenever the Luttinger parameters are exactly computed. We further characterize the transition by calculating the ground-state energy density, the field-field fluctuations and the density correlations.


title: Flexible and Low-Complexity Encoding and Decoding of Systematic Polar Codes

abstract: In this work, we present hardware and software implementations of flexible polar systematic encoders and decoders. The proposed implementations operate on polar codes of any length less than a maximum and of any rate. We describe the low-complexity, highly parallel, and flexible systematic-encoding algorithm that we use and prove its correctness. Our hardware implementation results show that the overhead of adding code rate and length flexibility is little, and the impact on operation latency minor compared to code-specific versions. Finally, the flexible software encoder and decoder implementations are also shown to be able to maintain high throughput and low latency.


title: Distance Preserving Graphs

abstract: Given a graph $G$ then a subgraph $H$ is $isometric$ if, for every pair of vertices $u,v$ of $H$, we have $d_H(u,v) = d_G(u,v)$. We say a graph $G$ is $distance\ preserving\ (dp)$ if it has an isometric subgraph of every possible order up to the order of $G$. We consider how to add a vertex to a dp graph so that the result is a dp graph. This condition implies that chordal graphs are dp. We also find a condition on the girth of $G$ which implies that it is not dp. In closing, we discuss other work and open problems concerning dp graphs.


title: Multichannel long-range Rydberg molecules

abstract: A generalized class of ultra-long-range Rydberg molecules is proposed which consist of a multichannel Rydberg atom whose outermost electron creates a chemical bond with a distant ground state atom. Such multichannel Rydberg molecules exhibit favorable properties for laser excitation, because states exist where the quantum defect varies strongly with the principal quantum number. The resulting occurrence of near degeneracies with states of high orbital angular momentum promotes the admixture of low $l$ into the high $l$ deeply bound `trilobite' molecule states, thereby circumventing the usual difficulty posed by electric dipole selection rules. Such states also can exhibit multi-scale binding possibilities that could present novel options for quantum manipulation.


title: Zero-one law for directional transience of one-dimensional random walks in dynamic random environments

abstract: We prove the trichotomy between transience to the right, transience to the left and recurrence of one-dimensional nearest-neighbour random walks in dynamic random environments under fairly general assumptions, namely: stationarity under space-time translations, ergodicity under spatial translations, and a mild ellipticity condition. In particular, the result applies to general uniformly elliptic models and also to a large class of non-uniformly elliptic cases that are i.i.d. in space and Markovian in time. An immediate consequence is the recurrence of models that are symmetric with respect to reflection through the origin.


title: Higgs potential and hidden light Higgs scenario in two Higgs doublet models

abstract: In two Higgs doublet models (2HDM), there exists an interesting possibility, the hidden light Higgs scenario, that the discovered SM-like Higgs boson is the heavier \textit{CP}-even Higgs boson $H^0$ and the lighter \textit{CP}-even $h^0$ has not been observed yet in any experiment. We study the current status of this scenario in Types I, II, X, and Y, through the scans of the 2HDM parameters with all relevant theoretical and experimental constraints. We employ not only the most up-to-date Higgs signal strength measurements with the feed-down effects, but also all the available LHC exclusion limits from heavy Higgs searches. Adjusting the heavier $H^0$ to the 125 GeV state while hiding the lighter $h^0$ from the LEP Higgs search prohibits the extreme decoupling limit: there exist upper bounds on the masses of the pseudoscalar $A^0$ and the charged Higgs $H^\pm$ below about 600 GeV. In addition, the $Z_2$ symmetry is shown to be a good approximate symmetry since the soft $Z_2$ symmetry breaking parameter $m_{12}^2$ should be less than about (45 GeV)$^2$. Most interestingly, a few parameters in the Higgs potential and the related Higgs triple and quartic couplings are shown to be meaningfully constrained by the current data. The double Higgs-strahlung process at an $e^+ e^-$ collider is also studied.


title: Comparison theorems for causal diamonds

abstract: We formulate certain inequalities for the geometric quantities characterizing causal diamonds in curved and Minkowski spacetimes. These inequalities involve the red-shift factor which, as we show explicitly in the spherically symmetric case, is monotonic in the radial direction and it takes its maximal value at the centre. As a byproduct of our discussion we re-derive Bishop's inequality without assuming the positivity of the spatial Ricci tensor. We then generalize our considerations to arbitrary, static and not necessarily spherically symmetric, asymptotically flat spacetimes. In the case of spacetimes with a horizon our generalization involves the so-called {\it domain of dependence}. The respective volume, expressed in terms of the duration measured by a distant observer compared with the volume of the domain in Minkowski spacetime, exhibits behaviours which differ if $d=4$ or $d>4$. This peculiarity of four dimensions is due to the logarithmic subleading term in the asymptotic expansion of the metric near infinity. In terms of the invariant duration measured by a co-moving observer associated with the diamond we establish an inequality which is universal for all $d$. We suggest some possible applications of our results including the comparison theorems for entanglement entropy, causal set theory and fundamental limits on computation.


title: On the conformal properties of topological terms in even dimensions

abstract: Conformal properties of the topological gravitational terms in $D=2$, $D=4$ and $D=6$ are discussed. It is shown that in the last two cases the integrands of these terms, when being settled into the dimension $D-1$ and multiplied by a scalar, become conformal invariant. Furthermore we present a simple covariant derivation of the Paneitz operator in $D=4$ and formulate two general conjectures concerning the conformal properties of topological structures in even dimensions.


title: Reactive Boundary Conditions as Limits of Interaction Potentials for Brownian and Langevin Dynamics

abstract: A popular approach to modeling bimolecular reactions between diffusing molecules is through the use of reactive boundary conditions. One common model is the Smoluchowski partial absorption condition, which uses a Robin boundary condition in the separation coordinate between two possible reactants. This boundary condition can be interpreted as an idealization of a reactive interaction potential model, in which a potential barrier must be surmounted before reactions can occur. In this work we show how the reactive boundary condition arises as the limit of an interaction potential encoding a steep barrier within a shrinking region in the particle separation, where molecules react instantly upon reaching the peak of the barrier. The limiting boundary condition is derived by the method of matched asymptotic expansions, and shown to depend critically on the relative rate of increase of the barrier height as the width of the potential is decreased. Limiting boundary conditions for the same interaction potential in both the overdamped Fokker-Planck equation (Brownian Dynamics), and the Kramers equation (Langevin Dynamics) are investigated. It is shown that different scalings are required in the two models to recover reactive boundary conditions that are consistent in the high friction limit where the Kramers equation solution converges to the solution of the Fokker-Planck equation.


title: Transits of Venus and Solar diameter measures from ground: method and results from Athens (2004) and Huairou (2012)

abstract: The variation of the solar diameter in time and in position angle has implications in astrophysics and in general relativity, as the long series of studies attest. The Transits of Venus in 2004 and 2012 have been carefully studied because of the rarity of the phenomenon and its historical importance due the AU measure and to the discovery of Venus atmosphere. The characterization of Venus atmosphere and the measure of the solar diameter to the milliarcsecond level of precision have been studied also from satellite images. The results of the solar diameter measurements made with the observations in Athens (2004) and at the Huairou Solar Observing Station in China (2012) are presented. The topic of the oblateness of the Sun at sunset and its intrinsic value is drafted to introduce the general public to the relativistic relevance of measuring the solar figure, in the occasion of the International Year of Light 2015.


title: Two characterizations of simple circulant tournaments

abstract: The \textit{acyclic disconnection} $\overrightarrow{\omega }(D)$ (resp. the \textit{directed triangle free disconnection } $\overrightarrow{\omega }_{3}(D)$) of a digraph $D$ is defined as the maximum possible number of connected components of the underlying graph of $D\setminus A(D^{\ast })$ where $D^{\ast }$ is an acyclic (resp. a directed triangle free) subdigraph of $D$. In this paper, we generalize some previous results and solve some problems posed by V. Neumann-Lara (The acyclic disconnection of a digraph, Discrete Math. 197/198 (1999), 617-632). Let $\overrightarrow{C}_{2n+1}(J)$ be a circulant tournament. We prove that $\overrightarrow{C}_{2n+1}(J)$ is $\overrightarrow{% \omega }$-keen and $\overrightarrow{\omega _{3}}$-keen, respectively, and $% \overrightarrow{\omega }(\overrightarrow{C}_{2n+1}(J))=\overrightarrow{% \omega }_{3}(\overrightarrow{C}_{2n+1}(J))=2$ for every $\overrightarrow{C}% _{2n+1}(J)$. Finally, it is showed that $\overrightarrow{\omega }_{3}(% \overrightarrow{C}_{2n+1}(J))=2$, $\overrightarrow{C}_{2n+1}(J)$ is simple and $J$ is aperiodic are equivalent propositions.


title: Scaling limit for a long-range divisible sandpile

abstract: We study the scaling limit of a divisible sandpile model associated to a truncated $\alpha$-stable random walk. We prove that the limiting distribution is related to an obstacle problem for a truncated fractional Laplacian. We also provide, as a fundamental tool, precise asymptotic expansions for the corresponding rescaled discrete Green's functions. In particular, the convergence rate of these Green's functions to its continuous counterpart is derived.


title: The Langlands-Shahidi method over function fields: Ramanujan Conjecture and Riemann Hypothesis for the unitary groups

abstract: We establish the Langlands-Shahidi method over a global field of characteristic p. We then focus on the unitary groups and prove global and local Langlands functoriality to general linear groups for generic representations. Main applications are to the Ramanujan Conjecture and Riemann Hypothesis.


title: Nonlinear absorption of high-intensity shortwave radiation in plasma within relativistic quantum theory

abstract: On the base of the quantum kinetic equation for density matrix in plasma at the stimulated bremsstrahlung of electrons on ions, the nonlinear absorption rate for high-intensity shortwave radiation in plasma has been obtained within relativistic quantum theory. Both classical Maxwellian and degenerate quantum plasma are considered for x-ray lasers of high intensities. Essentially different dependences of nonlinear absorption rate on polarization of strong laser radiation is stated.


title: The Hardy inequality and the heat flow in curved wedges

abstract: We show that the polynomial decay rate of the heat semigroup of the Dirichlet Laplacian in curved planar wedges equals the sum of the usual dimensional decay rate and a multiple of the reciprocal value of the opening angle. To prove the result, we develop the method of self-similar variables for the associated heat equation and study the asymptotic behaviour of the transformed non-autonomous parabolic problem for large times. We also establish an improved Hardy inequality for the Dirichlet Laplacian in non-trivially curved wedges and state a conjecture about an improved decay rate in this case.


title: Orbital and Pauli limiting effects in heavily doped Ba$_{1-x}$K$_x$Fe$_2$As$_2$

abstract: We investigated the thermodynamic properties of the Fe-based lightly disordered superconductor Ba$_{0.05}$K$_{0.95}$Fe$_2$As$_2$ in external magnetic field H applied along the FeAs layers (H//ab planes). The superconducting (SC) transition temperature for this doping level is T$_c$ = 6.6 K. Our analysis of the specific heat C(T,H) measured for T < T$_c$ implies a sign change of the superconducting order parameter across different Fermi pockets. We provide experimental evidence for the three components superconducting order parameter. We find that all three components have values which are comparable with the previously reported ones for the stochiometric compound KFe$_2$As$_2$. Our data for C(T,H) and resistivity rho(T,H) can be interpreted in favor of the dominant orbital contribution to the pair-breaking mechanism at low fields, while Pauli limiting effect dominates at high fields, giving rise to a gapless superconducting state with only the leading non-zero gap.


title: Magnetic reconnection in plasma under inertial confinement fusion conditions driven by heat flux effects in Ohm's law

abstract: In the interaction of high-power laser beams with solid density plasma there are a number of mechanisms that generate strong magnetic fields. Such fields subsequently inhibit or redirect electron flows, but can themselves be advected by heat fluxes, resulting in complex interplay between thermal transport and magnetic fields.We show that for heating by multiple laser spots reconnection of magnetic field lines can occur, mediated by these heat fluxes, using a fully implicit 2D Vlasov-Fokker-Planck code. Under such conditions, the reconnection rate is dictated by heat flows rather than Alfv\`enic flows. We find that this mechanism is only relevant in a high $\beta$ plasma. However, the Hall parameter $\omega_c \tau_{ei}$ can be large so that thermal transport is strongly modified by these magnetic fields, which can impact longer time scale temperature homogeneity and ion dynamics in the system.


title: Radiative corrections to the Triple Higgs Coupling in the Inert Higgs Doublet Model

abstract: We investigate the implication of the recent discovery of a Higgs-like particle in the first phase of the LHC Run 1 on the Inert Higgs Doublet Model (IHDM). The determination of the Higgs couplings to SM particles and its intrinsic properties will get improved during the new LHC Run 2 starting this year. The new LHC Run 2 would also shade some light on the triple Higgs coupling. Such measurement is very important in order to establish the details of the electroweak symmetry breaking mechanism. Given the importance of the Higgs couplings both at the LHC and $e^+e^-$ Linear Collider machines, accurate theoretical predictions are required. We study the radiative corrections to the triple Higgs coupling $hhh$ and to $hZZ$, $hWW$ couplings in the context of the IHDM. By combining several theoretical and experimental constraints on parameter space, we show that extra particles might modify the triple Higgs coupling near threshold regions. Finally, we discuss the effect of these corrections on the double Higgs production signal at the $e^+e^-$ LC and show that they can be rather important.


title: A survey on the kissing numbers

abstract: The maximum possible number of non-overlapping unit spheres that can touch a unit sphere in $n$ dimensions is called kissing number. The problem for finding kissing numbers is closely connected to the more general problems of finding bounds for spherical codes and sphere packings. We survey old and recent results on the kissing numbers keeping the generality of spherical codes.


title: On Recurrent Reachability for Continuous Linear Dynamical Systems

abstract: The continuous evolution of a wide variety of systems, including continuous-time Markov chains and linear hybrid automata, can be described in terms of linear differential equations. In this paper we study the decision problem of whether the solution $\boldsymbol{x}(t)$ of a system of linear differential equations $d\boldsymbol{x}/dt=A\boldsymbol{x}$ reaches a target halfspace infinitely often. This recurrent reachability problem can equivalently be formulated as the following Infinite Zeros Problem: does a real-valued function $f:\mathbb{R}_{\geq 0}\rightarrow\mathbb{R}$ satisfying a given linear differential equation have infinitely many zeros? Our main decidability result is that if the differential equation has order at most $7$, then the Infinite Zeros Problem is decidable. On the other hand, we show that a decision procedure for the Infinite Zeros Problem at order $9$ (and above) would entail a major breakthrough in Diophantine Approximation, specifically an algorithm for computing the Lagrange constants of arbitrary real algebraic numbers to arbitrary precision.


title: The First Fermi-LAT SNR Catalog SNR and Cosmic Ray Implications

abstract: Galactic cosmic ray (CRs) sources, classically proposed to be Supernova Remnants (SNRs), must meet the energetic particle content required by direct measurements of high energy CRs. Indirect gamma-ray measurements of SNRs with the Fermi Large Area Telescope (LAT) have now shown directly that at least three SNRs accelerate protons. With the first Fermi LAT SNR Catalog, we have systematically characterized the GeV gamma-rays emitted by 279 SNRs known primarily from radio surveys. We present these sources in a multiwavelength context, including studies of correlations between GeV and radio size, flux, and index, TeV index, and age and environment tracers, in order to better understand effects of evolution and environment on the GeV emission. We show that previously sufficient models of SNRs' GeV emission no longer adequately describe the data. To address the question of CR origins, we also examine the SNRs' maximal CR contribution assuming the GeV emission arises solely from proton interactions. Improved breadth and quality of multiwavelength data, including distances and local densities, and more, higher resolution gamma-ray data with correspondingly improved Galactic diffuse models will strengthen this constraint.


title: Idempotents in intensional type theory

abstract: We study idempotents in intensional Martin-L\"of type theory, and in particular the question of when and whether they split. We show that in the presence of propositional truncation and Voevodsky's univalence axiom, there exist idempotents that do not split; thus in plain MLTT not all idempotents can be proven to split. On the other hand, assuming only function extensionality, an idempotent can be split if and only if its witness of idempotency satisfies one extra coherence condition. Both proofs are inspired by parallel results of Lurie in higher category theory, showing that ideas from higher category theory and homotopy theory can have applications even in ordinary MLTT. Finally, we show that although the witness of idempotency can be recovered from a splitting, the one extra coherence condition cannot in general; and we construct "the type of fully coherent idempotents", by splitting an idempotent on the type of partially coherent ones. Our results have been formally verified in the proof assistant Coq.


title: On a monodromy theorem for sheaves of local fields and applications

abstract: We prove a monodromy theorem for local vector fields belonging to a sheaf satisfying the unique continuation property. In particular, in the case of admissible regular sheaves of local fields defined on a simply connected manifold, we obtain a global extension result for every local field of the sheaf. This generalizes previous works of Nomizu for semi-Riemannian Killing fields, of Ledger--Obata for conformal fields, and of Amores for fields preserving a $G$-structure of finite type. The result applies to Finsler or pseudo-Finsler Killing fields and, more generally, to affine fields of a spray. Some applications are discussed.


title: Optical Deformations in Solar Glass Filters for High Precision Astrometry

abstract: Measuring the solar diameter at all position angles gives the complete figure of the Sun. Their asphericities have implications in classical physics and general relativity, and the behavior of the optical systems used in the direct measurements is to be known accurately. A solar filter is a plane-parallel glass with given absorption, and here we study the departures from the parallelism of the faces of a crystal slab 5 mm thick, because of static deformations. These deformations are rescaled to the filter's dimensions. Related to the Solar Disk Sextant experiment and to the Reflecting Heliometer of Rio de Janeiro a simplified model of the influences of the inclination between the external and the internal surfaces of a glass solar filter, is discussed.


title: Spin-orbit interaction driven collective electron-hole excitations in a noncentrosymmetric nodal loop Weyl semimetal

abstract: NbP is one member of a new class of nodal loop semimetals characterized by the cooperative effects of spin-orbit coupling (SOC) and a lack of inversion center. Here transport and spectroscopic properties of NbP are evaluated using density functional theory methods. SOC together with the lack of inversion symmetry splits degeneracies, giving rise to "Russian doll nested" Fermi surfaces containing 4*10$^{-4}$ electron (hole) carriers/f.u. Due to the modest SOC strength in Nb, the Fermi surfaces map out the Weyl nodal loops. Calculated structure around T$^*$~100 K in transport properties reproduces well the observed transport behavior only when SOC is included, attesting to the precision of the (delicate) calculations and the stoichiometry of the samples. Low energy collective electron-hole excitations (plasmons) in the 20-60 meV range result from the nodal loop splitting.


title: Experimental analysis of data-driven control for a building heating system

abstract: Driven by the opportunity to harvest the flexibility related to building climate control for demand response applications, this work presents a data-driven control approach building upon recent advancements in reinforcement learning. More specifically, model assisted batch reinforcement learning is applied to the setting of building climate control subjected to a dynamic pricing. The underlying sequential decision making problem is cast on a markov decision problem, after which the control algorithm is detailed. In this work, fitted Q-iteration is used to construct a policy from a batch of experimental tuples. In those regions of the state space where the experimental sample density is low, virtual support samples are added using an artificial neural network. Finally, the resulting policy is shaped using domain knowledge. The control approach has been evaluated quantitatively using a simulation and qualitatively in a living lab. From the quantitative analysis it has been found that the control approach converges in approximately 20 days to obtain a control policy with a performance within 90% of the mathematical optimum. The experimental analysis confirms that within 10 to 20 days sensible policies are obtained that can be used for different outside temperature regimes.


title: Splitting of ISGMR strength in the light-mass nucleus $^{24}$Mg due to ground-state deformation

abstract: The isoscalar giant monopole resonance (ISGMR) strength distribution in $^{24}$Mg has been determined from background-free inelastic scattering of 386-MeV $\alpha$ particles at extreme forward angles, including 0$^{\circ}$. The ISGMR strength distribution has been observed for the first time to have a two-peak structure in a light-mass nucleus. This splitting of ISGMR strength is explained well by microscopic theory in terms of the prolate deformation of the ground state of $^{24}$Mg.


title: Gyromagnetic g_s factors of the spin-1/2 particles in the (1/2+,1/2-,3/2-) triad of the four-vector spinor, \psi_\mu, irreducibility, and linearity

abstract: We show that the spin (1/2-) particle from the (1/2,1)+(1,1/2) Lorentz irreducible sector of the four-vector spinor can not be described within a linear formalism but behaves as a genuinely quadratic fermion satisfying the generalized Feynman-Gell-Mann equation with a gyromagnetic factor of (-2/3). In contrast, spin (1/2 +) from the (1/2,0)+(0,1/2) sector is confirmed as a genuine linear Dirac fermion whose gyromagnetic factor takes the value of two units.We calculate Compton scatterings off each one of the two targets and obtain both times well behaved cross sections in the ultra relativistic limit and in accord with unitarity.


title: Neural CRF Parsing

abstract: This paper describes a parsing model that combines the exact dynamic programming of CRF parsing with the rich nonlinear featurization of neural net approaches. Our model is structurally a CRF that factors over anchored rule productions, but instead of linear potential functions based on sparse features, we use nonlinear potentials computed via a feedforward neural network. Because potentials are still local to anchored rules, structured inference (CKY) is unchanged from the sparse case. Computing gradients during learning involves backpropagating an error signal formed from standard CRF sufficient statistics (expected rule counts). Using only dense features, our neural CRF already exceeds a strong baseline CRF model (Hall et al., 2014). In combination with sparse features, our system achieves 91.1 F1 on section 23 of the Penn Treebank, and more generally outperforms the best prior single parser results on a range of languages.


title: On the number of open knight's tours

abstract: We review the state of the art in the problem of counting the number open knight tours, since the publication in internet of a computation of this quantity.


title: The APOSTLE project: Local Group kinematic mass constraints and simulation candidate selection

abstract: We use a large sample of isolated dark matter halo pairs drawn from cosmological N-body simulations to identify candidate systems whose kinematics match that of the Local Group of Galaxies (LG). We find, in agreement with the "timing argument" and earlier work, that the separation and approach velocity of the Milky Way (MW) and Andromeda (M31) galaxies favour a total mass for the pair of $\sim 5\times 10^{12} \,M_{\odot}$. A mass this large, however, is difficult to reconcile with the small relative tangential velocity of the pair, as well as with the small deceleration from the Hubble flow observed for the most distant LG members. Halo pairs that match these three criteria have average masses a factor of $\sim 2$ times smaller than suggested by the timing argument, but with large dispersion. Guided by these results, we have selected $12$ halo pairs with total mass in the range $1.6$-$3.6 \times 10^{12}\,M_{\odot}$ for the APOSTLE project (A Project Of Simulating The Local Environment), a suite of hydrodynamical resimulations at various numerical resolution levels (reaching up to $\sim10^{4}\,M_{\odot}$ per gas particle) that use the subgrid physics developed for the EAGLE project. These simulations reproduce, by construction, the main kinematics of the MW-M31 pair, and produce satellite populations whose overall number, luminosities, and kinematics are in good agreement with observations of the MW and M31 companions. The APOSTLE candidate systems thus provide an excellent testbed to confront directly many of the predictions of the $\Lambda$CDM cosmology with observations of our local Universe.


title: Three PT-symmetric Hamiltonians with completely different spectra

abstract: We discuss three Hamiltonians, each with a central-field part $H_{0}$ and a PT-symmetric perturbation $igz$. When $H_{0}$ is the isotropic Harmonic oscillator the spectrum is real for all $g$ because $H$ is isospectral to $H_{0}+g^{2}/2$. When $H_{0}$ is the Hydrogen atom then infinitely many eigenvalues are complex for all $g$. If the potential in $H_{0}$ is linear in the radial variable $r$ then the spectrum of $H$ exhibits real eigenvalues for $0<g<g_{c}$ and a PT phase transition at $g_c$.


title: Magnetar Driven Shock Breakout and Double Peaked Supernova Light Curves

abstract: The light curves of some luminous supernovae are suspected to be powered by the spindown energy of a rapidly rotating magnetar. Here we describe a possible signature of the central engine: a burst of shock breakout emission occurring several days after the supernova explosion. The energy input from the magnetar inflates a high-pressure bubble that drives a shock through the pre-exploded supernova ejecta. If the magnetar is powerful enough, that shock will near the ejecta surface and become radiative. At the time of shock breakout, the ejecta will have expanded to a large radius (~10^{14} cm) so that the radiation released is at optical/ultraviolet wavelengths (T ~ 20,000 K) and lasts for several days. The luminosity and timescale of this magnetar driven shock breakout are similar to the first peak observed recently in the double-peaked light curve of SNLSQ14BDQ. However, for a large region of model parameter space, the breakout emission is predicted to be dimmer than the diffusive luminosity from direct magnetar heating. A distinct double peaked light curve may therefore only be conspicuous if thermal heating from the magnetar is suppressed at early times. We describe how such a delay in heating may naturally result from inefficient dissipation and thermalization of the pulsar wind magnetic energy. Without such suppression, the breakout may only be noticeable as a small bump or kink in the early luminosity or color evolution, or as a small but abrupt rise in the photospheric velocity. A similar breakout signature may accompany other central engines in supernovae, such as a black hole accreting fallback material.


title: Astronomical Data Formats: What we have and how we got here

abstract: Despite almost all being acquired as photons, astronomical data from different instruments and at different stages in its life may exist in different formats to serve different purposes. Beyond the data itself, descriptive information is associated with it as metadata, either included in the data format or in a larger multi-format data structure. Those formats may be used for the acquisition, processing, exchange, and archiving of data. It has been useful to use similar formats, or even a single standard to ease interaction with data in its various stages using familiar tools. Knowledge of the evolution and advantages of present standards is useful before we discuss the future of how astronomical data is formatted. The evolution of the use of world coordinates in FITS is presented as an example.


title: The shape of the one-dimensional phylogenetic likelihood function

abstract: By fixing all parameters in a phylogenetic likelihood model except for one branch length, one obtains a one-dimensional likelihood function. In this work, we introduce a mathematical framework to characterize the shapes of such one-dimensional phylogenetic likelihood functions. This framework is based on analyses of algebraic structures on the space of all frequency patterns with respect to a polynomial representation of the likelihood functions. Using this framework, we provide conditions under which the one-dimensional phylogenetic likelihood functions are guaranteed to have at most one stationary point, and this point is the maximum likelihood branch length. These conditions are satisfied by common simple models including all binary models, the Jukes-Cantor model and the Felsenstein 1981 model. We then prove that for the simplest model that does not satisfy our conditions, namely, the Kimura 2-parameter model, the one-dimensional likelihood functions may have multiple stationary points. As a proof of concept, we construct a non-degenerate example in which the phylogenetic likelihood function has two local maxima and a local minimum. To construct such examples, we derive a general method of constructing a tree and sequence data with a specified frequency pattern at the root. We then extend the result to prove that the space of all rescaled and translated one-dimensional phylogenetic likelihood functions under the Kimura 2-parameter model is dense in the space of all non-negative continuous functions on $[0, \infty)$ with finite limits. These results indicate that one-dimensional likelihood functions under advanced evolutionary models can be more complex than it is typically assumed by phylogenetic inference algorithms; however, these complexities can be effectively captured by the Kimura 2-parameter model.


title: Deadline-aware Power Management in Data Centers

abstract: We study the dynamic power optimization problem in data centers. We formulate and solve the following offline problem: in which slot which server has to be assigned to which job; and in which slot which server has to be switched ON or OFF so that the total power is optimal for some time horizon. We show that the offline problem is a new version of generalized assignment problem including new constraints issuing from deadline characteristics of jobs and difference of activation energy of servers. We propose an online algorithm that solves the problem heuristically and compare it to randomized routing.


title: Effect of shear and magnetic field on the heat-transfer efficiency of convection in rotating spherical shells

abstract: We study rotating thermal convection in spherical shells. We base our analysis on a set of about 450 direct numerical simulations of the (magneto)hydrodynamic equations under the Boussinesq approximation. The Ekman number ranges from $10^{-3}$ to $10^{-5}$. The supercriticality of the convection reaches about 1000 in some models. Four sets of simulations are considered: non-magnetic simulations and dynamo simulations with either free-slip or no-slip flow boundary conditions. The non-magnetic setup with free-slip boundaries generates the strongest zonal flows. Both non-magnetic simulations with no-slip flow boundary conditions and self-consistent dynamos with free-slip boundaries have drastically reduced zonal-flows. Suppression of shear leads to a substantial gain in heat-transfer efficiency, increasing by a factor of 3 in some cases. Such efficiency enhancement occurs as long as the convection is significantly influenced by rotation. At higher convective driving the heat-transfer efficiency tends towards that of the classical non-rotating Rayleigh-B\'enard system. Analysis of the latitudinal distribution of heat flow at the outer boundary reveals that the shear is most effective at suppressing heat-transfer in the equatorial regions. Furthermore, we explore the influence of the magnetic field on the {\em non-zonal} flow components of the convection. For this we compare the heat-transfer efficiency of no-slip non-magnetic cases with that of the no-slip dynamo simulations. We find that at $E=10^{-5}$ magnetic field significantly affects the convection and a maximum gain of about 30\% (as compared to the non-magnetic case) in heat-transfer efficiency is obtained for an Elsasser number of about 3. Our analysis motivates us to speculate that convection in the polar regions in dynamos at $E=10^{-5}$ is probably in a `magnetostrophic' regime.

categories: astro-ph.EP astro-ph.SR physics.flu-dyn physics.geo-ph physics.plasm-ph