title: The Next Generation Virgo Cluster Survey. VI. The Kinematics of Ultra-compact Dwarfs and Globular Clusters in M87

abstract: The origin of ultra-compact dwarfs (UCDs)--objects larger and more massive than typical globular clusters (GCs), but more compact than typical dwarf galaxies--has been hotly debated in the 15 years since their discovery. Even whether UCDs should be considered galactic in origin, or simply the most extreme GCs, is not yet settled. We present the dynamical properties of 97 spectroscopically confirmed UCDs (rh >~10 pc) and 911 GCs associated with central cD galaxy of the Virgo cluster, M87. Our UCDs, of which 89% have M_star > ~2X10^6 M_sun and 92% are as blue as the classic blue GCs, nearly triple the sample of previous confirmed Virgo UCDs, providing by far the best opportunity for studying the global dynamics of a UCD system. We found that (1) UCDs have a surface number density profile that is shallower than that of the blue GCs in the inner ~ 70 kpc and as steep as that of the red GCs at larger radii; (2) UCDs exhibit a significantly stronger rotation than the GCs, and the blue GCs seem to have a velocity field that is more consistent with that of the surrounding dwarf ellipticals than with that of UCDs; (3) UCDs have a radially increasing orbital anisotropy profile, and are tangentially-biased at radii < ~ 40 kpc and radially-biased further out. In contrast, the blue GCs become more tangentially-biased at larger radii beyond ~ 40 kpc; (4) GCs with M_star > 2X10^6 M_sun have rotational properties indistinguishable from the less massive ones, suggesting that it is the size, instead of mass, that differentiates UCDs from GCs as kinematically distinct populations. We conclude that most UCDs in M87 are not consistent with being merely the most luminous and extended examples of otherwise normal GCs. The radially-biased orbital structure of UCDs at large radii is in general agreement with the "tidally threshed dwarf galaxy" scenario.


title: Complex organic molecules in organic-poor massive young stellar objects

abstract: Massive young stellar objects (MYSOs) with hot cores are classic sources of complex organic molecules. The origins of these molecules in such sources, as well as the small- and large-scale differentiation between nitrogen- and oxygen-bearing complex species, are poorly understood. We aim to use complex molecule abundances toward a chemically less explored class of MYSOs with weak hot organic emission lines to constrain the impact of hot molecular cores and initial ice conditions on the chemical composition toward MYSOs. We use the IRAM 30m and the Submillimeter Array to search for complex organic molecules over 8-16 GHz in the 1~mm atmospheric window toward three MYSOs with known ice abundances, but without luminous molecular hot cores. Complex molecules are detected toward all three sources at comparable abundances with respect to CH$_3$OH to classical hot core sources. The relative importance of CH$_3$CHO, CH$_3$CCH, CH$_3$OCH$_3$, CH$_3$CN, and HNCO differ between the organic-poor MYSOs and hot cores, however. Furthermore, the N-bearing molecules are generally concentrated toward the source centers, while most O- and C-bearing molecules are present both in the center and in the colder envelope. Gas-phase HNCO/CH$_3$OH ratios are tentatively correlated with the ratios of NH$_3$ ice over CH$_3$OH ice in the same lines of sight, which is consistent with new gas-grain model predictions. Hot cores are not required to form complex organic molecules, and source temperature and initial ice composition both seem to affect complex organic distributions toward MYSOs. To quantify the relative impact of temperature and initial conditions requires, however, a larger spatially resolved survey of MYSOs with ice detections.


title: Implicit large eddy simulations of anisotropic weakly compressible turbulence with application to core-collapse supernovae

abstract: (Abridged) In the implicit large eddy simulation (ILES) paradigm, the dissipative nature of high-resolution shock-capturing schemes is exploited to provide an implicit model of turbulence. Recent 3D simulations suggest that turbulence might play a crucial role in core-collapse supernova explosions, however the fidelity with which turbulence is simulated in these studies is unclear. Especially considering that the accuracy of ILES for the regime of interest in CCSN, weakly compressible and strongly anisotropic, has not been systematically assessed before. In this paper we assess the accuracy of ILES using numerical methods most commonly employed in computational astrophysics by means of a number of local simulations of driven, weakly compressible, anisotropic turbulence. We report a detailed analysis of the way in which the turbulent cascade is influenced by the numerics. Our results suggest that anisotropy and compressibility in CCSN turbulence have little effect on the turbulent kinetic energy spectrum and a Kolmogorov $k^{-5/3}$ scaling is obtained in the inertial range. We find that, on the one hand, the kinetic energy dissipation rate at large scales is correctly captured even at relatively low resolutions, suggesting that very high effective Reynolds number can be achieved at the largest scales of the simulation. On the other hand, the dynamics at intermediate scales appears to be completely dominated by the so-called bottleneck effect, \ie the pile up of kinetic energy close to the dissipation range due to the partial suppression of the energy cascade by numerical viscosity. An inertial range is not recovered until the point where relatively high resolution $\sim 512^3$, which would be difficult to realize in global simulations, is reached. We discuss the consequences for CCSN simulations.


title: On the characterization of the numbers $n$ such that any group of order $n$ has a given property $P$

abstract: One of the classical problems in group theory is determining the set of positive integers $n$ such that every group of order $n$ has a particular property $P$, such as cyclic or abelian. We first present the Sylow theorems and the idea of solvable groups, both of which will be invaluable in our analysis. We then gather various solutions to this problem for cyclic, abelian, nilpotent, and supersolvable groups, as well as groups with ordered Sylow towers. This work is an exposition of known results, but it is hoped that the reader will find useful the presentation in a single account of the various tools that have been used to solve this general problem. This article claims no originality, but is meant as a synthesis of related knowledge and resources.


title: Assembly history of subhalo populations in galactic and cluster sized dark haloes

abstract: We make use of two suits of ultra high resolution N-body simulations of individual dark matter haloes from the Phoenix and the Aquarius Projects to investigate systematics of assembly history of subhaloes in dark matter haloes differing by a factor of $1000$ in the halo mass. We have found that real progenitors which built up present day subhalo population are relatively more abundant for high mass haloes, in contrast to previous studies claiming a universal form independent of the host halo mass. That is mainly because of repeated counting of the 're-accreted' (progenitors passed through and were later re-accreted to the host more than once) and inclusion of the 'ejected' progenitor population(progenitors were accreted to the host in the past but no longer members at present day) in previous studies. The typical accretion time for all progenitors vary strongly with the host halo mass, which is typical about $z \sim 5$ for the galactic Aquarius and about $z \sim 3$ for the cluster sized Phoenix haloes. Once these progenitors start to orbit their parent haloes, they rapidly lose their original mass but not their identifiers, more than $55$ ($50$) percent of them survive to present day for the Phoenix(Aquarius) haloes. At given redshift, survival fraction of the accreted subhalo is independent of the parent halo mass, whilst the mass-loss of the subhalo is more efficient in high mass haloes. These systematics results in similarity and difference in the subhalo population in dark matter haloes of different masses at present day.


title: Kinematic and Spatial Substructure in NGC 2264

abstract: We present an expanded kinematic study of the young cluster NGC 2264 based upon optical radial velocities measured using multi-fiber echelle spectroscopy at the 6.5 meter MMT and Magellan telescopes. We report radial velocities for 695 stars, of which approximately 407 stars are confirmed or very likely members. Our results more than double the number of members with radial velocities from F{\H u}r{\'e}sz et al., resulting in a much better defined kinematic relationship between the stellar population and the associated molecular gas. In particular, we find that there is a significant subset of stars that are systematically blueshifted with respect to the molecular ($^{13}$CO) gas. The detection of Lithium absorption and/or infrared excesses in this blue-shifted population suggests that at least some of these stars are cluster members; we suggest some speculative scenarios to explain their kinematics. Our results also more clearly define the redshifted population of stars in the northern end of the cluster; we suggest that the stellar and gas kinematics of this region are the result of a bubble driven by the wind from O7 star S Mon. Our results emphasize the complexity of the spatial and kinematic structure of NGC 2264, important for eventually building up a comprehensive picture of cluster formation.


title: A note on choice and detect probabilities in the presence of choice bias

abstract: Recently we have presented the analytical relationship between choice probabilities, noise correlations and read-out weights in the classical feedforward decision-making framework (Haefner et al. 2013). The derivation assumed that behavioral reports are distributed evenly between the two possible choices. This assumption is often violated in empirical data - especially when computing so-called grand CPs combining data across stimulus conditions. Here, we extend our analytical results to situations when subjects show clear biases towards one choice over the other, e.g. in non-zero signal conditions. Importantly, this also extends our results from discrimination tasks to detection tasks and detect probabilities for which much empirical data is available. We find that CPs and DPs depend monotonously on the fraction, p, of choices assigned to the more likely option: CPs and DPs are smallest for p equal to 0.5 and increase as p increases, i.e. as the data deviates from the ideal, zero-signal, unbiased scenario. While this deviation is small, our results suggest a) an empirical test for the feedforward framework and b) a way in which to correct choice probability and detect probability measurements before combining different stimulus conditions to increase signal/noise.


title: Steady State and Modulated Heat Conduction in Layered Systems Predicted by the Phonon Boltzmann Transport Equation

abstract: Based on the phonon Boltzmann transport equation under the relaxation time approximation, analytical expressions for the temperature profiles of both steady state and modulated heat conduction inside a thin film deposited on a substrate are derived and analyzed. It is shown that both steady state and modulated components of the temperature depend strongly on the ratio between the film thickness and the average phonon mean free path, and they exhibit the diffusive behavior as predicted by the Fourier law of heat conduction when this ratio is much larger than the unity. In contrast, in the ballistic regime when this ratio is comparable to or smaller than the unity, the steady-state temperature tends to be independent of position, while the amplitude and the phase of the modulated temperature appear to be lower than those determined by the Fourier law. Furthermore, we derived an invariant of heat conduction and a simple formula for the cross-plane thermal conductivity of dielectric thin films, which could be a useful guide for understanding and optimizing the thermal performance of the layered systems. This work represents the Boltzmann transport equation-based extension of the Rosencwaig and Gerko work [J. Appl. Phys. 47, 64 (1976)], which is based on the Fourier law and has widely been used as the theoretical framework for the development of photoacoustic and photothermal techniques. This work might shed some light on developing a theoretical basis for the determination of the phonon MFP and relaxation time using ultrafast laser-based transient heating techniques.


title: The interaction between supermassive black holes and globular clusters

abstract: Almost all galaxies along the Hubble sequence host a compact massive object (CMO) in their center. The CMO can be either a supermassive black hole (SMBH) or a very dense stellar cluster, also known as nuclear star cluster (NSC). Generally, heavier galaxies (mass $\gtrsim 10^{11}\mathrm{M}_{\odot}$ ) host a central SMBH while lighter show a central NSC. Intermediate mass hosts, instead, contain both a NSC and a SMBH. One possible formation mechanisms of a NSC relies on the dry-merger (migratory) scenario, in which globular clusters (GCs) decay toward the center of the host galaxy and merge. In this framework, the absence of NSCs in high-mass galaxies can be imputed to destruction of the infalling GCs by the intense tidal field of the central SMBH. In this work, we report preliminary results of N-body simulations performed using our high-resolution, direct, code HiGPUs, to investigate the effects of a central SMBH on a single GC orbiting around it. By varying either the mass of the SMBH and the mass of the host galaxy, we derived an upper limit to the mass of the central SMBH, and thus to the mass of the host, above which the formation of a NSC is suppressed.


title: Influence of the electron-phonon interfacial conductance on the thermal transport at metal/dielectric interfaces

abstract: Thermal boundary conductance at a metal-dieletric interface is a quantity of prime importance for heat management at the nanoscale. While the boundary conductance is usually ascribed to the coupling between metal phonons and dielectric phonons, in this work we examine the influence of a direct coupling between the metal electrons and the dielectric phonons. The effect of electron- phonon processes is generally believed to be resistive, and tends to decrease the overall thermal boundary conductance as compared to the phonon-phonon conductance {\sigma}p . Here, we find that the effect of a direct coupling {\sigma}e is to enhance the effective thermal conductance, between the metal and the dielectric. Resistive effects turn out to be important only for thin films of metals having a low electron-phonon coupling strength. Two approaches are explored to reach these conclusions. First, we present an analytical solution of the two-temperature model to compute the effective conductance which account for all the relevant energy channels, as a function of {\sigma}e , {\sigma}p and the electron-phonon coupling factor G. Second, we use numerical resolution to examine the influence of {\sigma}e on two realistic cases: gold film on silicon or silica substrates. We point out the implications for the interpretation of time-resolved thermoreflectance experiments.


title: Looking for Population III stars with He II line intensity mapping

abstract: Constraining the properties of Population III (Pop III) stars will be very challenging because they reside in small galaxies at high redshift which will be difficult to directly detect. In this paper, we suggest that intensity mapping may be a promising method to study Pop III stars. Intensity mapping is a technique proposed to measure large-scale fluctuations of galaxy line emission in three dimensions without resolving individual sources. This technique is well suited for observing many faint galaxies because it can measure their cumulative emission even if they cannot be directly detected. We focus on intensity mapping of He II recombination lines. These lines are much stronger in Pop III stars than Pop II stars because the harder spectra of Pop III stars are expected to produce many He II ionizing photons. Measuring the He II 1640 \AA{} intensity mapping signal, along with the signals from other lines such as Ly$\alpha$, H$\alpha$, and metal lines, could give constraints on the initial mass function (IMF) and star formation rate density (SFRD) of Pop III stars as a function of redshift. To demonstrate the feasibility of these observations, we estimate the strength of the Pop III He II 1640 \AA{} intensity mapping signal from $z=10-20$. We show that at $z\approx10$, the signal could be measured accurately by two different hypothetical future instruments, one which cross-correlates He II 1640 \AA{} with CO(1-0) line emission from galaxies and the other with 21 cm emission from the intergalactic medium (IGM).


title: Metric compactification of infinite Sierpi\'nski carpet graphs

abstract: We associate, with every infinite word over a finite alphabet, an increasing sequence of rooted finite graphs, which provide a discrete approximation of the famous Sierpi\'nski carpet fractal. Each of these sequences converges, in the Gromov-Hausdorff topology, to an infinite rooted graph. We give an explicit description of the metric compactification of each of these limit graphs. In particular, we are able to classify Busemann and non-Busemann points of the metric boundary. It turns out that, with respect to the uniform Bernoulli measure on the set of words indexing the graphs, for almost all the infinite graphs, the boundary consists of four Busemann points and countably many non-Busemann points.


title: We Did See This Coming: Response to, We Should Have Seen This Coming, by D. Sam Schwarzkopf

abstract: We appreciate the effort by Schwarzkopf to examine alternative explanations for predictive anticipatory activity (PAA) or presentiment (for first response, see: Schwarzkopf 2014a; for additional response, see: Schwarzkopf 2014b, for original article, see: Mossbridge et al. 2014). These commentaries are a laudable effort to promote collegial discussion of the controversial claim of presentiment, whereby physiological measures preceding unpredictable emotional events differ from physiological measures preceding calm or neutral events (Mossbridge et al., 2012; Mossbridge et al., 2014). What is called truth at any given time in science has achieved that status through a continuous process of measurement and interpretation based on the current knowledge at hand. Here we address six points in his original commentary (Schwarzkopf 2014a), though our responses are informed by the points he made in his his supplementary commentary (Schwarzkopf 2014b). We hope our responses will help Schwarzkopf and others understand our interpretation of these data.


title: A 400 solar mass black hole in the Ultraluminous X-ray source M82 X-1 accreting close to its Eddington limit

abstract: M82 X-1, the brightest X-ray source in the galaxy M82, has been thought to be an intermediate-mass black hole (100 to 10,000 solar masses) because of its extremely high luminosity and variability characteristics, although some models suggest that its mass may be only about 20 solar masses. The previous mass estimates were based on scaling relations that use low-frequency characteristic timescales which have large intrinsic uncertainties. For stellar-mass black holes, we know that the high-frequency quasi-periodic oscillations (100-450 hertz) in the X-ray emission that occur in a 3:2 frequency ratio are stable and scale in frequency inversely with black hole mass with a reasonably small dispersion. The discovery of such stable oscillations thus potentially offers an alternative and less ambiguous means of mass determination for intermediate-mass black holes, but has hitherto not been realized. Here we report stable, twin-peak (3:2 frequency ratio) X-ray quasi-periodic oscillations from M82 X-1 at frequencies of 3.32$\pm$0.06 hertz and 5.07$\pm$0.06 hertz. Assuming that we can extrapolate the inverse-mass scaling that holds for stellar-mass black holes, we estimate the black hole mass of M82 X-1 to be 428$\pm$105 solar masses. In addition, we can estimate the mass using the relativistic precession model, from which we get a value of 415$\pm$63 solar masses.


title: The Inner Disk Structure, Disk-Planet Interactions, and Temporal Evolution in the Beta Pictoris System: A Two-Epoch HST/STIS Coronagraphic Study

abstract: We present deep HST/STIS coronagraphic images of the Beta Pic debris disk obtained at two epochs separated by 15 years. The new images and the re-reduction of the 1997 data provide the most sensitive and detailed views of the disk at optical wavelengths as well as the yet smallest inner working angle optical coronagraphic image of the disk. Our observations characterize the large-scale and inner-disk asymmetries and we identify multiple breaks in the disk radial surface brightness profile. We study in detail the radial and vertical disk structure and show that the disk is warped. We explore the disk at the location of the Beta Pic b super-jupiter and find that the disk surface brightness slope is continuous between 0.5 and 2.0 arcsec, arguing for no change at the separations where Beta Pic b orbits. The two epoch images constrain the disk surface brightness evolution on orbital and radiation pressure blow-out timescales. We place an upper limit of 3% on the disk surface brightness change between 3-5 arcsec, including the locations of the disk warp, and the CO and dust clumps. We discuss the new observations in the context of high-resolution multi-wavelength images and divide the disk asymmetries in two groups: axisymmetric and non-axisymmetric. The axisymmetric structures (warp, large-scale butterfly, etc.) are consistent with disk structure models that include interactions of a planetesimal belt and a non-coplanar giant planet. The non-axisymmetric features, however, require a different explanation.


title: Curvature properties of $4$-dimensional Riemannian manifolds with a circulant structure

abstract: We consider a $4$-dimensional Riemannian manifold $M$ equip\-ped with a circulant structure $q$, which is an isometry with respect to the metric $g$ and $q^{4}=\id$, $q^{2}\neq \pm \id$. For such a manifold $(M, g, q)$ we obtain some assertions for the sectional curvatures of $2$-planes. We construct an example of such a manifold on a Lie group and we find some of its geometric characteristics.


title: Lifting Commutation Relations in Cuntz Algebras

abstract: We examine splitting of the quotient map from the full free product $A*B$, or the unital free product $A*_{\mathbb C}B$, to the (maximal) tensor product $A\otimes B$, for unital C*-algebras $A$ and $B$. Such a splitting is very rare, but we show there is one if $A$ and $B$ are both the Cuntz algebra $O_2$ or $O_\infty$, and in a few other cases. The splitting is not explicit (and in principle probably cannot be). We also describe severe $K$-theoretic obstructions to a splitting.


title: Self-consistent hybrid functional for condensed systems

abstract: A self-consistent scheme for determining the optimal fraction of exact exchange for full-range hybrid functionals is presented and applied to the calculation of band gaps and dielectric constants of solids. The exchange-correlation functional is defined in a similar manner to the PBE0 functional, but the mixing parameter is set equal to the inverse macroscopic dielectric function and it is determined self-consistently by computing the optimal dielectric screening. We found excellent agreement with experiments for the properties of a broad class of systems, with band gaps ranging between 0.7 and 21.7 eV and dielectric constants within 1.23 and 15.9. We propose that the eigenvalues and eigenfunctions obtained with the present self-consistent hybrid scheme may be excellent inputs for G$_0$W$_0$ calculations.


title: Post-Selection and Post-Regularization Inference in Linear Models with Many Controls and Instruments

abstract: In this note, we offer an approach to estimating causal/structural parameters in the presence of many instruments and controls based on methods for estimating sparse high-dimensional models. We use these high-dimensional methods to select both which instruments and which control variables to use. The approach we take extends BCCH2012, which covers selection of instruments for IV models with a small number of controls, and extends BCH2014, which covers selection of controls in models where the variable of interest is exogenous conditional on observables, to accommodate both a large number of controls and a large number of instruments. We illustrate the approach with a simulation and an empirical example. Technical supporting material is available in a supplementary online appendix.


title: Correlation effects in the capacitance of a gated carbon nanotube

abstract: For a capacitor made of a semiconducting carbon nanotube (CNT) suspended above a metallic gate, Coulomb correlations between individual electrons can lead to a capacitance that is much larger than the geometric capacitance. We argue that when the average spacing $n^{-1}$ between electrons within the low density 1-dimensional electron gas (1DEG) in the CNT is larger than the physical separation $d$ between the CNT and the gate, the enhancement of capacitance is expected to be big. A recent Coulomb blockade experiment[1], however, has observed no obvious increase of capacitance even at very low electron density. We show that this smaller capacitance can be understood as the result of the confining potential produced by the potential difference between the source/drain electrodes and the gate, which compresses the 1DEG when the electron number decreases. We suggest that by profiling the potential with the help of multiple split gates, one can return to the case of a uniform 1DEG with anomalously large capacitance.


title: Subspaces with extra invariance nearest to observed data

abstract: Given an arbitrary finite set of data F= {f_1,..., f_m} in L2(Rd) we prove the existence and show how to construct a "small shift invariant space" that is "closest" to the data F over certain class of closed subspaces of L2(Rd). The approximating subspace is required to have extra-invariance properties, that is to be invariant under translations by a prefixed additive subgroup of Rd containing Zd. This is important for example in situations where we need to deal with jitter error of the data. Here small means that our solution subspace should be generated by the integer translates of a small number of generators. We give an expression for the error in terms of the data and construct a Parseval frame for the optimal space. We also consider the problem of approximating F from generalised Paley-Wiener spaces of Rd that are generated by the integer translates of a finite number of functions. That is finitely generated shift invariant spaces that are translation invariant. We characterise these spaces in terms of multi-tile sets of Rd, and show the connections with recent results on Riesz basis of exponentials on bounded sets of Rd. Finally we study the discrete case for our approximation problem.


title: Video Manipulation Techniques for the Protection of Privacy in Remote Presence Systems

abstract: Systems that give control of a mobile robot to a remote user raise privacy concerns about what the remote user can see and do through the robot. We aim to preserve some of that privacy by manipulating the video data that the remote user sees. Through two user studies, we explore the effectiveness of different video manipulation techniques at providing different types of privacy. We simultaneously examine task performance in the presence of privacy protection. In the first study, participants were asked to watch a video captured by a robot exploring an office environment and to complete a series of observational tasks under differing video manipulation conditions. Our results show that using manipulations of the video stream can lead to fewer privacy violations for different privacy types. Through a second user study, it was demonstrated that these privacy-protecting techniques were effective without diminishing the task performance of the remote user.


title: Propagation of Input Uncertainty in Presence of Model-Form Uncertainty: A Multi-fidelity Approach for CFD Applications

abstract: Proper quantification and propagation of uncertainties in computational simulations are of critical importance. This issue is especially challenging for CFD applications. A particular obstacle for uncertainty quantifications in CFD problems is the large model discrepancies associated with the CFD models used for uncertainty propagation. Neglecting or improperly representing the model discrepancies leads to inaccurate and distorted uncertainty distribution for the Quantities of Interest. High-fidelity models, being accurate yet expensive, can accommodate only a small ensemble of simulations and thus lead to large interpolation errors and/or sampling errors; low-fidelity models can propagate a large ensemble, but can introduce large modeling errors. In this work, we propose a multi-model strategy to account for the influences of model discrepancies in uncertainty propagation and to reduce their impact on the predictions. Specifically, we take advantage of CFD models of multiple fidelities to estimate the model discrepancies associated with the lower-fidelity model in the parameter space. A Gaussian process is adopted to construct the model discrepancy function, and a Bayesian approach is used to infer the discrepancies and corresponding uncertainties in the regions of the parameter space where the high-fidelity simulations are not performed. The proposed multi-model strategy combines information from models with different fidelities and computational costs, and is of particular relevance for CFD applications, where a hierarchy of models with a wide range of complexities exists. Several examples of relevance to CFD applications are performed to demonstrate the merits of the proposed strategy. Simulation results suggest that, by combining low- and high-fidelity models, the proposed approach produces better results than what either model can achieve individually.


title: Light-Element Abundances of Giant Stars in the Globular Cluster M71 (NGC 6838)

abstract: Aluminum is the heaviest light element displaying large star--to--star variations in Galactic globular clusters (GCs). This element may provide additional insight into the origin of the multiple populations now known to be a common place in GCs, and also the nature of the first-generation stars responsible for a cluster's chemical inhomogeneities. In a previous analysis we found that, unlike more metal-poor GCs, 47 Tuc did not exhibit a strong Na-Al correlation, which motivates a careful study of the similar metallicity but less massive GC M71. We present chemical abundances of O, Na, Al, and Fe for 33 giants in M71 using spectra obtained with the WIYN-Hydra spectrograph. Our spectroscopic analysis finds that, similar to 47 Tuc and in contrast with more metal-poor GCs, M71 stars do not exhibit a strong Na-Al correlation and span a relatively narrow range in [Al/Fe], which are characteristics that GC formation models must reproduce.


title: Annotating Cognates and Etymological Origin in Turkic Languages

abstract: Turkic languages exhibit extensive and diverse etymological relationships among lexical items. These relationships make the Turkic languages promising for exploring automated translation lexicon induction by leveraging cognate and other etymological information. However, due to the extent and diversity of the types of relationships between words, it is not clear how to annotate such information. In this paper, we present a methodology for annotating cognates and etymological origin in Turkic languages. Our method strives to balance the amount of research effort the annotator expends with the utility of the annotations for supporting research on improving automated translation lexicon induction.


title: Optical Tractor Beam with Chiral Light

abstract: We suggest a novel mechanism to induce the motion of a chiral material body towards an optical source. Our solution is based on the interference between a chiral light beam and its reflection on an opaque mirror. Surprisingly, it is theoretically shown that the electromagnetic response of the material may be tailored in such a way that independent of the specific body location with the respect to the mirror, it is always pushed upstream against the photon flow associated with the incoming wave. Moreover, it is proven that by controlling the handedness of the incoming light it may be possible to harness the sign of the optical force, switching from a pulling force to a pushing force.


title: Low-energy capture of asteroids onto KAM tori

abstract: We present a new method for engineering the artificial capture of asteroids. Based on theories of the chaos-assisted capture of natural satellites of the giant planets, we show how an unbound asteroid that passes close to a regular region of phase space can be easily moved onto the nearby KAM tori and essentially permanently captured with the Earth's Hill sphere without closing the zero velocity curves. The method has the advantages of a relatively low delta-v requirement and no need for control strategies. An illustration of the method is given for an example asteroid trajectory, demonstrating that it is a viable strategy for the final capture stage of asteroids in the Earth's neighbourhood.


title: The cavity method for analysis of large-scale penalized regression

abstract: Penalized regression methods aim to retrieve reliable predictors among a large set of putative ones from a limited amount of measurements. In particular, penalized regression with singular penalty functions is important for sparse reconstruction algorithms. For large-scale problems, these algorithms exhibit sharp phase transition boundaries where sparse retrieval breaks down. Large optimization problems associated with sparse reconstruction have been analyzed in the literature by setting up corresponding statistical mechanical models at a finite temperature. Using replica method for mean field approximation, and subsequently taking a zero temperature limit, this approach reproduces the algorithmic phase transition boundaries. Unfortunately, the replica trick and the non-trivial zero temperature limit obscure the underlying reasons for the failure of a sparse reconstruction algorithm, and of penalized regression methods, in general. In this paper, we employ the ``cavity method'' to give an alternative derivation of the mean field equations, working directly in the zero-temperature limit. This derivation provides insight into the origin of the different terms in the self-consistency conditions. The cavity method naturally involves a quantity, the average local susceptibility, whose behavior distinguishes different phases in this system. This susceptibility can be generalized for analysis of a broader class of sparse reconstruction algorithms.

categories: cs.IT cond-mat.dis-nn cond-mat.stat-mech math.IT


title: VUV-synchrotron absorption studies of N$_2$ and CO at 900 K

abstract: Photoabsorption spectra of N$_2$ and CO were recorded at 900K, using the vacuum-ultraviolet Fourier-transform spectrometer at the DESIRS beamline of synchrotron SOLEIL. These high-temperature and high-resolution measurements allow for precise determination of line wavelengths, oscillator strengths, and predissociative line broadening of highly-excited rotational states with $J$ up to about 50, and also vibrational hot bands. In CO, the perturbation of $A\,{}^1\Pi-X\,{}^1\Sigma^+$ vibrational bands $(0,0)$ and $(1,0)$ were studied, as well as the transitions to perturbing optically-forbidden states $e\,{}^3\Sigma^-$, $d\,{}^3\Delta$, $D\,{}^1\Delta$ and $a'\,{}^3\Sigma^+$. In N$_2$, we observed line shifts and broadening in several $b\,{}^1\Pi_u-X\,{}^1\Sigma^+_g$ bands due to unobserved forbidden states of ${}^3\Pi_u$ symmetry. The observed state interactions are deperturbed and, for N$_2$, used to validate a coupled-channels model of the interacting electronic states. This data is appropriate for use in astrophysical or (exo-)planetary atmospheric applications where high temperatures are important and in future spectroscopic models of these molecules.


title: Forward Delay-based Packet Scheduling Algorithm for Multipath TCP

abstract: Multipath TCP (MPTCP) is a transport layer protocol that allows network devices to transfer data over multiple concurrent paths, and hence, utilizes the network resources more effectively than does the traditional single-path TCP. However, as a reliable protocol, MPTCP still needs to deliver data packets (to the upper application) at the receiver in the same order they are transmitted at the sender. The out-of-order packet problem becomes more severe for MPTCP due to the heterogeneous nature of delay and bandwidth of each path. In this paper, we propose the forward-delay-based packet scheduling (FDPS) algorithm for MPTCP to address that problem. The main idea is that the sender dispatches packets via concurrent paths according to their estimated forward delay and throughput differences. Via simulations with various network conditions, the results show that our algorithm significantly maintains in-order arrival packets at the receiver compared with several previous algorithms.


title: The lower bound for the modulus of the derivatives and Jacobian of harmonic injective mappings

abstract: We give the lower bound for the modulus of the radial derivatives and Jacobian of harmonic injective mappings from the unit ball onto convex domain in plane and space. As an application we show co-Lipschitz property of some classes of qch mappings. We also review related results in planar case using some novelty.


title: Wave Function Collapse and the No-Superluminal-Signaling Principle

abstract: The assumption that wave function collapse is a real occurrence has very interesting consequences - both experimental and theoretical. Besides predicting observable deviations from linear evolution, it implies that these deviations must originate in nondeterministic effects at the elementary level in order to prevent superluminal signaling, as demonstrated by Gisin. This lack of determinism implies that information cannot be instantiated in a reproducible form in isolated microsystems (as illustrated by the No-cloning theorem). By stipulating that information is a reproducible and referential property of physical systems, one can formulate the no-signaling principle in strictly physical terms as a prohibition of the acquisition of information about spacelike-separated occurrences. This formulation provides a new perspective on the relationship between relativity and spacetime structure, and it imposes tight constraints on the way in which collapse effects are induced. These constraints indicate that wave function collapse results from (presumably small) nondeterministic deviations from linear evolution associated with nonlocally entangling interactions. This hypothesis can be formalized in a stochastic collapse equation and used to assess the feasibility of testing for collapse effects.


title: Expansion, Random Walks and Sieving in $SL_2 (\mathbb{F}_p [t])$

abstract: We construct new examples of expander Cayley graphs of finite groups, arising as congruence quotients of non-elementary subgroups of $SL_2 (\mathbb{F}_p [t])$ modulo certain square-free ideals. We describe some applications of our results to simple random walks on such subgroups, specifically giving bounds on the rate of escape of such walks from algebraic subvarieties, the set of squares and the set of elements with reducible characteristic polynomial in $SL_2 (\mathbb{F}_p [t])$.


title: Supremum distribution of Bessel process of drifting Brownian motion

abstract: Let (B^{(1)}_t ;B^{(2)}_t ;B^{(3)}_t + \mu t) be a three-dimensional Brownian motion with drift \mu, starting at the origin. Then X_t = ||(B^{(1)}_t ;B^{(2)}_t ;B^{(3)}_t +\mu t)||, its distance from the starting point, is a diffusion with many applications. We investigate the distribution of the supremum of (X_t), give an infinite-series formula for its density and an exact estimate by elementary functions.


title: Perturbative N=2 supersymmetric quantum mechanics and L-theory with complex coefficients

abstract: We construct L-theory with complex coefficients from the geometry of 1|2-dimensional perturbative mechanics. Methods of perturbative quantization lead to wrong-way maps that we identify with those coming from the MSO-orientation of L-theory tensored with the complex numbers.


title: No Return to Classical Reality

abstract: At a fundamental level, the classical picture of the world is dead, and has been dead now for almost a century. Pinning down exactly which quantum phenomena are responsible for this has proved to be a tricky and controversial question, but a lot of progress has been made in the past few decades. We now have a range of precise statements showing that whatever the ultimate laws of Nature are, they cannot be classical. In this article, we review results on the fundamental phenomena of quantum theory that cannot be understood in classical terms. We proceed by first granting quite a broad notion of classicality, describe a range of quantum phenomena (such as randomness, discreteness, the indistinguishability of states, measurement-uncertainty, measurement-disturbance, complementarity, noncommutativity, interference, the no-cloning theorem, and the collapse of the wave-packet) that do fall under its liberal scope, and then finally describe some aspects of quantum physics that can never admit a classical understanding -- the intrinsically quantum mechanical aspects of Nature. The most famous of these is Bell's theorem, but we also review two more recent results in this area. Firstly, Hardy's theorem shows that even a finite dimensional quantum system must contain an infinite amount of information, and secondly, the Pusey--Barrett--Rudolph theorem shows that the wave-function must be an objective property of an individual quantum system. Besides being of foundational interest, results of this sort now find surprising practical applications in areas such as quantum information science and the simulation of quantum systems.


title: Continuous versus discrete structures II -- Discrete Hamiltonian systems and Helmholtz conditions

abstract: We define discrete Hamiltonian systems in the framework of discrete embeddings. An explicit comparison with previous attempts is given. We then solve the discrete Helmholtz's inverse problem for the discrete calculus of variation in the Hamiltonian setting. Several applications are discussed.


title: Testing the Standard Model under the weight of heavy flavors

abstract: I review recently completed (since Lattice 2013) and ongoing lattice calculations in charm and bottom flavor physics. A comparison of the precision of lattice and experiment is made using both current experimental results and projected experimental precision in 2020. The combination of experiment and theory reveals several tensions between nature and the Standard Model. These tensions are reviewed in light of recent lattice results.


title: Adaptive behavior can produce maladaptive anxiety due to individual differences in experience

abstract: Normal anxiety is considered an adaptive response to the possible presence of danger, but is susceptible to dysregulation. Anxiety disorders are prevalent at high frequency in contemporary human societies, yet impose substantial disability upon their sufferers. This raises a puzzle: why has evolution left us vulnerable to anxiety disorders? We develop a signal detection model in which individuals must learn how to calibrate their anxiety responses: they need to learn which cues indicate danger in the environment. We derive the optimal strategy for doing so, and find that individuals face an inevitable exploration-exploitation tradeoff between obtaining a better estimate of the level of risk on one hand, and maximizing current payoffs on the other. Because of this tradeoff, a subset of the population can become trapped in a state of self-perpetuating over-sensitivity to threatening stimuli, even when individuals learn optimally. This phenomenon arises because when individuals become too cautious, they stop sampling the environment and fail to correct their misperceptions, whereas when individuals become too careless they continue to sample the environment and soon discover their mistakes. Thus, over-sensitivity to threats becomes common whereas under-sensitivity becomes rare. We suggest that this process may be involved in the development of excessive anxiety in humans.


title: Nanoconfined catalytic {\AA}ngstr\"om-size motors

abstract: Self-propelled chemically powered synthetic micron and nano-scale motors are being intensively studied because of the wide range of potential applications that exploit their directed motion. This paper considers even smaller {\AA}ngstr\"om-size synthetic motors. Such very small motors in bulk solution display effects arising from their self-propulsion. Recent experiments have shown that small-molecule catalysts and single enzyme molecules exhibit properties that have been attributed to their chemical activity. Molecular dynamics is used to investigate the properties of very small {\AA}ngstr\"om-size synthetic chemically powered sphere-dimer motors in a simple atomic-like solvent confined between walls separated by distances of tens of nanometers. Evidence for strong structural ordering of the motors between the walls, which reflects the finite size of solvent molecules and depends on solvent depletion forces, is provided. Dynamical properties, such as average motor velocity, orientational relaxation, and mean square displacement, are anisotropic and depend on the distance from the walls. This research provides information needed for potential applications that use molecular-scale motors in the complex confined geometries encountered in biology and the laboratory.


title: Detecting vortices in superconductors: Extracting one-dimensional topological singularities from a discretized complex scalar field

abstract: In type-II superconductors, the dynamics of superconducting vortices determine their transport properties. In the Ginzburg-Landau theory, vortices correspond to topological defects in the complex order parameter. Extracting their precise positions and motion from discretized numerical simulation data is an important, but challenging task. In the past, vortices have mostly been detected by analyzing the magnitude of the complex scalar field representing the order parameter and visualized by corresponding contour plots and isosurfaces. However, these methods, primarily used for small-scale simulations, blur the fine details of the vortices, scale poorly to large-scale simulations, and do not easily enable isolating and tracking individual vortices. Here we present a method for exactly finding the vortex core lines from a complex order parameter field. With this method, vortices can be easily described at a resolution even finer than the mesh itself. The precise determination of the vortex cores allows the interplay of the vortices inside a model superconductor to be visualized in higher resolution than has previously been possible. By representing the field as the set of vortices, this method also massively reduces the data footprint of the simulations and provides the data structures for further analysis and feature tracking.

categories: cond-mat.supr-con cond-mat.mes-hall math.NA physics.comp-ph


title: Compressive Sensing with Redundant Dictionaries and Structured Measurements

abstract: Consider the problem of recovering an unknown signal from undersampled measurements, given the knowledge that the signal has a sparse representation in a specified dictionary $D$. This problem is now understood to be well-posed and efficiently solvable under suitable assumptions on the measurements and dictionary, if the number of measurements scales roughly with the sparsity level. One sufficient condition for such is the $D$-restricted isometry property ($D$-RIP), which asks that the sampling matrix approximately preserve the norm of all signals which are sufficiently sparse in $D$. While many classes of random matrices are known to satisfy such conditions, such matrices are not representative of the structural constraints imposed by practical sensing systems. We close this gap in the theory by demonstrating that one can subsample a fixed orthogonal matrix in such a way that the $D$-RIP will hold, provided this basis is sufficiently incoherent with the sparsifying dictionary $D$. We also extend this analysis to allow for weighted sparse expansions. Consequently, we arrive at compressive sensing recovery guarantees for structured measurements and redundant dictionaries, opening the door to a wide array of practical applications.


title: Neural Implementation of Probabilistic Models of Cognition

abstract: Bayesian models of cognition hypothesize that human brains make sense of data by representing probability distributions and applying Bayes' rule to find the best explanation for available data. Understanding the neural mechanisms underlying probabilistic models remains important because Bayesian models provide a computational framework, rather than specifying mechanistic processes. Here, we propose a deterministic neural-network model which estimates and represents probability distributions from observable events --- a phenomenon related to the concept of probability matching. Our model learns to represent probabilities without receiving any representation of them from the external world, but rather by experiencing the occurrence patterns of individual events. Our neural implementation of probability matching is paired with a neural module applying Bayes' rule, forming a comprehensive neural scheme to simulate human Bayesian learning and inference. Our model also provides novel explanations of base-rate neglect, a notable deviation from Bayes.


title: Towards Deep Semantic Analysis Of Hashtags

abstract: Hashtags are semantico-syntactic constructs used across various social networking and microblogging platforms to enable users to start a topic specific discussion or classify a post into a desired category. Segmenting and linking the entities present within the hashtags could therefore help in better understanding and extraction of information shared across the social media. However, due to lack of space delimiters in the hashtags (e.g #nsavssnowden), the segmentation of hashtags into constituent entities ("NSA" and "Edward Snowden" in this case) is not a trivial task. Most of the current state-of-the-art social media analytics systems like Sentiment Analysis and Entity Linking tend to either ignore hashtags, or treat them as a single word. In this paper, we present a context aware approach to segment and link entities in the hashtags to a knowledge base (KB) entry, based on the context within the tweet. Our approach segments and links the entities in hashtags such that the coherence between hashtag semantics and the tweet is maximized. To the best of our knowledge, no existing study addresses the issue of linking entities in hashtags for extracting semantic information. We evaluate our method on two different datasets, and demonstrate the effectiveness of our technique in improving the overall entity linking in tweets via additional semantic information provided by segmenting and linking entities in a hashtag.


title: Measuring free energy in parallel tempering Monte Carlo

abstract: An efficient approach of measuring the absolute free energy in parallel tempering Monte Carlo using the exponential averaging method is discussed and the results are compared with those of population annealing Monte Carlo using the three-dimensional Edwards-Anderson Ising spin glass model as benchmark tests. Numerical results show that parallel tempering, even though uses a much less number of temperatures than population annealing, can nevertheless equally efficiently measure the absolute free energy by simulating each temperature for longer times.


title: Star formation in the first galaxies - III. Formation, evolution, and characteristics of the first stellar cluster

abstract: We simulate the formation of a low metallicity (0.01 Zsun) stellar cluster in a dwarf galaxy at redshift z~14. Beginning with cosmological initial conditions, the simulation utilizes adaptive mesh refinement and sink particles to follow the collapse and evolution of gas past the opacity limit for fragmentation, thus resolving the formation of individual protostellar cores. A time- and location-dependent protostellar radiation field, which heats the gas by absorption on dust, is computed by integration of protostellar evolutionary tracks with the MESA code. The simulation also includes a robust non-equilibrium chemical network that self-consistently treats gas thermodynamics and dust-gas coupling. The system is evolved for 18 kyr after the first protostellar source has formed. In this time span, 30 sink particles representing protostellar cores form with a total mass of 81 Msun. Their masses range from ~0.1 Msun to 14.4 Msun with a median mass ~0.5-1 Msun. Massive protostars grow by competitive accretion while lower-mass protostars are stunted in growth by close encounters and many-body ejections. In the regime explored here, the characteristic mass scale is determined by the temperature floor set by the cosmic microwave background and by the onset of efficient dust-gas coupling. It seems unlikely that host galaxies of the first bursts of metal-enriched star formation will be detectable with the James Webb Space Telescope or other next-generation infrared observatories. Instead, the most promising access route to the dawn of cosmic star formation may lie in the scrutiny of metal-poor, ancient stellar populations in the Galactic neighborhood. The observable targets that correspond to the system simulated here are ultra-faint dwarf satellite galaxies such as Bootes II, Segue I and II, and Willman I.


title: The Relativistic Green's function model and charged-current inclusive neutrino-nucleus scattering at T2K kinematics

abstract: We compare the results of the relativistic Green's function model with the experimental data of the charged-current inclusive differential neutrino-nucleus cross sections published by the T2K Collaboration. The model, which is able to describe both MINER$\nu$A and MiniBooNE charged-current quasielastic scattering data, underpredicts the inclusive T2K cross sections.


title: Quantifying Prosodic Variability in Middle English Alliterative Poetry

abstract: Interest in the mathematical structure of poetry dates back to at least the 19th century: after retiring from his mathematics position, J. J. Sylvester wrote a book on prosody called $\textit{The Laws of Verse}$. Today there is interest in the computer analysis of poems, and this paper discusses how a statistical approach can be applied to this task. Starting with the definition of what Middle English alliteration is, $\textit{Sir Gawain and the Green Knight}$ and William Langland's $\textit{Piers Plowman}$ are used to illustrate the methodology. Theory first developed for analyzing data from a Riemannian manifold turns out to be applicable to strings allowing one to compute a generalized mean and variance for textual data, which is applied to the poems above. The ratio of these two variances produces the analogue of the F test, and resampling allows p-values to be estimated. Consequently, this methodology provides a way to compare prosodic variability between two texts.


title: Creating, Automating, and Assessing Online Homework in Introductory Statistics and Mathematics Classes

abstract: Although textbook publishers offer course management systems, they do so to promote brand loyalty, and while an open source tool such as WeBWorK is promising, it requires administrative and IT buy-in. So supported in part by a College Access Challenge Grant from the Department of Education, we collaborated with other instructors to create online homework sets for three classes: Elementary Algebra, Intermediate Algebra, and Statistics for Behavioral Sciences I. After experimentation, some of these question pools are now created by Mathematica programs that can generate data sets from specified distributions, generate random polynomials that factor in a given way, create image files of histograms, scatterplots, and so forth. These programs produce files that can be read by the software package, Respondus, which then uploads the questions into Blackboard Learn, the course management system used by the Connecticut State University system. Finally, we summarize five classes worth of student performance data along with lessons learned while working on this project.


title: Sorting photon wave packets using temporal-mode interferometry based on multiple-stage quantum frequency conversion

abstract: All classical and quantum technologies that encode in and retrieve information from optical fields rely on the ability to selectively manipulate orthogonal field modes of light. Such manipulation can be achieved with high selectivity for polarization modes and transverse-spatial modes. For the time-frequency degree of freedom, this could efficiently be achieved for a limited choice of approximately orthogonal modes, i.e. non-overlapping bins in time or frequency. We recently proposed a method that surmounts the selectivity barrier for sorting arbitrary orthogonal temporal modes [Opt. Lett. {\bf 39}, 2924 (2014)] using cascaded interferometric quantum frequency conversion in nonlinear optical media. We call this method temporal-mode interferometry, as it has a close resemblance to the well-known separated-fields atomic interferometry method introduced by Ramsey. The method has important implications for quantum memories, quantum dense coding, quantum teleportation, and quantum key distribution. Here we explore the inner workings of the method in detail, and extend it to multiple stages with a concurrent asymptotic convergence of temporal-mode selectivity to unity. We also complete our analysis of pump-chirp compensation to counter pump-induced nonlinear phase-modulation in four-wave mixing implementations.


title: SgpDec: Cascade (De)Compositions of Finite Transformation Semigroups and Permutation Groups

abstract: We describe how the SgpDec computer algebra package can be used for composing and decomposing permutation groups and transformation semigroups hierarchically by directly constructing substructures of wreath products, the so called cascade products.


title: Expanding the Vocabulary of Multitouch Input using Magnetic Fingerprints

abstract: We present magnetic fingerprints; an input technique for mobile touchscreen devices that uses a small magnet attached to a user's fingernail in order to differentiate between a normal touch and a magnetic touch. The polarity of the magnet can be used to create different magnetic fingerprints where this technique takes advantage of the rich vocabulary offered by the use of multitouch input. User studies investigate the accuracy of magnetic fingerprint recognition in relation to magnet size, number of magnetic fingerprints used; and size of the touchscreen. Studies found our technique to be limited to using up to two fingerprints non-simultaneously, while achieving a high classification accuracy (95%) but it nearly triples the number of distinguishable multi touch events. Potential useful applications of this technique are presented.


title: On the use of shot noise for photon counting

abstract: Lieu et al. (2015) have recently claimed that it is possible to substantially improve the sensitivity of radio astronomical observations. In essence, their proposal is to make use of the intensity of the photon shot noise as a measure of the photon arrival rate. Lieu et al. (2015) provide a detailed quantum-mechanical calculation of a proposed measurement scheme that uses two detectors and conclude that this scheme avoids the sensitivity degradation that is associated with photon bunching. If correct, this result could have a profound impact on radio astronomy. Here I present a detailed analysis of the sensitivity attainable using shot-noise measurement schemes that use either one or two detectors, and demonstrate that neither scheme can avoid the photon bunching penalty. I perform both semiclassical and fully quantum calculations of the sensitivity, obtaining consistent results, and provide a formal proof of the equivalence of these two approaches. These direct calculations are furthermore shown to be consistent with an indirect argument based on a correlation method that establishes an independent limit to the sensitivity of shot-noise measurement schemes. Collectively, these results conclusively demonstrate that the photon bunching sensitivity penalty applies to shot noise measurement schemes just as it does to ordinary photon counting, in contradiction to the fundamental claim made by Lieu et al. (2015). The source of this contradiction is traced to a logical fallacy in their argument.


title: Quintessence and phantom emerging from the split-complex field and the split-quaternion field

abstract: Motivated by the mathematic theory of split-complex numbers (or hyperbolic numbers, also perplex numbers) and the split-quaternion numbers (or coquaternion numbers), we define the notion of split-complex scalar field and the split-quaternion scalar field. Then we explore the cosmic evolution of these scalar fields in the background of spatially flat Friedmann-Robertson-Walker Universe. We find that both the quintessence field and the phantom field could naturally emerge in these scalar fields. Introducing the metric of field space, these theories fall into a subclass of the multi-field theories which have been extensively studied in inflationary cosmology.


title: Regularized Principal Component Analysis for Spatial Data

abstract: In many atmospheric and earth sciences, it is of interest to identify dominant spatial patterns of variation based on data observed at $p$ locations and $n$ time points with the possibility that $p>n$. While principal component analysis (PCA) is commonly applied to find the dominant patterns, the eigenimages produced from PCA may exhibit patterns that are too noisy to be physically meaningful when $p$ is large relative to $n$. To obtain more precise estimates of eigenimages, we propose a regularization approach incorporating smoothness and sparseness of eigenimages, while accounting for their orthogonality. Our method allows data taken at irregularly spaced or sparse locations. In addition, the resulting optimization problem can be solved using the alternating direction method of multipliers, which is easy to implement, and applicable to a large spatial dataset. Furthermore, the estimated eigenfunctions provide a natural basis for representing the underlying spatial process in a spatial random-effects model, from which spatial covariance function estimation and spatial prediction can be efficiently performed using a regularized fixed-rank kriging method. Finally, the effectiveness of the proposed method is demonstrated by several numerical examples


title: Independence of Satellites of Torus Knots in the Smooth Concordance Group

abstract: The main goal of this article is to obtain a condition under which an infinite collection $\mathscr{F}$ of satellite knots (with companion a positive torus knot and pattern similar to the Whitehead link) freely generates a subgroup of infinite rank in the smooth concordance group. This goal is attained by examining both the instanton moduli space over a 4-manifold with tubular ends and the corresponding Chern-Simons invariant of the adequate 3-dimensional portion of the 4-manifold. More specifically, the result is derived from Furuta's criterion for the independence of Seifert fibred homology spheres in the homology cobordism group of oriented homology 3-spheres. Indeed, we first associate to $\mathscr{F}$ the corresponding collection of 2-fold covers of the 3-sphere branched over the elements of $\mathscr{F}$ and then introduce definite cobordisms from the aforementioned covers of the satellites to a number of Seifert fibered homology spheres. This allows us to apply Furuta's criterion and thus obtain a condition that guarantees the independence of the family $\mathscr{F}$ in the smooth concordance group.


title: The Interstellar Line of Sight to the Interacting Galaxy NGC 5195

abstract: We present moderately-high resolution echelle observations of the nucleus of NGC 5195, the line of sight to which samples intervening interstellar material associated with the outer spiral arm of M51. Our spectra reveal absorption from interstellar Na I, K I, Ca II, and CH+, and from a number of diffuse interstellar bands (DIBs), at a velocity close to that exhibited by H I 21 cm emission from M51 at the position of NGC 5195. The H I column density implied by the equivalent width of the 5780.5 DIB, based on the relationship between W(5780.5) and N(H I) derived for sight lines in the local Galactic interstellar medium, is consistent with the column density obtained from the integrated H I emission. The H2 column density predicted from the observed column density of K I, using the Galactic relationship between N(K I) and N(H2), is comparable to N(H I), suggesting a high molecular fraction (~0.65) for the M51 gas toward NGC 5195. The DIBs toward NGC 5195 are, on average, ~40% weaker than would be expected based on the K I column density, a further indication that the gas in this direction has a high molecular content. The M51 material is characterized also by a high N(Na I)/N(Ca II) ratio (>11), indicative of a high degree of Ca depletion, and a high W(5797.1)/W(5780.5) ratio (~1.6), suggestive of either a very weak ambient radiation field or a significantly shielded environment. A high N(CH+)/N(CH) ratio (>2.3) for the M51 material toward NGC 5195 may be the result of enhanced turbulence due to interactions between M51 and its companion.


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

abstract: The single-particle density is the most basic quantity that can be calculated from a given many-body wave function. It provides the probability to find a particle at a given position when the average over many realizations of an experiment is taken. However, the outcome of single experimental shots of ultracold atom experiments is determined by the $N$-particle probability density. This difference can lead to surprising results. For example, independent Bose-Einstein condensates (BECs) with definite particle numbers form interference fringes even though no fringes would be expected based on the single-particle density [1-4]. By drawing random deviates from the $N$-particle probability density single experimental shots can be simulated from first principles [1, 3, 5]. However, obtaining expressions for the $N$-particle probability density of realistic time-dependent many-body systems has so far been elusive. Here, we show how single experimental shots of general ultracold bosonic systems can be simulated based on numerical solutions of the many-body Schr\"odinger equation. We show how full counting distributions of observables involving any number of particles can be obtained and how correlation functions of any order can be evaluated. As examples we show the appearance of interference fringes in interacting independent BECs, fluctuations in the collisions of strongly attractive BECs, the appearance of randomly fluctuating vortices in rotating systems and the center of mass fluctuations of attractive BECs in a harmonic trap. The method described is broadly applicable to bosonic many-body systems whose phenomenology is driven by information beyond what is typically available in low-order correlation functions.


title: X-ray analysis of the proper motion and pulsar wind nebula for PSR J1741-2054

abstract: We obtained six observations of PSR J1741-2054 using the $Chandra$ ACIS-S detector totaling $\sim$300 ks. By registering this new epoch of observations to an archival observation taken 3.2 years earlier using X-ray point sources in the field of view, we have measured the pulsar proper motion at $\mu =109 \pm 10 {\rm mas\ yr}^{-1}$ in a direction consistent with the symmetry axis of the observed H$\alpha$ nebula. We investigated the inferred past trajectory of the pulsar but find no compelling association with OB associations in which the progenitor may have originated. We confirm previous measurements of the pulsar spectrum as an absorbed power law with photon index $\Gamma$=2.68$\pm$0.04, plus a blackbody with an emission radius of (4.5$^{+3.2}_{-2.5})d_{0.38}$ km, for a DM-estimated distance of $0.38d_{0.38}$ kpc and a temperature of $61.7\pm3.0$ eV. Emission from the compact nebula is well described by an absorbed power law model with a photon index of $\Gamma$ = 1.67$\pm$0.06, while the diffuse emission seen as a trail extending northeast of the pulsar shows no evidence of synchrotron cooling. We also applied image deconvolution techniques to search for small-scale structures in the immediate vicinity of the pulsar, but found no conclusive evidence for such structures.


title: Suppression of 3D-ordering by defects in the $S=1/2$ frustrated chain magnet LiCuVO$_4$

abstract: We report on a heat capacity study of high quality single crystal samples of \lcvo\ -- a frustrated spin $S=1/2$ chain system -- in magnetic field amounting to 3/4 of the saturation field. At low fields up to about 7~T, a linear temperature dependence of the specific heat, $C_p\propto T$, resulting from 1D magnetic correlations in the spin chains is followed upon cooling by a sharp lambda anomaly of the transition into a 3D helical phase with $C_p\propto T^3$ behavior at low temperature. The transition from a spin liquid into a spin-modulated (SM) phase at higher fields occurs via a hump-like anomaly which, as the temperature decreases further turns into a $C_p\propto T^2$ law distinctive for a quasi-2D system. We suggest an explanation for how nonmagnetic defects in the Cu$^{2+}$ chains can suppress 3D long-range ordering in the SM phase and leave it undisturbed in a helical phase.


title: Using Riemannian geometry for SSVEP-based Brain Computer Interface

abstract: Riemannian geometry has been applied to Brain Computer Interface (BCI) for brain signals classification yielding promising results. Studying electroencephalographic (EEG) signals from their associated covariance matrices allows a mitigation of common sources of variability (electronic, electrical, biological) by constructing a representation which is invariant to these perturbations. While working in Euclidean space with covariance matrices is known to be error-prone, one might take advantage of algorithmic advances in information geometry and matrix manifold to implement methods for Symmetric Positive-Definite (SPD) matrices. This paper proposes a comprehensive review of the actual tools of information geometry and how they could be applied on covariance matrices of EEG. In practice, covariance matrices should be estimated, thus a thorough study of all estimators is conducted on real EEG dataset. As a main contribution, this paper proposes an online implementation of a classifier in the Riemannian space and its subsequent assessment in Steady-State Visually Evoked Potential (SSVEP) experimentations.


title: Mixed Eigenvalues of {\LARGE $\pmb p$\,}-Laplacian

abstract: The mixed principal eigenvalue of $p\,$-Laplacian (equivalently, the optimal constant of weighted Hardy inequality in $L^p$ space) is studied in this paper. Several variational formulas for the eigenvalue are presented. As applications of the formulas, a criterion for the positivity of the eigenvalue is obtained. Furthermore, an approximating procedure and some explicit estimates are presented case by case. An example is included to illustrate the power of the results of the paper.


title: Thermoelectric properties of semiconductor nanowire networks

abstract: To examine thermoelectric (TE) properties of a semiconductor nanowire (NW) network, we propose a theoretical approach mapping the TE network on a two-port network. In contrast to a conventional single-port (i.e., resistor) network model, our model allows for large scale calculations showing convergence of TE figure of merit, $ZT$, with an increasing number of junctions. Using this model, numerical simulations are performed for the Bi$_2$Te$_3$ branched nanowire (BNW) and Cayley tree NW (CTNW) network. We find that the phonon scattering at the network junctions plays a dominant role in enhancing the network $ZT$. Specifically, disordered BNW and CTNW demonstrate an order of magnitude higher ZT enhancement compared to their ordered counterparts. Formation of preferential TE pathways in CTNW makes the network effectively behave as its BNW counterpart. We provide formalism for simulating large scale nanowire networks hinged upon experimentally measurable TE parameters of a single T-junction.


title: The optimal constant in Hardy-type inequalities

abstract: To estimate the optimal constant in Hardy-type inequalities, some variational formulas and approximating procedures are introduced. The known basic estimates are improved considerably. The results are illustrated by typical examples. It is shown that the sharp factor is meaningful for each finite interval and a classical sharp model is re-examined.


title: Long-lived Chaotic Orbital Evolution of Exoplanets in Mean Motion Resonances with Mutual Inclinations

abstract: We present N-body simulations of resonant planets with inclined orbits that show chaotically evolving eccentricities and inclinations that can persist for at least 10 Gyr. A wide range of behavior is possible, from fast, low amplitude variations to systems in which eccentricities reach 0.9999 and inclinations 179.9 degrees. While the orbital elements evolve chaotically, at least one resonant argument always librates. We show that the HD 73526, HD 45364 and HD 60532 systems may be in chaotically-evolving resonances. Chaotic evolution is apparent in the 2:1, 3:1 and 3:2 resonances, and for planetary masses from lunar- to Jupiter-mass. In some cases, orbital disruption occurs after several Gyr, implying the mechanism is not rigorously stable, just long-lived relative to the main sequence lifetimes of solar-type stars. Planet-planet scattering appears to yield planets in inclined resonances that evolve chaotically in about 0.5% of cases. These results suggest that 1) approximate methods for identifying unstable orbital architectures may have limited applicability, 2) the observed close-in exoplanets may be produced during the high eccentricity phases induced by inclined resonances, 3) those exoplanets' orbital planes may be misaligned with the host star's spin axis, 4) systems with resonances may be systematically younger than those without, 5) the distribution of period ratios of adjacent planets detected via transit may be skewed due to inclined resonances, and 6) potentially habitable planets in resonances may have dramatically different climatic evolution than the Earth. The GAIA spacecraft is capable of discovering giant planets in these types of orbits.


title: Efficiency fluctuations in quantum thermoelectric devices

abstract: We present a method, based on characterizing efficiency fluctuations, to asses the performance of nanoscale thermoelectric junctions. This method accounts for effects typically arising in small junctions, namely, stochasticity in the junction's performance, quantum effects, and nonequilibrium features preventing a linear response analysis. It is based on a nonequilibrium Green's function (NEGF) approach, which we use to derive the full counting statistics (FCS) for heat and work, and which in turn allows us to calculate the statistical properties of efficiency fluctuations. We simulate the latter for a variety of simple models where our method is exact. By analyzing the discrepancies with the semi-classical prediction of a quantum master equation (QME) approach, we emphasize the quantum nature of efficiency fluctuations for realistic junction parameters. We finally propose an approximate Gaussian method to express efficiency fluctuations in terms of nonequilibrium currents and noises which are experimentally measurable in molecular junctions.


title: Criteria for Discrete Spectrum of 1D Operators

abstract: For discrete spectrum of 1D second-order differential/difference operators (with or without potential (killing), with the maximal/minimal domain), a pair of unified dual criteria are presented in terms of two explicit measures and the harmonic function of the operators. Interes\-tingly, these criteria can be read out from the ones for the exponential convergence of four types of stability studied earlier, simply replacing the `finite supremum' by `vanishing at infinity'. Except a dual technique, the main tool used here is a transform in terms of the harmonic function, to which two new practical algorithms are introduced in the discrete context and two successive approximation schemes are reviewed in the continuous context. All of them are illustrated by examples. The main body of the paper is devoted to the hard part of the story, the easier part but powerful one is delayed to the end of the paper.


title: Quotient singularities, eta invariants, and self-dual metrics

abstract: There are three main components to this article: (i) A formula for the eta invariant of the signature complex for any finite subgroup of ${\rm{SO}}(4)$ acting freely on $S^3$ is given. An application of this is a non-existence result for Ricci-flat ALE metrics on certain spaces. (ii) A formula for the orbifold correction term that arises in the index of the self-dual deformation complex is proved for all finite subgroups of ${\rm{SO}}(4)$ which act freely on $S^3$. Some applications of this formula to the realm of self-dual and scalar-flat K\"ahler metrics are also discussed. (iii) Two infinite families of scalar-flat anti-self-dual ALE spaces with groups at infinity not contained in ${\rm{U}}(2)$ are constructed. Using these spaces, new examples of self-dual metrics on $n \# \mathbb{CP}^2$ are obtained for $n \geq 3$.


title: Successive Cancellation Decoding of Polar Codes using Stochastic Computing

abstract: Polar codes have emerged as the most favorable channel codes for their unique capacity-achieving property. To date, numerous works have been reported for efficient design of polar codes decoder. However, these prior efforts focused on design of polar decoders via deterministic computation, while the behavior of stochastic polar decoder, which can have potential advantages such as low complexity and strong error-resilience, has not been studied in existing literatures. This paper, for the first time, investigates polar decoding using stochastic logic. Specifically, the commonly-used successive cancellation (SC) algorithm is reformulated into the stochastic form. Several methods that can potentially improve decoding performance are discussed and analyzed. Simulation results show that a stochastic SC decoder can achieve similar error-correcting performance as its deterministic counterpart. This work can pave the way for future hardware design of stochastic polar codes decoders.


title: Footnotes to a paper by Domokos

abstract: The paper of Domokos is "Invariants of quivers and wreath products," and as in that paper we use Regev's double centralizer theorem for wreath products to study trace identities and invariant theory. In addition to extending to other algebras and groups, we also compute Poincare series and prove embedding theorems.


title: Blocks of defect of p-solvable groups

abstract: Let $p$ be a prime such that $p \geq 5$. Let $G$ be a finite $p$-solvable group and let $p^a$ be the largest power of $p$ dividing $\chi(1)$ for an irreducible character $\chi$ of $G$, we show that $|G:F(G)|_p \leq p^{5.5a}$. Let $G$ be a finite $p$-solvable group with trivial maximal normal solvable subgroup and we denote $|G|_p=p^n$, then $G$ contains a block of defect less than or equal to $\lfloor \frac {2n} {3} \rfloor$.


title: Valid plane trees: Combinatorial models for RNA secondary structures with Watson-Crick base pairs

abstract: The combinatorics of RNA plays a central role in biology. Mathematical biologists have several commonly-used models for RNA: words in a fixed alphabet (representing the primary sequence of nucleotides) and plane trees (representing the secondary structure, or folding of the RNA sequence). This paper considers an augmented version of the standard model of plane trees, one that incorporates some observed constraints on how the folding can occur. In particular we assume the alphabet consists of complementary pairs, for instance the Watson-Crick pairs A-U and C-G of RNA. Given a word in the alphabet, a valid plane tree is a tree for which, when the word is folded around the tree, each edge matches two complementary letters. Consider the graph whose vertices are valid plane trees for a fixed word and whose edges are given by Condon, Heitsch, and Hoos's local moves. We prove this graph is connected. We give an explicit algorithm to construct a valid plane tree from a primary sequence, assuming that at least one valid plane tree exists. The tree produced by our algorithm has other useful characterizations, including a uniqueness condition defined by local moves. We also study enumerative properties of valid plane trees, analyzing how the number of valid plane trees depends on the choice of sequence length and alphabet size. Finally we show that the proportion of words with at least one valid plane tree goes to zero as the word size increases. We also give some open questions.


title: Impact of the symmetry energy on nuclear pasta phases and crust-core transition in neutron stars

abstract: We study the impact of the symmetry energy on properties of nuclear pasta phases and crust-core transition in neutron stars. We perform a self-consistent Thomas--Fermi calculation employing the relativistic mean-field model. The properties of pasta phases presented in the inner crust of neutron stars are investigated and the crust-core transition is examined. It is found that the slope of the symmetry energy plays an important role in determining the pasta phase structure and the crust-core transition. The correlation between the symmetry energy slope and the crust-core transition density obtained in the Thomas--Fermi approximation is consistent with that predicted by the liquid-drop model.


title: Observation of topological transition of Fermi surface from a spindle-torus to a torus in large bulk Rashba spin-split BiTeCl

abstract: The recently observed large Rashba-type spin splitting in the BiTeX (X = I, Br, Cl) bulk states due to the absence of inversion asymmetry and large charge polarity enables observation of the transition in Fermi surface topology from spindle-torus to torus with varying the carrier density. These BiTeX systems with high spin-orbit energy scales offer an ideal platform for achieving practical spintronic applications and realizing non-trivial phenomena such as topological superconductivity and Majorana fermions. Here we use Shubnikov-de Haas oscillations to investigate the electronic structure of the bulk conduction band of BiTeCl single crystals with different carrier densities. We observe the topological transition of the Fermi surface (FS) from a spindle-torus to a torus. The Landau level fan diagram reveals the expected non-trivial {\pi} Berry phase for both the inner and outer FSs. Angle-dependent oscillation measurements reveal three-dimensional FS topology when the Fermi level lies in the vicinity of the Dirac point. All the observations are consistent with large Rashba spin-orbit splitting in the bulk conduction band.


title: Multi-phase IRC Botnet and Botnet Behavior Detection Model

abstract: Botnets are considered one of the most dangerous and serious security threats facing the networks and the Internet. Comparing with the other security threats, botnet members have the ability to be directed and controlled via C&C messages from the botmaster over common protocols such as IRC and HTTP, or even over covert and unknown applications. As for IRC botnets, general security instances like firewalls and IDSes do not provide by themselves a viable solution to prevent them completely. These devices could not differentiate well between the legitimate and malicious traffic of the IRC protocol. So, this paper is proposing an IDS-based and multi-phase IRC botnet and botnet behavior detection model based on C&C responses messages and malicious behaviors of the IRC bots inside the network environment. The proposed model has been evaluated on five network traffic traces from two different network environments (Virtual network and DARPA 2000 Windows NT Attack Data Set). The results show that the proposed model could detect all the infected IRC botnet member(s), state their current status of attack, filter their malicious IRC messages, pass the other normal IRC messages and detect the botnet behavior regardless of the botnet communication protocol with very low false positive rate. The proposed model has been compared with some of the existing and well-known approaches, including BotHunter, BotSniffer and Rishi regarding botnet characteristics taken in each approach. The comparison showed that the proposed model has made a progress on the comparative models by not to rely on a certain time window or specific bot signatures.


title: Certain homotopy properties related to $\text{map}(\Sigma^n \mathbb{C} P^2,S^m)$

abstract: For given spaces $X$ and $Y$, let $map(X,Y)$ and $map_\ast(X,Y)$ be the unbased and based mapping spaces from $X$ to $Y$, equipped with compact-open topology respectively. Then let $map(X,Y;f)$ and $map_\ast(X,Y;g)$ be the path component of $map(X,Y)$ containing $f$ and $map_\ast(X,Y)$ containing $g$, respectively. In this paper, we compute cohomotopy groups of suspended complex plane $\pi^{n+m}(\Sigma^n \mathbb{C} P^2)$ for $m=6,7$. Using these results, we classify path components of the spaces $map(\Sigma^n \mathbb{C} P^2,S^m)$ up to homotopy equivalent. We also determine the generalized Gottlieb groups $G_n(\mathbb{C} P^2,S^m)$. Finally, we compute homotopy groups of mapping spaces $map(\Sigma^n \mathbb{C}P^2,S^m;f)$ for all generators $[f]$ of $[\Sigma^n \mathbb{C} P^2,S^m]$, and Gottlieb groups of mapping components containing constant map $map(\Sigma^n \mathbb{C} P^2,S^m;0)$.


title: On non-Archimedean curves omitting few components and their arithmetic analogues

abstract: Let k be an algebraically closed field complete with respect to a non-Archimedean absolute value of arbitrary characteristic. Let D_1,...,D_n be effective nef divisors intersecting transversally in an n-dimensional nonsingular projective variety X. We study the degeneracy of non-Archimedean analytic maps from k into $X\setminus \cup_{i=1}^nD_i$ under various geometric conditions. When X is a rational ruled surface and D_1 and D_2 are ample, we obtain a necessary and sufficient condition such that there is no non-Archimedean analytic map from k into $X\setminus D_1 \cup D_2$. Using a dictionary between non-Archimedean Nevanlinna theory and Diophantine approximation, we also study arithmetic analogues of these problems, establishing results on integral points on these varieties over the integers or the ring of integers of an imaginary quadratic field.


title: Deterministic single soliton generation and compression in microring resonators avoiding the chaotic region

abstract: A path within the parameter space of detuning and pump power is demonstrated in order to obtain a single cavity soliton (CS) with certainty in SiN microring resonators in the anomalous dispersion regime. Once the single CS state is reached, it is possible to continue a path to compress it, broadening the corresponding single free spectral range (FSR) Kerr frequency comb. The first step to achieve this goal is to identify the stable regions in the parameter space via numerical simulations of the Lugiato-Lefever equation (LLE). Later, using this identification, we define a path from the stable modulation instability (SMI) region to the stable cavity solitons (SCS) region avoiding the chaotic and unstable regions.


title: Uniform Bounds for Weil-Petersson Curvatures

abstract: We find bounds for Weil-Petersson holomorphic sectional curvature, and the Weil-Petersson curvature operator in several regimes, that do not depend on the topology of the underlying surface. Among other results, we show that the minimal (most negative) eigenvalue of the curvature operator at any point in the Teichm\"uller space $\Teich(S_g)$ of a closed surface $S_g$ of genus $g$ is uniformly bounded away from zero. Restricting to a thick part of $\Teich(S_g)$, we show that the minimal eigenvalue is uniformly bounded below by an explicit constant which does not depend on the topology of the surface but only on the given bound on injectivity radius. We also show that the minimal Weil-Petersson holomorphic sectional curvature of a sufficiently thick hyperbolic surface is comparable to $-1$.


title: Tighter Estimates for epsilon-nets for Disks

abstract: The geometric hitting set problem is one of the basic geometric combinatorial optimization problems: given a set $P$ of points, and a set $\mathcal{D}$ of geometric objects in the plane, the goal is to compute a small-sized subset of $P$ that hits all objects in $\mathcal{D}$. In 1994, Bronniman and Goodrich made an important connection of this problem to the size of fundamental combinatorial structures called $\epsilon$-nets, showing that small-sized $\epsilon$-nets imply approximation algorithms with correspondingly small approximation ratios. Very recently, Agarwal and Pan showed that their scheme can be implemented in near-linear time for disks in the plane. Altogether this gives $O(1)$-factor approximation algorithms in $\tilde{O}(n)$ time for hitting sets for disks in the plane. This constant factor depends on the sizes of $\epsilon$-nets for disks; unfortunately, the current state-of-the-art bounds are large -- at least $24/\epsilon$ and most likely larger than $40/\epsilon$. Thus the approximation factor of the Agarwal and Pan algorithm ends up being more than $40$. The best lower-bound is $2/\epsilon$, which follows from the Pach-Woeginger construction for halfspaces in two dimensions. Thus there is a large gap between the best-known upper and lower bounds. Besides being of independent interest, finding precise bounds is important since this immediately implies an improved linear-time algorithm for the hitting-set problem. The main goal of this paper is to improve the upper-bound to $13.4/\epsilon$ for disks in the plane. The proof is constructive, giving a simple algorithm that uses only Delaunay triangulations. We have implemented the algorithm, which is available as a public open-source module. Experimental results show that the sizes of $\epsilon$-nets for a variety of data-sets is lower, around $9/\epsilon$.


title: Pfaffian Intersections and Multiplicity Cycles

abstract: We consider the problem of estimating the intersection multiplicity between an algebraic variety and a Pfaffian foliation, at every point of the variety. We show that this multiplicity can be majorized at every point $p$ by the local algebraic multiplicity at $p$ of a suitably constructed algebraic cycle. The construction is based on Gabrielov's complex analog of the Rolle-Khovanskii lemma. We illustrate the main result by deriving similar uniform estimates for the complexity of the Milnor fiber of a deformation (under a smoothness assumption) and for the order of contact between an algebraic hypersurface and an arbitrary non-singular one-dimensional foliation. We also use the main result to give an alternative geometric proof for a classical multiplicity estimate in the context of commutative group varieties.


title: Stability of filaments in star-forming clouds and the formation of prestellar cores in them

abstract: It is now widely accepted that dense filaments of molecular gas are integral to the process of stellar birth. While numerical simulations have succeeded in reproducing filamentary structure in turbulent gas and analytic calculations have predicted the formation of dense gas filaments via radial collapse, the exact process(es) that generate/s such filaments which then form prestellar cores within them, is unclear. In this work we therefore study numerically the formation of a dense filament using a relatively simple set-up of a uniform-density cylinder in pressure equilibrium with its confining medium. In particular, we examine if its propensity to form a dense filament and further, to the formation of prestellar cores within this filament bears on the gravitational state of the initial volume of gas. We report a radial collapse leading to the formation of a dense filamentary cloud is likely when the initial volume of gas is at least critically stable (characterised by the approximate equality between the mass line-density for this volume and its maximum value). Though self-gravitating, this volume of gas, however, is not seen to be in free-fall. This post-collapse filament then fragments along its length due to the growth of a Jeans-like instability to form prestellar cores like \emph{beads on a string}. We suggest, dense filaments in typical star-forming clouds classified as gravitationally super-critical under the assumption of : (i) isothermality when in fact, they are not, and (ii) extended radial profiles as against one that is pressure-truncated, thereby causing significant over-estimation of their mass line-density, are unlikely to experience gravitational free-fall. The radial density and temperature profile derived for this post-collapse filament is consistent with that deduced for typical filamentary clouds mapped in recent surveys of nearby star-forming regions.


title: Exposure Time Calculator for Immersion Grating Infrared Spectrograph: IGRINS

abstract: We present an exposure-time calculator (ETC) for the Immersion Grating Infrared Spectrograph (IGRINS). The signal and noise values are calculated by taking into account the telluric background emission and absorption, the emission and transmission of the telescope and instrument optics, and the dark current and read noise of the infrared detector arrays. For the atmospheric transmission, we apply models based on the amount of precipitable water vapor along the line of sight to the target. The ETC produces the expected signal-to-noise ratio (S/N) for each resolution element, given the exposure-time and number of exposures. In this paper, we compare the simulated continuum S/N for the early-type star HD 124683 and the late-type star GSS 32, and the simulated emission line S/N for the H2 rovibrational transitions from the Iris Nebula NGC 7023 with the observed IGRINS spectra. The simulated S/N from the ETC is overestimated by 40-50% for the sample continuum targets.


title: An elementary proof of convergence to the mean-field equations for an epidemic model

abstract: It is common to use a compartmental, fluid model described by a system of ordinary differential equations (ODEs) to model disease spread. In addition to their simplicity, these models are also the mean-field approximations of more accurate stochastic models of disease spread on contact networks. For the simplest case of a stochastic susceptible-infected-susceptible (SIS) process (infection with recovery) on a complete network, it has been shown that the fraction of infected nodes converges to the mean-field ODE as the number of nodes increases. However the proofs are not simple, requiring sophisticated probability, partial differential equations (PDE), or infinite systems of ODEs. We provide a short proof in this case for convergence in mean-square on finite time-intervals using a system of two ODEs and a moment inequality and also obtain the first lower bound on the expected fraction of infected nodes.


title: New Supersoft Supersymmetry Breaking Operators and a Solution to the $\mu$ Problem

abstract: We propose the framework, "generalized supersoft supersymmetry breaking". "Supersoft" models, with $D$-type supersymmetry breaking and heavy Dirac gauginos, are considerably less constrained by the LHC searches than the well studied MSSM. These models also ameliorate the supersymmetric flavor and $CP$ problems. However, previously considered mechanisms for obtaining a natural size Higgsino mass parameter (namely, $\mu$) in supersoft models have been relatively complicated and contrived. Obtaining a $125\gev$ for the mass of the lightest Higgs boson has also been difficult. Additional issues with the supersoft scenario arise from the fact that these models contain new scalars in the adjoint representation of the standard model, which may obtain negative squared-masses, breaking color and generating too large a $T$-parameter. In this work we introduce new operators into supersoft models which can potentially solve all these issues. A novel feature of this framework is that the new $\mu$-term can give unequal masses to the up and down type Higgs fields, and the Higgsinos can be much heavier than the Higgs boson without fine-tuning. However, unequal Higgs and Higgsino masses also remove some attractive features of supersoft susy.


title: Muonless Events in ICAL at INO

abstract: The primary physics signal events in the ICAL at INO are the ${\nu}_{\mu}$ charged current (CC) interactions with a well defined muon track. Apart from these events, ICAL can also detect other types of neutrino interactions, i.e. the electron neutrino charged current interactions and the neutral current events. It is possible to have a dataset containing mostly ${\nu}_e$CC events, by imposing appropriate selection cuts on the events. The ${\nu}_{\mu}$ CC and the neutral current events form the background to these events. This study uses the Monte Carlo generated neutrino events, to design the necessary selection cuts to obtain a ${\nu}_e$ CC rich dataset. An optimized set of constraints are developed which balance the need for improving the purity of the sample and having a large enough event sample. Depending on the constraints used, one can obtain a neutrino data sample, with the purity of ${\nu}_e$ events varying between 55% to 70%.