title: Hypersurfaces quartiques de dimension 3 : non rationalit\'e stable

abstract: Inspir\'es par un argument de C. Voisin, nous montrons l'existence d'hypersurfaces quartiques lisses dans ${\bf P}^4_{\mathbb C}$ qui ne sont pas stablement rationnelles, plus pr\'ecis\'ement dont le groupe de Chow de degr\'e z\'ero n'est pas universellement \'egal \`a $\mathbb Z$. --- There are (many) smooth quartic hypersurfaces in ${\bf P}^4_{\mathbb C}$ which are not stably rational. More precisely, their degree zero Chow group is not universally equal to $\mathbb Z$. The proof uses a variation of a specialisation method due to C. Voisin.


title: The Giant Gemini GMOS survey of z>4.4 quasars - I. Measuring the mean free path across cosmic time

abstract: We have obtained spectra of 163 quasars at $z_\mathrm{em}>4.4$ with the Gemini Multi Object Spectrometers on the Gemini North and South telescopes, the largest publicly available sample of high-quality, low-resolution spectra at these redshifts. From this homogeneous data set, we generated stacked quasar spectra in three redshift intervals at $z\sim 5$. We have modelled the flux below the rest-frame Lyman limit ($\lambda_\mathrm{r}<912$\AA) to assess the mean free path $\lambda_\mathrm{mfp}^{912}$ of the intergalactic medium to HI-ionizing radiation. At mean redshifts $z_\mathrm{q}=4.56$, 4.86 and 5.16, we measure $\lambda_\mathrm{mfp}^{912}=(22.2\pm 2.3, 15.1\pm 1.8, 10.3\pm 1.6)h_{70}^{-1}$ proper Mpc with uncertainties dominated by sample variance. Combining our results with $\lambda_\mathrm{mfp}^{912}$ measurements from lower redshifts, the data are well modelled by a simple power-law $\lambda_\mathrm{mfp}^{912}=A[(1+z)/5]^\eta$ with $A=(37\pm 2)h_{70}^{-1}$ Mpc and $\eta = -5.4\pm 0.4$ between $z=2.3$ and $z=5.5$. This rapid evolution requires a physical mechanism -- beyond cosmological expansion -- which reduces the cosmic effective Lyman limit opacity. We speculate that the majority of HI Lyman limit opacity manifests in gas outside galactic dark matter haloes, tracing large-scale structures (e.g. filaments) whose average density (and consequently neutral fraction) decreases with cosmic time. Our measurements of the strongly redshift-dependent mean free path shortly after the completion of HI reionization serve as a valuable boundary condition for numerical models thereof. Having measured $\lambda_\mathrm{mfp}^{912}\approx 10$ Mpc at $z=5.2$, we confirm that the intergalactic medium is highly ionized by that epoch and that the redshift evolution of the mean free path does not show a break that would indicate a recent end to HI reionization.


title: Getting a Grip on the Transverse Motion in a Zeeman Decelerator

abstract: Zeeman deceleration is an experimental technique in which inhomogeneous, time-dependent magnetic fields generated inside an array of solenoid coils are used to manipulate the velocity of a supersonic beam. A 12-stage Zeeman decelerator has been built and characterized using hydrogen atoms as a test system. The instrument has several original features including the possibility to replace each deceleration coil individually. In this article, we give a detailed description of the experimental setup, and illustrate its performance. We demonstrate that the overall acceptance in a Zeeman decelerator can be significantly increased with only minor changes to the setup itself. This is achieved by applying a rather low, anti-parallel magnetic field in one of the solenoid coils that forms a temporally varying quadrupole field, and improves particle confinement in the transverse direction. The results are reproduced by three-dimensional numerical particle trajectory simulations thus allowing for a rigorous analysis of the experimental data. The findings suggest the use of a modified coil configuration to improve transverse focusing during the deceleration process.


title: Algebraic vs. topological vector bundles on spheres

abstract: We study the problem of when a topological vector bundle on a smooth complex affine variety admits an algebraic structure. We prove that all rank $2$ topological complex vector bundles on smooth affine quadrics of dimension $11$ over the complex numbers admit algebraic structures.


title: Conservative collision prediction and avoidance for stochastic trajectories in continuous time and space

abstract: Existing work in multi-agent collision prediction and avoidance typically assumes discrete-time trajectories with Gaussian uncertainty or that are completely deterministic. We propose an approach that allows detection of collisions even between continuous, stochastic trajectories with the only restriction that means and variances can be computed. To this end, we employ probabilistic bounds to derive criterion functions whose negative sign provably is indicative of probable collisions. For criterion functions that are Lipschitz, an algorithm is provided to rapidly find negative values or prove their absence. We propose an iterative policy-search approach that avoids prior discretisations and yields collision-free trajectories with adjustably high certainty. We test our method with both fixed-priority and auction-based protocols for coordinating the iterative planning process. Results are provided in collision-avoidance simulations of feedback controlled plants.


title: Runge-Kutta Central Discontinuous Galerkin BGK Method for the Navier-Stokes Equations

abstract: In this paper, we propose a Runge-Kutta (RK) central discontinuous Galerkin (CDG) gas-kinetic BGK method for the Navier-Stokes equations. The proposed method is based on the CDG method defined on two sets of overlapping meshes to avoid discontinuous solutions at cell interfaces, as well as the gas-kinetic BGK model to evaluate fluxes for both convection and diffusion terms. Redundant representation of the numerical solution in the CDG method offers great convenience in the design of gas-kinetic BGK fluxes. Specifically, the evaluation of fluxes at cell interfaces of one set of computational mesh is right inside the cells of the staggered mesh, hence the corresponding particle distribution function for flux evaluation is much simpler than that in existing gas-kinetic BGK methods. As a central scheme, the proposed CDG-BGK has doubled the memory requirement as the corresponding DG scheme; on the other hand, {for the convection part,} the CFL time step constraint of the CDG method for numerical stability is relatively large compared with that for the DG method. Numerical boundary conditions have to be treated with special care. Numerical examples for 1D and 2D viscous flow simulations are presented to validate the accuracy and robustness of the proposed RK CDG-BGK method.


title: Application of Pseudo-Transient Continuation Method in Dynamic Stability Analysis

abstract: In this paper, pseudo-transient continuation method has been modified and implemented in power system long-term stability analysis. This method is a middle ground between integration and steady state calculation, thus is a good compromise between accuracy and efficiency. Pseudo-transient continuation method can be applied in the long-term stability model directly to accelerate simulation speed and can also be implemented in the QSS model to overcome numerical difficulties. Numerical examples show that pseudo-transient continuation method can provide correct approximations for the long-term stability model in terms of trajectories and stability assessment.


title: Maximizing the Signal-to-Alias Ratio in Non-Uniform Filter Banks for Acoustic Echo Cancellation

abstract: A new method for designing non-uniform filter-banks for acoustic echo cancellation is proposed. In the method, the analysis prototype filter design is framed as a convex optimization problem that maximizes the signal-to-alias ratio (SAR) in the analysis banks. Since each sub-band has a different bandwidth, the contribution to the overall SAR from each analysis bank is taken into account during optimization. To increase the degrees of freedom during optimization, no constraints are imposed on the phase or group delay of the filters; at the same time, low delay is achieved by ensuring that the resulting filters are minimum phase. Experimental results show that the filter bank designed using the proposed method results in a sub-band adaptive filter with a much better echo return loss enhancement (ERLE) when compared with existing design methods.


title: Collapsing rotating shells in Myers-Perry-AdS$_5$ spacetime: A perturbative approach

abstract: We study gravitational perturbations sourced by a rotating test shell collapsing into five-dimensional Myers-Perry black holes in anti-de Sitter (AdS). Our attention is restricted to the case in which the two possible angular momenta of the geometry are set equal. In this situation the background is cohomogeneity-1, which amounts to a crucial technical simplification. It is found that the linearized Einstein equations are consistent only when the test shell is corotating with the spacetime. However, it is argued that this is a consequence of the matter on the shell being described by dust or, more precisely, noninteracting test particles. We compute the mass and angular momenta of the perturbed spacetime using a counterterm subtraction method, for which we provide an explicit formula that has not appeared previously in the literature. The results are in agreement with the expected expressions for energy and angular momenta of geodesic particles in AdS$_5$.


title: An Analytic Model for Buoyancy Resonances in Protoplanetary Disks

abstract: Zhu, Stone, and Rafikov (2012) found in 3D shearing box simulations a new form of planet-disk interaction that they attributed to a vertical buoyancy resonance in the disk. We describe an analytic linear model for this interaction. We adopt a simplified model involving azimuthal forcing that produces the resonance and permits an analytic description of its structure. We derive an analytic expression for the buoyancy torque and show that the vertical torque distribution agrees well with results of Athena simulations and a Fourier method for linear numerical calculations carried out with the same forcing. The buoyancy resonance differs from the classic Lindblad and corotation resonances in that the resonance lies along tilted planes. Its width depends on damping effects and is independent of the gas sound speed. The resonance does not excite propagating waves. At a given large azimuthal wavenumber k_y > 1/h (for disk thickness h), the buoyancy resonance exerts a torque over a region that lies radially closer to the corotation radius than the Lindblad resonance. Because the torque is localized to the region of excitation, it is potentially subject to the effects of nonlinear saturation. In addition, the torque can be reduced by the effects of radiative heat transfer between the resonant region and its surroundings. For each azimuthal wavenumber, the resonance establishes a large scale density wave pattern in a plane within the disk.


title: General random walk in a random environment defined on Galton-Watson trees

abstract: We consider the motion of a particle on a Galton Watson tree, when the probabilities of jumping from a vertex to any one of its neighbours is determined by a random process. Given the tree, positive weights are assigned to the edges in such a way that, viewed along any line of descent, they evolve as a random process. In order to introduce our method for proving transience or recurrence, we first suppose that the weights are i.i.d., reproving a result of Lyons and Pemantle. We then extend the argument to allow a Markovian environment, and finally to a random walk on a Markovian environment that changes the environment. Our approach involves studying the typical behaviour of processes on fixed lines of descent, which we then show determines the behaviour of the process on the whole tree.


title: Fine-grained Patches for Java Software Upgrades

abstract: We present a novel methodology for deriving fine-grained patches of Java software. We consider an abstract-syntax tree (AST) representation of Java classes compiled to the Java Virtual Machine (JVM) format, and a difference analysis over the AST representation to derive patches. The AST representation defines an appropriate abstraction level for analyzing differences, yielding compact patches that correlate modularly to actual source code changes. The approach contrasts to other common, coarse-grained approaches, like plain binary differences, which may easily lead to disproportionately large patches. We present the main traits of the methodology, a prototype tool called aspa that implements it, and a case-study analysis on the use of aspa to derive patches for the Java 2 SE API. The case-study results illustrate that aspa patches have a significantly smaller size than patches derived by binary differencing tools.


title: Existence of solutions for a system of coupled nonlinear stationary bi-harmonic Schr\"{o}dinger equations

abstract: We obtain existence and multiplicity results for the solutions of a class of coupled semilinear bi-harmonic Schr\"{o}dinger equations. Actually, using the classical Mountain Pass Theorem and minimization techniques, we prove the existence of critical points of the associated functional constrained on the \emph{Nehari manifold}. Furthermore, we show that using the so-called \emph{fibering method} and the \emph{Lusternik-Schnirel'man theory} there exist infinitely many solutions, actually a countable family of critical points, for such a semiliner bi-hamonic Schr\"{o}dinger system under study in this work.


title: The Dynamic Age of Centaurus A

abstract: In this paper I present dynamic models of the radio source Centaurus A, and critique possible models of in situ particle reacceleration (ISR) within the radio lobes. The radio and gamma-ray data require neither homogeneous plasma nor quasi-equipartition between plasma and magnetic field; inhomogeneous models containing both high-field and low-field regions are equally likely. Cen A cannot be as young as the radiative lifetimes of its relativistic electrons, which range from a few to several tens of Myr. Two classes of dynamic models -- flow driven and magnetically driven -- are consistent with current observations; each requires Cen A to be on the order of a Gyr old. Thus, ongoing ISR must be occurring within the radio source. Alfven-wave ISR is probably occurring throughout the source, and may be responsible for maintaining the gamma-ray-loud electrons. It is likely to be supplemented by shock or reconnection ISR which maintains the radio-loud electrons in high-field regions.


title: Observation of suppression of light scattering induced by dipole-dipole interactions in a cold atomic ensemble

abstract: We study the emergence of collective scattering in the presence of dipole-dipole interactions when we illuminate a cold cloud of rubidium atoms with a near-resonant and weak intensity laser. The size of the atomic sample is comparable to the wavelength of light. When we gradually increase the atom number from 1 to 450, we observe a broadening of the line, a small red shift and, consistently with these, a strong suppression of the scattered light with respect to the noninteracting atom case. Numerical simulations, which include the internal atomic level structure, agree with the data.


title: On Topological Structure of the First Non-abelian Cohomology of Topological Groups

abstract: Let $G$, $R$ and $A$ be topological groups. Suppose that $G$ and $R$ act continuously on $A$, and $G$ acts continuously on $R$. In this paper, we define a partially crossed topological $G-R$-bimodule $(A,\mu)$, where $\mu:A\rightarrow R$ is a continuous homomorphism. Let $Der_{c}(G,(A,\mu))$ be the set of all $(\alpha,r)$ such that $\alpha:G\rightarrow A$ is a continuous crossed homomorphism and $\mu\alpha(g)=r^{g}r^{-1}$. We introduce a topology on $Der_{c}(G,(A,\mu))$. We show that $Der_{c}(G,(A,\mu))$ is a topological group, wherever $G$ and $R$ are locally compact. We define the first cohomology, $H^{1}(G,(A,\mu))$, of $G$ with coefficients in $(A,\mu)$ as a quotient space of $Der_{c}(G,(A,\mu))$. Also, we state conditions under which $H^{1}(G,(A,\mu))$ is a topological group. Finally, we show that under what conditions $H^{1}(G,(A,\mu))$ is one of the following: $k$-space, discrete, locally compact and compact.


title: The looping rate and sandpile density of planar graphs

abstract: We give a simple formula for the looping rate of loop-erased random walk on a finite planar graph. The looping rate is closely related to the expected amount of sand in a recurrent sandpile on the graph. The looping rate formula is well-suited to taking limits where the graph tends to an infinite lattice, and we use it to give an elementary derivation of the (previously computed) looping rate and sandpile densities of the square, triangular, and honeycomb lattices, and compute (for the first time) the looping rate and sandpile densities of many other lattices, such as the kagome lattice, the dice lattice, and the truncated hexagonal lattice (for which the values are all rational), and the square-octagon lattice (for which it is transcendental).


title: A Free-Form Lensing Grid Solution for A1689 with New Mutiple Images

abstract: Hubble Space Telescope imaging of the galaxy cluster Abell 1689 has revealed an exceptional number of strongly lensed multiply-imaged galaxies, including high-redshift candidates. Previous studies have used this data to obtain the most detailed dark matter reconstructions of any galaxy cluster to date, resolving substructures ~25 kpc across. We examine Abell 1689 (hereafter, A1689) non-parametrically, combining strongly lensed images and weak distortions from wider field Subaru imaging, and we incorporate member galaxies to improve the lens solution. Strongly lensed galaxies are often locally affected by member galaxies, however, these perturbations cannot be recovered in grid based reconstructions because the lensing information is too sparse to resolve member galaxies. By adding luminosity-scaled member galaxy deflections to our smooth grid we can derive meaningful solutions with sufficient accuracy to permit the identification of our own strongly lensed images, so our model becomes self consistent. We identify 11 new multiply lensed system candidates and clarify previously ambiguous cases, in the deepest optical and NIR data to date from Hubble and Subaru. Our improved spatial resolution brings up new features not seen when the weak and strong lensing effects are used separately, including clumps and filamentary dark matter around the main halo. Our treatment means we can obtain an objective mass ratio between the cluster and galaxy components, for examining the extent of tidal stripping of the luminous member galaxies. We find a typical mass-to-light ratios of M/L_B = 21 inside the r<1 arcminute region that drops to M/L_B = 17 inside the r<40 arcsecond region. Our model independence means we can objectively evaluate the competitiveness of stacking cluster lenses for defining the geometric lensing-distance-redshift relation in a model independent way.


title: Reconstructing the world trade multiplex: the role of intensive and extensive biases

abstract: In economic and financial networks, the strength of each node has always an important economic meaning, such as the size of supply and demand, import and export, or financial exposure. Constructing null models of networks matching the observed strengths of all nodes is crucial in order to either detect interesting deviations of an empirical network from economically meaningful benchmarks or reconstruct the most likely structure of an economic network when the latter is unknown. However, several studies have proved that real economic networks and multiplexes are topologically very different from configurations inferred only from node strengths. Here we provide a detailed analysis of the World Trade Multiplex by comparing it to an enhanced null model that simultaneously reproduces the strength and the degree of each node. We study several temporal snapshots and almost one hundred layers (commodity classes) of the multiplex and find that the observed properties are systematically well reproduced by our model. Our formalism allows us to introduce the (static) concept of extensive and intensive bias, defined as a measurable tendency of the network to prefer either the formation of extra links or the reinforcement of link weights, with respect to a reference case where only strengths are enforced. Our findings complement the existing economic literature on (dynamic) intensive and extensive trade margins. More in general, they show that real-world multiplexes can be strongly shaped by layer-specific local constraints.


title: A propositional system induced by Japaridze's approach to IF logic

abstract: Cirquent calculus is a new proof-theoretic and semantic approach introduced for the needs of computability logic by G.Japaridze, who also showed that, through cirquent calculus, one can capture, refine and generalize independence-friendly (IF) logic. Specifically, the approach allows us to account for independence from propositional connectives in the same spirit as the traditional IF logic accounts for independence from quantifiers. Japaridze's treatment of IF logic, however, was purely semantical, and no deductive system was proposed. The present paper constructs a formal system sound and complete w.r.t. the propositional fragment of Japaridze's cirquent-based semantics for IF logic. Such a system can thus be considered an axiomatization of purely propositional IF logic in its full generality.


title: Nonlocal $s$-minimal surfaces and Lawson cones

abstract: The nonlocal $s$-fractional minimal surface equation for $\Sigma= \partial E$ where $E$ is an open set in $R^N$ is given by $$ H_\Sigma^ s (p) := \int_{R^N} \frac {\chi_E(x) - \chi_{E^c}(x)} {|x-p|^{N+s}}\, dx \ =\ 0 \quad \text{for all } p\in \Sigma. $$ Here $0<s<1$, $\chi$ designates characteristic function, and the integral is understood in the principal value sense. The classical notion of minimal surface is recovered by letting $s\to 1$. In this paper we exhibit the first concrete examples (beyond the plane) of nonlocal $s-$minimal surfaces. When $s$ is close to $1$, we first construct a connected embedded $s$-minimal surface of revolution in $R^3$, the {\bf nonlocal catenoid}, an analog of the standard catenoid $|x_3| = \log (r + \sqrt{r^2 -1})$. Rather than eventual logarithmic growth, this surface becomes asymptotic to the cone $|x_3|= r\sqrt{1-s}$. We also find a two-sheet embedded $s$-minimal surface asymptotic to the same cone, an analog to the simple union of two parallel planes. On the other hand, for any $0<s<1$, $n,m\ge 1$, $s-$minimal Lawson cones $|v|=\alpha|u|$, $(u,v)\in R^n\times R^m$, are found to exist. In sharp contrast with the classical case, we prove their stability for small $s$ and $n+m=7$, which suggests that unlike the classical theory (or the case $s$ close to 1), the regularity of $s$-area minimizing surfaces may not hold true in dimension $7$.


title: Chiral Langevin theory for non-Abelian plasmas

abstract: Charged plasmas with chirality imbalance are unstable and tend to reduce the imbalance. This chiral plasma instability is, however, not captured in (anomalous) hydrodynamics for high-temperature non-Abelian plasmas. We derive a Langevin-type classical effective theory with anomalous parity-violating effects for non-Abelian plasmas that describes the chiral plasma instability at the magnetic scale. We show that the time scale of the instability is of order $[g^4 T \ln(1/g)]^{-1}$ at weak coupling.


title: Perturbation Theory for Parent Hamiltonians of Matrix Product States

abstract: This article investigates the stability of the ground state subspace of a canonical parent Hamiltonian of a Matrix product state against local perturbations. We prove that the spectral gap of such a Hamiltonian remains stable under weak local perturbations even in the thermodynamic limit, where the entire perturbation might not be bounded. Our discussion is based on preceding work by D.A. Yarotsky that develops a perturbation theory for relatively bounded quantum perturbation of classical Hamiltonians. We exploit a renormalization procedure, which on large scale transforms the parent Hamiltonian of a Matrix product state into a classical Hamiltonian plus some perturbation. We can thus extend D.A. Yarotsky's results to provide a perturbation theory for parent Hamiltonians of Matrix product states and recover some of the findings of independent contributions by S. Michalakis et al. [arXiv:1109.1588] and J. I. Cirac et al. [arXiv:1306.4003].

categories: math-ph cond-mat.stat-mech math.MP quant-ph


title: Some remarks on Hodge symmetry

abstract: I make some remarks on Hodge symmetry, and prove for instance that if $k$ is a perfect field of characteristic $p>0$ and $X/k$ smooth, proper and Hodge-Witt scheme, and Hodge de Rham sequence of $X$ degenerates at $E_1$ and $X$ has torsion-free crystalline cohomology, then Hodge symmetry holds for $X$.


title: Damage processes in thermoviscoelastic materials with damage-dependent thermal expansion coefficients

abstract: In this paper we prove existence of global in time weak solutions for a highly nonlinear PDE system arising in the context of damage phenomena in thermoviscoelastic materials. The main novelty of the present contribution with respect to the ones already present in the literature consists in the possibility of taking into account a damage-dependent thermal expansion coefficient. This term implies the presence of nonlinear couplings in the PDE system, which make the analysis more challenging.


title: A reclaimer scheduling problem arising in coal stockyard management

abstract: We study a number of variants of an abstract scheduling problem inspired by the scheduling of reclaimers in the stockyard of a coal export terminal. We analyze the complexity of each of the variants, providing complexity proofs for some and polynomial algorithms for others. For one, especially interesting variant, we also develop a constant factor approximation algorithm.


title: Network robustness assessed within a dual connectivity perspective

abstract: Network robustness against attacks has been widely studied in fields as diverse as the Internet, power grids and human societies. Typically, in these studies, robustness is assessed only in terms of the connectivity of the nodes unaffected by the attack. Here we put forward the idea that the connectivity of the affected nodes can play a crucial role in properly evaluating the overall network robustness and its future recovery from the attack. Specifically, we propose a dual perspective approach wherein at any instant in the network evolution under attack, two distinct networks are defined: (i) the Active Network (AN) composed of the unaffected nodes and (ii) the Idle Network (IN) composed of the affected nodes. The proposed robustness metric considers both the efficiency of destroying the AN and the efficiency of building-up the IN. We show, via analysis of both prototype networks and real world data, that trade-offs between the efficiency of Active and Idle network dynamics give rise to surprising crossovers and re-ranking of different attack strategies, pointing to significant implications for decision making.


title: Effect of Deck Deterioration on Overall System Behavior, Resilience and Remaining Life of Composite Steel Girder Bridges

abstract: During past few decades, several studies have been conducted to characterize the performance of in-service girder-type bridge superstructures under operating conditions. Few of these efforts have focused on evaluating the actual response of the bridge systems, especially beyond the elastic limit of their behavior, and correlating the impact of damage to the overall system behavior. In practice, most of the in-service bridge superstructures behave elastically under the routine daily traffic; however, existing damage and deteriorating conditions would significantly influence different aspects of the structural performance including reserve capacity, resilience and remaining service-life. The main purpose of this study is to evaluate the response of composite steel girder bridges under the effect of subsurface delamination in the reinforced concrete deck. Commercial finite element computer software, ANSYS, was implemented to perform a nonlinear analysis on a representative single-span simply supported bridge superstructure. The system failure characteristics were captured in the numerical models by incorporating different sources of material non-linearities including cracking/crushing in the concrete and plasticity in steel components. Upon validation, non-linear behavior of the system with both intact and degraded configurations was used to evaluate the impact of integrated damage mechanism on the overall system performance. Reserve capacity of this bridge superstructure was also determined with respect to the nominal element-level design capacity. As vision to the future path, this framework can be implemented to evaluate the performance of other in-service bridges degraded under the effect of different damage scenarios, thus providing a mechanism to determine a measure of capacity, resilience and remaining service-life.


title: Understanding Persistence: A 3D Trap Map of an H2RG Imaging Sensor

abstract: Several theories exist to explain persistence, most of which revolve around the distribution of traps. We aim to simulate persistence and illustrate this complex issue with a 3D trap map. For this experiment, we vary the detector voltage bias to simulate the change in the depletion region that occurs when the detector is exposed to light. This allows us to measure the distribution of traps in the depletion region. This paper will explore the results from this experiment and discuss the implications.


title: Learning Natural Coding Conventions

abstract: Every programmer has a characteristic style, ranging from preferences about identifier naming to preferences about object relationships and design patterns. Coding conventions define a consistent syntactic style, fostering readability and hence maintainability. When collaborating, programmers strive to obey a project's coding conventions. However, one third of reviews of changes contain feedback about coding conventions, indicating that programmers do not always follow them and that project members care deeply about adherence. Unfortunately, programmers are often unaware of coding conventions because inferring them requires a global view, one that aggregates the many local decisions programmers make and identifies emergent consensus on style. We present NATURALIZE, a framework that learns the style of a codebase, and suggests revisions to improve stylistic consistency. NATURALIZE builds on recent work in applying statistical natural language processing to source code. We apply NATURALIZE to suggest natural identifier names and formatting conventions. We present four tools focused on ensuring natural code during development and release management, including code review. NATURALIZE achieves 94% accuracy in its top suggestions for identifier names and can even transfer knowledge about conventions across projects, leveraging a corpus of 10,968 open source projects. We used NATURALIZE to generate 18 patches for 5 open source projects: 14 were accepted.


title: An Improved Algorithm for Fixed-Hub Single Allocation Problem

abstract: This paper discusses the fixed-hub single allocation problem (FHSAP). In this problem, a network consists of hub nodes and terminal nodes. Hubs are fixed and fully connected; each terminal node is connected to a single hub which routes all its traffic. The goal is to minimize the cost of routing the traffic in the network. In this paper, we propose a linear programming (LP)-based rounding algorithm. The algorithm is based on two ideas. First, we modify the LP relaxation formulation introduced in Ernst and Krishnamoorthy (1996, 1999) by incorporating a set of validity constraints. Then, after obtaining a fractional solution to the LP relaxation, we make use of a geometric rounding algorithm to obtain an integral solution. We show that by incorporating the validity constraints, the strengthened LP often provides much tighter upper bounds than the previous methods with a little more computational effort, and the solution obtained often has a much smaller gap with the optimal solution. We also formulate a robust version of the FHSAP and show that it can guard against data uncertainty with little cost.


title: Borel structurability on the 2-shift of a countable group

abstract: We show that for any infinite countable group $G$ and for any free Borel action $G \curvearrowright X$ there exists an equivariant class-bijective Borel map from $X$ to the free part $\mathrm{Free}(2^G)$ of the $2$-shift $G \curvearrowright 2^G$. This implies that any Borel structurability which holds for the equivalence relation generated by $G \curvearrowright \mathrm{Free}(2^G)$ must hold a fortiori for all equivalence relations coming from free Borel actions of $G$. A related consequence is that the Borel chromatic number of $\mathrm{Free}(2^G)$ is the maximum among Borel chromatic numbers of free actions of $G$. This answers a question of Marks. Our construction is flexible and, using an appropriate notion of genericity, we are able to show that in fact the generic $G$-equivariant map to $2^G$ lands in the free part. As a corollary we obtain that for every $\epsilon > 0$, every free pmp action of $G$ has a free factor which admits a $2$-piece generating partition with Shannon entropy less than $\epsilon$. This generalizes a result of Danilenko and Park.


title: Ultrahigh-Speed Rotating Nanoelectromechanical System (NEMS) Devices Assembled from Nanoscale Building Blocks

abstract: The development of rotary nanomotors is crucial for further advancing NEMS technology. In this work, we report innovative design, assembly, and rotation of ordered arrays of nanomotors. The nanomotors were bottom-up assembled from nanoscale building blocks with nanowires as rotors, patterned nanomagnets as bearings, and quadrupole microelectrodes as stators. Arrays of nanomotors can be synchronously rotated with controlled angle, speed (to at least 18,000 rpm), and chirality by electric fields. Using analytical modeling, we revealed the fundamental nanoscale electrical, mechanical, and magnetic interactions in the nanomotor system, which excellently agrees with experimental results and provides critical understanding for designing future metallic NEMS devices. The nanowire motors were applied for controlled biomolecule release, and for the first time demonstrated releasing rate of biochemicals on nanoparticles can be tuned by mechanical rotation. The innovations reported in this research, from concept, design, actuation, to application are relevant to NEMS, nanomedicine, microfluidics, and lab-on-chip architectures.


title: Two mod-p Johnson filtrations

abstract: We consider two mod-p central series of the free group given by Stallings and Zassenhaus. Applying these series to definitions of Dennis Johnson's filtration of the mapping class group we obtain two mod-p Johnson filtrations. Further, we adapt the definition of the Johnson homomorphisms to obtain mod-p Johnson homomorphisms. We calculate the image of the first of these homomorphisms. We give generators for the kernels of these homomorphisms as well. We restrict the range of our mod-p Johnson homomorphisms using work of Morita. We finally prove the announced result of Perron that a rational homology 3-sphere may be given as a Heegaard splitting with gluing map coming from certain members of our mod-p Johnson filtrations.


title: Invariantly Admissible Policy Iteration for a Class of Nonlinear Optimal Control Problems

abstract: In this paper, we propose a generalized successive approximation method (SAM), called invariantly admissible policy iteration (PI), for finding the solution to a class of input-affine nonlinear optimal control problems by iterations. Unlike the existing SAM, the proposed method updates the domain of the next policy and value function for admissibility (and invariance). In the existing SAM, the admissibility of the generated policies are guaranteed under the two implicit assumptions regarding Lyapunov's theorem and invariance, both of which are presented and discussed in this paper and are generally not true. On the contrary, the proposed invariantly admissible PI guarantees the admissibility in a more refined manner, without such assumptions. The admissibility and invariance of the updated region, with respect to the corresponding policies, are mathematically prove under the specific invariant admissible update rule. We also provide monotonic decreasing and uniform convergence properties of the sequence of value functions under certain conditions. Finally, numerical simulations are presented to illustrate the proposed PI method and its effectiveness.


title: Validating a Time-Dependent Wave-Turbulence-Driven Model of the Solar Wind

abstract: Although the mechanisms responsible for heating the Sun's corona and accelerating the solar wind are still being actively investigated, it is largely accepted that photospheric motions provide the energy source and that the magnetic field must play a key role in the process. \citet{2010ApJ...708L.116V} presented a model for heating and accelerating the solar wind based on the turbulent dissipation of Alfv\'en waves. We first use a time-dependent model of the solar wind to reproduce one of \citeauthor{2010ApJ...708L.116V}'s solutions; then we extend its application to the case when the energy equation includes thermal conduction and radiation losses, and the upper chromosphere is part of the computational domain. Using this model, we explore parameter space and describe the characteristics of a fast-solar-wind solution. We discuss how this formulation may be applied to a 3D MHD model of the coron a and solar wind \citep{2009ApJ...690..902L}.


title: Computations of Nonlinear Propagation of Sound Emitted from High Speed Mixing Layers

abstract: Non-linear sound propagation is investigated computationally by simulating compressible time-developing mixing layers using the Large Eddy Simulation (LES) approach and solving the viscous Burgers Equation. The mixing layers are of convective Mach numbers of 0.4, 0.8 and 1.2. The LES results agree qualitatively with known flow behavior. Mach waves are observed in the near sound field of the supersonic mixing layer computed by the LES. These waves show steepening typical to non-linear propagation. Further calculations using the Burgers equation support this finding, where the initial wave slope has a role in kicking them. No visible non-linear propagation effects were found for the subsonic mixing layers. The effects of geometrical spreading and viscosity are also considered.


title: A Typical Medium Dynamical Cluster Approximation for the Study of Anderson Localization in Three Dimensions

abstract: We develop a systematic typical medium dynamical cluster approximation that provides a proper description of the Anderson localization transition in three dimensions (3D). Our method successfully captures the localization phenomenon both in the low and large disorder regimes, and allows us to study the localization in different momenta cells, which renders the discovery that the Anderson localization transition occurs in a cell-selective fashion. As a function of cluster size, our method systematically recovers the re-entrance behavior of the mobility edge and obtains the correct critical disorder strength for Anderson localization in 3D.


title: Effect of dispersion interactions on the properties of LiF in condensed phases

abstract: Classical molecular dynamics simulations are performed on LiF in the framework of the polarizable ion model. The overlap-repulsion and polarization terms of the interaction potential are derived on a purely non empirical, first-principles basis. For the dispersion, three cases are considered: a first one in which the dispersion parameters are set to zero and two others in which they are included, with different parameterizations. Various thermodynamic, structural and dynamic properties are calculated for the solid and liquid phases. The melting temperature is also obtained by direct coexistence simulations of the liquid and solid phases. Dispersion interactions appear to have an important effect on the density of both phases and on the melting point, although the liquid properties are not affected when simulations are performed in the NVT ensemble at the experimental density.


title: Tunable Band Gap and Anisotropic Optical Response in Few-layer Black Phosphorus

abstract: We report the quasiparticle band gap, excitons, and highly anisotropic optical responses of few-layer black phosphorous (phosphorene). It is shown that these new materials exhibit unique many-electron effects; the electronic structures are dispersive essentially along one dimension, leading to particularly enhanced self-energy corrections and excitonic effects. Additionally, within a wide energy range, including infrared light and part of visible light, few-layer black phosphorous absorbs light polarized along the structure's armchair direction and is transparent to light polarized along the zigzag direction, making them viable linear polarizers for applications. Finally, the number of phosphorene layers included in the stack controls the material's band gap, optical absorption spectrum, and anisotropic polarization energy-window across a wide range.


title: Is gauge-invariant complete decomposition of the nucleon spin possible ?

abstract: Is gauge-invariant complete decomposition of the nucleon spin possible? Although it is a difficult theoretical question which has not reached a complete consensus yet, a general agreement now is that there are at least two physically inequivalent gauge-invariant decompositions (I) and (II) of the nucleon. %The one is a nontrivial gauge-invariant %generalization of the Jaffe-Manohar decomposition. %The other is an extension of the Ji decomposition, which allows %a gauge-invariant decomposition of the total gluon angular %momentum into the intrinsic spin and orbital parts. In these two decompositions, the intrinsic spin parts of quarks and gluons are just common. What discriminate these two decompositions are the orbital angular momentum parts. The orbital angular momenta of quarks and gluons appearing in the decomposition (I) are the so-called "mechanical" orbital angular momenta, while those appearing in the decomposition (II) are the generalized (gauge-invariant) "canonical" ones. By this reason, these decompositions are also called the "mechanical" and "canonical" decompositions of the nucleon spin, respectively. A crucially important question is which decomposition is more favorable from the observational viewpoint. The main objective of this concise review is to try to answer this question with careful consideration of recent intensive researches on this problem.


title: On the Hardness of Signaling

abstract: There has been a recent surge of interest in the role of information in strategic interactions. Much of this work seeks to understand how the realized equilibrium of a game is influenced by uncertainty in the environment and the information available to players in the game. Lurking beneath this literature is a fundamental, yet largely unexplored, algorithmic question: how should a "market maker" who is privy to additional information, and equipped with a specified objective, inform the players in the game? This is an informational analogue of the mechanism design question, and views the information structure of a game as a mathematical object to be designed, rather than an exogenous variable. We initiate a complexity-theoretic examination of the design of optimal information structures in general Bayesian games, a task often referred to as signaling. We focus on one of the simplest instantiations of the signaling question: Bayesian zero-sum games, and a principal who must choose an information structure maximizing the equilibrium payoff of one of the players. In this setting, we show that optimal signaling is computationally intractable, and in some cases hard to approximate, assuming that it is hard to recover a planted clique from an Erdos-Renyi random graph. This is despite the fact that equilibria in these games are computable in polynomial time, and therefore suggests that the hardness of optimal signaling is a distinct phenomenon from the hardness of equilibrium computation. Necessitated by the non-local nature of information structures, en-route to our results we prove an "amplification lemma" for the planted clique problem which may be of independent interest.


title: Electron g-2 in Light-Front Quantization

abstract: Basis Light-front Quantization has been proposed as a nonperturbative framework for solving quantum field theory. We apply this approach to Quantum Electrodynamics and explicitly solve for the light-front wave function of a physical electron. Based on the resulting light-front wave function, we evaluate the electron anomalous magnetic moment. Nonperturbative mass renormalization is performed. Upon extrapolation to the infinite basis limit our numerical results agree with the Schwinger result obtained in perturbation theory to an accuracy of 0.06%.


title: Spontaneous Breaking of Scale Invariance in a d=3 U(N) Model with Chern-Simons Gauge Field

abstract: We study spontaneous breaking of scale invariance in the large N limit of three dimensional $U(N)_\kappa$ Chern-Simons theories coupled to a scalar field in the fundamental representation. When a $\lambda_6(\phi^\dagger\cdot\phi)^3$ self interaction term is added to the action we find a massive phase at a certain critical value for a combination of the $\lambda_6$ and 't Hooft's $\lambda=N/\kappa$ couplings. This model attracted recent attention since at finite $\kappa$ it contains a singlet sector which is conjectured to be dual to Vasiliev's higher spin gravity on $AdS_4$. Our paper concentrates on the massive phase of the 3d boundary theory. We discuss the advantage of introducing masses in the boundary theory through spontaneous breaking of scale invariance.


title: Adjoinable Homology

abstract: The notion of a duality between two derived functors as well as an extension theorem for derived functors to larger categories in which they need not be defined is introduced. These ideas are then applied to extend and study the coext functors to an arbitrary coalgebra. A new homology theory theory is then built therefrom and is shown to exhibit certain duality relations to the Hochschild cohomology of certain coalgebras. Lastly, a certain exceptional type of coalgebra is introduced and it is used to make explicit connections between this new homology theory and the continuous cohomology of this exceptional algebra's pro-finite dual algebra.


title: An SDP Approach For Solving Quadratic Fractional Programming Problems

abstract: This paper considers a fractional programming problem (P) which minimizes a ratio of quadratic functions subject to a two-sided quadratic constraint. As is well-known, the fractional objective function can be replaced by a parametric family of quadratic functions, which makes (P) highly related to, but more difficult than a single quadratic programming problem subject to a similar constraint set. The task is to find the optimal parameter $\lambda^*$ and then look for the optimal solution if $\lambda^*$ is attained. Contrasted with the classical Dinkelbach method that iterates over the parameter, we propose a suitable constraint qualification under which a new version of the S-lemma with an equality can be proved so as to compute $\lambda^*$ directly via an exact SDP relaxation. When the constraint set of (P) is degenerated to become an one-sided inequality, the same SDP approach can be applied to solve (P) {\it without any condition}. We observe that the difference between a two-sided problem and an one-sided problem lies in the fact that the S-lemma with an equality does not have a natural Slater point to hold, which makes the former essentially more difficult than the latter. This work does not, either, assume the existence of a positive-definite linear combination of the quadratic terms (also known as the dual Slater condition, or a positive-definite matrix pencil), our result thus provides a novel extension to the so-called "hard case" of the generalized trust region subproblem subject to the upper and the lower level set of a quadratic function.


title: Sums of Squares on the Hypercube

abstract: Let X be a finite set of points in R^n. A polynomial p nonnegative on X can be written as a sum of squares of rational functions modulo the vanishing ideal I(X). From the point of view of applications, such as polynomial optimization, we are interested in rational function representations of small degree. We derive a general upper bound in terms of the Hilbert function of X, and we show that this upper bound is tight for the case of quadratic functions on the hypercube C={0,1}^n, a very well studied case in combinatorial optimization. Using the lower bounds for C we construct a family of globally nonnegative quartic polynomials, which are not sums of squares of rational functions of small degree. To our knowledge this is the first construction for Hilbert's 17th problem of a family of polynomials of bounded degree which need increasing degrees in rational function representations as the number of variables n goes to infinity. We note that representation theory of the symmetric group S_n play a crucial role in our proofs of the lower bounds.


title: Moduli spaces of SO(8) instantons on smooth ALE spaces as Higgs branches of 4d N=2 supersymmetric theories

abstract: The worldvolume theory of D3-branes probing four D7-branes and an O7-plane on $\mathbb{C}^2/\mathbb{Z}_2$ is given by a supersymmetric USp x USp gauge theory. We demonstrate that, at least for a particular choice of the holonomy at infinity, we can go to a dual description of the gauge theory, where we can add a Fayet-Iliopoulos term describing the blowing-up of the orbifold to the smooth ALE space. This allows us to express the moduli space of SO(8) instantons on the smooth ALE space as a hyperk\"ahler quotient of a flat space times the Higgs branch of a class S theory. We also discuss a generalization to $\mathbb{C}^2/\mathbb{Z}_{2n}$, and speculate how to extend the analysis to bigger SO groups and to ALE spaces of other types.


title: Reciprocity sheaves

abstract: We start developing a notion of reciprocity sheaves, generalizing Voevodsky's homotopy invariant presheaves with transfers which were used in the construction of his triangulated categories of motives. We hope reciprocity sheaves will eventually lead to the definition of a larger triangulated category of motivic nature, encompassing non homotopy invariant phenomena.


title: A Liouville theorem for minimizers with finite potential energy for the vectorial Allen-Cahn equation

abstract: We prove that if a globally minimizing solution to the vectorial Allen-Cahn equation has finite potential energy, then it is a constant.


title: Higgs bundles and local systems on Riemann surfaces

abstract: Lecture notes from the Third International School on Geometry and Physics at the Centre de Recerca Matematica in Barcelona, March 26--30, 2012.


title: Electroweak baryogenesis in the MSSM with vector-like superfields

abstract: Introducing heavy particles with strong couplings to the Higgs field can strengthen electroweak phase transition, through the entropy release mechanism from both bosons and fermions. We analyze the possibility of electroweak baryogenesis in the MSSM with new vector-like superfields. The new vector-like particles belong to the representation 5+5bar+10+10bar of SU(5). By analyzing in detail the effective potential at finite temperature, we show that a strongly first order electroweak phase transition in this model is ruled out by a combination of 125 GeV Higgs requirement, the bound for exotic quarks, the gluon fusion Higgs production rate and the Higgs diphoton decay rate as well as the electroweak precision measurement.


title: Method of Studying $\Lambda_b^0$ decays with one missing particle

abstract: A new technique is discussed that can be applied to $\Lambda_b^0$ baryon decays where decays with one missing particle can be discerned from background and their branching fractions determined, along with other properties of the decays. Applications include measurements of the CKM elements $|V_{ub}|$ and $|V_{cb}|$, selected charmless decays, and detection of any exotic objects coupling to $b\to s$ decays, such as the inflaton.


title: Stress relaxation models with polyconvex entropy in Lagrangean and Eulerian coordinates

abstract: The embedding of the equations of polyconvex elastodynamics to an augmented symmetric hyperbolic system provides in conjunction with the relative entropy method a robust stability framework for approximate solutions. We devise here a model of stress relaxation motivated by the format of the enlargement process which formally approximates the equations of polyconvex elastodynamics. The model is endowed with an entropy function which is not convex but rather of polyconvex type. Using the relative entropy we prove a stability estimate and convergence of the stress relaxation model to polyconvex elastodynamics in the smooth regime. As an application, we show that models of pressure relaxation for real gases in Eulerian coordinates fit into the proposed framework.


title: Satellite structures in the spectral functions of semiconductor quantum well two-dimensional electron gases: a GW plus cumulant study

abstract: We present theoretical calculations for the spectral functions and single-particle densities of states of the two-dimensional electron gas in semiconductor quantum wells at different electron densities using the GW+cumulant method. We compare our results to GW only calculations and find significant differences in the description of the satellites between the two theories: while GW theory predicts the existence of a plasmaron excitation, no such excitation is found in GW+cumulant theory. We compare our results to experimental tunneling spectra from semiconductor quantum wells and find good agreement for the satellite properties.


title: Gru\v{s}in operators, Riesz transforms and nilpotent Lie groups

abstract: We establish that the Riesz transforms of all orders corresponding to the Gru\v{s}in operator $H_N=-\nabla_{x}^2-|x|^{2N}\,\nabla_{y}^2$, and the first-order operators $(\nabla_{x},x^\nu\,\nabla_{y})$ where $x\in \Ri^n$, $y\in\Ri^m$, $N\in\Ni_+$, and $\nu\in\{1,\ldots,n\}^N$, are bounded on $L_p(\Ri^{n+m})$ for all $p\in\langle1,\infty\rangle$ and are also weak-type $(1,1)$. Moreover, the transforms of order less than or equal to $N+1$ corresponding to $H_N$ and the operators $(\nabla_{x}, |x|^N\nabla_{y})$ are bounded on $L_p(\Ri^{n+m})$ for all $p\in\langle1,\infty\rangle$. But all transforms of order $N+2$ are bounded if and only if $p\in\langle1,n\rangle$. The proofs are based on the observation that the $(\nabla_{x},x^\nu\,\nabla_{y})$ generate a finite-dimensional nilpotent Lie algebra, the corresponding connected, simply connected, nilpotent Lie group is isometrically represented on the spaces $L_p(\Ri^{n+m})$ and $H_N$ is the corresponding sublaplacian


title: On non-Archimedean recurrence equations and their applications

abstract: In the present paper we study stability of recurrence equations (which in particular case contain a dynamics of rational functions) generated by contractive functions defined on an arbitrary non-Archimedean algebra. Moreover, multirecurrence equations are considered. We also investigate reverse recurrence equations which have application in the study of $p$-adic Gibbs measures. Note that our results also provide the existence of unique solutions of nonlinear functional equations as well.


title: Nonadiabatic Dynamics of a Dissipative Two-level System

abstract: We study the dynamics of a two-level system described by a slowly varying Hamiltonian and weakly coupled to the Ohmic environment. We follow the Bloch--Redfield perturbative approach to include the effect of the environment on qubit evolution and take into account modification of the spectrum and matrix elements of qubit transitions due to time-dependence of the Hamiltonian. This formalism is applied to two problems. (1) We consider a qubit, or a spin-1/2, in a rotating magnetic field. We show that once the rotation starts, the spin has a component perpendicular to the rotation plane of the field that initially wiggles and eventually settles to the value proportional to the product of angular rotation velocity of the field and the Berry curvature. (2) We re-examine the Landau--Zener transition for a system coupled to environment at arbitrary temperature. We show that as temperature increases, the thermal excitation and relaxation become leading processes responsible for transition between states of the system. We also apply the Lindblad master equations to these two problems and compare results with those obtained from the Bloch--Redfield equations.


title: Exploring gender differences on general and specific computer self-efficacy in mobile learning adoption

abstract: Reasons for contradictory findings regarding the gender moderate effect on computer self-efficacy in the adoption of e-learning/mobile learning are limited. Recognizing the multilevel nature of the computer self-efficacy (CSE), this study attempts to explore gender differences in the adoption of mobile learning, by extending the Technology Acceptance Model (TAM) with general and specific CSE. Data collected from 137 university students were tested against the research model using the structural equation modeling approach. The results suggest that there are significant gender differences in perceptions of general CSE, perceived ease of use and behavioral intention to use but no significant differences in specific CSE, perceived usefulness. Additionally, the findings reveal that specific CSE is more salient than general CSE in influencing perceived ease of use while general CSE seems to be the salient factor on perceived usefulness for both female and male combined. Moreover, general CSE was salient to determine the behavioral intention to use indirectly for female despite lower perception of general CSE than male's, and specific CSE exhibited stronger indirect effect on behavioral intention to use than general CSE for female despite similar perception of specific CSE as males'. These findings provide important implications for mobile learning adoption and usage.


title: Non-supersymmetric Microstates of the MSW System

abstract: We present an analysis parallel to that of Giusto, Ross, and Saxena (arXiv:0708.3845) and construct a discrete family of non-supersymmetric microstate geometries of the Maldacena-Strominger-Witten system. The supergravity configuration in which we look for the smooth microstates is constructed using SO(4,4) dualities applied to an appropriate seed solution. The SO(4,4) approach offers certain technical advantages. Our microstate solutions are smooth in five dimensions, as opposed to all previously known non-supersymmetric microstates with AdS3 cores, which are smooth only in six dimensions. The decoupled geometries for our microstates are related to global AdS3 x S2 by spectral flows.


title: On homomorphisms from Ringel-Hall algebras to quantum cluster algebras

abstract: In \cite{rupel3},the authors defined algebra homomorphisms from the dual Ringel-Hall algebra of certain hereditary abelian category $\mathcal{A}$ to an appropriate $q$-polynomial algebra. In the case that $\mathcal{A}$ is the representation category of an acyclic quiver, we give an alternative proof by using the cluster multiplication formulas in \cite{DX}. Moreover, if the underlying graph of $Q$ is bipartite and the matrix $B$ associated to the quiver $Q$ is of full rank, we show that the image of the algebra homomorphisms is in the corresponding quantum cluster algebra.


title: Quantum Yields from Stationary States: Cis-Trans Isomerization of Model Retinal

abstract: Cis-trans isomerization in retinal, the first step in vision, is often computationally studied from a time dependent viewpoint. Motivation for such studies lies in coherent pulsed laser experiments that explore the isomerization dynamics. However, such biological processes take place naturally in the presence of incoherent light, which excites a non-evolving mixture of stationary states. Here the isomerization problem is considered from the latter viewpoint and applied to a standard two-state, two-mode linear vibronic coupling model of retinal that explicitly includes a conical intersection between the ground and first excited electronic states. The calculated quantum yield at 500 nm agrees well with both the previous time-dependent calculations of Hahn and Stock (0.63) and with experiment ($0.65\pm0.01$), as does its wavelength dependence. Significantly, the effects of environmental relaxation on the quantum yield in this well-established model are found to be negligible. The results make clear the connection of the photoisomerization quantum yield to properties of stationary eigenstates, providing alternate insights into conditions for yield optimization.


title: The greybody factor for scalar fields in the Schwarzschild spacetime with an $f(R)$ global monopole

abstract: The greybody factor of massless scalar fields in the four-dimensional Schwarzschild spacetime involving an $f(R)$ global monopole is derived. We show how the monopole parameter and the deviation from the standard general relativity adjust the greybody factor. We also demonstrate that the effects from the global monopole and $f(R)$ gravity theory are manifest in the energy emission rate and the generalized absorption cross section of the scalar fields.


title: The iPod Model

abstract: We introduce a Voter Model variant, inspired by social evolution of musical preferences. In our model, agents have preferences over a set of songs and upon meeting update their own preferences incrementally towards those of the other agents they meet. Using the spectral gap of an associated Markov chain, we give a geometry dependent result on the asymptotic consensus time of the model.


title: Mass-Galaxy offsets in Abell 3827, 2218 and 1689: intrinsic properties or line-of-sight substructures?

abstract: We have made mass maps of three strong-lensing clusters, Abell 3827, Abell 2218 and Abell 1689, in order to test for mass-light offsets. The technique used is GRALE, which enables lens reconstruction with minimal assumptions, and specifically with no information about the cluster light being given. In the first two of these clusters, we find local mass peaks in the central regions that are displaced from the nearby galaxies by a few to several kpc. These offsets {\em could\/} be due to line of sight structure unrelated to the clusters, but that is very unlikely, given the typical levels of chance line-of-sight coincidences in $\Lambda CDM$ simulations --- for Abell 3827 and Abell 2218 the offsets appear to be intrinsic. In the case of Abell 1689, we see no significant offsets in the central region, but we do detect a possible line of sight structure: it appears only when sources at $z\ga 3$ are used for reconstructing the mass. We discuss possible origins of the mass-galaxy offsets in Abell 3827 and Abell 2218: these include pure gravitational effects like dynamical friction, but also non-standard mechanisms like self-interacting dark-matter.


title: Surface atomic diffusion processes observed at milliseconds time resolution using environmental TEM

abstract: Significant progress has been made in spatial resolution using environmental transmission electron microscopes (ETEM), which now enables atomic resolution visualization of structural transformation under variable temperature and gas environments close to materials real operational conditions. Structural transformations are observed by recording images or diffraction patterns at various time intervals using a video camera or by taking snap shots using electron pulses. While time resolution at 15 ns has been reported using pulsed electron beams, the time interval that can be recorded by this technique is currently very limited. For longer recording, however, time resolution inside ETEM has been limited by electron cameras to ~1/30 seconds for a long time. Using the recently developed direct electron detection technology, we have significantly improved the time resolution of ETEM to 2.5 ms (milliseconds) for full frame or 0.625 ms for 0.25 frames.


title: Detection of bridge emission above 50 GeV from the Crab pulsar with the MAGIC telescopes

abstract: The Crab pulsar is the only astronomical pulsed source detected at very high energy (VHE, E>100GeV) gamma-rays. The emission mechanism of VHE pulsation is not yet fully understood, although several theoretical models have been proposed. In order to test the new models, we measured the light curve and the spectra of the Crab pulsar with high precision by means of deep observations. We analyzed 135 hours of selected MAGIC data taken between 2009 and 2013 in stereoscopic mode. In order to discuss the spectral shape in connection with lower energies, 4.6 years of {\it Fermi}-LAT data were also analyzed. The known two pulses per period were detected with a significance of $8.0 \sigma$ and $12.6 \sigma$. In addition, significant emission was found between the two pulses with $6.2 \sigma$. We discovered the bridge emission above 50 GeV between the two main pulses. This emission can not be explained with the existing theories. These data can be used for testing new theoretical models.


title: Functional Models and Minimal Contractive Liftings

abstract: Based on a careful analysis of functional models for contractive multi-analytic operators we establish a one-to-one correspondence between unitary equivalence classes of minimal contractive liftings of a row contraction and injective symbols of contractive multi-analytic operators. This allows an effective construction and classification of all such liftings with given defects. Popescu's theory of characteristic functions of completely non-coisometric row contractions is obtained as a special case satisfying a Szeg\"{o} condition. In another special case of single contractions and defects equal to $1$ all non-zero Schur functions on the unit disk appear in the classification. It is also shown that the process of constructing liftings iteratively reflects itself in a factorization of the corresponding symbols.


title: Blow-up formulae of high genus Gromov-Witten invariants in dimensional six

abstract: Using the degeneration formula and absolute/relative correspondence, one studied the change of Gromov-Witten invariants under blow-up for six dimensional symplectic manifolds and obtained closed blow-up formulae for high genus Gromov-Witten invariants. Our formulae also imply some relations among generalized BPS numbers introduced by Pandharipande.


title: MAP, MAC, and Vortex-rings Configurations in the Weinberg-Salam Model

abstract: We report on the presence of new axially symmetric monopoles, antimonopoles and vortex-rings solutions of the SU(2)$\times$U(1) Weinberg-Salam model of electromagnetic and weak interactions. When the $\phi$-winding number $n=1$, and 2, the configurations are monopole-antimonopole pair (MAP) and monopole-antimonopole chain (MAC) with poles of alternating sign magnetic charge arranged along the $z$-axis. Vortex-rings start to appear from the MAP and MAC configurations when the winding number $n=3$. The MAP configurations possess zero net magnetic charge whereas the MAC configurations possess net magnetic charge of $4\pi n/e$. In the MAP configurations, the monopole-antimonopole pair is bounded by the ${\cal Z}^0$ field flux string and there is an electromagnetic current loop encircling it. The monopole and antimonopole possess magnetic charges $\pm\frac{4\pi n}{e}\sin^2\theta_W$ respectively. In the MAC configurations there is no string connecting the monopole and the adjacent antimonopole and they possess magnetic charges $\pm\frac{4\pi n}{e}$ respectively. The MAC configurations possess infinite total energy and zero magnetic dipole moment whereas the MAP configurations which are actually sphalerons possess finite total energy and magnetic dipole moment. The configurations were investigated for varying values of Higgs self-coupling constant $0\leq \lambda\leq 40$ at Weinberg angle $\theta_W=\frac{\pi}{4}$.


title: On the origins and the historical roots of the Higgs boson research from a bibliometric perspective

abstract: Subject of our present paper is the analysis of the origins or historical roots of the Higgs boson research from a bibliometric perspective, using a segmented regression analysis in a reference publication year spectroscopy (RPYS). Our analysis is based on the references cited in the Higgs boson publications published since 1974. The objective of our analysis consists of identifying concrete individual publications in the Higgs boson research context to which the scientific community frequently had referred to. As a consequence, we are interested in seminal works which contributed to a high extent to the discovery of the Higgs boson. Our results show that researchers in the Higgs boson field preferably refer to more recently published papers - particular papers published since the beginning of the sixties. For example, our analysis reveals seven major contributions which appeared within the sixties: Englert and Brout (1964), Higgs (1964, 2 papers), and Guralnik et al. (1964) on the Higgs mechanism as well as Glashow (1961), Weinberg (1967), and Salam (1968) on the unification of weak and electromagnetic interaction. Even if the Nobel Prize award highlights the outstanding importance of the work of Peter Higgs and Francois Englert, bibliometrics offer the additional possibility of getting hints to other publications in this research field (especially to historical publications), which are of vital importance from the expert point of view.


title: Asymptotic of 'rigid-body' motions for nonlinear dynamics: physical insight and methodologies

abstract: The purpose of the present work is to show that an adequate basis for understanding the essentially nonlinear phenomena must also be essentially nonlinear however still simple enough to play the role of a basis. It is shown that such types of 'elementary' nonlinear models can be revealed by tracking the hidden links between analytical tools of analyses and subgroups of the rigid-body motions or, in other terms, rigid Euclidean transformation. While the subgroup of rotations is linked with linear and weakly nonlinear vibrations, the translations with reflections can be viewed as a geometrical core of the strongly nonlinear dynamics associated with the so-called vibro-impact behaviors. It is shown that the corresponding analytical approach develops through non-smooth temporal substitutions generated by the impact models.


title: Information Theory of Matrix Completion

abstract: Matrix completion is a fundamental problem that comes up in a variety of applications like the Netflix problem, collaborative filtering, computer vision, and crowdsourcing. The goal of the problem is to recover a k-by-n unknown matrix from a subset of its noiseless (or noisy) entries. We define an information-theoretic notion of completion capacity C that quantifies the maximum number of entries that one observation of an entry can resolve. This number provides the minimum number m of entries required for reliable reconstruction: m=kn/C. Translating the problem into a distributed joint source-channel coding problem with encoder restriction, we characterize the completion capacity for a wide class of stochastic models of the unknown matrix and the observation process. Our achievability proof is inspired by that of the Slepian-Wolf theorem. For an arbitrary stochastic matrix, we derive an upper bound on the completion capacity.


title: Oscillatory instability and fluid patterns in low-Prandtl-number Rayleigh-B\'{e}nard convection with uniform rotation

abstract: We present the results of direct numerical simulations of flow patterns in a low-Prandtl-number ($Pr = 0.1$) fluid above the onset of oscillatory convection in a Rayleigh-B\'{e}nard system rotating uniformly about a vertical axis. Simulations were carried out in a periodic box with thermally conducting and stress-free top and bottom surfaces. We considered a rectangular box ($L_x \times L_y \times 1$) and a wide range of Taylor numbers ($750 \le Ta \le 5000$) for the purpose. The horizontal aspect ratio $\eta = L_y/L_x$ of the box was varied from $0.5$ to $10$. The primary instability appeared in the form of two-dimensional standing waves for shorter boxes ($0.5 \le \eta < 1$ and $1 < \eta < 2$). The flow patterns observed in boxes with $\eta = 1$ and $\eta = 2$ were different from those with $\eta < 1$ and $1 < \eta < 2$. We observed a competition between two sets of mutually perpendicular rolls at the primary instability in a square cell ($\eta = 1$) for $Ta < 2700$, but observed a set of parallel rolls in the form of standing waves for $Ta \geq 2700$. The three-dimensional convection was quasiperiodic or chaotic for $750 \le Ta < 2700$, and then bifurcated into a two-dimensional periodic flow for $Ta \ge 2700$. The convective structures consisted of the appearance and disappearance of straight rolls, rhombic patterns, and wavy rolls inclined at an angle $\phi = \frac{\pi}{2} - \arctan{(\eta^{-1})}$ with the straight rolls.


title: Mass spectra and decay properties of $D_s$ Meson in a relativistic Dirac formalism

abstract: The mass spectra of $D_s$ meson is obtained in the framework of relativistic independent quark model using Martin like potential for the quark confinement. The predicted excited states are in good agreement with the experimental results as well as with the lattice and other theoretical predictions. The spectroscopic parameters are employed further to compute the decay constant, electromagnetic transition and leptonic decay widths. The present result for its decay constant, $f_P$ (252.82 MeV) is in excellent agreement with the value 252.6 $\pm$ 11.1 MeV reported by CLEO-c and the predicted branching ratios for $(D_s \rightarrow \tau \bar{\nu}_\tau, \mu \bar{\nu}_\mu)$ ($5.706 \times 10^{-2}, 5.812 \times 10^{-3}$) are in close agreement with the PDG values ($ (5.43 \pm 0.31)\times 10^{-2}, (5.90 \pm 0.33)\times 10^{-3} $) respectively.


title: A question of Doctor Malte Wandel on automorphisms of the punctural Hilbert schemes of K3 surfaces

abstract: We present a sufficient condition for the punctural Hilbert scheme of length two of a K3 surface with finite automorphism group to have automorphism group of infinite order in geometric terms (Theorem 2.1). We then give concrete examples (Theorem 1.2). We also discuss about Mori dream space (MDS) structures under an extermal crepant resolution (Theorems 1.2, 5.2, 4.1) from the viewpoint of automorphisms. These affirmatively answer a question of Doctor Malte Wandel.


title: Realization of white light cavity for ulra-short optical signal storage and Processing

abstract: We have demonstrated that chirped Bragg gratings can compensate for a chromatic dispersion by reflecting different wavelengths at different location along the axis of the gratings. In this paper, we make use of such a dispersion-compensating property to create a white light cavity (WLC) effect that cancels the round trip phase shift due to propagation as the optical frequency is tuned. A pair of LCBGs is used to make the cavity. To fulfill WLC condition, the dispersion profile of the LCBGs is tailored by adjusting parameters such as chirping rate, modulation depth and the length of the grating region. Numerical simulation shows WLC bandwidth is ~10GHz.


title: On the effect of intermittency of turbulence on the parabolic relation between skewness and kurtosis in magnetized plasmas

abstract: This paper is aimed to contribute to the scientific discussions that have been triggered by the experimental observation of a quadratic relation between the kurtosis and skewness of turbulent fluctuations present in fusion plasmas and other nonlinear physical systems. In this paper we offer a general statistical model which attributes the observed $K = aS^2 + b$ relation to the varying intermittency of the experimental signals. The model is a two random variable model constructed to catch the essential intermittent feature of the real signal. One of the variables is the amplitude of the underlying intermittent event (e.g. turbulent structure) while the other is connected to the intermittency level of the system. This simple model can attribute physical meaning to the $a$ and $b$ coefficients, as they characterize the spatio-temporal statistics of intermittent events. By constructing a particle-conserving Gaussian model for the underlying coherent structures the experimentally measured $a$ and $b$ coefficients could be adequately reproduced.


title: Some Centrally Symmetric Manifolds

abstract: We show existence of centrally symmetric maps on surfaces all of whose faces are quadrangles and pentagons for each orientable genus $g \geq 0$. We also show existence of centrally symmetric maps on surfaces all of whose faces are hexagons for each orientable genus $g = 2k-1$, $k\in \mathbb{N}$. We enumerate centrally symmetric triangulated manifolds of dimensions 2 and 3 with few vertices.


title: Harnack Estimates for Nonlinear Backward Heat Equations in Geometric Flows

abstract: Let $M$ be a closed Riemannian manifold with a family of Riemannian metrics $g_{ij}(t)$ evolving by a geometric flow $\partial_{t}g_{ij} = -2{S}_{ij}$, where $S_{ij}(t)$ is a family of smooth symmetric two-tensors. We derive several differential Harnack estimates for positive solutions to the nonlinear backward heat-type equation \begin{eqnarray*} \frac{\partial f}{\partial t} = -{\Delta}f + \gamma f\log f +aSf \end{eqnarray*} where $a$ and $\gamma$ are constants and $S=g^{ij}S_{ij}$ is the trace of $S_{ij}$. Our abstract formulation provides a unified framework for some known results proved by various authors, and moreover lead to new Harnack inequalities for a variety of geometric flows.


title: On the Phase Space of Block-Hiding Strategies in Bitcoin-like networks

abstract: We calculate the probability of success of block-hiding mining strategies in Bitcoin-like networks. These strategies involve building a secret branch of the block-tree and publishing it opportunistically, aiming to replace the top of the main branch and rip the reward associated with the secretly mined blocks. We identify two types of block-hiding strategies and chart the parameter space where those are more beneficial than the standard mining strategy described in Nakamoto's paper. Our analysis suggests a generalization of the notion of the relative hashing power as a measure for a miner's influence on the network. Block-hiding strategies are beneficial only when this measure of influence exceeds a certain threshold.


title: Time-Optimal Universal Control of Two-Level Systems under Strong Driving

abstract: We experimentally study the time-optimal construction of arbitrary single-qubit rotations under a single strong driving field of finite amplitude. Using radiation-dressed states of nitrogen vacancy centers in diamond, we realize a strongly-driven two-level system and achieve driving frequencies four times larger than its Larmor frequency. We implement time optimal universal rotations on this system, characterize their performance using quantum process tomography, and demonstrate a dual-axis ac magnetometry sequence with pulses at sub-Larmor time scales. Our results pave the way for applying fast qubit control and high-density pulse schemes in the fields of quantum information processing and quantum metrology.


title: Signifying quantum benchmarks for qubit teleportation and secure communication using Einstein-Podolsky-Rosen steering inequalities

abstract: The demonstration of quantum teleportation of a photonic qubit from Alice to Bob usually relies on data conditioned on detection at Bob's location. I show that Bohm's Einstein-Podolsky-Rosen (EPR) paradox can be used to verify that the quantum benchmark for qubit teleportation has been reached, without postselection. This is possible for scenarios insensitive to losses at the generation station, and with efficiencies of $\eta_{B}>1/3$ for the teleportation process. The benchmark is obtained, if it is shown that Bob can {}"steer" Alice's record of the qubit as stored by Charlie. EPR steering inequalities involving $m$ measurement settings can also be used to confirm quantum teleportation, for efficiencies $\eta_{B}>1/m$, if one assumes trusted detectors for Charlie and Alice. Using proofs of monogamy, I show that two-setting EPR steering inequalities can signify secure teleportation of the qubit state.


title: Harnack Estimates for Nonlinear Heat Equations with Potentials in Geometric Flows

abstract: Let $M$ be a closed Riemannian manifold with a family of Riemannian metrics $g_{ij}(t)$ evolving by geometric flow $\partial_{t}g_{ij} = -2{S}_{ij}$, where $S_{ij}(t)$ is a family of smooth symmetric two-tensors on $M$. In this paper we derive differential Harnack estimates for positive solutions to the nonlinear heat equation with potential: \begin{eqnarray*} \frac{\partial f}{\partial t} = {\Delta}f + \gamma (t) f\log f +aSf, \end{eqnarray*} where $\gamma (t)$ is a continuous function on $t$, $a$ is a constant and $S=g^{ij}S_{ij}$ is the trace of $S_{ij}$. Our Harnack estimates include many known results as special cases, and moreover lead to new Harnack inequalities for a variety geometric flows.


title: Prediction of the thermophysical properties of molten salt fast reactor fuel from first-principles

abstract: Molten fluorides are known to show favorable thermophysical properties which make them good candidate coolants for nuclear fission reactors. Here we investigate the special case of mixtures of lithium fluoride and thorium fluoride, which act both as coolant and fuel in the molten salt fast reactor concept. By using ab initio parameterized polarizable force fields, we show that it is possible to calculate the whole set of properties (density, thermal expansion, heat capacity, viscosity and thermal conductivity) which are necessary for assessing the heat transfer performance of the melt over the whole range of compositions and temperatures. We then deduce from our calculations several figures of merit which are important in helping the optimization of the design of molten salt fast reactors.


title: Downlink and Uplink Energy Minimization Through User Association and Beamforming in Cloud RAN

abstract: The cloud radio access network (C-RAN) concept, in which densely deployed access points (APs) are empowered by cloud computing to cooperatively support mobile users (MUs), to improve mobile data rates, has been recently proposed. However, the high density of active ("on") APs results in severe interference and also inefficient energy consumption. Moreover, the growing popularity of highly interactive applications with stringent uplink (UL) requirements, e.g. network gaming and real-time broadcasting by wireless users, means that the UL transmission is becoming more crucial and requires special attention. Therefore in this paper, we propose a joint downlink (DL) and UL MU-AP association and beamforming design to coordinate interference in the C-RAN for energy minimization, a problem which is shown to be NP hard. Due to the new consideration of UL transmission, it is shown that the two state-of-the-art approaches for finding computationally efficient solutions of joint MU-AP association and beamforming considering only the DL, i.e., group-sparse optimization and relaxed-integer programming, cannot be modified in a straightforward way to solve our problem. Leveraging on the celebrated UL-DL duality result, we show that by establishing a virtual DL transmission for the original UL transmission, the joint DL and UL optimization problem can be converted to an equivalent DL problem in C-RAN with two inter-related subproblems for the original and virtual DL transmissions, respectively. Based on this transformation, two efficient algorithms for joint DL and UL MU-AP association and beamforming design are proposed, whose performances are evaluated and compared with other benchmarking schemes through extensive simulations.


title: Modeling Covariate Effects in Group Independent Component Analysis with Applications to Functional Magnetic Resonance Imaging

abstract: Independent component analysis (ICA) is a powerful computational tool for separating independent source signals from their linear mixtures. ICA has been widely applied in neuroimaging studies to identify and characterize underlying brain functional networks. An important goal in such studies is to assess the effects of subjects' clinical and demographic covariates on the spatial distributions of the functional networks. Currently, covariate effects are not incorporated in existing group ICA decomposition methods. Hence, they can only be evaluated through ad-hoc approaches which may not be accurate in many cases. In this paper, we propose a hierarchical covariate ICA model that provides a formal statistical framework for estimating and testing covariate effects in ICA decomposition. A maximum likelihood method is proposed for estimating the covariate ICA model. We develop two expectation-maximization (EM) algorithms to obtain maximum likelihood estimates. The first is an exact EM algorithm, which has analytically tractable E-step and M-step. Additionally, we propose a subspace-based approximate EM, which can significantly reduce computational time while still retain high model-fitting accuracy. Furthermore, to test covariate effects on the functional networks, we develop a voxel-wise approximate inference procedure which eliminates the needs of computationally expensive covariance estimation. The performance of the proposed methods is evaluated via simulation studies. The application is illustrated through an fMRI study of Zen meditation.


title: Boltzmann-equation approach to anomalous transport in a Weyl metal

abstract: Weyl metal is regarded as a platform toward interacting topological states of matter, where its topological structure gives rise to anomalous transport phenomena, referred to as chiral magnetic effect and "negative" magneto-resistivity, the origin of which is chiral anomaly. Recently, the negative magneto-resistivity has been observed with the signature of weak anti-localization at $x = 3 \sim 4 ~ \%$ in Bi$_{1-x}$Sb$_{x}$, where magnetic field is applied in parallel with electric field. Based on the Boltzmann-equation approach, we find the negative magneto-resistivity in the presence of weak anti-localization. An essential ingredient is to introduce the topological structure of chiral anomaly into the Boltzmann-equation approach, resorting to semi-classical equations of motion with Berry curvature.


title: Astrometric asteroid masses: a simultaneous determination

abstract: Using over 89 million astrometric observations of 349,737 numbered minor planets, an attempt was made to determine the masses of 230 of them by simultaneously solving for corrections to all orbital elements and the masses. For 132 minor planets an acceptable result was obtained, 50 of which appear to be new.


title: Constraints on the skewness coefficient of symmetric nuclear matter within the nonlinear relativistic mean field model

abstract: Within the nonlinear relativistic mean field (NL-RMF) model, we show that both the pressure of symmetric nuclear matter at supra-saturation densities and the maximum mass of neutron stars are sensitive to the skewness coefficient $J_0$ of symmetric nuclear matter. Using experimental constraints on the pressure of symmetric nuclear matter at supra-saturation densities from flow data in heavy ion collisions and the astrophysical observation of a large mass neutron star PSR J0348+0432, with the former favoring a smaller $J_0$ while the latter a larger $J_0$, we extract a constraint of $-494 \mathrm{MeV}\leq J_0\leq -10 \mathrm{MeV}$ based on the NL-RMF model. This constraint is compared with the results obtained in other analyses.


title: Searching for tau neutrinos with Cherenkov telescopes

abstract: Cherenkov telescopes have the capability of detecting high energy tau neutrinos in the energy range of 1--1000 PeV by searching for very inclined showers. If a tau lepton, produced by a tau neutrino, escapes from the Earth or a mountain, it will decay and initiate a shower in the air which can be detected by an air shower fluorescence or Cherenkov telescope. In this paper, we present detailed Monte Carlo simulations of corresponding event rates for the VERITAS and two proposed Cherenkov Telescope Array sites: Meteor Crater and Yavapai Ranch, which use representative AGN neutrino flux models and take into account topographic conditions of the detector sites. The calculated neutrino sensitivities depend on the observation time and the shape of the energy spectrum, but in some cases are comparable or even better than corresponding neutrino sensitivities of the IceCube detector. For VERITAS and the considered Cherenkov Telescope Array sites the expected neutrino sensitivities are up to factor 3 higher than for the MAGIC site because of the presence of surrounding mountains.


title: A comparison theorem for backward SPDEs with jumps

abstract: In this paper we obtain a comparison theorem for backward stochastic partial differential equation (SPDEs) with jumps. We apply it to introduce space-dependent convex risk measures as a model for risk in large systems of interacting components.


title: Irreducible many-body correlations in topologically ordered systems

abstract: Topologically ordered systems exhibit large-scale correlation in their ground states, which may be characterized by quantities such as topological entanglement entropy. We propose that the concept of irreducible many-body correlation, the correlation that cannot be implied by all local correlations, may also be used as a signature of topological order. In a topologically ordered system, we demonstrate that for a part of the system with holes, the reduced density matrix exhibits irreducible many-body correlation which becomes reducible when the holes are removed. The appearance of these irreducible correlations then represents a key feature of topological phase. We analyze the many-body correlation structures in the ground state of the toric code model in an external magnetic field, and show that the topological phase transition is signaled by the irreducible many-body correlations.


title: Precoding by Priority: A UEP Scheme for RaptorQ Codes

abstract: Raptor codes are the first class of fountain codes with linear time encoding and decoding. These codes are recommended in standards such as Third Generation Partnership Project (3GPP) and digital video broadcasting. RaptorQ codes are an extension to Raptor codes, having better coding efficiency and flexibility. Standard Raptor and RaptorQ codes are systematic with equal error protection of the data. However, in many applications such as MPEG transmission, there is a need for Unequal Error Protection (UEP): namely, some data symbols require higher error correction capabilities compared to others. We propose an approach that we call Priority Based Precode Ratio (PBPR) to achieve UEP for systematic RaptorQ and Raptor codes. Our UEP assumes that all symbols in a source block belong to the same importance class. The UEP is achieved by changing the number of precode symbols depending on the priority of the information symbols in the source block. PBPR provides UEP with the same number of decoding overhead symbols for source blocks with different importance classes. We demonstrate consistent improvements in the error correction capability of higher importance class compared to the lower importance class across the entire range of channel erasure probabilities. We also show that PBPR does not result in a significant increase in decoding and encoding times compared to the standard implementation.


title: Effects of Easy Hybrid Parallelization with CUDA for Numerical-Atomic-Orbital Density Functional Theory Calculation

abstract: We modified a MPI-friendly density functional theory (DFT) source code within hybrid parallelization including CUDA. Our objective is to find out how simple conversions within the hybrid parallelization with mid-range GPUs affect DFT code not originally suitable to CUDA. We settled several rules of hybrid parallelization for numerical-atomic-orbital (NAO) DFT codes. The test was performed on a magnetite material system with OpenMX code by utilizing a hardware system containing 2 Xeon E5606 CPUs and 2 Quadro 4000 GPUs. 3-way hybrid routines obtained a speedup of 7.55 while 2-way hybrid speedup by 10.94. GPUs with CUDA complement the efficiency of OpenMP and compensate CPUs' excessive competition within MPI.


title: Sensitivity relations for the Mayer problem with differential inclusions

abstract: In optimal control,sensitivity relations are usually understood as inclusions that identify the pair formed by the dual arc and the Hamiltonian, evaluated along the associated minimizing trajectory, as a suitable generalized gradient of the value function. In this paper, sensitivity relations are obtained for the Mayer problem associated with the differential inclusion $\dot x\in F(x)$ and applied to derive optimality conditions. Our first application concerns the maximum principle and consists in showing that a dual arc can be constructed for every element of the superdifferential of the final cost. As our second application, with every nonzero limiting gradient of the value function at some point $(t,x)$ we associate a family of optimal trajectories at $(t,x)$ with the property that families corresponding to distinct limiting gradients have empty intersection.


title: Quantum flag manifolds as quotients of degenerate quantized universal enveloping algebras

abstract: Let $\mathfrak{g}$ be a semi-simple Lie algebra with fixed root system, and $U_q(\mathfrak{g})$ the quantization of its universal enveloping algebra. Let $\mathcal{S}$ be a subset of the simple roots of $\mathfrak{g}$. We show that the defining relations for $U_q(\mathfrak{g})$ can be slightly modified in such a way that the resulting algebra $U_q(\mathfrak{g};\mathcal{S})$ allows a homomorphism onto (an extension of) the algebra $\mathrm{Pol}(\mathbb{G}_q/\mathbb{K}_{\mathcal{S},q})$ of functions on the quantum flag manifold $\mathbb{G}_q/\mathbb{K}_{\mathcal{S},q}$ corresponding to $\mathcal{S}$. Moreover, this homomorphism is equivariant with respect to a natural adjoint action of $U_q(\mathfrak{g})$ on $U_q(\mathfrak{g};\mathcal{S})$ and the standard action of $U_q(\mathfrak{g})$ on $Pol(\mathbb{G}_q/\mathbb{K}_{\mathcal{S},q})$.


title: Nonlinear response for external field and perturbation in the Vlasov system

abstract: A nonlinear response theory is provided by use of the transient linearization method in the spatially one-dimensional Vlasov systems. The theory inclusively gives responses to external fields and to perturbations for initial stationary states, and is applicable even to the critical point of a second order phase transition. We apply the theory to the Hamiltonian mean-field model, a toy model of a ferromagnetic body, and investigate the critical exponent associated with the response to the external field at the critical point in particular. The obtained critical exponent is nonclassical value 3/2, while the classical value is 3. However, interestingly, one scaling relation holds with another nonclassical critical exponent of susceptibility in the isolated Vlasov systems. Validity of the theory is numerically confirmed by directly simulating temporal evolutions of the Vlasov equation.


title: Possible Evidence for Metal Accretion onto the Surfaces of Metal-Poor Main-Sequence Stars

abstract: The entire evolution of the Milky Way, including its mass-assembly and star-formation history, is imprinted onto the chemo-dynamical distribution function of its member stars, f(x, v, [X/H]), in the multi-dimensional phase space spanned by position, velocity, and elemental abundance ratios. In particular, the chemo-dynamical distribution functions for low-mass stars (e.g., G- or K-type dwarfs) are precious tracers of the earliest stages of the Milky Way's formation, since their main-sequence lifetimes approach or exceed the age of the universe. A basic tenet of essentially all previous analyses is that the stellar metallicity, usually parametrized as [Fe/H], is conserved over time for main-sequence stars (at least those that have not been polluted due to mass transfer from binary companions). If this holds true, any correlations between metallicity and kinematics for long-lived main-sequence stars of different masses, effective temperatures, or spectral types must strictly be the same, since they reflect the same mass-assembly and star-formation histories. By analyzing a sample of nearby metal-poor halo and thick-disk stars on the main sequence, taken from Data Release 8 of the Sloan Digital Sky Survey, we find that the median metallicity of G-type dwarfs is systematically higher (by about 0.2 dex) than that of K-type dwarfs having the same median rotational velocity about the Galactic center. If it can be confirmed, this finding may invalidate the long-accepted assumption that the atmospheric metallicities of long-lived stars are conserved over time.


title: A Finite-Volume Method for Nonlinear Nonlocal Equations with a Gradient Flow Structure

abstract: We propose a positivity preserving entropy decreasing finite volume scheme for nonlinear nonlocal equations with a gradient flow structure. These properties allow for accurate computations of stationary states and long-time asymptotics demonstrated by suitably chosen test cases in which these features of the scheme are essential. The proposed scheme is able to cope with non-smooth stationary states, different time scales including metastability, as well as concentrations and self-similar behavior induced by singular nonlocal kernels. We use the scheme to explore properties of these equations beyond their present theoretical knowledge.