title: Optimal Subsampling Approaches for Large Sample Linear Regression

abstract: A significant hurdle for analyzing large sample data is the lack of effective statistical computing and inference methods. An emerging powerful approach for analyzing large sample data is subsampling, by which one takes a random subsample from the original full sample and uses it as a surrogate for subsequent computation and estimation. In this paper, we study subsampling methods under two scenarios: approximating the full sample ordinary least-square (OLS) estimator and estimating the coefficients in linear regression. We present two algorithms, weighted estimation algorithm and unweighted estimation algorithm, and analyze asymptotic behaviors of their resulting subsample estimators under general conditions. For the weighted estimation algorithm, we propose a criterion for selecting the optimal sampling probability by making use of the asymptotic results. On the basis of the criterion, we provide two novel subsampling methods, the optimal subsampling and the predictor- length subsampling methods. The predictor-length subsampling method is based on the L2 norm of predictors rather than leverage scores. Its computational cost is scalable. For unweighted estimation algorithm, we show that its resulting subsample estimator is not consistent to the full sample OLS estimator. However, it has better performance than the weighted estimation algorithm for estimating the coefficients. Simulation studies and a real data example are used to demonstrate the effectiveness of our proposed subsampling methods.


title: Models and Representations for Fractal Social Organizations

abstract: Our subject is oriented towards investigation of potential ways of societal organization, that allow for collective intelligent organization and management of resources. The main objective of such organizations is the exploration of the social energy from the existing societies. We conjecture that an organizational model that fulfills the mentioned requirements is the Fractal Social Organization (FSO). Our goal is to prove and verify the effectiveness of this model by performing various simulations using the NetLogo environment, a tool that allows agent-based rapid prototyping. We begin by simulating trivial real life activities that demonstrate the main properties of the core unit of the FSO, namely the SoC. Further, more complex scenarios involving various nested SoCs are simulated. Two main simulation models are presented, allowing us to obtain preliminary results using the FSO concepts as a potential solution. In the first simulation model we demonstrate that by the use of FSO properties the individuals may benefit by receiving more qualitative healthcare services, while in the second simulation model we show how it might be possible to improve the fall detection systems by the use of FSO mechanism.


title: Revealed Preference at Scale: Learning Personalized Preferences from Assortment Choices

abstract: We consider the problem of learning the preferences of a heterogeneous population by observing choices from an assortment of products, ads, or other offerings. Our observation model takes a form common in assortment planning applications: each arriving customer is offered an assortment consisting of a subset of all possible offerings; we observe only the assortment and the customer's single choice. In this paper we propose a mixture choice model with a natural underlying low-dimensional structure, and show how to estimate its parameters. In our model, the preferences of each customer or segment follow a separate parametric choice model, but the underlying structure of these parameters over all the models has low dimension. We show that a nuclear-norm regularized maximum likelihood estimator can learn the preferences of all customers using a number of observations much smaller than the number of item-customer combinations. This result shows the potential for structural assumptions to speed up learning and improve revenues in assortment planning and customization. We provide a specialized factored gradient descent algorithm and study the success of the approach empirically.


title: Non-abelian tensor square of finite-by-nilpotent groups

abstract: Let $G$ be a group. We denote by $\nu(G)$ an extension of the non-abelian tensor square $G \otimes G$ by $G \times G$. We prove that if $G$ is finite-by-nilpotent, then the non-abelian tensor square $G \otimes G$ is finite-by-nilpotent. Moreover, $\nu(G)$ is nilpotent-by-finite (Theorem A). Also we characterize BFC-groups in terms of $\nu(G)$ (Theorem B).


title: Face numbers of manifolds with boundary

abstract: We study face numbers of simplicial complexes that triangulate manifolds (or even normal pseudomanifolds) with boundary. Specifically, we establish a sharp lower bound on the number of interior edges of a simplicial normal pseudomanifold with boundary in terms of the number of interior vertices and relative Betti numbers. Moreover, for triangulations of manifolds with boundary all of whose vertex links have the weak Lefschetz property, we extend this result to sharp lower bounds on the number of higher-dimensional interior faces. Along the way we develop a version of Bagchi and Datta's $\sigma$- and $\mu$-numbers for the case of relative simplicial complexes and prove stronger versions of the above statements with the Betti numbers replaced by the $\mu$-numbers. Our results provide natural generalizations of known theorems and conjectures for closed manifolds and appear to be new even for the case of a ball.


title: Quasi-normal acoustic oscillations in the transonic Bondi flow

abstract: In recent work, we analyzed the dynamics of spherical and nonspherical acoustic perturbations of the Michel flow, describing the steady radial accretion of a relativistic perfect fluid into a nonrotating black hole. We showed that such perturbations undergo quasi-normal oscillations and computed the corresponding complex frequencies as a function of the black hole mass M and the radius r_c of the sonic horizon. It was found that when r_c is much larger than the Schwarzschild radius r_H = 2GM/c^2 of the black hole, these frequencies scale like the surface gravity of the analogue black hole associated with the acoustic metric. In this work, we analyze the Newtonian limit of the Michel solution and its acoustic perturbations. In this limit, the flow outside the sonic horizon reduces to the transonic Bondi flow, and the acoustic metric reduces to the one introduced by Unruh in the context of experimental black hole evaporation. We show that for the transonic Bondi flow, Unruh's acoustic metric describes an analogue black hole and compute the associated quasi-normal frequencies. We prove that they do indeed scale like the surface gravity of the acoustic black hole, thus providing an explanation for our previous results in the relativistic setting.


title: The effect of randomness for dependency map on the robustness of interdependent lattices

abstract: For interdependent networks with identity dependency map, percolation is exactly the same with that on a single network and follows a second-order phase transition, while for random dependency, percolation follows a first-order phase transition. In real networks, the dependency relations between networks are neither identical nor completely random. Thus in this paper, we study the influence of randomness for dependency maps on the robustness of interdependent lattice networks. We introduce approximate entropy($ApEn$) as the measure of randomness of the dependency maps. We find that there is critical $ApEn_c$ below which the percolation is continuous, but for larger $ApEn$, it is a first-order transition. With the increment of $ApEn$, the $p_c$ increases until $ApEn$ reaching ${ApEn}_c'$ and then remains almost constant. The time scale of the system shows rich properties as $ApEn$ increases. Our results uncover that randomness is one of the important factors that lead to cascading failures of spatially interdependent networks.


title: High-speed optical correlation-domain reflectometry without using acousto-optic modulator

abstract: To achieve a distributed reflectivity measurement along an optical fiber, we develop a simplified cost-effective configuration of optical correlation- (or coherence-) domain reflectometry based on a synthesized optical coherence function by sinusoidal modulation. By excluding conventional optical heterodyne detection (practically, without using an acousto-optic modulator) and by exploiting the foot of the Fresnel reflection spectrum, the electrical bandwidth required for signal processing is lowered down to several megahertz. We evaluate the basic system performance and demonstrate its high-speed operation (10 ms for one scan) by tracking a moving reflection point in real time.


title: MRL order, log-concavity and an application to peacocks

abstract: We provide an equivalent log-concavity condition to the mean residual life (MRL) ordering for real-valued processes. This result, combined with classical properties of total positivity of order 2, allows to exhibit new families of integrable processes which increase in the MRL order (MRL processes). Note that MRL processes with constant mean are peacocks to which the Az\'ema-Yor (Skorokhod embedding) algorithm yields an explicit associated martingale.


title: Normalized incomplete beta function: log-concavity in parameters and other properties

abstract: The normalized incomplete beta function can be defined either as cumulative distribution function of beta density or as the Gauss hypergeometric function with one of the upper parameters equal to unity. Logarithmic concavity/convexity of this function in parameters was established by Finner and Roters in 1997. Their proof is indirect and rather difficult; it is based on generalized reproductive property of certain more general distributions. These authors remark that these results "seems to be very hard to obtain by usual analytic methods". In the first part of this paper we provide such proof based on standard tools of analysis. In the second part we go one step further and investigate the sign of generalized Tur\'{a}n determinants formed by shifts of the normalized incomplete beta function. Under some additional restrictions we demonstrate that these coefficients are of the same sign. We further conjecture that such restrictions can be removed without altering the results. Our method of proof also leads to various companion results which may be of independent interest. In particular, we establish linearization formulas and two-sided bounds for the above mentioned Tur\'{a}n determinants. Further, we find two combinatorial style identities for finite sums which we believe to be new.


title: Detecting Community Structures in Hi-C Genomic Data

abstract: Community detection (CD) algorithms are applied to Hi-C data to discover new communities of loci in the 3D conformation of human and mouse DNA. We find that CD has some distinct advantages over pre-existing methods: (1) it is capable of finding a variable number of communities, (2) it can detect communities of DNA loci either adjacent or distant in the 1D sequence, and (3) it allows us to obtain a principled value of k, the number of communities present. Forcing k = 2, our method recovers earlier findings of Lieberman-Aiden, et al. (2009), but letting k be a parameter, our method obtains as optimal value k = 6, discovering new candidate communities. In addition to discovering large communities that partition entire chromosomes, we also show that CD can detect small-scale topologically associating domains (TADs) such as those found in Dixon, et al. (2012). CD thus provides a natural and flexible statistical framework for understanding the folding structure of DNA at multiple scales in Hi-C data.


title: Impact of ionization of ferrocene: EOES of alpha- and beta- electrons and the fingerprint orbital of ferrocenium

abstract: Ionization process of ferrocene (Fc) to produce ferrocenium cation (Fc+) has been debated as much as the eclipsed and staggered ferrocene conformers. The present quantum mechanical study reveals that removal of an electron does not apparently affect the geometry and symmetry of the cation, as the geometric changes are < 2% with respect to neutral Fc, but produce the fingerprint orbital 8a1 of Fc+. The excess orbital energy spectrum (EOES) of the \alpha- and \beta-electrons reveals that the electrons originated from the transition metal Fe in both core and valence shells experience significant energy changes in the cation with respect to the neutral ferrocene counterparts, indicating that the Fe-electrons correlate stronger than electrons from other atoms such as carbons in ferrocene. The EOES also exhibits that the orbital energies of the \alpha-electrons in ferrocenium change more significantly than the \beta-electrons after one \beta-electron being ionized with respect to ferrocene. The most significant changes upon ionization are dominated by the Fe-electrons in non-degenerate signature orbitals. That is, 4a1 which is dominated by the Fe 3s orbital, 3a2 by the Fe 3pz orbital and the 8a1 orbital (alpha only) by the Fe 3dz2 (one beta-electron is removed from this orbital). The \Delta DFT and \Delta SCF calculations yield the IP of 6.90 eV and 6.85 eV, respectively, in excellent agreement to the recent measurements of 6.9 \pn 0.1 eV, suggesting that significant relaxation energy exists whereas electron correlation energies are largely cancelled out. A further dual space analysis (DSA) identifies that the momentum profiles of such the singly occupied fingerprint 8a1 orbital indeed experiences most significant changes in all orbitals of Fc+, again indicates that the electron removal of ferrocene ionization is not from the highest occupied molecular orbital.


title: An application of multivariate total positivity to peacocks

abstract: We use multivariate total positivity theory to exhibit new families of peacocks. As the authors of \cite{HPRY}, our guiding example is the result of Carr-Ewald-Xiao \cite{CEX}. We shall introduce the notion of strong conditional monotonicity. This concept is strictly more restrictive than the conditional monotonicity as defined in \cite{HPRY} (see also \cite{Be}, \cite{BPR1} and \cite{ShS1}). There are many random vectors which are strongly conditionally monotone (SCM). Indeed, we shall prove that multivariate totally positive of order 2 (MTP$_2$) random vectors are SCM. As a consequence, stochastic processes with MTP$_2$ finite-dimensional marginals are SCM. This family includes processes with independent and log-concave increments, and one-dimensional diffusions which have absolutely continuous transition kernels.


title: Pole placement design for quantum systems via coherent observers

abstract: We previously extended Luenberger's approach for observer design to the quantum case, and developed a class of coherent observers which tracks linear quantum stochastic systems in the sense of mean values. In light of the fact that the Luenberger observer is commonly and successfully applied in classical control, it is interesting to investigate the role of coherent observers in quantum feedback. As the first step in exploring observer-based coherent control, in this paper we study pole-placement techniques for quantum systems using coherent observers, and in such a fashion, poles of a closed-loop quantum system can be relocated at any desired locations. In comparison to classical feedback control design incorporating the Luenberger observer, here direct coupling between a quantum plant and the observer-based controller are allowed to enable a greater degree of freedom for the design of controller parameters. A separation principle is presented, and we show how to design the observer and feedback independently to be consistent with the laws of quantum mechanics. The proposed scheme is applicable to coherent feedback control of quantum systems, especially when the transient dynamic response is of interest, and this issue is illustrated in an example.


title: Asymptotic convergence rates for coordinate descent in polyhedral sets

abstract: We consider a family of parallel methods for constrained optimization based on projected gradient descents along individual coordinate directions. In the case of polyhedral feasible sets, local convergence towards a regular solution occurs unconstrained in a reduced space, allowing for the computation of tight asymptotic convergence rates by sensitivity analysis, this even when global convergence rates are unavailable or too conservative. We derive linear asymptotic rates of convergence in polyhedra for variants of the coordinate descent approach, including cyclic, synchronous, and random modes of implementation. Our results find application in stochastic optimization, and with recently proposed optimization algorithms based on Taylor approximations of the Newton step.


title: Quenching of inclusive and tagged b-jets at the LHC

abstract: We present theoretical predictions for the nuclear-induced attenuation of the differential cross sections for inclusive and tagged b-jet production in heavy ion collisions at the LHC. We find that for inclusive b-jet production at high transverse momentum the mass effects are negligible, and that the attenuation is comparable to the one observed for light jets. On the other hand, for isolated-photon and B-meson-tagged b-jets the sample of events with heavy quarks produced at the early stages of the collision is greatly enhanced. Thus, these tagged b-jet final-states have a much more direct connection to the physics of b-quark energy loss. We present theoretical predictions for the quenching of such tagged b-jet events at the LHC and the QGP-induced modification of the related momentum imbalance and asymmetry. We demonstrate that these tagged processes can be used to accurately study the physics of heavy quark production and propagation in dense QCD matter.


title: Controllability function as time of motion. I

abstract: The admissible positional control problem for the canonical system with geometrical restrictions on the control is considered. The investigation is performed with the help of the controllability function method. We obtain controllability functions which are the time of motion from an arbitrary initial condition to the origin. We also reveal a set of controls which solves this problem.


title: Cosmic ray penetration in diffuse clouds

abstract: Cosmic rays are a fundamental source of ionization for molecular and diffuse clouds, influencing their chemical, thermal, and dynamical evolution. The amount of cosmic rays inside a cloud also determines the $\gamma$-ray flux produced by hadronic collisions between cosmic rays and cloud material. We study the spectrum of cosmic rays inside and outside of a diffuse cloud, by solving the stationary transport equation for cosmic rays including diffusion, advection and energy losses due to ionization of neutral hydrogen atoms. We found that the cosmic ray spectrum inside a diffuse cloud differs from the one in the interstellar medium (ISM) for energies smaller than $E_{br}\approx 100$ MeV, irrespective of the model details. Below $E_{br}$, the spectrum is harder (softer) than that in the ISM if the latter is a power law $\propto p^{-s}$ with $s$ larger (smaller) than $\sim0.42$. As a consequence also the ionization rate due to CRs is strongly affected. Assuming an average Galactic spectrum similar to the one inferred from AMS-2 and Voyager 1 data, we discuss the resulting ionization rate in a typical diffuse cloud.


title: Determining the open pit to underground transition: A new method

abstract: Many Ore Reserves are harvested by a combination of open pit and underground mining methods. In these cases there is often material that could be mined by either method, and a choice has to be made. The area containing this material is referred to as the transition zone. Deciding where to finish the open pit and start the underground is referred to as the transition problem and it has received some attention in the literature since the 1980s. In this paper we provide a review of existing approaches to the transition problem encompassing: graph-theory based optimisation employing an opportunity cost approach; heuristics and integer programming. We also present a novel opportunity cost approach, allowing it to take into account a crown pillar, and show how the new approach can be best applied through the unconventional application of an existing mine optimisation tool.


title: Long-lived universal resonant Bose gases

abstract: Quantum simulations based on near-resonance Bose gases are limited by their short lifetimes due to severe atom losses. In addition to this, the recently predicted thermodynamical instability adds another constraint on accessing the resonant Bose gases. In this article, we offer a potential solution by proposing long-lived resonant Bose gases in both two and three dimensions, where the conventional few-body losses are strongly suppressed. We show that the thermodynamical properties as well as the lifetimes of these strongly interacting systems are universal, and independent of short-range physics.


title: Null Geodesics in a Magnetically Charged Stringy Black Hole Spacetime

abstract: We study the geodesic motion of massless test particles in the background of a magnetic charged black hole spacetime in four dimensions in dilaton-Maxwell gravity. The behaviour of effective potential in view of the different values of black hole parameters is analysed in the equatorial plane. The possible orbits for null geodesics are also discussed in detail in view of the different values of the impact parameter. We have also calculated the frequency shift of photons in this spacetime. The results obtained are then compared with those for the electrically charged stringy black hole spacetime and the Schwarzschild black hole spacetime. It is observed that there exists no stable circular orbit outside the event horizon for massless test particles.


title: Dense free sets

abstract: We prove the existence of infinite dense free sets (in the usual topology) for set mappings on the reals, under reasonable assumptions.


title: Growth and electronic structure of epitaxial single-layer WS$_2$ on Au(111)

abstract: Large-area single-layer WS$_2$ is grown epitaxially on Au(111) using evaporation of W atoms in a low pressure H$_2$S atmosphere. It is characterized by means of scanning tunneling microscopy, low-energy electron diffraction and core-level spectroscopy. Its electronic band structure is determined by angle-resolved photoemission spectroscopy. The valence band maximum at $\bar{K}$ is found to be significantly higher than at $\bar{\Gamma}$. The observed dispersion around $\bar{K}$ is in good agreement with density functional theory calculations for a free-standing monolayer, whereas the bands at $\bar{\Gamma}$ are found to be hybridized with states originating from the Au substrate. Strong spin-orbit coupling leads to a large spin-splitting of the bands in the neighborhood of the $\bar{K}$ points, with a maximum splitting of 419(11)~meV. The valence band dispersion around $\bar{K}$ is found to be highly anisotropic with spin-branch dependent effective hole masses of $0.40(02)m_e$ and $0.57(09)m_e$ for the upper and lower split valence band, respectively. The large size of the spin-splitting and the low effective mass of the valence band maximum make single-layer WS$_2$ a promising alternative to the widely studied MoS$_2$ for applications in electronics, spintronics and valleytronics.


title: SZ effects in the Magneticum Pathfinder Simulation: Comparison with the Planck, SPT, and ACT results

abstract: We calculate the one-point probability density distribution functions (PDF) and the power spectra of the thermal and kinetic Sunyaev-Zeldovich (tSZ and kSZ) effects and the mean Compton Y parameter using the Magneticum Pathfinder simulations, state-of-the-art cosmological hydrodynamical simulations of a large cosmological volume of (896 Mpc/h)^3. These simulations follow in detail the thermal and chemical evolution of the intracluster medium as well as the evolution of super-massive black holes and their associated feedback processes. We construct full-sky maps of tSZ and kSZ from the light-cones out to z=0.17, and one realization of 8.8x8.8 degree wide, deep light-cone out to z=5.2. The local universe at z<0.027 is simulated by a constrained realisation. The tail of the one-point PDF of tSZ from the deep light-cone follows a power-law shape with an index of -3.2. Once convolved with the effective beam of Planck, it agrees with the PDF measured by Planck. The predicted tSZ power spectrum agrees with that of the Planck data at all multipoles up to l~1000, once the calculations are scaled to the Planck 2015 cosmological parameters with \Omega_m=0.308 and \sigma_8=0.8149. Consistent with the results in the literature, however, we continue to find the tSZ power spectrum at l=3000 that is significantly larger than that estimated from the high-resolution ground-based data. The simulation predicts the mean fluctuating Compton Y value of <Y>=1.18x10^{-6} for \Omega_m=0.272 and \sigma_8=0.809. Nearly half (~ 5x10^{-7}) of the signal comes from halos below a virial mass of 10^{13}M_\odot/h. Scaling this to the Planck 2015 parameters, we find <Y>=1.57x10^{-6}. The PDF and the power spectrum of kSZ from our simulation agree broadly with the previous work.


title: Dark matter implications of the WMAP-Planck Haze

abstract: Gamma rays and microwave observations of the Galactic Center and surrounding areas indicate the presence of anomalous emission, whose origin remains ambiguous. The possibility of dark matter (DM) annihilation explaining both signals through prompt emission at gamma-rays and secondary emission at microwave frequencies from interactions of high-energy electrons produced in annihilation with the Galactic magnetic fields has attracted much interest in recent years. We investigate the DM interpretation of the Galactic Center gamma-ray excess by searching for the associated synchrotron in the WMAP-Planck data. Considering various magnetic field and cosmic-ray propagation models, we predict the synchrotron emission due to DM annihilation in our Galaxy, and compare it with the WMAP-Planck data at 23-70GHz. In addition to standard microwave foregrounds, we separately model the microwave counterpart to the Fermi Bubbles and the signal due to DM, and use component separation techniques to extract the signal associated with each template from the total emission. We confirm the presence of the Haze at the level of 7% of the total sky intensity at 23GHz in our chosen region of interest, with a harder spectrum $I \sim \nu^{-0.8}$ than the synchrotron from regular cosmic-ray electrons. The data do not show a strong preference towards fitting the Haze by either the Bubbles or DM emission only. Inclusion of both components provides a better fit with a DM contribution to the Haze emission of 20% at 23GHz, however, due to significant uncertainties in foreground modeling, we do not consider this a clear detection of a DM signal. We set robust upper limits on the annihilation cross section by ignoring foregrounds, and also report best-fit DM annihilation parameters obtained from a complete template analysis. We conclude that the WMAP-Planck data are consistent with a DM interpretation of the gamma-ray excess.

categories: astro-ph.CO astro-ph.GA astro-ph.HE hep-ph


title: Are weak measurements really necessary for Leggett-Garg type measurements?

abstract: Leggett-Garg inequalities are an important milestone in our quest to bridge the classical-quantum divide. An experimental investigation of these inequalities requires the so called \emph{non-invasive measurements}(NIM). It has become popular to invoke weak measurements as the means of realising NIMs to very good approximation, because of their allegedly low disturbance of systems under measurement. In this note, this is shown to be a myth; it is shown, by simple estimates of errors, that for comparable levels of statistical errors, even the strong or projective measurements can be used. In fact, it is shown that resource-wise, strong measurements are even preferable.


title: Joint Channel and Queue Aware Scheduling for Wireless Links with Multiple Fading States

abstract: In this work, we address the delay optimal scheduling problem for wireless transmission with fixed modulation over multi-state fading channels. We propose a stochastic scheduling policy which schedules the source to transmit with probability jointly based on the buffer and channel states, with an average power constraint at the transmitter. Our objective is to minimize the average queueing delay by choosing the optimal transmission probabilities. Using Markov chain modeling, we formulate a power-constrained delay minimization problem, and then transform it into a Linear Programming (LP) one. By analyzing its property, we can derive the optimal threshold-based scheduling policy together with the corresponding transmission probabilities. Our theoretical analysis is corroborated by simulation results.


title: Superfluid density of a pseudogapped superconductor near SIT

abstract: We analyze critical behavior of superfluid density $\rho_s$ in strongly disordered superconductors near superconductor-insulator transition and compare it with the behavior of the spectral gap $\Delta$ for collective excitations. We show that in contrast to conventional superconductors, the superconductors with preformed pairs display unusual scaling relation $\rho_s \propto \Delta^2$ close to superconductor-insulator transition. This relation have been reported in very recent experiments.


title: Dark matter as a cancer hazard

abstract: We comment on the paper "Dark Matter collisions with the Human Body" by K. Freese and C. Savage (Phys. Lett. B 717, 25 (2012) [arXiv:1204.1339]) and describe a dark matter model for which the results of the previous paper do not quite apply. Within this mirror dark matter model, potentially hazardous objects, mirror micrometeorites, can exist and may lead to diseases triggered by multiple mutations, such as cancer, though with very low probability.


title: Multiscale modeling of microscale fiber reinforced composites with nano-engineered interphases

abstract: This study is focused on the mechanical properties and stress transfer behavior of multiscale composites containing nano- and micro-scale reinforcements. The distinctive feature of construction of this composite is such that the carbon nanostructures (CNS) are dispersed in the matrix around the continuous microscale fiber to modify microfiber-matrix interfacial adhesion. Such CNS are considered to be made of aligned CNTs (A-CNTs). Accordingly, multiscale models are developed for such hybrid composites. First, molecular dynamics simulations in conjunction with the Mori-Tanaka method are used to determine the effective elastic properties of nano-engineered interphase layer composed of CNS and epoxy. Subsequently, a micromechanical pull-out model for a continuous fiber multi-scale composite is developed, and stress transfer behavior is studied for different orientations of CNS considering their perfect and imperfect interfacial bonding conditions with the surrounding epoxy. Such interface condition was modeled using the linear spring layer model with a continuous traction but a displacement jump. The current pull-out model accounts for the radial as well as the axial deformations of different orthotropic constituent phases of the multiscale composite. The results from the developed pull-out model are compared with those of the finite element analyses and are found to be in good agreement. Our results reveal that the stress transfer characteristics of the multiscale composite are significantly improved by controlling the CNT morphology around the fiber, particularly, when they are aligned along the axial direction of the microscale fiber. The results also show that the CNS-epoxy interface weakening significantly influences the radial stress along the length of the microscale fiber.


title: A Tight SDP Relaxation for MIQCQP Problems in Power Systems Based on Disjunctive Programming

abstract: An optimization problem considering AC power flow constraints and integer decision variables can usually be posed as a mixed-integer quadratically constrained quadratic program (MIQCQP) problem. In this paper, first, a set of valid linear equalities are applied to strengthen the semidefinite program (SDP) relaxation of the MIQCQP problem without significantly increasing the problem dimension so that an enhanced mixed-integer SDP (MISDP) relaxation, which is a mixed-integer convex problem, is obtained. Then, the enhanced MISDP relaxation is reformulated as a disjunctive programming (DP) problem which is tighter than the former one, since the disjunctions are designed to capture the disjunctive nature of the terms in the rank-1 constraint about the integral variables. The DP relaxation is then equivalently convert-ed back into a MISDP problem the feasible set of whose continu-ous relaxation is the convex hull of feasible region of the DP prob-lem. Finally, globally optimal solution of the DP problem which is the tightest relaxation for the MIQCQP proposed in the paper is obtained by solving the resulting MISDP problem using a branch-and-bound (B&B) algorithm. Computational efficiency of the B&B algorithm is expected to be high since feasible set of the continuous relaxation of a MISDP sub-problem is the convex hull of that of the corresponding DP sub-problem. To further reduce the dimension of the resulting MISDP problem, a compact for-mulation of this problem is proposed considering the sparsity. An optimal placement problem of smart PV inverter in distribution systems integrated with high penetration of PV, which is an MIQCQP problem, is studied as an example. The proposed ap-proach is tested on an IEEE distribution system. The results show that it can effectively improve the tightness and feasibility of the SDP relaxation.


title: Fast Gaussian Process Regression for Big Data

abstract: Gaussian Processes are widely used for regression tasks. A known limitation in the application of Gaussian Processes to regression tasks is that the computation of the solution requires performing a matrix inversion. The solution also requires the storage of a large matrix in memory. These factors restrict the application of Gaussian Process regression to small and moderate size data sets. We present an algorithm that combines estimates from models developed using subsets of the data obtained in a manner similar to the bootstrap. The sample size is a critical parameter for this algorithm. Guidelines for reasonable choices of algorithm parameters, based on detailed experimental study, are provided. Various techniques have been proposed to scale Gaussian Processes to large scale regression tasks. The most appropriate choice depends on the problem context. The proposed method is most appropriate for problems where an additive model works well and the response depends on a small number of features. The minimax rate of convergence for such problems is attractive and we can build effective models with a small subset of the data. The Stochastic Variational Gaussian Process and the Sparse Gaussian Process are also appropriate choices for such problems. These methods pick a subset of data based on theoretical considerations. The proposed algorithm uses bagging and random sampling. Results from experiments conducted as part of this study indicate that the algorithm presented in this work can be as effective as these methods. Model stacking can be used to combine the model developed with the proposed method with models from other methods for large scale regression such as Gradient Boosted Trees. This can yield performance gains.


title: Extended Lyman alpha haloes around individual high-redshift galaxies revealed by MUSE

abstract: We report the detection of extended Ly alpha emission around individual star-forming galaxies at redshifts z = 3-6 in an ultradeep exposure of the Hubble Deep Field South obtained with MUSE on the ESO-VLT. The data reach a limiting surface brightness (1sigma) of ~1 x 10^-19 erg s^-1 cm^-2 arcsec^-2 in azimuthally averaged radial profiles, an order of magnitude improvement over previous narrowband imaging. Our sample consists of 26 spectroscopically confirmed Ly alpha-emitting, but mostly continuum-faint (m_AB >~ 27) galaxies. In most objects the Ly alpha emission is considerably more extended than the UV continuum light. While 5 of the faintest galaxies in the sample show no significantly detected Ly alpha haloes, the derived upper limits suggest that this is just due to insufficient S/N. Ly alpha haloes therefore appear to be (nearly) ubiquitous even for low-mass (~10^8-10^9 M_sun) star-forming galaxies at z>3. We decompose the Ly alpha emission of each object into a compact `continuum-like' and an extended halo component, and infer sizes and luminosities of the haloes. The extended Ly alpha emission approximately follows an exponential surface brightness distribution with a scale length of a few kpc. While these haloes are thus quite modest in terms of their absolute sizes, they are larger by a factor of 5-15 than the corresponding rest-frame UV continuum sources as seen by HST. They are also much more extended, by a factor ~5, than Ly alpha haloes around low-redshift star-forming galaxies. Between ~40% and >90% of the observed Ly alpha flux comes from the extended halo component, with no obvious correlation of this fraction with either the absolute or the relative size of the Ly alpha halo. Our observations provide direct insights into the spatial distribution of at least partly neutral gas residing in the circumgalactic medium of low to intermediate mass galaxies at z > 3.


title: Distributed PROMPT-LTL Synthesis

abstract: We consider the synthesis of distributed implementations for specifications in Prompt Linear Temporal Logic (PROMPT-LTL), which extends LTL by temporal operators equipped with parameters that bound their scope. For single process synthesis it is well-established that such parametric extensions do not increase worst-case complexities. For synchronous systems, we show that, despite being more powerful, the distributed realizability problem for PROMPT-LTL is not harder than its LTL counterpart. For asynchronous systems we have to consider an assume-guarantee synthesis problem, as we have to express scheduling assumptions. As asynchronous distributed synthesis is already undecidable for LTL, we give a semi-decision procedure for the PROMPT-LTL assume-guarantee synthesis problem based on bounded synthesis.


title: Modified TOV in gravity's rainbow: properties of neutron stars and dynamical stability conditions

abstract: In this paper, we consider a spherical symmetric metric to extract the hydrostatic equilibrium equation of stars in $(3+1)-$dimensional gravity's rainbow in the presence of cosmological constant. Then, we generalize the hydrostatic equilibrium equation to $d$-dimensions and obtain the hydrostatic equilibrium equation for this gravity. Also, we obtain the maximum mass of neutron star using the modern equations of state of neutron star matter derived from the microscopic calculations. It is notable that, in this paper, we consider the effects of rainbow functions on the diagrams related to the mass-central mass density ($M$-$\rho _{c}$) relation and also the mass-radius ($M$-$R$) relation of neutron star. We also study the effects of rainbow functions on the other properties of neutron star such as the Schwarzschild radius, average density, strength of gravity and gravitational redshift. Then, we apply the cosmological constant to this theory to obtain the diagrams of $M$-$\rho _{c}$ (or $M$-$R$) and other properties of these stars. Next, we investigate the dynamical stability condition for these stars in gravity's rainbow and show that these stars have dynamical stability. We also obtain a relation between mass of neutron stars and Planck mass. In addition, we compare obtained results of this theory with the observational data.


title: Detection of CI line emission from the detached CO shell of the AGB star R Sculptoris

abstract: Stars on the asymptotic giant branch (AGB) lose substantial amounts of matter, to the extent that they are important for the chemical evolution of, and dust production in, the universe. The mass loss is believed to increase gradually with age on the AGB, but it may also occur in the form of bursts, possibly related to the thermal pulsing phenomenon. Detached, geometrically thin, CO shells around carbon stars are good signposts of brief and intense mass ejection. We aim to put further constraints on the physical properties of detached CO shells around AGB stars. The photodissociation of CO and other carbon-bearing species in the shells leads to the possibility of detecting lines from neutral carbon. We have therefore searched for the CI($^3P_1-\,^3P_0$) line at 492 GHz towards two carbon stars, S Sct and R Scl, with detached CO shells of different ages, about 8000 and 2300 years, respectively. The CI($^3P_1-\,^3P_0$) line was detected towards R Scl. The line intensity is dominated by emission from the detached shell. The detection is at a level consistent with the neutral carbon coming from the full photodissociation of all species except CO, and with only limited photoionisation of carbon. The best fit to the observed $^{12}$CO and $^{13}$CO line intensities, assuming a homogeneous shell, is obtained for a shell mass of about 0.002 $M_\odot$, a temperature of about 100 K, and a CO abundance with respect to H$_2$ of 10$^{-3}$. The estimated CI/CO abundance ratio is about 0.3 for the best-fit model. However, a number of arguments point in the direction of a clumpy medium, and a viable interpretation of the data within such a context is provided.


title: Algebraic K-theory of quotient stacks

abstract: We prove some fundamental results like localization, excision, Nisnevich descent and the Mayer-Vietoris property for equivariant regular blow-up for the equivariant K-theory of schemes with an affine group scheme action. We also show that the equivariant K-theory with finite coefficients is invariant under equivariant vector bundle morphisms. We show that the equivariant homotopy K-theory is invariant under equivariant vector bundle morphisms and satisfies all the above properties, including nil-invariance.


title: Nonequilibrium work energy relation for non-Hamiltonian dynamics

abstract: Recent years have witnessed major advances in our understanding of nonequilibrium processes. The Jarzynski equality, for example, provides a link between equilibrium free energy differences and finite-time, nonequilibrium dynamics. We propose a generalization of this relation to non-Hamiltonian dynamics, relevant for active matter systems, continuous feedback, and computer simulation. Surprisingly, this relation allows us to calculate the free energy difference between the desired initial and final equilibrium states using arbitrary dynamics. As a practical matter, this dissociation between the dynamics and the initial and final states promises to facilitate a range of techniques for free energy estimation in a single, universal expression.


title: Iterated limits for aggregation of randomized INAR(1) processes with Poisson innovations

abstract: We discuss joint temporal and contemporaneous aggregation of $N$ independent copies of strictly stationary INteger-valued AutoRegressive processes of order 1 (INAR(1)) with random coefficient $\alpha\in(0,1)$ and with idiosyncratic Poisson innovations. Assuming that $\alpha$ has a density function of the form $\psi(x)(1 - x)^\beta$, $x\in(0,1)$, with $\lim_{x\uparrow 1}\psi(x) = \psi_1 \in(0,\infty)$, different limits of appropriately centered and scaled aggregated partial sums are shown to exist for $\beta\in(-1,0)$, $\beta = 0$, $\beta\in(0,1)$ or $\beta\in(1,\infty)$, when taking first the limit as $N\to\infty$ and then the time scale $n\to\infty$, or vice versa. In fact, we give a partial solution to an open problem of Pilipauskaite and Surgailis (2014) by replacing the random-coefficient AR(1) process with a certain randomized INAR(1) process.


title: Direct Urca processes involving singlet proton superfluidity in neutron star cooling

abstract: A detailed description of the baryon direct Urca processes A: $n\rightarrow p+e+\bar{\nu}_{e}$, B: $\Lambda\rightarrow p+e+\bar{\nu}_{e}$, C: $\Xi^{-}\rightarrow\Lambda+e+\bar{\nu}_{e}$ related to the neutron star cooling is given in the relativistic mean field approximation.The contributions of the reactions B and C on the neutrino luminosity are calculated by means of the relativistic expressions of the neutrino energy losses.Our results show that the total neutrino luminosities of the reactions A, B, C within the mass range 1.603-2.067$M_{\odot}$ (1.515-1.840$M_{\odot}$ for TM1 model) for GM1 model are larger than the corresponding values for neutron stars in npe$\mu$ matter. Although the hyperon direct Urca processes B and C reduce the neutrino emissivity of the reaction A, it illustrates the reactions B and C still make the total neutrino luminosity enhancement in the above mentioned areas.Furthermore, when we only consider the $^{1}S_{0}$ proton superfluidity in neutron star cooling, we find that although the neutrino emissivity of the reactions A and B is suppressed with the appearance of $^{1}S_{0}$ proton superfluidity, the total contribution of the reactions A, B, C can still quicken a massive neutron star cooling.These results could be used to help prove appearing hyperons in PSR J1614-2230 and J0348+0432 from neutron star cooling perspective.


title: Deconfinement Phase Transition with External Magnetic Field in Friedberg-Lee Model

abstract: The deconfinement phase transition with external magnetic field is investigated in the Friedberg-Lee model. In the frame of functional renormalization group, we extend the often used potential expansion method for continuous phase transitions to the first-order phase transition in the model. By solving the flow equations we find that, the magnetic field displays a catalysis effect and it becomes more difficult to break through the confinement in hot and dense medium.


title: Polariton resonances for ultra-strong coupling cavity optomechanics in GaAs/AlAs multiple quantum wells

abstract: Polariton-mediated light-sound interaction is investigated through resonant Brillouin experiments in GaAs/AlAs multiple-quantum wells. Photoelastic coupling enhancement at exciton-polariton resonance reaches 105 at 30 K as compared to a typical bulk solid room temperature transparency value. When applied to GaAs based cavity optomechanical nanodevices, this result opens the path to huge displacement sensitivities and to novel ultrastrong-coupling cavity phenomena with optomechanical couplings g0 in the range of 100 THz.


title: Identifying Biochemical Reaction Networks From Heterogeneous Datasets

abstract: In this paper, we propose a new method to identify biochemical reaction networks (i.e. both reactions and kinetic parameters) from heterogeneous datasets. Such datasets can contain (a) data from several replicates of an experiment performed on a biological system; (b) data measured from a biochemical network subjected to different experimental conditions, for example, changes/perturbations in biological inductions, temperature, gene knock-out, gene over-expression, etc. Simultaneous integration of various datasets to perform system identification has the potential to avoid non-identifiability issues typically arising when only single datasets are used.


title: Phase space factors and half-life predictions for Majoron emitting $\beta^-\beta^-$ decay

abstract: A complete calculation of phase space factors (PSF) for Majoron emitting $0\nu\beta^-\beta^-$ decay modes is presented. The calculation makes use of exact Dirac wave functions with finite nuclear size and electron screening and includes life-times, single electron spectra, summed electron spectra, and angular electron correlations. Combining these results with recent interacting boson nuclear matrix elements (NME) we make half-life predictions for the the ordinary Majoron decay (spectral index $n$=1). Furthermore, comparing theoretical predictions with the obtained experimental lower bounds for this decay mode we are able to set limits on the effective Majoron-neutrino coupling constant $\langle g_{ee}^M\rangle$.


title: Decoupling with random diagonal unitaries

abstract: We investigate decoupling, one of the most important primitives in quantum Shannon theory, by replacing the uniformly distributed random unitaries commonly used to achieve the protocol, with repeated applications of random unitaries diagonal in the Pauli-$Z$ and -$X$ bases. This strategy was recently shown to achieve an approximate unitary $2$-design after a number of repetitions of the process, which implies that the strategy gradually achieves decoupling. Here, we prove that even fewer repetitions of the process achieve decoupling at the same rate as that with the uniform ones, showing that rather imprecise approximations of unitary $2$-designs are sufficient for decoupling. We also briefly discuss efficient implementations of them and implications of our decoupling theorem to coherent state merging and relative thermalisation.


title: The Cotton tensor and Chern-Simons invariants in dimension $3$: an introduction

abstract: Chern-Simons invariants of closed oriented Riemannian $3$-manifolds are introduced and studied from the basics. Their first-order variation is the Cotton tensor. The properties of the Cotton tensor: symmetry, conformal covariance, trace- and divergence-freedom, are recovered as corollaries of the Chern-Simons invariant. We prove that the Cotton tensor is the obstruction to local conformal flatness in dimension $3$. Finally we discuss the link between Chern-Simons invariants, the eta invariant, and the central value of the Selberg zeta function. This is a survey paper without original results, based on lecture notes for a doctoral course at the University of Bucharest, March 2014.


title: Neutrinoless double positron decay and positron emitting electron capture in the interacting boson model

abstract: Neutrinoless double-$\beta$ decay is of fundamental importance for determining the neutrino mass. Although double electron ($\beta^-\beta^-$) decay is the most promising mode, in very recent years interest in double positron ($\beta^+\beta^+$) decay, positron emitting electron capture ($EC\beta^+$), and double electron capture ($ECEC$) has been renewed. We present here results of a calculation of nuclear matrix elements for neutrinoless double-$\beta^+$ decay and positron emitting electron capture within the framework of the microscopic interacting boson model (IBM-2) for $^{58}$Ni, $^{64}$Zn, $^{78}$Kr, $^{96}$Ru, $^{106}$Cd, $^{124}$Xe, $^{130}$Ba, and $^{136}$Ce decay. By combining these with a calculation of phase space factors we calculate expected half-lives.


title: Response to a local quench of a system near many body localization transition

abstract: We consider a one dimensional spin $1/2$ chain with Heisenberg interaction in a disordered parallel magnetic field. This system is known to exhibit the many body localization (MBL) transition at critical strength of disorder. We analyze the response of the chain when additional perpendicular magnetic field is applied to an individual spin and propose a method for accurate determination of the mobility edge via local spin measurements. We further demonstrate that the exponential decrease of the spin response with the distance between perturbed spin and measured spin can be used to characterize the localization length in the MBL phase.


title: Phonon effects on x-ray absorption and nuclear magnetic resonance spectroscopies

abstract: In material sciences, spectroscopic approaches combining ab initio calculations with experiments are commonly used to accurately analyze the experimental spectral data. Most state-of-the-art first-principle calculations are usually performed assuming an equilibrium static lattice. Yet, nuclear motion affects spectra even when reduced to the zero-point motion at 0 K. We propose a framework based on Density-Functional Theory that includes quantum thermal fluctuations in theoretical X- ray Absorption Near-Edge Structure (XANES) and solid-state Nuclear Magnetic Resonance (NMR) spectroscopies and allows to well describe temperature effects observed experimentally. Within the Born-Oppenheimer and quasi-harmonic approximations, we incorporate the nuclear motion by generating several non-equilibrium configurations from the dynamical matrix. The averaged calculated XANES and NMR spectral data have been compared to experiments in MgO, proof-of-principle compound. The good agreement obtained between experiments and calculations validates the developed approach, which suggests that calculating the XANES spectra at finite temperature by averaging individual non-equilibrium configurations is a suitable approximation. This study high- lights the relevance of phonon renormalization and the relative contributions of thermal expansion and nuclear dynamics on NMR and XANES spectra on a wide range of temperatures.


title: Evolving Social Networks via Friend Recommendations

abstract: A social network grows over a period of time with the formation of new connections and relations. In recent years we have witnessed a massive growth of online social networks like Facebook, Twitter etc. So it has become a problem of extreme importance to know the destiny of these networks. Thus predicting the evolution of a social network is a question of extreme importance. A good model for evolution of a social network can help in understanding the properties responsible for the changes occurring in a network structure. In this paper we propose such a model for evolution of social networks. We model the social network as an undirected graph where nodes represent people and edges represent the friendship between them. We define the evolution process as a set of rules which resembles very closely to how a social network grows in real life. We simulate the evolution process and show, how starting from an initial network, a network evolves using this model. We also discuss how our model can be used to model various complex social networks other than online social networks like political networks, various organizations etc..


title: The role of water in the behavior of wood

abstract: Wood, due to its biological origin, has the capacity to interact with water. Sorption/desorption of moisture is accompanied with swelling/shrinkage and softening/hardening of its stiffness. The correct prediction of the behavior of wood components undergoing environmental loading requires that the moisture behavior and mechanical behavior of wood are considered in a coupled manner. We propose a comprehensive framework using a fully coupled poromechanical approach, where its multiscale implementation provides the capacity to take into account, directly, the exact geometry of the wood cellular structure, using computational homogenization. A hierarchical model is used to take into account the subcellular composite-like organization of the material. Such advanced modeling requires high resolution experimental data for the appropriate determination of inputs and for its validation. High-resolution x-ray tomography, digital image correlation, and neutron imaging are presented as valuable methods to provide the required information.


title: Perspectives of optical lattices with state-dependent tunneling in approaching quantum magnetism in the presence of the external harmonic trapping potential

abstract: We study theoretically potential advantages of two-component mixtures in optical lattices with state-dependent tunneling for approaching long-range-order phases and detecting easy-axis antiferromagnetic correlations. While we do not find additional advantages of mixtures with large hopping imbalance for approaching quantum magnetism in a harmonic trap, it is shown that a nonzero difference in hopping amplitudes remains highly important for a proper symmetry breaking in the pseudospin space for the single-site-resolution imaging and can be advantageously used for a significant increase of the signal-to-noise ratio and thus detecting long-range easy-axis antiferromagnetic correlations in the corresponding experiments.


title: Numerical study of variable lung ventilation strategies

abstract: Mechanical ventilation is used for patients with a variety of lung diseases. Traditionally, ventilators have been designed to monotonously deliver equal sized breaths. While it may seem intuitive that lungs may benefit from unvarying and stable ventilation pressure strategy, recently it has been reported that variable lung ventilation is advantageous. In this study, we analyze the mean tidal volume in response to different `variable ventilation pressure' strategies. We found that uniformly distributed variability in pressure gives the best tidal volume as compared to that of normal, scale- free, log normal and linear distributions.


title: Interacting galaxies in the nearby Universe: only moderate increase of star formation

abstract: We investigate the influence of interactions on the star formation by studying a sample of almost 1500 of the nearest galaxies, all within a distance of ~45 Mpc. We define the massive star formation rate (SFR), as measured from far-IR emission, and the specific star formation rate (SSFR), which is the former quantity normalised by the stellar mass of the galaxy, and explore their distribution with morphological type and with stellar mass. We then calculate the relative enhancement of these quantities for each galaxy by normalising them by the median SFR and SSFR values of individual control populations of similar non-interacting galaxies. We find that both SFR and SSFR are enhanced in interacting galaxies, and more so as the degree of interaction is higher. The increase is, however, moderate, reaching a maximum of a factor of 1.9 for the highest degree of interaction (mergers). The SFR and SSFR are enhanced statistically in the population, but in many individual interacting galaxies they are not enhanced at all. We discuss how those galaxies with the largest SFR and/or SSFR enhancement can be defined as starbursts. This study is based on a representative sample of nearby galaxies, including many low-mass and dwarf/irregular galaxies, and we argue that it should be used to place constraints on studies based on samples of galaxies at larger distances, beyond the local Universe.


title: Minimal control power of the controlled teleportation

abstract: We generalize the control power of a perfect controlled teleportation of an entangled three-qubit pure state, suggested by Li and Ghose [Phys. Rev. A {\bf 90}, 052305 (2014)], to the control power of a general controlled teleportation of a multiqubit pure state. Thus, we define the minimal control power, and calculate the values of the minimal control power for a class of general three-qubit Greenberger-Horne-Zeilinger (GHZ) states and the three-qubit W class whose states have zero three-tangles. Moreover, we show that the standard three-qubit GHZ state and the standard three-qubit W state have the maximal values of the minimal control power for the two classes, respectively. This means that the minimal control power can be interpreted as not only an operational quantity of a three-qubit quantum communication but also a degree of three-qubit entanglement. In addition, we calculate the values of the minimal control power for general n-qubit GHZ states and the n-qubit W-type states.


title: Density Induced Phases in Active Nematic

abstract: We introduce a minimal model for a collection of self-propelled apolar active particles, also called as `active nematic', on a two-dimensional substrate and study the order-disorder transition with the variation of density. The particles interact with their neighbours within the framework of the Lebwohl-Lasher model and move asymmetrically, along their orientation, to unoccupied nearest neighbour lattice sites. At a density lower than the equilibrium isotropic-nematic transition density, the active nematic shows a first order transition from the isotropic state to a banded state. The banded state extends over a range of density, and the scalar order parameter of the system shows a plateau like behaviour, similar to that of the magnetic systems. In the large density limit the active nematic shows a bistable behaviour between a homogeneous ordered state with global ordering and an inhomogeneous mixed state with local ordering. The study of the above phases with density variation is scant and gives significant insight of complex behaviours of many biological systems.


title: Computing the Kummer function U(a,b,z) for small values of the arguments

abstract: We describe methods for computing the Kummer function $U(a,b,z)$ for small values of $z$, with special attention to small values of $b$. For these values of $b$ the connection formula that represents $U(a,b,z)$ as a linear combination of two ${}_1F_1$-functions needs a limiting procedure. We use the power series of the ${}_1F_1$-functions and consider the terms for which this limiting procedure is needed. We give recursion relations for higher terms in the expansion, and we consider the derivative $U^\prime(a,b,z)$ as well. We also discuss the performance for small $\vert z\vert$ of an asymptotic approximation of the Kummer function in terms of modified Bessel functions.


title: Weak Infeasibility in Second Order Cone Programming

abstract: The objective of this work is to study weak infeasibility in second order cone programming. For this purpose, we consider a relaxation sequence of feasibility problems that mostly preserve the feasibility status of the original problem. This is used to show that for a given weakly infeasible problem at most $m$ directions are needed to approach the cone, where $m$ is the number of Lorentz cones. We also tackle a closely related question and show that given a bounded optimization problem satisfying Slater's condition, we may transform it into another problem that has the same optimal value but it is ensured to attain it. From solutions to the new problem, we discuss how to obtain solution to the original problem which are arbitrarily close to optimality. Finally, we discuss how to obtain finite certificate of weak infeasibility by combining our own techniques with facial reduction. The analysis is similar in spirit to previous work by the authors on SDPs, but a different approach is required to obtain tighter bounds.


title: A Review of Proximity Effect Correction in Electron-beam Lithography

abstract: I review the work of proximity effect correction (PEC) in electron-beam (e-beam) lithography with emphasis on dose modification and shape modification PEC techniques.


title: Mannheim Curves in the three-dimensional Sphere

abstract: Mannheim curves are defined for immersed curves in 3-dimensional sphere S^3 . The definition is given by considering the geodesics of S^3. First, two special geodesics, called principal normal geodesic and binormal geodesic, of S^3 are defined by using Frenet vectors of a curve immersed in S^3. Later, the curve alpha is called a Mannheim curve if there exits another curve beta in S^3 such that the principal normal geodesics of beta coincide with the binormal geodesics of S^3 . It is obtained that if alpha and beta form a Mannheim pair then there exist a constant lambda that is not equal 0 and a non-constant function Mu such that Lambda.(kappa_alpha)+M(Tau_alpha)=1 where kappa_alpha, Tau_alpha are the curvatures of alpha. Moreover, the relation between a Mannheim curve immersed in S^3 and a generalized Mannheim curve in E^4 is obtained and a table containing comparison of Bertrand and Mannheim curves in S^3 is introduced.


title: The GALEX Ultraviolet Virgo Cluster Survey (GUViCS). VI: The UV luminosity function of the Virgo cluster and its surrounding regions

abstract: We use the GALEX data of the GUViCS survey to construct the NUV luminosity function of the Virgo cluster over ~ 300 deg.2, an area covering the cluster and its surrounding regions up to ~ 1.8 virial radii. The NUV luminosity function is also determined for galaxies of different morphological type and NUV-i colour, and for the different substructures within the cluster. These luminosity functions are robust vs. statistical corrections since based on a sample of 833 galaxies mainly identified as cluster members with spectroscopic redshift (808) or high-quality optical scaling relations (10). We fit these luminosity functions with a Schechter function, and compare the fitted parameters with those determined for other nearby clusters and for the field. The faint end slope of the Virgo NUV luminosity function (alpha = -1.19), here sampled down to ~ NUV = -11.5 mag, is significantly flatter than the one measured in other nearby clusters and similar to the field one. Similarly M* = -17.56 is one-to-two magnitudes fainter than measured in Coma, A1367, the Shapley supercluster, and the field. These differences seem due to the quite uncertain statistical corrections and the small range in absolute magnitude sampled in these clusters. We do not observe strong systematic differences in the overall NUV luminosity function of the core of the cluster with respect to that of its periphery. We notice, however, that the relative contribution of red-to-blue galaxies at the faint end is inverted, with red quiescent objects dominating the core of the cluster and star forming galaxies dominating beyond one virial radius. These observational evidences are discussed in the framework of galaxy evolution in dense environments.


title: Generalized Emphatic Temporal Difference Learning: Bias-Variance Analysis

abstract: We consider the off-policy evaluation problem in Markov decision processes with function approximation. We propose a generalization of the recently introduced \emph{emphatic temporal differences} (ETD) algorithm \citep{SuttonMW15}, which encompasses the original ETD($\lambda$), as well as several other off-policy evaluation algorithms as special cases. We call this framework \ETD, where our introduced parameter $\beta$ controls the decay rate of an importance-sampling term. We study conditions under which the projected fixed-point equation underlying \ETD\ involves a contraction operator, allowing us to present the first asymptotic error bounds (bias) for \ETD. Our results show that the original ETD algorithm always involves a contraction operator, and its bias is bounded. Moreover, by controlling $\beta$, our proposed generalization allows trading-off bias for variance reduction, thereby achieving a lower total error.


title: Taming the ReLU with Parallel Dither in a Deep Neural Network

abstract: Rectified Linear Units (ReLU) seem to have displaced traditional 'smooth' nonlinearities as activation-function-du-jour in many - but not all - deep neural network (DNN) applications. However, nobody seems to know why. In this article, we argue that ReLU are useful because they are ideal demodulators - this helps them perform fast abstract learning. However, this fast learning comes at the expense of serious nonlinear distortion products - decoy features. We show that Parallel Dither acts to suppress the decoy features, preventing overfitting and leaving the true features cleanly demodulated for rapid, reliable learning.


title: Wavy way to the Kerr metric and the quantum nature of its ring singularity

abstract: From inherent non-linearity two gravitational waves, unless they are unidirectional, fail to satisfy a superposition law. They collide to develop a new spacetime carrying the imprints of the incoming waves. Same behaviour is valid also for any massless lightlike field. As a result of the violent collision process either a naked singularity or a Cauchy horizon (CH) develops. It was shown by Chandrasekhar and Xanthopoulos (CX) that a particular class of colliding gravitational waves (CGW) spacetime is locally isometric to the Kerr metric for rotating black holes. This relation came to be known as the CX duality. Such a duality can be exploited as an alternative derivation for the Kerr metric as we do herein. Not each case gives rise to a CH but those which do are transient to a black hole state provided stability requirements are met. These classical considerations can be borrowed to shed light on black hole formation in high energy collisions. Their questionable stability and many other sophisticated agenda, we admit that await for a full - fledged quantum gravity. Yet, to add an element of novelty, a quantum probe is sent in the plane $\theta =\pi /2$ to the naked ring singularity of Kerr which develops for the overspinning case $(a>M)$ to test it from a quantum picture. We show that the spatial operator of the reduced Klein-Gordon equation has a unique self-adjoint extension. As a result, the classical Kerr`s ring singularity is healed and becomes quantum regular. Our poetic message of the paper is summarized as Let there be light that collide with might to disperse the night and create holes that are white


title: A theory of Galois descent for finite inseparable extensions

abstract: We present a generalization of Galois descent to finite modular normal field extension $L/K$, using the Heerma-Galois group $Aut(L[\bar{X}]/K[\bar{X}])$ where $L[\bar{X}]=L[X]/(X^{p^e})$ and $e$ is the exponent of $L$ over $K$.


title: Performance Analysis and Prediction in Educational Data Mining: A Research Travelogue

abstract: In this era of computerization, education has also revamped itself and is not limited to old lecture method. The regular quest is on to find out new ways to make it more effective and efficient for students. Nowadays, lots of data is collected in educational databases, but it remains unutilized. In order to get required benefits from such a big data, powerful tools are required. Data mining is an emerging powerful tool for analysis and prediction. It is successfully applied in the area of fraud detection, advertising, marketing, loan assessment and prediction. But, it is in nascent stage in the field of education. Considerable amount of work is done in this direction, but still there are many untouched areas. Moreover, there is no unified approach among these researches. This paper presents a comprehensive survey, a travelogue (2002-2014) towards educational data mining and its scope in future.


title: Some Theorems for Feed Forward Neural Networks

abstract: In this paper we introduce a new method which employs the concept of "Orientation Vectors" to train a feed forward neural network and suitable for problems where large dimensions are involved and the clusters are characteristically sparse. The new method is not NP hard as the problem size increases. We `derive' the method by starting from Kolmogrov's method and then relax some of the stringent conditions. We show for most classification problems three layers are sufficient and the network size depends on the number of clusters. We prove as the number of clusters increase from N to N+dN the number of processing elements in the first layer only increases by d(logN), and are proportional to the number of classes, and the method is not NP hard. Many examples are solved to demonstrate that the method of Orientation Vectors requires much less computational effort than Radial Basis Function methods and other techniques wherein distance computations are required, in fact the present method increases logarithmically with problem size compared to the Radial Basis Function method and the other methods which depend on distance computations e.g statistical methods where probabilistic distances are calculated. A practical method of applying the concept of Occum's razor to choose between two architectures which solve the same classification problem has been described. The ramifications of the above findings on the field of Deep Learning have also been briefly investigated and we have found that it directly leads to the existence of certain types of NN architectures which can be used as a "mapping engine", which has the property of "invertibility", thus improving the prospect of their deployment for solving problems involving Deep Learning and hierarchical classification. The latter possibility has a lot of future scope in the areas of machine learning and cloud computing.


title: Boundary controllability for a one-dimensional heat equation with a singular inverse-square potential

abstract: We analyze controllability properties for the one-dimensional heat equation with singular inverse-square potential $$ u_t-u_{xx}-\frac{\mu}{x^2}u=0,\;\;\; (x,t)\in(0,1)\times(0,T).$$ For any $\mu<1/4$, we prove that the equation is null controllable through a boundary control $f\in H^1(0,T)$ acting at the singularity point $x=0$. This result is obtained employing the moment method by Fattorini and Russell.


title: A pharmacokinetic -- viral kinetic model describes the effect of alisporivir monotherapy or in combination with peg-IFN on 2 hepatitis C virologic response

abstract: Alisporivir is a cyclophilin inhibitor with demonstrated in vitro and in vivo activity against hepatitis C 11 virus (HCV). We estimated antiviral effectiveness of alisporivir alone or in combination with 12 pegylated-Inteferon (peg-IFN) in 88 patients infected with different HCV genotypes treated for four 13 weeks. The pharmacokinetics of both drugs were modeled and used as driving functions for the viral 14 kinetic model. Genotype was found to significantly affect pegylated-Inteferon effectiveness ($\epsilon$= 86.3% 15 and 99.1% in genotype-1/4 and genotype-2/3, respectively, p\textless{}10 -7) and infected cells loss rate ($\delta$= 16 0.22 vs 0.39 day -1 in genotype-1/4 and genotype-2/3, respectively, p\textless{}10 -6). Alisporivir effectiveness 17 was not significantly different across genotype and was high for doses $\ge$600 mg QD. We simulated 18 virologic responses with other alisporivir dosing regimens in HCV genotype-2/3 patients using the 19 model. Our predictions consistently matched the observed responses, demonstrating that this model 20 could be a useful tool for anticipating virologic response and optimize alisporivir-based therapies.


title: A magnetic Impurity in a Weyl semimetal

abstract: We utilize the variational method to study the Kondo screening of a spin-$1/2$ magnetic impurity in a three-dimensional (3D) Weyl semimetal with two Weyl nodes along the $k_z$-axis. The model reduces to a 3D Dirac semimetal when the separation of the two Weyl nodes vanishes. When the chemical potential lies at the nodal point, $\mu=0$, the impurity spin is screened only if the coupling between the impurity and the conduction electron exceeds a critical value. For finite but small $\mu$, the impurity spin is weakly bound due to the low density of state, which is proportional to $\mu^2$, contrary to that in a 2D Dirac metal such as graphene and 2D helical metal where the density of states is proportional to $|\mu|$. The spin-spin correlation function $J_{uv}(\mathbf{r})$ between the spin $v$-component of the magnetic impurity at the origin and the spin $u$-component of a conduction electron at spatial point $\mathbf{r}$, is found to be strongly anisotropic due to the spin-orbit coupling, and it decays in the power-law. The main difference of the Kondo screening in 3D Weyl semimetals and in Dirac semimetals is in the spin $x$- ($y$-) component of the correlation function in the spatial direction of the $z$-axis.


title: Efficient Task Collaboration with Execution Uncertainty

abstract: We study a general task allocation problem, involving multiple agents that collaboratively accomplish tasks and where agents may fail to successfully complete the tasks assigned to them (known as execution uncertainty). The goal is to choose an allocation that maximises social welfare while taking their execution uncertainty into account. We show that this can be achieved by using the post-execution verification (PEV)-based mechanism if and only if agents' valuations satisfy a multilinearity condition. We then consider a more complex setting where an agent's execution uncertainty is not completely predictable by the agent alone but aggregated from all agents' private opinions (known as trust). We show that PEV-based mechanism with trust is still truthfully implementable if and only if the trust aggregation is multilinear.


title: Ground state patterns and phase transition of spin-1 Bose-Einstein condensates via \Gamma-convergence theory

abstract: We develop an analytic theory for the ground state patterns and their phase transitions for spin-1 Bose-Einstein condensates on a bounded domain in the presence of a uniform magnetic field. Within the Thomas-Fermi approximation, these ground state patterns are composed of four basic states: magnetic state, nematic state, two-component state and three-component state, separated by interfaces. A complete phase diagram of the ground state patterns are found analytically with different quadratic Zeeman energy q and total magnetization M for both ferromagnetic and antiferromagnetic systems. Using the \Gamma-convergence technique, it is found that the semi-classical limits of these ground states minimize an energy functional which consists of interior interface energy plus a boundary contact energy. As a consequence, the interface between two different basic states has constant mean curvature, and the contact angle between the interface and the boundary obeys Young's relation.


title: Continuous utility factor in segregation models

abstract: We consider the constrained Schelling model of social segregation in which the utility factor of agents strictly increases and non-local jumps of the agents are allowed. In the present study, the utility factor u is defined in a way such that it can take continuous values and depends on the tolerance threshold as well as the fraction of unlike neighbours. Two models are proposed: in model A the jump probability is determined by the sign of u only which makes it equivalent to the discrete model. In model B the actual values of u are considered. Model A and model B are shown to differ drastically as far as segregation behaviour and phase transitions are concerned. In model A, although segregation can be achieved, the cluster sizes are rather small. Also, a frozen state is obtained in which steady states comprise of many unsatisfied agents. In model B, segregated states with much larger cluster sizes are obtained. The correlation function is calculated to show quantitatively that larger clusters occur in model B. Moreover for model B, no frozen states exist even for very low dilution and small tolerance parameter. This is in contrast to the unconstrained discrete model considered earlier where agents can move even when utility remains same. In addition, we also consider a few other dynamical aspects which have not been studied in segregation models earlier.


title: Optimization of the Superconducting Linear Magnetic Bearing of a Maglev Vehicle

abstract: Considering the need for cost/performance prediction and optimization of superconducting maglev vehicles, we develop and validate here a 3D finite element model to simulate superconducting linear magnetic bearings. Then we reduce the 3D model to a 2D model in order to decrease the computing time. This allows us to perform in a reasonable time a stochastic optimization considering the superconductor properties and the vehicle operation. We look for the permanent magnet guideway geometry that minimizes the cost and maximizes the lateral force during a displacement sequence, with a constraint on the minimum levitation force. The displacement sequence reproduces a regular maglev vehicle operation with both vertical and lateral movements. For the sake of comparison, our reference is the SupraTrans prototype bearing. The results of the optimization suggest that the bearing cost could be substantially reduced, while keeping the same performances as the initial design. Alternatively, the performances could be significantly improved for the same original cost.


title: Nanoscale nonlinear effects in Erbium-implanted Yttrium Orthosilicate

abstract: Doping of substrates at desired locations is a key technology for spin-based quantum memory devices. Focused ion beam implantation is well-suited for this task due to its high spacial resolution. In this work, we investigate ion-beam implanted erbium ensembles in Yttrium Orthosilicate crystals by means of confocal photoluminescence spectroscopy. The sample temperature and the post-implantation annealing step strongly reverberate in the properties of the implanted ions. We find that hot implantation leads to a higher activation rate of the ions. At high enough fluences, the relation between the fluence and final concentration of ions becomes non-linear. Two models are developed explaining the observed behaviour.


title: Accelerated Distance Computation with Encoding Tree for High Dimensional Data

abstract: We propose a novel distance to calculate distance between high dimensional vector pairs, utilizing vector quantization generated encodings. Vector quantization based methods are successful in handling large scale high dimensional data. These methods compress vectors into short encodings, and allow efficient distance computation between an uncompressed vector and compressed dataset without decompressing explicitly. However for large datasets, these distance computing methods perform excessive computations. We avoid excessive computations by storing the encodings on an Encoding Tree(E-Tree), interestingly the memory consumption is also lowered. We also propose Encoding Forest(E-Forest) to further lower the computation cost. E-Tree and E-Forest is compatible with various existing quantization-based methods. We show by experiments our methods speed-up distance computing for high dimensional data drastically, and various existing algorithms can benefit from our methods.


title: Uniqueness of W*-tensor Products

abstract: In contrast to C*-algebras, distinct C*-norms on the algebraic tensor product of two W*-algebras produce isomorphic W*-tensor products


title: The Bulk and The Tail of Minimal Absent Words in Genome Sequences

abstract: Minimal absent words (MAW) of a genomic sequence are subsequences that are absent themselves but the subwords of which are all present in the sequence. The characteristic distribution of genomic MAWs as a function of their length has been observed to be qualitatively similar for all living organisms, the bulk being rather short, and only relatively few being long. It has been an open issue whether the reason behind this phenomenon is statistical or reflects a biological mechanism, and what biological information is contained in absent words. In this work we demonstrate that the bulk can be described by a probabilistic model of sampling words from random sequences, while the tail of long MAWs is of biological origin. We introduce the novel concept of a core of a minimal absent word, which are sequences present in the genome and closest to a given MAW. We show that in bacteria and yeast the cores of the longest MAWs, which exist in two or more copies, are located in highly conserved regions the most prominent example being ribosomal RNAs (rRNAs). We also show that while the distribution of the cores of long MAWs is roughly uniform over these genomes on a coarse-grained level, on a more detailed level it is strongly enhanced in 3' untranslated regions (UTRs) and, to a lesser extent, also in 5' UTRs. This indicates that MAWs and associated MAW cores correspond to fine-tuned evolutionary relationships, and suggest that they can be more widely used as markers for genomic complexity.


title: Grand canonical ensemble, multi-particle wave functions, scattering data, and lattice field theories

abstract: We show that information about scattering data of a quantum field theory can be obtained from studying the system at finite density and low temperatures. In particular we consider models formulated on the lattice which can be exactly dualized to theories of conserved charge fluxes on lattice links. Apart from eliminating the complex action problem at nonzero chemical potential mu, these dualizations allow for a particle world line interpretation of the dual fluxes from which one can extract data about the 2-particle wave function. As an example we perform dual Monte Carlo simulations of the 2-dimensional O(3) model at nonzero mu and finite volume, whose non-perturbative spectrum consists of a massive triplet of particles. At nonzero mu particles are induced in the system, which at sufficiently low temperature give rise to sectors of fixed particle number. We show that the scattering phase shifts can be obtained either from the critical chemical potential values separating the sectors or directly from the wave function in the 2-particle sector. We find that both methods give excellent agreement with the exact result. We discuss the applicability and generality of the new approaches.


title: Measurement of jet charge in dijet events from sqrt(s)=8 TeV pp collisions with the ATLAS detector

abstract: The momentum-weighted sum of the charges of tracks associated to a jet is sensitive to the charge of the initiating quark or gluon. This paper presents a measurement of the distribution of momentum-weighted sums, called jet charge, in dijet events using 20.3 fb$^{-1}$ of data recorded with the ATLAS detector at $\sqrt{s}=8$ TeV in $pp$ collisions at the LHC. The jet charge distribution is unfolded to remove distortions from detector effects and the resulting particle-level distribution is compared with several models. The $p_T$-dependence of the jet charge distribution average and standard deviation are compared to predictions obtained with several LO and NLO parton distribution functions. The data are also compared to different Monte Carlo simulations of QCD dijet production using various settings of the free parameters within these models. The chosen value of the strong coupling constant used to calculate gluon radiation is found to have a significant impact on the predicted jet charge. There is evidence for a $p_{T}$-dependence of the jet charge distribution for a given jet flavor. In agreement with perturbative QCD predictions, the data show that the average jet charge of quark-initiated jets decreases in magnitude as the energy of the jet increases.


title: The H.E.S.S. II GRB Program

abstract: Gamma-ray bursts (GRBs) are some of the most energetic and exotic events in the Universe, however their behaviour at the highest energies (>10 GeV) is largely unknown. Although the Fermi-LAT space telescope has detected several GRBs in this energy range, it is limited by the relatively small collection area of the instrument. The H.E.S.S. experiment has now entered its second phase by adding a fifth telescope of 600 m$^{2}$ mirror area to the centre of the array. This new telescope increases the energy range of the array, allowing it to probe the sub-100 GeV range while maintaining the large collection area of ground based gamma-ray observatories, essential to probing short-term variability at these energies. We will present a description of the GRB observation scheme used by the H.E.S.S. experiment, summarising the behaviour and performance of the rapid GRB repointing system, the conditions under which potential GRB repointings are made and the data analysis scheme used for these observations.


title: A balanced rail-to-rail all digital comparator using only standard cells

abstract: An all-digital comparator with full input range is presented. It outperforms the nowaday all-digital comparators with its large rail-to-rail input range. This is achieved by the proposed Yin-yang balance mechanism between the two logic gates: NAND3 and OAI (Or-And-Invert). The important design considerations to achieve this balance are presented, such as the driving strength manipulation and the use of pre-distortion technique. Constructed only by commercially available digital standard cells, the layout of the proposed comparator is generated automatically by standard digital Place & Route routine within several minutes. The Verilog code for the proposed circuit is given, and the circuit is successfully implemented in 130nm CMOS technology with the power consumption of 0.176mW at the clock of 330MHz.


title: Was the soft X-ray flare in NGC 3599 due to an AGN disc instability or a delayed tidal disruption event?

abstract: We present unpublished data from a tidal disruption candidate in NGC 3599 which show that the galaxy was already X-ray bright 18 months before the measurement which led to its classification. This removes the possibility that the flare was caused by a classical, fast-rising, short-peaked, tidal disruption event. Recent relativistic simulations indicate that the majority of disruptions will actually take months or years to rise to a peak, which will then be maintained for longer than previously thought. NGC 3599 could be one of the first identified examples of such an event. The optical spectra of NGC 3599 indicate that it is a low-luminosity Seyfert/LINER with L_bol~10^40 ergs/s The flare may alternatively be explained by a thermal instability in the accretion disc, which propagates through the inner region at the sound speed, causing an increase of the disc scale height and local accretion rate. This can explain the <9 years rise time of the flare. If this mechanism is correct then the flare may repeat on a timescale of several decades as the inner disc is emptied and refilled.


title: HCLAE: High Capacity Locally Aggregating Encodings for Approximate Nearest Neighbor Search

abstract: Vector quantization-based approaches are successful to solve Approximate Nearest Neighbor (ANN) problems which are critical to many applications. The idea is to generate effective encodings to allow fast distance approximation. We propose quantization-based methods should partition the data space finely and exhibit locality of the dataset to allow efficient non-exhaustive search. In this paper, we introduce the concept of High Capacity Locality Aggregating Encodings (HCLAE) to this end, and propose Dictionary Annealing (DA) to learn HCLAE by a simulated annealing procedure. The quantization error is lower than other state-of-the-art. The algorithms of DA can be easily extended to an online learning scheme, allowing effective handle of large scale data. Further, we propose Aggregating-Tree (A-Tree), a non-exhaustive search method using HCLAE to perform efficient ANN-Search. A-Tree achieves magnitudes of speed-up on ANN-Search tasks, compared to the state-of-the-art.


title: Improved Residual Vector Quantization for High-dimensional Approximate Nearest Neighbor Search

abstract: Quantization methods have been introduced to perform large scale approximate nearest search tasks. Residual Vector Quantization (RVQ) is one of the effective quantization methods. RVQ uses a multi-stage codebook learning scheme to lower the quantization error stage by stage. However, there are two major limitations for RVQ when applied to on high-dimensional approximate nearest neighbor search: 1. The performance gain diminishes quickly with added stages. 2. Encoding a vector with RVQ is actually NP-hard. In this paper, we propose an improved residual vector quantization (IRVQ) method, our IRVQ learns codebook with a hybrid method of subspace clustering and warm-started k-means on each stage to prevent performance gain from dropping, and uses a multi-path encoding scheme to encode a vector with lower distortion. Experimental results on the benchmark datasets show that our method gives substantially improves RVQ and delivers better performance compared to the state-of-the-art.


title: A Class of Prediction-Correction Methods for Time-Varying Convex Optimization

abstract: This paper considers unconstrained convex optimization problems with time-varying objective functions. We propose algorithms with a discrete time-sampling scheme to find and track the solution trajectory based on prediction and correction steps, while sampling the problem data at a constant rate of $1/h$, where $h$ is the length of the sampling interval. The prediction step is derived by analyzing the iso-residual dynamics of the optimality conditions. The correction step adjusts for the distance between the current prediction and the optimizer at each time step, and consists either of one or multiple gradient steps or Newton steps, which respectively correspond to the gradient trajectory tracking (GTT) or Newton trajectory tracking (NTT) algorithms. Under suitable conditions, we establish that the asymptotic error incurred by both proposed methods behaves as $O(h^2)$, and in some cases as $O(h^4)$, which outperforms the state-of-the-art error bound of $O(h)$ for correction-only methods in the gradient-correction step. Moreover, when the characteristics of the objective function variation are not available, we propose approximate gradient and Newton tracking algorithms (AGT and ANT, respectively) that still attain these asymptotical error bounds. Numerical simulations demonstrate the practical utility of the proposed methods and that they improve upon existing techniques by several orders of magnitude.


title: A Reliable and Cost-Efficient Auto-Scaling System for Web Applications Using Heterogeneous Spot Instances

abstract: Cloud providers sell their idle capacity on markets through an auction-like mechanism to increase their return on investment. The instances sold in this way are called spot instances. In spite that spot instances are usually 90% cheaper than on-demand instances, they can be terminated by provider when their bidding prices are lower than market prices. Thus, they are largely used to provision fault-tolerant applications only. In this paper, we explore how to utilize spot instances to provision web applications, which are usually considered availability-critical. The idea is to take advantage of differences in price among various types of spot instances to reach both high availability and significant cost saving. We first propose a fault-tolerant model for web applications provisioned by spot instances. Based on that, we devise novel auto-scaling polices for hourly billed cloud markets. We implemented the proposed model and policies both on a simulation testbed for repeatable validation and Amazon EC2. The experiments on the simulation testbed and the real platform against the benchmarks show that the proposed approach can greatly reduce resource cost and still achieve satisfactory Quality of Service (QoS) in terms of response time and availability.


title: A Paley-like graph in characteristic two

abstract: The Paley graph is a well-known self-complementary pseudo-random graph, defined over a finite field of odd order. We describe an attempt at an analogous construction using fields of even order. Some properties of the graph are noted, such as the existence of a Hamiltonian decomposition.


title: Singularity analysis for heavy-tailed random variables

abstract: We propose a novel complex-analytic method for sums of i.i.d. random variables that are heavy-tailed and integer-valued. The method combines singularity analysis, Lindel\"of integrals, and bivariate saddle points. As an application, we prove three theorems on precise large and moderate deviations which provide a local variant of a result by S. V. Nagaev (1973). The theorems generalize five theorems by A. V. Nagaev (1968) on stretched exponential laws $p(k) = c\exp( -k^\alpha)$ and apply to logarithmic hazard functions $c\exp( - (\log k)^\beta)$, $\beta>2$; they cover the big jump domain as well as the small steps domain. The analytic proof is complemented by clear probabilistic heuristics. Critical sequences are determined with a non-convex variational problem.


title: Notions of maximality for integral lattice-free polyhedra: the case of dimension three

abstract: Lattice-free sets (convex subsets of $\mathbb{R}^d$ without interior integer points) and their applications for cutting-plane methods in mixed-integer optimization have been studied in recent literature. Notably, the family of all integral lattice-free polyhedra which are not properly contained in another integral lattice-free polyhedron has been of particular interest. We call these polyhedra $\mathbb{Z}^d$-maximal. It is known that, for fixed $d$, the family $\mathbb{Z}^d$-maximal integral lattice-free polyhedra is finite up to unimodular equivalence. In view of possible applications in cutting-plane theory, one would like to have a classification of this family. However, this turns out to be a challenging task already for small dimensions. In contrast, the subfamily of all integral lattice-free polyhedra which are not properly contained in any other lattice-free set, which we call $\mathbb{R}^d$-maximal lattice-free polyhedra, allow a rather simple geometric characterization. Hence, the question was raised for which dimensions the notions of $\mathbb{Z}^d$-maximality and $\mathbb{R}^d$-maximality are equivalent. This was known to be the case for dimensions one and two. On the other hand, Nill and Ziegler (2011) showed that for dimension $d \ge 4$, there exist polyhedra which are $\mathbb{Z}^d$-maximal but not $\mathbb{R}^d$-maximal. In this article, we consider the remaining case $d = 3$ and prove that for integral polyhedra the notions of $\mathbb{R}^3$-maximality and $\mathbb{Z}^3$-maximality are equivalent. As a consequence, the classification of all $\mathbb{R}^3$-maximal integral polyhedra by Averkov, Wagner and Weismantel (2011) contains all $\mathbb{Z}^3$-maximal integral polyhedra.


title: Spin-ice behavior of three-dimensional inverse opal-like magnetic structures: micromagnetic simulations

abstract: We perform micromagnetic simulations of the magnetization distribution in inverse opal-like structures (IOLS) made from ferromagnetic materials (nickel and cobalt). It is shown that the unit cell of these complex structures, whose characteristic length is approximately 700 nm, can be divided into a set of structural elements some of which behave like Ising-like objects. A spin-ice behavior of IOLS is observed in a broad range of external magnetic fields. Numerical results describe successfully the experimental hysteresis curves of the magnetization in Ni- and Co-based IOLS. We conclude that ferromagnetic IOLS can be considered as the first realization of three-dimensional artificial spin ice.


title: Constructive Solutions to the Riemann--Hilbert Problem and Middle Convolution

abstract: In this paper we present a general scheme to generate constructive solutions to the Riemann--Hilbert problem via middle convolution and illustrate this approach for a Fuchsian system with four singular points.