title: Minimal genus for 4-manifolds with $b^+=1$

abstract: We derive an adjunction inequality for any smooth, closed, connected, oriented 4-manifold $X$ with $b^+=1$. This inequality depends only on the cohomology algebra and generalizes the inequality of Strle in the case of $b_1=0$. We demonstrate that the inequality is especially powerful when $2\tilde \chi+3\sigma\geq 0$, where $\tilde \chi$ is the modified Euler number taking account of the cup product on $H^1$.


title: On the calculation of soft phase space integral

abstract: The recent discovery of the Higgs boson at the LHC attracts much attention to the precise calculation of its production cross section in quantum chromodynamics. In this work, we discuss the calculation of soft triple-emission phase space integral, which is an essential ingredient in the recently calculated soft-virtual corrections to Higgs boson production at next-to-next-to-next-to-leading order. The main techniques used this calculation are method of differential equation for Feynman integral, and integration of harmonic polylogarithms.


title: Topological properties of a Valence-Bond-Solid

abstract: We present a projector quantum Monte Carlo study of the topological properties of the valence-bond-solid ground state in the $J$-$Q_3$ spin model on the square lattice. The winding number is a topological number counting the number of domain walls in the system and is a good quantum number in the thermodynamic limit. We study the finite-size behaviour and obtain the domain wall energy density for a topological nontrivial valence-bond-solid state.


title: Far-from-Equilibrium Distribution from Near-Steady-State Work Fluctuations

abstract: A longstanding goal of nonequilibrium statistical mechanics has been to extend the conceptual power of the Boltzmann distribution to driven systems. We report some new progress towards this goal. Instead of writing the nonequilibrium steady-state distribution in terms of perturbations around thermal equilibrium, we start from the linearized driven dynamics of observables about their stable fixed point, and expand in the strength of the nonlinearities encountered during typical fluctuations away from the fixed point. The first terms in this expansion retain the simplicity of known expansions about equilibrium, but can correctly describe the statistics of a certain class of systems even under strong driving. We illustrate this approach by comparison with a numerical simulation of a sheared Brownian colloid, where we find that the first two terms in our expansion are sufficient to account for the shear thinning behavior at high shear rates.


title: Measurement theory in local quantum physics

abstract: In this paper, we aim to establish foundations of measurement theory in local quantum physics. For this purpose, we discuss a representation theory of completely positive (CP) instruments on arbitrary von Neumann algebras. We introduce a condition called the normal extension property (NEP) and establish a one-to-one correspondence between CP instruments with the NEP and statistical equivalence classes of measuring processes. We show that every CP instrument on an atomic von Neumann algebra has the NEP, extending the well-known result for type I factors. Moreover, we show that every CP instrument on an injective von Neumann algebra is approximated by CP instruments with the NEP. The concept of posterior states is also discussed to show that the NEP is equivalent to the existence of a strongly measurable family of posterior states for every normal state. Two examples of CP instruments without the NEP are obtained from this result. It is thus concluded that in local quantum physics not every CP instrument represents a measuring process, but in most of physically relevant cases every CP instrument can be realized by a measuring process within arbitrary error limits, as every approximately finite dimensional (AFD) von Neumann algebra on a separable Hilbert space is injective. To conclude the paper, the concept of local measurement in algebraic quantum field theory is examined in our framework. In the setting of the Doplicher-Haag-Roberts and Doplicher-Roberts (DHR-DR) theory describing local excitations, we show that an instrument on a local algebra can be extended to a local instrument on the global algebra if and only if it is a CP instrument with the NEP, provided that the split property holds for the net of local algebras.


title: Prediction of Superconductivity in Potassium-Doped Benzene

abstract: To explore underlying mechanism for the superconducting phase in recent discovered aromatic hydrocarbons, we carry out the first-principles calculations on benzene, the basic and the simplest unit of the series and examine the structural and phase stability when doped by potassium, K$_x$C$_6$H$_6$ ($x=1,2,3$). We find that K$_2$C$_6$H$_6$ with the space group of $Pbca$ is the most stable phase with superconducting transition temperature around 6.2 K. Moreover, we argue that all existing hydrocarbons should have a unified superconducting phase in the same temperature range of 5$-$7 K, when doped by two potassium atoms. Our results indicate that the electron-phonon interaction is enough to account for the superconductivity of this unified superconducting phase.


title: Quantum delocalization of protons in the hydrogen bond network of an enzyme active site

abstract: Enzymes utilize protein architectures to create highly specialized structural motifs that can greatly enhance the rates of complex chemical transformations. Here we use experiments, combined with ab initio simulations that exactly include nuclear quantum effects, to show that a triad of strongly hydrogen bonded tyrosine residues within the active site of the enzyme ketosteroid isomerase (KSI) facilitates quantum proton delocalization. This delocalization dramatically stabilizes the deprotonation of an active site tyrosine residue, resulting in a very large isotope effect on its acidity. When an intermediate analog is docked, it is incorporated into the hydrogen bond network, giving rise to extended quantum proton delocalization in the active site. These results shed light on the role of nuclear quantum effects in the hydrogen bond network that stabilizes the reactive intermediate of KSI, and the behavior of protons in biological systems containing strong hydrogen bonds.


title: Combined 3D PET and Optical Projection Tomography Techniques for Plant Root Phenotyping

abstract: New imaging techniques are in great demand for investigating underground plant roots systems which play an important role in crop production. Compared with other non-destructive imaging modalities, PET can image plant roots in natural soil and produce dynamic 3D functional images which reveal the temporal dynamics of plant-environment interactions. In this study, we combined PET with optical projection tomography (OPT) to evaluate its potential for plant root phenotyping. We used a dedicated high resolution plant PET imager that has a 14 cm transaxial and 10 cm axial field of views, and multi-bed imaging capability. The image resolution is around 1.25 mm using ML-EM reconstruction algorithm. B73 inbred maize seeds were germinated and then grown in a sealed jar with transparent gel-based media. PET scanning started on the day when the first green leaf appeared, and was carried out once a day for 5 days. Each morning, around 10 mCi of 11CO2 was administrated into a custom built plant labeling chamber. After 10 minutes, residual activity was flushed out with fresh air before a 2-h PET scan started. For the OPT imaging, the jar was placed inside an acrylic cubic container filled with water, illuminated with a uniform surface light source, and imaged by a DSLR camera from 72 angles to acquire optical images for OPT reconstruction. The same plant was imaged 3 times a day by the OPT system. Plant roots growth is measured from the optical images. Co-registered PET and optical images indicate that most of the hot spots appeared in later time points of the PET images correspond to the most actively growing root tips. The strong linear correlation between 11C allocation at root tips measured by PET and eventual root growth measured by OPT suggests that we can use PET as a phenotyping tool to measure how a plant makes subterranean carbon allocation decisions in different environmental scenarios.


title: New chaos indicators for systems with extremely small Lyapunov exponents

abstract: We propose new chaos indicators for systems with extremely small positive Lyapunov exponents. These chaos indicators can firstly detect a sharp transition between the Arnold diffusion regime and the Chirikov diffusion regime of the Froeschl\'e map and secondly detect chaoticity in systems with zero Lyapunov exponent such as the Boole transformation and the $S$-unimodal function to characterize sub-exponential diffusions.


title: Domain of attraction for maps tangent to the identity in $\mathbb{C}^2$ with characteristic direction of higher degree

abstract: We study holomorphic fixed point germs in two complex variables that are tangent to the identity and have a degenerate characteristic direction. We show that if that characteristic direction is also a characteristic direction for higher degree terms, is non-degenerate for a higher degree term, and satisfies some additional properties, then there is a domain of attraction on which points converge to the origin along that direction.


title: Weak-star convergence and a polynomial approximation problem

abstract: Let $E$ be an arbitrary subset of the unit circle $T$ and let $f$ be a function defined on $E$. When there exist polynomials $P_n$ which are uniformly bounded by a number $M > 0$ on $T$ and converge (pointwise) to $f$ at each point of $E$? We give a necessary and sufficient description of such functions $f$. The necessity part of our result, in fact, is a classical theorem of S.Ya. Khavinson, while the proof of sufficiency uses the method that has been recently applied in particular in the author's solution of an approximation problem proposed by L. Zalcman.


title: Mechanically-driven spreading of bacterial populations

abstract: The effect of mechanical interactions between cells in the spreading of bacterial populations was investigated in one-dimensional space. A continuum-mechanics approach, comprising cell migration, proliferation, and exclusion processes, was employed to elucidate the dynamics. The consequent nonlinear reaction-diffusion-like equation describes the constitution dynamics of a bacterial population. In this model, bacterial cells were treated as rod-like particles that interact with each other through hard-core repulsion, which introduces the exclusion effect that causes bacterial populations to migrate quickly and at high density. The propagation of bacterial density as a traveling wave front over extended times was also analysed. The analytical and numerical solutions revealed that the front speed was enhanced by the exclusion process, which depended upon the cell-packing fraction. Finally, we qualitatively compared our theoretical results with experimental evidence.


title: On the zero-free polynomial approximation problem

abstract: Let $E$ be a compact set in $\mathbb C$ with connected complement, and let $A(E)$ be the class of all complex continuous function on $E$ that are analytic in the interior $E^0$ of $E$. Let $f \in A(E)$ be zero free on $E^0$. By Mergelyan's theorem $f$ can be uniformly approximated on $E$ by polynomials, but is it possible to realize such approximation by polynomials that are zero-free on $E$? This natural question has been proposed by J. Andersson and P. Gauthier. So far it has been settled for some particular sets $E$. The present paper describes classes of functions for which zero free approximation is possible on an arbitrary $E$.


title: On Berman-Gibbs stability and K-stability of $\mathbb{Q}$-Fano varieties

abstract: The notion of Berman-Gibbs stability was originally introduced by Robert Berman for $\mathbb{Q}$-Fano varieties $X$. We show that the pair $(X, -K_X)$ is K-stable (resp. K-semistable) provided that $X$ is Berman-Gibbs stable (resp. semistable).


title: Normality of Orthogonal and Sympletic Nilpotent Orbit Closures in Positive Characteristic

abstract: In this note we investigate the normality of closures of orthogonal and symplectic nilpotent orbits in positive characteristic. We prove that the closure of such a nilpotent orbit is normal provided that neither type d nor type e minimal irreducible degeneration occurs in the closure, and conversely if the closure is normal, then any type e minimal irreducible degeneration does not occur in it. Here, the minimal irreducible degenerations of a nilpotent orbit are introduced by W. Hesselink in [6] (or see [11] from which we take Table 1 for the complete list of all minimal irreducible degenerations). Our result is a weak version in positive characteristic of [11, Theorem 16.2(ii)], one of the main results of [11] over complex numbers.


title: Kondo lattice and antiferromagnetic behavior in quaternary CeTAl$_4$Si$_2$ (T~=~Rh, Ir) single crystals

abstract: We report the synthesis and the magnetic properties of single crystalline CeRhAl$_4$Si$_2$ and CeIrAl$_4$Si$_2$ and their non magnetic La-analogs. The single crystals of these quaternary compounds were grown using Al-Si binary eutectic as flux. The anisotropic magnetic properties of the cerium compounds were explored in detail by means of magnetic susceptibility, isothermal magnetization, electrical resistivity, magnetoresistivity and heat capacity measurements. Both CeRhAl$_4$Si$_2$ and CeIrAl$_4$Si$_2$ undergo two antiferromagnetic transitions, first from the paramagnetic to an antiferromagnetic state at $T_{\rm N1}$~=~12.6~K and 15.5~K, followed by a second transition at lower temperatures $T_{\rm N2}$~=~9.4~K and 13.8~K, respectively. The paramagnetic susceptibility is highly anisotropic and its temperature dependence in the magnetically ordered state suggests the $c$-axis to be the relatively easy axis of magnetization. Concomitantly, isothermal magnetization at 2~K along the $c$-axis shows a sharp spin-flop transition accompanied by a sizeable hysteresis, while it varies nearly linearly with field along the [100] direction up to the highest field 14~T, of our measurement. The electrical resistivity provides evidence of the Kondo interaction in both compounds, inferred from its $-lnT$ behavior in the paramagnetic region. The heat capacity data confirm the bulk nature of the two magnetic transitions in each compound, and further confirm the presence of Kondo interaction by a reduced value of the entropy associated with the magnetic ordering. From the heat capacity data below 1~K, the coefficient of the linear term in the electronic heat capacity, $\gamma$, is inferred to be 195.6 and 49.4~mJ/mol K$^2$ in CeRhAl$_4$Si$_2$ and CeIrAl$_4$Si$_2$, respectively classifying these materials as moderate heavy fermion compounds.


title: A Tutorial on Inverse Problems for Anomalous Diffusion Processes

abstract: Over the last two decades, anomalous diffusion processes in which the mean squares variance grows slower or faster than that in a Gaussian process have found many applications. At a macroscopic level, these processes are adequately described by fractional differential equations, which involves fractional derivatives in time or/and space. The fractional derivatives describe either history mechanism or long range interactions of particle motions at a microscopic level. The new physics can change dramatically the behavior of the forward problems. Naturally one expects that the new physics will impact related inverse problems in terms of uniqueness, stability, and degree of ill-posedness. The last aspect is especially important from a practical point of view, i.e., stably reconstructing the quantities of interest. In this paper, we employ a formal analytic and numerical way to examine the degree of ill-posedness of several "classical" inverse problems for fractional differential equations involving a Djrbashian-Caputo fractional derivative in either time or space, which represent the fractional analogues of that for classical integral order differential equations. We discuss four inverse problems, i.e., backward fractional diffusion, sideways problem, inverse source problem and inverse potential problem for time fractional diffusion, and inverse Sturm-Liouville problem, Cauchy problem, backward fractional diffusion and sideways problem for space fractional diffusion. It is found that contrary to the wide belief, the influence of anomalous diffusion on the degree of ill-posedness is not definitive: it can either significantly improve or worsen the conditioning of related inverse problems, depending crucially on the specific type of given data and quantity of interest. Further, the study exhibits distinct new features of "fractional" inverse problems.


title: Torsional Alfven Waves in Solar Magnetic Flux Tubes of Axial Symmetry

abstract: Aims: Propagation and energy transfer of torsional Alfv\'en waves in solar magnetic flux tubes of axial symmetry is studied. Methods: An analytical model of a solar magnetic flux tube of axial symmetry is developed by specifying a magnetic flux and deriving general analytical formulae for the equilibrium mass density and a gas pressure. The main advantage of this model is that it can be easily adopted to any axisymmetric magnetic structure. The model is used to simulate numerically the propagation of nonlinear Alfv\'en waves in such 2D flux tubes of axial symmetry embedded in the solar atmosphere. The waves are excited by a localized pulse in the azimuthal component of velocity and launched at the top of the solar photosphere, and they propagate through the solar chromosphere, transition region, and into the solar corona. Results: The results of our numerical simulations reveal a complex scenario of twisted magnetic field lines and flows associated with torsional Alfv\'en waves as well as energy transfer to the magnetoacoustic waves that are triggered by the Alfv\'en waves and are akin to the vertical jet flows. Alfv\'en waves experience about 5 % amplitude reflection at the transition region. Magnetic (velocity) field perturbations experience attenuation (growth) with height is agreement with analytical findings. Kinetic energy of magnetoacoustic waves consists of 25 % of the total energy of Alfv\'en waves. The energy transfer may lead to localized mass transport in the form of vertical jets, as well as to localized heating as slow magnetoacoustic waves are prone to dissipation in the inner corona.


title: An analysis of the L1 Scheme for the subdiffusion equation with nonsmooth data

abstract: The subdiffusion equation with a Caputo fractional derivative of order $\alpha\in(0,1)$ in time arises in a wide variety of practical applications, and it is often adopted to model anomalous subdiffusion processes in heterogeneous media. The L1 scheme is one of the most popular and successful numerical methods for discretizing the Caputo fractional derivative in time. The scheme was analyzed earlier independently by Lin and Xu (2007) and Sun and Wu (2006), and an $O(\tau^{2-\alpha})$ convergence rate was established, under the assumption that the solution is twice continuously differentiable in time. However, in view of the smoothing property of the subdiffusion equation, this regularity condition is restrictive, since it does not hold even for the homogeneous problem with a smooth initial data. In this work, we revisit the error analysis of the scheme, and establish an $O(\tau)$ convergence rate for both smooth and nonsmooth initial data. The analysis is valid for more general sectorial operators. In particular, the L1 scheme is applied to one-dimensional space-time fractional diffusion equations, which involves also a Riemann-Liouville derivative of order $\beta\in(3/2,2)$ in space, and error estimates are provided for the fully discrete scheme. Numerical experiments are provided to verify the sharpness of the error estimates, and robustness of the scheme with respect to data regularity.


title: Spectral Dimension of kappa-deformed space-time

abstract: We investigate the spectral dimension of $\kappa$-space-time using the $\kappa$-deformed diffusion equation. The deformed equation is constructed for two different choices of Laplacians in $n$-dimensional, $\kappa$-deformed Euclidean space-time. We use an approach where the deformed Laplacians are expressed in the commutative space-time itself. Using the perturbative solutions to diffusion equations, we calculate the spectral dimension of $\kappa$-deformed space-time and show that it decreases as the probe length decreases. By introducing a bound on the deformation parameter, spectral dimension is guaranteed to be positive definite. We find that, for one of the choices of the Laplacian, the non-commutative correction to the spectral dimension depends on the topological dimension of the space-time whereas for the other, it is independent of the topological dimension. We have also analysed the dimensional flow for the case where the probe particle has a finite extension, unlike a point particle.


title: Speculative Approximations for Terascale Analytics

abstract: Model calibration is a major challenge faced by the plethora of statistical analytics packages that are increasingly used in Big Data applications. Identifying the optimal model parameters is a time-consuming process that has to be executed from scratch for every dataset/model combination even by experienced data scientists. We argue that the incapacity to evaluate multiple parameter configurations simultaneously and the lack of support to quickly identify sub-optimal configurations are the principal causes. In this paper, we develop two database-inspired techniques for efficient model calibration. Speculative parameter testing applies advanced parallel multi-query processing methods to evaluate several configurations concurrently. The number of configurations is determined adaptively at runtime, while the configurations themselves are extracted from a distribution that is continuously learned following a Bayesian process. Online aggregation is applied to identify sub-optimal configurations early in the processing by incrementally sampling the training dataset and estimating the objective function corresponding to each configuration. We design concurrent online aggregation estimators and define halting conditions to accurately and timely stop the execution. We apply the proposed techniques to distributed gradient descent optimization -- batch and incremental -- for support vector machines and logistic regression models. We implement the resulting solutions in GLADE PF-OLA -- a state-of-the-art Big Data analytics system -- and evaluate their performance over terascale-size synthetic and real datasets. The results confirm that as many as 32 configurations can be evaluated concurrently almost as fast as one, while sub-optimal configurations are detected accurately in as little as a $1/20^{\text{th}}$ fraction of the time.


title: Comment on Quantum paraelectric glass state in SrCu3Ti4O12 [Appl. Phys. Lett.104, 262905 (2014)]

abstract: Recently, Kumar et al. [APL 104, 262905 (2014)] have reported quantum paraelectric glass state in SrCu3Ti4O12. Their analysis is based on the misconception about the effect of quantum fluctuations on paraelectric properties. There are very serious issues with their analysis of permittivity vs T data where authors tried to prove the presence of quantum paraelectricity in SrCu3Ti4O12. In the comment we have clearly affirm that the analysis presented by the authors do not prove the presence of quantum paraelectric behavior and there by the quantum paraelectric glass state. The quantum fluctuations and the frequency dependence of permittivity, showing relaxor like behavior cannot exist simultaneously.


title: Magnetic moment of the positronium ion and virial relations

abstract: We derive a hypervirial relation for the positronium ion, a three-body bound state of two electrons and a positron. It relates expectation values of three operators and reconciles two recently published, seemingly inequivalent formulas for the magnetic moment of the ion. As a by-product, the precision of the leading binding correction is improved.


title: Instabilities and shape variation phase transitions in tubular lipid membranes

abstract: Changes of external parameters in proximity of critical point can increase thermal fluctuations of tubular lipid membrane (TLM) and result in variation of the membrane shape. The phase transitions in the system are shown to be controlled by a single effective parameter, which depends on the pressure difference between inner and outer regions of membrane and the applied stretching force. We determine an interval of the parameter values corresponding to the stability region of the cylindrical shape of TLM and investigate the behavior of the system in the vicinity of critical instabilities, where the cylindrical shape of membrane becomes unstable with respect to thermal fluctuations. The applied boundary conditions strongly influence the behavior of TLM. For example, small negative effective parameter corresponds to chiral shape of TLM only in the case of periodic boundary conditions. We also discuss other three types of phase transitions emerging in the system.


title: Cascade Model of an Anomaly in Blazar Spectra at Very High Energy

abstract: It is well known that the effect of gamma-ray absorption on extragalactic background light (EBL) is weakly expressed in the spectra of some blazars. It is shown that a secondary component generated by electromagnetic cascades might considerably decrease the statistical significance of this anomaly. Observational results indicate the existence of the cascade component in the spectra of extragalactic gamma-ray sources, thus supporting the proposed model.


title: The collision between positronium (Ps) and muonium (Mu)

abstract: The collision between a positronium (Ps) and a muonium (Mu) is studied for the first time using the static-exchange model and considering the system as a four-center Coulomb problem in the center of mass frame. An exact analysis is made to find the s-wave elastic phase-shifts, the scattering-lengths for both singlet and triplet channels, the integrated/total elastic cross section and the quenching cross section due to ortho to para conversion of Ps and the conversion ratio.


title: Single electron tunneling with "slow" insulators

abstract: Usual paradigm in the theory of electron transport is related to the fact that the dielectric permittivity of the insulator is assumed to be constant, no time dispersion. We take into account the "slow" polarization dynamics of the dielectric layers in the tunnel barriers in the fluctuating electric fields induced by single-electron tunneling events and study transport in the single electron transistor (SET). Here "slow" dielectric implies slow compared to the characteristic time scales of the SET charging-discharging effects. We show that for strong enough polarizability, such that the induced charge on the island is comparable with the elementary charge, the transport properties of the SET substantially deviate from the known results of transport theory of SET. In particular, the coulomb blockade is more pronounced at finite temperature, the conductance peaks change their shape and the current-voltage characteristics show the memory-effect (hysteresis). However, in contrast to SETs with ferroelectric tunnel junctions, here the periodicity of the conductance in the gate voltage is not broken, instead the period strongly depends on the polarizability of the gate-dielectric. We uncover the fine structure of the hysteresis-effect where the "large" hysteresis loop may include a number of "smaller" loops. Also we predict the memory effect in the current-voltage characteristics $I(V)$, with $I(V)\neq -I(-V)$.


title: Existence and uniqueness of weak solutions of the compressible spherically symmetric Navier-Stokes equations

abstract: One of the most influential fundamental tools in harmonic analysis is Riesz transform. It maps $L^p$ functions to $L^p$ functions for any $p\in (1,\infty)$ which plays an important role in singular operators. As an application in fluid dynamics, the norm equivalence between $\|\nabla u\|_{L^p}$ and $\|\mbox{div} u\|_{L^p}+\|\mbox{curl} u\|_{L^p}$ is well established for $p\in (1,\infty)$. However, since Riesz operators sent bounded functions only to BMO functions, there is no hope to bound $\|\nabla u\|_{L^\infty}$ in terms of $\|\mbox{div} u\|_{L^\infty}+\|\mbox{curl} u\|_{L^\infty}$. As pointed out by Hoff[{\it SIAM J. Math. Anal.} {\bf 37}(2006), No. 6, 1742-1760], this is the main obstacle to obtain uniqueness of weak solutions for isentropic compressible flows. Fortunately, based on new observations, see Lemma \ref{Riesz}, we derive an exact estimate for $\|\nabla u\|_{L^\infty}\le (2+1/N)\|\mbox{div} u\|_{L^\infty}$ for any N-dimensional radially symmetric vector functions $u$. As a direct application, we give an affirmative answer to the open problem of uniqueness of some weak solutions to the compressible spherically symmetric flows in a bounded ball.


title: Communication-Efficient Distributed Optimization of Self-Concordant Empirical Loss

abstract: We consider distributed convex optimization problems originated from sample average approximation of stochastic optimization, or empirical risk minimization in machine learning. We assume that each machine in the distributed computing system has access to a local empirical loss function, constructed with i.i.d. data sampled from a common distribution. We propose a communication-efficient distributed algorithm to minimize the overall empirical loss, which is the average of the local empirical losses. The algorithm is based on an inexact damped Newton method, where the inexact Newton steps are computed by a distributed preconditioned conjugate gradient method. We analyze its iteration complexity and communication efficiency for minimizing self-concordant empirical loss functions, and discuss the results for distributed ridge regression, logistic regression and binary classification with a smoothed hinge loss. In a standard setting for supervised learning, the required number of communication rounds of the algorithm does not increase with the sample size, and only grows slowly with the number of machines.


title: Bayesian Design of Experiments using Approximate Coordinate Exchange

abstract: The construction of decision-theoretic Bayesian designs for realistically-complex nonlinear models is computationally challenging, as it requires the optimization of analytically intractable expected utility functions over high-dimensional design spaces. We provide the most general solution to date for this problem through a novel approximate coordinate exchange algorithm. This methodology uses a Gaussian process emulator to approximate the expected utility as a function of a single design coordinate in a series of conditional optimization steps. It has flexibility to address problems for any choice of utility function and for a wide range of statistical models with different numbers of variables, numbers of runs and randomization restrictions. In contrast to existing approaches to Bayesian design, the method can find multi-variable designs in large numbers of runs without resorting to asymptotic approximations to the posterior distribution or expected utility. The methodology is demonstrated on a variety of challenging examples of practical importance, including design for pharmacokinetic models and design for mixed models with discrete data. For many of these models, Bayesian designs are not currently available. Comparisons are made to results from the literature, and to designs obtained from asymptotic approximations.


title: On the complexity of finite valued functions

abstract: The essential variables in a finite function $f$ are defined as variables which occur in $f$ and weigh with the values of that function. The number of essential variables is an important measure of complexity for discrete functions. When replacing some variables in a function with constants the resulting functions are called subfunctions, and when replacing all essential variables in a function with constants we obtain an implementation of this function. Such an implementation corresponds with a path in an ordered decision diagram (ODD) of the function which connects the root with a leaf of the diagram. The sets of essential variables in subfunctions of $f$ are called separable in $f$. In this paper we study several properties of separable sets of variables in functions which directly impact on the number of implementations and subfunctions in these functions. We define equivalence relations which classify the functions of $k$-valued logic into classes with same number of implementations, subfunctions or separable sets. These relations induce three transformation groups which are compared with the lattice of all subgroups of restricted affine group (RAG). This allows us to solve several important computational and combinatorial problems.


title: Morse bifurcations of transition states in bimolecular reactions

abstract: The transition states and dividing surfaces used to find rate constants for bimolecular reactions are shown to undergo qualitative changes, known as Morse bifurcations, and to exist for a large range of energies, not just immediately above the critical energy for first connection between reactants and products. Specifically, we consider capture between two molecules and the associated transition states for the case of non-zero angular momentum and general attitudes. The capture between an atom and a diatom, and then a general molecule are presented, providing concrete examples of Morse bifurcations of transition states and dividing surfaces. The reduction of the $n$-body systems representing the reactions is discussed and reviewed with comments on the difficulties associated with choosing appropriate charts and the global geometry of the reduced spaces.

categories: physics.chem-ph math-ph math.DS math.MP nlin.CD


title: Fast Generation of Random Spanning Trees and the Effective Resistance Metric

abstract: We present a new algorithm for generating a uniformly random spanning tree in an undirected graph. Our algorithm samples such a tree in expected $\tilde{O}(m^{4/3})$ time. This improves over the best previously known bound of $\min(\tilde{O}(m\sqrt{n}),O(n^{\omega}))$ -- that follows from the work of Kelner and M\k{a}dry [FOCS'09] and of Colbourn et al. [J. Algorithms'96] -- whenever the input graph is sufficiently sparse. At a high level, our result stems from carefully exploiting the interplay of random spanning trees, random walks, and the notion of effective resistance, as well as from devising a way to algorithmically relate these concepts to the combinatorial structure of the graph. This involves, in particular, establishing a new connection between the effective resistance metric and the cut structure of the underlying graph.


title: Quantum transport in 3D Weyl semimetals: Is there a metal-insulator transition?

abstract: We calculate the transport properties of three-dimensional Weyl fermions in a disordered environment. The resulting conductivity depends only on the Fermi energy and the scattering rate. First we study the conductivity at the spectral node for a fixed scattering rate and obtain a continuous transition from an insulator at weak disorder to a metal at stronger disorder. In the self-consistent Born approximation the scattering rate depends on the Fermi energy. Then it is crucial that the limits of the conductivity for a vanishing Fermi energy and a vanishing scattering rate do not commute. As a result, there is also metallic behavior in the phase with vanishing scattering rate and only a quantum critical point remains as an insulating state.


title: Criterion of Hurwitz equivalence for quasipositive factorizations of 3-braids

abstract: We consider the Hurwitz action on quasipositive factorizations of a 3-braid. In a previous paper, for any given 3-braid we described a certain finite set which contains at least one representative of each orbit. Here we give an algorithm to decide if two elements of this finite set belong to the same orbit.


title: Single and double polarization observables in timelike Compton scattering off proton

abstract: We study the exclusive photoproduction of a lepton pair off the proton with the aim of studying the proton quark structure via the Generalized Parton Distributions (GPD) formalism. After deriving the amplitudes of the processes contributing to the $\gamma P\to P' e^+e^-$, the Timelike Compton Scattering and the Bethe-Heitler process, we calculate all unpolarized, single- and double- beam-target spin observables in the valence region in terms of GPDs.


title: Takayasu cofibrations revisited

abstract: This note gives a new construction of Takayasu's cofibrations.


title: Intra-cavity gain shaping of mode-locked Ti:Sapphire laser oscillations

abstract: The gain properties of an oscillator strongly affect its behavior. When the gain is homogeneous, different modes compete for gain resources in a `winner takes all' manner, whereas with inhomogeneous gain, modes can coexist if they utilize different gain resources. We demonstrate precise control over the mode competition in a mode locked Ti:sapphire oscillator by manipulation and spectral shaping of the gain properties, thus steering the competition towards a desired, otherwise inaccessible, oscillation. Specifically, by adding a small amount of spectrally shaped inhomogeneous gain to the standard homogeneous gain oscillator, we selectively enhance a desired two-color oscillation, which is inherently unstable to mode competition and could not exist in a purely homogeneous gain oscillator. By tuning the parameters of the additional inhomogeneous gain we flexibly control the center wavelengths, relative intensities and widths of the two colors.


title: On the spot-futures no-arbitrage relations in commodity markets

abstract: In commodity markets the convergence of futures towards spot prices, at the expiration of the contract, is usually justified by no-arbitrage arguments. In this article, we propose an alternative approach that relies on the expected profit maximization problem of an agent, producing and storing a commodity while trading in the associated futures contracts. In this framework, the relation between the spot and the futures prices holds through the well-posedness of the maximization problem. We show that the futures price can still be seen as the risk-neutral expectation of the spot price at maturity and we propose an explicit formula for the forward volatility. Moreover, we provide an heuristic analysis of the optimal solution for the production/storage/trading problem, in a Markovian setting. This approach is particularly interesting in the case of energy commodities, like electricity: this framework indeed remains suitable for commodities characterized by storability constraints, when standard no-arbitrage arguments cannot be safely applied.


title: Orthogonal jumps of wavefunction in white-noise potentials

abstract: We investigate the evolution of the quantum state for a free particle placed into a random external potential of white-noise type. The master equation for the density matrix is derived by means of path integral method. We propose an equivalent stochastic process where the wavefunction satisfies a nonlinear Schr\"odinger equation except for random moments at which it shows orthogonal jumps. The relation of our work to the usual theory of quantum noise and damping is briefly discussed.


title: Two examples of eigenforms of the Hodge-de Rham Laplacian

abstract: In the present paper, we consider the Hodge-de Rham Laplacian that acts on conformal Killing and projective Killing one-forms of a compact Riemannian manifold.


title: The theta-exact Seiberg-Witten maps at the e3 order

abstract: We study wo distinct theta-exact Seiberg-Witten (SW) map expansions, (I) and (II) respectively, up to the e3 order for the gauge parameter, gauge field and the gauge field strengths of the noncommutative U*(1) gauge theory on the Moyal space. We derive explicitly the closed form expression for the SW map ambiguity between the two and observe the emergence of several new totally commutative generalized star products. We also identify the additional gauge freedoms within each of the e3 order field strength expansions and define corresponding sets of deformation/ratio/weight parameters, (kappa,kappa_i) and (kappa,kappa'_i), for these two SW maps respectively.


title: Star Formation in Disk Galaxies. III. Does stellar feedback result in cloud death?

abstract: Stellar feedback, star formation and gravitational interactions are major controlling forces in the evolution of Giant Molecular Clouds (GMCs). To explore their relative roles, we examine the properties and evolution of GMCs forming in an isolated galactic disk simulation that includes both localised thermal feedback and photoelectric heating. The results are compared with the three previous simulations in this series which consists of a model with no star formation, star formation but no form of feedback and star formation with photoelectric heating in a set with steadily increasing physical effects. We find that the addition of localised thermal feedback greatly suppresses star formation but does not destroy the surrounding GMC, giving cloud properties closely resembling the run in which no stellar physics is included. The outflows from the feedback reduce the mass of the cloud but do not destroy it, allowing the cloud to survive its stellar children. This suggests that weak thermal feedback such as the lower bound expected for supernova may play a relatively minor role in the galactic structure of quiescent Milky Way-type galaxies, compared to gravitational interactions and disk shear.


title: A search for optical counterparts of the complex Vela X system

abstract: The pulsar wind nebula (PWN) associated with the Vela pulsar is a bright source in the radio, X-ray and gamma-ray bands, but not in the optical. This source is very near, lying at a distance of 290 pc, as inferred from the radio and optical parallax measurements of the pulsar. Knowledge of the brightness and structure of the Vela PWN in optical is important in order to constrain the underlying particle spectrum (and possibly the B-field properties and particle losses) associated with this extended source. We use results from the Digital Sky Survey, as well as results obtained using the SAAO 1.0 m telescope equipped with an imaging CCD (STE4) and BV filters, in an attempt to measure optical radiation from Vela X. To enlarge our field of view, we constructed a mosaic consisting of 3 x 3 frames around the pulsar position. We present spectral measurements from the High Energy Stereoscopic System (H.E.S.S.), Fermi Large Area Telescope (LAT), ASCA, Hubble Space Telescope (HST), Very Large Telescope (VLT), New Technology Telescope (NTT), and Wilkinson Microwave Anisotropy Probe (WMAP), in addition to our optical results. Using these data, we investigate whether or not the radio synchrotron component can be smoothly extrapolated to the optical band. This would constrain the electron population to consist of either a single or multiple components, representing a significant advancement in our understanding of this complex multiwavelength source.


title: Extremal functions for Trudinger-Moser inequalities of Adimurthi-Druet type in dimension two

abstract: Combining Carleson-Chang's result with blow-up analysis, we prove existence of extremal functions for certain Trudinger-Moser inequalities in dimension two. This kind of inequality was originally proposed by Adimurthi and O. Druet, extended by the author to high dimensional case and Riemannian surface case, generalized by C. Tintarev to wider cases including singular form and by M. de Souza and J. M. do \'O to the whole Euclidean space R^2. In addition to the Euclidean case, we also consider the Riemannian surface case. The results in the current paper complement that of L. Carleson and A. Chang, M. Struwe, M. Flucher, K. Lin, and Adimurthi-Druet, our previous ones, and part of C. Tintarev.


title: Spline Galerkin methods for the Sherman-Lauricella equation on contours with corners

abstract: Spline Galerkin approximation methods for the Sherman-Lauricella integral equation on simple closed piecewise smooth contours are studied, and necessary and sufficient conditions for their stability are obtained. It is shown that the method under consideration is stable if and only if certain operators associated with the corner points of the contour are invertible. Numerical experiments demonstrate a good convergence of the spline Galerkin methods and validate theoretical results. Moreover, it is shown that if all corners of the contour have opening angles located in interval $(0.1\pi, 1.9\pi)$, then the corresponding Galerkin method based on splines of order $0$, $1$ and $2$ is always stable. These results are in strong contrast with the behaviour of the Nystr\"om method, which has a number of instability angles in the interval mentioned.


title: Observing Radio Pulsars in the Galactic Centre with the Square Kilometre Array

abstract: The discovery and timing of radio pulsars within the Galactic centre is a fundamental aspect of the SKA Science Case, responding to the topic of "Strong Field Tests of Gravity with Pulsars and Black Holes" (Kramer et al. 2004; Cordes et al. 2004). Pulsars have in many ways proven to be excellent tools for testing the General theory of Relativity and alternative gravity theories (see Wex (2014) for a recent review). Timing a pulsar in orbit around a companion, provides a unique way of probing the relativistic dynamics and spacetime of such a system. The strictest tests of gravity, in strong field conditions, are expected to come from a pulsar orbiting a black hole. In this sense, a pulsar in a close orbit ($P_{\rm orb}$ < 1 yr) around our nearest supermassive black hole candidate, Sagittarius A* - at a distance of ~8.3 kpc in the Galactic centre (Gillessen et al. 2009a) - would be the ideal tool. Given the size of the orbit and the relativistic effects associated with it, even a slowly spinning pulsar would allow the black hole spacetime to be explored in great detail (Liu et al. 2012). For example, measurement of the frame dragging caused by the rotation of the supermassive black hole, would allow a test of the "cosmic censorship conjecture." The "no-hair theorem" can be tested by measuring the quadrupole moment of the black hole. These are two of the prime examples for the fundamental studies of gravity one could do with a pulsar around Sagittarius A*. As will be shown here, SKA1-MID and ultimately the SKA will provide the opportunity to begin to find and time the pulsars in this extreme environment.


title: CMB Cold Spot from Inflationary Feature Scattering

abstract: We propose a "feature-scattering" mechanism to explain the cosmic microwave background cold spot seen from WMAP and Planck maps. If there are hidden features in the potential of multi-field inflation, the inflationary trajectory can be scattered by such features. The scattering is controlled by the amount of isocurvature fluctuations, and thus can be considered as a mechanism to convert isocurvature fluctuations into curvature fluctuations. This mechanism predicts localized cold spots (instead of hot ones) on the CMB. In addition, it may also bridge a connection between the cold spot and a dip on the CMB power spectrum at $\ell \sim 20$.


title: A categorification of twisted Heisenberg algebras

abstract: We categorify a quantized Heisenberg algebra associated to a finite subgroup of SL(2,C).


title: Model dependence of the multi-transonic behavior, stability properties and corresponding acoustic geometry for accretion onto a spinning black hole

abstract: Multi-transonic accretion for a spinning black hole has been compared among different disc geometries within post Newtonian pseudo potential framework. The variation of stationary shock characteristics with black hole spin has been studied in details for all the disc models and compared for adiabatic as well as for isothermal scenario. The variations of surface gravity with spin for all these cases have also been investigated.


title: Axiomatization of geometry employing group actions

abstract: The aim of this paper is to develop a new axiomatization of planar geometry by reinterpreting the original axioms of Euclid. The basic concept is still that of a line segment but its equivalent notion of betweenness is viewed as a topological, not a metric concept. That leads quickly to the notion of connectedness without any need to dwell on the definition of topology. In our approach line segments must be connected. Lines and planes are unified via the concept of separation: lines are separated into two components by each point, planes contain lines that separate them into two components as well. We add a subgroup of bijections preserving line segments and establishing unique isomorphism of basic geometrical sets, and the axiomatic structure is complete. Of fundamental importance is the Fixed Point Theorem that allows for creation of the concepts of length and congruency of line segments. The resulting structure is much more in sync with modern science than other axiomatic approaches to planar geometry. For instance, it leads naturally to the Erlangen Program in geometry. Our Conditions of Homogeneity and Rigidity have two interpretations. In physics, they correspond to the basic tenet that independent observers should arrive at the same measurement and are related to boosts in special relativity. In geometry, they mean uniqueness of congruence for certain geometrical figures. Another thread of the paper is the introduction of boundary at infinity, an important concept of modern mathematics, and linking of Pasch Axiom to endowing boundaries at infinity with a natural relation of betweenness. That way spherical geometry can be viewed as geometry of boundaries at infinity.


title: Muon spin rotation and neutron scattering study of the non-centrosymmetric tetragonal compound CeAuAl3

abstract: We have investigated the non-centrosymmetric tetragonal heavy-fermion compound CeAuAl3 using muon spin rotation (muSR), neutron diffraction (ND) and inelastic neutron scattering (INS) measurements. We have also revisited the magnetic, transport and thermal properties. The magnetic susceptibility reveals an antiferromagnetic transition at 1.1 K with a possibility of another magnetic transition near 0.18 K. The heat capacity shows a sharp lambda-type anomaly at 1.1 K in zero-filed, which broadens and moves to higher temperature in applied magnetic field. Our zero-field muSR and ND measurements confirm the existence of a long-range magnetic ground state below 1.2 K. Further the ND study reveals an incommensurate magnetic ordering with a magnetic propagation vector k = (0, 0, 0.52) and a spiral structure of Ce moments coupled ferromagnetically within the ab-plane. Our INS study reveals the presence of two well-defined crystal electric field (CEF) excitations at 5.1 meV and 24.6 meV in the paramagnetic phase of CeAuAl3 which can be explained on the basis of the CEF theory. Furthermore, low energy quasi-elastic excitations show a Gaussian line shape below 30 K compared to a Lorentzian line shape above 30 K, indicating a slowdown of spin fluctuation below 30 K. We have estimated a Kondo temperature of TK=3.5 K from the quasi-elastic linewidth, which is in good agreement with that estimated from the heat capacity. This study also indicates the absence of any CEF-phonon coupling unlike that observed in isostructural CeCuAl3. The CEF parameters, energy level scheme and their wave functions obtained from the analysis of INS data explain satisfactorily the single crystal susceptibility in the presence of two-ion anisotropic exchange interaction in CeAuAl3.


title: Consistent Classification Algorithms for Multi-class Non-Decomposable Performance Metrics

abstract: We study consistency of learning algorithms for a multi-class performance metric that is a non-decomposable function of the confusion matrix of a classifier and cannot be expressed as a sum of losses on individual data points; examples of such performance metrics include the macro F-measure popular in information retrieval and the G-mean metric used in class-imbalanced problems. While there has been much work in recent years in understanding the consistency properties of learning algorithms for `binary' non-decomposable metrics, little is known either about the form of the optimal classifier for a general multi-class non-decomposable metric, or about how these learning algorithms generalize to the multi-class case. In this paper, we provide a unified framework for analysing a multi-class non-decomposable performance metric, where the problem of finding the optimal classifier for the performance metric is viewed as an optimization problem over the space of all confusion matrices achievable under the given distribution. Using this framework, we show that (under a continuous distribution) the optimal classifier for a multi-class performance metric can be obtained as the solution of a cost-sensitive classification problem, thus generalizing several previous results on specific binary non-decomposable metrics. We then design a consistent learning algorithm for concave multi-class performance metrics that proceeds via a sequence of cost-sensitive classification problems, and can be seen as applying the conditional gradient (CG) optimization method over the space of feasible confusion matrices. To our knowledge, this is the first efficient learning algorithm (whose running time is polynomial in the number of classes) that is consistent for a large family of multi-class non-decomposable metrics. Our consistency proof uses a novel technique based on the convergence analysis of the CG method.


title: LP approximations to mixed-integer polynomial optimization problems

abstract: We present a class of linear programming approximations for constrained optimization problems. In the case of mixed-integer polynomial optimization problems, if the intersection graph of the constraints has bounded tree-width our construction yields a class of linear size formulations that attain any desired tolerance. As a result, we obtain an approximation scheme for the "AC-OPF" problem on graphs with bounded tree-width. We also describe a more general construction for pure binary optimization problems where individual constraints are available through a membership oracle; if the intersection graph for the constraints has bounded tree-width our construction is of linear size and exact. This improves on a number of results in the literature, both from the perspective of formulation size and generality.


title: The role of pion exchange in $\eta$ meson photo production on the deuteron

abstract: Incoherent $\eta$ photoproduction on the deuteron is studied with the main emphasis on the role of final state interactions. In addition to the previously studied mechanisms of $NN$ and $\eta N$ rescatterings, the role of an intermediate pion exchange is considered in detail, where first a pion is photoproduced on one nucleon and then rescatters into an $\eta$ meson on the other, the spectator nucleon. It is found, that the role of this pion mediated contribution is comparable in size to that of $\eta N$ rescattering. Results for total and semi-inclusive differential cross sections and associated polarization observables are presented. In particular polarization observables show a significant sensititvity to final state interactions.


title: Chiral Thirring-Wess model with Faddeevian Regularization

abstract: Replacing vector type of interaction of the Thirring-Wess model by the chiral type a new model is presented which is termed here as chiral Thirring-Wess model. Ambiguity parameters of regularization is so chosen that the model falls into the Faddeevian class. The resulting Faddeevian class of model in general do not possess Lorentz invariance. However we can exploit the arbitrariness admissible in the ambiguity parameters to relate the quantum mechanically generated ambiguity parameters with the classical parameter involved in the masslike term of the gauge field which helps to maintain physical Lorentz invariance instead of the absence of manifestly lorentz covariance of the model. The the phase space structure and the theoretical spectrum of this class of model has been determined through Dirac's method of quantization of constraint system.


title: Motivation and challenge to capture both large scale and local transport in next generation accretion theory

abstract: Accretion disc theory is less developed than stellar evolution theory although a similarly mature phenomenological picture is ultimately desired. While the interplay of theory and numerical simulations has amplified community awareness of the role of magnetic fields in angular momentum transport, there remains a long term challenge to incorporate insight gained from simulations back into improving practical models for comparison with observations. Here we emphasize the need to incorporate the role of non-local transport more precisely. To show where large scale transport would fit into the theoretical framework and how it is currently missing, we review why the wonderfully practical approach of Shakura-Sunyaev (1973,SS73) is necessarily a mean field theory, and one which does not include large scale transport. Observations of coronae and jets combined with the interpretation of results even from shearing box simulations of the magnetorotational instability (MRI) suggest that a significant fraction of disc transport is indeed non-local. We show that the Maxwell stresses in saturation are dominated by large scale contributions and the physics of MRI transport is not fully captured by a viscosity. We also clarify the standard physical interpretation of the MRi as it applies to shearing boxes. Computational limitations have so far focused most attention toward local simulations but the next generation of global simulations should help to inform improved mean field theories. Mean field accretion theory and mean field dynamo theory should in fact be unified into a single theory that predicts the time evolution of spectra and luminosity from separate disc, corona, and outflow contributions. Finally, we note that any mean field theory has a finite predictive precision that needs to be quantified when comparing the predictions to observations.


title: Computation of the multi-chord distribution of convex and concave polygons

abstract: Analytical expressions for the distribution of the length of chords corresponding to the affine invariant measure on the set of chords are given for convex polygons. These analytical expressions are a computational improvement over other expressions published in 2011. The correlation function of convex polygons can be computed from the results obtained in this work, because it is determined by the distribution of chords. An analytical expression for the multi-chord distribution of the length of chords corresponding to the affine invariant measure on the set of chords is found for non convex polygons. In addition we give an algorithm to find this multi-chord distribution which, for many concave polygons, is computationally more efficient than the said analytical expression. The results also apply to non simply connected polygons.


title: Continuous-time limit of dynamic games with incomplete information and a more informed player

abstract: We study a two-player, zero-sum, dynamic game with incomplete information where one of the players is more informed than his opponent. We analyze the limit value as the players play more and more frequently. The more informed player observes the realization of a Markov process (X,Y) on which the payoffs depend, while the less informed player only observes Y and his opponent's actions. We show the existence of a limit value as the time span between two consecutive stages goes to zero. This value is characterized through an auxiliary optimization problem and as the unique viscosity solution of a second order Hamilton-Jacobi equation with convexity constraints.


title: The Milnor-Thurston determinant and the Ruelle transfer operator

abstract: The topological entropy $h_{\rm top}$ of a continuous piecewise monotone interval map measures the exponential growth in the number of monotonicity intervals for iterates of the map. Milnor and Thurston showed that $\exp(-h_{\rm top})$ is the smallest zero of an analytic function, now coined the Milnor-Thurston determinant, that keeps track of relative positions of forward orbits of critical points. On the other hand $\exp(h_{\rm top})$ equals the spectral radius of a Ruelle transfer operator $L$, associated with the map. Iterates of $L$ keep track of inverse orbits of the map. For no obvious reason, a Fredholm determinant for the transfer operator has not only the same leading zero as the M-T determinant but all peripheral (those lying in the unit disk) zeros are the same. The purpose of this note is to show that on a suitable function space, the dual of the Ruelle transfer operator has a regularized determinant, identical to the Milnor-Thurston determinant, hereby providing a natural explanation for the above puzzle.


title: Lifting curves simply

abstract: We provide linear lower bounds for $f_\rho(L)$, the smallest integer so that every curve on a fixed hyperbolic surface $(S,\rho)$ of length at most $L$ lifts to a simple curve on a cover of degree at most $f_\rho(L)$. This bound is independent of hyperbolic structure $\rho$, and improves on a recent bound of Gupta-Kapovich. When $(S,\rho)$ is without punctures, using work of Patel we conclude asymptotically linear growth of $f_\rho$. When $(S,\rho)$ has a puncture, we obtain exponential lower bounds for $f_\rho$.


title: Constraining the Radiation and Plasma Environment of the Kepler Circumbinary Habitable Zone Planets

abstract: The discovery of many planets using the Kepler telescope includes ten planets orbiting eight binary stars. Three binaries, Kepler-16, Kepler-47, and Kepler-453, have at least one planet in the circumbinary habitable-zone (BHZ). We constrain the level of high-energy radiation and the plasma environment in the BHZ of these systems. With this aim, BHZ limits in these Kepler binaries are calculated as a function of time, and the habitability lifetimes are estimated for hypothetical terrestrial planets and/or moons within the BHZ. With the time-dependent BHZ limits established, a self-consistent model is developed describing the evolution of stellar activity and radiation properties as proxies for stellar aggression toward planetary atmospheres. Modeling binary stellar rotation evolution, including the effect of tidal interaction between stars in binaries is key to establishing the environment around these systems. We find that Kepler-16 and its binary analogs provide a plasma environment favorable for the survival of atmospheres of putative Mars-sized planets and exomoons. Tides have modified the rotation of the stars in Kepler-47 making its radiation environment less harsh in comparison to the solar system. This is a good example of the mechanism first proposed by Mason et al. Kepler-453 has an environment similar to that of the solar system with slightly better than Earth radiation conditions at the inner edge of the BHZ. These results can be reproduced and even reparametrized as stellar evolution and binary tidal models progress, using our online tool http://bhmcalc.net.


title: Complete Homology over associative rings

abstract: We compare two generalizations of Tate homology: stable homology and the J-completion of Tor, also known as complete homology. For finitely generated modules, we show that the two theories agree over Artin algebras and over commutative noetherian rings that are Gorenstein, or local and complete.


title: On the Gromov width of polygon spaces

abstract: For generic $r=(r_1,\ldots,r_n) \in \mathbb{R}^n_+$ the space $\mathcal{M}(r)$ of $n$--gons in $\mathbb{R}^3$ with edges of lengths $r$ is a smooth, symplectic manifold. We investigate its Gromov width and prove that the expression $$2\pi \min \{2 r_j, (\sum_{i \neq j} r_i) - r_j\,\,|\, j=1,\ldots,n\}$$ is the Gromov width of all (smooth) $5$--gon spaces and of $6$--gon spaces, under some condition on $r \in \mathbb{R}^6_+$. The same formula constitutes a lower bound for all (smooth) spaces of $6$--gons. Moreover, we prove that the Gromov width of $\mathcal{M}(r)$ is given by the above expression when $\mathcal{M}(r)$ is symplectomorphic to $\mathbb{C}\mathbb{P}^{n-3}$, for any $n \geq 4$.


title: Sequence Modeling using Gated Recurrent Neural Networks

abstract: In this paper, we have used Recurrent Neural Networks to capture and model human motion data and generate motions by prediction of the next immediate data point at each time-step. Our RNN is armed with recently proposed Gated Recurrent Units which has shown promising results in some sequence modeling problems such as Machine Translation and Speech Synthesis. We demonstrate that this model is able to capture long-term dependencies in data and generate realistic motions.


title: An almost flat manifold with a cyclic or quaternionic holonomy group bounds

abstract: A long-standing conjecture of Farrell and Zdravkovska and independently S.~T.~Yau states that every almost flat manifold is the boundary of a compact manifold. This paper gives a simple proof of this conjecture when the holonomy group is cyclic or quaternionic. The proof is based on the interaction between flat bundles and involutions.


title: A Liouville type result for bounded, entire solutions to a class of variational semilinear elliptic systems

abstract: We prove a Liouville type result for bounded, entire solutions to a class of variational semilinear elliptic systems, based on the growth of their potential energy over balls with growing radius. Important special cases to which our result applies are the Ginzburg-Landau system and systems that arise in the study of multi-phase transitions.


title: An Event-Driven Approach for Studying Gene Block Evolution in Bacteria

abstract: Motivation: Gene blocks are genes co-located on the chromosome. In many cases, genes blocks are conserved between bacterial species, sometimes as operons, when genes are co-transcribed. The conservation is rarely absolute: gene loss, gain, duplication, block splitting, and block fusion are frequently observed. An open question in bacterial molecular evolution is that of the formation and breakup of gene blocks, for which several models have been proposed. These models, however, are not generally applicable to all types of gene blocks, and consequently cannot be used to broadly compare and study gene block evolution. To address this problem we introduce an event-based method for tracking gene block evolution in bacteria. Results: We show here that the evolution of gene blocks in proteobacteria can be described by a small set of events. Those include the insertion of genes into, or the splitting of genes out of a gene block, gene loss, and gene duplication. We show how the event-based method of gene block evolution allows us to determine the evolutionary rate, and to trace the ancestral states of their formation. We conclude that the event-based method can be used to help us understand the formation of these important bacterial genomic structures. Availability: The software is available under GPLv3 license on http://github.com/reamdc1/gene_block_evolution.git Supplementary online material: http://iddo-friedberg.net/operon-evolution Contact: Iddo Friedberg i.friedberg@miamioh.edu


title: Inferring Latent States and Refining Force Estimates via Hierarchical Dirichlet Process Modeling in Single Particle Tracking Experiments

abstract: Optical microscopy provides rich spatio-temporal information characterizing in vivo molecular motion. However, effective forces and other parameters used to summarize molecular motion change over time in live cells due to latent state changes, e.g., changes induced by dynamic micro-environments, photobleaching, and other heterogeneity inherent in biological processes. This study focuses on techniques for analyzing Single Particle Tracking (SPT) data experiencing abrupt state changes. We demonstrate the approach on GFP tagged chromatids experiencing metaphase in yeast cells and probe the effective forces resulting from dynamic interactions that reflect the sum of a number of physical phenomena. State changes are induced by factors such as microtubule dynamics exerting force through the centromere, thermal polymer fluctuations, etc. Simulations are used to demonstrate the relevance of the approach in more general SPT data analyses. Refined force estimates are obtained by adopting and modifying a nonparametric Bayesian modeling technique, the Hierarchical Dirichlet Process Switching Linear Dynamical System (HDP-SLDS), for SPT applications. The HDP-SLDS method shows promise in systematically identifying dynamical regime changes induced by unobserved state changes when the number of underlying states is unknown in advance (a common problem in SPT applications). We expand on the relevance of the HDP-SLDS approach, review the relevant background of Hierarchical Dirichlet Processes, show how to map discrete time HDP-SLDS models to classic SPT models, and discuss limitations of the approach. In addition, we demonstrate new computational techniques for tuning hyperparameters and for checking the statistical consistency of model assumptions directly against individual experimental trajectories; the techniques circumvent the need for "ground-truth" and subjective information.


title: Contact Representations of Graphs in 3D

abstract: We study contact representations of graphs in which vertices are represented by axis-aligned polyhedra in 3D and edges are realized by non-zero area common boundaries between corresponding polyhedra. We show that for every 3-connected planar graph, there exists a simultaneous representation of the graph and its dual with 3D boxes. We give a linear-time algorithm for constructing such a representation. This result extends the existing primal-dual contact representations of planar graphs in 2D using circles and triangles. While contact graphs in 2D directly correspond to planar graphs, we next study representations of non-planar graphs in 3D. In particular we consider representations of optimal 1-planar graphs. A graph is 1-planar if there exists a drawing in the plane where each edge is crossed at most once, and an optimal n-vertex 1-planar graph has the maximum (4n - 8) number of edges. We describe a linear-time algorithm for representing optimal 1-planar graphs without separating 4-cycles with 3D boxes. However, not every optimal 1-planar graph admits a representation with boxes. Hence, we consider contact representations with the next simplest axis-aligned 3D object, L-shaped polyhedra. We provide a quadratic-time algorithm for representing optimal 1-planar graph with L-shaped polyhedra.


title: Massive MIMO and Waveform Design for 5th Generation Wireless Communication Systems

abstract: This article reviews existing related work and identifies the main challenges in the key 5G area at the intersection of waveform design and large-scale multiple antenna systems, also known as Massive MIMO. The property of self-equalization is introduced for Filter Bank Multicarrier (FBMC)-based Massive MIMO, which can reduce the number of subcarriers required by the system. It is also shown that the blind channel tracking property of FBMC can be used to address pilot contamination -- one of the main limiting factors of Massive MIMO systems. Our findings shed light into and motivate for an entirely new research line towards a better understanding of waveform design with emphasis on FBMC-based Massive MIMO networks.


title: Phase structure and Higgs boson mass in a Higgs-Yukawa model with a dimension-6 operator

abstract: We investigate the impact of a $\lambda_6 \varphi^6$ term included in a chirally invariant lattice Higgs-Yukawa model. Such a term could emerge from BSM physics at some larger energy scale. We map out the phase structure of the Higgs-Yukawa model with positive $\lambda_6$ and negative quartic self coupling of the scalar fields. To this end, we evaluate the constraint effective potential in lattice perturbation theory and also determine the magnetization of the model via numerical simulations which allow us to reach also non-perturbative values of the couplings. As a result, we find a complex phase structure with first and second order phase transitions identified through the magnetization. Further we analyze the effect of such a $\varphi^6$ term on the lower Higgs boson mass bound to see, whether the standard model lower mass bound can be altered.


title: Coderivative characterizations of maximal monotonicity for set-valued mappings

abstract: This paper concerns generalized differential characterizations of maximal monotone set-valued mappings. Using advanced tools of variational analysis, we establish coderivative criteria for maximal monotonicity of set-valued mappings, which seem to be the first infinitesimal characterizations of maximal monotonicity outside the single-valued case. We also present second-order necessary and sufficient conditions for lower-${\mathcal C}^2$ functions to be convex and strongly convex. Examples are provided to illustrate the obtained results and the imposed assumptions.


title: Non-diagonal metric on a product riemanniann manifold

abstract: In this paper, We construct the symmetric tensor field $G_{f_1f_2}$ and $h_{f_1f_2}$ on a product manifold and we give conditions under which $G_{f_1f_2}$ becomes a metric tensor, theses tensors fields will be called the generalized warped product, and then we develop an expression of curvature for the connection of the generalized warped product in relation to those corresponding analogues of its base and fiber and warping functions. By constructing a frame field in $M_1\times_{f_1f_2}M_2$ with respect to the Riemannian metric $G_{f_1f_2}$ and $h_{f_1f_2}$, then we calculate the Laplacian$-$Beltrami operator of a function on a generalized warped product which may be expressed in terms of the local restrictions of the functions to the base and fiber. Finally, we conclude some interesting relationships between the geometry of the couples $(M_1,g_1)$ and $(M_2,g_2)$ and that of $(M_1\times M_2,h_{f_1f_2})$.


title: Formulation of the relativistic heat equation and the relativistic kinetic Fokker-Planck equation using GENERIC

abstract: In this paper, we formulate the relativistic heat equation and the relativistic kinetic Fokker-Planck equations into the GENERIC (General Equation for Non-Equilibrium Reversible-Irreversible Coupling) framework. We also show that the relativistic Maxwellian distribution is the stationary solution of the latter. The GENERIC formulation provides an alternative justification that the two equations are meaningful relativistic generalizations of their nonrelativistic counterparts.


title: Alignment of metabolic trajectories with application to metabonomic toxicology

abstract: Geometry of the metabolic trajectories is characteristic of the biological response (Keun, Ebbels et al. 2004). Yet, due to unavoidable inter-individual variations, the exact trajectories characterising the biological responses differ. We examined whether the differences seen between metabolic trajectories of a specific treatment, correspond to the variations seen in the other biological manifestations of the same treatment. Differences in trajectories were measured via alignment procedures which introduced and implemented in this study. Our study revealed strong correlation between the scales of the aligned trajectories of metabolic responses and the severity of the hepatocelluar lesions induced after administration of hydrazine. Thus the results confirm that aligned trajectories are characteristic of a specific treatment. They then can be used for comparison with other treatment specific or unknown metabolic trajectories and can have many metabonomic applications such as preclinical toxicological screening


title: QANUS: An Open-source Question-Answering Platform

abstract: In this paper, we motivate the need for a publicly available, generic software framework for question-answering (QA) systems. We present an open-source QA framework QANUS which researchers can leverage on to build new QA systems easily and rapidly. The framework implements much of the code that will otherwise have been repeated across different QA systems. To demonstrate the utility and practicality of the framework, we further present a fully functioning factoid QA system QA-SYS built on top of QANUS.


title: Statistical consistency and asymptotic normality for high-dimensional robust M-estimators

abstract: We study theoretical properties of regularized robust M-estimators, applicable when data are drawn from a sparse high-dimensional linear model and contaminated by heavy-tailed distributions and/or outliers in the additive errors and covariates. We first establish a form of local statistical consistency for the penalized regression estimators under fairly mild conditions on the error distribution: When the derivative of the loss function is bounded and satisfies a local restricted curvature condition, all stationary points within a constant radius of the true regression vector converge at the minimax rate enjoyed by the Lasso with sub-Gaussian errors. When an appropriate nonconvex regularizer is used in place of an l_1-penalty, we show that such stationary points are in fact unique and equal to the local oracle solution with the correct support---hence, results on asymptotic normality in the low-dimensional case carry over immediately to the high-dimensional setting. This has important implications for the efficiency of regularized nonconvex M-estimators when the errors are heavy-tailed. Our analysis of the local curvature of the loss function also has useful consequences for optimization when the robust regression function and/or regularizer is nonconvex and the objective function possesses stationary points outside the local region. We show that as long as a composite gradient descent algorithm is initialized within a constant radius of the true regression vector, successive iterates will converge at a linear rate to a stationary point within the local region. Furthermore, the global optimum of a convex regularized robust regression function may be used to obtain a suitable initialization. The result is a novel two-step procedure that uses a convex M-estimator to achieve consistency and a nonconvex M-estimator to increase efficiency.


title: Dynamic rotor mode in antiferromagnetic nanoparticles

abstract: We present experimental, numerical, and theoretical evidence for a new mode of antiferromagnetic dynamics in nanoparticles. Elastic neutron scattering experiments on 8 nm particles of hematite display a loss of diffraction intensity with temperature, the intensity vanishing around 150 K. However, the signal from inelastic neutron scattering remains above that temperature, indicating a magnetic system in constant motion. In addition, the precession frequency of the inelastic magnetic signal shows an increase above 100 K. Numerical Langevin simulations of spin dynamics reproduce all measured neutron data and reveal that thermally activated spin canting gives rise to a new type of coherent magnetic precession mode. This "rotor" mode can be seen as a high-temperature version of superparamagnetism and is driven by exchange interactions between the two magnetic sublattices. The frequency of the rotor mode behaves in fair agreement with a simple analytical model, based on a high temperature approximation of the generally accepted Hamiltonian of the system. The extracted model parameters, as the magnetic interaction and the axial anisotropy, are in excellent agreement with results from Mossbauer spectroscopy.


title: Inclusive semileptonic B decays and |Vcb|

abstract: The magnitude of the CKM matrix elements Vcb can be extracted from inclusive semileptonic B decays in a model independent way pioneered by Kolya Uraltsev and collaborators. I review here the present status and latest developments in this field.


title: Filaments of the radio cosmic web: opportunities and challenges for SKA

abstract: The detection of the diffuse gas component of the cosmic web remains a formidable challenge. In this work we study synchrotron emission from the cosmic web with simulated SKA1 observations, which can represent an fundamental probe of the warm-hot intergalactic medium. We investigate radio emission originated by relativistic electrons accelerated by shocks surrounding cosmic filaments, assuming diffusive shock acceleration and as a function of the (unknown) large-scale magnetic fields. The detection of the brightest parts of large ($>10 \rm Mpc$) filaments of the cosmic web should be within reach of the SKA1-LOW, if the magnetic field is at the level of a $\sim 10$ percent equipartition with the thermal gas, corresponding to $\sim 0.1 \mu G$ for the most massive filaments in simulations. In the course of a 2-years survey with SKA1-LOW, this will enable a first detection of the "tip of the iceberg" of the radio cosmic web, and allow for the use of the SKA as a powerful tool to study the origin of cosmic magnetism in large-scale structures. On the other hand, the SKA1-MID and SKA1-SUR seem less suited for this science case at low redshift ($z \leq 0.4$), owing to the missing short baselines and the consequent lack of signal from the large-scale brightness fluctuations associated with the filaments. In this case only very long exposures ($\sim 1000$ hr) may enable the detection of $\sim 1-2$ filament for field of view in the SKA1-SUR PAF Band1.


title: Phenomenological modelling for Time-Resolved Electron Paramagnetic Resonance in radical-triplet system

abstract: The spin dynamics of radical-triplet system (RTS) has been calculated by using the Lindblad formalism within the theory of open quantum system. The single-radical-triplet system (SRTS) is considered here for single-qubit quantum gate operations while double-radical-triplet system (DRTS) for two-qubit operations. The environment effects taken into account include the spin-lattice relaxation of the triplet exciton and radical spin-$\frac{1}{2}$, the inter-system crossing process that induces the transition from singlet excited state to the triplet ground state, and the rather slow relaxation process from the triplet ground state back down to the singlet ground state. These calculations shown that the line shape broadening is strongly related to the exchange interaction between triplet and exciton, which can be understood as a spontaneous magnetic field created by the triplet renormalises the original spin-$\frac{1}{2}$ electron spin resonance spectra. This work will provide key information about the spin dynamics for building optically-controlled molecular quantum gate out of radical-bearing molecules. Moreover, this has generated the further theoretical question on how the mixture of fermion and boson behaves.


title: The Influence of Mergers on Scatter and Evolution in Sunyaev-Zel'dovich Effect Scaling Relations

abstract: The Sunyaev-Zel'dovich effect (SZE) observable-mass (Y-M) scaling relation is a promising technique for obtaining mass estimates for large samples of galaxy clusters and holds a key to studying the nature of dark matter and dark energy. However, cosmological inference based on SZE cluster surveys is limited by our incomplete knowledge of bias, scatter, and evolution in the Y-M relation. In this work, we investigate the effects of galaxy cluster mergers on the scaling relation using the Omega500 high-resolution cosmological hydrodynamic simulation. We show that the non-thermal pressure associated with merger-induced gas motions contributes significantly to the bias, scatter, and evolution of the scaling relation. After the merger, the kinetic energy of merging systems is slowly converted into thermal energy through dissipation of turbulent gas motions, which causes the thermal SZE signal to increase steadily with time. This post-merger evolution is one of the primary source of bias and scatter in the Y-M relation. However, we show that when the missing non-thermal energy is accounted for, the resulting relation exhibits little to no redshift evolution and the scatter around the scaling relation is ~20-30 % smaller than that of the thermal SZE signal alone. Our work opens up a possibility to further improve the current robust mass proxy, Y, by accounting for the missing non-thermal pressure component. We discuss future prospect of measuring internal gas motions in galaxy clusters and its implication for cluster-based cosmological tests.


title: Contact Representations of Sparse Planar Graphs

abstract: We study representations of graphs by contacts of circular arcs, CCA-representations for short, where the vertices are interior-disjoint circular arcs in the plane and each edge is realized by an endpoint of one arc touching the interior of another. A graph is (2,k)-sparse if every s-vertex subgraph has at most 2s - k edges, and (2, k)-tight if in addition it has exactly 2n - k edges, where n is the number of vertices. Every graph with a CCA- representation is planar and (2, 0)-sparse, and it follows from known results on contacts of line segments that for k >= 3 every (2, k)-sparse graph has a CCA-representation. Hence the question of CCA-representability is open for (2, k)-sparse graphs with 0 <= k <= 2. We partially answer this question by computing CCA-representations for several subclasses of planar (2,0)-sparse graphs. In particular, we show that every plane (2, 2)-sparse graph has a CCA-representation, and that any plane (2, 1)-tight graph or (2, 0)-tight graph dual to a (2, 3)-tight graph or (2, 4)-tight graph has a CCA-representation. Next, we study CCA-representations in which each arc has an empty convex hull. We characterize the plane graphs that have such a representation, based on the existence of a special orientation of the graph edges. Using this characterization, we show that every plane graph of maximum degree 4 has such a representation, but that finding such a representation for a plane (2, 0)-tight graph with maximum degree 5 is an NP-complete problem. Finally, we describe a simple algorithm for representing plane (2, 0)-sparse graphs with wedges, where each vertex is represented with a sequence of two circular arcs (straight-line segments).


title: Electronic structure and magnetic properties of $\mathrm{CaCr}\mathrm{O}_3$: The interplay between spin- and orbital-orderings

abstract: The electronic structure and magnetic properties of CaCr$\mathrm{O}_3$ have been calculated by two methods, including hybrid-exchange density-function theory and density-functional theory + $U$. The computed densities of states from both of these methods are in a qualitative agreement with the previous x-ray spectroscopy. On the other hand, the opening of the band gap separates them apart. hybrid-exchange density-functional theory always gives a finite band gap, down to $\sim 1.2$ eV from HSE06 functional, whereas by tuning the Hubbard-$U$ parameter down to $0.5$ eV, a conducting state with AFM-C (defined in the text) spin configuration can be achieved. From hybrid density-functional theory, the computed nearest-neighbouring exchange interaction along the $c$-axis and in the $ab$-plane are $\sim 4$ meV and $\sim 6$ meV (anti-ferromagnetic), respectively, which are qualitatively in agreement with the previous magnetic measurements. These anti-ferromagnetic exchange interaction, together with the in-plane anti-ferro-orbital ordering will induce a spin-orbital frustration, which could play a role for the abnormal electronic properties in CaCrO$_3$. In hybrid-exchange density-functional theory, an abrupt reduction ($\sim 0.2$ eV) of the majority-spin band gap of the ferromagnetic state between 60 K and 100 K has been found as lowering temperature, which shows a strong link to the previous optical conductivity measurements in [A. C. Komarek, et. al., Phys. Rev. B \textbf{84}, 125114 (2011)]. In sharp contrast, the density-functional theory + $U$ methods predicted AFM-C state as the lowest AFM state for the crystal structure measured below 90 K, above which AFM-A is however the lowest. The closely related concepts including electron-hole liquid and surface-plasmon-mediating spin-spin interactions have been discussed as well.


title: A robust sub-linear time R-FFAST algorithm for computing a sparse DFT

abstract: The Fast Fourier Transform (FFT) is the most efficiently known way to compute the Discrete Fourier Transform (DFT) of an arbitrary n-length signal, and has a computational complexity of O(n log n). If the DFT X of the signal x has only k non-zero coefficients (where k < n), can we do better? In [1], we addressed this question and presented a novel FFAST (Fast Fourier Aliasing-based Sparse Transform) algorithm that cleverly induces sparse graph alias codes in the DFT domain, via a Chinese-Remainder-Theorem (CRT)-guided sub-sampling operation of the time-domain samples. The resulting sparse graph alias codes are then exploited to devise a fast and iterative onion-peeling style decoder that computes an n length DFT of a signal using only O(k) time-domain samples and O(klog k) computations. The FFAST algorithm is applicable whenever k is sub-linear in n (i.e. k = o(n)), but is obviously most attractive when k is much smaller than n. In this paper, we adapt the FFAST framework of [1] to the case where the time-domain samples are corrupted by a white Gaussian noise. In particular, we show that the extended noise robust algorithm R-FFAST computes an n-length k-sparse DFT X using O(klog ^3 n) noise-corrupted time-domain samples, in O(klog^4n) computations, i.e., sub-linear time complexity. While our theoretical results are for signals with a uniformly random support of the non-zero DFT coefficients and additive white Gaussian noise, we provide simulation results which demonstrates that the R-FFAST algorithm performs well even for signals like MR images, that have an approximately sparse Fourier spectrum with a non-uniform support for the dominant DFT coefficients.


title: Unravelling the origin of large-scale magnetic fields in galaxy clusters and beyond through Faraday Rotation Measures with the SKA

abstract: We investigate the possibility for the SKA to detect and study the magnetic fields in galaxy clusters and in the less dense environments surrounding them using Faraday Rotation Measures. To this end, we produce 3-dimensional magnetic field models for galaxy clusters of different masses and in different stages of their evolution, and derive mock rotation measure observations of background radiogalaxies. According to our results, already in phase I, we will be able to infer the magnetic field properties in galaxy clusters as a function of the cluster mass, down to $10^{13}$ solar-masses. Moreover, using cosmological simulations to model the gas density, we have computed the expected rotation measure through shock-fronts that occur in the intra-cluster medium during cluster mergers. The enhancement in the rotation measure due to the density jump will permit to constraint the magnetic field strength and structure after the shock passage. SKA observations of polarised sources located behind galaxy clusters will answer several questions about the magnetic field strength and structure in galaxy clusters, and its evolution with cosmic time.


title: A hybrid-exchange density-functional theory study of the electronic structure of $\mathrm{MnV}_2\mathrm{O}_4$: Exotic orbital ordering in the cubic structure

abstract: The electronic structures of the cubic and tetragonal $\mathrm{MnV}_2\mathrm{O}_4$ have been studied by using hybrid-exchange density functional theory. The computed electronic structure of the tetragonal phase shows an anti-ferro orbital ordering on V sites and a ferrimagnetic ground state (the spins on V and Mn are anti-aligned). These results are in a good agreement with the previous theoretical result obtained from the local-density approximation+$U$ methods [S. Sarkar, et. al., Phys. Rev. Lett. 102, 216405 (2009)]. Moreover, the electronic structure, especially the projected density of states of the cubic phase has been predicted with a good agreement with the recent soft x-ray spectroscopy experiment. Similar to the tetragonal phase, the spins on V and Mn in the cubic structure favour a ferrimagnetic configuration. Most interesting is that the computed charge densities of the spin-carrying orbitals on V in the cubic phase show an exotic orbital ordering, i.e., a ferro-orbital ordering along [110] but an anti-ferro-orbital ordering along [$\overline{1}$10].


title: A Semi-linear Energy Critical Wave Equation With Applications

abstract: In this work we consider a semi-linear energy critical wave equation in ${\mathbb R}^d$ ($3\leq d \leq 5$) \[ \partial_t^2 u - \Delta u = \pm \phi(x) |u|^{4/(d-2)} u, \qquad (x,t)\in {\mathbb R}^d \times {\mathbb R} \] with initial data $(u, \partial_t u)|_{t=0} = (u_0,u_1) \in \dot{H}^1 \times L^2 ({\mathbb R}^d)$. Here the function $\phi \in C({\mathbb R}^d; (0,1])$ converges to zero as $|x| \rightarrow \infty$. We follow the same compactness-rigidity argument as Kenig and Merle applied on the Cauchy problem of the equation \[ \partial_t^2 u - \Delta u = |u|^{4/(d-2)} u \] and obtain a similar result when $\phi$ satisfies some technical conditions. In the defocusing case we prove that the solution scatters for any initial data in the energy space $\dot{H}^1 \times L^2$. While in the focusing case we can determine the global behaviour of the solutions, either scattering or finite-time blow-up, according to their initial data when the energy is smaller than a certain threshold.


title: A New Sparse Matrix Vector Multiplication GPU Algorithm Designed for Finite Element Problems

abstract: Recently, graphics processors (GPUs) have been increasingly leveraged in a variety of scientific computing applications. However, architectural differences between CPUs and GPUs necessitate the development of algorithms that take advantage of GPU hardware. As sparse matrix vector multiplication (SPMV) operations are commonly used in finite element analysis, a new SPMV algorithm and several variations are developed for unstructured finite element meshes on GPUs. The effective bandwidth of current GPU algorithms and the newly proposed algorithms are measured and analyzed for 15 sparse matrices of varying sizes and varying sparsity structures. The effects of optimization and differences between the new GPU algorithm and its variants are then subsequently studied. Lastly, both new and current SPMV GPU algorithms are utilized in the GPU CG Solver in GPU finite element simulations of the heart. These results are then compared against parallel PETSc finite element implementation results. The effective bandwidth tests indicate that the new algorithms compare very favorably with current algorithms for a wide variety of sparse matrices and can yield very notable benefits. GPU finite element simulation results demonstrate the benefit of using GPUs for finite element analysis, and also show that the proposed algorithms can yield speedup factors up to 12-fold for real finite element applications.


title: The Host Galaxies and Narrow Line Regions of Four Double-Peaked [OIII] AGN

abstract: Major gas-rich mergers of galaxies are expected to play an important role in triggering and fuelling luminous AGN. We present deep multi-band (u/r/z) imaging and long slit spectroscopy of four double-peaked [OIII] emitting AGN, a class of objects associated with either kcp-separated binary AGN or final stage major mergers, though AGN with complex narrow-line regions are known contaminants. Such objects are of interest since they represent the onset of AGN activity during the merger process. Three of the objects studied have been confirmed as major mergers using near-infrared imaging, one is a confirmed X-ray binary AGN. All AGN are luminous and have redshifts of 0.1 < z < 0.4. Deep r-band images show that a majority (3/4) of the sources have disturbed host morphologies and tidal features, while the remaining source is morphologically undisturbed down to low surface brightness limits. The lack of morphological disturbances in this galaxy despite the fact that is is a close binary AGN suggests that the merger of a binary black hole can take longer than ~1 Gyr. The narrow line regions (NLRs) have large sizes (10 kpc < r < 100 kpc) and consist of compact clumps with considerable relative velocities (~ 200-650 km/s). We detect broad, predominantly blue, wings with velocities up to ~1500 km/s in [OIII], indicative of powerful outflows. The outflows are compact (<5 kpc) and co-spatial with nuclear regions showing considerable reddening, consistent with enhanced star formation. One source shows an offset between gas and stellar kinematics, consistent with either a bipolar flow or a counter-rotating gas disk. We are not able to unambiguously identify the sources as binary AGN using our data, X-ray or radio data is required for an unambiguous identification. However, the data still yield interesting results for merger triggering of AGN and time-scales of binary black hole mergers.


title: Decomposition approaches to integration without a measure

abstract: Extending the idea of Even and Lehrer [3], we discuss a general approach to integration based on a given decomposition system equipped with a weighting function, and a decomposition of the integrated function. We distinguish two type of decompositions: sub-decomposition based integrals (in economics linked with optimization problems to maximize the possible profit) and super-decomposition based integrals (linked with costs minimization). We provide several examples (both theoretical and realistic) to stress that our approach generalizes that of Even and Lehrer [3] and also covers problems of linear programming and combinatorial optimization. Finally, we introduce some new types of integrals related to optimization tasks.