title: Isotropic matroids I: Multimatroids and neighborhoods

abstract: Several properties of the isotropic matroid of a looped simple graph are presented. Results include a characterization of the multimatroids that are associated with isotropic matroids and several ways in which the isotropic matroid of G incorporates information about graphs locally equivalent to G. Specific results of the latter type include a characterization of graphs that are locally equivalent to bipartite graphs, a direct proof that two forests are isomorphic if and only if their isotropic matroids are isomorphic, and a way to express local equivalence indirectly, using only edge pivots.


title: Spatial autocorrelation approaches to testing residuals from least squares regression

abstract: In statistics, the Durbin-Watson test is always employed to detect the presence of serial correlation of residuals from a least squares regression analysis. However, the Durbin-Watson statistic is only suitable for ordered time or spatial series. If the variables comprise cross-sectional data coming from spatial random sampling, the Durbin-Watson will be ineffectual because the value of Durbin-Watson's statistic depends on the sequences of data point arrangement. Based on the ideas from spatial autocorrelation, this paper presents two new statistics for testing serial correlation of residuals from least squares regression based on spatial samples. By analogy with the new form of Moran's index, an autocorrelation coefficient is defined with a standardized residual vector and a normalized spatial weight matrix. Then on the analogy of the Durbin-Watson statistic, a serial correlation index is constructed. As a case, the two statistics are applied to the spatial sample of 29 China's regions. These results show that the new spatial autocorrelation model can be used to test the serial correlation of residuals from regression analysis. In practice, the new statistics can make up for the deficiency of the Durbin-Watson test.


title: Expressing Forms as a sum of Pfaffians

abstract: Let A= (a_{ij}) be a symmetric non-negative integer 2k x 2k matrix. A is homogeneous if a_{ij} + a_{kl}=a_{il} + a_{kj} for any choice of the four indexes. Let A be a homogeneous matrix and let F be a general form in C[x_1, \dots x_n] with 2deg(F) = trace(A). We look for the least integer, s(A), so that F= pfaff(M_1) + \cdots + pfaff(M_{s(A)}), where the M_i's are 2k x 2k skew-symmetric matrices of forms with degree matrix A. We consider this problem for n= 4 and we prove that s(A) < k+1 for all A.


title: Expressing a General Form as a Sum of Determinants

abstract: Let A= (a_{ij}) be a non-negative integer k x k matrix. A is a homogeneous matrix if a_{ij} + a_{kl}=a_{il} + a_{kj} for any choice of the four indexes. We ask: If A is a homogeneous matrix and if F is a form in C[x_1, \dots x_n] with deg(F) = trace(A), what is the least integer, s(A), so that F = det M_1 + ... + det M_{s(A)}, where the M_i's are k x k matrices of forms with degree matrix A? We consider this problem for n>3 and we prove that s(A) is at most k^{n-3} and s(A) <k^{n-3} in infinitely many cases. However s(A) = k^{n-3} when the entries of A are large with respect to k.


title: Characterizing common cause closedness of quantum probability theories

abstract: We prove new results on common cause closedness of quantum probability spaces, where by a quantum probability space is meant the projection lattice of a non-commutative von Neumann algebra together with a countably additive probability measure on the lattice. Common cause closedness is the feature that for every correlation between a pair of commuting projections there exists in the lattice a third projection commuting with both of the correlated projections and which is a Reichenbachian common cause of the correlation. The main result we prove is that a quantum probability space is common cause closed if and only if it has at most one measure theoretic atom. This result improves earlier ones published in Z. GyenisZ and M. Redei Erkenntnis 79 (2014) 435-451. The result is discussed from the perspective of status of the Common Cause Principle. Open problems on common cause closedness of general probability spaces $(\mathcal{L},\phi)$ are formulated, where $\mathcal{L}$ is an orthomodular bounded lattice and $\phi$ is a probability measure on $\mathcal{L}$.


title: Oscillatory integrals related to Carleson's theorem: fractional monomials

abstract: Stein and Wainger proved the $L^p$ bounds of the polynomial Carleson operator for all integer-power polynomials without linear term. In the present paper, we partially generalise this result to all fractional monomials in dimension one. Moreover, the connections with Carleson's theorem and the Hilbert transform along vector fields or (variable) curves are also discussed in details.


title: Central Mass Profiles of the Nearby Cool-core Galaxy Clusters Hydra A and A478

abstract: We perform a weak-lensing study of the nearby cool-core galaxy clusters, Hydra A ($z=0.0538$) and A478 ($z=0.0881$), of which brightest cluster galaxies (BCGs) host powerful activities of active galactic nuclei (AGNs). For each cluster, the observed tangential shear profile is well described either by a single Navarro--Frenk--White model or a two-component model including the BCG as an unresolved point mass. For A478, we determine the BCG and its host-halo masses from a joint fit to weak-lensing and stellar photometry measurements. We find that the choice of initial mass functions (IMFs) can introduce a factor of two uncertainty in the BCG mass, whereas the BCG host halo mass is well constrained by data. We perform a joint analysis of weak-lensing and stellar kinematics data available for the Hydra A cluster, which allows us to constrain the central mass profile without assuming specific IMFs. We find that the central mass profile ($r<300$kpc) determined from the joint analysis is in excellent agreement with those from independent measurements, including dynamical masses estimated from the cold gas disk component, X-ray hydrostatic total mass estimates, and the central stellar mass estimated based on the Salpeter IMF. The observed dark-matter fraction around the BCG for Hydra A is found to be smaller than those predicted by adiabatic contraction models, suggesting the importance of other physical processes, such as the the AGN feedback and/or dissipationless mergers.


title: Exponential Stabilization of Nonholonomic Systems by means of Oscillating Controls

abstract: This paper is devoted to the stabilization problem for nonlinear driftless control systems by means of a time-varying feedback control. It is assumed that the vector fields of the system together with their first order Lie brackets span the whole tangent space at the equilibrium. A family of trigonometric open-loop controls is constructed to approximate the gradient flow associated with a Lyapunov function. These controls are applied for the derivation of a time-varying feedback law under the sampling strategy. By using Lyapunov's direct method, we prove that the controller proposed ensures exponential stability of the equilibrium. As an example, this control design procedure is applied to stabilize the Brockett integrator.


title: Possible Excitonic Phase of Graphite in the Quantum Limit State

abstract: The in-plane resistivity, Hall resistivity and magnetization of graphite were investigated in pulsed magnetic fields applied along the \textit{c}-axis. The Hall resistivity approaches zero at around 53 T where the in-plane and out-of-plane resistivities steeply decrease. The differential magnetization also shows an anomaly at around this field with a similar amplitude compared to that of de Haas-van Alphen oscillations at lower fields. This transition field appears insensitive to disorder, but reduces with doping holes. These results suggest the realization of the quantum limit states above 53 T. As a plausible explanation for the observed gapped out-of-plane conduction above 53 T, the emergence of the excitonic BCS-like state in graphite is proposed.


title: Fluctuations of the Self-Normalized Sum in the Curie-Weiss Model of SOC

abstract: We extend the main theorem of arXiv:1301.6911 about the fluctuations in the Curie-Weiss model of SOC. We present a short proof using the Hubbard-Stratonovich transformation with the self-normalized sum of the random variables.


title: Restoration of supersymmetry in two-dimensional SYM with sixteen supercharges on the lattice

abstract: We perform lattice simulations of two-dimensional supersymmetric Yang-Mills theory with sixteen supercharges with a lattice action which has two exact supercharges (Sugino lattice action). According to the gauge/gravity duality, the theory at finite temperature is expected to be well described by the corresponding black 1-branes, at low temperature in the large N limit. We aim to confirm the duality conjecture by comparing the lattice results with the theoretical predictions obtained in the gravity side. In this article, at the beginning of this study, we examine the supersymmetric Ward-Takahashi identity to test whether the lattice action reproduces the correct continuum theory. Numerical results of the SUSY WTI strongly suggest us that any cut-off effects, which break supersymmetry, disappear in the continuum limit. In addition, we study the issue of degenerate vacua and find that the admissiblilty condition or any other constraints of the link fields which guarantee the unique vacuum are not always needed.


title: Tunable mechanical and thermal properties of ZnS/CdS core/shell nanowires

abstract: Using all atom molecular dynamics (MD) simulations, we have studied the mechanical properties of ZnS/CdS core/shell nanowires. Our results show that the coating of a few atomic layer CdS shell on the ZnS nanowire leads to a significant change in the stiffness of the core/shell nanowires compared to the stiffness of pure ZnS nanowires. The binding energy between the core and shell region decreases due to the lattice mismatch at the core-shell interface. This reduction in binding energy plays an important role in determining the stiffness of a core/shell nanowire. We have also investigated the effects of the shell on the thermal conductivity and melting behavior of the nanowires.


title: On totally decomposable algebras with involution in characteristic two

abstract: A necessary and sufficient condition for a central simple algebra with involution over a field of characteristic two to be decomposable as a tensor product of quaternion algebras with involution, in terms of its Frobenius subalgebras, is given. It is also proved that a bilinear Pfister form, recently introduced by A. Dolphin, can classify totally decomposable central simple algebras of orthogonal type.


title: Robust criticality of Ising model on rewired directed networks

abstract: We show that preferential rewiring, which is supposed to mimick the behaviour of financial agents, changes a directed-network Ising ferromagnet with a single critical point into a model with robust critical behaviour. For the non-rewired random graph version, due to a constant number of links out-going from each site, we write a simple mean-field-like equation describing the behaviour of magnetization; we argue that it is exact and support the claim with extensive Monte Carlo simulations. For the rewired version, this equation is obeyed only at low temperatures. At higher temperatures, rewiring leads to strong heterogeneities, which apparently invalidates mean-field arguments and induces large fluctuations and divergent susceptibility. Such behaviour is traced back to the formation of a relatively small core of agents which influence the entire system.


title: Entropy degeneration of convex projective surfaces

abstract: We show that the volume entropy of the Hilbert metric on a closed convex projective surface tends to zero as the corresponding Pick differential tends to infinity. The proof is based on the theorem, due to Benoist and Hulin, that the Hilbert metric and Blaschke metric are comparable.


title: Production of heavy neutron-rich nuclei in transfer reactions within the dinuclear system model

abstract: The dynamics of nucleon transfer processes in heavy-ion collisions is investigated within the dinuclear system model. The production cross sections of nuclei in the reactions $^{136}$Xe+$^{208}$Pb and $^{238}$U+$^{248}$Cm are calculated and the calculations are in good agreement with the experimental data. The transfer cross sections for the $^{58}$Ni+$^{208}$Pb reaction are calculated and compared with the experimental data. We predict the production cross sections of neutron-rich nuclei $^{165-168}$Eu, $^{169-173}$Tb, $^{173-178}$Ho, and $^{181-185}$Yb based on the reaction $^{176}$Yb+$^{238}$U. It can be seen that the production cross sections of the neutron-rich nuclei $^{165}$Eu, $^{169}$Tb, $^{173}$Ho, and $^{181}$Yb are 2.84 $\mu$b, 6.90 $\mu$b, 46.24 $\mu$b, and 53.61 $\mu$b, respectively, which could be synthesized in experiment.


title: Making Availability as a Service in the Clouds

abstract: Cloud computing has achieved great success in modern IT industry as an excellent computing paradigm due to its flexible management and elastic resource sharing. To date, cloud computing takes an irrepalceable position in our socioeconomic system and influences almost every aspect of our daily life. However, it is still in its infancy, many problems still exist.Besides the hotly-debated security problem, availability is also an urgent issue.With the limited power of availability mechanisms provided in present cloud platform, we can hardly get detailed availability information of current applications such as the root causes of availability problem,mean time to failure, etc. Thus a new mechanism based on deep avaliability analysis is neccessary and benificial.Following the prevalent terminology 'XaaS',this paper proposes a new win-win concept for cloud users and providers in term of 'Availability as a Service' (abbreviated as 'AaaS').The aim of 'AaaS' is to provide comprehensive and aimspecific runtime avaliabilty analysis services for cloud users by integrating plent of data-driven and modeldriven approaches. To illustrate this concept, we realize a prototype named 'EagleEye' with all features of 'AaaS'. By subscribing corresponding services in 'EagleEye', cloud users could get specific availability information of their applications deployed in cloud platform. We envision this new kind of service will be merged into the cloud management mechanism in the near future.


title: A Finite Element Framework for Some Mimetic Finite Difference Discretizations

abstract: In this work we derive equivalence relations between mimetic finite difference schemes on simplicial grids and modified N\'ed\'elec-Raviart-Thomas finite element methods for model problems in $\mathbf{H}(\operatorname{\mathbf{curl}})$ and $H(\operatorname{div})$. This provides a simple and transparent way to analyze such mimetic finite difference discretizations using the well-known results from finite element theory. The finite element framework that we develop is also crucial for the design of efficient multigrid methods for mimetic finite difference discretizations, since it allows us to use canonical inter-grid transfer operators arising from the finite element framework. We provide special Local Fourier Analysis and numerical results to demonstrate the efficiency of such multigrid methods.


title: Bridging Social Media via Distant Supervision

abstract: Microblog classification has received a lot of attention in recent years. Different classification tasks have been investigated, most of them focusing on classifying microblogs into a small number of classes (five or less) using a training set of manually annotated tweets. Unfortunately, labelling data is tedious and expensive, and finding tweets that cover all the classes of interest is not always straightforward, especially when some of the classes do not frequently arise in practice. In this paper we study an approach to tweet classification based on distant supervision, whereby we automatically transfer labels from one social medium to another for a single-label multi-class classification task. In particular, we apply YouTube video classes to tweets linking to these videos. This provides for free a virtually unlimited number of labelled instances that can be used as training data. The classification experiments we have run show that training a tweet classifier via these automatically labelled data achieves substantially better performance than training the same classifier with a limited amount of manually labelled data; this is advantageous, given that the automatically labelled data come at no cost. Further investigation of our approach shows its robustness when applied with different numbers of classes and across different languages.


title: Propagation of moments and uniqueness of weak solutions to the Vlasov-Poisson-Fokker-Planck system

abstract: In this paper, we prove the uniqueness of weak solutions to the Vlasov-Poisson-Fokker-Planck system in $C([0,T]; L^p)$, by assuming the solution has a local bounded density which tends to infinite with a "reasonable" rate as $t\to 0$. And particularly as a corollary, we get the uniqueness of weak solutions with initial data $f_0$ satisfying $f_0|v|^2\in L^1$, which solves the uniqueness of weak solutions with finite energy. In addition, we prove that the moments with respect to the velocity propagate for any order higher than 2.


title: An Improved Pseudo-Polynomial Upper Bound for the Value Problem and Optimal Strategy Synthesis in Mean Payoff Games

abstract: In this work we offer an $O(|V|^2 |E|\, W)$ pseudo-polynomial time deterministic algorithm for solving the Value Problem and Optimal Strategy Synthesis in Mean Payoff Games. This improves by a factor $\log(|V|\, W)$ the best previously known pseudo-polynomial time upper bound due to Brim,~\etal The improvement hinges on a suitable characterization of values, and a description of optimal positional strategies, in terms of reweighted Energy Games and Small Energy-Progress Measures.


title: Rates of convergence of rho-estimators for sets of densities satisfying shape constraints

abstract: The purpose of this paper is to pursue our study of rho-estimators built from i.i.d. observations that we defined in Baraud et al. (2014). For a \rho-estimator based on some model S (which means that the estimator belongs to S) and a true distribution of the observations that also belongs to S, the risk (with squared Hellinger loss) is bounded by a quantity which can be viewed as a dimension function of the model and is often related to the "metric dimension" of this model, as defined in Birg\'e (2006). This is a minimax point of view and it is well-known that it is pessimistic. Typically, the bound is accurate for most points in the model but may be very pessimistic when the true distribution belongs to some specific part of it. This is the situation that we want to investigate here. For some models, like the set of decreasing densities on [0,1], there exist specific points in the model that we shall call "extremal" and for which the risk is substantially smaller than the typical risk. Moreover, the risk at a non-extremal point of the model can be bounded by the sum of the risk bound at a well-chosen extremal point plus the square of its distance to this point. This implies that if the true density is close enough to an extremal point, the risk at this point may be smaller than the minimax risk on the model and this actually remains true even if the true density does not belong to the model. The result is based on some refined bounds on the suprema of empirical processes that are established in Baraud (2016).


title: Totally-Reflective Genera of Integral Lattices

abstract: In this paper we give a complete classification of totally-reflective, primitive genera in dimension 3 and 4. Our method breaks up into two parts. The first part consists of classifying the square free, totally-reflective, primitive genera by calculating strong bounds on the prime factors of the determinant of genera of positive definite quadratic forms (lattices) with this property. We achieve these bounds by combining the Minkowski-Siegel mass formula with the combinatorial classification of reflective lattices accomplished by Scharlau \& Blaschke. In a second part, we use a lattice transformation that goes back to Watson, to generate all totally-reflective, primitive genera when starting with the square-free case.


title: Relative Stability and Local Curvature Analysis in Carbon Nanotori

abstract: We introduce a concise formalism to characterize nanometer-sized tori based on carbon nanotubes and to determine their stability by combining {\em ab initio} density functional calculations with a continuum elasticity theory approach that requires only shape information. We find that the high strain energy in nanotori containing only hexagonal rings is significantly reduced in nanotori containing also other polygons. Our approach allows to determine local curvature and link it to local strain energy, which is correlated with local stability and chemical reactivity.


title: Search for a Heavy Neutral Particle Decaying to $e\mu$, $e\tau$, or $\mu\tau$ in $pp$ Collisions at $\sqrt{s}=8$ TeV with the ATLAS Detector

abstract: This Letter presents a search for a heavy neutral particle decaying into an opposite-sign different-flavor dilepton pair, $e^\pm \mu^\mp$, $e^\pm \tau^\mp$, or $\mu^\pm \tau^\mp$ using 20.3 fb$^{-1}$ of $pp$ collision data at $\sqrt{s} = 8$ TeV collected by the ATLAS detector at the LHC. The numbers of observed candidate events are compatible with the Standard Model expectations. Limits are set on the cross section of new phenomena in two scenarios: the production of $\tilde{\nu}_{\tau}$ in $R$-parity-violating supersymmetric models and the production of a lepton-flavor-violating $Z'$ vector boson.


title: $X(3915)$ as a $D_s \bar D_s$ bound state

abstract: We suggest that the observed properties of the charmoniumlike resonance $X(3915)$ can possibly be explained if it is an $S-$wave molecular bound state of $D_s \bar D_s$ meson pair with binding energy about 18\,MeV. In particular, the decays of $X(3915)$ to pairs of non-strange $D$ mesons are suppressed and proceed at a rate comparable to that of the decay $X(3915) \to \omega J/\psi$, whose branching fraction can be as large as about 0.3. Other major types of decay of $X(3915)$ with a comparable (or slightly smaller) rate are the transition $X(3915) \to \eta \eta_c$ and the decays into light hadrons due to the annihilation of the $c \bar c$ quark pair. The existence of a bound state should lead to an enhancement in the spectrum of the invariant mass for the $D_s \bar D_s$ near threshold in $B$ decays, e.g. in $B \to K D_s \bar D_s$ which enhancement can be tested experimentally.


title: Sequence determines degree of knottedness in a coarse-grained protein model

abstract: Knots are abundant in globular homopolymers but rare in globular proteins. To shed new light on this long-standing conundrum, we study the influence of sequence on the formation of knots in proteins under native conditions within the framework of the hydrophobic-polar (HP) lattice protein model. By employing large scale Wang-Landau simulations combined with suitable Monte Carlo trial moves we show that, even though knots are still abundant on average, sequence introduces large variability in the degree of self-entanglements. Moreover, we are able to design sequences which are either almost always or almost never knotted. Our findings serve as proof of concept that the introduction of just one additional degree of freedom per monomer (in our case sequence) facilitates evolution towards a protein universe in which knots are rare.


title: Versatile approach to access the low temperature thermodynamics of lattice polymers and proteins

abstract: We show that Wang-Landau sampling, combined with suitable Monte Carlo trial moves, provides a powerful method for both the ground state search and the determination of the density of states for the hydrophobic-polar (HP) protein model and the interacting self-avoiding walk (ISAW) model for homopolymers. We obtained accurate estimates of thermodynamic quantities for HP sequences with $>100$ monomers and for ISAWs up to $>500$ monomers. Our procedure possesses an intrinsic simplicity and overcomes the limitations inherent in more tailored approaches making it interesting for a broad range of protein and polymer models.


title: Relativistic modelling of stable anisotropic super-dense star

abstract: In the present article we have obtained new set of exact solutions of Einstein field equations for anisotropic fluid spheres by using the Herrera et al.[1] algorithm. The anisotropic fluid solution so obtained join continuously to Schwarzschild exterior solution across the pressure free boundary. It is observed that most of the new anisotropic solutions are well behaved and utilized to construct the super-dense star models such as neutron star and pulsars.


title: Precise mass and radius measurements for the components of the bright solar-type eclipsing binary star V1094 Tau

abstract: V1094 Tau is bright eclipsing binary star with an orbital period close to 9 days containing two stars similar to the Sun. Our aim is to test models of Sun-like stars using precise and accurate mass and radius measurements for both stars in V1094 Tau. We present new spectroscopy of V1094 Tau which we use to estimate the effective temperatures of both stars and to refine their spectroscopic orbits. We also present new, high-quality photometry covering both eclipses of V1094 Tau in the Stroemgren uvby system and in the Johnson V-band. The masses, radii and effective temperatures of the stars in V1094 Tau are found to be M$_A$ = 1.0965 $\pm$ 0.0040 M$_{\odot}$, R$_A$ = 1.4109 $\pm$ 0.0058 R$_{\odot}$, T$_{\rm eff,A}$ = 5850 $\pm$ 100 K, and M$_B$ = 1.0120 $\pm$ 0.0028 M$_{\odot}$, R$_B$ = 1.1063 $\pm$ 0.0066 R$_{\odot}$, T$_{\rm eff,B}$ = 5700 $\pm$ 100 K. An analysis of the times of mid-eclipse and the radial velocity data reveals apsidal motion with a period of 14500 $\pm$ 3700 years. The observed masses, radii and effective temperatures are consistent with stellar models for an age $\approx$ 6 Gyr if the stars are assumed to have a metallicity similar to the Sun. This estimate is in reasonable agreement with our estimate of the metallicity derived using Stroemgren photometry and treating the binary as a single star ([Fe/H] $= -0.09 \pm 0.11$). The rotation velocities of the stars suggest that V1094 Tau is close to the limit at which tidal interactions between the stars force them to rotate pseudo-synchronously with the orbital motion.


title: Hint at an axion-like particle from the redshift dependence of blazar spectra

abstract: We consider the largest observed sample including all intermediate-frequency peaked (IBL) and high-frequency peaked (HBL) flaring blazars above 100 GeV up to redshift $z = 0.6$. We show that the best-fit regression line of the emitted spectral indices $\Gamma_{\rm em} (z)$ is a concave parabola decreasing as $z$ increases, thereby implying a statistical correlation between the $\{\Gamma_{\rm em} (z) \}$ distribution and $z$. This result contradicts our expectation that such a distribution should be $z$-independent. We argue that the above correlation does not arise from any selection bias. We show that our expectation naturally emerges provided that axion-like particles (ALPs) are put into the game. Moreover, ALPs can also explain why flat spectrum radio quasars emit up to 400 GeV, in sharp contradiction with conventional physics. So, the combination of the two very different but consistent results -- taken at face value -- leads to a hint at an ALP with mass $m = {\cal O} (10^{-10} \, {\rm eV})$ and two-photon coupling in the range $2.94 \times 10^{- 12} \, {\rm GeV}^{- 1} < g_{a \gamma \gamma} < 0.66 \times 10^{- 10} \, {\rm GeV}^{- 1}$. As a bonus, the Universe would become considerably more transparent above energies $E \gtrsim 1 \, {\rm TeV}$ than dictated by conventional physics. Our prediction can be checked not only by the new generation of observatories like CTA, HAWC, GAMMA-400, LHAASO, TAIGA-HiSCORE and HERD, but also thanks to the planned laboratory experiments ALPS II (upgraded), STAX, IAXO and with other techniques now being developed by Avignone and collaborators.


title: A mean value inequality for the generalized self-expander type submanifolds and its application

abstract: In this paper we get a version of mean value inequality for generalized self-expander type submanifolds in Euclidean space. As the application, we prove that if mean curvature flow $M(t)$ on the self-expander in Euclidean space subconverges to an $n$-rectifiable varifold $T$ in weak sense for $t$ goes to the singular time, then $T$ must be the cone.


title: Stochastic Stability Analysis of Discrete Time System Using Lyapunov Measure

abstract: In this paper, we study the stability problem of a stochastic, nonlinear, discrete-time system. We introduce a linear transfer operator-based Lyapunov measure as a new tool for stability verification of stochastic systems. Weaker set-theoretic notion of almost everywhere stochastic stability is introduced and verified, using Lyapunov measure-based stochastic stability theorems. Furthermore, connection between Lyapunov functions, a popular tool for stochastic stability verification, and Lyapunov measures is established. Using the duality property between the linear transfer Perron-Frobenius and Koopman operators, we show the Lyapunov measure and Lyapunov function used for the verification of stochastic stability are dual to each other. Set-oriented numerical methods are proposed for the finite dimensional approximation of the Perron-Frobenius operator; hence, Lyapunov measure is proposed. Stability results in finite dimensional approximation space are also presented. Finite dimensional approximation is shown to introduce further weaker notion of stability referred to as coarse stochastic stability. The results in this paper extend our earlier work on the use of Lyapunov measures for almost everywhere stability verification of deterministic dynamical systems ("Lyapunov Measure for Almost Everywhere Stability", {\it IEEE Trans. on Automatic Control}, Vol. 53, No. 1, Feb. 2008).


title: Time-Dependent Meson Melting in External Magnetic Field

abstract: The dynamics of a probe D7-brane in an asymptotically AdS-Vaidya background has been investigated in the presence of an external magnetic field. Holographically, this is dual to the dynamical meson melting in the N = 2 super Yang-Milles theory. If the final temperature of the system is large enough, the probe D7-brane will dynamically cross the horizon (black hole embedding). By turning on the external magnetic field and raising it sufficiently, the final embedding of the corresponding D7-brane changes to Minkowski embedding. In the field theory side, this means that the mesons which melt due to the raise in the temperature, will form bound states again by applying an external magnetic field. We will also show that the evolution of the system to its final equilibrium state is postponed due to the presence of the magnetic field.


title: Unconventional Ballooning Structures for Toroidal Drift Waves

abstract: With strong gradients in the pedestal of high confinement mode (H-mode) fusion plasmas, gyrokinetic simulations are carried out for the trapped electron and ion temperature gradient modes. A broad class of unconventional mode structures is found to localize at arbitrary poloidal positions or with multiple peaks. It is found that these unconventional ballooning structures are associated with different eigen states for the most unstable mode. At weak gradient (low confinement mode or L-mode), the most unstable mode is usually in the ground eigen state, which corresponds to a conventional ballooning mode structure peaking in the outboard mid-plane of tokamaks. However, at strong gradient (H-mode), the most unstable mode is usually not the ground eigen state and the ballooning mode structure becomes unconventional. This result implies that the pedestal of H-mode could have better confinement than L-mode.


title: Schwinger-Keldysh canonical formalism for electronic Raman scattering

abstract: Inelastic low-energy Raman and high-energy X-ray scatterings have made great progress in instrumentation to investigate the strong electronic correlations in matter. However, theoretical study of the relevant scattering spectrum is still a challenge. In this article, we present a Schwinger-Keldysh canonical perturbation formalism for the electronic Raman scattering, where all the resonant, non-resonant and mixed responses are considered uniformly. We show how to use this formalism to evaluate the cross section of the electronic Raman scattering off an one-band superconductor. All the two-photon scattering processes from electrons, the non-resonant charge density response, the elastic Rayleigh scattering, the fluorescence, the intrinsic energy-shift Raman scattering and the mixed response, are included. In the mean-field superconducting state, Cooper pairs contribute only to the non-resonant response. All the other responses are dominated by the single-particle excitations and are strongly suppressed due to the opening of the superconducting gap. Our formalism for the electronic Raman scattering can be easily extended to study the high-energy resonant inelastic X-ray scattering.


title: Double-nucleus elliptical MCG-01-12-005 in an X-ray emitting cluster of galaxies

abstract: The scenario of galaxy formation is believed to follow a structure that builds up from the bottom, with large galaxies being formed by several merging episodes of smaller ones. In this scenario a number of galaxies can be expected to be seen in the merging phase, with their external regions already mixed, while their nuclei, with stronger self-gravitation, are still recognizable as such. During a photometric monitoring of AGNs in the field of a long-exposure INTEGRAL pointing, we serendipitously found an elliptical galaxy in the center of the X-ray cluster (EXO 0422-086) with two nuclei. We performed surface photometry on our images and those of the SDSS archive and obtained slit spectra of both nuclei. Aperture photometry of the two stellar-like nuclei showed very similar colors in the SDSS image and in our Johnson BVRI images, which is typical of an elliptical galaxy nucleus. The spectra of the nuclei showed the typical absorption lines of an elliptical galaxy without appreciable emission lines. The redshifts derived from each nucleus were equal and fully consistent with the literature value (0.0397). We can therefore exclude the possibility that one of the nuclei is a foreground star or a background AGN and consider this elliptical galaxy as a bona fide example of a galaxy merger.


title: Extended Conway-Maxwell-Poisson distribution and its properties and applications

abstract: A new three parameter natural extension of the Conway-Maxwell-Poisson (COM-Poisson) distribution is proposed. This distribution includes the recently proposed COM-Poisson type negative binomial (COM-NB) distribution [Chakraborty, S. and Ong, S. H. (2014): A COM-type Generalization of the Negative Binomial Distribution, Accepted in Communications in Statistics-Theory and Methods] and the generalized COM-Poisson (GCOMP) distribution [Imoto, T. :(2014) A generalized Conway-Maxwell-Poisson distribution which includes the negative binomial distribution, Applied Mathematics and Computation, 247, 824-834]. The proposed distribution is derived from a queuing system with state dependent arrival and service rates and also from an exponential combination of negative binomial and COM-Poisson distribution. Some distributional, reliability and stochastic ordering properties are investigated. Computational asymptotic approximations, different characterizations, parameter estimation and data fitting example also discussed.


title: Pattern Recognition of Bearing Faults using Smoother Statistical Features

abstract: A pattern recognition (PR) based diagnostic scheme is presented to identify bearing faults, using time domain features. Vibration data is acquired from faulty bearings using a test rig. The features are extracted from the data, and processed prior to utilize in the PR process. The processing involves smoothing of feature distributions. This reduces the undesired impact of vibration randomness on the PR process, and thus enhances the diagnostic accuracy of the model.


title: Associative production of $\Upsilon$ and open charm at LHC

abstract: The yield of $\Upsilon$ associated with open charm has been estimated with different approaches. The crucial differences between SPS and DPS predictions are discussed.


title: Enhancement of the CP-odd effect in the nuclear electric dipole moment of $^6$Li

abstract: We calculate for the first time the electric dipole moment (EDM) of the $^6$Li nucleus within the alpha + p + n three-body cluster model using the Gaussian expansion method, assuming the one meson exchange P, CP-odd nuclear forces. It is found that the EDM of the $^6$Li is 2 times more sensitive on the isovector pion exchange P, CP-odd nuclear force than the deuteron EDM, due to the CP-odd interaction between the nucleons and the alpha cluster. The $^9$Be EDM is also calculated in the same framework as an alpha + alpha + n three-body system. We also test the ab initio calculation of the EDM of the deuteron, $^3$H and $^3$He nuclei using the realistic Argonne $v18$ nuclear force. In the ab initio calculations, good agreements with previous studies are obtained. We finally discuss the prospects for the new physics beyond the standard model.


title: Equipped graded graphs, projective limits of simplices, and their boundaries

abstract: In this paper, we develop a theory of equipped graded graphs (or Bratteli diagrams) and an alternative theory of projective limits of finite-dimensional simplices. An equipment is an additional structure on the graph, namely, a system of "cotransition" probabilities on the set of its paths. The main problem is to describe all probability measures on the path space of a graph with given cotransition probabilities; it goes back to the problem, posed by E.~B.~Dynkin in the 1960s, of describing exit and entrance boundaries for Markov chains. The most important example is the problem of describing all central measures, to which one can reduce the problems of describing states on AF-algebras or characters on locally finite groups. We suggest an unification of the whole theory, an interpretation of the notions of Martin, Choquet, and Dynkin boundaries in terms of equipped graded graphs and in terms of the theory of projective limits of simplices. In the last section, we study the new notion of "standardness" of projective limits of simplices and of equipped Bratteli diagrams, as well as the notion of "lacunarization."


title: Security in the Multi-Dimensional Fibonacci Protocol

abstract: Security against simple eavesdropping attacks is demonstrated for a recently proposed quantum key distribution protocol which uses the Fibonacci recursion relation to enable high-capacity key generation with entangled photon pairs. No transmitted pairs need to be discarded in reconciliation; the only pairs not used for key generation are those used for security-checking. Although the proposed approach does not allow eavesdropper-induced errors to be detected on single trials, it can nevertheless reveal the eavesdropper's action on the quantum channel by detecting changes in the distribution of outcome probabilities over multiple trials, and can do so as well as the BB84 protocol. The mutual information shared by the participants is calculated and used to show that a secret key can always be distilled.


title: Visual, user-interactive generation of bond networks in 3D particle configurations

abstract: We present a new program able to perform visual structural analysis on 3D particle systems called PASYVAT (PArticle SYstem Visual Analysis Tool). More specifically, it can select multiple interparticle distance ranges from a radial distribution function (RDF) plot and display them in 3D as bonds between the particles falling within the selected distance range, thus generating a network of bonds. This software can be used with any data set representing a system of points or other objects having a welldefined center of mass or geometric center in 3D space. In this article we describe the program and its internal structure, with emphasis on its applicability in the study of certain particle configurations, obtained from classical molecular dynamics simulation in condensed matter physics.


title: Stability of toroidal magnetic fields in stellar interiors

abstract: We present 3D MHD simulations of purely toroidal and mixed poloidal-toroidal magnetic field configurations to study the behavior of the Tayler instability. For the first time the simultaneous action of rotation and magnetic diffusion are taken into account and the effects of a poloidal field on the dynamic evolution of unstable toroidal magnetic fields is included. In the absence of diffusion, fast rotation (rotation rate compared to Alfv\'en frequency) is able to suppress the instability when the rotation and magnetic axes are aligned and when the radial field strength gradient p < 1.5. When diffusion is included, this system turns unstable for diffusion dominated and marginally diffusive dominated regions. If the magnetic and rotation axes are perpendicular to each other the stabilizing effect induced by the Coriolis force is scale dependent and decreases with increasing wavenumber. In toroidal fields with radial field gradients bigger than p > 1.5, rapid rotation does not suppress the instability but instead introduces a damping factor to the growth rate in agreement with the analytic predictions. For the mixed poloidal-toroidal fields we find an unstable axisymmetric mode, not predicted analytically, right at the stability threshold for the non-axisymmetric modes; it has been argued that an axisymmetric mode is necessary for the closure of the Tayler-Spruit dynamo loop.


title: Giant comets and mass extinctions of life

abstract: I find evidence for clustering in age of well-dated impact craters over the last 500 Myr. At least nine impact episodes are identified, with durations whose upper limits are set by the dating accuracy of the craters. Their amplitudes and frequency are inconsistent with an origin in asteroid breakups or Oort cloud disturbances, but are consistent with the arrival and disintegration in near-Earth orbits of rare, giant comets, mainly in transit from the Centaur population into the Jupiter family and Encke regions. About 1 in 10 Centaurs in Chiron-like orbits enter Earth-crossing epochs, usually repeatedly, each such epoch being generally of a few thousand years duration. On time-scales of geological interest, debris from their breakup may increase the mass of the near-Earth interplanetary environment by two or three orders of magnitude, yielding repeated episodes of bombardment and stratospheric dusting. I find a strong correlation between these bombardment episodes and major biostratigraphic and geological boundaries, and propose that episodes of extinction are most effectively driven by prolonged encounters with meteoroid streams during bombardment episodes. Possible mechanisms are discussed.


title: Svestka's Research: Then and Now

abstract: Zdenek Svestka's research work influenced many fields of solar physics, especially in the area of flare research. In this article I take five of the areas that particularly interested him and assess them in a "then and now" style. His insights in each case were quite sound, although of course in the modern era we have learned things that he could not readily have envisioned. His own views about his research life have been published recently in this journal, to which he contributed so much, and his memoir contains much additional scientific and personal information (Svestka, 2010).


title: Reentrant Superfluidity and Pair Density Wave in Single Component Dipolar Fermi Gases

abstract: We study the superfluidity of single component dipolar Fermi gases in three dimensions within a pairing fluctuation theory. The transition temperature $T_{c}$ for the dominant $p_z$ wave superfluidity exhibits a remarkable re-entrant behavior as a function of the pairing strength induced by the dipole-dipole interaction (DDI), which leads to an anisotropic pair dispersion. The anisotropy and the long range nature of the DDI cause $T_c$ to vanish for a narrow range of intermediate interaction strengths, where a pair density wave state emerges as the ground state. The superfluid density and thermodynamics below $T_{c}$, along with the density profiles in a harmonic trap, are investigated as well, throughout the BCS-BEC crossover. Implications for experiments are discussed.


title: Generalised Stability of ultraproducts of finite residue rings

abstract: We study ultraproducts of finite residue rings $\underset{n\in\mathbb{N}}{\prod}\mathbb{Z}/n\mathbb{Z} \diagup \mathcal{U} $ where $\mathcal{U}$ is a non-principal ultrafilter. We find sufficient conditions of the ultrafilter $\mathcal{U}$ to determine if the resulting ultraproduct $\underset{n\in\mathbb{N}}{\prod}\mathbb{Z}/n\mathbb{Z} \diagup \mathcal{U}$ has simple, NIP, NTP2 but not simple nor NIP, or TP2 theory, noting that all these four cases occur.


title: Estimation of the time of change in panel data

abstract: We consider the problem of estimating the common time of a change in the mean parameters of panel data when dependence is allowed between the panels in the form of a common factor. A CUSUM type estimator is proposed, and we establish first and second order asymptotics that can be used to derive consistent confidence intervals for the time of change. Our results improve upon existing theory in two primary directions. Firstly, the conditions we impose on the model errors only pertain to the order of their long run moments, and hence our results hold for nearly all stationary time series models of interest, including nonlinear time series like the ARCH and GARCH processes. Secondly, we study how the asymptotic distribution and norming sequences of the estimator depend on the magnitude of the changes in each panel and the common factor loadings. The performance of our results in finite samples is demonstrated with a Monte Carlo simulation study, and we consider applications to two real data sets: the exchange rates of 23 currencies with respect to the US dollar, and the GDP per capita in 113 countries.


title: Exotic attractors of the non-equilibrium Rabi-Hubbard model

abstract: We explore the phase diagram of the dissipative Rabi-Hubbard model, as could be realized by a Raman-pumping scheme applied to a coupled cavity array. There exist various exotic attractors, including ferroelectric, antiferroelectric, and inccomensurate fixed points, as well as regions of persistent oscillations. Many of these features can be understood analytically by truncating to the two lowest lying states of the Rabi model on each site. We also show that these features survive beyond mean-field, using Matrix Product Operator simulations.


title: The evolution of the X-ray phase lags during the outbursts of the black hole candidate GX 339-4

abstract: Owing to the frequency and reproducibility of its outbursts, the black-hole candidate GX 339-4 has become the standard against which the outbursts of other black-hole candidate are matched up. Here we present the first systematic study of the evolution of the X-ray lags of the broad-band variability component (0.008-5 Hz) in GX 339-4 as a function of the position of the source in the hardness-intensity diagram. The hard photons always lag the soft ones, consistent with previous results. In the low-hard state the lags correlate with X-ray intensity, and as the source starts the transition to the intermediate/soft states, the lags first increase faster, and then appear to reach a maximum, although the exact evolution depends on the outburst and the energy band used to calculate the lags. The time of the maximum of the lags appears to coincide with a sudden drop of the Optical/NIR flux, the fractional RMS amplitude of the broadband component in the power spectrum, and the appearance of a thermal component in the X-ray spectra, strongly suggesting that the lags can be very useful to understand the physical changes that GX 339-4 undergoes during an outburst. We find strong evidence for a connection between the evolution of the cut-off energy of the hard component in the energy spectrum and the phase lags, suggesting that the average magnitude of the lags is correlated with the properties of the corona/jet rather than those of the disc. Finally, we show that the lags in GX 339-4 evolve in a similar manner to those of the black-hole candidate Cygnus X-1, suggesting similar phenomena could be observable in other black-hole systems.


title: On the monodromy of almost toric fibrations on the complex projective plane

abstract: The monodromy of almost toric Lagrangian fibrations on the complex projective plane is described.


title: Center vortices as composites of monopole fluxes

abstract: We study the relation between the flux of a center vortex obtained from the center vortex model and the flux formed between monopoles obtained from the Abelian gauge fixing method. Motivated by the Monte Carlo simulations which have shown that almost all monopoles are sitting on the top of vortices, we construct the fluxes of center vortices for $SU(2)$ and $SU(3)$ gauge groups using fractional fluxes of monopoles. Then, we compute the potentials in the fundamental representation induced by center vortices and fractional fluxes of monopoles. We show that by combining the fractional fluxes of monopoles one can produce the center vortex fluxes for $SU(3)$ gauge group in a "center vortex model". Comparing the potentials, we conclude that the fractional fluxes of monopoles attract each other.


title: Optimal risk allocation in a market with non-convex preferences

abstract: The aims of this study are twofold. First, we consider an optimal risk allocation problem with non-convex preferences. By establishing an infimal representation for distortion risk measures, we give some necessary and sufficient conditions for the existence of optimal and asymptotic optimal allocations. We will show that, similar to a market with convex preferences, in a non-convex framework with distortion risk measures the boundedness of the optimal risk allocation problem depends only on the preferences. Second, we consider the same optimal allocation problem by adding a further assumption that allocations are co-monotone. We characterize the co-monotone optimal risk allocations within which we prove the "marginal risk allocations" take only the values zero or one. Remarkably, we can separate the role of the market preferences and the total risk in our representation.


title: Exact Controllability of the Distributed System Governed by the Wave Equation with Memory

abstract: We will consider exact controllability of the distributed system governed by the wave equation with memory. It will be proved that this mechanical system can be driven to rest in finite time, the absolute value of the distributed control function being bounded. In this case the memory kernel is a linear combination of exponentials.


title: Enabling entanglement distillation via optomechanics

abstract: Quantum networking based on optical Gaussian states, although promising in terms of scalability, is hindered by the fact that their entanglement cannot be distilled via Gaussian operations. We show that optomechanics, integrable (on-chip) availability, and particularly the scope to measure the mechanical degree of freedom, can address this problem. Here, one of the optical modes of a two-mode squeezed vacuum is injected into a single-sided Fabry-P\'{e}rot cavity and non-linearly coupled to a mechanical oscillator. Afterward, the position of the oscillator is measured using pulsed optomechanics and homodyne detection. We show that this measurement can supply non-Gaussian entangled states frequently enough to enable scalable entanglement distillation. Moreover, it can conditionally increase the initial entanglement under an optimal radiation-pressure interaction strength, which corresponds to an effective unsharp measurement of the photon number inside the cavity. We show how the resulting entanglement enhancement can be verified by using a standard teleportation procedure.


title: Point Charges and Polygonal Linkages

abstract: We investigate the critical points of Coulomb potential of point charges placed at the vertices of a planar polygonal linkage. It is shown that, for a collection of positive charges on a pentagonal linkage, there is a unique critical point in the set of convex configurations which is the point of absolute minimum. This enables us to prove that two controlling charges are sufficient to navigate between any two convex configurations of a pentagonal linkage.


title: Frequency assortativity can induce chaos in oscillator networks

abstract: We investigate the effect of preferentially connecting oscillators with similar frequency to each other in networks of coupled phase oscillators (i.e., frequency assortativity). Using the network Kuramoto model as an example, we find that frequency assortativity can induce chaos in the macroscopic dynamics. By applying a mean-field approximation in combination with the dimension reduction method of Ott and Antonsen, we show that the dynamics can be described by a low dimensional system of equations. We use the reduced system to characterize the macroscopic chaos using Lyapunov exponents, bifurcation diagrams, and time-delay embeddings. Finally, we show that the emergence of chaos stems from the formation of multiple groups of synchronized oscillators, i.e., meta-oscillators.


title: Severi Varieties and Brill-Noether theory of curves on abelian surfaces

abstract: Severi varieties and Brill-Noether theory of curves on K3 surfaces are well understood. Yet, quite little is known for curves on abelian surfaces. Given a general abelian surface $S$ with polarization $L$ of type $(1,n)$, we prove nonemptiness and regularity of the Severi variety parametrizing $\delta$-nodal curves in the linear system $|L|$ for $0\leq \delta\leq n-1=p-2$ (here $p$ is the arithmetic genus of any curve in $|L|$). We also show that a general genus $g$ curve having as nodal model a hyperplane section of some $(1,n)$-polarized abelian surface admits only finitely many such models up to translation; moreover, any such model lies on finitely many $(1,n)$-polarized abelian surfaces. Under certain assumptions, a conjecture of Dedieu and Sernesi is proved concerning the possibility of deforming a genus $g$ curve in $S$ equigenerically to a nodal curve. The rest of the paper deals with the Brill-Noether theory of curves in $|L|$. It turns out that a general curve in $|L|$ is Brill-Noether general. However, as soon as the Brill-Noether number is negative and some other inequalities are satisfied, the locus $|L|^r_d$ of smooth curves in $|L|$ possessing a $g^r_d$ is nonempty and has a component of the expected dimension. As an application, we obtain the existence of a component of the Brill-Noether locus $\mathcal{M}^r_{p,d}$ having the expected codimension in the moduli space of curves $\mathcal{M}_p$. For $r=1$, the results are generalized to nodal curves.


title: Re-examining Archie's law: Conductance description by tortuosity and constriction

abstract: In this article we investigate the electrical conductance of an insulating porous medium (e.g., a sedimentary rock) filled with an electrolyte (e.g., brine), usually described using the Archie cementation exponent. We show how the electrical conductance depends on changes in the drift velocity and the length of the electric field lines, in addition to the porosity and the conductance of the electrolyte. We characterized the length of the electric field lines by a tortuosity and the changes in drift velocity by a constriction factor. Both the tortuosity and the constriction factor are descriptors of the pore microstructure. We define a conductance reduction factor to measure the local contributions of the pore microstructure to the global conductance. It is shown that the global conductance reduction factor is the product of the tortuosity squared divided by the constriction factor, thereby proving that the combined effect of tortuosity and constriction, in addition to the porosity and conductance of the electrolyte, fully describes the effective electrical conductance of a porous medium. We show that our tortuosity, constriction factor, and conductance reduction factor reproduce the electrical conductance for idealized porous media. They are also applied to Bentheimer sandstone, where we describe a microstructure-related correlation between porosity and conductivity using both the global conductance reduction factor and the distinct contributions from tortuosity and constriction. Overall, this work shows how the empirical Archie cementation exponent can be substituted by more descriptive, physical parameters, either by the global conductance reduction factor or by tortuosity and constriction.


title: A Consistency Relation for the Observed Galaxy Bispectrum and the Local non-Gaussianity from Relativistic Corrections

abstract: We obtain a consistency relation for the observed three-point correlator of galaxies. It includes relativistic effects and it is valid in the squeezed limit. Furthermore, the consistency relation is non-perturbative and can be used at arbitrarily small scales for the short modes. Our results are also useful to compute the non-linear relativistic corrections which induce a signal in the observations that might be misinterpreted as primordial non-Gaussianity with a local shape. We estimate the effective local non-Gaussianity parameter from the relativistic corrections. The exact value depends on the redshift and the magnification bias. At redshift of $ z = 1$, in the absence of magnification bias, we get $\,\, f^{\rm GR}_{\rm NL} \simeq - 3.7 $.


title: Ordered Median Hub Location Problems with Capacity Constraints

abstract: The Single Allocation Ordered Median Hub Location problem is a recent hub model introduced in Puerto et al. (2011) that provides a unifying analysis of a wide class of hub location mod- els. In this paper, we deal with the capacitated version of this problem, presenting two formulations as well as some preprocessing phases for fixing variables. In addition, a strengthening of one of these formulations is also studied through the use of some fami- lies of valid inequalities. A battery of test problems with data taken from the AP library are solved where it is shown that the running times have been significantly reduced with the improvements presented in the paper.


title: From physical principles to relativistic classical Hamiltonian and Lagrangian particle mechanics

abstract: We show that classical particle mechanics (Hamiltonian and Lagrangian consistent with relativistic electromagnetism) can be derived from three fundamental assumptions: infinite reducibility, deterministic and reversible evolution, and kinematic equivalence. The core idea is that deterministic and reversible systems preserve the cardinality of a set of states, which puts considerable constraints on the equations of motion. This perspective links different concepts from different branches of math and physics (e.g. cardinality of a set, cotangent bundle for phase space, Hamiltonian flow, locally Minkowskian space-time manifold), providing new insights. The derivation strives to use definitions and mathematical concepts compatible with future extensions to field theories and quantum mechanics.


title: A criterion for the existence of zero modes for the Pauli operator with fastly decaying fields

abstract: We consider the Pauli operator in $\mathbb R^3$ for magnetic fields in $L^{3/2}$ that decay at infinity as $|x|^{-2-\beta}$ with $\beta > 0$. In this case we are able to prove that the existence of a zero mode for this operator is equivalent to a quantity $\delta(\mathbf B)$, defined below, being equal to zero. Complementing a result from [Balinsky, Evans, Lewis (2001)], this implies that for the class of magnetic fields considered, Sobolev, Hardy and CLR inequalities hold whenever the magnetic field has no zero mode.


title: Energetically stable discretizations for charge carrier transport and electrokinetic models

abstract: A finite element discretization using a method of lines approached is proposed for approximately solving the Poisson-Nernst-Planck (PNP) equations. This discretization scheme enforces positivity of the computed solutions, corresponding to particle density functions, and a discrete energy estimate is established that resembles the familiar energy law for the PNP system. This energy estimate is extended to finite element solutions to an electrokinetic model, which couples the PNP system with the Navier-Stokes equations. Numerical experiments are conducted to validate convergence of the computed solution and verify the discrete energy estimate.


title: Local $L_\infty$-estimates, weak Harnack inequality, and stochastic continuity of solutions of SPDEs

abstract: We consider stochastic partial differential equations under minimal assumptions: the coefficients are merely bounded and measurable and satisfy the stochastic parabolicity condition. In particular, the diffusion term is allowed to be scaling-critical. We derive local supremum estimates with a stochastic adaptation of De Giorgi's iteration and establish a weak Harnack inequality for the solutions. The latter is then used to obtain pointwise almost sure continuity.


title: Guarding Networks Through Heterogeneous Mobile Guards

abstract: In this article, the issue of guarding multi-agent systems against a sequence of intruder attacks through mobile heterogeneous guards (guards with different ranges) is discussed. The article makes use of graph theoretic abstractions of such systems in which agents are the nodes of a graph and edges represent interconnections between agents. Guards represent specialized mobile agents on specific nodes with capabilities to successfully detect and respond to an attack within their guarding range. Using this abstraction, the article addresses the problem in the context of eternal security problem in graphs. Eternal security refers to securing all the nodes in a graph against an infinite sequence of intruder attacks by a certain minimum number of guards. This paper makes use of heterogeneous guards and addresses all the components of the eternal security problem including the number of guards, their deployment and movement strategies. In the proposed solution, a graph is decomposed into clusters and a guard with appropriate range is then assigned to each cluster. These guards ensure that all nodes within their corresponding cluster are being protected at all times, thereby achieving the eternal security in the graph.


title: Statistical Estimation and Clustering of Group-invariant Orientation Parameters

abstract: We treat the problem of estimation of orientation parameters whose values are invariant to transformations from a spherical symmetry group. Previous work has shown that any such group-invariant distribution must satisfy a restricted finite mixture representation, which allows the orientation parameter to be estimated using an Expectation Maximization (EM) maximum likelihood (ML) estimation algorithm. In this paper, we introduce two parametric models for this spherical symmetry group estimation problem: 1) the hyperbolic Von Mises Fisher (VMF) mixture distribution and 2) the Watson mixture distribution. We also introduce a new EM-ML algorithm for clustering samples that come from mixtures of group-invariant distributions with different parameters. We apply the models to the problem of mean crystal orientation estimation under the spherically symmetric group associated with the crystal form, e.g., cubic or octahedral or hexahedral. Simulations and experiments establish the advantages of the extended EM-VMF and EM-Watson estimators for data acquired by Electron Backscatter Diffraction (EBSD) microscopy of a polycrystalline Nickel alloy sample.


title: Simulation of Genetic Algorithm: Traffic Light Efficiency

abstract: Traffic is a problem in many urban areas worldwide. Traffic flow is dictated by certain devices such as traffic lights. The traffic lights signal when each lane is able to pass through the intersection. Often, static schedules interfere with ideal traffic flow. The purpose of this project was to find a way to make intersections controlled with traffic lights more efficient. This goal was accomplished through the creation of a genetic algorithm, which enhances an input algorithm through genetic principles to produce the fittest algorithm. The program was comprised of two major elements: coding in Java and coding in Simulation of Urban Mobility (SUMO), which is an environment that simulates real traffic. The Java code called upon the SUMO simulation via a command prompt which ran the simulation, received the output, altered the algorithm, and looped. The SUMO component initialized a simulation in which a 1 x 1 street layout was created, each intersection with its own traffic light. Each loop enhanced the input algorithm by altering the scheduling string (dictates the light changes). After the looped simulations were executed, the data was then analyzed. This was accomplished by creating an algorithm based upon regular practice, timed traffic lights, and comparing the output which was comprised of the total time it took for all vehicles to exit the system and the average time it took each individual vehicle to exit the system. These different variables: the time it took the average vehicle to exit the system and total time for all vehicles to exit the system, where then graphed together to provide a visual aid. The genetic algorithm did improve traffic light and traffic flow efficiency in comparison to traditional scheduling methods.


title: Structural Cohesion: Visualization and Heuristics for Fast Computation

abstract: The structural cohesion model is a powerful theoretical conception of cohesion in social groups, but its diffusion in empirical literature has been hampered by operationalization and computational problems. In this paper we start from the classic definition of structural cohesion as the minimum number of actors who need to be removed in a network in order to disconnect it, and extend it by using average node connectivity as a finer grained measure of cohesion. We present useful heuristics for computing structural cohesion that allow a speed-up of one order of magnitude over the algorithms currently available. We analyze three large collaboration networks (co-maintenance of Debian packages, co-authorship in Nuclear Theory and High-Energy Theory) and show how our approach can help researchers measure structural cohesion in relatively large networks. We also introduce a novel graphical representation of the structural cohesion analysis to quickly spot differences across networks.


title: Solar Sector Structure

abstract: The interplanetary magnetic field near has a characteristic "sector" structure that reflects its polarity relative to the solar direction. Typically we observe large-scale coherence in these directions, with two or four "away" or "towards" sectors per solar rotation, from any platform in deep space and near the ecliptic plane. In a simple picture, this morphology simply reflects the idea that the sources of the interplanetary field lie mainly in or near the Sun, and that the solar-wind flow enforces a radial component in this field. Although defined strictly via the interplanetary field near one AU, recent evidence confirms that this pattern also appears clearly at the level of the photosphere, with signatures including not only the large-scale structures (e.g., the streamers) but also highly concentrated fields such as those found in sunspots and even solar flares. This association with small-scale fields strengthens at the Hale sector boundary, defining the Hale boundary as the one for which the polarity switch matches that of the leading-to-following polarity alternation in the sunspots of a given hemisphere.


title: Generation of longitudinal electric current by transversal electromagnetic field in Maxwellian plasmas

abstract: The analysis of nonlinear interaction of transversal electromagnetic field with Maxwellian collisionless classical and quntum plasmas is carried out. Formulas for calculation electric current in Maxwellian collisionless classical and quntum plasmas are deduced. It has appeared, that the nonlinearity account leads to occurrence of the longitudinal electric current directed along a wave vector. This second current is orthogonal to the known transversal current, received at the classical linear analysis. Graphic comparison of density of electric current for classical Maxwellian plasmas and Fermi---Dirac plasmas (plasmas with any degree of degeneration of electronic gas) is carried out. Graphic comparison of density of electric current for classical and quantum Maxwellian plasmas is carried out. Also comparison of dependence of density of electric current of quantum Maxwellian plasmas from dimensionless wave number at various values of dimensionless frequency of oscillations of electromagnetic field is carried out.


title: Entanglement Entropy Near Kondo-Destruction Quantum Critical Points

abstract: We study the impurity entanglement entropy $S_e$ in quantum impurity models that feature a Kondo-destruction quantum critical point (QCP) arising from a pseudogap in the conduction-band density of states or from coupling to a bosonic bath. On the local-moment (Kondo-destroyed) side of the QCP, the entanglement entropy contains a critical component that can be related to the order parameter characterizing the quantum phase transition. In Kondo models describing a spin-$\Simp$, $S_e$ assumes its maximal value of $\ln(2\Simp+1)$ at the QCP and throughout the Kondo phase, independent of features such as particle-hole symmetry and under- or over-screening. In Anderson models, $S_e$ is nonuniversal at the QCP, and at particle-hole symmetry, rises monotonically on passage from the local-moment phase to the Kondo phase; breaking this symmetry can lead to a cusp peak in $S_e$ due to a divergent charge susceptibility at the QCP. Implications of these results for quantum critical systems and quantum dots are discussed.


title: Verbal covering properties of topological spaces

abstract: For any topological space $X$ we study the relation between the universal uniformity $\mathcal U_X$, the universal quasi-uniformity $q\mathcal U_X$ and the universal pre-uniformity $p\mathcal U_X$ on $X$. For a pre-uniformity $\mathcal U$ on a set $X$ and a word $v$ in the two-letter alphabet $\{+,-\}$ we define the verbal power $\mathcal U^v$ of $\mathcal U$ and study its boundedness numbers $\ell(\mathcal U^v)$ and $\bar \ell(\mathcal U^v)$. The boundedness numbers of the (Boolean operations over) the verbal powers of the canonical pre-uniformities $p\mathcal U_X$, $q\mathcal U_X$ and $\mathcal U_X$ yield new cardinal characteristics $\ell^v(X)$, $\bar \ell^v(X)$, $q\ell^v(X)$, $q\bar \ell^v(X)$, $u\ell(X)$ of a topological space $X$, which generalize all known cardinal topological invariants related to (star)-covering properties. We study the relation of the new cardinal invariants $\ell^v$, $\bar \ell^v$ to classical cardinal topological invariants such as Lindel\"of number $\ell$, density $d$, and spread $s$. The simplest new verbal cardinal invariant is the foredensity $\ell^-(X)$ defined for a topological space $X$ as the smallest cardinal $\kappa$ such that for any neighborhood assignment $(O_x)_{x\in X}$ there is a subset $A\subset X$ of cardinality $|A|\le\kappa$ that meets each neighborhood $O_x$, $x\in X$. It is clear that $\ell^-(X)\le d(X)\le \ell^-(X)\cdot \chi(X)$. We shall prove that $\ell^-(X)=d(X)$ if $|X|<\aleph_\omega$. On the other hand, for every singular cardinal $\kappa$ (with $\kappa\le 2^{2^{cf(\kappa)}}$) we construct a (totally disconnected) $T_1$-space $X$ such that $\ell^-(X)=cf(\kappa)<\kappa=|X|=d(X)$.


title: From symplectic groupoids to double structures

abstract: These notes are an introduction to symplectic groupoids and the double structures associated with them. The treatment is intended to lie about midway between the original account of Coste, Dazord and Weinstein, which relied on effective use of the symplectic structures, and the account in my 2005 book, which showed, on the level of Poisson groupoids, that the basic results of the theory follow from `categorical' compatibility conditions between the associated Lie algebroid and Lie groupoid structures. The reader needs to know only the most basic ideas of symplectic geometry, Lie groups, and vector bundles. These notes aim to introduce the reader to certain important, and relatively new, ideas quickly; accordingly they omit much standard material. All the main results here are known, but the approach has a few new features.


title: Classification of isoparametric hypersurfaces in spheres with $(g,m)=(6,1)$

abstract: We classify isoparametric hypersurfaces in spheres with $(g,m)=(6,1)$ and thereby reprove a result of Dorfmeister and Neher.


title: A String-Inspired Model for the Low-$\ell$ CMB

abstract: We present a semi--analytic exploration of some low--$\ell$ angular power spectra inspired by "Brane Supersymmetry Breaking". This mechanism splits Bose and Fermi excitations in String Theory, leaving behind an exponential potential that is just too steep for the inflaton to emerge from the initial singularity while descending it. As a result, the scalar generically bounces against the exponential wall, which typically introduces an infrared depression and a pre--inflationary peak in the power spectrum of scalar perturbations. We elaborate on a possible link between this phenomenon and the low--$\ell$ CMB. For the first 32 multipoles, combining the hard exponential with a milder one leading to $n_s\simeq 0.96$ and with a small gaussian bump we have attained a reduction of $\chi^{\,2}$ to about 46% of the standard $\Lambda$CDM setting, with both WMAP9 and PLANCK 2013 data. This result corresponds to a $\chi^{\,2}/DOF$ of about 0.45, to be compared with a $\Lambda$CDM value of about 0.85. The preferred choices combine naturally quadrupole depression, a first peak around $\ell=5$ and a wide minimum around $\ell=20$. We have also gathered some evidence that similar spectra emerge if the hard exponential is combined with more realistic models of inflation. A problem of the preferred examples is their slow convergence to an almost scale--invariant profile.


title: Separation of transport lifetimes in SrTiO3-based two-dimensional electron liquids

abstract: Deviations from Landau Fermi liquid behavior are ubiquitous features of the normal state of unconventional superconductors. Despite several decades of investigation, the underlying mechanisms of these properties are still not completely understood. In this work, we show that two-dimensional electron liquids at SrTiO3/RTiO3 (R = Gd or Sm) interfaces reveal strikingly similar physics. Analysis of Hall and resistivity data show a clear separation of transport and Hall scattering rates, also known as "two-lifetime" behavior. This framework gives a remarkably simple and general description of the temperature dependence of the Hall coefficient. Distinct transport lifetimes accurately describe the transport phenomena irrespective of the nature of incipient magnetic ordering, the degree of disorder, confinement, or the emergence of non-Fermi liquid behavior. The Hall scattering rate diverges at a critical quantum well thickness, coinciding with a quantum phase transition. Collectively, these results introduce new constraints on the existing microscopic theories of lifetime separation and point to the need for unified understanding.


title: Orthogonal polynomial projection error measured in Sobolev norms in the unit disk

abstract: We study approximation properties of weighted $L^2$-orthogonal projectors onto the space of polynomials of degree less than or equal to $N$ on the unit disk where the weight is of the generalized Gegenbauer form $x \mapsto (1-|x|^2)^\alpha$. The approximation properties are measured in Sobolev-type norms involving canonical weak derivatives, all measured in the same weighted $L^2$ norm. Our basic tool consists in the analysis of orthogonal expansions with respect to Zernike polynomials. The sharpness of the main result is proved in some cases and otherwise strongly hinted at by reported numerical tests. A number of auxiliary results of independent interest are obtained including some properties of the uniformly and non-uniformly weighted Sobolev spaces involved, a Markov-type inequality, connection coefficients between Zernike polynomials and relations between the Fourier-Zernike expansions of a function and its derivatives.


title: The complexity of computing the minimum rank of a sign pattern matrix

abstract: We show that computing the minimum rank of a sign pattern matrix is NP hard. Our proof is based on a simple but useful connection between minimum ranks of sign pattern matrices and the stretchability problem for pseudolines arrangements. In fact, our hardness result shows that it is already hard to determine if the minimum rank of a sign pattern matrix is $\leq 3$. We complement this by giving a polynomial time algorithm for determining if a given sign pattern matrix has minimum rank $\leq 2$. Our result answers one of the open problems from Linial et al. [Combinatorica, 27(4):439--463, 2007].


title: Magnetic structure of the antiferromagnetic half-Heusler compound NdBiPt

abstract: We present results of single crystal neutron diffraction experiments on the rare-earth, half-Heusler antiferromagnet (AFM) NdBiPt. This compound exhibits an AFM phase transition at $T_{\mathrm N}=2.18$~K with an ordered moment of $1.78(9)$~$\mu_{\mathrm B}$ per Nd atom. The magnetic moments are aligned along the $[001]$-direction, arranged in a type-I AFM structure with ferromagnetic planes, alternating antiferromagnetically along a propagation vector $\tau$ of $(100)$. The $R$BiPt ($R$= Ce-Lu) family of materials has been proposed as candidates of a new family of antiferromagnetic topological insulators (AFTI) with magnetic space group that corresponds to a type-II AFM structure where ferromagnetic sheets are stacked along the space diagonal. The resolved structure makes it unlikely, that NdBiPt qualifies as an AFTI.


title: Convergence of the conical Ricci flow on S2 to a soliton

abstract: In our previous work [PSSW], we showed that the Ricci flow on S^2 whose initial metric has conical singularities \sum_{j=1}^k \beta_j[p_j] converges to a constant curvature metric with conic singularities (in the stable and semi-stable cases) or to a gradient shrinking soliton with conical singularities (in the unstable case). The purpose of this note is to show that in the unstable case, that is, the case where \beta_k>\beta_k'=\s_{j<k}\beta_j, that the limiting metric is the unique shrinking soliton with cone singularity \beta_k[p_\infty]+\beta_k'[q_\infty]. This verifies the prediction made in [PSSW].


title: Dephasing due to quasiparticle tunneling in fluxonium qubits: a phenomenological approach

abstract: The fluxonium qubit has arisen as one of the most promising candidate devices for implementing quantum information in superconducting devices, since it is both insensitive to charge noise (like flux qubits) and insensitive to flux noise (like charge qubits). Here, we investigate the stability of the quantum information to quasiparticle tunneling through a Josephson junction. Microscopically, this dephasing is due to the dependence of the quasiparticle transmission probability on the qubit state. We argue that on a phenomenological level the dephasing mechanism can be understood as originating from heat currents, which are flowing in the device due to possible effective temperature gradients, and their sensitivity to the qubit state. The emerging dephasing time is found to be insensitive to the number of junctions with which the superinductance of the fluxonium qubit is realised. Furthermore, we find that the dephasing time increases quadratically with the shunt-inductance of the circuit which highlights the stability of the device to this dephasing mechanism.


title: Two radio-emission mechanisms in PSR J0901-4624

abstract: We have detected sporadic, bright, short-duration radio pulses from PSR J0901$-$4624. These pulses are emitted simultaneously with persistent, periodic emission that dominates the flux density when averaging over many periods of the pulsar. The bright pulses have energies that are consistent with a power-law distribution. The integrated profile of PSR J0901$-$4624 is highly polarized and shows four distinct components. The bright pulses appear to originate near the magnetic pole of the pulsar and have polarization properties unlike that of the underlying emission at the same pulse phase. We conclude that the bright pulses represent a secondary giant-micropulse emission process, possibly from a different region in the pulsar magnetosphere.


abstract: We prove that on any compact complex manifold one can find Gauduchon metrics with prescribed volume form. This is equivalent to prescribing the Chern-Ricci curvature of the metrics, and thus solves a conjecture of Gauduchon from 1984.


title: Dynamic choosability of triangle-free graphs and sparse random graphs

abstract: The \textit{$r$-dynamic choosability} of a graph $G$, written ${\rm ch}_r(G)$, is the least $k$ such that whenever each vertex is assigned a list of at least $k$ colors a proper coloring can be chosen from the lists so that every vertex $v$ has at least $\min\{d_G(v),r\}$ neighbors of distinct colors. Let ${\rm ch}(G)$ denote the choice number of $G$. In this paper, we prove ${\rm ch}_r(G)\leq (1+o(1)){\rm ch}(G)$ when $\frac{\Delta(G)}{\delta(G)}$ is bounded. We also show that there exists a constant $C$ such that for the random graph $G=G(n,p)$ with $\frac{2}{n}<p\leq \frac{1}{2}$, it holds that ${\rm ch}_2(G)\leq {\rm ch}(G) + C$, asymptotically almost surely. Also if $G$ is triangle-free regualr graph, then ${\rm ch}_2(G)\leq {\rm ch}(G)+86$ holds.


title: Semiparametric Bernstein-von Mises Theorem: Second Order Studies

abstract: The major goal of this paper is to study the second order frequentist properties of the marginal posterior distribution of the parametric component in semiparametric Bayesian models, in particular, a second order semiparametric Bernstein-von Mises (BvM) Theorem. Our first contribution is to discover an interesting interference phenomenon between Bayesian estimation and frequentist inferential accuracy: more accurate Bayesian estimation on the nuisance function leads to higher frequentist inferential accuracy on the parametric component. As the second contribution, we propose a new class of dependent priors under which Bayesian inference procedures for the parametric component are not only efficient but also adaptive (w.r.t. the smoothness of nonparametric component) up to the second order frequentist validity. However, commonly used independent priors may even fail to produce a desirable root-n contraction rate for the parametric component in this adaptive case unless some stringent assumption is imposed. Three important classes of semiparametric models are examined, and extensive simulations are also provided.


title: Static and dynamic XY-like short-range order in a frustrated magnet with exchange disorder

abstract: A single crystal of the Co2+ based pyrochlore NaCaCo2F7 was studied by inelastic neutron scattering. This frustrated magnet with quenched exchange disorder remains in a strongly correlated paramagnetic state down to one 60th of the Curie-Weiss temperature. Below T_f = 2.4 K, diffuse elastic scattering develops and comprises 30 +/- 10% of the total magnetic scattering, as expected for J_{eff} = 1/2 moments frozen on a time scale that exceeds \hbar/\delta E=3.8 ps. The diffuse scattering is consistent with short range XY antiferromagnetism with a correlation length of 16 \AA. The momentum (Q) dependence of the inelastic intensity indicates relaxing XY-like antiferromagnetic clusters at energies below ~ 5.5 meV, and collinear antiferromagnetic fluctuations above this energy. The relevant XY configurations form a continuous manifold of symmetry-related states. Contrary to well-known models that produce this continuous manifold, order-by-disorder does not select an ordered state in NaCaCo2F7 despite evidence for weak (~12 %) exchange disorder. Instead, NaCaCo2F7 freezes into short range ordered clusters that span this manifold.


title: Fabrication of (111)-oriented Ca0.5Sr0.5IrO3/SrTiO3 superlattices; a designed playground for honeycomb physics

abstract: We report the fabrication of (111)-oriented superlattice structures with alternating 2m-layers (m = 1, 2, and 3) of Ca0.5Sr0.5IrO3 perovskite and two layers of SrTiO3 perovskite on SrTiO3(111) substrates. In the case of m = 1 bilayer films, the Ir sub-lattice is a buckled honeycomb, where a topological state may be anticipated. The successful growth of superlattice structures on an atomic level along the [111] direction was clearly demonstrated by superlattice reflections in x-ray diffraction patterns and by atomically-resolved transmission electron microscope images. The ground states of the superlattice films were found to be magnetic insulators, which may suggest the importance of electron correlations in Ir perovskites in addition to the much discussed topological effects.


title: A multi-band, multi-level, multi-electron model for efficient FDTD simulations of electromagnetic interactions with semiconductor quantum wells

abstract: We report a new computational model for simulations of electromagnetic interactions with semiconductor quantum well(s) (SQW) in complex electromagnetic geometries using the finite difference time domain (FDTD) method. The presented model is based on an approach of spanning a large number of electron transverse momentum states in each SQW sub-band (multi-band) with a small number of discrete multi-electron states (multi-level, multi-electron). This enables accurate and efficient two dimensional (2-D) and 3-D simulations of nanophotonic devices with SQW active media. The model includes the following features: (1) Optically induced interband transitions between various SQW conduction and heavy-hole or light-hole sub-bands are considered. (2) Novel intra sub-band and inter sub-band transition terms are derived to thermalize the electron and hole occupational distributions to the correct Fermi-Dirac distributions. (3) The terms in (2) result in an explicit update scheme which circumvents numerically cumbersome iterative procedures. This significantly augments computational efficiency. (4) Explicit update terms to account for carrier leakage to unconfined states are derived which thermalize the bulk and SQW populations to a common quasi-equilibrium Fermi-Dirac distribution. (5) Auger recombination and intervalence band absorption are included. The model is validated by comparisons to analytic band filling calculations, simulations of SQW optical gain spectra and photonic crystal lasers.


title: Supercontinuants

abstract: Morier-Genoud, Ovsienko and Tabachnikov introduced supersymmetric frieze patterns (see arXiv:1501.07476). This note gives a solution to Problem 1 from that article: determine the formula for the entries of a superfrieze.


title: Plasma diagnostic potential of 2p4f in N$^+$ -- accurate wavelengths and oscillator strengths

abstract: Radiative emission lines from nitrogen and its ions are often observed in nebulae spectra, where the N$^{2+}$ abundance can be inferred from lines of the 2p4f configuration. In addition, intensity ratios between lines of the 2p3p -- 2p3s and 2p4f -- 2p3d transition arrays can serve as temperature diagnostics. To aid abundance determinations and plasma diagnostics, wavelengths and oscillator strengths were calculated with high-precision for electric-dipole (E1) transitions from levels in the 2p4f configuration of N$^{+}$. Electron correlation and relativistic effects, including the Breit interaction, were systematically taken into account within the framework of the multiconfiguration Dirac-Hartree-Fock (MCDHF) method. Except for the 2p4f - 2p4d transitions with quite large wavelengths and the two-electron-one-photon 2p4f -2s2p$^3$ transitions, the uncertainties of the present calculations were controlled to within 3% and 5% for wavelengths and oscillator strengths, respectively. We also compared our results with other theoretical and experimental values when available. Discrepancies were found between our calculations and previous calculations due to the neglect of relativistic effects in the latter.


title: Cumulative Conditional Expectation Index

abstract: In this paper we study the cumulative conditional expectation function (CCEF) in the copula context. It is shown how to compute CCEF in terms of the cumulative copula function, this natural representation allows to deduce some useful properties, for instance with applications to convex combination of copulas. We introduce approximations of CCEF based on Bernstein polynomial copulas. We introduce estimators for CCEF, which were constructed through Bernstein polynomial estimators for copulas. The estimators are asymptotically normal and biased for CCEF.


title: A Residual Based Sparse Approximate Inverse Preconditioning Procedure for Large Sparse Linear Systems

abstract: The SPAI algorithm, a sparse approximate inverse preconditioning technique for large sparse linear systems, proposed by Grote and Huckle [SIAM J. Sci. Comput., 18 (1997), pp.~838--853.], is based on the F-norm minimization and computes a sparse approximate inverse $M$ of a large sparse matrix $A$ adaptively. However, SPAI may be costly to seek the most profitable indices at each loop and $M$ may be ineffective for preconditioning. In this paper, we propose a residual based sparse approximate inverse preconditioning procedure (RSAI), which, unlike SPAI, is based on only the {\em dominant} rather than all information on the current residual and augments sparsity patterns adaptively during the loops. RSAI is less costly to seek indices and is more effective to capture a good approximate sparsity pattern of $A^{-1}$ than SPAI. To control the sparsity of $M$ and reduce computational cost, we develop a practical RSAI($tol$) algorithm that drops small nonzero entries adaptively during the process. Numerical experiments are reported to demonstrate that RSAI($tol$) is at least competitive with SPAI and can be considerably more efficient and effective than SPAI. They also indicate that RSAI($tol$) is comparable to the PSAI($tol$) algorithm proposed by one of the authors in 2009.


title: Computer Assisted Parallel Program Generation

abstract: Parallel computation is widely employed in scientific researches, engineering activities and product development. Parallel program writing itself is not always a simple task depending on problems solved. Large-scale scientific computing, huge data analyses and precise visualizations, for example, would require parallel computations, and the parallel computing needs the parallelization techniques. In this Chapter a parallel program generation support is discussed, and a computer-assisted parallel program generation system P-NCAS is introduced. Computer assisted problem solving is one of key methods to promote innovations in science and engineering, and contributes to enrich our society and our life toward a programming-free environment in computing science. Problem solving environments (PSE) research activities had started to enhance the programming power in 1970's. The P-NCAS is one of the PSEs; The PSE concept provides an integrated human-friendly computational software and hardware system to solve a target class of problems


title: The Simulation Powers and Limitations of Higher Temperature Hierarchical Self-Assembly Systems

abstract: In this paper, we extend existing results about simulation and intrinsic universality in a model of tile-based self-assembly. Namely, we work within the 2-Handed Assembly Model (2HAM), which is a model of self-assembly in which assemblies are formed by square tiles that are allowed to combine, using glues along their edges, individually or as pairs of arbitrarily large assemblies in a hierarchical manner, and we explore the abilities of these systems to simulate each other when the simulating systems have a higher "temperature" parameter, which is a system wide threshold dictating how many glue bonds must be formed between two assemblies to allow them to combine. It has previously been shown that systems with lower temperatures cannot simulate arbitrary systems with higher temperatures, and also that systems at some higher temperatures can simulate those at particular lower temperatures, creating an infinite set of infinite hierarchies of 2HAM systems with strictly increasing simulation power within each hierarchy. These previous results relied on two different definitions of simulation, one (strong simulation) seemingly more restrictive than the other (standard simulation), but which have previously not been proven to be distinct. Here we prove distinctions between them by first fully characterizing the set of pairs of temperatures such that the high temperature systems are intrinsically universal for the lower temperature systems (i.e. one tile set at the higher temperature can simulate any at the lower) using strong simulation. This includes the first impossibility result for simulation downward in temperature. We then show that lower temperature systems which cannot be simulated by higher temperature systems using the strong definition, can in fact be simulated using the standard definition, proving the distinction between the types of simulation.


title: Superradiance for atoms trapped along a photonic crystal waveguide

abstract: We report observations of superradiance for atoms trapped in the near field of a photonic crystal waveguide (PCW). By fabricating the PCW with a band edge near the D$_1$ transition of atomic cesium, strong interaction is achieved between trapped atoms and guided-mode photons. Following short-pulse excitation, we record the decay of guided-mode emission and find a superradiant emission rate scaling as $\bar{\Gamma}_{\rm SR}\propto\bar{N}\cdot\Gamma_{\rm 1D}$ for average atom number $0.19 \lesssim \bar{N} \lesssim 2.6$ atoms, where $\Gamma_{\rm 1D}/\Gamma_0 =1.1\pm0.1$ is the peak single-atom radiative decay rate into the PCW guided mode and $\Gamma_{0}$ is the Einstein-$A$ coefficient for free space. These advances provide new tools for investigations of photon-mediated atom-atom interactions in the many-body regime.


title: Full-order fluctuation-dissipation relation for a class of non-equilibrium steady states

abstract: Acceleration of relaxation toward a fixed stationary distribution via violation of detailed balance was reported in the context of a Markov chain Monte Carlo method recently. Inspired by this result, systematic methods to violate detailed balance in Langevin dynamics were formulated by using exponential and rotational nonconservative forces. In the present paper, we accentuate that such specific nonconservative forces relate to the large deviation of total heat in an equilibrium state. The response to these nonconservative forces can be described by the intrinsic fluctuation of the total heat in the equilibrium state. Consequently, the fluctuation-dissipation relation for nonequilibrium steady states is derived without recourse to a linear response approximation.


title: Interaction between an emerging flux region and a pre-existing fan-spine dome observed by IRIS and SDO

abstract: We present multi-wavelength observations of a fan-spine dome in the active region NOAA 11996 with the \textit{Interface Region Imaging Spectrograph} (IRIS) and the Atmospheric Imaging Assembly on board the \textit{Solar Dynamics Observatory} (SDO) on March 9, 2014. The destruction of the fan-spine topology owing to the interaction between its magnetic fields and an nearby emerging flux region (EFR) is firstly observed. The line-of-sight magnetograms from the Helioseismic and Magnetic Imager on board the SDO reveal that the dome is located on the mixed magnetic fields, with its rim rooted in the redundant positive polarity surrounding the minority parasitic negative fields. The fan surface of the dome consists of a filament system and recurring jets are observed along its spine. The jet occurring around 13:54 UT is accompanied with a quasi-circular ribbon that brightens in the clockwise direction along the bottom rim of the dome, which may indicate an occurrence of slipping reconnection in the fan-spine topology. The EFR emerges continuously and meets with the magnetic fields of the dome. Magnetic cancellations take place between the emerging negative polarity and the outer positive polarity of the dome's fields, which lead to the rise of the loop connecting the EFR and brightenings related to the dome. A single Gaussian fit to the profiles of the IRIS SI IV 1394 \AA\ line is used in the analysis. It appears that there are two rising components along the slit, except for the rise in the line-of-sight direction. The cancellation process repeats again and again. Eventually the fan-spine dome is destroyed and a new connectivity is formed. We suggest that magnetic reconnection between the EFR and the magnetic fields of the fan-spine dome in the process is responsible for the destruction of the dome.