title: Prevalence problem in the set of quadratic stochastic operators acting on L1

abstract: This paper is devoted to the study of the problem of prevalence in the class of quadratic stochastic operators acting on the L1 space for the uniform topology. We obtain that the set of norm quasi-mixing quadratic stochastic operators is a dense and open set in the topology induced by a very natural metric. This shows the typical long-term behaviour of iterates of quadratic stochastic operators.


title: Injectable Spontaneous Generation of Tremendous Self-Fueled Liquid Metal Droplet Motors in a Moment

abstract: Micro motors that could run in liquid environment is very important for a variety of practices such as serving as pipeline robot, soft machine, drug delivery, or microfluidics system etc. However, fabrication of such tiny motors is generally rather time and cost consumptive and has been a tough issue due to involve too many complicated procedures and tools. Here, we show a straightforward injectable way for spontaneously generating autonomously running soft motors in large quantity. A basic fabrication strategy thus enabled is established and illustrated. It was found that, injecting the GaIn alloy pre-fueled with aluminum into electrolyte would automatically split in seconds into tremendous droplet motors swiftly running here and there. The driving force originated from the galvanic cell reaction among alloy, aluminum and surrounding electrolyte which offers interior electricity and hydrogen gas as motion power. This finding opens the possibility to develop injectable tiny-robots, droplet machines or microfluidic elements. It also raised important scientific issues regarding characterizing the complicated fluid mechanics stimulated by the quick running of the soft metal droplet and the gases it generated during the traveling.


title: Fast Estimation of the Median Covariation Matrix with Application to Online Robust Principal Components Analysis

abstract: The geometric median covariation matrix is a robust multivariate indicator of dispersion which can be extended without any difficulty to functional data. We define estimators, based on recursive algorithms, that can be simply updated at each new observation and are able to deal rapidly with large samples of high dimensional data without being obliged to store all the data in memory. Asymptotic convergence properties of the recursive algorithms are studied under weak conditions. The computation of the principal components can also be performed online and this approach can be useful for online outlier detection. A simulation study clearly shows that this robust indicator is a competitive alternative to minimum covariance determinant when the dimension of the data is small and robust principal components analysis based on projection pursuit and spherical projections for high dimension data. An illustration on a large sample and high dimensional dataset consisting of individual TV audiences measured at a minute scale over a period of 24 hours confirms the interest of considering the robust principal components analysis based on the median covariation matrix. All studied algorithms are available in the R package Gmedian on CRAN.


title: Scattering for the focusing $L^{2}$ -supercritical and $\dot H^2$-subcritical biharmonic NLS Equations

abstract: We consider the focusing $\dot H^{s_c}$-critical biharmonic Schr\"odinger equation, and prove a global wellposedness and scattering result for the radial data $u_0\in H^2(\mathbb R^N)$ satisfying $ M(u_0)^{\frac{2-s_c}{s_c}}E(u_0)<M(Q)^{\frac{2-s_c}{s_c}}E(Q) $ and $ \|u_{0}\|^{\frac{2-s_c}{s_c}}_{2}\|\Delta u_{0}\|_{2}<\|Q\|^{\frac{2-s_c}{s_c}}_{2}\|\Delta Q\|_{2}, $ where $s_c\in(0,2)$ and $Q$ is the ground state of $\Delta^2Q+(2-s_c)Q-|Q|^{p-1}Q=0$.


title: Mission Analysis For the Ion Beam Deflection of Fictitious Asteroid 2015PDC

abstract: A realistic mission scenario for the deflection of fictitious asteroid 2015PDC is investigated that makes use of the ion beam shepherd concept as primary deflection technique. The article deals with the design of a low thrust rendezvous trajectory to the asteroid, the estimation of the propagated covariance ellipsoid and the outcome of a slow-push deflection starting from three worst case scenarios (impacts in New Delhi, Dhaka and Teheran). Displacing the impact point towards very low populated areas, as opposed to full deflection, is found to be the simplest and most effective mitigation approach. Mission design, technical and political aspects are discussed.


title: Symmetry consideration in the problem of wave modes of thin viscous liquid layer flow

abstract: The equations in conservative form for nonlinear waves modeling on a liquid film flowing down a vertical plane have been investigated. It has been found that in the computational domain extended along the transverse axis the equations with boundary conditions are invariant under parity transformation. It is numerically shown that for moderate Reynolds numbers the steady-state travelling solutions of the equations have the detected symmetry. It is demonstrated that using this symmetry for the numerical solution of the problem by Galerkin methods significantly increases the efficiency of calculations.


title: High Density Noise Removal by Cascading Algorithms

abstract: An advanced non-linear cascading filter algorithm for the removal of high density salt and pepper noise from the digital images is proposed. The proposed method consists of two stages. The first stage Decision base Median Filter (DMF) acts as the preliminary noise removal algorithm. The second stage is either Modified Decision Base Partial Trimmed Global Mean Filter (MDBPTGMF) or Modified Decision Based Unsymmetric Trimmed Median Filter (MDBUTMF) which is used to remove the remaining noise and enhance the image quality. The DMF algorithm performs well at low noise density but it fails to remove the noise at medium and high level. The MDBPTGMF and MDUTMF have excellent performance at low, medium and high noise density but these reduce the image quality and blur the image at high noise level. So the basic idea behind this paper is to combine the advantages of the filters used in both the stages to remove the Salt and Pepper noise and enhance the image quality at all the noise density level. The proposed method is tested against different gray scale images and it gives better Mean Absolute Error (MAE), Peak Signal to Noise Ratio (PSNR) and Image Enhancement Factor (IEF) than the Adaptive Median Filter (AMF), Decision Base Unsymmetric Trimmed Median Filter (DBUTMF), Modified Decision Base Unsymmetric Trimmed Median Filter (MDBUTMF) and Decision Base Partial Trimmed Global Mean Filter (DBPTGMF).


title: Emergence of coexisting coherence and incoherence by an external forcing

abstract: A common external forcing can cause a saddle-node bifurcation in an ensemble of identical Duffing oscillators by breaking the symmetry of the individual bistable (double-well) unit. The strength of the forcing determines the separation between the saddle and node, which in turn dictates different dynamical transitions depending on the distribution of the initial states of the oscillators. In particular, chimera-like states appear in the vicinity of the saddle-node bifurcation for which theoretical explanation is provided from the stability of slow-scale dynamics of the original system of equations. Further, as a consequence, it is shown that even a linear nearest neighbor coupling can lead to the manifestation of the chimera states in an ensemble of identical Duffing oscillators in the presence of the common external forcing.


title: Optimal Phonon-to-Spin Mapping in a system of a trapped ion

abstract: We propose a protocol for measurement of the phonon number distribution of a harmonic oscillator based on selective mapping to a discrete spin-1/2 degree of freedom. We consider a system of a harmonically trapped ion, where a transition between two long lived states can be driven with resolved motional sidebands. The required unitary transforms are generated by amplitude-modulated polychromatic radiation fields, where the time-domain ramps are obtained from numerical optimization by application of the Chopped RAndom Basis (CRAB) algorithm. We provide a detailed analysis of the scaling behavior of the attainable fidelities and required times for the mapping transform with respect to the size of the Hilbert space. As one application we show how the mapping can be employed as a building block for experiments which require measurement of the work distribution of a quantum process.


title: QAM Adaptive Measurements Feedback Quantum Receiver Performance

abstract: We theoretically study the quantum receivers with adaptive measurements feedback for discriminating quadrature amplitude modulation (QAM) coherent states in terms of average symbol error rate. For rectangular 16-QAM signal set, with different stages of adaptive measurements, the effects of realistic imperfection parameters including the sub-unity quantum efficiency and the dark counts of on-off detectors, as well as the transmittance of beam splitters and the mode mismatch factor between the signal and local oscillating fields on the symbol error rate are separately investigated through Monte Carlo simulations. Using photon-number-resolving detectors (PNRD) instead of on-off detectors, all the effects on the symbol error rate due to the above four imperfections can be suppressed in a certain degree. The finite resolution and PNR capability of PNRDs are also considered. We find that for currently available technology, the receiver shows a reasonable gain from the standard quantum limit (SQL) with moderate stages.


title: Collapse of spin-orbit coupled Bose-Einstein condensates

abstract: A finite-size quasi two-dimensional Bose-Einstein condensate collapses if the attraction between atoms is sufficiently strong. Here we present a theory of collapse for condensates with the interatomic attraction and spin-orbit coupling. We consider two realizations of spin-orbit coupling: the axial Rashba coupling and balanced, effectively one-dimensional, Rashba-Dresselhaus one. In both cases spin-dependent "anomalous" velocity, proportional to the spin-orbit coupling strength, plays a crucial role. For the Rashba coupling, this velocity forms a centrifugal component in the density flux opposite to that arising due to the attraction between particles and prevents the collapse at a sufficiently strong coupling. For the balanced Rashba-Dresselhaus coupling, the spin-dependent velocity can spatially split the initial state in one dimension and form spin-projected wavepackets, reducing the total condensate density. Depending on the spin-orbit coupling strength, interatomic attraction, and the initial state, this splitting either prevents the collapse or modifies the collapse process. These results show that the collapse can be controlled by a spin-orbit coupling, thus, extending the domain of existence of condensates of attracting atoms.


title: Explicit Model Checking of Very Large MDP using Partitioning and Secondary Storage

abstract: The applicability of model checking is hindered by the state space explosion problem in combination with limited amounts of main memory. To extend its reach, the large available capacities of secondary storage such as hard disks can be exploited. Due to the specific performance characteristics of secondary storage technologies, specialised algorithms are required. In this paper, we present a technique to use secondary storage for probabilistic model checking of Markov decision processes. It combines state space exploration based on partitioning with a block-iterative variant of value iteration over the same partitions for the analysis of probabilistic reachability and expected-reward properties. A sparse matrix-like representation is used to store partitions on secondary storage in a compact format. All file accesses are sequential, and compression can be used without affecting runtime. The technique has been implemented within the Modest Toolset. We evaluate its performance on several benchmark models of up to 3.5 billion states. In the analysis of time-bounded properties on real-time models, our method neutralises the state space explosion induced by the time bound in its entirety.


title: Coherence measures and optimal conversion for coherent states

abstract: We discuss a general strategy to construct coherence measures. One can build an important class of coherence measures which cover the relative entropy measure for pure states, the $l_1$-norm measure for pure states and the $\alpha$-entropy measure. The optimal conversion of coherent states under incoherent operations is presented which sheds some light on the coherence of a single copy of a pure state.


title: Appearance-Based Gaze Estimation in the Wild

abstract: Appearance-based gaze estimation is believed to work well in real-world settings, but existing datasets have been collected under controlled laboratory conditions and methods have been not evaluated across multiple datasets. In this work we study appearance-based gaze estimation in the wild. We present the MPIIGaze dataset that contains 213,659 images we collected from 15 participants during natural everyday laptop use over more than three months. Our dataset is significantly more variable than existing ones with respect to appearance and illumination. We also present a method for in-the-wild appearance-based gaze estimation using multimodal convolutional neural networks that significantly outperforms state-of-the art methods in the most challenging cross-dataset evaluation. We present an extensive evaluation of several state-of-the-art image-based gaze estimation algorithms on three current datasets, including our own. This evaluation provides clear insights and allows us to identify key research challenges of gaze estimation in the wild.


title: Orbital Feshbach Resonance in Alkali-Earth Atoms

abstract: For a mixture of alkali-earth atomic gas in the long-lived excited state ${}^3P_0$ and ground state ${}^1S_0$, in addition to nuclear spin, another "orbital" index is introduced to distinguish these two internal states. In this letter we propose a mechanism to induce Feshbach resonance between two atoms with different orbital and nuclear spin quantum numbers. Two essential ingredients are inter-orbital spin-exchanging scattering and orbital dependence of the Land\'e g-factors. Here the orbital degrees of freedom plays similar role as electron spin degree of freedom in magnetic Feshbach resonance in alkali-metal atoms. This resonance is particularly accessible for ${}^{173}$Yb system. The BCS-BEC crossover in this system requires two fermion pairing order parameters, and displays significant difference comparing to that in alkali-metal system.


title: Analytic Solutions to 3-D Finite Deformation Problems Governed by St Venant-Kirchhoff Material

abstract: This paper presents a detailed study on analytical solutions to a general nonlinear boundary-value problem in finite deformation theory. Based on canonical duality theory and the associated pure complementary energy principle in nonlinear elasticity proposed by Gao in 1999, we show that the general nonlinear partial differential equation for deformation is actually equivalent to an algebraic (tensor) equation in stress space. For St Venant-Kirchhoff materials, this coupled cubic algebraic equation can be solved principally to obtain all possible solutions. Our results show that for any given external source field such that the statically admissible first Piola-Kirchhoff stress field has no-zero eigenvalues, the problem has a unique global minimal solution, which is corresponding to a positive-definite second Piola-Kirchhoff stress S, and at most eight local solutions corresponding to negative-definite S. Additionally, the problem could have 15 unstable solutions corresponding to indefinite S. This paper demonstrates that the canonical duality theory and the pure complementary energy principle play fundamental roles in nonconvex analysis and finite deformation theory.


title: A Delay Aware Routing Protocol for Wireless Sensor Networks

abstract: Wireless Sensor Networks (WSNs) consist of sensor nodes which can be deployed for various operations such as agriculture and environmental sensing, wild life monitoring, health care, military surveillance, industrial control, home automation, security etc. Quality of Service (QoS) is an important issue in wireless sensor networks (WSNs) and providing QoS support in WSNs is an emerging area of research. Due to resource constraints nature of sensor networks like processing power, memory, bandwidth, energy etc. providing QoS support in WSNs is a challenging task. Delay is an important QoS parameter for forwarding data in a time constraint WSNs environment. In this paper we propose a delay aware routing protocol for transmission of time critical event information to the Sink of WSNs. The performance of the proposed protocol is evaluated by NS2 simulations under different scenarios.


title: Interplay between superconductivity and chiral symmetry breaking in a (2+1)-dimensional model with compactified spatial coordinate

abstract: In this paper a (2+1)-dimensional model with four-fermion interactions is investigated in the case when one spatial coordinate is compactified and the space topology takes the form of an infinite cylinder, $R^1\otimes S^1$. It is supposed that the system is embedded into real three-dimensional space and that a magnetic flux $\Phi$ crosses the transverse section of the cylinder. The model includes four-fermion interactions both in the fermion-antifermion (or chiral) and fermion-fermion (or superconducting) channels. We then study phase transitions in dependence on the chemical potential $\mu$ and the flux $\Phi$ in the leading order of the large-$N$ expansion technique, where $N$ is the number of fermion fields. It is demonstrated that for arbitrary relations between coupling constants in the chiral and superconducting channels, superconductivity appears in the system at rather high values of $\mu$ (the length $L$ of the circumference $S^1$ is fixed). Moreover, it is shown that at sufficiently small values of $\mu$ the growth of the magnetic flux $\Phi$ leads to a periodical reentrance of the chiral symmetry breaking or superconducting phase, depending on the values of $\mu$ and coupling constants.


title: Entropic Topological Invariants in Three Dimensions

abstract: We evaluate the entanglement entropy of exactly solvable Hamiltonians corresponding to general families of three-dimensional topological models. We show that the modification to the entropic area law due to three-dimensional topological properties is richer than the two-dimensional case. In addition to the reduction of the entropy caused by non-zero vacuum expectation value of contractible loop operators a new topological invariant appears that increases the entropy if the model consists of non-trivially braiding anyons. As a result the three-dimensional topological entanglement entropy provides only partial information about the two entropic topological invariants.


title: Far field imaging of a dielectric inclusion

abstract: A non-iterative topological sensitivity framework for guaranteed far field detection of a dielectric inclusion is presented. The cases of single and multiple measurements of the electric far field scattering amplitude at a fixed frequency are taken into account. The performance of the algorithm is analyzed theoretically in terms of resolution, stability, and signal-to-noise ratio.


title: Quick sensitivity analysis for incremental data modification and its application to leave-one-out CV in linear classification problems

abstract: We introduce a novel sensitivity analysis framework for large scale classification problems that can be used when a small number of instances are incrementally added or removed. For quickly updating the classifier in such a situation, incremental learning algorithms have been intensively studied in the literature. Although they are much more efficient than solving the optimization problem from scratch, their computational complexity yet depends on the entire training set size. It means that, if the original training set is large, completely solving an incremental learning problem might be still rather expensive. To circumvent this computational issue, we propose a novel framework that allows us to make an inference about the updated classifier without actually re-optimizing it. Specifically, the proposed framework can quickly provide a lower and an upper bounds of a quantity on the unknown updated classifier. The main advantage of the proposed framework is that the computational cost of computing these bounds depends only on the number of updated instances. This property is quite advantageous in a typical sensitivity analysis task where only a small number of instances are updated. In this paper we demonstrate that the proposed framework is applicable to various practical sensitivity analysis tasks, and the bounds provided by the framework are often sufficiently tight for making desired inferences.


title: Vertical transport in InAs/GaSb superlattice: role of ionized impurity and interface roughness scatterings

abstract: We report on the vertical electron mobility versus temperature by applying the interface roughness scattering and ionized impurity scattering in InAs/GaSb superlattices. Using the Finite difference K.P method, we calculated the band structure of InAs/Gasb superlattices and then studied the transport properties of these systems. Several structural parameters such as layer thicknesses, interface roughness height, correlation length and ion density have been investigated and characterized that how the vertical mobility change with varying these parameters. Theoretical modeling results show that these two scattering mechanism are important in lower temperatures and thin layer systems.


title: Regularized system identification using orthonormal basis functions

abstract: Most of existing results on regularized system identification focus on regularized impulse response estimation. Since the impulse response model is a special case of orthonormal basis functions, it is interesting to consider if it is possible to tackle the regularized system identification using more compact orthonormal basis functions. In this paper, we explore two possibilities. First, we construct reproducing kernel Hilbert space of impulse responses by orthonormal basis functions and then use the induced reproducing kernel for the regularized impulse response estimation. Second, we extend the regularization method from impulse response estimation to the more general orthonormal basis functions estimation. For both cases, the poles of the basis functions are treated as hyperparameters and estimated by empirical Bayes method. Then we further show that the former is a special case of the latter, and more specifically, the former is equivalent to ridge regression of the coefficients of the orthonormal basis functions.


title: A new perspective on Einstein's philosophy of cosmology

abstract: The recent discovery that Einstein once attempted - and quickly abandoned - a steady-state model of the expanding universe sheds new light on his philosophical journey from static to dynamic cosmologies.


title: Magnetic response of a disordered binary ferromagnetic alloy to an oscillating magnetic field

abstract: By means of Monte Carlo simulation with local spin update Metropolis algorithm, we have elucidated non-equilibrium phase transition properties and stationary-state treatment of a disordered binary ferromagnetic alloy of the type $A_{p}B_{1-p}$ on a square lattice. After a detailed analysis, we have found that the system shows many interesting and unusual thermal and magnetic behaviors, for instance, the locations of dynamic phase transition points change significantly depending upon amplitude and period of the external magnetic field as well as upon the active concentration of $A-$ type components. Much effort has also been dedicated to clarify the hysteresis tools, such as coercivity, dynamic loop area as well as dynamic correlations between time dependent magnetizations and external time dependent applied field as a functions of period and amplitude of field as well as active concentration of of $A-$ type components, and outstanding physical findings have been reported in order to better understand the dynamic process underlying present system.


title: Discovery of the $D$-basis in binary tables based on hypergraph dualization

abstract: Discovery of (strong) association rules, or implications, is an important task in data management, and it finds application in artificial intelligence, data mining and the semantic web. We introduce a novel approach for the discovery of a specific set of implications, called the $D$-basis, that provides a representation for a reduced binary table, based on the structure of its Galois lattice. At the core of the method are the $D$-relation defined in the lattice theory framework, and the hypergraph dualization algorithm that allows us to effectively produce the set of transversals for a given Sperner hypergraph. The latter algorithm, first developed by specialists from Rutgers Center for Operations Research, has already found numerous applications in solving optimization problems in data base theory, artificial intelligence and game theory. One application of the method is for analysis of gene expression data related to a particular phenotypic variable, and some initial testing is done for the data provided by the University of Hawaii Cancer Center.


title: The Hardness of Subgraph Isomorphism

abstract: Subgraph Isomorphism is a very basic graph problem, where given two graphs $G$ and $H$ one is to check whether $G$ is a subgraph of $H$. Despite its simple definition, the Subgraph Isomorphism problem turns out to be very broad, as it generalizes problems such as Clique, $r$-Coloring, Hamiltonicity, Set Packing and Bandwidth. However, for all of the mentioned problems $2^{\mathcal{O}(n)}$ time algorithms exist, so a natural and frequently asked question in the past was whether there exists a $2^{\mathcal{O}(n)}$ time algorithm for Subgraph Isomorphism. In the monograph of Fomin and Kratsch [Springer'10] this question is highlighted as an open problem, among few others. Our main result is a reduction from 3-SAT, producing a subexponential number of sublinear instances of the Subgraph Isomorphism problem. In particular, our reduction implies a $2^{\Omega(n \sqrt{\log n})}$ lower bound for Subgraph Isomorphism under the Exponential Time Hypothesis. This shows that there exist classes of graphs that are strictly harder to embed than cliques or Hamiltonian cycles. The core of our reduction consists of two steps. First, we preprocess and pack variables and clauses of a 3-SAT formula into groups of logarithmic size. However, the grouping is not arbitrary, since as a result we obtain only a limited interaction between the groups. In the second step, we overcome the technical hardness of encoding evaluations as permutations by a simple, yet fruitful scheme of guessing the sizes of preimages of an arbitrary mapping, reducing the case of arbitrary mapping to bijections. In fact, when applying this step to a recent independent result of Fomin et al.[arXiv:1502.05447 (2015)], who showed hardness of Graph Homomorphism, we can transfer their hardness result to Subgraph Isomorphism, implying a nearly tight lower bound of $2^{\Omega(n \log n / \log \log n)}$.


title: Role of Non-quantized Fluxes in Coulombic and Casimir Scaling Regimes of the Thick Center Vortex Potentials

abstract: We discuss a non-quantized Gaussian flux limited to vary in a specific region of space which describes the Coulombic potential as well as confinement at intermediate distances for color potentials in various representations. Our calculations show if uncorrelated vortices are not quantized with center elements, the N-ality dependence of asymptotic string tensions is lost while the ratio of the Coulombic strengths and intermediate string tensions are in agreement with Casimir scaling. Using both uncorrelated quantized and non-quantized fluxes in the potential between static color sources results in the correct N-ality dependence at large distances and very good agreement with Casimir scaling at short and intermediate distances.


title: Data Science and Ebola

abstract: Data Science---Today, everybody and everything produces data. People produce large amounts of data in social networks and in commercial transactions. Medical, corporate, and government databases continue to grow. Sensors continue to get cheaper and are increasingly connected, creating an Internet of Things, and generating even more data. In every discipline, large, diverse, and rich data sets are emerging, from astrophysics, to the life sciences, to the behavioral sciences, to finance and commerce, to the humanities and to the arts. In every discipline people want to organize, analyze, optimize and understand their data to answer questions and to deepen insights. The science that is transforming this ocean of data into a sea of knowledge is called data science. This lecture will discuss how data science has changed the way in which one of the most visible challenges to public health is handled, the 2014 Ebola outbreak in West Africa.


title: Magnetic Nanorods Confined in a Lamellar Lyotropic Phase

abstract: The dilute lamellar phase of the nonionic surfactant C$_{12}$EO$_5$ was doped with goethite (iron oxide) nanorods up to a fraction of 5 vol. %. The interaction between the inclusions and the host phase was studied by polarized optical microscopy (with or without an applied magnetic field) and by small-angle X-ray scattering. We find that, when the orientation of the nanorods is modified using the magnetic field, the texture of the lamellar phase changes accordingly; one can thus induce a homeotropic-planar reorientation transition. On the other hand, the lamellar phase induces an attractive interaction between the nanorods. In more concentrated lamellar phases (under stronger confinement) the particles form aggregates. This behavior is not encountered for a similar system doped with spherical particles, emphasizing the role of particle shape in the interaction between doping particles and the host phase.


title: Jacobi stability analysis of the Lorenz system

abstract: We perform the study of the stability of the Lorenz system by using the Jacobi stability analysis, or the Kosambi-Cartan-Chern (KCC) theory. The Lorenz model plays an important role for understanding hydrodynamic instabilities and the nature of the turbulence, also representing a non-trivial testing object for studying non-linear effects. The KCC theory represents a powerful mathematical method for the analysis of dynamical systems. In this approach we describe the evolution of the Lorenz system in geometric terms, by considering it as a geodesic in a Finsler space. By associating a non-linear connection and a Berwald type connection, five geometrical invariants are obtained, with the second invariant giving the Jacobi stability of the system. The Jacobi (in)stability is a natural generalization of the (in)stability of the geodesic flow on a differentiable manifold endowed with a metric (Riemannian or Finslerian) to the non-metric setting. In order to apply the KCC theory we reformulate the Lorenz system as a set of two second order non-linear differential equations. The geometric invariants associated to this system (nonlinear and Berwald connections), and the deviation curvature tensor, as well as its eigenvalues, are explicitly obtained. The Jacobi stability of the equilibrium points of the Lorenz system is studied, and the condition of the stability of the equilibrium points is obtained. Finally, we consider the time evolution of the components of the deviation vector near the equilibrium points.


title: Numerical methods and comparison for the Dirac equation in the nonrelativistic limit regime

abstract: We analyze rigorously error estimates and compare numerically spatial/temporal resolution of various numerical methods for the discretization of the Dirac equation in the nonrelativistic limit regime, involving a small dimensionless parameter $0<\varepsilon\ll 1$ which is inversely proportional to the speed of light. In this limit regime, the solution is highly oscillatory in time, i.e. there are propagating waves with wavelength $O(\varepsilon^2)$ and $O(1)$ in time and space, respectively. We begin with several frequently used finite difference time domain (FDTD) methods and obtain rigorously their error estimates in the nonrelativistic limit regime by paying particular attention to how error bounds depend explicitly on mesh size $h$ and time step $\tau$ as well as the small parameter $\varepsilon$. Based on the error bounds, in order to obtain `correct' numerical solutions in the nonrelativistic limit regime, i.e. $0<\varepsilon\ll 1$, the FDTD methods share the same $\varepsilon$-scalability on time step: $\tau=O(\varepsilon^3)$. Then we propose and analyze two numerical methods for the discretization of the Dirac equation by using the Fourier spectral discretization for spatial derivatives combined with the exponential wave integrator and time-splitting technique for temporal derivatives, respectively. Rigorous error bounds for the two numerical methods show that their $\varepsilon$-scalability on time step is improved to $\tau=O(\varepsilon^2)$ when $0<\varepsilon\ll 1$. Extensive numerical results are reported to support our error estimates.


title: Quantitative Analysis of Whether Machine Intelligence Can Surpass Human Intelligence

abstract: Whether the machine intelligence can surpass the human intelligence is a controversial issue. On the basis of traditional IQ, this article presents the Universal IQ test method suitable for both the machine intelligence and the human intelligence. With the method, machine and human intelligences were divided into 4 major categories and 15 subcategories. A total of 50 search engines across the world and 150 persons at different ages were subject to the relevant test. And then, the Universal IQ ranking list of 2014 for the test objects was obtained.


title: Interaction of Alamethicin Pores in DMPC Bilayers

abstract: We have investigated the x-ray scattering signal of highly aligned multilayers of the zwitterionic lipid 1,2-dimyristoyl-sn-glycero-3-phosphatidylcholine containing pores formed by the antimicrobial peptide alamethicin as a function of the peptide/lipid ratio. We are able to obtain information on the structure factor of the pore fluid, which then yields the interaction potential between pores in the plane of the bilayers. Aside from a hard core with a radius corresponding to the geometric radius of the pore, we find a repulsive lipid-mediated interaction with a range of $\simeq 30$ {\AA} and a contact value of 2.4 $k_BT$. This result is in qualitative agreement with recent theoretical models.


title: Electrophoretic Mobility of Polyelectrolytes within a Confining Well

abstract: We present a numerical study of polyelectrolytes electrophoresing in free solution while squeezed by an axisymmetric confinement force transverse to their net displacement. Hybrid multi-particle collision dynamics and molecular dynamics simulations with mean-field finite Debye layers show that even though the polyelectrolyte chains remain "free-draining", their electrophoretic mobility increases with confinement in nanoconfining potential wells. The primary mechanism leading to the increase in mobility above the free-solution value, despite long-range hydrodynamic screening by counterion layers, is the orientation of polymer segments within Debye layers. The observed length-dependence of the electrophoretic mobility arises due to secondary effects of counterion condensation related to confinement compactification.


title: The derived non-commutative Poisson bracket on Koszul Calabi-Yau algebras

abstract: Let $A$ be a Koszul (or more generally, $N$-Koszul) Calabi-Yau algebra. Inspired by the works of Kontsevich, Ginzburg and Van den Bergh, we show that there is a derived non-commutative Poisson structure on $A$, which induces a graded Lie algebra structure on the cyclic homology of $A$; moreover, we show that the Hochschild homology of $A$ is a Lie module over the cyclic homology and the Connes long exact sequence is in fact a sequence of Lie modules. Finally, we show that the Leibniz-Loday bracket associated to the derived non-commutative Poisson structure on $A$ is naturally mapped to the Gerstenhaber bracket on the Hochschild cohomology of its Koszul dual algebra and hence on that of $A$ itself. Relations with some other brackets in literature are also discussed and several examples are given in detail.


title: Fundamental Relations Between Reactive and Proactive Relay-Selection Strategies

abstract: Two major relay-selection strategies widely applied in cooperative decode-and-forward (DF) relaying networks, namely, reactive relay selection (RRS) and proactive relay selection (PRS), are generally looked upon as independent and studied separately. In this paper, RRS and PRS are proven to be equivalent with respect to the end-to-end outage probability from the first principle, i.e. their respective relay-selection criteria. On the other hand, RRS is shown to be superior to PRS with respect to the end-to-end symbol error rate. Afterwards, a case study of a general DF relaying system, subject to co-channel interferences and additive white Gaussian noise at both the relaying nodes and the destination, is performed to explicitly illustrate the aforementioned outage equivalence. These fundamental relations provide intuitive yet insightful performance benchmarks for comparing various applications of these two relay-selection strategies.


title: Block-Maxima of Vines

abstract: We examine the dependence structure of finite block-maxima of multivariate distributions. We provide a closed form expression for the copula density of the vector of the block-maxima. Further, we show how partial derivatives of three-dimensional vine copulas can be obtained by only one-dimensional integration. Combining these results allows the numerical treatment of the block-maxima of any three-dimensional vine copula for finite block-sizes. We look at certain vine copula specifications and examine how the density of the block-maxima behaves for different block-sizes. Additionally, a real data example from hydrology is considered. In extreme-value theory for multivariate normal distributions, a certain scaling of each variable and the correlation matrix is necessary to obtain a non-trivial limiting distribution when the block-size goes to infinity. This scaling is applied to different three-dimensional vine copula specifications.


title: On Entropy Bumps for Calder\'on-Zygmund Operators

abstract: We study two weight inequalities in the recent innovative language of `entropy' due to Treil-Volberg. The inequalities are extended to $ L ^{p}$, for $ 1< p \neq 2 < \infty $, with new short proofs. A result proved is as follows. Let $ \varepsilon $ be a monotonic increasing function on $ (1, \infty)$ which satisfy $ \int _{1} ^{\infty} \frac {dt} {\varepsilon (t) t} = 1$. Let $ \sigma $ and $ w$ be two weights on $ \mathbb R ^{d}$. If this supremum is finite, for a choice of $ 1< p < \infty $, $$ \sup _{Q} \biggl[ \frac {\sigma (Q)} {\lvert Q\rvert} \biggr]^{p-1} \frac {\int _{Q} M (\sigma \chi_{Q})} {\sigma (Q)} \cdot \frac {w (Q)} {\lvert Q\rvert}\biggl[ \frac {\int _{Q} M (w \chi_{Q})} {w (Q)}\biggr]^{p-1} < \infty, $$ then any Calder\'on-Zygmund operator $ T$ satisfies the bound $ \lVert T _{\sigma} f \rVert _{L ^{p} (w)} \lesssim \lVert f\rVert _{L ^{p} (\sigma)} $.


title: Stellar Kinematics and Metallicities in the Ultra-Faint Dwarf Galaxy Reticulum II

abstract: We present Magellan/M2FS, VLT/GIRAFFE, and Gemini South/GMOS spectroscopy of the newly discovered Milky Way satellite Reticulum II. Based on the spectra of 25 Ret II member stars selected from Dark Energy Survey imaging, we measure a mean heliocentric velocity of 62.8 +/- 0.5 km/s and a velocity dispersion of 3.3 +/- 0.7 km/s. The mass-to-light ratio of Ret II within its half-light radius is 470 +/- 210 Msun/Lsun, demonstrating that it is a strongly dark matter-dominated system. Despite its spatial proximity to the Magellanic Clouds, the radial velocity of Ret II differs from that of the LMC and SMC by 199 and 83 km/s, respectively, suggesting that it is not gravitationally bound to the Magellanic system. The likely member stars of Ret II span 1.3 dex in metallicity, with a dispersion of 0.28 +/- 0.09 dex, and we identify several extremely metal-poor stars with [Fe/H] < -3. In combination with its luminosity, size, and ellipticity, these results confirm that Ret II is an ultra-faint dwarf galaxy. With a mean metallicity of [Fe/H] = -2.65 +/- 0.07, Ret II matches Segue~1 as the most metal-poor galaxy known. Although Ret II is the third-closest dwarf galaxy to the Milky Way, the line-of-sight integral of the dark matter density squared is log J = 18.8 +/- 0.6 Gev^2/cm^5 within 0.2 degrees, indicating that the predicted gamma-ray flux from dark matter annihilation in Ret II is lower than that of several other dwarf galaxies.


title: Error estimates for a numerical approximation to the compressible barotropic Navier-Stokes equations

abstract: We present here a general method based on the investigation of the relative energy of the system, that provides an unconditional error estimate for the approximate solution of the barotropic Navier Stokes equations obtained by time and space discretization. We use this methodology to derive an error estimate for a specific DG/finite element scheme for which the convergence was proved in [27]. This is an extended version of the paper submitted to IMAJNA.


title: A Regularized Newton Method for Computing Ground States of Bose-Einstein condensates

abstract: In this paper, we propose a regularized Newton method for computing ground states of Bose-Einstein condensates (BECs), which can be formulated as an energy minimization problem with a spherical constraint. The energy functional and constraint are discretized by either the finite difference, or sine or Fourier pseudospectral discretization schemes and thus the original infinite dimensional nonconvex minimization problem is approximated by a finite dimensional constrained nonconvex minimization problem. Then an initial solution is first constructed by using a feasible gradient type method, which is an explicit scheme and maintains the spherical constraint automatically. To accelerate the convergence of the gradient type method, we approximate the energy functional by its second-order Taylor expansion with a regularized term at each Newton iteration and adopt a cascadic multigrid technique for selecting initial data. It leads to a standard trust-region subproblem and we solve it again by the feasible gradient type method. The convergence of the regularized Newton method is established by adjusting the regularization parameter as the standard trust-region strategy. Extensive numerical experiments on challenging examples, including a BEC in three dimensions with an optical lattice potential and rotating BECs in two dimensions with rapid rotation and strongly repulsive interaction, show that our method is efficient, accurate and robust.


title: A short proof of the equivalence of left and right convergence for sparse graphs

abstract: There are several notions of convergence for sequences of bounded degree graphs. One such notion is left convergence, which is based on counting neighborhood distributions. Another notion is right convergence, based on counting homomorphisms to a target (weighted) graph. Borgs, Chayes, Kahn and Lov\'asz showed that a sequence of bounded degree graphs is left convergent if and only if it is right convergent for certain target graphs $H$ with all weights (including loops) close to $1$. We give a short alternative proof of this statement. In particular, for each bounded degree graph $G$ we associate functions $f_{G,k}$ for every positive integer $k$, and we show that left convergence of a sequence of graphs is equivalent to the convergence of the partial derivatives of each of these functions at the origin, while right convergence is equivalent to pointwise convergence. Using the bound on the maximum degree of the graphs, we can uniformly bound the partial derivatives at the origin, and show that the Taylor series converges uniformly on a domain independent of the graph, which implies the equivalence.


title: Re-examining sin(2beta) and Delta m(d) from evolution of B(d) mesons with decoherence

abstract: In the time evolution of neutral meson systems, a perfect quantum coherence is usually assumed. The important quantities of the B(d) system, such as sin (2beta) and Delta m(d), are determined under this assumption. However, the meson system interacts with its environment. This interaction can lead to decoherence in the mesons even before they decay. In our formalism this decoherence is modelled by a single parameter lambda. It is desirable to re-examine the procedures of determination of sin(2beta) and Delta m(d) in meson systems with decoherence. We find that the present values of these two quantities are modulated by lambda. Re-analysis of B(d) data from B-factories and LHCb can lead to a clean determination of lambda, sin(2beta) and Delta m(d).


title: Spin-charge separated quasiparticles in one dimensional quantum fluids

abstract: We revisit the problem of dynamical response in spin-charge separated one dimensional quantum fluids. In the framework of Luttinger liquid theory, the dynamical response is formulated in terms of noninteracting bosonic collective excitations carrying either charge or spin. We argue that, as a result of spectral nonlinearity, long-lived excitations are best understood in terms of generally strongly interacting fermionic holons and spinons. This has far reaching ramifications for the construction of mobile impurity models used to determine threshold singularities in response functions. We formulate and solve the appropriate mobile impurity model describing the spinon threshold in the single-particle Green's function. Our formulation further raises the question whether it is possible to realize a model of noninteracting fermionic holons and spinons in microscopic lattice models of interacting spinful fermions. We investigate this issue in some detail by means of density matrix renormalization group (DMRG) computations.


title: One-dimensional chain "melting" in incommensurate potassium

abstract: Between 19 and 54 GPa, potassium has a complex composite incommensurate host-guest structure which undergoes two intraphase transitions over this pressure range. The temperature dependence of these host-guest phases is further complicated by the onset of an order-disorder transition in their guest chains. Here, we report single crystal, quasi-single crystal, and powder synchrotron X-ray diffraction measurements of this order-disorder phenomenon in incommensurate potassium to 47 GPa and 750 K. The so-called chain "melting" transition is clearly visible over a 22 GPa pressure range, and there are significant changes in the slope of the phase boundary which divides the ordered and disordered phases, one of which results from the intraphase transitions in the guest structure.


title: Pricing and Risk Management with High-Dimensional Quasi Monte Carlo and Global Sensitivity Analysis

abstract: We review and apply Quasi Monte Carlo (QMC) and Global Sensitivity Analysis (GSA) techniques to pricing and risk management (greeks) of representative financial instruments of increasing complexity. We compare QMC vs standard Monte Carlo (MC) results in great detail, using high-dimensional Sobol' low discrepancy sequences, different discretization methods, and specific analyses of convergence, performance, speed up, stability, and error optimization for finite differences greeks. We find that our QMC outperforms MC in most cases, including the highest-dimensional simulations and greeks calculations, showing faster and more stable convergence to exact or almost exact results. Using GSA, we are able to fully explain our findings in terms of reduced effective dimension of our QMC simulation, allowed in most cases, but not always, by Brownian bridge discretization. We conclude that, beyond pricing, QMC is a very promising technique also for computing risk figures, greeks in particular, as it allows to reduce the computational effort of high-dimensional Monte Carlo simulations typical of modern risk management.


title: Accurate and efficient numerical methods for computing ground states and dynamics of dipolar Bose-Einstein condensates via the nonuniform FFT

abstract: In this paper, we propose efficient and accurate numerical methods for computing the ground state and dynamics of the dipolar Bose-Einstein condensates utilising a newly developed dipole-dipole interaction (DDI) solver that is implemented with the non-uniform fast Fourier transform (NUFFT) algorithm. We begin with the three-dimensional (3D) Gross-Pitaevskii equation (GPE) with a DDI term and present the corresponding two-dimensional (2D) model under a strongly anisotropic confining potential. Different from existing methods, the NUFFT based DDI solver removes the singularity by adopting the spherical/polar coordinates in Fourier space in 3D/2D, respectively, thus it can achieve spectral accuracy in space and simultaneously maintain high efficiency by making full use of FFT and NUFFT whenever it is necessary and/or needed. Then, we incorporate this solver into existing successful methods for computing the ground state and dynamics of GPE with a DDI for dipolar BEC. Extensive numerical comparisons with existing methods are carried out for computing the DDI, ground states and dynamics of the dipolar BEC. Numerical results show that our new methods outperform existing methods in terms of both accuracy and efficiency.


title: Recent advances in percolation theory and its applications

abstract: Percolation is the simplest fundamental model in statistical mechanics that exhibits phase transitions signaled by the emergence of a giant connected component. Despite its very simple rules, percolation theory has successfully been applied to describe a large variety of natural, technological and social systems. Percolation models serve as important universality classes in critical phenomena characterized by a set of critical exponents which correspond to a rich fractal and scaling structure of their geometric features. In this review we will first outline the basic features of the ordinary model and take a glimpse at a number of selective variations and modifications of the original model. Directed percolation process will be also discussed as a prototype of systems displaying a nonequilibrium phase transition. After a short review on SLE, we will provide an overview on existence of the scaling limit and conformal invariance of the critical percolation. We will also establish a connection with the magnetic models. Recent applications of the percolation theory in natural and artificial landscapes are also reviewed.


title: On absorption in semigroups and $n$-ary semigroups

abstract: The notion of absorption was developed a few years ago by Barto and Kozik and immediately found many applications, particularly in topics related to the constraint satisfaction problem. We investigate the behavior of absorption in semigroups and n-ary semigroups (that is, algebras with one n-ary associative operation). In the case of semigroups, we give a simple necessary and sufficient condition for a semigroup to be absorbed by its subsemigroup. We then proceed to n-ary semigroups, where we conjecture an analogue of this necessary and sufficient condition, and prove that the conjectured condition is indeed necessary and sufficient for B to absorb A (where A is an n-ary semigroup and B is its n-ary subsemigroup) in the following three cases: when A is commutative, when |A-B|=1 and when A is an idempotent ternary semigroup.


title: The Dark Energy Camera

abstract: The Dark Energy Camera is a new imager with a 2.2-degree diameter field of view mounted at the prime focus of the Victor M. Blanco 4-meter telescope on Cerro Tololo near La Serena, Chile. The camera was designed and constructed by the Dark Energy Survey Collaboration, and meets or exceeds the stringent requirements designed for the wide-field and supernova surveys for which the collaboration uses it. The camera consists of a five element optical corrector, seven filters, a shutter with a 60 cm aperture, and a CCD focal plane of 250 micron thick fully-depleted CCDs cooled inside a vacuum Dewar. The 570 Mpixel focal plane comprises 62 2kx4k CCDs for imaging and 12 2kx2k CCDs for guiding and focus. The CCDs have 15 microns x15 microns pixels with a plate scale of 0.263 arc sec per pixel. A hexapod system provides state-of-the-art focus and alignment capability. The camera is read out in 20 seconds with 6-9 electrons readout noise. This paper provides a technical description of the camera's engineering, construction, installation, and current status.


title: Work measurement in an optomechanical quantum heat engine

abstract: We analyze theoretically the measurement of the mean output work and its fluctuations in a recently proposed optomechanical quantum heat engine [K. Zhang {\it et al.} Phys. Rev. Lett. {\bf112}, 150602 (2014)]. After showing that this work can be evaluated by a continuous measurements of the intracavity photon number we discuss both dispersive and absorptive measurement schemes and analyze their back-action effects on the efficiency of the engine. Both measurements are found to reduce the efficiency of the engine, but their back-action is both qualitatively and quantitatively different. For dispersive measurements the efficiency decreases as a result of the mixing of photonic and phononic excitations, while for absorptive measurements, its reduction results from photon losses due to the interaction with the quantum probe.


title: Gradual Training Method for Denoising Auto Encoders

abstract: Stacked denoising auto encoders (DAEs) are well known to learn useful deep representations, which can be used to improve supervised training by initializing a deep network. We investigate a training scheme of a deep DAE, where DAE layers are gradually added and keep adapting as additional layers are added. We show that in the regime of mid-sized datasets, this gradual training provides a small but consistent improvement over stacked training in both reconstruction quality and classification error over stacked training on MNIST and CIFAR datasets.


title: A singular limit problem for rotating capillary fluids with variable rotation axis

abstract: In the present paper we study a singular perturbation problem for a Navier-Stokes-Korteweg model with Coriolis force. Namely, we perform the incompressible and fast rotation asymptotics simultaneously, while we keep the capillarity coefficient constant in order to capture surface tension effects in the limit. We consider here the case of variable rotation axis: we prove the convergence to a linear parabolic-type equation with variable coefficients. The proof of the result relies on compensated compactness arguments. Besides, we look for minimal regularity assumptions on the variations of the axis.


title: Large Prandtl Number Asymptotics in Randomly Forced Turbulent Convection

abstract: We establish the convergence of statistically invariant states for the stochastic Boussinesq Equations in the infinite Prandtl number limit and in particular demonstrate the convergence of the Nusselt number (a measure of heat transport in the fluid). This is a singular parameter limit significant in mantle convection and for gasses under high pressure. The equations are subject to a both temperature gradient on the boundary and internal heating in the bulk driven by a stochastic, white in time, gaussian forcing. Here, the stochastic source terms have a strong physical motivation for example as a model of radiogenic heating. Our approach uses mixing properties of the formal limit system to reduce the convergence of invariant states to an analysis of the finite time asymptotics of solutions and parameter-uniform moment bounds. Here, it is notable that there is a phase space mismatch between the finite Prandtl system and the limit equation, and we implement methods to lift both finite and infinite time convergence results to an extended phase space which includes velocity fields. For the infinite Prandtl stochastic Boussinesq equations, we show that the associated invariant measure is unique and that the dual Markovian dynamics are contractive in an appropriate Kantorovich-Wasserstein metric. We then address the convergence of solutions on finite time intervals, which is still a singular perturbation. In the process we derive well-posed equations which accurately approximate the dynamics up to the initial time when the Prandtl number is large.


title: The Roles of Reconnected Flux and Overlying Fields in CME Speeds

abstract: Researchers have reported (i) correlations of coronal mass ejection (CME) speeds and the total photospheric magnetic flux swept out by flare ribbons in flare-associated eruptive events, and, separately, (ii) correlations of CME speeds and more rapid decay, with height, of magnetic fields in potential coronal models above eruption sites. Here, we compare the roles of both ribbon fluxes and the decay rates of overlying fields in a set of 16 eruptive events. We confirm previous results that higher CME speeds are associated with both larger ribbon fluxes and more rapidly decaying overlying fields. We find the association with ribbon fluxes to be weaker than a previous report, but stronger than the dependence on the decay rate of overlying fields. Since the photospheric ribbon flux is thought to approximate the amount of coronal magnetic flux reconnected during the event, the correlation of speeds with ribbon fluxes suggests that reconnection plays some role in accelerating CMEs. One possibility is that reconnected fields that wrap around the rising ejection produce an increased outward hoop force, thereby increasing CME acceleration. The correlation of CME speeds with more rapidly decaying overlying fields might be caused by greater downward magnetic tension in stronger overlying fields, which could act as a source of drag on rising ejections.


title: Skew ray tracing in a step-index optical fiber using Geometric Algebra

abstract: We used Geometric Algebra to compute the paths of skew rays in a cylindrical, step-index multimode optical fiber. To do this, we used the vector addition form for the law of propagation, the exponential of an imaginary vector form for the law of refraction, and the juxtaposed vector product form for the law of reflection. In particular, the exponential forms of the vector rotations enables us to take advantage of the addition or subtraction of exponential arguments of two rotated vectors in the derivation of the ray tracing invariants in cylindrical and spherical coordinates. We showed that the light rays inside the optical fiber trace a polygonal helical path characterized by three invariants that relate successive reflections inside the fiber: the ray path distance, the difference in axial distances, and the difference in the azimuthal angles. We also rederived the known generalized formula for the numerical aperture for skew rays, which simplifies to the standard form for meridional rays.


title: Quantum transport in ultracold atoms

abstract: Ultracold atoms confined by engineered magnetic or optical potentials are ideal systems for studying phenomena otherwise difficult to realize or probe in the solid state because their atomic interaction strength, number of species, density, and geometry can be independently controlled. This review focuses on quantum transport phenomena in atomic gases that mirror and oftentimes either better elucidate or show fundamental differences with those observed in mesoscopic and nanoscopic systems. We discuss significant progress in performing transport experiments in atomic gases, contrast similarities and differences between transport in cold atoms and in condensed matter systems, and survey inspiring theoretical predictions that are difficult to verify in conventional setups. These results further demonstrate the versatility offered by atomic systems in the study of nonequilibrium phenomena and their promise for novel applications.


title: Superconducting Circuitry for Quantum Electromechanical Systems

abstract: Superconducting systems have a long history of use in experiments that push the frontiers of mechanical sensing. This includes both applied and fundamental research, which at present day ranges from quantum computing research and efforts to explore Planck-scale physics to fundamental studies on the nature of motion and the quantum limits on our ability to measure it. In this paper, we first provide a short history of the role of superconducting circuitry and devices in mechanical sensing, focusing primarily on efforts in the last decade to push the study of quantum mechanics to include motion on the scale of human-made structures. This background sets the stage for the remainder of the paper, which focuses on the development of quantum electromechanical systems (QEMS) that incorporate superconducting quantum bits (qubits), superconducting transmission line resonators and flexural nanomechanical elements. In addition to providing the motivation and relevant background on the physical behavior of these systems, we discuss our recent efforts to develop a particular type of QEMS that is based upon the Cooper-pair box (CPB) and superconducting coplanar waveguide (CPW) cavities, a system which has the potential to serve as a testbed for studying the quantum properties of motion in engineered systems.


title: Counting designs

abstract: We give estimates on the number of combinatorial designs, which prove (and generalise) a conjecture of Wilson from 1974 on the number of Steiner Triple Systems. This paper also serves as an expository treatment of our recently developed method of Randomised Algebraic Construction: we give a simpler proof of a special case of our result on clique decompositions of hypergraphs, namely triangle decompositions of quasirandom graphs.


title: Extending four dimensional Ricci flows with bounded scalar curvature

abstract: We consider smooth solutions (M,g(t)), 0 <= t <T, to Ricci flow on compact, connected, four dimensional manifolds without boundary. We assume that the scalar curvature is bounded uniformly, and that T is finite. In this case, we show that the metric space (M,d(t)) associated to (M,g(t)) converges uniformly in the C^0 sense to (X,d), as t approaches T, where (X,d) is a C^0 Riemannian orbifold with at most finitely many orbifold points. Estimates on the rate of convergence near and away from the orbifold points are given. We also show that it is possible to continue the flow past (X,d) using the orbifold Ricci flow.


title: Minimum Distance Approach to Inference with Many Instruments

abstract: I analyze a linear instrumental variables model with a single endogenous regressor and many instruments. I use invariance arguments to construct a new minimum distance objective function. With respect to a particular weight matrix, the minimum distance estimator is equivalent to the random effects estimator of Chamberlain and Imbens (2004), and the estimator of the coefficient on the endogenous regressor coincides with the limited information maximum likelihood estimator. This weight matrix is inefficient unless the errors are normal, and I construct a new, more efficient estimator based on the optimal weight matrix. Finally, I show that when the minimum distance objective function does not impose a proportionality restriction on the reduced-form coefficients, the resulting estimator corresponds to a version of the bias-corrected two-stage least squares estimator. I use the objective function to construct confidence intervals that remain valid when the proportionality restriction is violated.


title: Physical Aspects of Unitary evolution of Bianchi-I Quantum Cosmological Model

abstract: In this work, we study some physical aspects of unitary evolution of Bianchi-I model. In particular, we study the behavior of the volume and the scale factor as a function of time for the Bianchi-I universe with ultra-relativistic fluid ($\alpha=1$). The expectation value of volume is shown not to hit any singularity. We elucidate on the anisotropic nature of the solution and physically interpret the wavefunction as a superposition of collapsing universe and expanding universe mimicking Hartle-Hawking type wavefunction. The same analysis has been done for $\alpha\neq 1$ as well. The work also comments on how to obtain Vilenkin type wavefunction representing an expanding universe without that being superposed with a solution representing a collapsing universe. We also point at a possible phase transition phenomenon for $\alpha\neq1$ case with existence of a critical parameter.


title: A pairwise likelihood approach to simultaneous clustering and dimensional reduction of ordinal data

abstract: The literature on clustering for continuous data is rich and wide; differently, that one developed for categorical data is still limited. In some cases, the problem is made more difficult by the presence of noise variables/dimensions that do not contain information about the clustering structure and could mask it. The aim of this paper is to propose a model for simultaneous clustering and dimensionality reduction of ordered categorical data able to detect the discriminative dimensions discarding the noise ones. Following the underlying response variable approach, the observed variables are considered as a discretization of underlying first-order latent continuous variables distributed as a Gaussian mixture. To recognize discriminative and noise dimensions, these variables are considered to be linear combinations of two independent sets of second-order latent variables where only one contains the information about the cluster structure while the other contains noise dimensions. The model specification involves multidimensional integrals that make the maximum likelihood estimation cumbersome and in some cases infeasible. To overcome this issue the parameter estimation is carried out through an EM-like algorithm maximizing a pairwise log-likelihood. Examples of application of the model on real and simulated data are performed to show the effectiveness of the proposal.


title: Representing numeric data in 32 bits while preserving 64-bit precision

abstract: Data files often consist of numbers having only a few significant decimal digits, whose information content would allow storage in only 32 bits. However, we may require that arithmetic operations involving these numbers be done with 64-bit floating-point precision, which precludes simply representing the data as 32-bit floating-point values. Decimal floating point gives a compact and exact representation, but requires conversion with a slow division operation before it can be used. Here, I show that interesting subsets of 64-bit floating-point values can be compactly and exactly represented by the 32 bits consisting of the sign, exponent, and high-order part of the mantissa, with the lower-order 32 bits of the mantissa filled in by table lookup, indexed by bits from the part of the mantissa retained, and possibly from the exponent. For example, decimal data with 4 or fewer digits to the left of the decimal point and 2 or fewer digits to the right of the decimal point can be represented in this way using the lower-order 5 bits of the retained part of the mantissa as the index. Data consisting of 6 decimal digits with the decimal point in any of the 7 positions before or after one of the digits can also be represented this way, and decoded using 19 bits from the mantissa and exponent as the index. Encoding with such a scheme is a simple copy of half the 64-bit value, followed if necessary by verification that the value can be represented, by checking that it decodes correctly. Decoding requires only extraction of index bits and a table lookup. Lookup in a small table will usually reference cache; even with larger tables, decoding is still faster than conversion from decimal floating point with a division operation. I discuss how such schemes perform on recent computer systems, and how they might be used to automatically compress large arrays in interpretive languages such as R.


title: Solution of the Hyperon Puzzle within a Relativistic Mean-Field Model

abstract: The equation of state of cold baryonic matter is studied within a relativistic mean-field model with hadron masses and coupling constants depending on the scalar field. All hadron masses undergo a universal scaling, whereas the coupling constants are scaled differently. The appearance of hyperons in dense neutron star interiors is accounted for, however the equation of state remains sufficiently stiff if a reduction of the $\phi$ meson mass is included. Our equation of state matches well the constraints known from analyses of the astrophysical data and the particle production in heavy-ion collisions.


title: Perturbative study of the QCD phase diagram for heavy quarks at nonzero chemical potential

abstract: We investigate the phase diagram of QCD with heavy quarks at finite temperature and chemical potential in the context of background field methods. In particular, we use a massive extension of the Landau-DeWitt gauge which is motivated by previous studies of the deconfinement phase transition in pure Yang-Mills theories. We show that a simple one-loop calculation is able to capture the richness of the phase diagram in the heavy quark region, both at real and imaginary chemical potential. Moreover, dimensionless ratios of quantities directly measurable in numerical simulations are in good agreement with lattice results.


title: Tidal Truncation of Inclined Circumstellar and Circumbinary Discs in Young Stellar Binaries

abstract: Recent observations have shown that circumstellar and circumbinary discs in young stellar binaries are often misaligned with respect to the binary orbital plane. We analyze the tidal truncation of such misaligned discs due to torques applied to the disc at the Lindblad resonances from the tidal forcings of the binary. We consider eccentric binaries with arbitrary binary-disc inclination angles. We determine the dependence of the tidal forcing strengths on the binary parameters and show that they are complicated non-monotonic functions of eccentricity and inclination. We adopt a truncation criterion determined by the balance between resonant torque and viscous torque, and use it to calculate the outer radii of circumstellar discs and the inner radii of circumbinary discs. Misaligned circumstellar discs have systematically larger outer radii than aligned discs, and are likely to fill their Roche lobes if inclined by more than $45^\circ - 90^\circ$, depending on the binary mass ratio and disc viscosity parameter. Misaligned circumbinary discs generally have smaller inner radii than aligned discs, but the details depend sensitively on the binary and disc parameters.


title: Global regularity criterion for the 3D Navier-Stokes equations with large data

abstract: In this paper, we study the global regularity of strong solution to the Cauchy problem of 3D incompressible Navier-Stokes equations with large data and non-zero force. We prove that the strong solution exists globally for $\nabla u\in L^{\alpha}\left( 0,T,L^{2}\left( \Omega\right) \right) $ with $2\leq\alpha\leq4$.


title: Propagation peculiarities of mean field massive gravity

abstract: Massive gravity (mGR) describes a dynamical "metric" on a fiducial, background one. We investigate fluctuations of the dynamics about mGR solutions, that is about its "mean field theory". Analyzing mean field massive gravity propagation characteristics is not only equivalent to studying those of the full non-linear theory, but also in direct correspondence with earlier analyses of charged higher spin systems, the oldest example being the charged, massive spin 3/2 Rarita--Schwinger (RS) theory. The fiducial and mGR mean field background metrics in the mean field mGR model correspond to the RS Minkowski metric and external EM field. The common implications in both systems are that hyperbolicity holds only in a weak background-mean-field limit, immediately ruling both theories out as fundamental theories. Although both can still be considered as predictive effective models in the weak regime, their lower helicities exhibit superluminal behavior: lower helicity gravitons are superluminal as compared to photons propagating on either the fiducial or background metric. This "crystal-like" phenomenon of differing helicities having differing propagation speeds in both mean field mGR and mGR is a peculiar feature of these models.


title: The Donaldson-Thomas theory of $K3\times E$ via the topological vertex

abstract: Oberdieck and Pandharipande conjectured that the curve counting invariants of $S\times E$, the product of a $K3$ surface and an elliptic curve, is given by minus the reciprocal of the Igusa cusp form of weight 10. For a fixed primitive curve class in $S$ of square $2h-2$, their conjecture predicts that the corresponding partition functions are given by meromorphic Jacobi forms of weight $-10$ and index $h-1$. We calculate the partition functions for primitive classes of square -2 and of square 0. Our computation uses reduced Donaldson-Thomas invariants which are defined as the Behrend function weighted Euler characteristics of the quotient of the Hilbert scheme of curves in $S\times E$ by the action of $E$. Our technique is a mixture of motivic and toric methods (developed with Martijn Kool) which allows us to express the partition functions in terms of the topological vertex and subsequently in terms of Jacobi forms. We compute the partition functions for both Behrend function weighted Euler characteristics and for unweighted Euler characteristics. The results for the Behrend function weighted case depends on Conjecture 18 from https://arxiv.org/abs/1608.07369


title: Antenna Array Signal Processing for Quaternion-Valued Wireless Communication Systems

abstract: Quaternion-valued wireless communication systems have been studied in the past. Although progress has been made in this promising area, a crucial missing link is lack of effective and efficient quaternion-valued signal processing algorithms for channel equalisation and beamforming. With most recent developments in quaternion-valued signal processing, in this work, we fill the gap to solve the problem and further derive the quaternion-valued Wiener solution for block-based calculation.


title: A New Companion to Capparelli's Identities

abstract: We discuss a new companion to Capparelli's identities. Capparelli's identities for m=1,2 state that the number of partitions of $n$ into distinct parts not congruent to m, -m modulo $6$ is equal to the number of partitions of n into distinct parts not equal to m, where the difference between parts is greater than or equal to 4, unless consecutive parts are either both consecutive multiples of 3 or add up to to a multiple of 6. In this paper we show that the set of partitions of n into distinct parts where the odd-indexed parts are not congruent to m modulo 3, the even-indexed parts are not congruent to -m modulo 3, and 3l+1 and 3l+2 do not appear together as consecutive parts for any integer l has the same number of elements as the above mentioned Capparelli's partitions of n. In this study we also extend the work of Alladi, Andrews and Gordon by providing a complete set of generating functions for the refined Capparelli partitions, and conjecture some combinatorial inequalities.


title: Compressed Sensing Recovery via Nonconvex Shrinkage Penalties

abstract: The $\ell^0$ minimization of compressed sensing is often relaxed to $\ell^1$, which yields easy computation using the shrinkage mapping known as soft thresholding, and can be shown to recover the original solution under certain hypotheses. Recent work has derived a general class of shrinkages and associated nonconvex penalties that better approximate the original $\ell^0$ penalty and empirically can recover the original solution from fewer measurements. We specifically examine p-shrinkage and firm thresholding. In this work, we prove that given data and a measurement matrix from a broad class of matrices, one can choose parameters for these classes of shrinkages to guarantee exact recovery of the sparsest solution. We further prove convergence of the algorithm iterative p-shrinkage (IPS) for solving one such relaxed problem.


title: The Constrained Krasnosel'skii Formula for Parabolic Differential Inclusions

abstract: We consider a constrained evolution inclusions of parabolic type \eqref{inkluzja-rozn} involving an $m$-dissipative linear operator and the source term of multivalued type in a Banach space and topological properties of the solution map. We show a relation between the constrained fixed point index of the Krasnosel'skii--Poincar\'{e} operator of translation along trajectories associated with \eqref{inkluzja-rozn} and the appropriately defined constrained degree of $A + F\le 0 , \cdot \pr $ of the right-hand side in \eqref{inkluzja-rozn}. Our results extend those of \cite{cw} and \cite{gab-krysz}.


title: On the vanishing of local cohomology of the absolute integral closure in positive characteristic

abstract: The aim of this paper is to extend the main result of C. Huneke and G. Lyubeznik in [Adv. Math. 210 (2007), 498--504] to the class of rings that are images of Cohen-Macaulay local rings. Namely, let $R$ be a local Noetherian domain of positive characteristic that is an image of a Cohen-Macaulay local ring. We prove that all local cohomology of $R$ (below the dimension) maps to zero in a finite extension of the ring. As a direct consequence we obtain that the absolute integral closure of $R$ is a big Cohen-Macaulay algebra. Since every excellent local ring is an image of a Cohen-Macaulay local ring, this result is a generalization of the main result of M. Hochster and Huneke in [Ann. of Math. 135 (1992), 45--79] with a simpler proof.


title: Strategic Interaction Among Different Entities in Internet of Things

abstract: The economic model of the Internet of Things (IoT) consists of end users, advertisers and three different kinds of providers--IoT service provider (IoTSP), Wireless service provider (WSP) and cloud service provider (CSP). We investigate three different kinds of interactions among the providers. First, we consider that the IoTSP prices a bundled service to the end-users, and the WSP and CSP pay the IoTSP (push model). Next, we consider the model where the end-users independently pay the each provider (pull model). Finally, we consider a hybrid model of the above two where the IoTSP and WSP quote their prices to the end-users, but the CSP quotes its price to the IoTSP. We characterize and quantify the impact of the advertisement revenue on the equilibrium pricing strategy and payoff of providers, and corresponding demands of end users in each of the above interaction models. Our analysis reveals that the demand of end-users, and the payoffs of the providers are non decreasing functions of the advertisement revenue. For sufficiently high advertisement revenue, the IoTSP will offer its service free of cost in each interaction model. However, the payoffs of the providers, and the demand of end-users vary across different interaction models. Our analysis shows that the demand of end-users, and the payoff of the WSP are the highest in the pull (push, resp.) model in the low (high, resp.) advertisement revenue regime. The payoff of the IoTSP is always higher in the pull model irrespective of the advertisement revenue. The payoff of the CSP is the highest in the hybrid model in the low advertisement revenue regime. However, in the high advertisement revenue regime the payoff of the CSP in the hybrid model or in the push model can be higher depending on the equilibrium chosen in the push model.


title: The Origin of Inference: Ediacaran Ecology and the Evolution of Bayesian Brains

abstract: The evolution of spiking neurons and nervous systems in the late Ediacaran period simultaneously with the evolution of carnivores around 550 million years ago can be explained by the need for accurately timed decisions under an imminent threat of being eaten. A simple model shows that threshold triggering devices, spiking neurons, are utility-maximizing decision-makers for the timing of escape reflexes given the sensory cues available to Ediacaran animals at the onset of carnivory. Decisions are suboptimal for very weak stimuli, providing selection pressure for secondary processing of primary spike train data. A simple network can make approximately Bayes optimal decisions given stochastic spike trains. Decisions that are arbitrarily close to Bayes optimal can be obtained by enlarging this network. A subnetwork that computes the Bayesian posterior density of the critical state variable - distance between predator and prey - emerges as a core component of the decision-making mechanism. This is a neural analog of a Bayesian particle filter with cerebellar-like architecture. The model explains fundamental properties of neurons and nervous systems in modern animals and makes testable predictions.


title: Spherical quadrilaterals with three non-integer angles

abstract: We classify spherical quadrilaterals up to isometry in the case when one inner angle is a multiple of pi while the other three are not. This is equivalent to classification of Heun's equations with real parameters and one apparent singularity such that the monodromy consists of unitary transformations.


title: Lattice mismatch induced ripples and wrinkles in planar graphene/boron nitride superlattices

abstract: A continuum theory to describe periodic ripple formation in planar graphene/boron nitride superlattices is formulated. Due to the lattice mismatch between the two materials, it is shown that flat superlattices are unstable with respect to ripple formation of appropriate wavelengths. A competition between bending energy and transverse stretching energy gives rise to an optimal ripple wavelength that depends on the superlattice pitch. The optimal wavelengths predicted by the continuum theory are in good agreement with atomic scale total energy calculations previously reported in Nandwana and Ertekin, Nano Lett. 15 1468 (2015).


title: Knowledge reduction of dynamic covering decision information systems with varying attribute values

abstract: Knowledge reduction of dynamic covering information systems involves with the time in practical situations. In this paper, we provide incremental approaches to computing the type-1 and type-2 characteristic matrices of dynamic coverings because of varying attribute values. Then we present incremental algorithms of constructing the second and sixth approximations of sets by using characteristic matrices. We employ experimental results to illustrate that the incremental approaches are effective to calculate approximations of sets in dynamic covering information systems. Finally, we perform knowledge reduction of dynamic covering information systems with the incremental approaches.


title: Generalized Correntropy for Robust Adaptive Filtering

abstract: As a robust nonlinear similarity measure in kernel space, correntropy has received increasing attention in domains of machine learning and signal processing. In particular, the maximum correntropy criterion (MCC) has recently been successfully applied in robust regression and filtering. The default kernel function in correntropy is the Gaussian kernel, which is, of course, not always the best choice. In this work, we propose a generalized correntropy that adopts the generalized Gaussian density (GGD) function as the kernel (not necessarily a Mercer kernel), and present some important properties. We further propose the generalized maximum correntropy criterion (GMCC), and apply it to adaptive filtering. An adaptive algorithm, called the GMCC algorithm, is derived, and the mean square convergence performance is studied. We show that the proposed algorithm is very stable and can achieve zero probability of divergence (POD). Simulation results confirm the theoretical expectations and demonstrate the desirable performance of the new algorithm.


title: Binary astrometric microlensing with Gaia

abstract: We investigate whether Gaia can specify the binary fractions of massive stellar populations in the Galactic disk through astrometric microlensing. Furthermore, we study if some information about their mass distributions can be inferred via this method. In this regard, we simulate the binary astrometric microlensing events due to massive stellar populations according to the Gaia observing strategy by considering (a) stellar-mass black holes, (b) neutron stars, (c) white dwarfs and (d) main-sequence stars as microlenses. The Gaia efficiency for detecting the binary signatures in binary astrometric microlensing events is $\sim 10-20$ per cent. By calculating the optical depth due to the mentioned stellar populations, the number of the binary astrometric microlensing events being observed with Gaia with detectable binary signatures, for the binary fraction about 0.1, is estimated as 6, 11, 77 and 1316 respectively. Consequently, Gaia can potentially specify the binary fractions of these massive stellar populations. However, the binary fraction of black holes measured with this method has the large uncertainty owing to a low number of the estimated events. Knowing the binary fractions in massive stellar populations helps for studying the gravitational waves. Moreover, we investigate the number of massive microlenses which Gaia specifies their masses through astrometric microlensing of single lenses toward the Galactic bulge. The resulted efficiencies of measuring the mass of mentioned populations are 9.8, 2.9, 1.2 and 0.8 per cent respectively. The number of their astrometric microlensing events being observed in the Gaia era in which the lens mass can be inferred with the relative error less than 0.5 toward the Galactic bulge is estimated as 45, 34, 76 and 786 respectively.


title: A treatment procedure for VLT/SINFONI data cubes: application to NGC 5643

abstract: In this second paper of a series, we present a treatment procedure for data cubes obtained with the Spectrograph for Integral Field Observations in the Near Infrared of the Very Large Telescope. We verified that the treatment procedure improves significantly the quality of the images of the data cubes, allowing a more detailed analysis. The images of the Br$\gamma$ and H$_2 \lambda 21218$ emission lines from the treated data cube of the nuclear region of NGC 5643 reveal the existence of ionized and molecular-gas clouds around the nucleus, which cannot be seen clearly in the images from the non-treated data cube of this galaxy. The ionized-gas clouds represent the narrow-line region, in the form of a bicone. We observe a good correspondence between the positions of the ionized-gas clouds in the Br$\gamma$ image and in an [O III] image, obtained with the Hubble Space Telescope, of the nuclear region of this galaxy convolved with an estimate of the point-spread function of the data cube of NGC 5643. The morphologies of the ionized and molecular gas seem to be compatible with the existence of a molecular torus/disc that collimates the active galactic nucleus (AGN) emission. The molecular gas may also flow along this torus/disc, feeding the AGN. This scenario is compatible with the unified model for AGNs.


title: Quadratic unitary Cayley graphs of finite commutative rings

abstract: The purpose of this paper is to study spectral properties of a family of Cayley graphs on finite commutative rings. Let $R$ be such a ring and $R^\times$ its set of units. Let $Q_R=\{u^2: u\in R^\times\}$ and $T_R=Q_R\cup(-Q_R)$. We define the quadratic unitary Cayley graph of $R$, denoted by $\mathcal{G}_R$, to be the Cayley graph on the additive group of $R$ with respect to $T_R$; that is, $\mathcal{G}_R$ has vertex set $R$ such that $x, y \in R$ are adjacent if and only if $x-y\in T_R$. It is well known that any finite commutative ring $R$ can be decomposed as $R=R_1\times R_2\times\cdots\times R_s$, where each $R_i$ is a local ring with maximal ideal $M_i$. Let $R_0$ be a local ring with maximal ideal $M_0$ such that $|R_0|/|M_0| \equiv 3\,(\mod\,4)$. We determine the spectra of $\mathcal{G}_R$ and $\mathcal{G}_{R_0\times R}$ under the condition that $|R_i|/|M_i|\equiv 1\,(\mod\,4)$ for $1 \le i \le s$. We compute the energies and spectral moments of such quadratic unitary Cayley graphs, and determine when such a graph is hyperenergetic or Ramanujan.


title: Optimal Multiple Testing Under a Gaussian Prior on the Effect Sizes

abstract: We develop a new method for frequentist multiple testing with Bayesian prior information. Our procedure finds a new set of optimal p-value weights called the Bayes weights. Prior information is relevant to many multiple testing problems. Existing methods assume fixed, known effect sizes available from previous studies. However, the case of uncertain information is usually the norm. For a Gaussian prior on effect sizes, we show that finding the optimal weights is a non-convex problem. Despite the non-convexity, we give an efficient algorithm that solves this problem nearly exactly. We show that our method can discover new loci in genome-wide association studies. On several data sets it compares favorably to other methods. Open source code is available.


title: Automated Transient Identification in the Dark Energy Survey

abstract: We describe an algorithm for identifying point-source transients and moving objects on reference-subtracted optical images containing artifacts of processing and instrumentation. The algorithm makes use of the supervised machine learning technique known as Random Forest. We present results from its use in the Dark Energy Survey Supernova program (DES-SN), where it was trained using a sample of 898,963 signal and background events generated by the transient detection pipeline. After reprocessing the data collected during the first DES-SN observing season (Sep. 2013 through Feb. 2014) using the algorithm, the number of transient candidates eligible for human scanning decreased by a factor of 13.4, while only 1 percent of the artificial Type Ia supernovae (SNe) injected into search images to monitor survey efficiency were lost, most of which were very faint events. Here we characterize the algorithm's performance in detail, and we discuss how it can inform pipeline design decisions for future time-domain imaging surveys, such as the Large Synoptic Survey Telescope and the Zwicky Transient Facility. An implementation of the algorithm and the training data used in this paper are available at http://portal.nersc.gov/project/dessn/autoscan.


title: Uniform hyperbolicity for curve graphs of non-orientable surfaces

abstract: Hensel-Przytycki-Webb proved that all curve graphs of orientable surfaces are 17-hyperbolic. In this paper, we show that curve graphs of non-orientable surfaces are 17-hyperbolic by applying Hensel-Przytycki-Webb's argument. We also show that arc graphs of non-orientable surfaces are 7-hyperbolic, and arc-curve graphs of (non-)orientable surfaces are 9-hyperbolic.


title: ExoMol molecular line lists: IX The spectrum of AlO

abstract: Accurate line lists are calculated for aluminium monoxide covering the pure rotation, rotation-vibration and electronic (B -- X blue-green and A -- X infrared bands) spectrum. Line lists are presented for the main isotopologue, $^{27}$Al$^{16}$O, as well as for $^{27}$Al$^{17}$O, $^{27}$Al$^{18}$O and $^{26}$Al$^{16}$O. These line lists are suitable for high temperatures (up to 8000 K) including those relevant to exoplanetary atmospheres and cool stars. A combination of empirical and \textit{ab initio} methods is used: the potential energy curves were previously determined to high accuracy by fitting to extensive data from analysis of laboratory spectra; a high quality {\it ab initio} dipole moment curve is calculated using quadruple zeta basis set and the multi-reference configuration interaction (MRCI) method. Partition functions plus full line lists of transitions are made available in an electronic form as supplementary data to this article and at \url{www.exomol.com}.


title: On Elliptic Systems involving critical Hardy-Sobolev exponents (Part II)

abstract: This paper is the second part of a work devoted to the study of elliptic systems involving multiple Hardy-Sobolev critical exponents: $$\begin{cases} -\Delta u-\lambda \frac{|u|^{2^*(s_1)-2}u}{|x|^{s_1}}=\kappa\alpha \frac{1}{|x|^{s_2}}|u|^{\alpha-2}u|v|^\beta\quad &\hbox{in}\;\Omega,\\ -\Delta v-\mu \frac{|v|^{2^*(s_1)-2}v}{|x|^{s_1}}=\kappa\beta \frac{1}{|x|^{s_2}}|u|^{\alpha}|v|^{\beta-2}v\quad &\hbox{in}\;\Omega,\\ \kappa>0,(u,v)\in \mathscr{D}:=D_{0}^{1,2}(\Omega)\times D_{0}^{1,2}(\Omega), \end{cases}$$ where $s_1\neq s_2\in (0,2), \alpha>1,\beta>1, \lambda>0,\mu>0,\kappa>0, \alpha+\beta=2^*(s_2)$. Here, $2^*(s):=\frac{2(N-s)}{N-2}$ is the critical Hardy-Sobolev exponent. When $\Omega$ is a cone (especially $\Omega=\R_+^N$ or $\Omega=\R^N$), we study the existence of positive ground state solution.


title: Evolution of Kinetic and Magnetic Energy in Intra Cluster Media

abstract: Intra Cluster Media (ICMs) located at galaxy clusters is in the state of hot, tenuous, magnetized, and highly ionized X-ray emitting plasmas. This overall collisionless, viscous, and conductive magnetohydrodynamic (MHD) turbulence in ICM is simulated using hyper and physical magnetic diffusivity. The results show that fluctuating random plasma motion amplifies the magnetic field, which cascades toward the diffusivity scale passing through the viscous scale. The kinetic eddies in the subviscous scale are driven and constrained by the magnetic tension which finally gets balanced with the highly damping effect of the kinetic eddies. However, the saturated kinetic energy spectrum is deeper than that of the incompressible or compressible hydrodynamics fluid. To explain this unusual field profile we set up two simultaneous differential equations for the kinetic and magnetic energy spectrum using an Eddy Damped Quasi Normal Markovianized (EDQNM) approximation. The analytic solution tells us that the magnetic energy in addition to the viscous damping effect constrains the plasma motion leading to the power spectra: kinetic energy spectrum $E_V^k\sim k^{-3}$ and corresponding representative magnetic energy spectrum $E_M^k\sim k^{-1/2}$. Also the comparison of simulation results with different resolutions and magnetic diffusivities implies the role of small scale magnetic energy in dynamo.


title: The Archimedean Projection Property

abstract: Let $H$ be a hypersurface in $\mathbb R^n$ and let $\pi$ be an orthogonal projection in $\mathbb R^n$ restricted to $H$. We say that $H$ satisfies the $Archimedean$ $projection$ $property$ corresponding to $\pi$ if there exists a constant $C$ such that $Vol(\pi^{-1}(U)) = C \cdot Vol(U)$ for every measurable $U$ in the range of $\pi$. It is well-known that the $(n-1)$-dimensional sphere, as a hypersurface in $\mathbb R^n$, satisfies the Archimedean projection property corresponding to any codimension 2 orthogonal projection in $\mathbb R^n$, the range of any such projection being an $(n-2)$-dimensional ball. Here we construct new hypersurfaces that satisfy Archimedean projection properties. Our construction works for any projection codimension $k$, $2 \leq k \leq n - 1$, and it allows us to specify a wide variety of desired projection ranges $\Omega^{n-k} \subset \mathbb R^{n-k}$. Letting $\Omega^{n-k}$ be an $(n-k)$-dimensional ball for each $k$, it produces a new family of smooth, compact hypersurfaces in $\mathbb R^n$ satisfying codimension $k$ Archimedean projection properties that includes, in the special case $k = 2$, the $(n-1)$-dimensional spheres.


title: Unwinding of circular helicoidal molecules versus size

abstract: The thermodynamical stability of a set of circular double helical molecules is analyzed by path integral techniques. The minicircles differ only in \textit{i)} the radius and \textit{ii)} the number of base pairs ($N$) arranged along the molecule axis. Instead, the rise distance is kept constant. For any molecule size, the computational method simulates a broad ensemble of possible helicoidal configurations while the partition function is a sum over the path trajectories describing the base pair fluctuational states. The stablest helical repeat of every minicircle is determined by free energy minimization. We find that, for molecules with $N$ larger than $100$, the helical repeat grows linearly with the size and the twist number is constant. On the other hand, by reducing the size below $100$ base pairs, the double helices sharply unwind and the twist number drops to one for $N=\,20$. This is predicted as the minimum size for the existence of helicoidal molecules in the closed form. The helix unwinding appears as a strategy to release the bending stress associated to the circularization of the molecules.


title: Polarization in exoplanetary systems caused by transits, grazing transits, and starspots

abstract: We present results of numerical simulations of flux and linear polarization variations in transiting exoplanetary systems, caused by the host star disk symmetry breaking. We consider different configurations of planetary transits depending on orbital parameters. Starspot contribution to the polarized signal is also estimated. Applying the method to known systems and simulating observational conditions, a number of targets is selected where transit polarization effects could be detected. We investigate several principal benefits of the transit polarimetry, particularly, for determining orbital spatial orientation and distinguishing between grazing and near-grazing planets. Simulations show that polarization parameters are also sensitive to starspots, and they can be used to determine spot positions and sizes.