title: A linear programming based approach for determining maximal closest reference set in DEA

abstract: Identification of the reference set for each decision making unit (DMU) is a main concern in the data envelopment analysis (DEA). All of the methods developed to date have been focused on finding the furthest reference DMUs. In this paper, we introduce the new notion of maximal closest reference set (MCRS) containing the maximum number of closest reference DMUs to the assessed DMU. Then, we develop a linear programming (LP) model for determining the MCRS for each inefficient DMU. We illustrate our method through a numerical example.


title: Determining Maximal Reference Set in Data Envelopment Analysis

abstract: In data envelopment analysis (DEA), the occurrence of multiple reference sets is a crucial issue in identifying all the reference DMUs to a given decision making unit (DMU). To resolve this difficulty, we introduce the useful notion of maximal reference set (MRS) which contains all the reference DMUs. Based on both primal and dual formulations, we then propose a new mixed integer linear programming (MILP) based approach for determining the MRS. The proposed approach is more general than the existing ones, and has several desirable properties: (i) having a unified formulation, (ii) adaptable with different DEA models, and (iii) compatible with both primal and dual forms of the DEA models. Furthermore, our approach is made computationally effective by establishing a LP-based procedure to treat with the primal-based MILP.


title: The dynamics of eruptive prominences

abstract: This chapter discusses the dynamical properties of eruptive prominences in relation to coronal mass ejections (CMEs). The fact that eruptive prominences are a part of CMEs is emphasized in terms of their physical association and kinematics. The continued propagation of prominence material into the heliosphere is illustrated using in-situ observations. The solar-cycle variation of eruptive prominence locations is discussed with a particular emphasis on the rush-to-the-pole (RTTP) phenomenon. One of the consequences of the RTTP phenomenon is polar CMEs, which are shown to be similar to the low-latitude CMEs. This similarity is important because it provides important clues to the mechanism by which CMEs erupt. The nonradial motion of CMEs is discussed, including the deflection by coronal holes that have important space weather consequences. Finally, the implications of the presented observations for the modeling CME modeling are outlined.


title: Nonlinear drift-diffusion model of gating in the fast Cl channel

abstract: The dynamics of the open or closed state region of an ion channel may be described by a probability density $p(x,t)$ which satisfies a Fokker-Planck equation. The closed state dwell-time distribution $f_c(t)$ derived from the Fokker-Planck equation with a nonlinear diffusion coefficient $D(x) \propto \exp(-\gamma x)$, $\gamma > 0$ and a linear ramp potential $U_c(x)$, is in good agreement with experimental data and it may be shown analytically that if $\gamma$ is sufficiently large, $f_c(t) \propto t^{-2 - \nu}$ for intermediate times, where $\nu = U_c^{\prime}/\gamma \approx -0.3$ for a fast Cl channel. The solution of a master equation which approximates the Fokker-Planck equation exhibits an oscillation superimposed on the power law trend and can account for an empirical rate-amplitude correlation that applies to several ion channels.


title: Values of symmetric polynomials and a truncated analogue of the Riemann zeta function

abstract: For each positive integer n, we determine the set of symmetric functions f for which the congruence f(p/1,p/2,...,p/(p-1)) \equiv 0 mod p^n holds for all sufficiently large primes p. Our determination is conditional on a conjecture regarding the modulo p independence of Bernoulli numbers. In a recent work the author introduced a new truncated analogue of the multiple zeta function and investigated a class of relations among values of this function at positive integers. The question answered in the present work is equivalent to the determination of the relations satisfied by values of the corresponding analogue of the ordinary Riemann zeta function.


title: Universality conjecture and results for a model of several coupled positive-definite matrices

abstract: The paper contains two main parts: in the first part, we analyze the general case of $p\geq 2$ matrices coupled in a chain subject to Cauchy interaction. Similarly to the Itzykson-Zuber interaction model, the eigenvalues of the Cauchy chain form a multi level determinantal point process. We first compute all correlations functions in terms of Cauchy biorthogonal polynomials and locate them as specific entries of a $(p+1)\times (p+1)$ matrix valued solution of a Riemann-Hilbert problem. In the second part, we fix the external potentials as classical Laguerre weights. We then derive strong asymptotics for the Cauchy biorthogonal polynomials when the support of the equilibrium measures contains the origin. As a result, we obtain a new family of universality classes for multi-level random determinantal point fields which include the Bessel$_\nu$ universality for $1$-level and the Meijer-$G$ universality for $2$-level. Our analysis uses the Deift-Zhou nonlinear steepest descent method and the explicit construction of a $(p+1)\times (p+1)$ origin parametrix in terms of Meijer G-functions. The solution of the full Riemann-Hilbert problem is derived rigorously only for $p=3$ but the general framework of the proof can be extended to the Cauchy chain of arbitrary length $p$.


title: From Rigid Supersymmetry to Twisted Holomorphic Theories

abstract: We study N=1 field theories with a U(1)_R symmetry on compact four-manifolds M. Supersymmetry requires M to be a complex manifold. The supersymmetric theory on M can be described in terms of conventional fields coupled to background supergravity, or in terms of twisted fields adapted to the complex geometry of M. Many properties of the theory that are difficult to see in one formulation are simpler in the other one. We use the twisted description to study the dependence of the partition function Z_M on the geometry of M, as well as coupling constants and background gauge fields, recovering and extending previous results. We also indicate how to generalize our analysis to three-dimensional N=2 theories with a U(1)_R symmetry. In this case supersymmetry requires M to carry a transversely holomorphic foliation, which endows it with a near-perfect analogue of complex geometry. Finally, we present new explicit formulas for the dependence of Z_M on the choice of U(1)_R symmetry in four and three dimensions, and illustrate them for complex manifolds diffeomorphic to S^3 x S^1, as well as general squashed three-spheres.


title: Analyzing the concept of super-efficiency in data envelopment analysis: A directional distance function approach

abstract: Based on the framework of the directional distance function, we conduct a systematic analysis on the measurement of super-efficiency in order to achieve two main objectives. Our primary purpose is developing two generalized directional measures of super-efficiency that completely resolve the crucial infeasibility issue, commonly arisen in the traditional super-efficiency measures. The secondary goal is to demonstrate that our directional super-efficiency models encompass the conventional ones as special cases. The proposed measures are advantageous because they circumvent biases in super-efficiency estimation due to input and output slacks. They are general, and satisfy several desirable properties, such as always feasibility, monotonicity, unit independence, and translation invariance.


title: Dark Sky Simulations: Early Data Release

abstract: The Dark Sky Simulations are an ongoing series of cosmological N-body simulations designed to provide a quantitative and accessible model of the evolution of the large-scale Universe. Such models are essential for many aspects of the study of dark matter and dark energy, since we lack a sufficiently accurate analytic model of non-linear gravitational clustering. In July 2014, we made available to the general community our early data release, consisting of over 55 Terabytes of simulation data products, including our largest simulation to date, which used $1.07 \times 10^{12}~(10240^3)$ particles in a volume $8h^{-1}\mathrm{Gpc}$ across. Our simulations were performed with 2HOT, a purely tree-based adaptive N-body method, running on 200,000 processors of the Titan supercomputer, with data analysis enabled by yt. We provide an overview of the derived halo catalogs, mass function, power spectra and light cone data. We show self-consistency in the mass function and mass power spectrum at the 1% level over a range of more than 1000 in particle mass. We also present a novel method to distribute and access very large datasets, based on an abstraction of the World Wide Web (WWW) as a file system, remote memory-mapped file access semantics, and a space-filling curve index. This method has been implemented for our data release, and provides a means to not only query stored results such as halo catalogs, but also to design and deploy new analysis techniques on large distributed datasets.


title: BRST invariant effective action of shadow fields, conformal fields, and AdS/CFT

abstract: Totally symmetric arbitrary spin massless and massive fields in AdS space are studied. For such fields, we obtain Lagrangians which are invariant under global BRST transformations. The Lagrangians are used for computation of partition functions and effective actions. We demonstrate that BRST invariant bulk action for massless field evaluated on the solution of Dirichlet problem for gauge massless fields and Faddeev-Popov fields leads to BRST invariant effective action for canonical shadow gauge fields and shadow Faddeev-Popov fields, while the BRST invariant bulk action for massive field evaluated on the solution of Dirichlet problem for gauge massive fields and Faddeev-Popov fields leads to BRST invariant effective action for anomalous shadow gauge fields and shadow Faddeev-Popov fields. The leading logarithmic divergence of the regularized effective action for the canonical shadow field leads to simple BRST invariant action of conformal field. We demonstrate that the Nakanishi-Laudrup fields entering the BRST invariant Lagrangian of conformal field can geometrically be interpreted as boundary values of massless AdS fields.


title: Compressed sensing for longitudinal MRI: An adaptive-weighted approach

abstract: Purpose: Repeated brain MRI scans are performed in many clinical scenarios, such as follow up of patients with tumors and therapy response assessment. In this paper, the authors show an approach to utilize former scans of the patient for the acceleration of repeated MRI scans. Methods: The proposed approach utilizes the possible similarity of the repeated scans in longitudinal MRI studies. Since similarity is not guaranteed, sampling and reconstruction are adjusted during acquisition to match the actual similarity between the scans. The baseline MR scan is utilized both in the sampling stage, via adaptive sampling, and in the reconstruction stage, with weighted reconstruction. In adaptive sampling, k-space sampling locations are optimized during acquisition. Weighted reconstruction uses the locations of the nonzero coefficients in the sparse domains as a prior in the recovery process. The approach was tested on 2D and 3D MRI scans of patients with brain tumors. Results: The longitudinal adaptive CS MRI (LACS-MRI) scheme provides reconstruction quality which outperforms other CS-based approaches for rapid MRI. Examples are shown on patients with brain tumors and demonstrate improved spatial resolution. Compared with data sampled at Nyquist rate, LACS-MRI exhibits Signal-to-Error Ratio (SER) of 24.8dB with undersampling factor of 16.6 in 3D MRI. Conclusions: The authors have presented a novel method for image reconstruction utilizing similarity of scans in longitudinal MRI studies, where possible. The proposed approach can play a major part and significantly reduce scanning time in many applications that consist of disease follow-up and monitoring of longitudinal changes in brain MRI.


title: To Principles of Quantum Mechanics Development

abstract: New insight to the principles of the quantum physics development is given. The correct ways for the construction of new versions of quantum mechanics on the second main postulate base are discussed. The conclusion on the status of the second main postulate is given. Its formulation in all textbooks has to be represented in the form of statement, since the hypothesis of Schr\"odinger on the existance of the field scalar function, being to be observable quantity, just charge density, is strictly proved for the case of EM-field, the role of which is argued to be decisive for the dynamics of the atomic systems. It is shown, that the field scalar function, being to be the function the only of coordinates and time, actually describes the state of the system. The second main postulate in Schr\"odinger formulation is mathematically strictly grounded, but in the popular probabilistic form used in modern textbooks on quantum theory it cannot be proved. The probabilistic theatise, proposed by Born is true in a number of special cases, quite correctly indicated by Dirac. The possible ways of the development of quantum theory, based on clear understanding of the origin of corpuscular-wave dualism are analysed.


title: Cosmological perturbations in warm-tachyon inflationary universe model with viscous pressure

abstract: We study the warm-tachyon inflationary universe model with viscous pressure in high-dissipation regime. General conditions which are required for this model to be realizable are derived in the slow-roll approximation. We present analytic expressions for density perturbation and amplitude of tensor perturbation in longitudinal gauge. Expressions of tensor-to-scalar ratio, scalar spectral index and its running are obtained. We develop our model by using exponential potential, the characteristics of this model are calculated for two specific cases in great details: 1- Dissipative parameter $\Gamma$ and bulk viscous parameter $\zeta$ are constant parameters. 2- Dissipative parameter is a function of tachyon field $\phi$ and bulk viscous parameter is a function of matter-radiation mixture energy density $\rho$. The parameters of the model are restricted by recent observational data from the nine-year Wilkinson microwave anisotropy probe (WMAP9), Planck and BICEP2 data.


title: Is a Single Photon Always Circularly Polarized? A Proposed Experiment using a Superconducting Microcalorimeter Photon Detector

abstract: A single photon is well known to have spin S = hbar, which would correspond to circular polarization, and all quantum transitions with photon absorption or emission correspond to DeltaS = +/-hbar. However, it is also widely believed that a single photon may be linearly polarized, which would correspond to a state with S = 0. Indeed, linearly polarized single photons are central to most quantum entanglement experiments. On the contrary, it has recently been suggested (based on a realistic spin-quantized wave picture of quantum states) that a linearly polarized photon state must be a superposition of a pair of circularly polarized photons, each with S = +/-hbar. This question cannot be resolved using a conventional photon detector, which generally cannot distinguish one photon from two simultaneous photons. However, it can be addressed using a superconducting microcalorimeter detector with sub-eV energy resolution and high quantum efficiency (QE). A careful experiment demonstrating this photon pairing could place in question some of the paradoxical central foundations of modern quantum theory, including quantum entanglement and nonlocality.


title: Strategies for Probing Non-Minimal Dark Sectors at Colliders: The Interplay Between Cuts and Kinematic Distributions

abstract: In this paper, we examine the strategies and prospects for distinguishing between traditional dark-matter models and models with non-minimal dark sectors --- including models of Dynamical Dark Matter (DDM) --- at hadron colliders. For concreteness, we focus on events with two hadronic jets and large missing transverse energy at the Large Hadron Collider (LHC). As we discuss, simple "bump-hunting" searches are not sufficient; probing non-minimal dark sectors typically requires an analysis of the actual shapes of the distributions of relevant kinematic variables. We therefore begin by identifying those kinematic variables whose distributions are particularly suited to this task. However, as we demonstrate, this then leads to a number of additional subtleties, since cuts imposed on the data for the purpose of background reduction can at the same time have the unintended consequence of distorting these distributions in unexpected ways, thereby obscuring signals of new physics. We therefore proceed to study the correlations between several of the most popular relevant kinematic variables currently on the market, and investigate how imposing cuts on one or more of these variables can impact the distributions of others. Finally, we combine our results in order to assess the prospects for distinguishing non-minimal dark sectors in this channel at the upgraded LHC.


title: Unbroken $SU(2)$ at a 100 TeV collider

abstract: A future 100 TeV pp collider will explore energies much higher than the scale of electroweak (EW) symmetry breaking. In this paper we study some of the phenomenological consequences of this fact, concentrating on enhanced bremsstrahlung of EW gauge bosons. We survey a handful of possible new physics experimental searches one can pursue at a 100 TeV collider using this phenomenon. The most dramatic effect is the non-negligible radiation of EW gauge bosons from neutrinos, making them partly visible objects. The presence of collinear EW radiation allows for the full reconstruction of neutrinos under certain circumstances. We also show that the presence of EW radiation allows one to distinguish the $SU(2)$ quantum numbers of various new physics particles. We consider examples of two completely different new physics paradigms, additional gauge groups and SUSY, where the bremsstrahlung radiation of $W$ and $Z$ from $W'$s, $Z'$s or stops allows one to determine the couplings and the mixing angles of the new particles (respectively). Finally, we show how the emission of $W$s and $Z$s from high $p_T$ Higgs bosons can be used to test the couplings of new physics to the Higgs boson.


title: Determining stellar atmospheric parameters and chemical abundances of FGK stars with iSpec

abstract: Context. An increasing number of high-resolution stellar spectra is available today thanks to many past and ongoing extensive spectroscopic surveys. Consequently, the scientific community needs automatic procedures to derive atmospheric parameters and individual element abundances. Aims. Based on the widely known SPECTRUM code by R. O. Gray, we developed an integrated spectroscopic software framework suitable for the determination of atmospheric parameters (i.e., effective temperature, surface gravity, metallicity) and individual chemical abundances. The code, named iSpec and freely distributed, is written mainly in Python and can be used on different platforms. Methods. iSpec can derive atmospheric parameters by using the synthetic spectral fitting technique and the equivalent width method. We validated the performance of both approaches by developing two different pipelines and analyzing the Gaia FGK benchmark stars spectral library. The analysis was complemented with several tests designed to assess other aspects, such as the interpolation of model atmospheres and the performance with lower quality spectra. Results. We provide a code ready to perform automatic stellar spectral analysis. We successfully assessed the results obtained for FGK stars with high-resolution and high signal-to-noise spectra.


title: Discovery of an Overdensity of Lyman-alpha Emitters Around a $\mathrm{z}\sim4$ QSO with the Large Binocular Telescope

abstract: Measurements of QSO clustering in the SDSS show that $\mathrm{z}>4$ QSOs are some of the most highly biased objects in the universe. Their large correlation lengths of $r_0 \sim 20h^{-1}$Mpc are comparable to the most massive clusters of galaxies in the universe today and suggest that these QSOs may mark the locations of massive cluster progenitors at high redshift. We report the discovery of an overdensity of LBGs around QSO SDSSJ114514.18+394715.9 as part of our survey to identify Lyman-Break galaxies (LBGs) around luminous $\mathrm{z}\sim4$ QSOs. In this field three of the eight LBGs with secure redshifts are consistent with the redshift of the QSO. We find that the likelihood that this is merely an apparent overdensity due to the chance selection of field galaxies is only 0.02%, based on comparisons to simulations and our modeled selection efficiency. Overall, our survey finds four of the 15 LBGs with secure redshifts are consistent with the redshifts of their respective QSOs, which is consistent with luminous QSOs residing in larger haloes.


title: The Most Distant Stars in the Milky Way

abstract: We report on the discovery of the most distant Milky Way (MW) stars known to date: ULAS J001535.72$+$015549.6 and ULAS J074417.48$+$253233.0. These stars were selected as M giant candidates based on their infrared and optical colors and lack of proper motions. We spectroscopically confirmed them as outer halo giants using the MMT/Red Channel spectrograph. Both stars have large estimated distances, with ULAS J001535.72$+$015549.6 at $274 \pm 74$ kpc and ULAS J074417.48$+$253233.0 at 238 $\pm$ 64 kpc, making them the first MW stars discovered beyond 200 kpc. ULAS J001535.72$+$015549.6 and ULAS J074417.48$+$253233.0 are both moving away from the Galactic center at $52 \pm 10$ km s$^{-1}$ and $24 \pm 10$ km s$^{-1}$, respectively. Using their distances and kinematics, we considered possible origins such as: tidal stripping from a dwarf galaxy, ejection from the MW's disk, or membership in an undetected dwarf galaxy. These M giants, along with two inner halo giants that were also confirmed during this campaign, are the first to map largely unexplored regions of our Galaxy's outer halo.


title: Transcendence and CM on Borcea-Voisin towers of Calabi-Yau manifolds

abstract: This paper is a sequel to a paper by the author and Marvin D. Tretkoff (reference [49]), in which we showed the validity of a special case of a conjecture of Green, Griffths and Kerr for certain families of Calabi-Yau manifolds over Hermitian symmetric domains. These results are analogues of a celebrated theorem of Th. Schneider on the transcendence of values of the elliptic modular function, and its generalization to the context of abelian varieties by the author, Shiga and Wolfart. In the present paper, we apply related techniques to many of the examples of families of Calabi-Yau varieties with dense sets of CM fibers in the work of Rohde, and in particular to Borcea-Voisin towers. Our results fit into a broader context of transcendence theory for variations of Hodge structure of higher weight.


title: The Vast Polar Structure of the Milky Way Attains New Members

abstract: The satellite galaxies of the Milky Way (MW) align with and preferentially orbit in a vast polar structure (VPOS), which also contains globular clusters, stellar and gaseous streams. Similar alignments have been discovered around several other host galaxies. We test whether recently discovered objects in the MW halo, the satellite galaxy/globular cluster transition object named PSO J174.0675-10.8774 or Crater and three stellar streams, are part of the VPOS. Crater is situated close to the VPOS. Incorporating the new object in the VPOS-plane fit slightly improves the alignment of the plane with other features such as the Magellanic stream and the average orbital plane of the satellites co-orbiting in the VPOS. We predict Crater's proper motion by assuming that it, too, orbits in the VPOS. One of the three streams aligns well with the VPOS. Surprisingly, it appears to lie in the exact same orbital plane as the Palomar 5 stream and shares its distance, suggesting a direct connection between the two. The stream also crosses close to the Fornax dwarf galaxy and is oriented approximately along the galaxy's direction of motion. The two other streams cannot align closely with the VPOS because they were discovered in the direction of M31/M33, which is outside of the satellite structure. The VPOS thus attains two new members. This further emphasizes that the highly anisotropic and correlated distribution of satellite objects requires an explanation beyond the suggestion that the MW satellite system is an extreme statistical outlier of a $\Lambda$CDM sub-halo system.


title: On the parallel spectrum in MHD Turbulence

abstract: Anisotropy of MHD turbulence has been studied extensively for many years, most prominently by measurements in the solar wind and high resolution simulations. The spectrum parallel to the local magnetic field was observed to be steeper than perpendicular spectrum, typically $k^{-2}$, consistent with the widely accepted Goldreich-Sridhar (1995) model. In this Letter I looked deeper into the nature of the relation between parallel and perpendicular spectra and argue that this $k^{-2}$ scaling has the same origin as the $\omega^{-2}$ scaling of Lagrangian frequency spectrum in strong hydrodynamic turbulence. This follows from the fact that Alfv\'en waves propagate along magnetic field lines. It is now became clear that the observed anisotropy can be argued without invocation of the "critical balance" argument and is more robust that was previously thought. The relation between parallel (Lagrangian) and perpendicular (Eulerian) spectra is inevitable consequence of strong turbulence of Alfven waves, rather than a conjecture based on the uncertainty relation. I tested this using high-resolution simulations of MHD turbulence, in particular I verified that the cutoff of the parallel spectrum scales as Kolmogorov timescale, not lengthscale.


title: The Gentle Growth of Galaxies at High Redshifts in Overdense Environments

abstract: We have explored prevailing modes of galaxy growth for redshifts z ~ 6-14, comparing substantially overdense and normal regions of the universe, using high-resolution zoom-in cosmological simulations. Such rare overdense regions have been projected to host high-z quasars. We demonstrate that galaxies in such environments grow predominantly by a smooth accretion from cosmological filaments which dominates the mass input from major, intermediate and minor mergers. We find that by z ~6, the accumulated galaxy mass fraction from mergers falls short by a factor of 10 of the cumulative accretion mass for galaxies in the overdense regions, and by a factor of 5 in the normal environments. Moreover, the rate of the stellar mass input from mergers also lies below that of an in-situ star formation (SF) rate. The fraction of stellar masses in galaxies contributed by mergers in overdense regions is ~12%, and ~33% in the normal regions, at these redshifts. Our median SF rates for ~few X 10^9 Mo galaxies agrees well with the recently estimated rates for z ~ 7 galaxies from Spitzer's SURF-UP survey. Finally, we find that the main difference between the normal and overdense regions lies in the amplified growth of massive galaxies in massive dark matter halos. This leads to the formation of >= 10^{10} Mo galaxies due to the ~100-fold increase in mass during the above time period. Such galaxies are basically absent in the normal regions at these redshifts.


title: Non-Linear Dynamics of Accretion Disks with Stochastic Viscosity

abstract: We present a non-linear numerical model for a geometrically thin accretion disk with the addition of stochastic non-linear fluctuations in the viscous parameter. These numerical realizations attempt to study the stochastic effects on the disk angular momentum transport. We show that this simple model is capable of reproducing several observed phenomenologies of accretion driven systems. The most notable of these is the observed linear rms-flux relationship in the disk luminosity. This feature is not formally captured by the linearized disk equations used in previous work. A Fourier analysis of the dissipation and mass accretion rates across disk radii show coherence for frequencies below the local viscous frequency. This is consistent with the coherence behavior observed in astrophysical sources such as Cygnus X-1.


title: Ward Identities for Hall Transport

abstract: We derive quantum field theory Ward identities based on linear area preserving and conformal transformations in 2+1 dimensions. The identities relate Hall viscosities, Hall conductivities and the angular momentum. They apply both for relativistic and non relativistic systems, at zero and at finite temperature. We consider systems with or without translation invariance, and introduce an external magnetic field and viscous drag terms. A special case of the identities yields the well known relation between the Hall conductivity and half the angular momentum density.


title: A near-infrared SETI experiment: instrument overview

abstract: We are designing and constructing a new SETI (Search for Extraterrestrial Intelligence) instrument to search for direct evidence of interstellar communications via pulsed laser signals at near-infrared wavelengths. The new instrument design builds upon our past optical SETI experiences, and is the first step toward a new, more versatile and sophisticated generation of very fast optical and near-infrared pulse search devices. We present our instrumental design by giving an overview of the opto-mechanical design, detector selection and characterization, signal processing, and integration procedure. This project makes use of near-infrared (950-1650 nm) discrete amplification Avalanche Photodiodes (APD) that have greater than 1 GHz bandwidths with low noise characteristics and moderate gain (~10^4). We have investigated the use of single versus multiple detectors in our instrument (see Maire et al., this conference), and have optimized the system to have both high sensitivity and low false coincidence rates. Our design is optimized for use behind a 1m telescope and includes an optical camera for acquisition and guiding. A goal is to make our instrument relatively economical and easy to duplicate. We describe our observational setup and our initial search strategies for SETI targets, and for potential interesting compact astrophysical objects.


title: A near-infrared SETI experiment: probability distribution of false coincidences

abstract: A Search for Extraterrestrial Life (SETI), based on the possibility of interstellar communication via laser signals, is being designed to extend the search into the near-infrared spectral region (Wright et al, this conference). The dedicated near-infrared (900 to 1700 nm) instrument takes advantage of a new generation of avalanche photodiodes (APD), based on internal discrete amplification. These discrete APD (DAPD) detectors have a high speed response ($>$ 1 GHz) and gain comparable to photomultiplier tubes, while also achieving significantly lower noise than previous APDs. We are investigating the use of DAPD detectors in this new astronomical instrument for a SETI search and transient source observations. We investigated experimentally the advantages of using a multiple detector device operating in parallel to remove spurious signals. We present the detector characterization and performance of the instrument in terms of false positive detection rates both theoretically and empirically through lab measurements. We discuss the required criteria that will be needed for laser light pulse detection in our experiment. These criteria are defined to optimize the trade between high detection efficiency and low false positive coincident signals, which can be produced by detector dark noise, background light, cosmic rays, and astronomical sources. We investigate experimentally how false coincidence rates depend on the number of detectors in parallel, and on the signal pulse height and width. We also look into the corresponding threshold to each of the signals to optimize the sensitivity while also reducing the false coincidence rates. Lastly, we discuss the analytical solution used to predict the probability of laser pulse detection with multiple detectors.


title: Formation of irregular and runaway moons/exomoons through moon-moon scattering

abstract: Gas giant planets in the Solar system host large satellite systems with multiple regular and irregular moons. Regular moons revolve around their host planet in circular, low inclination short period orbits, and are thought to form in-situ through coagulation processes. In contrast, irregular moons have highly inclined (and even retrograde), typically eccentric and long period orbits around their host planet. Irregular moons are therefore often thought to be unbound objects in helio-centric orbits that were later captured to their current orbits around the planet. Here we study the a different possibility in which regular moons form in-situ, and dynamically evolve through mutual moon-moon scattering. We find that such evolution can excite the satellites into high eccentricities and inclinations. We find that moons are either ejected from the host planet to become runaway moons, or stay bound and become prograde orbiting irregular moons with inclined and eccentric orbits around their host planet. Ejected moons, unbound to the planet, can later be temporarily re-captured by the host planet even at retrograde orbits. Such moons are eventually re-ejected from the system or collide with the planet, at least in the absence of dissipative processes (e.g. collisions with existing bound moons, a debris disk or through tidal interactions with the host planet), not currently modeled. Uncaptured runaway moons may eventually be ejected from the Solar system, or be captured into stable helio-centric orbits and contribute to the populations of asteroidal or trans-Neptunian objects. Such scenarios are potentially relevant both for the gas-giant satellites in the Solar system and for the dynamical evolution of exomoons.


title: Student Use of a Single Lecture Video in a Flipped Introductory Mechanics Course

abstract: In the Fall of 2013, Georgia Tech offered a 'flipped' calculus-based introductory mechanics class as an alternative to the traditional large-enrollment lecture class. This class flipped instruction by introducing new material outside of the classroom through pre-recorded, lecture videos. Video lectures constituted students' initial introduction to course material. We analyze how students engaged with online lecture videos via 'clickstream' data, consisting of time-stamped interactions (plays, pauses, seeks, etc.) with the online video player. Analysis of these events has shown that students may be focusing on elements of the video that facilitate a 'correct' solution.


title: B-modes and the Nature of Inflation

abstract: Observations of the cosmic microwave background do not yet determine whether inflation was driven by a slowly-rolling scalar field or involved another physical mechanism. In this paper we discuss the prospects of using the power spectra of scalar and tensor modes to probe the nature of inflation. We focus on the leading modification to the slow-roll dynamics, which entails a sound speed $c_s$ for the scalar fluctuations. We derive analytically a lower bound on $c_s$ in terms of a given tensor-to-scalar ratio $r$, taking into account the difference in the freeze-out times between the scalar and tensor modes. We find that any detection of primordial B-modes with $r > 0.01$ implies a lower bound on $c_s$ that is stronger than the bound derived from the absence of non-Gaussianity in the Planck data. For $r \gtrsim 0.1$, the bound would be tantalizingly close to a critical value for the sound speed, $(c_s)_\star = 0.47$ (corresponding to $(f_{\rm NL}^{\rm equil})_\star = -0.93$), which we show serves as a threshold for non-trivial dynamics beyond slow-roll. We also discuss how an order-one level of equilateral non-Gaussianity is a natural observational target for other extensions of the canonical paradigm.


title: Abiotic Ozone and Oxygen in Atmospheres Similar to Prebiotic Earth

abstract: The search for life on planets outside our solar system will use spectroscopic identification of atmospheric biosignatures. The most robust remotely-detectable potential biosignature is considered to be the detection of oxygen (O_2) or ozone (O_3) simultaneous to methane (CH_4) at levels indicating fluxes from the planetary surface in excess of those that could be produced abiotically. Here, we use an altitude-dependent photochemical model with the enhanced lower boundary conditions necessary to carefully explore abiotic O_2 and O_3 production on lifeless planets with a wide variety of volcanic gas fluxes and stellar energy distributions. On some of these worlds, we predict limited O_2 and O_3 build up, caused by fast chemical production of these gases. This results in detectable abiotic O_3 and CH_4 features in the UV-visible, but no detectable abiotic O_2 features. Thus, simultaneous detection of O_3 and CH_4 by a UV-visible mission is not a strong biosignature without proper contextual information. Discrimination between biological and abiotic sources of O_2 and O_3 is possible through analysis of the stellar and atmospheric context - particularly redox state and O atom inventory - of the planet in question. Specifically, understanding the spectral characteristics of the star and obtaining a broad wavelength range for planetary spectra should allow more robust identification of false positives for life. This highlights the importance of wide spectral coverage for future exoplanet characterization missions. Specifically, discrimination between true- and false-positives may require spectral observations that extend into infrared wavelengths, and provide contextual information on the planet's atmospheric chemistry.


title: Scattering of Dark Particles with Light Mediators

abstract: We present a treatment of the high energy scattering of dark Dirac fermions from nuclei, mediated by the exchange of a light vector boson. The dark fermions are produced by proton-nucleus interactions in a fixed target and, after traversing shielding that screens out strongly interacting products, appear similarly to neutrino neutral current scattering in a detector. Using the Fermilab experiment E613 as an example, we place limits on a secluded dark matter scenario. Visible scattering in the detector includes both the familiar regime of large momentum transfer to the nucleus ($Q^2$) described by deeply inelastic scattering, as well as small $Q^2$ kinematics described by the exchanged vector mediator fluctuating into a quark-antiquark pair whose interaction with the nucleus is described by a saturation model. We find that the improved description of the low $Q^2$ scattering leads to important corrections, resulting in more robust constraints in a regime where a description entirely in terms of deeply inelastic scattering cannot be trusted.


title: Valley Zeeman Effect in Elementary Optical Excitations of a Monolayer WSe2

abstract: A monolayer of a transition metal dichalcogenide (TMD) such as WSe$_2$ is a two-dimensional (2D) direct band-gap valley-semiconductor having an effective Honeycomb lattice structure with broken inversion symmetry. The inequivalent valleys in the Brillouin zone could be selectively addressed using circularly-polarized light fields, suggesting the possibility for magneto-optical measurement and manipulation of the valley pseudospin degree of freedom. Here we report such experiments that demonstrate the valley Zeeman effect -- strongly anisotropic lifting of the degeneracy of the valley pseudospin degree of freedom using an external magnetic field. While the valley-splitting measured using the exciton transition is consistent with the difference of the conduction and valence band orbital magnetic moments, the trion transition exhibits an unexpectedly large valley Zeeman effect which cannot be understood using an independent electron-hole picture. Instead, we find an explanation using the recently predicted large Berry curvature and the associated magnetic moment for the electron-hole exchange interaction modified trion dispersion. Our results raise the possibility of observing optical excitation induced valley Hall effect in monolayer TMDs or topological states of photons strongly coupled to trion excitations in a microcavity.


title: An XMM-Newton View of the Radio Galaxy 3C 411

abstract: We present the first high signal-to-noise XMM-Newton observations of the broad-line radio galaxy 3C 411. After fitting various spectral models, an absorbed double power-law continuum and a blurred relativistic disk reflection model (kdblur) are found to be equally plausible descriptions of the data. While the softer power-law component ($\Gamma$=2.11) of the double power-law model is entirely consistent with that found in Seyfert galaxies (and hence likely originates from a disk corona), the additional power law component is very hard ($\Gamma$=1.05); amongst the AGN zoo, only flat-spectrum radio quasars have such hard spectra. Together with the very flat radio-spectrum displayed by this source, we suggest that it should instead be classified as a FSRQ. This leads to potential discrepancies regarding the jet inclination angle, with the radio morphology suggesting a large jet inclination but the FSRQ classification suggesting small inclinations. The kdblur model predicts an inner disk radius of at most 20 r$_g$ and relativistic reflection.


title: The Complexity of Primes in Computable UFDs

abstract: In many simple integral domains, such as $\mathbb{Z}$ or $\mathbb{Z}[i]$, there is a straightforward procedure to determine if an element is prime by simply reducing to a direct check of finitely many potential divisors. Despite the fact that such a naive approach does not immediately translate to integral domains like $\mathbb{Z}[x]$ or the ring of integers in an algebraic number field, there still exist computational procedures that work to determine the prime elements in these cases. In contrast, we will show how to computably extend $\mathbb{Z}$ in such a way that we can control the ordinary integer primes in any $\Pi_2^0$ way, all while maintaining unique factorization. As a corollary, we establish the existence of a computable UFD such that the set of primes is $\Pi_2^0$-complete in every computable presentation.


title: Super-intuition and correlations with the future in Quantum Consciousness

abstract: The role of quantum information is discussed in the framework of Quantum Concept of Consciousness (QCC), based on the many-worlds interpretation of quantum mechanics (Everett interpretation). Within QCC the phenomenon of super-intuition is analyzed, which explains in particular the great scientific insights and realizes a sort of a "mental time machine". The recently expressed opinion that emergence of super-intuition requires transferring quantum information, which is banned by the impossibility of "quantum cloning", is critically considered. It is shown that in fact under QCC the emergence of the phenomenon of super-intuition requires not transferring quantum information, but only correlating various fragments of it. Actual examples of scientific insights that have been documented by prominent scientists, suggest that the marker of the correlation in this case is a strong positive emotion, which can be called cognitive euphoria.


title: On the Spectral Efficiency of Full-Duplex Small Cell Wireless Systems

abstract: We investigate the spectral efficiency of full-duplex small cell wireless systems, in which a full-duplex capable base station (BS) is designed to send/receive data to/from multiple halfduplex users on the same system resources. The major hurdle for designing such systems is due to the self-interference at the BS and co-channel interference among users. Hence, we consider a joint beamformer design to maximize the spectral efficiency subject to certain power constraints. The design problem is first formulated as a rank-constrained optimization one, and the rank relaxation method is then applied. However the relaxed problem is still nonconvex, and thus optimal solutions are hard to find. Herein, we propose two provably convergent algorithms to obtain suboptimal solutions. Based on the concept of the difference of convex functions programming, we approximate the design problem by a determinant maximization program in each iteration of the first algorithm. The second method is built upon the sequential parametric convex approximation method, which allows us to transform the relaxed problem into a semidefinite program in each iteration. Extensive numerical experiments under small cell setups illustrate that the full-duplex system with the proposed algorithms can achieve a large gain over the half-duplex one.


title: Torification of diagonalizable group actions on toroidal schemes

abstract: We study actions of diagonalizable groups on toroidal schemes (i.e. logarithmically regular logarithmic schemes). In particular, we show that for so-called toroidal actions the quotient is again a toroidal scheme. Our main result constructs for an arbitrary action a canonical torification - making the action toridal after an equivariant blowings up. This extends earlier results of Abramovich-de Jong, Abramovich-Karu-Matsuki-W{\l}odarczyk, and Gabber in various aspects.


title: A Statistical Modeling Approach to Computer-Aided Quantification of Dental Biofilm

abstract: Biofilm is a formation of microbial material on tooth substrata. Several methods to quantify dental biofilm coverage have recently been reported in the literature, but at best they provide a semi-automated approach to quantification with significant input from a human grader that comes with the graders bias of what are foreground, background, biofilm, and tooth. Additionally, human assessment indices limit the resolution of the quantification scale; most commercial scales use five levels of quantification for biofilm coverage (0%, 25%, 50%, 75%, and 100%). On the other hand, current state-of-the-art techniques in automatic plaque quantification fail to make their way into practical applications owing to their inability to incorporate human input to handle misclassifications. This paper proposes a new interactive method for biofilm quantification in Quantitative light-induced fluorescence (QLF) images of canine teeth that is independent of the perceptual bias of the grader. The method partitions a QLF image into segments of uniform texture and intensity called superpixels; every superpixel is statistically modeled as a realization of a single 2D Gaussian Markov random field (GMRF) whose parameters are estimated; the superpixel is then assigned to one of three classes (background, biofilm, tooth substratum) based on the training set of data. The quantification results show a high degree of consistency and precision. At the same time, the proposed method gives pathologists full control to post-process the automatic quantification by flipping misclassified superpixels to a different state (background, tooth, biofilm) with a single click, providing greater usability than simply marking the boundaries of biofilm and tooth as done by current state-of-the-art methods.


title: Exponential self-similar mixing and loss of regularity for continuity equations

abstract: We consider the mixing behaviour of the solutions of the continuity equation associated with a divergence-free velocity field. In this announcement we sketch two explicit examples of exponential decay of the mixing scale of the solution, in case of Sobolev velocity fields, thus showing the optimality of known lower bounds. We also describe how to use such examples to construct solutions to the continuity equation with Sobolev but non-Lipschitz velocity field exhibiting instantaneous loss of any fractional Sobolev regularity.


title: Lessons from the Milky Way: the Kapteyn Universe

abstract: Jacobus Cornelius Kapteyn (1851-1922) presented a model for the distribution of stars in space together with a dynamical interpretation in terms of an equilibrium between the gravitational field of the stars and their random motion and rotation. In the vertical direction Kapteyn's results are substantially correct. Usually the Kapteyn Universe is described as being flawed due to neglect of interstellar absorption. Kapteyn was led to adopt this on the basis of widely accepted evidence by Shapley on an absence of reddening of stars in globular clusters. But another, equally important misconception was Kapteyn's interpretation of the two Star Streams as manifestations of two groups of stars rotating around a center in opposite directions. This was supported by the observation of very different mixes in stellar types in the two streams. Had Kapteyn adopted the absorption as he himself had determined it he would not have been able to arrive at a consistent picture.


title: Effective field theory for the bulk-edge correspondence in a two-dimensional Z_2 topological insulator with Rashba interactions

abstract: We determine the effective field theory in (2+1)-dimensional space and time that captures the long wave length and low-energy limit of fermions hopping on a honeycomb lattice at half-filling when both a dominant intrinsic and subdominant Rashba spin-orbit couplings are present. This effective field theory for a Z_2 topological insulator (the Kane-Mele model at vanishing uniform and staggered chemical potentials) is a perturbation around a double Chern-Simons theory, with the U(1) gauge invariance associated to spin conservation explicitly broken due to the Rashba spin orbit coupling. Nonetheless, we find that the effective field theory has a BRST symmetry that allows us to construct the bulk-edge correspondence.


title: Vortex solutions of an Abelian Higgs model with visible and hidden sectors

abstract: We study vortex solutions in a theory with dynamics governed by two weakly coupled Abelian Higgs models, describing a hidden sector and a visible sector. We analyze the radial dependence of the axially symmetric solutions constructed numerically and discuss the stability of vortex configurations for different values of the model parameters, studying in detail vortex decay into lower energy configurations. We find that even in a weak coupling regime vortex solutions strongly depend on the parameters of both the visible and hidden sectors. We also discuss on qualitative grounds possible implications of the existence of a hidden sector in connection with superconductivity.


title: Nonparametric empirical Bayes and maximum likelihood estimation for high-dimensional data analysis

abstract: Nonparametric empirical Bayes methods provide a flexible and attractive approach to high-dimensional data analysis. One particularly elegant empirical Bayes methodology, involving the Kiefer-Wolfowitz nonparametric maximum likelihood estimator (NPMLE) for mixture models, has been known for decades. However, implementation and theoretical analysis of the Kiefer-Wolfowitz NPMLE are notoriously difficult. A fast algorithm was recently proposed that makes NPMLE-based procedures feasible for use in large-scale problems, but the algorithm calculates only an approximation to the NPMLE. In this paper we make two contributions. First, we provide upper bounds on the convergence rate of the approximate NPMLE's statistical error, which have the same order as the best known bounds for the true NPMLE. This suggests that the approximate NPMLE is just as effective as the true NPMLE for statistical applications. Second, we illustrate the promise of NPMLE procedures in a high-dimensional binary classification problem. We propose a new procedure and show that it vastly outperforms existing methods in experiments with simulated data. In real data analyses involving cancer survival and gene expression data, we show that it is very competitive with several recently proposed methods for regularized linear discriminant analysis, another popular approach to high-dimensional classification.


title: Parallel MATLAB Techniques

abstract: In this chapter, we show why parallel MATLAB is useful, provide a comparison of the different parallel MATLAB choices, and describe a number of applications in Signal and Image Processing: Audio Signal Processing, Synthetic Aperture Radar (SAR) Processing and Superconducting Quantum Interference Filters (SQIFs). Each of these applications have been parallelized using different methods (Task parallel and Data parallel techniques). The applications presented may be considered representative of type of problems faced by signal and image processing researchers. This chapter will also strive to serve as a guide to new signal and image processing parallel programmers, by suggesting a parallelization strategy that can be employed when developing a general parallel algorithm. The objective of this chapter is to help signal and image processing algorithm developers understand the advantages of using parallel MATLAB to tackle larger problems while staying within the powerful environment of MATLAB.


title: Nonlinear Bias of Cosmological Halo Formation in the Early Universe

abstract: We present estimates of the nonlinear bias of cosmological halo formation, spanning a wide range in the halo mass from $\sim 10^{5} M_\odot$ to $\sim 10^{12} M_\odot$, based upon both a suite of high-resolution cosmological N-body simulations and theoretical predictions. The halo bias is expressed in terms of the mean bias and stochasticity as a function of local overdensity ($\delta$), under different filtering scales, which is realized as the density of individual cells in uniform grids. The sampled overdensities span a range wide enough to provide the fully nonlinear bias effect on the formation of haloes. A strong correlation between $\delta$ and halo population overdensity $\delta_h$ is found, along with sizable stochasticity. We find that the empirical mean halo bias matches, with good accuracy, the prediction by the peak-background split method based on the excursion set formalism, as long as the empirical, globally-averaged halo mass function is used. Consequently, this bias formalism is insensitive to uncertainties caused by varying halo identification schemes, and can be applied generically. We also find that the probability distribution function of biased halo numbers has wider distribution than the pure Poisson shot noise, which is attributed to the sub-cell scale halo correlation. We explicitly calculate this correlation function and show that both overdense and underdense regions have positive correlation, leading to stochasticity larger than the Poisson shot noise in the range of haloes and halo-collapse epochs we study.


title: Hopf bifurcation from fronts in the Cahn-Hilliard equation

abstract: We study Hopf bifurcation from traveling-front solutions in the Cahn-Hilliard equation. The primary front is induced by a moving source term. Models of this form have been used to study a variety of physical phenomena, including pattern formation in chemical deposition and precipitation processes. Technically, we study bifurcation in the presence of essential spectrum. We contribute a simple and direct functional analytic method and determine bifurcation coefficients explicitly. Our approach uses exponential weights to recover Fredholm properties and spectral flow ideas to compute Fredholm indices. Simple mass conservation helps compensate for negative indices. We also construct an explicit, prototypical example, prove the existence of a bifurcating front, and determine the direction of bifurcation.


title: Critical points of the linear entropy for pure L-qubit states

abstract: We present a substancially improved version of the method proposed in Sawicki et al (2012, 2014) for finding critical points of the linear entropy for L-qubit system. The new approach is based on the corespondance between momentum maps for abelian and non-abelian groups, as described in Kirwan (1984). The proposed method can be implemented numerically much easier than the previous one.


title: Approximately Stable, School Optimal, and Student-Truthful Many-to-One Matchings (via Differential Privacy)

abstract: We present a mechanism for computing asymptotically stable school optimal matchings, while guaranteeing that it is an asymptotic dominant strategy for every student to report their true preferences to the mechanism. Our main tool in this endeavor is differential privacy: we give an algorithm that coordinates a stable matching using differentially private signals, which lead to our truthfulness guarantee. This is the first setting in which it is known how to achieve nontrivial truthfulness guarantees for students when computing school optimal matchings, assuming worst- case preferences (for schools and students) in large markets.


title: Private Pareto Optimal Exchange

abstract: We consider the problem of implementing an individually rational, asymptotically Pareto optimal allocation in a barter-exchange economy where agents are endowed with goods and have preferences over the goods of others, but may not use money as a medium of exchange. Because one of the most important instantiations of such economies is kidney exchange -- where the "input"to the problem consists of sensitive patient medical records -- we ask to what extent such exchanges can be carried out while providing formal privacy guarantees to the participants. We show that individually rational allocations cannot achieve any non-trivial approximation to Pareto optimality if carried out under the constraint of differential privacy -- or even the relaxation of \emph{joint} differential privacy, under which it is known that asymptotically optimal allocations can be computed in two-sided markets, where there is a distinction between buyers and sellers and we are concerned only with privacy of the buyers~\citep{Matching}. We therefore consider a further relaxation that we call \emph{marginal} differential privacy -- which promises, informally, that the privacy of every agent $i$ is protected from every other agent $j \neq i$ so long as $j$ does not collude or share allocation information with other agents. We show that, under marginal differential privacy, it is possible to compute an individually rational and asymptotically Pareto optimal allocation in such exchange economies.


title: A Bellman View of Jesse Livermore

abstract: Richard Bellman's Principle of Optimality, formulated in 1957, is the heart of dynamic programming, the mathematical discipline which studies the optimal solution of multi-period decision problems. In this paper, we look at the main trading principles of Jesse Livermore, the legendary stock operator whose method was published in 1923, from a Bellman point of view.


title: Thermal Properties of Asymmetric Nuclear Matter

abstract: The thermal properties of asymmetric nuclear matter are investigated in a relativistic mean- field approach. We start from free space NN-interactions and derive in-medium self-energies by Dirac-Brueckner theory. By the DDRH procedure we derive in a self-consistent approach density- dependent meson-baryon vertices. At the mean-field level, we include isoscalar and isovector scalar and vector interactions. The nuclear equation of state is investigated for a large range of total baryon densities up to the neutron star regime, the full range of asymmetries from symmetric nuclear matter to pure neutron matter, and temperatures up to T~100 MeV. The isovector-scalar self-energies are found to modify strongly the thermal properties of asymmetric nuclear matter. A striking result is the change of phase transitions when isovector-scalar self-energies are included.


title: Quadrature rules for $L^1$-weighted norms of orthogonal polynomials

abstract: In this paper we obtain $L^1$-weighted norms of classical orthogonal polynomials (Hermite, Laguerre and Jacobi polynomials) in terms of the zeros of these orthogonal polynomials; these expressions are usually known as quadrature rules. In particular these new formulae are useful to calculate directly some positive defined integrals as several examples show.


title: Magnetic Control of Valley Pseudospin in Monolayer WSe2

abstract: Local energy extrema of the bands in momentum space, or valleys, can endow electrons in solids with pseudo-spin in addition to real spin. In transition metal dichalcogenides this valley pseudo-spin, like real spin, is associated with a magnetic moment which underlies the valley-dependent circular dichroism that allows optical generation of valley polarization, intervalley quantum coherence, and the valley Hall effect. However, magnetic manipulation of valley pseudospin via this magnetic moment, analogous to what is possible with real spin, has not been shown before. Here we report observation of the valley Zeeman splitting and magnetic tuning of polarization and coherence of the excitonic valley pseudospin, by performing polarization-resolved magneto-photoluminescence on monolayer WSe2. Our measurements reveal both the atomic orbital and lattice contributions to the valley orbital magnetic moment; demonstrate the deviation of the band edges in the valleys from an exact massive Dirac fermion model; and reveal a striking difference between the magnetic responses of neutral and charged valley excitons which is explained by renormalization of the excitonic spectrum due to strong exchange interactions.


title: Learning Probabilistic Programs

abstract: We develop a technique for generalising from data in which models are samplers represented as program text. We establish encouraging empirical results that suggest that Markov chain Monte Carlo probabilistic programming inference techniques coupled with higher-order probabilistic programming languages are now sufficiently powerful to enable successful inference of this kind in nontrivial domains. We also introduce a new notion of probabilistic program compilation and show how the same machinery might be used in the future to compile probabilistic programs for efficient reusable predictive inference.


title: Interaction Driven Subgap Spin Exciton in the Kondo Insulator SmB6

abstract: Using inelastic neutron scattering, we map a 14 meV coherent resonant mode in the topological Kondo insulator SmB6 and describe its relation to the low energy insulating band structure. The resonant intensity is confined to the X and R high symmetry points, repeating outside the first Brillouin zone and dispersing less than 2 meV, with a 5d-like magnetic form factor. We present a slave-boson treatment of the Anderson Hamiltonian with a third neighbor dominated hybridized band structure. This approach produces a spin exciton below the charge gap with features that are consistent with the observed neutron scattering. We find that maxima in the wave vector dependence of the inelastic neutron scattering indicate band inversion.


title: Probing Satellite Quenching With Galaxy Clustering

abstract: Satellites within simulated massive clusters are significantly spatially correlated with each other, even when those satellites are not gravitationally bound to each other. This correlation is produced by satellites that entered their hosts relatively recently, and is undetectable for satellites that have resided in their hosts for multiple dynamical timescales. Therefore, a measurement of clustering statistics of cluster satellites may be used to determine the typical accretion redshifts of those satellites into their observed hosts. We argue that such measurements may be used to determine the fraction of satellite galaxies that were quenched by their current hosts, thereby discriminating among models for quenching of star formation in satellite galaxies.


title: Classifying Fonts and Calligraphy Styles Using Complex Wavelet Transform

abstract: Recognizing fonts has become an important task in document analysis, due to the increasing number of available digital documents in different fonts and emphases. A generic font-recognition system independent of language, script and content is desirable for processing various types of documents. At the same time, categorizing calligraphy styles in handwritten manuscripts is important for palaeographic analysis, but has not been studied sufficiently in the literature. We address the font-recognition problem as analysis and categorization of textures. We extract features using complex wavelet transform and use support vector machines for classification. Extensive experimental evaluations on different datasets in four languages and comparisons with state-of-the-art studies show that our proposed method achieves higher recognition accuracy while being computationally simpler. Furthermore, on a new dataset generated from Ottoman manuscripts, we show that the proposed method can also be used for categorizing Ottoman calligraphy with high accuracy.


title: Logic and linear algebra: an introduction

abstract: We give an introduction to logic tailored for algebraists, explaining how proofs in linear logic can be viewed as algorithms for constructing morphisms in symmetric closed monoidal categories with additional structure. This is made explicit by showing how to represent proofs in linear logic as linear maps between vector spaces. The interesting part of this vector space semantics is based on the cofree cocommutative coalgebra of Sweedler.


title: The Eccentric Collective BFKL Pomeron

abstract: We apply the flow analysis for multi-particle correlations used in heavy ion collisions to multi-particle production from a Pomeron. We show that the n'th order angular harmonic arising from an m particle correlation v_n[m] satisfies v_n[m] ~ v_n[p] for n $\ge$ 1. We discuss some implications of this for the Color Glass Condensate description of high energy hadronic collisions.


title: ALMA Observations of the T Tauri Binary System AS 205: Evidence for Molecular Winds and/or Binary Interactions

abstract: In this study, we present high-resolution millimeter observations of the dust and gas disk of the T Tauri star AS 205 N and its companion, AS 205 S, obtained with the Atacama Large Millimeter Array. The gas disk around AS 205 N, for which infrared emission spectroscopy demonstrates significant deviations from Keplerian motion that has been interpreted as evidence for a disk wind (Pontoppidan et al. 2011; Bast et al. 2011), also displays significant deviations from Keplerian disk emission in the observations presented here. Detections near both AS 205 N and S are obtained in 1.3 mm continuum, 12CO 2-1, 13CO 2-1 and C18O 2-1. The 12CO emission is extended up to 2 arcsec from AS 205N, and both 12CO and 13CO display deviations from Keplerian rotation at all angular scales. Two possible explanations for these observations hold up best to close scrutiny - tidal interaction with AS 205 S or disk winds (or a combination of the two), and we discuss these possibilities in some detail.


title: Radio-Continuum Study of the Nearby Sculptor Group Galaxies. Part 3: NGC 7793 at lambda=12.2, 6 and 3 cm

abstract: We re-examine a series of archived centimetre radio-continuum observations (lambda=16, 6 and 3 cm) focusing on NGC7793 using the Australia Telescope Compact Array. These new images are both very sensitive (rms=0.011 mJy/beam) and feature reasonably high angular resolution (down to < 3"). Using these images, a total of 76 discrete radio sources are identified, of which 57 have been classified. We also studied the radio component of the micro-quasar NGC7793-S26 which shows two distinct regions of somewhat steep spectral index between -0.3 and -0.7.


title: Genus-2 Jacobians with torsion points of large order

abstract: We produce new explicit examples of genus-2 curves over the rational numbers whose Jacobian varieties have rational torsion points of large order. In particular, we produce a family of genus-2 curves over Q whose Jacobians have a rational point of order 48, parametrized by a rank-2 elliptic curve over Q, and we exhibit a single genus-2 curve over Q whose Jacobian has a rational point of order 70, the largest order known. We also give new examples of genus-2 Jacobians with rational points of order 27, 28, and 39. Most of our examples are produced by `gluing' two elliptic curves together along their n-torsion subgroups, where n is either 2 or 3. The 2-gluing examples arise from techniques developed by the author in joint work with Lepr\'evost and Poonen 15 years ago. The 3-gluing examples are made possible by an algorithm for explicit 3-gluing over non-algebraically closed fields recently developed by the author in joint work with Br\"oker, Lauter, and Stevenhagen.


title: Systematic pathway to ${\mathcal{PT}}$ symmetry breaking in scattering systems

abstract: Recently [Phys. Rev. Lett. {\bf 106}, 093902 (2011)] it has been shown that $\mathcal{PT}$-symmetric scattering systems with balanced gain and loss, undergo a transition from $\mathcal{PT}$-symmetric scattering eigenstates, which are norm preserving, to symmetry broken pairs of eigenstates exhibiting net amplification and loss. In the present work we derive the existence of an invariant non-local current which can be directly associated with the observed transition playing the role of an "order parameter". The use of this current for the description of the $\mathcal{PT}$-symmetry breaking allows the extension of the known phase diagram to higher dimensions incorporating scattering states which are not eigenstates of the scattering matrix.


title: The Error Term in the Sato-Tate Conjecture

abstract: Let $f(z)=\sum_{n=1}^\infty a(n)e^{2\pi i nz}\in S_k^{new}(\Gamma_0(N))$ be a newform of even weight $k\geq2$ that does not have complex multiplication. Then $a(n)\in\mathbb{R}$ for all $n$, so for any prime $p$, there exists $\theta_p\in[0,\pi]$ such that $a(p)=2p^{(k-1)/2}\cos(\theta_p)$. Let $\pi(x)=\#\{p\leq x\}$. For a given subinterval $I\subset[0,\pi]$, the now-proven Sato-Tate Conjecture tells us that as $x\to\infty$, \[ \#\{p\leq x:\theta_p\in I\}\sim \mu_{ST}(I)\pi(x),\quad \mu_{ST}(I)=\int_{I} \frac{2}{\pi}\sin^2(\theta)~d\theta. \] Let $\epsilon>0$. Assuming that the symmetric power $L$-functions of $f$ are automorphic, we prove that as $x\to\infty$, \[ \#\{p\leq x:\theta_p\in I\}=\mu_{ST}(I)\pi(x)+O\left(\frac{x}{(\log x)^{9/8-\epsilon}}\right), \] where the implied constant is effectively computable and depends only on $k,N,$ and $\epsilon$.


title: Beyond Disagreement-based Agnostic Active Learning

abstract: We study agnostic active learning, where the goal is to learn a classifier in a pre-specified hypothesis class interactively with as few label queries as possible, while making no assumptions on the true function generating the labels. The main algorithms for this problem are {\em{disagreement-based active learning}}, which has a high label requirement, and {\em{margin-based active learning}}, which only applies to fairly restricted settings. A major challenge is to find an algorithm which achieves better label complexity, is consistent in an agnostic setting, and applies to general classification problems. In this paper, we provide such an algorithm. Our solution is based on two novel contributions -- a reduction from consistent active learning to confidence-rated prediction with guaranteed error, and a novel confidence-rated predictor.


title: Reconstructing quantum states from local data

abstract: We consider the problem of reconstructing global quantum states from local data. Because the reconstruction problem has many solutions in general, we consider the reconstructed state of maximum global entropy consistent with the local data. We show that unique ground states of local Hamiltonians are exactly reconstructed as the maximal entropy state. More generally, we show that if the state in question is a ground state of a local Hamiltonian with a degenerate space of locally indistinguishable ground states, then the maximal entropy state is close to the ground state projector. We also show that local reconstruction is possible for thermal states of local Hamiltonians. Finally, we discuss a procedure to certify that the reconstructed state is close to the true global state. We call the entropy of our reconstructed maximum entropy state the "reconstruction entropy", and we discuss its relation to emergent geometry in the context of holographic duality.


title: Moduli spaces of graded representations of finite dimensional algebras

abstract: Let $\Lambda$ be a basic finite dimensional algebra over an algebraically closed field, presented as a path algebra modulo relations; further, assume that $\Lambda$ is graded by lengths of paths. The paper addresses the classifiability, via moduli spaces, of classes of graded $\Lambda$-modules with fixed dimension $d$ and fixed top $T$. It is shown that such moduli spaces exist far more frequently than they do for ungraded modules. In the local case (i.e., when $T$ is simple), the graded $d$-dimensional $\Lambda$-modules with top $T$ always possess a fine moduli space which classifies these modules up to graded-isomorphism; moreover, this moduli space is a projective variety with a distinguished affine cover that can be constructed from quiver and relations of $\Lambda$. When $T$ is not simple, existence of a coarse moduli space for the graded $d$-dimensional $\Lambda$-modules with top $T$ forces these modules to be direct sums of local modules; under the latter condition, a finite collection of isomorphism invariants of the modules in question yields a partition into subclasses, each of which has a fine moduli space (again projective) parametrizing the corresponding graded-isomorphism classes.


title: High resolution numerical-relativity simulations for the merger of binary magnetized neutron stars

abstract: We perform high-resolution magnetohydrodynamics simulations of binary neutron star mergers in numerical relativity on the Japanese supercomputer K. The neutron stars and merger remnants are covered by a grid spacing of 70\,m, which yields the highest-resolution results among those derived so far. By an in-depth resolution study, we clarify several amplification mechanisms of magnetic fields during the binary neutron star merger for the first time. First, the Kelvin-Helmholtz instability developed in the shear layer at the onset of the merger significantly amplifies the magnetic fields. A hypermassive neutron star (HMNS) formed after the merger is then subject to the nonaxisymmetric magnetorotational instability, which amplifies the magnetic field in the HMNS. These two amplification mechanisms cannot be found with insufficient-resolution runs. We also show that the HMNS eventually collapses to a black hole surrounded by an accretion torus which is strongly magnetized at birth.

categories: astro-ph.HE astro-ph.CO astro-ph.SR gr-qc


title: Stacks of algebras and their homology

abstract: For any increasing function $f: {\Bbb N} \rightarrow {\Bbb N}_{\ge 2}$ which takes only finitely many distinct values, a connected finite dimensional algebra $\Lambda$ is constructed, with the property that $\text{fin.dim}_n\, \Lambda = f(n)$ for all $n$; here $\text{fin.dim}_n\, \Lambda$ is the $n$-generated finitistic dimension of $\Lambda$. The stacking technique developed for this construction of homological examples permits strong control over the higher syzygies of $\Lambda$-modules in terms of the algebras serving as layers.


title: Learning Privately with Labeled and Unlabeled Examples

abstract: A private learner is an algorithm that given a sample of labeled individual examples outputs a generalizing hypothesis while preserving the privacy of each individual. In 2008, Kasiviswanathan et al. (FOCS 2008) gave a generic construction of private learners, in which the sample complexity is (generally) higher than what is needed for non-private learners. This gap in the sample complexity was then further studied in several followup papers, showing that (at least in some cases) this gap is unavoidable. Moreover, those papers considered ways to overcome the gap, by relaxing either the privacy or the learning guarantees of the learner. We suggest an alternative approach, inspired by the (non-private) models of semi-supervised learning and active-learning, where the focus is on the sample complexity of labeled examples whereas unlabeled examples are of a significantly lower cost. We consider private semi-supervised learners that operate on a random sample, where only a (hopefully small) portion of this sample is labeled. The learners have no control over which of the sample elements are labeled. Our main result is that the labeled sample complexity of private learners is characterized by the VC dimension. We present two generic constructions of private semi-supervised learners. The first construction is of learners where the labeled sample complexity is proportional to the VC dimension of the concept class, however, the unlabeled sample complexity of the algorithm is as big as the representation length of domain elements. Our second construction presents a new technique for decreasing the labeled sample complexity of a given private learner, while roughly maintaining its unlabeled sample complexity. In addition, we show that in some settings the labeled sample complexity does not depend on the privacy parameters of the learner.


title: Using Realistic MHD Simulations for Modeling and Interpretation of Quiet-Sun Observations with the Solar Dynamics Observatory Helioseismic and Magnetic Imager

abstract: The solar atmosphere is extremely dynamic, and many important phenomena develop on small scales that are unresolved in observations with the Helioseismic and Magnetic Imager (HMI) instrument on the Solar Dynamics Observatory (SDO). For correct calibration and interpretation of the observations, it is very important to investigate the effects of small-scale structures and dynamics on the HMI observables, such as Doppler shift, continuum intensity, spectral line depth, and width. We use 3D radiative hydrodynamics simulations of the upper turbulent convective layer and the atmosphere of the Sun, and a spectro-polarimetric radiative transfer code to study observational characteristics of the Fe I 6173A line observed by HMI in quiet-Sun regions. We use the modeling results to investigate the sensitivity of the line Doppler shift to plasma velocity, and also sensitivities of the line parameters to plasma temperature and density, and determine effective line formation heights for observations of solar regions located at different distances from the disc center. These estimates are important for the interpretation of helioseismology measurements. In addition, we consider various center-to-limb effects, such as convective blue-shift, variations of helioseismic travel-times, and the 'concave' Sun effect, and show that the simulations can qualitatively reproduce the observed phenomena, indicating that these effects are related to a complex interaction of the solar dynamics and radiative transfer.


title: Classifying representations by way of Grassmannians

abstract: Let $\Lambda$ be a finite dimensional algebra over an algebraically closed field. Criteria are given which characterize existence of a fine or coarse moduli space classifying, up to isomorphism, the representations of $\Lambda$ with fixed dimension $d$ and fixed squarefree top $T$. Next to providing a complete theoretical picture, some of these equivalent conditions are readily checkable from quiver and relations. In case of existence of a moduli space -- unexpectedly frequent in light of the stringency of fine classification -- this space is always projective and, in fact, arises as a closed subvariety ${\mathfrak{Grass}}^T_d$ of a classical Grassmannian. Even when the full moduli problem fails to be solvable, the variety ${\mathfrak{Grass}}^T_d$ is seen to have distinctive properties recommending it as a substitute for a moduli space. As an application, a characterization of the algebras having only finitely many representations with fixed simple top is obtained; in this case of `finite local representation type at a given simple $T$', the radical layering $\bigl( J^lM/ J^{l+1}M \bigr)_{l \ge 0}$ is shown to be a classifying invariant for the modules with top $T$. This relies on the following general fact obtained as a byproduct: Proper degenerations of a local module $M$ never have the same radical layering as $M$.


title: Top-stable degenerations of finite dimensional representations I

abstract: Given a finite dimensional representation $M$ of a finite dimensional algebra, two hierarchies of degenerations of $M$ are analyzed in the context of their natural orders: the poset of those degenerations of $M$ which share the top $M/JM$ with $M$ - here $J$ denotes the radical of the algebra - and the sub-poset of those which share the full radical layering $\bigl(J^lM/J^{l+1}M\bigr)_{l \ge 0}$ with $M$. In particular, the article addresses existence of proper top-stable or layer-stable degenerations - more generally, it addresses the sizes of the corresponding posets including bounds on the lengths of saturated chains - as well as structure and classification.


title: The Occurrence of Wide-Orbit Planets in Binary Star Systems

abstract: The occurrence of planets in binary star systems has been investigated via a variety of techniques that sample a wide range of semi-major axes, but with a preponderance of such results applicable to planets with semi-major axes less than a few AU. We utilize a new method -- the presence or absence of heavy elements in the atmospheres of white dwarf stars -- to elucidate the frequency in main sequence binary star systems of planets with semi-major axes greater than a few AU. We consider only binaries where a putative planetary system orbits one member (no circumbinary planets). For main sequence binaries where the primary star is of spectral type A or F, data in the published literature suggests that the existence of a secondary star with a semi-major axis less than about 1000 AU suppresses the formation and/or long term stability of an extended planetary system around the primary. For these spectral types and initial semi-major axis $\geq$1000 AU, extended planetary systems appear to be as common around stars in binary systems as they are around single stars.


title: Information Geometry for Husimi-Temperley Model

abstract: We examine phase transition of the Husimi-Temperley model in terms of information geometry. For this purpose, we introduce the Fisher metric defined by the density matrix of the model. We find that the metric becomes hyperbolic at the critical point with respect to the energy scale. Then, the metric is invariant under the scale transformation. We also find that the equation of states is naturally derived from a necessary condition for the entropy operator that is a building block of the metric. Based on these findings, we conclude that the geometric quantities clearly detect the phase transition of the model.


title: Transport and Capacitance properties of Charge Density Wave in few layer 2H-TaS2 Devices

abstract: We carefully investigated the transport and capacitance properties of few layer charge density wave (CDW) 2H-TaS2 devices. The CDW transition temperature and the threshold voltage vary from device to device, which is attributed to the interlayer interaction and inhomogeneous local defects of these micro-devices based on few layer 2H-TaS2 flakes. Semiconductivity rather than metallic property of 2H-TaS2 devices was observed in our experiment at low temperature. The temperature dependence of the relative threshold voltage can be scaled to (1- T / Tr )^0.5+delta with delta=0.08 for the different measured devices with presence of the CDWs. The conductance-voltage and capacity-voltage measurements were performed simultaneously. At very low ac active voltage, we found that the hysteresis loops of these two measurements exactly match each other. Our results point out that the capacity-voltage measurements can also be used to define the threshold depinning voltage of the CDW, which give us a new method to investigate the CDWs.


title: Elastic Properties of Chemical-Vapor-Deposited Monolayer MoS2, WS2, and Their Bilayer Heterostructures

abstract: Elastic properties of materials are an important factor in their integration in applications. Chemical vapor deposited (CVD) monolayer semiconductors are proposed as key components in industrial-scale flexible devices and building blocks of 2D van der Waals heterostructures. However, their mechanical and elastic properties have not been fully characterized. Here we report high 2D elastic moduli of CVD monolayer MoS2 and WS2 (~ 170 N/m), which is very close to the value of exfoliated MoS2 monolayers and almost half the value of the strongest material, graphene. The 2D moduli of their bilayer heterostructures are lower than the sum of 2D modulus of each layer, but comparable to the corresponding bilayer homostructure, implying similar interactions between the hetero monolayers as between homo monolayers. These results not only provide deep insight to understanding interlayer interactions in 2D van der Waals structures, but also potentially allow engineering of their elastic properties as desired.


title: Thermal design and performance of the balloon-borne large aperture submillimeter telescope for polarimetry BLASTPol

abstract: We present the thermal model of the Balloon-borne Large-Aperture Submillimeter Telescope for Polarimetry (BLASTPol). This instrument was successfully flown in two circumpolar flights from McMurdo, Antarctica in 2010 and 2012. During these two flights, BLASTPol obtained unprecedented information about the magnetic field in molecular clouds through the measurement of the polarized thermal emission of interstellar dust grains. The thermal design of the experiment addresses the stability and control of the payload necessary for this kind of measurement. We describe the thermal modeling of the payload including the sun-shielding strategy. We present the in-flight thermal performance of the instrument and compare the predictions of the model with the temperatures registered during the flight. We describe the difficulties of modeling the thermal behavior of the balloon-borne platform and establish landmarks that can be used in the design of future balloon-borne instruments.


title: Nuclear spin content and constraints on exotic spin-dependent couplings

abstract: There are numerous recent and ongoing experiments employing a variety of atomic species to search for couplings of atomic spins to exotic fields. In order to meaningfully compare these experimental results, the coupling of the exotic field to the atomic spin must be interpreted in terms of the coupling to electron, proton, and neutron spins. Traditionally, constraints from atomic experiments on exotic couplings to neutron and proton spins have been derived using the single-particle Schmidt model for nuclear spin. In this model, particular atomic species are sensitive to either neutron or proton spin couplings, but not both. More recently, semi-empirical models employing nuclear magnetic moment data have been used to derive new constraints for non-valence nucleons. However, comparison of such semi-empirical models to detailed large-scale nuclear shell model calculations and analysis of known physical effects in nuclei show that existing semi-empirical models cannot reliably be used to predict the spin polarization of non-valence nucleons. The results of our re-analysis of nuclear spin content are applied to searches for exotic long-range monopole-dipole and dipole-dipole couplings of nuclei leading to significant revisions of some published constraints.


title: Truncated path algebras are homologically transparent

abstract: It is shown that path algebras modulo relations of the form $\Lambda = KQ/I$, where $Q$ is a quiver, $K$ a coefficient field, and $I \subseteq KQ$ the ideal generated by all paths of a given length, can be readily analyzed homologically, while displaying a wealth of phenomena. In particular, the syzygies of their modules, and hence their finitistic dimensions, allow for smooth descriptions in terms of $Q$ and the Loewy length of $\Lambda$. The same is true for the distributions of projective dimensions attained on the irreducible components of the standard parametrizing varieties for the modules of fixed $K$- dimension.


title: Generic representation theory of quivers with relations

abstract: The irreducible components of varieties parametrizing the finite dimensional representations of a finite dimensional algebra $\Lambda$ are explored, with regard to both their geometry and the structure of the modules they encode. Provided that the base field is algebraically closed and has infinite transcendence degree over its prime field, we establish the existence and uniqueness (not up to isomorphism, but in a strong sense to be specified) of a generic module for any irreducible component $\mathcal C$, that is, of a module which displays all categorically defined generic properties of the modules parametrized by $\mathcal C$; the crucial point of the existence statement - a priori almost obvious - lies in the description of such a module in a format accessible to representation-theoretic techniques. Our approach to generic modules over path algebras modulo relations, by way of minimal projective resolutions, is largely constructive. It is explicit for large classes of algebras of wild type. We follow with an investigation of the properties of such generic modules in terms of quiver and relations. The sharpest specific results on all fronts are obtained for truncated path algebras, that is, for path algebras of quivers modulo ideals generated by all paths of a fixed length; this class of algebras extends the comparatively thoroughly studied hereditary case, displaying many novel features.


title: Private Learning and Sanitization: Pure vs. Approximate Differential Privacy

abstract: We compare the sample complexity of private learning [Kasiviswanathan et al. 2008] and sanitization~[Blum et al. 2008] under pure $\epsilon$-differential privacy [Dwork et al. TCC 2006] and approximate $(\epsilon,\delta)$-differential privacy [Dwork et al. Eurocrypt 2006]. We show that the sample complexity of these tasks under approximate differential privacy can be significantly lower than that under pure differential privacy. We define a family of optimization problems, which we call Quasi-Concave Promise Problems, that generalizes some of our considered tasks. We observe that a quasi-concave promise problem can be privately approximated using a solution to a smaller instance of a quasi-concave promise problem. This allows us to construct an efficient recursive algorithm solving such problems privately. Specifically, we construct private learners for point functions, threshold functions, and axis-aligned rectangles in high dimension. Similarly, we construct sanitizers for point functions and threshold functions. We also examine the sample complexity of label-private learners, a relaxation of private learning where the learner is required to only protect the privacy of the labels in the sample. We show that the VC dimension completely characterizes the sample complexity of such learners, that is, the sample complexity of learning with label privacy is equal (up to constants) to learning without privacy.


title: A hierarchy of parametrizing varieties for representations

abstract: The primary purpose is to introduce and explore projective varieties, $\text{GRASS}_{\bf d}(\Lambda)$, parametrizing the full collection of those modules over a finite dimensional algebra $\Lambda$ which have dimension vector $\bf d$. These varieties extend the smaller varieties previously studied by the author; namely, the projective varieties encoding those modules with dimension vector $\bf d$ which, in addition, have a preassigned top or radical layering. Each of the $\text{GRASS}_{\bf d}(\Lambda)$ is again partitioned by the action of a linear algebraic group, and covered by certain representation-theoretically defined affine subvarieties which are stable under the unipotent radical of the acting group. A special case of the pertinent theorem served as a cornerstone in the work on generic representations by Babson, Thomas, and the author. Moreover, applications are given to the study of degenerations.


title: A New Optimal Stepsize For Approximate Dynamic Programming

abstract: Approximate dynamic programming (ADP) has proven itself in a wide range of applications spanning large-scale transportation problems, health care, revenue management, and energy systems. The design of effective ADP algorithms has many dimensions, but one crucial factor is the stepsize rule used to update a value function approximation. Many operations research applications are computationally intensive, and it is important to obtain good results quickly. Furthermore, the most popular stepsize formulas use tunable parameters and can produce very poor results if tuned improperly. We derive a new stepsize rule that optimizes the prediction error in order to improve the short-term performance of an ADP algorithm. With only one, relatively insensitive tunable parameter, the new rule adapts to the level of noise in the problem and produces faster convergence in numerical experiments.


title: Impacts of Regional Trade Agreements(RTAs) on Food Security: A Case of ASEAN Free Trade Agreement

abstract: Discriminatory trade liberalization policies are becoming more popular among world economies. Countries are motivated to enter for regional trade agreements to capture faster economic growth for alleviating poverty. In developing economies like most of the member countries of the Association of South East Asian Nations (ASEAN), a sizeable portion of people are suffering from poverty by exposing them to food insecurity. Low level of income and low productivity of agricultural sector have augmented the severity of food insecurity of those people. Discriminatory trade liberalization policies are expected to reduce poverty and strengthen the food security. The objective of this paper is to examine the effect of ASEAN Free Trade Agreement (AFTA) on food security of its member countries. The multiple regression analysis in panel data was employed to disentangle the impacts of trade liberalization on food securit y with use of regional trade agreement dummy variable. The finding of the study supports that AFTA has influenced positively on food security of its member nations. After the formation of AFTA, the level of per-capita daily dietary energy supply of the member countries has been increased moderately over time.


title: Study of Simulation Method of Time Evolution of Atomic and Molecular Systems by Quantum Electrodynamics

abstract: We discuss a method to follow step-by-step time evolution of atomic and molecular systems based on QED (Quantum Electrodynamics). Our strategy includes expanding the electron field operator by localized wavepackets to define creation and annihilation operators and following the time evolution using the equations of motion of the field operator in the Heisenberg picture. We first derive a time evolution equation for the excitation operator, the product of two creation or annihilation operators, which is necessary for constructing operators of physical quantities such as the electronic charge density operator. We then describe our approximation methods to obtain time differential equations of the electronic density matrix, which is defined as the expectation value of the excitation operator. By solving the equations numerically, we show "electron-positron oscillations", the fluctuations originated from virtual electron-positron pair creations and annihilations, appear in the charge density of a hydrogen atom and molecule. We also show that the period of the electron-positron oscillations becomes shorter by including the self-energy process, in which the electron emits a photon and then absorbs it again, and it can be interpreted as the increase in the electron mass due to the self-energy.


title: Smaller SDP for SOS Decomposition

abstract: A popular numerical method to compute SOS (sum of squares of polynomials) decompositions for polynomials is to transform the problem into semi-definite programming (SDP) problems and then solve them by SDP solvers. In this paper, we focus on reducing the sizes of inputs to SDP solvers to improve the efficiency and reliability of those SDP based methods. Two types of polynomials, convex cover polynomials and split polynomials, are defined. A convex cover polynomial or a split polynomial can be decomposed into several smaller sub-polynomials such that the original polynomial is SOS if and only if the sub-polynomials are all SOS. Thus the original SOS problem can be decomposed equivalently into smaller sub-problems. It is proved that convex cover polynomials are split polynomials and it is quite possible that sparse polynomials with many variables are split polynomials, which can be efficiently detected in practice. Some necessary conditions for polynomials to be SOS are also given, which can help refute quickly those polynomials which have no SOS representations so that SDP solvers are not called in this case. All the new results lead to a new SDP based method to compute SOS decompositions, which improves this kind of methods by passing smaller inputs to SDP solvers in some cases. Experiments show that the number of monomials obtained by our program is often smaller than that by other SDP based software, especially for polynomials with many variables and high degrees. Numerical results on various tests are reported to show the performance of our program.