title: Standing Shock Instability in Advection-Dominated Accretion Flows

abstract: Depending on the values of the energy and angular momentum per unit mass in the gas supplied at large radii, inviscid advection-dominated accretion flows can display velocity profiles with either preshock deceleration or preshock acceleration. Nakayama has shown that these two types of flow configurations are expected to have different stability properties. By employing the Chevalier \& Imamura linearization method and the Nakayama instability boundary conditions, we discover that there are regions of parameter space where disks/shocks with outflows can be stable or unstable. In regions of instability, we find that preshock deceleration is always unstable to the zeroth mode with zero frequency of oscillation, but is always stable to the fundamental mode and overtones. Furthermore, we also find that preshock acceleration is always unstable to the zeroth mode and that the fundamental mode and overtones become increasingly less stable as the shock location moves away from the horizon when the disk half-height expands above $\sim 12$ gravitational radii at the shock radius. In regions of stability, we demonstrate the zeroth mode to be stable for the velocity profiles that exhibit preshock acceleration and deceleration. Moreover, for models that are linearly unstable, our model suggests the possible existence of quasi-periodic oscillations (QPOs) with ratios 2:3 and 3:5. These ratios are believed to occur in stellar and supermassive black hole candidates, for example, in GRS 1915+105 and Sgr A*, respectively. We expect that similar QPO ratios also exist in regions of stable shocks.


title: Design of Massive-MIMO-NOMA with Limited Feedback

abstract: In this letter, a low-feedback non-orthogonal multiple access (NOMA) scheme using massive multiple-input multiple-output (MIMO) transmission is proposed. In particular, the proposed scheme can decompose a massive-MIMO-NOMA system into multiple separated single-input single-output NOMA channels, and analytical results are developed to evaluate the performance of the proposed scheme for two scenarios, with perfect user ordering and with one-bit feedback, respectively.


title: Production of heavy Higgs bosons and decay into top quarks at the LHC

abstract: We investigate the production of heavy, neutral Higgs boson resonances and their decays to top-quark top-antiquark ($t\bar{t}$) pairs at the Large Hadron Collider (LHC) at next-to-leading order (NLO) in the strong coupling of quantum chromodynamics (QCD). The NLO corrections to heavy Higgs boson production and the Higgs-QCD interference are calculated in the large $m_t$ limit with an effective K-factor rescaling. The nonresonant $t\bar{t}$ background is taken into account at NLO QCD including weak-interaction corrections. In order to consistently determine the total decay widths of the heavy Higgs bosons, we consider for definiteness the type-II two-Higgs-doublet extension of the standard model and choose three parameter scenarios that entail two heavy neutral Higgs bosons with masses above the $t\bar{t}$ threshold and unsuppressed Yukawa couplings to top quarks. For these three scenarios we compute, for the LHC operating at 13 TeV, the $t\bar{t}$ cross section and the distributions of the $t\bar{t}$ invariant mass, of the transverse top-quark momentum and rapidity, and of the cosine of the Collins-Soper angle with and without the two heavy Higgs resonances. For selected $M_{t\bar{t}}$ bins we estimate the significances for detecting a heavy Higgs signal in the $t\bar{t}$ dileptonic and lepton plus jets decay channels.


title: Model-Driven Automatic Tiling with Cache Associativity Lattices

abstract: Traditional compiler optimization theory distinguishes three separate classes of cache miss -- Cold, Conflict and Capacity. Tiling for cache is typically guided by capacity miss counts. Models of cache function have not been effectively used to guide cache tiling optimizations due to model error and expense. Instead, heuristic or empirical approaches are used to select tilings. We argue that conflict misses, traditionally neglected or seen as a small constant effect, are the only fundamentally important cache miss category, that they form a solid basis by which caches can become modellable, and that models leaning on cache associatvity analysis can be used to generate cache performant tilings. We develop a mathematical framework that expresses potential and actual cache misses in associative caches using Associativity Lattices. We show these lattices to possess two theoretical advantages over rectangular tiles -- volume maximization and miss regularity. We also show that to generate such lattice tiles requires, unlike rectangular tiling, no explicit, expensive lattice point counting. We also describe an implementation of our lattice tiling approach, show that it can be used to give speedups of over 10x versus unoptimized code, and despite currently only tiling for one level of cache, can already be competitive with the aggressive compiler optimizations used in general purposes compares such as GCC and Intel's ICC. We also show that the tiling approach can lead to reasonable automatic parallelism when compared to existing auto-threading compilers.


title: A new free-floating planet in the Upper Scorpius association

abstract: We report on a deep photometric survey covering an area of 1.17 deg$^2$ in the young Upper Scorpius stellar association using VIMOS $Iz$ and UKIDSS $ZJHK$ data taking several years apart. The search for the least massive population of Upper Scorpius ($\sim$5-10 Myr, 145 pc) is performed on the basis of various optical and infrared color-color and color-magnitude diagrams, including WISE photometry, in the magnitude interval $J$=14.5-19 mag (completeness), which corresponds to substellar masses from 0.028 through 0.004 M$_\odot$ at the age and distance of Upper Scorpius. We also present the proper motion analysis of the photometric candidates, finding that two objects successfully pass all photometric and astrometric criteria for membership in the young stellar association. One of them, UScoJ155150.2$-$213457, is a new discovery. We obtained low resolution, near-infrared spectroscopy ($R\sim$450, 0.85--2.35 $\mu$m) of this new finding using the FIRE instrument. We confirmed its low-gravity atmosphere expected for an Upper Scorpius member (weak alkaline lines, strong VO absorption, peaked $H$-band pseudocontinuum). By comparison with spectroscopic standards, we derive a spectral type of L6$\pm$1, and estimate a mass of $\approx$0.008-0.010 M$_\odot$ for UScoJ155150.2$-$213457. The colors and spectral slope of this object resemble those of other young, cool members of Upper Scorpius and $\sigma$ Orionis ($\sim$3 Myr) and field, high gravity dwarfs of related classification in contrast with the very red indices of field, low gravity, L-type dwarfs of intermediate age. UScoJ155150.2$-$213457, which does not show infrared flux excesses up to 4.5 $\mu$m, becomes one of the least massive and latest type objects known in the entire Upper Scorpius stellar association.


title: Maximum lifetime broadcasting problem in sensor networks

abstract: We solve the maximum lifetime problem for the point-to-point and point-to-multipoint broadcast data transmission in one dimensional regular sensor network. Based on the analytical solution of the problem for one dimension we propose an algorithm solving the maximum lifetime broadcasting problem for point-to-point data transmission for any dimension.


title: Fano resonance scattering in waveguide with an impedance boundary condition

abstract: Sound propagation in a waveguide lined with one section of locally reactive material is studied by resonance scattering approach. The objective is to understand the effects of mode coupling in the lined section on the transmission. It is shown that a transmission zero is present in the vicinity of a resonance peak when a numerically real resonance frequency of the open lined section (opened to infinities through the rigid parts of the waveguide) is crossed. The transmission zero and immediate resonance peak form a Fano resonance, it has been explained as an interaction between a resonance and the non-resonant background. The real resonance frequency and its corresponding trapped mode are formed by the interferences (couplings) between two neighbor modes with complex resonance frequencies. It is also linked to the avoided crossing of eigenvalues and the exceptional point. The scattering matrix is expressed in terms of a matrix $\mathsf{H_{eff}}$ which describes approximately the complex resonances in the open lined section. With the aid of the eigenvalues and eigenfunctions of matrix $\mathsf{H_{eff}}$, the traditional acoustic resonance scattering formula can be extended to describe the coupling effects between the open lined section and the rigid parts of the waveguide.


title: An alternative aggregate preference ranking algorithm to assess environmental effects on macrobenthic abundance in coastal water

abstract: Coastal marine waters are ranked among the most important aquatic ecosystems on earth in terms of ecological and economic significance. Since the Industrial Revolution, human activities have drastically changed coastal marine ecosystems. The development of rules and regulations to protect these ecosystems against human activities requires availability of environmental assessment standard. This necessitates the identification of the key parameters that reflect condition of the coastal water ecosystem. Macrobenthic assemblages are recognized to rapidly respond to changes in the quality of water or habitat. Therefore, it would be useful to study the population of macrobenthos and assess the influential factors on the growth of this species. This study is categorized as multidisciplinary approach which contains two perspectives, ecological and mathematical. In the ecological section, the effect of the water quality parameters (e.g. pH, temperature, dissolved oxygen and salinity) and the sediment characteristics on the macrobenthic abundance is studied. A total of 432 samples were collected and analyzed from four touristic costal locations (at various distances form the coast) of Penang National Park to investigate the spatial change of macrobenthic assemblage. In a mathematical perspective, this paper pursues a new algorithm based on the performance evaluation methods. For this purpose, first, Data Envelopment Analysis (DEA) is employed to evaluate a group of Decision Making Units (DMUs). Consequently, inputs of the mentioned DMUs are considered as alternatives (or candidates), and using a modified DEA model that is categorized as aggregating preference ranking method, the influence of inputs in efficiency of DMUs is investigated.


title: The Universal Transverse Momentum Dependent Soft Function at NNLO

abstract: All (un)polarized transverse momentum dependent functions (TMDs), both distribution and fragmentation functions, are defined with the same universal soft function, which cancels spurious rapidity divergences within an individual TMD and renders them well-defined hadronic quantities. Moreover, it is independent of the kinematics, whether it is Drell-Yan, deep inelastic scattering or $e^+e^-\rightarrow 2$ hadrons. In this paper we provide this soft function at next-to-next-to-leading order (NNLO), necessary for the calculation of all TMDs at the same order, and to perform the resummation of large logarithms at next-to-next-to-next-to-leading-logarithmic accuracy. From the results we obtain the $D$ function at NNLO, which governs the evolution of all TMDs. This work represents the first independent and direct calculation of this quantity. Given the all order relation through a Casimir scaling between the soft function relevant for gluon TMDs and the one for quark TMDs, we also obtain the first at NNLO. The used regularization method to deal with the rapidity divergences is discussed as well.


title: Improved determination of hadron matrix elements using the variational method

abstract: The extraction of hadron form factors in lattice QCD using the standard two- and three-point correlator functions has its limitations. One of the most commonly studied sources of systematic error is excited state contamination, which occurs when correlators are contaminated with results from higher energy excitations. We apply the variational method to calculate the axial vector current gA and compare the results to the more commonly used summation and two-exponential fit methods. The results demonstrate that the variational approach offers a more efficient and robust method for the determination of nucleon matrix elements.


title: Twist deformations leading to kappa-Poincare Hopf algebra and their application to physics

abstract: We consider two twist operators that lead to kappa-Poincare Hopf algebra, the first being an Abelian one and the second corresponding to a light-like kappa-deformation of Poincare algebra. The advantage of the second one is that it is expressed solely in terms of Poincare generators. In contrast to this, the Abelian twist goes out of the boundaries of Poincare algebra and runs into envelope of the general linear algebra. Some of the physical applications of these two different twist operators are considered. In particular, we use the Abelian twist to construct the statistics flip operator compatible with the action of deformed symmetry group. Furthermore, we use the light-like twist operator to define a star product and subsequently to formulate a free scalar field theory compatible with kappa-Poincare Hopf algebra and appropriate for considering the interacting phi^4 scalar field model on kappa-deformed space.


title: High-responsivity vertical-illumination Si/Ge uni-traveling-carrier photodiodes based on silicon-on-insulator substrate

abstract: Si/Ge uni-traveling carrier photodiodes exhibit higher output current when space-charge effects are overcome and thermal effects are suppressed, which is highly beneficial for increasing the dynamic range of various microwave photonic systems and simplifying high-bit-rate digital receivers in different applications. From the point of view of packaging, detectors with vertical-illumination configuration can be easily handled by pick-and-place tools and are a popular choice for making photo-receiver modules. However, vertical-illumination Si/Ge uni-traveling carrier (UTC) devices suffer from inter-constraint between high speed and high responsivity. Here, we report a high responsivity vertical-illumination Si/Ge UTC photodiode based on a silicon-on-insulator substrate. The maximum absorption efficiency of the devices was 2.4 times greater than the silicon substrate owing to constructive interference. The Si/Ge UTC photodiode was successfully fabricated and had a dominant responsivity at 1550 nm of 0.18 A/W, a 50% improvement even with a 25% thinner Ge absorption layer.


title: Visco-thermal effects in acoustic metamaterials: from total transmission to total reflection and high absorption

abstract: We theoretically and experimentally investigate visco-thermal effects on the acoustic propagation through metamaterials consisting of rigid slabs with subwavelength slits embedded in air. We demonstrate that this unavoidable loss mechanism is not merely a refinement, but it plays a dominant role in the actual acoustic response of the structure. Specifically, in the case of very narrow slits, the visco-thermal losses avoid completely the excitation of Fabry-Perot resonances, leading to 100% reflection. This is exactly opposite to the perfect transmission predicted in the idealised lossless case. Moreover, for a wide range of geometrical parameters, there exists an optimum slit width at which the energy dissipated in the structure can be as high as 50%. This work provides a clear evidence that visco-thermal effects are necessary to describe realistically the acoustic response of locally resonant metamaterials.


title: High-gradient High-charge CW Superconducting RF gun with CsK2Sb photocathode

abstract: High-gradient CW photo-injectors operating at high accelerating gradients promise to revolutionize many sciences and applications. They can establish the basis for super-bright monochromatic X-ray free-electron lasers, super-bright hadron beams, nuclear- waste transmutation or a new generation of microchip production. In this letter we report on our operation of a superconducting RF electron gun with a record-high accelerating gradient at the CsK2Sb photocathode (i.e. ~ 20 MV/m) generating a record-high bunch charge (i.e., 3 nC). We briefly describe the system and then detail our experimental results. This achievement opens new era in generating high-power electron beams with a very high brightness.


title: Multi-domain electromagnetic absorption of triangular quantum rings

abstract: We present a theoretical study of the unielectronic energy spectra, electron localization, and optical absorption of triangular core-shell quantum rings. We show how these properties depend on geometric details of the triangle, such as side thickness or corners' symmetry. For equilateral triangles, the lowest six energy states (including spin) are grouped in an energy shell, are localized only around corner areas, and are separated by a large energy gap from the states with higher energy which are localized on the sides of the triangle. The energy levels strongly depend on the aspect ratio of the triangle sides, i.e., thickness/length ratio, in such a way that the energy differences are not monotonous functions of this ratio. In particular, the energy gap between the group of states localized in corners and the states localized on the sides strongly decreases with increasing the side thickness, and then slightly increases for thicker samples. With increasing the thickness the low-energy shell remains distinct but the spatial distribution of these states spreads. The behavior of the energy levels and localization leads to a thickness dependent absorption spectrum where one transition may be tuned in the THz domain and a second transition can be tuned from THz to the infrared range of electromagnetic spectrum. We show how these features may be further controlled with an external magnetic field. In this work the electron-electron Coulomb repulsion is neglected.


title: Lyman-\alpha{} Emitters in the context of hierarchical galaxy formation: predictions for VLT/MUSE surveys

abstract: The VLT Multi Unit Spectroscopic Explorer (MUSE) integral-field spectrograph can detect Ly\alpha{} emitters (LAE) in the redshift range $2.8 \lesssim z \lesssim 6.7$ in a homogeneous way. Ongoing MUSE surveys will notably probe faint Ly\alpha{} sources that are usually missed by current narrow-band surveys. We provide quantitative predictions for a typical wedding-cake observing strategy with MUSE based on mock catalogs generated with a semi-analytic model of galaxy formation coupled to numerical Ly\alpha{} radiation transfer models in gas outflows. We expect $\approx$ 1500 bright LAEs ($F_{Ly\alpha}$ $\gtrsim$ $10^{-17}$ erg s$^{-1}$ cm$^{-2}$) in a typical Shallow Field (SF) survey carried over $\approx$ 100 arcmin$^2$, and $\approx$ 2,000 sources as faint as $10^{-18}$ erg s$^{-1}$ cm$^{-2}$ in a Medium-Deep Field (MDF) survey over 10 arcmin$^2$. In a typical Deep Field (DF) survey of 1 arcmin$^2$, we predict that $\approx$ 500 extremely faint LAEs ($F_{Ly\alpha}$ $\gtrsim$ $4 \times 10^{-19}$ erg s$^{-1}$ cm$^{-2}$) will be found. Our results suggest that faint Ly\alpha{} sources contribute significantly to the cosmic Ly\alpha{} luminosity and SFR budget. While the host halos of bright LAEs at z $\approx$ 3 and 6 have descendants with median masses of $2 \times 10^{12}$ and $5 \times 10^{13}$ $M_{\odot}$ respectively, the faintest sources detectable by MUSE at these redshifts are predicted to reside in halos which evolve into typical sub-$L^{*}$ and $L^{*}$ galaxy halos at z = 0. We expect typical DF and MDF surveys to uncover the building blocks of Milky Way-like objects, even probing the bulk of the stellar mass content of LAEs located in their progenitor halos at z $\approx$ 3.


title: Rhodium Nanoparticles for Ultraviolet Plasmonics

abstract: The non-oxidizing catalytic noble metal rhodium is introduced for ultraviolet plasmonics. Planar tripods of 8 nm Rh nanoparticles, synthesized by a modified polyol reduction method, have a calculated local surface plasmon resonance near 330 nm. By attaching p-aminothiophenol, local field-enhanced Raman spectra and accelerated photo-damage were observed under near-resonant ultraviolet illumination, while charge transfer simultaneously increased fluorescence for up to 13 minutes. The combined local field enhancement and charge transfer demonstrate essential steps toward plasmonically-enhanced ultraviolet photocatalysis.


title: Geometric analysis of pathways dynamics: application to versatility of TGF-{\beta} receptors

abstract: We propose a new geometric approach to describe the qualitative dynamics of chemical reactions networks. By this method we identify metastable regimes, defined as low dimensional regions of the phase space close to which the dynamics is much slower compared to the rest of the phase space. Given the network topology and the orders of magnitude of kinetic parameters, the number of such metastable regimes is finite. The dynamics of the network can be described as a sequence of jumps from one metastable regime to another. We show that a geometrically computed connectivity graph restricts the set of possible jumps. We also provide finite state machine (Markov chain) models for such dynamic changes. Applied to signal transduction models, our approach unravels dynamical and functional capacities of signaling pathways, as well as parameters responsible for specificity of the pathway response. In particular, for a model of TGF$\beta$ signalling, we find that the ratio of TGFBR1 to TGFBR2 concentrations can be used to discriminate between metastable regimes. Using expression data from the NCI60 panel of human tumor cell lines, we show that aggressive and non-aggressive tumour cell lines function in different metastable regimes and can be distinguished by measuring the relative concentrations of receptors of the two types.


title: Semiparametric Estimation of a CES Demand System with Observed and Unobserved Product Characteristics

abstract: We develop a characteristics based demand estimation framework for the Marshallian demand system obtained by solving a budget-constrained constant elasticity of substitution (CES) utility maximization problem. From our Marshallian CES demand system, we derive the same market share equation of Berry (1994); Berry, Levinsohn, and Pakes (1995)'s characteristics based logit demand system. Our CES demand estimation framework can accommodate zero predicted and observed market shares by conceptually separating the whether-to-buy decision and how-much-to-buy decision. Furthermore, the estimator we suggest allows a tractable semiparametric estimation strategy that is flexible regarding the distribution of unobservable product characteristics. We apply our framework to scanner data on cola sales, where we show estimated demand curves can be upward sloping if zero market shares are not accommodated properly.


title: Rayleigh Scattering in the Atmosphere of the Warm Exo-Neptune GJ 3470b

abstract: GJ 3470b is a warm Neptune-size planet transiting a M dwarf star. Like the handful of other small exoplanets for which transmission spectroscopy has been obtained, GJ 3470b exhibits a flat spectrum in the near- and mid-infrared. Recently, a tentative detection of Rayleigh scattering in its atmosphere has been reported. This signal manifests itself as an observed increase of the planetary radius as a function of decreasing wavelength in the visible. We set out to verify this detection and observed several transits of this planet with the LCOGT network and the Kuiper telescope in four different bands (Sloan g', Sloan i', Harris B and Harris V). Our analysis reveals a strong Rayleigh scattering slope, thus confirming previous results. This makes GJ 3470b the smallest known exoplanet with a detection of Rayleigh scattering. We find that the most plausible scenario is a hydrogen/helium-dominated atmosphere covered by clouds which obscure absorption features in the infrared and hazes which give rise to scattering in the visible. Our results demonstrate the feasibility of exoplanet atmospheric characterization from the ground, even with meter-class telescopes.


title: Gravitational contraction versus Supernova driving and the origin of the velocity dispersion-size relation in molecular clouds

abstract: Molecular cloud observations show that clouds have non-thermal velocity dispersions that scale with the cloud size as $\sigma\propto R^{1/2}$ at constant surface density, and for varying surface density scale with both the cloud`s size and surface density, $\sigma^2 \propto R \Sigma$. The energy source driving these chaotic motions remains poorly understood. We describe the velocity dispersions observed in a cloud population formed in a kiloparsec-scale numerical simulation of a magnetized, supernova-driven, self-gravitating, interstellar medium, including diffuse heating and radiative cooling. We compare the relationships between velocity dispersion, size, and surface density measured in the simulated cloud population to those found in observations of Galactic molecular clouds. We find that external supernova explosions can not drive turbulent motions of the observed magnitudes within dense clouds. On the other hand, self-gravity also induces non-thermal motions as gravitationally bound clouds begin to collapse in our model, and by doing so their internal velocity dispersions recover the observed relations. Energy conservation suggests that the observed behavior is consistent with the kinetic energy being proportional to the gravitational energy. However, the clouds in our model show no sign of reaching a stable equilibrium state at any time, even for strongly magnetized clouds. We conclude that gravitationally bound molecular clouds are always in a state of gravitational collapse and their properties are a natural result of this chaotic collapse. In order to agree with observed star formation efficiencies, this process must be terminated by the early destruction of the clouds, presumably from internal stellar feedback.


title: Electrostatic energy barriers from dielectric membranes upon approach of translocating DNA molecules

abstract: We probe the electrostatic cost associated with the approach phase of DNA translocation events. Within an analytical theory at the Debye-Huckel level, we calculate the electrostatic free energy of a rigid DNA molecule interacting with a dielectric membrane. For carbon or silicon based low permittivity neutral membranes, the DNA molecule experiences a repulsive energy barrier between 10 kBT and 100 kBT. In the case of engineered membranes with high dielectric permittivity, the membrane surface attracts the DNA with an energy of the same magnitude. Both the repulsive and attractive interactions result from image-charge effects and their magnitude survive even for the thinnest graphene-based membranes of size d~6 A. For weakly charged membranes, the electrostatic free energy is always attractive at large separation distances but switches to repulsive close to the membrane surface. We also characterise the polymer length dependence of the interaction energy. For specific values of the membrane charge density, low permittivity membranes repel short polymers but attract long polymers. Our results can be used to control the strong electrostatic free energy of DNA-membrane interactions prior to translocation events by chemical engineering of the relevant system parameters.


title: blavaan: Bayesian structural equation models via parameter expansion

abstract: This article describes blavaan, an R package for estimating Bayesian structural equation models (SEMs) via JAGS and for summarizing the results. It also describes a novel parameter expansion approach for estimating specific types of models with residual covariances, which facilitates estimation of these models in JAGS. The methodology and software are intended to provide users with a general means of estimating Bayesian SEMs, both classical and novel, in a straightforward fashion. Users can estimate Bayesian versions of classical SEMs with lavaan syntax, they can obtain state-of-the-art Bayesian fit measures associated with the models, and they can export JAGS code to modify the SEMs as desired. These features and more are illustrated by example, and the parameter expansion approach is explained in detail.


title: H\"older estimates and large time behavior for a nonlocal doubly nonlinear evolution

abstract: The nonlinear and nonlocal PDE $$ |v_t|^{p-2}v_t+(-\Delta_p)^sv=0 \, , $$ where $$ (-\Delta_p)^s v\, (x,t)=2 \,\text{PV} \int_{\mathbb{R}^n}\frac{|v(x,t)-v(x+y,t)|^{p-2}(v(x,t)-v(x+y,t))}{|y|^{n+sp}}\, dy, $$ has the interesting feature that an associated Rayleigh quotient is non-increasing in time along solutions. We prove the existence of a weak solution of the corresponding initial value problem which is also unique as a viscosity solution. Moreover, we provide H\"older estimates for viscosity solutions and relate the asymptotic behavior of solutions to the eigenvalue problem for the fractional $p$-Laplacian.


title: Rescue of endemic states in interconnected networks with adaptive coupling

abstract: We study the Susceptible-Infected-Susceptible model of epidemic spreading on two layers of networks interconnected by adaptive links, which are rewired at random to avoid contacts between infected and susceptible nodes at the interlayer. We find that the rewiring reduces the effective connectivity for the transmission of the disease between layers, and may even totally decouple the networks. Weak endemic states, in which the epidemics spreads only if the two layers are interconnected, show a transition from the endemic to the healthy phase when the rewiring overcomes a threshold value that depends on the infection rate, the strength of the coupling and the mean connectivity of the networks. In the strong endemic scenario, in which the epidemics is able to spread on each separate network, the prevalence in each layer decreases when increasing the rewiring, arriving to single network values only in the limit of infinitely fast rewiring. We also find that finite-size effects are amplified by the rewiring, as there is a finite probability that the epidemics stays confined in only one network during its lifetime.


title: Identifying the Absorption Bump with Deep Learning

abstract: The pervasive interstellar dust grains provide significant insights to understand the formation and evolution of the stars, planetary systems, and the galaxies, and may harbor the building blocks of life. One of the most effective way to analyze the dust is via their interaction with the light from background sources. The observed extinction curves and spectral features carry the size and composition information of dust. The broad absorption bump at 2175 Angstrom is the most prominent feature in the extinction curves. Traditionally, statistical methods are applied to detect the existence of the absorption bump. These methods require heavy preprocessing and the co-existence of other reference features to alleviate the influence from the noises. In this paper, we apply Deep Learning techniques to detect the broad absorption bump. We demonstrate the key steps for training the selected models and their results. The success of Deep Learning based method inspires us to generalize a common methodology for broader science discovery problems. We present our on-going work to build the DeepDis system for such kind of applications.


title: Globular Cluster Populations: Results Including S$^4$G Late-Type Galaxies

abstract: Using 3.6 and 4.5$\mu$m images of 73 late-type, edge-on galaxies from the S$^4$G survey, we compare the richness of the globular cluster populations of these galaxies to those of early type galaxies that we measured previously. In general, the galaxies presented here fill in the distribution for galaxies with lower stellar mass, M$_*$, specifically $\log({\rm M}_*/{\rm M}_\odot) < 10$, overlap the results for early-type galaxies of similar masses, and, by doing so, strengthen the case for a dependence of the number of globular clusters per $10^9\ {\rm M}_\odot$ of galaxy stellar mass, T$_{\rm N}$, on M$_*$. For $8.5 < \log ({\rm M}_*/{\rm M}_\odot) < 10.5$ we find the relationship can be satisfactorily described as T$_{\rm N} = ({\rm M}_*/10^{6.7})^{-0.56}$ when M$_*$ is expressed in solar masses. The functional form of the relationship is only weakly constrained and extrapolation outside this range is not advised. Our late-type galaxies, in contrast to our early-types, do not show the tendency for low mass galaxies to split into two T$_{\rm N}$ families. Using these results and a galaxy stellar mass function from the literature, we calculate that in a volume limited, local Universe sample, clusters are most likely to be found around fairly massive galaxies (M$_* \sim 10^{10.8}$ M$_\odot$) and present a fitting function for the volume number density of clusters as a function of parent galaxy stellar mass. We find no correlation between T$_{\rm N}$ and large-scale environment, but do find a tendency for galaxies of fixed M$_*$ to have larger T$_{\rm N}$ if they have converted a larger proportion of their baryons into stars.


title: Deformation of an Elastic Beam due to Viscous Flow in an Embedded Parallel Channel Network

abstract: Elastic deformation due to embedded fluidic networks is currently studied in the context of soft-actuators and soft-robotic applications. In this work, we analyze interaction between the elastic deflection of a slender beam and viscous flow within a long serpentine channel, embedded in the elastic beam. The channel is positioned asymmetrically with regard to the midplane of the beam, and thus pressure within the channel creates a local moment deforming the beam. We focus on creeping flows and small deflections of the elastic beam and obtain, in leading order, a fourth-order partial integro-differential equation governing the time-dependent deflection field. This relation enables the design of complex time-dependent deformation patterns of beams with embedded channel networks, including inertia-like standing and moving wave solutions in configurations with negligible inertia.


title: Bounded Stability in Networked Systems with Parameter Mismatch and Adaptive Decentralized Estimation

abstract: Here, we study the ultimately bounded stability of network of mismatched systems using Lyapunov direct method. The upper bound on the error of oscillators from the center of the neighborhood is derived. Then the performance of an adaptive compensation via decentralized control is analyzed. Finally, the analytical results for a network of globally connected Lorenz oscillators are verified.


title: Co-modularity and Detection of Co-communities

abstract: This paper introduces the notion of co-modularity, to co-cluster observations of bipartite networks into co-communities. The task of co-clustering is to group together nodes of one type with nodes of another type, according to the interactions that are the most similar. The novel measure of co-modularity is introduced to assess the strength of co-communities, as well as to arrange the representation of nodes and clusters for visualisation, and to define an objective function for optimisation. We demonstrate the usefulness of our proposed methodology on simulated data, and with examples from genomics and consumer-product reviews.


title: A Block Regression Model for Short-Term Mobile Traffic Forecasting

abstract: Accurate mobile traffic forecast is important for efficient network planning and operations. However, existing traffic forecasting models have high complexity, making the forecasting process slow and costly. In this paper, we analyze some characteristics of mobile traffic such as periodicity, spatial similarity and short term relativity. Based on these characteristics, we propose a \emph{Block Regression} ({BR}) model for mobile traffic forecasting. This model employs seasonal differentiation so as to take into account of the temporally repetitive nature of mobile traffic. One of the key features of our {BR} model lies in its low complexity since it constructs a single model for all base stations. We evaluate the accuracy of {BR} model based on real traffic data and compare it with the existing models. Results show that our {BR} model offers equal accuracy to the existing models but has much less complexity.


title: Local existence of solutions to the Euler-Poisson system, including densities without compact support

abstract: Local existence and well posedness for a class of solutions for the Euler Poisson system is shown. These solutions have a density $\rho$ which either falls off at infinity or has compact support. The solutions have finite mass, finite energy functional and include the static spherical solutions for the adiabatic constant $\gamma=\frac{6}{5}$. The result is achieved by using weighted Sobolev spaces of fractional order and a new non linear estimate which allows to estimate the physical density by the regularised non linear matter variable. Gamblin also has studied this setting but using very different functional spaces. However we believe that the functional setting we use is more appropriate to describe a physical isolated body and more suitable to study the Newtionan limit.


title: Model-based Dashboards for Customer Analytics

abstract: Automating the customer analytics process is crucial for companies that manage distinct customer bases. In such data-rich and dynamic environments, visualization plays a key role in understanding events of interest. These ideas have led to the popularity of analytics dashboards, yet academic research has paid scant attention to these managerial needs. We develop a probabilistic, nonparametric framework for understanding and predicting individual-level spending using Gaussian process priors over latent functions that describe customer spending along calendar time, interpurchase time, and customer lifetime dimensions. These curves form a dashboard that provides a visual model-based representation of purchasing dynamics that is easily comprehensible. The model flexibly and automatically captures the form and duration of the impact of events that influence spend propensity, even when such events are unknown a-priori. We illustrate the use of our Gaussian Process Propensity Model (GPPM) on data from two popular mobile games. We show that the GPPM generalizes hazard and buy-till-you-die models by incorporating calendar time dynamics while simultaneously accounting for recency and lifetime effects. It therefore provides insights about spending propensity beyond those available from these models. Finally, we show that the GPPM outperforms these benchmarks both in fitting and forecasting real and simulated spend data.


title: A pragmatic approach to sovereignty on Mars

abstract: Rising interest in Mars colonization from both private and public sectors necessitates a renewed discussion about sovereignty in space. The non-appropriation principle of the Outer Space Treaty currently prohibits any sovereign claims to celestial bodies, but it remains unclear how this principle should be applied to the peaceful colonization of Mars. Here we develop a pragmatic approach to guide the settlement of Mars, which is based upon a "bounded first possession" model with mandatory planetary parks. Scientists, experts, and leaders will establish planetary park locations and regulations through worldwide community solicitation in order to protect sites of scientific, aesthetic, historical, cultural, environmental, spiritual value. Colonization parties may occupy limited plots of martian land and may claim exclusive economic rights within this zone, while still refraining from any claims to sovereignty. All colonists remain under the legal jurisdiction of their host nation, with conflicts to be resolved diplomatically or through a temporary tribunal system composed of representatives from other Mars colonies. We also propose the formation of a Mars Secretariat as an administrative body with limited power to facilitate communication among parties. Our model for Mars colonization remains consistent with the Outer Space Treaty, but we also recommend revisiting or amending the non-appropriation and province of mankind principles to resolve the ambiguity of how nations, corporations, and individuals may utilize the resources of space.


title: Learning Structured Inference Neural Networks with Label Relations

abstract: Images of scenes have various objects as well as abundant attributes, and diverse levels of visual categorization are possible. A natural image could be assigned with fine-grained labels that describe major components, coarse-grained labels that depict high level abstraction or a set of labels that reveal attributes. Such categorization at different concept layers can be modeled with label graphs encoding label information. In this paper, we exploit this rich information with a state-of-art deep learning framework, and propose a generic structured model that leverages diverse label relations to improve image classification performance. Our approach employs a novel stacked label prediction neural network, capturing both inter-level and intra-level label semantics. We evaluate our method on benchmark image datasets, and empirical results illustrate the efficacy of our model.


title: Two-photon interference from a bright single photon source at telecom wavelengths

abstract: Long-distance quantum communication relies on the ability to efficiently generate and prepare single photons at telecom wavelengths. In many applications these photons must also be indistinguishable such that they exhibit interference on a beamsplitter, which implements effective photon-photon interactions. However, deterministic generation of indistinguishable single photons with high brightness remains a challenging problem. We demonstrate two-photon interference at telecom wavelengths using an InAs/InP quantum dot in a nanophotonic cavity. The cavity enhances the quantum dot emission, resulting in a nearly Gaussian transverse mode profile with high out-coupling efficiency exceeding 36% after multi-photon correction. We also observe Purcell enhanced spontaneous emission rate up to 4. Using this source, we generate linearly polarized, high purity single photons at 1.3 micron wavelength and demonstrate the indistinguishable nature of the emission using a two-photon interference measurement, which exhibits indistinguishable visibilities of 18% without post-selection and 67% with post-selection. Our results provide a promising approach to generate bright, deterministic single photons at telecom wavelength for applications in quantum networking and quantum communication.


title: Topic Model Based Behaviour Modeling and Clustering Analysis for Wireless Network Users

abstract: User behaviour analysis based on traffic log in wireless networks can be beneficial to many fields in real life: not only for commercial purposes, but also for improving network service quality and social management. We cluster users into groups marked by the most frequently visited websites to find their preferences. In this paper, we propose a user behaviour model based on Topic Model from document classification problems. We use the logarithmic TF-IDF (term frequency - inverse document frequency) weighing to form a high-dimensional sparse feature matrix. Then we apply LSA (Latent semantic analysis) to deduce the latent topic distribution and generate a low-dimensional dense feature matrix. K-means++, which is a classic clustering algorithm, is then applied to the dense feature matrix and several interpretable user clusters are found. Moreover, by combining the clustering results with additional demographical information, including age, gender, and financial information, we are able to uncover more realistic implications from the clustering results.


title: Cyclus Archetypes

abstract: The current state of nuclear fuel cycle simulation exists in highly customized form. Satisfying a wide range of users requires model modularity within such a tool. Cyclus is a fuel cycle simulator specifically designed to combat the lack of adaptability of previous generations of simulators. This is accomplished through an agent-based infrastructure and treating time discretely. The Cyclus kernel was developed to allow for models, called archetypes, of differing fidelity and function depending on need of the users. To take advantage of this flexibility, a user must write an archetype for their desired simulation if it does not yet exist within the Cyclus ecosystem. At this stage, a user graduates to the title of archetype developer. Without automation, archetype development is difficult for the uninitiated. This paper presents the framework developed for simplifying the writing of archetypes: the Cyclus preprocessor, or cycpp. cycpp addresses the computer science and software development aspects of archetype development that can be addressed algorithmically, allowing the developer to focus on modeling the physics, social policies, and economics. cycpp passes through the code three times to perform the following tasks: normalizing the code via the C preprocessor, accumulation of notations, and code generation. Not only does this reduce the amount of code a developer must write by approximately an order of magnitude, but the archetypes are automatically validated.


title: External Photoevaporation of the Solar Nebula II: Effects on Disk Structure and Evolution with Non-Uniform Turbulent Viscosity due to the Magnetorotational Instability

abstract: The structure and evolution of protoplanetary disks, especially the radial flows of gas through them, are sensitive to a number of factors. One that has been considered only occasionally in the literature is external photoevaporation by far-ultraviolet (FUV) radiation from nearby, massive stars, despite the fact that nearly half of all disks will experience photoevaporation. Another effect apparently not considered in the literature is a spatially and temporally varying value of $\alpha$ in the disk [where the turbulent viscosity $\nu$ is $\alpha$ times the sound speed C times the disk scale height H]. Here we use the formulation of Bai \& Stone (2011) to relate $\alpha$ to the ionization fraction in the disk, assuming turbulent transport of angular momentum is due to the magnetorotational instability. We find that disk evolution is most sensitive to the surface area of dust. Typically $\alpha \lesssim 10^{-5}$ in the inner disk ($< 2$ AU), rising to $\sim 10^{-1}$ beyond 20 AU. This drastically alters the structure of the disk and the flow of mass through it: while the outer disk rapidly viscously spreads, the inner disk hardly evolves; this leads to a steep surface density profile with a slope < p > $\approx$ 2 - 5 in the 5-30 AU region) that is made steeper by external photoevaporation. We also find that the combination of variable $\alpha$ and external photoevaporation eventually causes gas as close as 3 AU, previously accreting inward, to be drawn outward to the photoevaporated outer edge of the disk. These effects have drastic consequences for planet formation and volatile transport in protoplanetary disks.


title: Origin of the long-term modulation of radio emission of LS I +61 303

abstract: One of the most unusual aspects of the X-ray binary LSI +61 303 is that at each orbit (P1=26.4960 \pm 0.0028 d) one radio outburst occurs whose amplitude is modulated with Plong, a long-term period of more than 4 yr. It is still not clear whether the compact object of the system or the companion Be star is responsible for the long-term modulation. We study here the stability of Plong. Such a stability is expected if Plong is due to periodic (P2) Doppler boosting of periodic (P1) ejections from the accreting compact object of the system. On the contrary it is not expected if Plong is related to variations in the mass loss of the companion Be star. We built a database of 36.8 yr of radio observations of LSI +61 303 covering more than 8 long-term cycles. We performed timing and correlation analysis. In addition to the two dominant features at P1 and P2, the timing analysis gives a feature at Plong=1628 \pm 48 days. The determined value of Plong agrees with the beat of the two dominant features, i.e. Pbeat=1/(\nu1 -\nu2)=1626 \pm 68 d. The correlation coefficient of the radio data oscillates at multiples of Pbeat. Cycles in varying Be stars change in length and disappear after 2-3 cycles following the well-studied case of the binary system zeta Tau. On the contrary, in LSI +61 303 the long-term period is quite stable and repeats itself over the available 8 cycles. The long-term modulation in LSI +61 303 accurately reflects the beat of periodical Doppler boosting (induced by precession) with the periodicity of the ejecta. The peak of the long-term modulation occurs at the coincidence of the maximum number of ejected particles with the maximum Doppler boosting of their emission; this coincidence occurs every 1/(\nu1 - \nu2) and creates the long-term modulation observed in LSI +61 303.


title: Predicting distributions with Linearizing Belief Networks

abstract: Conditional belief networks introduce stochastic binary variables in neural networks. Contrary to a classical neural network, a belief network can predict more than the expected value of the output $Y$ given the input $X$. It can predict a distribution of outputs $Y$ which is useful when an input can admit multiple outputs whose average is not necessarily a valid answer. Such networks are particularly relevant to inverse problems such as image prediction for denoising, or text to speech. However, traditional sigmoid belief networks are hard to train and are not suited to continuous problems. This work introduces a new family of networks called linearizing belief nets or LBNs. A LBN decomposes into a deep linear network where each linear unit can be turned on or off by non-deterministic binary latent units. It is a universal approximator of real-valued conditional distributions and can be trained using gradient descent. Moreover, the linear pathways efficiently propagate continuous information and they act as multiplicative skip-connections that help optimization by removing gradient diffusion. This yields a model which trains efficiently and improves the state-of-the-art on image denoising and facial expression generation with the Toronto faces dataset.


title: Characterization of steady solutions to the 2D Euler equation

abstract: Steady fluid flows have very special topology. In this paper we describe necessary and sufficient conditions on the vorticity function of a 2D ideal flow on a surface with or without boundary, for which there exists a steady flow among isovorticed fields. For this we introduce the notion of an antiderivative (or circulation function) on a measured graph, the Reeb graph associated to the vorticity function on the surface, while the criterion is related to the total negativity of this antiderivative. It turns out that given topology of the vorticity function, the set of coadjoint orbits of the symplectomorphism group admitting steady flows with this topology forms a convex polytope. As a byproduct of the proposed construction, we also describe a complete list of Casimirs for the 2D Euler hydrodynamics: we define generalized enstrophies which, along with circulations, form a complete set of invariants for coadjoint orbits of area-preserving diffeomorphisms on a surface.


title: Finite Type Invariants of w-Knotted Objects IV: Some Computations

abstract: In the previous three papers in this series, [WKO1]-[WKO3] (arXiv:1405.1956, arXiv:1405.1955, and to appear), Z. Dancso and I studied a certain theory of "homomorphic expansions" of "w-knotted objects", a certain class of knotted objects in 4-dimensional space. When all layers of interpretation are stripped off, what remains is a study of a certain number of equations written in a family of spaces $A^w$, closely related to degree-completed free Lie algebras and to degree-completed spaces of cyclic words. The purpose of this paper is to introduce mathematical and computational tools that enable explicit computations (up to a certain degree) in these $A^w$ spaces and to use these tools to solve the said equations and verify some properties of their solutions, and as a consequence, to carry out the computation (up to a certain degree) of certain knot-theoretic invariants discussed in [WKO1]-[WKO3] and in my related paper [KBH] (arXiv:1308.1721).


title: MOEA/D-GM: Using probabilistic graphical models in MOEA/D for solving combinatorial optimization problems

abstract: Evolutionary algorithms based on modeling the statistical dependencies (interactions) between the variables have been proposed to solve a wide range of complex problems. These algorithms learn and sample probabilistic graphical models able to encode and exploit the regularities of the problem. This paper investigates the effect of using probabilistic modeling techniques as a way to enhance the behavior of MOEA/D framework. MOEA/D is a decomposition based evolutionary algorithm that decomposes a multi-objective optimization problem (MOP) in a number of scalar single-objective subproblems and optimizes them in a collaborative manner. MOEA/D framework has been widely used to solve several MOPs. The proposed algorithm, MOEA/D using probabilistic Graphical Models (MOEA/D-GM) is able to instantiate both univariate and multi-variate probabilistic models for each subproblem. To validate the introduced framework algorithm, an experimental study is conducted on a multi-objective version of the deceptive function Trap5. The results show that the variant of the framework (MOEA/D-Tree), where tree models are learned from the matrices of the mutual information between the variables, is able to capture the structure of the problem. MOEA/D-Tree is able to achieve significantly better results than both MOEA/D using genetic operators and MOEA/D using univariate probability models, in terms of the approximation to the true Pareto front.


title: Observation of Quantum Interferences via Light Induced Conical Intersections in Diatomic Molecules

abstract: We observe energy-dependent angle-resolved diffraction patterns in protons from strong-field dissociation of the molecular hydrogen ion H$_2^+$. The interference is a characteristic of dissociation around a laser-induced conical intersection (LICI), which is a point of contact between two surfaces in the dressed 2-dimensional Born-Oppenheimer potential energy landscape of a diatomic molecule in a strong laser field. The interference magnitude and angular period depend strongly on the energy difference between the initial state and the LICI, consistent with coherent diffraction around a cone-shaped potential barrier whose width and thickness depend on the relative energy of the initial state and the cone apex. These findings are supported by numerical solutions of the time-dependent Schr\"{o}dinger equation for similar experimental conditions.


title: On the construction of charged operators inside an eternal black hole

abstract: We revisit the holographic construction of (approximately) local bulk operators inside an eternal AdS black hole in terms of operators in the boundary CFTs. If the bulk operator carries charge, the construction must involve a qualitatively new object: a Wilson line that stretches between the two boundaries of the eternal black hole. This operator - more precisely, its zero mode - cannot be expressed in terms of the boundary currents and only exists in entangled states dual to two-sided geometries, which suggests that it is a state-dependent operator. We determine the action of the Wilson line on the relevant subspaces of the total Hilbert space, and show that it behaves as a local operator from the point of view of either CFT. For the case of three bulk dimensions, we give explicit expressions for the charged bulk field and the Wilson line. Furthermore, we show that when acting on the thermofield double state, the Wilson line may be extracted from a limit of certain standard CFT operator expressions. We also comment on the relationship between the Wilson line and previously discussed mirror operators in the eternal black hole.


title: Quantum modularity and complex Chern-Simons theory

abstract: The Quantum Modularity Conjecture of Zagier predicts the existence of a formal power series with arithmetically interesting coefficients that appears in the asymptotics of the Kashaev invariant at each root of unity. Our goal is to construct a power series from a Neumann-Zagier datum (i.e., an ideal triangulation of the knot complement and a geometric solution to the gluing equations) and a complex root of unity $\zeta$. We prove that the coefficients of our series lie in the trace field of the knot, adjoined a complex root of unity. We conjecture that our series are those that appear in the Quantum Modularity Conjecture and confirm that they match the numerical asymptotics of the Kashaev invariant (at various roots of unity) computed by Zagier and the first author. Our construction is motivated by the analysis of singular limits in Chern-Simons theory with gauge group $SL(2,C)$ at fixed level $k$, where $\zeta^k=1$.


title: Marginalized Two Part Models for Generalized Gamma Family of Distributions

abstract: Positive continuous outcomes with a point mass at zero are prevalent in biomedical research. To model the point mass at zero and to provide marginalized covariate effect estimates, marginalized two part models (MTP) have been developed for outcomes with lognormal and log skew normal distributions. In this paper, we propose MTP models for outcomes from a generalized gamma (GG) family of distributions. In the proposed MTP-GG model, the conditional mean from a two-part model with a three-parameter GG distribution is parameterized to provide regression coefficients that have marginal interpretation. MTP-gamma and MTP-Weibull are developed as special cases of MTP-GG. We derive marginal covariate effect estimators from each model and through simulations assess their finite sample operating characteristics in terms of bias, standard errors, 95% coverage, and rate of convergence. We illustrate the models using data sets from The Medical Expenditure Survey (MEPS) and from a randomized trial of addictive disorders and we provide SAS code for implementation. The simulation results show that when the response distribution is unknown or mis-specified, which is usually the case in real data sets, the MTP-GG is preferable over other models.


title: Dissociation and ionization equilibria of deuterium fluid over a wide range of temperatures and densities

abstract: We investigate the dissociation and ionization equilibria of deuterium fluid over a wide range of temperatures and densities. The partition functions for molecular and atomic species are evaluated, in a statistical-mechanically consistent way, implementing recent developments in the literature and taking high-density effects into account. A new chemical model (free energy function) is introduced in which the fluid is considered as a mixture of diatomic molecules, atoms, ions and free electrons. Intensive short range hard core repulsion is taken into account together with partial degeneracy of free electrons and Coulomb interactions among charged particles. Samples of computational results are presented as a set of isotherms for the degree of ionization, dissociated fraction of molecules, pressure, and specific internal energy for a wide range of densities and temperatures. Predictions from the present model calculations show an improved and sensible physical behavior compared to other results in the literature.


title: Two laws of large numbers for sublinear expectations

abstract: In this paper, we consider the sublinear expectation on bounded random variables. With the notion of uncorrelatedness for random variables under the sublinear expectation, a weak law of large numbers is obtained. With the notion of independence for random variable sequences and regular property for sublinear expectations, we get a strong one. These results are helpful for the application of the law of large numbers in engineering and statistics such as to assess the level of product quality and to explain the relationship between the statistical distribution and the sample distribution.


title: Tables of equation-of-state, thermodynamic properties, and shock Hugoniot for hot dense fluid deuterium

abstract: We present computational results and tables of the equation-of-state, thermodynamic properties, and shock Hugoniot for hot dense fluid deuterium. The present results are generated using a recently developed chemical model that takes into account different high density effects such as Coulomb interactions among charged particles, partial degeneracy, and intensive short range hard core repulsion. Internal partition functions are evaluated in a statistical-mechanically consistent way implementing recent developments in the literature. The shock Hugoniot curve derived from the present tables is in reasonable overall agreement with the Hugoniot derived from the Nova-laser shock wave experiments on liquid deuterium, showing that deuterium has a significantly higher compressibility than predicted by the SESAME tables or by Path Integral Monte Carlo (PIMC) calculations. Computational results are presented as surface plots for the dissociated fraction, degree of ionization, pressure, and specific internal energy for densities ranging from 0.0001 to 40 g/cm3 and temperatures from 2000 to ~ 106 K. Tables for values of the above mentioned quantities in addition to the specific heat at constant pressure, cp, ratio of specific heats, cp/cv, sound speed and Hugoniot curve (for a specific initial state) are presented for practical use.


title: On the Star Formation Law for Spiral and Irregular Galaxies

abstract: A dynamical model for star formation on a galactic scale is proposed in which the interstellar medium is constantly condensing to star-forming clouds on the dynamical time of the average midplane density, and the clouds are constantly being disrupted on the dynamical time scale appropriate for their higher density. In this model, the areal star formation rate scales with the 1.5 power of the total gas column density throughout the main regions of spiral galaxies, and with a steeper power, 2, in the far outer regions and in dwarf irregular galaxies because of the flaring disks. At the same time, there is a molecular star formation law that is linear in the main and outer parts of disks and in dIrrs because the duration of individual structures in the molecular phase is also the dynamical time scale, canceling the additional 0.5 power of surface density. The total gas consumption time scales directly with the midplane dynamical time, quenching star formation in the inner regions if there is no accretion, and sustaining star formation for ~100 Gyr or more in the outer regions with no qualitative change in gas stability or molecular cloud properties. The ULIRG track follows from high densities in galaxy collisions.


title: Local entropy as a measure for sampling solutions in Constraint Satisfaction Problems

abstract: We introduce a novel Entropy-driven Monte Carlo (EdMC) strategy to efficiently sample solutions of random Constraint Satisfaction Problems (CSPs). First, we extend a recent result that, using a large-deviation analysis, shows that the geometry of the space of solutions of the Binary Perceptron Learning Problem (a prototypical CSP), contains regions of very high-density of solutions. Despite being sub-dominant, these regions can be found by optimizing a local entropy measure. Building on these results, we construct a fast solver that relies exclusively on a local entropy estimate, and can be applied to general CSPs. We describe its performance not only for the Perceptron Learning Problem but also for the random $K$-Satisfiabilty Problem (another prototypical CSP with a radically different structure), and show numerically that a simple zero-temperature Metropolis search in the smooth local entropy landscape can reach sub-dominant clusters of optimal solutions in a small number of steps, while standard Simulated Annealing either requires extremely long cooling procedures or just fails. We also discuss how the EdMC can heuristically be made even more efficient for the cases we studied.


title: Competitive Multi-scale Convolution

abstract: In this paper, we introduce a new deep convolutional neural network (ConvNet) module that promotes competition among a set of multi-scale convolutional filters. This new module is inspired by the inception module, where we replace the original collaborative pooling stage (consisting of a concatenation of the multi-scale filter outputs) by a competitive pooling represented by a maxout activation unit. This extension has the following two objectives: 1) the selection of the maximum response among the multi-scale filters prevents filter co-adaptation and allows the formation of multiple sub-networks within the same model, which has been shown to facilitate the training of complex learning problems; and 2) the maxout unit reduces the dimensionality of the outputs from the multi-scale filters. We show that the use of our proposed module in typical deep ConvNets produces classification results that are either better than or comparable to the state of the art on the following benchmark datasets: MNIST, CIFAR-10, CIFAR-100 and SVHN.


title: Theoretical study of transport properties of B40 and its endohedral borospherenes in single-molecule junctions

abstract: C60 fullerene has been studied extensively, as it is considered to be a good candidate for building single-molecule junctions. Here, we theoretically demonstrate that the conductance of single-molecule junctions based on a newly discovered molecule, borospherene (B40), is comparable to that for the C60-based junction with its more delocalized {\pi} electrons. The charge injection efficiency in the B40-based junction is improved, as up to 7 atoms in direct contact with the electrode are possible in the Au-B40-Au junction. Interestingly, a higher number of atoms in direct contact with the electrode does not result in a higher number of conduction channels because of the unique chemical bonding in the B40 molecule, without two-center two- electron bonds. The transport properties of Au-B40-Au junctions can be proved by doping. With a Ca, Sr, or Y atom encapsulated into the B40 cage, the conductance at zero bias increases significantly. Moreover, our calculations show that the lowest unoccupied molecular orbital dominates the low-bias transport, as the thermopower in these junctions is negative. Our study indicates that B40 is an attractive new platform for designing highly conductive single-molecule junctions for future molecular circuits.


title: Discrete one-dimensional coverage process on a renewal process

abstract: We consider the {following} coverage model on $\mathbb{N}$. For each site $i\in \mathbb{N}$ we associate a pair $(\xi_i, R_i)$ where $\{\xi_0, \xi_1, \ldots \}$ is a 1-dimensional {undelayed} discrete renewal point process and $\{R_0,R_1,\ldots\}$ is an i.i.d. sequence of $\mathbb{N}$-valued random variables. At each site where $\xi_i=1$ we start an interval of length $R_i$. Coverage occurs if every site of $\mathbb{N}$ is covered by some interval. We obtain sharp conditions for both, positive and null probability of coverage. As corollaries, we extend results of the literature of rumor processes and discrete one-dimensional Boolean percolation.


title: Strong enhancement of critical current density in both low & high fields and flux pinning mechanism under hydrostatic pressure in optimally doped (Ba,K)Fe2As2 single crystals

abstract: Strong pinning depends on the pinning force strength and number density of effective defects. Using hydrostatic pressure method, we demonstrate that hydrostatic pressure up to 1.2 GPa can significantly enhance flux pinning or Jc by a factor of up to 5 especially in both low and high fields in optimally doped Ba0.6K0.4Fe2As2 crystals. Our analysis on the flux pining mechanism indicate that both pinning centre number density (Np) and pinning force (Fp) are greatly increased by the pressure and contribute to strong pinning.


title: A Mathematical Framework Exhibiting the Emergence of Dynamic Expansion of Task Repertoire in \emph{Pheidole dentata}

abstract: The division of labor (DOL) and task allocation among groups of ants living in a colony is thought to be highly efficient, and key to the robust survival of a colony. A great deal of experimental and theoretical work has been done toward gaining a clear understanding of the evolution of, and underlying mechanisms of these phenomena. Much of this research has utilized mathematical modeling. Here we continue this tradition by developing a mathematical model for a particular aspect of task allocation, known as age-related repertoire expansion, that has been observed in the minor workers of the ant species \emph{Pheidole dentata}. In fact, we present a relatively broad mathematical modeling framework based on the dynamics of the frequency with which members of specific age groups carry out distinct tasks. We apply our modeling approach to a specific task allocation scenario, and compare our theoretical results with experimental data. It is observed that the model predicts perceived behavior, and provides a possible explanation for the aforementioned experimental results.


title: A Quick Test on Rotation Period Clustering for the Small Members of the Koronis Family

abstract: Rotation period clustering in prograde/retrograde rotators might indicate the preliminary indication of the Slivan state in the Koronis family as a result of the YORP effect. We follow the general scenario of dispersion in semimajor axis of the asteroid family members to separate prograde and retrograde rotators in the Koronis family. From the available rotation periods obtained from PTF/iPTF, we can not found the rotation period clustering of objects with H grater than 12 mag in the Koronis family. This could be the result of the intermittent collisional process of small asteroids (D less than ~20 km) which leads to astray Yarkovsky drifting. Measurement of the pole orientations of our sample will verify our preliminary result and validate our method.


title: Net2Net: Accelerating Learning via Knowledge Transfer

abstract: We introduce techniques for rapidly transferring the information stored in one neural net into another neural net. The main purpose is to accelerate the training of a significantly larger neural net. During real-world workflows, one often trains very many different neural networks during the experimentation and design process. This is a wasteful process in which each new model is trained from scratch. Our Net2Net technique accelerates the experimentation process by instantaneously transferring the knowledge from a previous network to each new deeper or wider network. Our techniques are based on the concept of function-preserving transformations between neural network specifications. This differs from previous approaches to pre-training that altered the function represented by a neural net when adding layers to it. Using our knowledge transfer mechanism to add depth to Inception modules, we demonstrate a new state of the art accuracy rating on the ImageNet dataset.


title: Gravitational wave production by Hawking radiation from rotating primordial black holes

abstract: In this paper we analyze in detail a rarely discussed question of gravity wave production from evaporating primordial black holes. These black holes emit gravitons which are, at classical level, registered as gravity waves. We use the latest constraints on their abundance, and calculate the power emitted in gravitons at the time of their evaporation. We then solve the coupled system of equations that gives us the evolution of the frequency and amplitude of gravity waves during the expansion of the universe. The spectrum of gravitational waves that can be detected today depends on multiple factors: fraction of the total energy density which was occupied by primordial black holes, the epoch in which they were formed, and quantities like their mass and angular momentum. We conclude that very small primordial black holes which evaporate before the big-bang nucleosynthesis emit gravitons whose spectral energy fraction today can be as large as $10^{-7.5}$. On the other hand, those which are massive enough so that they still exist now can yield a signal as high as $10^{-6.5}$. However, typical frequencies of the gravity waves from primordial black holes are still too high to be observed with the current and near future gravity wave observations.


title: A New Smooth Approximation to the Zero One Loss with a Probabilistic Interpretation

abstract: We examine a new form of smooth approximation to the zero one loss in which learning is performed using a reformulation of the widely used logistic function. Our approach is based on using the posterior mean of a novel generalized Beta-Bernoulli formulation. This leads to a generalized logistic function that approximates the zero one loss, but retains a probabilistic formulation conferring a number of useful properties. The approach is easily generalized to kernel logistic regression and easily integrated into methods for structured prediction. We present experiments in which we learn such models using an optimization method consisting of a combination of gradient descent and coordinate descent using localized grid search so as to escape from local minima. Our experiments indicate that optimization quality is improved when learning meta-parameters are themselves optimized using a validation set. Our experiments show improved performance relative to widely used logistic and hinge loss methods on a wide variety of problems ranging from standard UC Irvine and libSVM evaluation datasets to product review predictions and a visual information extraction task. We observe that the approach: 1) is more robust to outliers compared to the logistic and hinge losses; 2) outperforms comparable logistic and max margin models on larger scale benchmark problems; 3) when combined with Gaussian- Laplacian mixture prior on parameters the kernelized version of our formulation yields sparser solutions than Support Vector Machine classifiers; and 4) when integrated into a probabilistic structured prediction technique our approach provides more accurate probabilities yielding improved inference and increasing information extraction performance.


title: Adversarial Autoencoders

abstract: In this paper, we propose the "adversarial autoencoder" (AAE), which is a probabilistic autoencoder that uses the recently proposed generative adversarial networks (GAN) to perform variational inference by matching the aggregated posterior of the hidden code vector of the autoencoder with an arbitrary prior distribution. Matching the aggregated posterior to the prior ensures that generating from any part of prior space results in meaningful samples. As a result, the decoder of the adversarial autoencoder learns a deep generative model that maps the imposed prior to the data distribution. We show how the adversarial autoencoder can be used in applications such as semi-supervised classification, disentangling style and content of images, unsupervised clustering, dimensionality reduction and data visualization. We performed experiments on MNIST, Street View House Numbers and Toronto Face datasets and show that adversarial autoencoders achieve competitive results in generative modeling and semi-supervised classification tasks.


title: ABC Implies There are Infinitely Many non-Fibonacci-Wieferich Primes - An Application of ABC Conjecture over Number Fields

abstract: In this paper, we define $X$-base Fibonacci-Wieferich prime which is a generalized Wieferich prime where $X$ is a finite set of algebraic numbers. We are going to show that there are infinitely many non-$X$-base Fibonacci-Wieferich primes assuming the $abc$-conjecture of Masser-Oesterl\'{e}-Szpiro for number fields. We also provide a new conjecture concerning the rank of free part of abelian group generated by all elements in $X$, and we will use the arithmetic point of view and geometric point of view to give heuristic.


title: The Invisible Hand of Dynamic Market Pricing

abstract: Walrasian prices, if they exist, have the property that one can assign every buyer some bundle in her demand set, such that the resulting assignment will maximize social welfare. Unfortunately, this assumes carefully breaking ties amongst different bundles in the buyer demand set. Presumably, the shopkeeper cleverly convinces the buyer to break ties in a manner consistent with maximizing social welfare. Lacking such a shopkeeper, if buyers arrive sequentially and simply choose some arbitrary bundle in their demand set, the social welfare may be arbitrarily bad. In the context of matching markets, we show how to compute dynamic prices, based upon the current inventory, that guarantee that social welfare is maximized. Such prices are set without knowing the identity of the next buyer to arrive. We also show that this is impossible in general (e.g., for coverage valuations), but consider other scenarios where this can be done. We further extend our results to Bayesian and bounded rationality models.


title: On the relation between a graph code and a graph state

abstract: A graph state and a graph code respectively are defined based on a mathematical simple graph. In this work, we examine a relation between a graph state and a graph code both obtained from the same graph, and show that a graph state is a superposition of logical qubits of the related graph code. By using the relation, we first discuss that a local complementation which has been used for a graph state can be useful for searching locally equivalent stabilizer codes, and second provide a method to find a stabilizer group of a graph code.


title: Hexagonal Rare-Earth Manganites as Promising Photovoltaics and Light Polarizers

abstract: Ferroelectric materials possess a spontaneous electric polarization and may be utilized in various technological applications ranging from non-volatile memories to solar cells and light polarizers. Recently, hexagonal rare-earth manganites, h-RMnO$_3$ (R is a rare-earth ion) have attracted considerable interest due to their intricate multiferroic properties and improper ferroelectricity characterized by a sizable remnant polarization and high Curie temperature. Here, we demonstrate that these compounds can serve as very efficient photovoltaic materials and, in addition, possess remarkable optical anisotropy properties. Using first-principles methods based on density-functional theory and considering h-TbMnO$_3$ as a representative manganite, we predict a strong light absorption of this material in the solar spectrum range, resulting in the maximum light-to-electricity energy conversion efficiency up to 33%. We also predict an extraordinary optical linear dichroism and linear birefringence properties of h-TbMnO$_3$ in a broad range of optical frequencies. These results uncover the unexplored potential of hexagonal rare-earth manganites to serve as photovoltaics in solar cells and as absorptive and birefringent light polarizers.


title: Optimal CCD readout by digital correlated double sampling

abstract: Digital correlated double sampling (DCDS), a readout technique for charge-coupled devices (CCD), is gaining popularity in astronomical applications. By using an oversampling ADC and a digital filter, a DCDS system can achieve a better performance than traditional analogue readout techniques at the expense of a more complex system analysis. Several attempts to analyse and optimize a DCDS system have been reported, but most of the work presented in the literature has been experimental. Some approximate analytical tools have been presented for independent parameters of the system, but the overall performance and trade-offs have not been yet modelled. Furthermore, there is disagreement among experimental results that cannot be explained by the analytical tools available. In this work, a theoretical analysis of a generic DCDS readout system is presented, including key aspects such as the signal conditioning stage, the ADC resolution, the sampling frequency and the digital filter implementation. By using a time-domain noise model, the effect of the digital filter is properly modelled as a discrete-time process, thus avoiding the imprecision of continuous-time approximations that have been used so far. As a result, an accurate, closed-form expression for the signal-to-noise ratio at the output of the readout system is reached. This expression can be easily optimized in order to meet a set of specifications for a given CCD, thus providing a systematic design methodology for an optimal readout system. Simulated results are presented to validate the theory, obtained with both time- and frequency-domain noise generation models for completeness.


title: Tree-Guided MCMC Inference for Normalized Random Measure Mixture Models

abstract: Normalized random measures (NRMs) provide a broad class of discrete random measures that are often used as priors for Bayesian nonparametric models. Dirichlet process is a well-known example of NRMs. Most of posterior inference methods for NRM mixture models rely on MCMC methods since they are easy to implement and their convergence is well studied. However, MCMC often suffers from slow convergence when the acceptance rate is low. Tree-based inference is an alternative deterministic posterior inference method, where Bayesian hierarchical clustering (BHC) or incremental Bayesian hierarchical clustering (IBHC) have been developed for DP or NRM mixture (NRMM) models, respectively. Although IBHC is a promising method for posterior inference for NRMM models due to its efficiency and applicability to online inference, its convergence is not guaranteed since it uses heuristics that simply selects the best solution after multiple trials are made. In this paper, we present a hybrid inference algorithm for NRMM models, which combines the merits of both MCMC and IBHC. Trees built by IBHC outlines partitions of data, which guides Metropolis-Hastings procedure to employ appropriate proposals. Inheriting the nature of MCMC, our tree-guided MCMC (tgMCMC) is guaranteed to converge, and enjoys the fast convergence thanks to the effective proposals guided by trees. Experiments on both synthetic and real-world datasets demonstrate the benefit of our method.


title: General de Finetti type theorems in noncommutative probability

abstract: We prove general de Finetti type theorems for classical and free independence. The de Finetti type theorems work for all non-easy quantum groups, which generalize a recent work of Banica, Curran and Speicher. We determine maximal distributional symmetries which means the corresponding de Finetti type theorem fails if a sequence of random variables satisfy more symmetry relations other than the maximal one. In addition, we define Boolean quantum semigroups in analogous to the easy quantum groups, by universal conditions on matrix coordinate generators and an orthogonal projection. Then, we show a general de Finetti type theorem for Boolean independence.


title: Systematic study of $\alpha$ decay half-lives for even-even nuclei within a two-potential approach

abstract: $\alpha$ decay is a common and important process for natural radioactivity of heavy and superheavy nuclei. The $\alpha$ decay half-lives for even-even nuclei from Z=62 to Z=118 are systematically researched based on the two-potential approach with a quasi-stationary state approximation. To describe the deviations between experimental half-lives and calculated results due to the nuclear shell structure, a hindrance factor related with $\alpha$ particle preformation probability is introduced. Our results can well reproduce the experimental data equally to the density-dependent cluster model and the generalized liquid drop model. We also study the isospin effect of nuclear potential in this work. Considering the isospin effect the calculated results improved about 7.3$\%$.


title: Why are deep nets reversible: A simple theory, with implications for training

abstract: Generative models for deep learning are promising both to improve understanding of the model, and yield training methods requiring fewer labeled samples. Recent works use generative model approaches to produce the deep net's input given the value of a hidden layer several levels above. However, there is no accompanying "proof of correctness" for the generative model, showing that the feedforward deep net is the correct inference method for recovering the hidden layer given the input. Furthermore, these models are complicated. The current paper takes a more theoretical tack. It presents a very simple generative model for RELU deep nets, with the following characteristics: (i) The generative model is just the reverse of the feedforward net: if the forward transformation at a layer is $A$ then the reverse transformation is $A^T$. (This can be seen as an explanation of the old weight tying idea for denoising autoencoders.) (ii) Its correctness can be proven under a clean theoretical assumption: the edge weights in real-life deep nets behave like random numbers. Under this assumption ---which is experimentally tested on real-life nets like AlexNet--- it is formally proved that feed forward net is a correct inference method for recovering the hidden layer. The generative model suggests a simple modification for training: use the generative model to produce synthetic data with labels and include it in the training set. Experiments are shown to support this theory of random-like deep nets; and that it helps the training.


title: Explaining the reportedly over-massive black holes in early-type galaxies with intermediate-scale discs

abstract: The classification "early-type" galaxy includes both elliptically- and lenticular-shaped galaxies. Theoretically, the spheroid-to-disc flux ratio of an early-type galaxy can assume any positive value, but in practice studies often consider only spheroid/disc decompositions in which the disc neatly dominates over the spheroid at large galaxy radii, creating an inner "bulge" as observed in most spiral galaxies. Here we show that decompositions in which the disc remains embedded within the spheroid, labelled by some as "unphysical", correctly reproduce both the photometric and kinematic properties of early-type galaxies with intermediate-scale discs. Intermediate-scale discs have often been confused with large-scale discs and incorrectly modelled as such; when this happens, the spheroid luminosity is considerably underestimated. This has recently led to some surprising conclusions, such as the claim that a number of galaxies with intermediate-scale discs (Mrk 1216, NGC 1277, NGC 1271, and NGC 1332) host a central black hole whose mass is abnormally large compared to expectations from the (underestimated) spheroid luminosity. We show that when these galaxies are correctly modelled, they no longer appear as extreme outliers in the (black hole mass)-(spheroid mass) diagram. This not only nullifies the need for invoking different evolutionary scenarios for these galaxies but it strengthens the significance of the observed (black hole mass)-(spheroid mass) correlation and confirms its importance as a fundamental ingredient for theoretical and semi-analytic models used to describe the coevolution of spheroids and their central supermassive black holes.


title: Axino dark matter and baryon number asymmetry production by the Q-ball decay in gauge mediation

abstract: We investigate the Q-ball decay into the axino dark matter in the gauge-mediated supersymmetry breaking. In our scenario, the Q ball decays mainly into nucleons and partially into axinos to account respectively for the baryon asymmetry and the dark matter of the universe. The Q ball decays well before the big bang nucleosynthesis so that it is not affected by the decay. We show the region of the parameters which realizes this scenario.


title: Time Independent Universal Computing with Spin Chains: Quantum Plinko Machine

abstract: We present a scheme for universal quantum computing using XY Heisenberg spin chains. Information is encoded into packets propagating down these chains, and they interact with each other to perform universal quantum computation. A circuit using g gate blocks on m qubits can be encoded into chains of length $O(g^{3+\delta} m^{3+\delta})$ for all $\delta >0$ with vanishingly small error.


title: A short proof of the converse to a theorem of Steinhaus

abstract: A result of H. Steinhaus states that any positive Lebesgue measurable set has a property that its difference set contains an open interval around the origin. Y. V. Mospan proved that this result is the characterization of absolutely continuous measure. In this note we give a short proof of it.


title: Minimalist coupled evolution model for stellar x-ray activity, rotation, mass loss, and magnetic field

abstract: Late-type main sequence stars exhibit an x-ray to bolometric flux ratio that depends on ${\tilde Ro}$, the ratio of rotation period to convective turnover time, as ${\tilde Ro}^{-\zeta}$ with $2\le \zeta \le 3$ for ${\tilde Ro} >0.13$, but saturates with $|\zeta| <0.2$ for ${\tilde Ro} < 0.13$. Saturated stars are younger than unsaturated stars and show a broader spread of rotation rates and x-ray activity. The unsaturated stars have magnetic fields and rotation speeds that scale roughly with the square root of their age, though possibly flattening for stars older than the sun. The connection between faster rotators, stronger fields, and higher activity has been established observationally, but a theory for the unified time-evolution of x-ray luminosity, rotation, magnetic field and mass loss that captures the above trends has been lacking. Here we derive a minimalist holistic framework for the time evolution of these quantities built from combining a Parker wind with new ingredients: (1) explicit sourcing of both the thermal energy launching the wind and the x-ray luminosity via dynamo produced magnetic fields; (2) explicit coupling of x-ray activity and mass loss saturation to dynamo saturation (via magnetic helicity build-up and convection eddy shredding); (3) use of coronal equilibrium to determine how magnetic energy is divided into wind and x-ray contributions. For solar-type stars younger than the sun, we infer conduction to be a subdominant power loss compared to x-rays and wind. For older stars, conduction is more important, possibly quenching the wind and reducing angular momentum loss. We focus on the time evolution for stars younger than the sun, highlighting what is possible for further generalizations. Overall, the approach shows promise toward a unified explanation of all of the aforementioned observational trends.


title: Learning Discriminative Representations for Semantic Cross Media Retrieval

abstract: Heterogeneous gap among different modalities emerges as one of the critical issues in modern AI problems. Unlike traditional uni-modal cases, where raw features are extracted and directly measured, the heterogeneous nature of cross modal tasks requires the intrinsic semantic representation to be compared in a unified framework. This paper studies the learning of different representations that can be retrieved across different modality contents. A novel approach for mining cross-modal representations is proposed by incorporating explicit linear semantic projecting in Hilbert space. The insight is that the discriminative structures of different modality data can be linearly represented in appropriate high dimension Hilbert spaces, where linear operations can be used to approximate nonlinear decisions in the original spaces. As a result, an efficient linear semantic down mapping is jointly learned for multimodal data, leading to a common space where they can be compared. The mechanism of "feature up-lifting and down-projecting" works seamlessly as a whole, which accomplishes crossmodal retrieval tasks very well. The proposed method, named as shared discriminative semantic representation learning (\textbf{SDSRL}), is tested on two public multimodal dataset for both within- and inter- modal retrieval. The experiments demonstrate that it outperforms several state-of-the-art methods in most scenarios.


title: Bayesian hypothesis testing for one bit compressed sensing with sensing matrix perturbation

abstract: This letter proposes a low-computational Bayesian algorithm for noisy sparse recovery in the context of one bit compressed sensing with sensing matrix perturbation. The proposed algorithm which is called BHT-MLE comprises a sparse support detector and an amplitude estimator. The support detector utilizes Bayesian hypothesis test, while the amplitude estimator uses an ML estimator which is obtained by solving a convex optimization problem. Simulation results show that BHT-MLE algorithm offers more reconstruction accuracy than that of an ML estimator (MLE) at a low computational cost.


title: Comparison of the analytical approximation formula and Newton's method for solving a class of nonlinear Black-Scholes parabolic equations

abstract: Market illiquidity, feedback effects, presence of transaction costs, risk from unprotected portfolio and other nonlinear effects in PDE based option pricing models can be described by solutions to the generalized Black-Scholes parabolic equation with a diffusion term nonlinearly depending on the option price itself. Different linearization techniques such as Newton's method and analytic asymptotic approximation formula are adopted and compared for a wide class of nonlinear Black-Scholes equations including, in particular, the market illiquidity model and the risk-adjusted pricing model. Accuracy and time complexity of both numerical methods are compared. Furthermore, market quotes data was used to calibrate model parameters.


title: Discovering Underlying Plans Based on Distributed Representations of Actions

abstract: Plan recognition aims to discover target plans (i.e., sequences of actions) behind observed actions, with history plan libraries or domain models in hand. Previous approaches either discover plans by maximally "matching" observed actions to plan libraries, assuming target plans are from plan libraries, or infer plans by executing domain models to best explain the observed actions, assuming complete domain models are available. In real world applications, however, target plans are often not from plan libraries and complete domain models are often not available, since building complete sets of plans and complete domain models are often difficult or expensive. In this paper we view plan libraries as corpora and learn vector representations of actions using the corpora; we then discover target plans based on the vector representations. Our approach is capable of discovering underlying plans that are not from plan libraries, without requiring domain models provided. We empirically demonstrate the effectiveness of our approach by comparing its performance to traditional plan recognition approaches in three planning domains.


title: Spin current and magnetization induced by circularly polarized synchrotron radiation in magnetically-doped topological insulator Bi$_{1.37}$V$_{0.03}$Sb$_{0.6}$Te$_2$Se

abstract: We propose a hole-induced mechanism of spin-polarized current generation by circularly polarized synchrotron radiation and corresponding induced magnetization in magnetically-doped topological insulators Bi$_{1.37}$V$_{0.03}$Sb$_{0.6}$Te$_2$Se. Considered spin-polarized current is generated due to the spin-dependent depopulation of the Dirac cone topological surface states at the Fermi level and subsequent compensation of the generated holes. We have found experimentally and theoretically a relation between the generated spin-polarized current and the shift of the electrochemical potential. The out-of-plane magnetization induced by circularly polarized synchrotron radiation and its inversion with switching the direction of circular polarization were experimentally shown and theoretically confirmed.


title: Gravitational acceleration and tidal effects in spherical-symmetric density profiles

abstract: Pure power-law density profiles, $\rho(r)\propto r^{b-3}$, are classified in connection with the following reference cases: (i) isodensity, $b=3$, $\rho=$ const; (ii) isogravity, $b=2$, $g=$ const; (iii) isothermal, $b=1$, $v=[GM(r)/r]^{1/2}=$ const; (iv) isomass, $b=0$, $M=$ const. A restricted number of different families of density profiles including, in addition, cored power-law, generalized power-law, polytropes, are studied in detail with regard to both one-component and two-component systems. Considerable effort is devoted to the existence of an extremum point (maximum absolute value) in the gravitational acceleration within the matter distribution. Predicted velocity curves are compared to the data inferred from observations. Tidal effects on an inner subsystem are investigated and an application is made to globular clusters within the Galaxy. To this aim, the tidal radius is defined by balancing the opposite gravitational forces from the Galaxy and the selected cluster on a special point of the cluster boundary, lying between related centres of mass. The position of 17 globular clusters with respect to the stability region, where the tidal radius exceeds the observed radius, is shown for assigned dark-to-visible mass ratios and density profiles, among those considered, which are currently used for the description of galaxies and/or dark matter haloes.


title: Constructive stability and stabilizability of positive linear discrete-time switching systems

abstract: We describe a new class of positive linear discrete-time switching systems for which the problems of stability or stabilizability can be resolved constructively. This class generalizes the class of systems with independently switching state vector components. The distinctive feature of this class is that, for any system from this class, its components or blocks can be arbitrarily connected in parallel or in series without loss of the `constructive resolvability' property. It is shown also that, for such systems, it is possible to build constructively the individual positive trajectories with the greatest or the lowest rate of convergence to the zero.