title: Optimal control of fast and high-fidelity quantum gates with electron and nuclear spins of a nitrogen-vacancy center in diamond

abstract: A negatively charged nitrogen vacancy (NV) center in diamond has been recognized as a good solid-state qubit. A system consisting of the electronic spin of the NV center and hyperfine-coupled nitrogen and additionally nearby carbon nuclear spins can form a quantum register of several qubits for quantum information processing or as a node in a quantum repeater. Several impressive experiments on the hybrid electron and nuclear spin register have been reported, but fidelities achieved so far are not yet at or below the thresholds required for fault-tolerant quantum computation (FTQC). Using quantum optimal control theory based on the Krotov method, we show here that fast and high-fidelity single-qubit and two-qubit gates in the universal quantum gate set for FTQC, taking into account the effects of the leakage state, nearby noise qubits and distant bath spins, can be achieved with errors less than those required by the threshold theorem of FTQC.


title: Optical Tuning of Exciton and Trion Emissions in Monolayer Phosphorene

abstract: Monolayer phosphorene provides a unique two-dimensional (2D) platform to investigate the fundamental dynamics of excitons and trions (charged excitons) in reduced dimensions. However, owing to its high instability, unambiguous identification of monolayer phosphorene has been elusive. Consequently, many important fundamental properties, such as exciton dynamics, remain underexplored. We report a rapid, noninvasive, and highly accurate approach based on optical interferometry to determine the layer number of phosphorene, and confirm the results with reliable photoluminescence measurements. Furthermore, we successfully probed the dynamics of excitons and trions in monolayer phosphorene by controlling the photo-carrier injection in a relatively low excitation power range. Based on our measured optical gap and the previously measured electronic energy gap, we determined the exciton binding energy to be ~0.3 eV for the monolayer phosphorene on SiO2/Si substrate, which agrees well with theoretical predictions. A huge trion binding energy of ~100 meV was first observed in monolayer phosphorene, which is around five times higher than that in transition metal dichalcogenide (TMD) monolayer semiconductor, such as MoS2. The carrier lifetime of exciton emission in monolayer phosphorene was measured to be ~220 ps, which is comparable to those in other 2D TMD semiconductors. Our results open new avenues for exploring fundamental phenomena and novel optoelectronic applications using monolayer phosphorene.


title: Throughput Optimal and Fast Near-Optimal Scheduling with Heterogeneously Delayed Network-State Information (Extended Version)

abstract: We consider the problem of distributed scheduling in wireless networks where heterogeneously delayed information about queue lengths and channel states of all links are available at all the transmitters. In an earlier work (by Reddy et al. in Queueing Systems, 2012), a throughput optimal scheduling policy (which we refer to henceforth as the R policy) for this setting was proposed. We study the R policy, and examine its two drawbacks -- (i) its huge computational complexity, and (ii) its non-optimal average per-packet queueing delay. We show that the R policy unnecessarily constrains itself to work with information that is more delayed than that afforded by the system. We propose a new policy that fully exploits the commonly available information, thereby greatly improving upon the computational complexity and the delay performance of the R policy. We show that our policy is throughput optimal. Our main contribution in this work is the design of two fast and near-throughput-optimal policies for this setting, whose explicit throughput and runtime performances we characterize analytically. While the R policy takes a few milliseconds to several tens of seconds to compute the schedule once (for varying number of links in the network), the running times of the proposed near-throughput-optimal algorithms range from a few microseconds to only a few hundred microseconds, and are thus suitable for practical implementation in networks with heterogeneously delayed information.


title: Thickness Measurements from Single X-ray Phase-contrast Speckle Projection

abstract: We propose a one-shot thickness measurement method for sponge-like structures using a propagation-based X-ray phase-contrast imaging (P-PCI) method. In P-PCI, the air-material interface refracts the incident X-ray. Refracted many times along their paths by such a structure, incident X-rays propagate randomly within a small divergent angle range, resulting in a speckle pattern in the captured image. We found structure thickness and contrast of a phase-contrast projection are directly related in images. This relationship can be described by a natural logarithm equation. Thus, from the one phase-contrast view, depth information can be retrieved from its contrast. Our preliminary biological experiments indicate promise in its application to measurements requiring in vivo and ongoing assessment of lung tumor progression.


title: New universality class in percolation on multifractal scale-free planar stochastic lattice

abstract: We investigate site percolation on a weighted planar stochastic lattice (WPSL) which is a multifractal and whose dual is a scale-free network. Percolation is typically characterized by percolation threshold $p_c$ and by a set of critical exponents $\beta$, $\gamma$, $\nu$ which describe the critical behavior of percolation probability $P(p)\sim (p_c-p)^\beta$, mean cluster size $S\sim (p_c-p)^{-\gamma}$ and the correlation length $\xi\sim (p_c-p)^{-\nu}$. Besides, the exponent $\tau$ characterizes the cluster size distribution function $n_s(p_c)\sim s^{-\tau}$ and the fractal dimension $d_f$ the spanning cluster. We obtain an exact value for $p_c$ and for all these exponents. Our results suggest that the percolation on WPSL belong to a new universality class as its exponents do not share the same value as for all the existing planar lattices.


title: Fundamental rate-loss tradeoff for optical quantum key distribution

abstract: Since 1984, various optical quantum key distribution (QKD) protocols have been proposed and examined. In all of them, the rate of secret key generation decays exponentially with distance. A natural and fundamental question is then whether there are yet-to-be discovered optical QKD protocols (without quantum repeaters) that could circumvent this rate-distance tradeoff. This paper provides a major step towards answering this question. We show that the secret-key-agreement capacity of a lossy and noisy optical channel assisted by unlimited two-way public classical communication is limited by an upper bound that is solely a function of the channel loss, regardless of how much optical power the protocol may use. Our result has major implications for understanding the secret-key-agreement capacity of optical channels---a long-standing open problem in optical quantum information theory---and strongly suggests a real need for quantum repeaters to perform QKD at high rates over long distances.


title: On the Stability of Online Language Features: How Much Text do you Need to know a Person?

abstract: In recent years, numerous studies have inferred personality and other traits from people's online writing. While these studies are encouraging, more information is needed in order to use these techniques with confidence. How do linguistic features vary across different online media, and how much text is required to have a representative sample for a person? In this paper, we examine several large sets of online, user-generated text, drawn from Twitter, email, blogs, and online discussion forums. We examine and compare population-wide results for the linguistic measure LIWC, and the inferred traits of Big5 Personality and Basic Human Values. We also empirically measure the stability of these traits across different sized samples for each individual. Our results highlight the importance of tuning models to each online medium, and include guidelines for the minimum amount of text required for a representative result.


title: Pure sea-quark contributions to the magnetic form factors of $\Sigma$ baryons

abstract: We propose the pure sea-quark contributions to the magnetic form factors of $\Sigma$ baryons, $G_{\Sigma^-}^u$ and $G_{\Sigma^+}^d$, as priority observables for the examination of sea-quark contributions to baryon structure, both in present lattice QCD simulations and possible future experimental measurement. $G_{\Sigma^-}^u$, the $u$-quark contribution to the magnetic form factor of $\Sigma^-$, and $G_{\Sigma^+}^d$, the $d$-quark contribution to the magnetic form factor of $\Sigma^+$, are similar to the strange quark contribution to the magnetic form factor of the nucleon, but promise to be larger by an order of magnitude. We explore the size of this quantity within chiral effective field theory, including both octet and decuplet intermediate states. The finite range regularization approach is applied to deal with ultraviolet divergences. Drawing on an established connection between quenched and full QCD, this approach makes it possible to predict the sea quark contribution to the magnetic form factor purely from the meson loop. In the familiar convention where the quark charge is set to unity $G_{\Sigma^-}^u = G_{\Sigma^+}^d$. We find a value of $-0.38^{+0.16}_{-0.17}\ \mu_N$, which is about seven times larger than the strange magnetic moment of the nucleon found in the same approach. Including quark charge factors, the $u$-quark contribution to the $\Sigma^-$ magnetic moment exceeds the strange quark contribution to the nucleon magnetic moment by a factor of 14.


title: Tuning the electronic structure of Sr2IrO4 thin films by bulk electronic doping using molecular beam epitaxy

abstract: By means of oxide molecular beam epitaxy with shutter-growth mode, we have fabricated a series of electron-doped (Sr1-xLax)2IrO4(001)(x = 0, 0.05, 0.1 and 0.15) single crystalline thin films and then investigated the doping dependence of electronic structure utilizing in-situ angle-resolved photoemission spectroscopy. We find that with increasing doping proportion, the Fermi levels of samples progressively shift upward. Prominently, an extra electron pocket crossing the Fermi level around the M point has been evidently observed in 15 % nominal doping sample. Moreover, bulk-sensitive transport measurements confirm that doping effectively suppresses the insulating state with respect to the as-grown Sr2IrO4, though doped samples still remain insulating at low temperatures due to the localization effect possibly stemming from disorders including oxygen deficiencies. Our work provides another feasible doping method to tune electronic structure of Sr2IrO4.


title: Social Trust Prediction via Max-norm Constrained 1-bit Matrix Completion

abstract: Social trust prediction addresses the significant problem of exploring interactions among users in social networks. Naturally, this problem can be formulated in the matrix completion framework, with each entry indicating the trustness or distrustness. However, there are two challenges for the social trust problem: 1) the observed data are with sign (1-bit) measurements; 2) they are typically sampled non-uniformly. Most of the previous matrix completion methods do not well handle the two issues. Motivated by the recent progress of max-norm, we propose to solve the problem with a 1-bit max-norm constrained formulation. Since max-norm is not easy to optimize, we utilize a reformulation of max-norm which facilitates an efficient projected gradient decent algorithm. We demonstrate the superiority of our formulation on two benchmark datasets.


title: Resource Allocation with Reverse Pricing for Communication Networks

abstract: Reverse pricing has been recognized as an effective tool to handle demand uncertainty in the travel industry (e.g., airlines and hotels). To investigate its viability for communication networks, we study the practical limitations of (operator-driven) time-dependent pricing that has been recently introduced, taking into account demand uncertainty. Compared to (operator-driven) time-dependent pricing, we show that the proposed pricing scheme can achieve "triple-win" solutions: an increase in the total average revenue of the operator; higher average resource utilization efficiency; and an increment in the total average payoff of the users. Our findings provide a new outlook on resource allocation, and design guidelines for adopting the reverse pricing scheme.


title: Periodicity for the Hadamard walk on cycles

abstract: The present paper treats the period T_N of the Hadamard walk on a cycle C_N with N vertices. Dukes (2014) considered the periodicity of more general quantum walks on C_N and showed T_2 =2, T_4=8, T_8=24 for the Hadamard walk case. We prove that the Hadamard walk does not have any period except for his case, i.e., N=2, 4, 8. Our method is based on a path counting and cyclotomic polynomials which is different from his approach based on the property of eigenvalues for unitary matrix that determines the evolution of the walk.


title: Testing the performance of technical trading rules in the Chinese market

abstract: Technical trading rules have a long history of being used by practitioners in financial markets. Their profitable ability and efficiency of technical trading rules are yet controversial. In this paper, we test the performance of more than seven thousands traditional technical trading rules on the Shanghai Securities Composite Index (SSCI) from May 21, 1992 through June 30, 2013 and Shanghai Shenzhen 300 Index (SHSZ 300) from April 8, 2005 through June 30, 2013 to check whether an effective trading strategy could be found by using the performance measurements based on the return and Sharpe ratio. To correct for the influence of the data-snooping effect, we adopt the Superior Predictive Ability test to evaluate if there exists a trading rule that can significantly outperform the benchmark. The result shows that for SSCI, technical trading rules offer significant profitability, while for SHSZ 300, this ability is lost. We further partition the SSCI into two sub-series and find that the efficiency of technical trading in sub-series, which have exactly the same spanning period as that of SHSZ 300, is severely weakened. By testing the trading rules on both indexes with a five-year moving window, we find that the financial bubble from 2005 to 2007 greatly improve the effectiveness of technical trading rules. This is consistent with the predictive ability of technical trading rules which appears when the market is less efficient.


title: Erdos-Feller-Kolmogorov-Petrowsky law of the iterated logarithm for self-normalized martingales: a game-theoretic approach

abstract: We prove an Erdos-Feller-Kolmogorov-Petrowsky law of the iterated logarithm for self-normalized martingales. Our proof is given in the framework of the game-theoretic probability of Shafer and Vovk. As many other game-theoretic proofs, our proof is self-contained and explicit.


title: Key Science Observations of AGNs with KaVA Array

abstract: KaVA (KVN and VERA Array) is a new combined VLBI array with KVN (Korean VLBI Network) and VERA (VLBI Exploration of Radio Astrometry). First, we briefly review the imaging capabilities of KaVA array which actually achieves more than three times better dynamic range than that achieved by VERA alone. The KaVA images clearly show detailed structures of extended radio jets in AGNs. Next, we represent the key science program to be led by KaVA AGN sub working group. We will conduct the monitoring observations of Sgr A* and M87 because of the largeness of their central super-massive black hole angular sizes. The main science goals of the program are (i) testing magnetically-driven-jet paradigm by mapping velocity fields of the M87 jet, and (ii) obtaining tight constraints on physical properties of radio emitting region in Sgr A*.


title: Thin and Thick Strip Passage Times for L\'evy Flights and L\'evy Processes

abstract: We review some of the theory relevant to passage times of one-dimensional L\'evy processes out of bounded regions, highlighting results that are useful in physical phenomena modelled by heavy-tailed L\'evy flights. The process is hypothesised to describe the motion of a particle on the line, starting at $0$, and exiting either a fixed interval $[-r, r]$, $r > 0$, or a time-dependent, expanding, set of intervals of the form $[-r t^\kappa, r t^\kappa]$, $r > 0$, $\kappa > 0$. Asymptotic behaviour of the exit time may be as $r \downarrow 0$ or as $r \to \infty$, but particular emphasis is placed herein on "small time" approximations, corresponding to exits from or transmissions through thin strips. Applications occur for example in the transmission of photons through moderately doped thin or thick wafers by means of "photon recycling", and in atmospheric radiation modelling.


title: Seeking instructional specificity: an example from analogical instruction

abstract: Broad instructional methods like interactive engagement have been shown to be effective, but such general characterization provides little guidance on the details of how to structure the instructional materials. In this study, we seek instructional specificity by comparing two ways of using an analogy to learn a target physical principle: (i) applying the analogy to the target physical domain on a Case-by-Case basis and (ii) using the analogy to create a General Rule in the target physical domain. In the discussion sections of a large, introductory physics course (N = 231), students who sought a General Rule were better able to discover and apply a correct physics principle than students who analyzed the examples Case-by-Case. The difference persisted at a reduced level after subsequent direct instruction. We argue that students who performed Case-by-Case analyses are more likely to focus on idiosyncratic problem-specific features rather than the deep structural features. This study provides an example of investigating how the specific structure of instructional materials can be consequential for what is learned.


title: Shape optimization for surface functionals in Navier--Stokes flow using a phase field approach

abstract: We consider shape and topology optimization for fluids which are governed by the Navier--Stokes equations. Shapes are modelled with the help of a phase field approach and the solid body is relaxed to be a porous medium. The phase field method uses a Ginzburg--Landau functional in order to approximate a perimeter penalization. We focus on surface functionals and carefully introduce a new modelling variant, show existence of minimizers and derive first order necessary conditions. These conditions are related to classical shape derivatives by identifying the sharp interface limit with the help of formally matched asymptotic expansions. Finally, we present numerical computations based on a Cahn--Hilliard type gradient descent which demonstrate that the method can be used to solve shape optimization problems for fluids with the help of the new approach.


title: Driven Anomalous Diffusion: An example from Polymer Stretching

abstract: The way tension propagates along a chain is a key to govern many of anomalous dynamics in macromolecular systems. After introducing the weak and the strong force regimes of the tension propagation, we focus on the latter, in which the dynamical fluctuations of a segment in a long polymer during its stretching process is investigated. We show that the response, i.e., average drift, is anomalous, which is characterized by the nonlinear memory kernel, and its relation to the fluctuation is nontrivial. These features are discussed on the basis of the generalized Langevin equation, in which the role of the temporal change in spring constant due to the stress hardening is pinpointed. We carried out the molecular dynamics simulation, which supports our theory.


title: Universal coarsening dynamics of a quenched ferromagnetic spin-1 condensate

abstract: We demonstrate that a quasi-two-dimensional spin-1 condensate quenched to a ferromagnetic phase undergoes universal coarsening in its late time dynamics. The quench can be implemented by a sudden change in the applied magnetic field and, depending on the final value, the ferromagnetic phase has easy-axis (Ising) or easy-plane (XY) symmetry, with different dynamical critical exponents. Our results for the easy-plane phase reveal a fractal domain structure and the crucial role of polar-core spin vortices in the coarsening dynamics.


title: Pump Electron-Positron Pairs from Potential Well

abstract: In this paper we show that electron-positron pairs can be pumped inexhaustibly with a constant production rate from the one dimensional potential well with oscillating depth or width. Bound states embedded in the the Dirac sea can be pulled out and pushed to the positive continuum, and become scattering states. Pauli block, which dominant the saturation of pair creation in the static super-critical potential well, can be broken by the ejection of electrons. We find that the width oscillating mode is more efficient that the depth oscillating mode. In the adiabatic limit, pair number as a function of upper boundary of the oscillating, will reveal the diving of the bound states.


title: Twisted cubics on singular cubic fourfolds - On Starr's fibration

abstract: We study lines and twisted cubics on cubic fourfolds with simple isolated singularities. We show that the Hilbert scheme compactification of the total space of Starr's fibration on the space of twisted cubics on a cubic fourfold with simple isolated singularities and not containing a plane admits a contraction to a singular projective symplectic variety of dimension eight. The latter admits a crepant resolution which is a deformation of the symplectic eightfold constructed from the space of twisted cubics on a smooth cubic fourfold. Thereby we give another proof of the fact that the latter is a deformation of the Hilbert scheme of four points on a K3 surface. We also show that the variety of lines on a cubic fourfold with simple singularities is a singular symplectic variety birational to the Hilbert scheme of two points on a K3 surface.


title: Thermodynamics of Error Correction

abstract: Information processing at the molecular scale is limited by thermal fluctuations. This can cause undesired consequences in copying information since thermal noise can lead to errors that can compromise the functionality of the copy. For example, a high error rate during DNA duplication can lead to cell death. Given the importance of accurate copying at the molecular scale, it is fundamental to understand its thermodynamic features. In this paper, we derive a universal expression for the copy error as a function of entropy production and {\cred work dissipated by the system during wrong incorporations}. Its derivation is based on the second law of thermodynamics, hence its validity is independent of the details of the molecular machinery, be it any polymerase or artificial copying device. Using this expression, we find that information can be copied in three different regimes. In two of them, work is dissipated to either increase or decrease the error. In the third regime, the protocol extracts work while correcting errors, reminiscent of a Maxwell demon. As a case study, we apply our framework to study a copy protocol assisted by kinetic proofreading, and show that it can operate in any of these three regimes. We finally show that, for any effective proofreading scheme, error reduction is limited by the chemical driving of the proofreading reaction.


title: Location of eigenvalues for the wave equation with dissipative boundary conditions

abstract: We examine the location of the eigenvalues of the generator $G$ of a semi-group $V(t) = e^{tG},\: t \geq 0,$ related to the wave equation in an unbounded domain $\Omega \subset {\mathbb R}^d$ with dissipative boundary condition $\partial_{\nu}u - \gamma(x) \partial_t u = 0$ on $\Gamma = \partial \Omega.$ We study two cases: $(A): \: 0 < \gamma(x) < 1,\: \forall x \in \Gamma$ and $(B):\: 1 < \gamma(x), \: \forall x \in \Gamma.$ We prove that for every $0 < \epsilon \ll 1,$ the eigenvalues of $G$ in the case $(A)$ lie in the region $\Lambda_{\epsilon} = \{z \in {\mathbb C}:\: |\Re z | \leq C_{\epsilon} (|\Im z|^{\frac{1}{2} + \epsilon} + 1), \: \Re z < 0\},$ while in the case $(B)$ for every $0 < \epsilon \ll 1$ and every $N \in {\mathbb N}$ the eigenvalues lie in $\Lambda_{\epsilon} \cup {\mathcal R}_N,$ where ${\mathcal R}_N = \{z \in {\mathbb C}:\: |\Im z| \leq C_N (|\Re z| + 1)^{-N},\: \Re z < 0\}.$


title: Satisfiability of Modal Inclusion Logic: Lax and Strict Semantics

abstract: We investigate the computational complexity of the satisfiability problem of modal inclusion logic. We distinguish two variants of the problem: one for the strict and another one for the lax semantics. Both problems turn out to be EXPTIME-complete on general structures. Finally, we show how for a specific class of structures NEXPTIME-completeness for these problems under strict semantics can be achieved.


title: On the paper "statistical approximation by positive linear operators"

abstract: We give an example to show that the main result of [1] is incorrect.


title: Non-equilibrium quasi-long-range order of driven random field O(N) model

abstract: We investigate three-dimensional O(N) spin models driven with a uniform velocity over a random field. Within a spin-wave approximation, it is shown that in the strong driving regime the model with N=2 exhibits a quasi-long-range order in which the spatial correlation function decays in a power-law form. Furthermore, for the cases that N=2 and 3, we numerically demonstrate a non-equilibrium phase transition between the quasi-long-range order phase and the disordered phase, which turns out to resemble the Kosterlitz-Thouless transition in the two-dimensional pure XY model in equilibrium.


title: Giant anomalous self-steepening in photonic crystal waveguides

abstract: Self-steepening of optical pulses arises due the dispersive contribution of the $\chi^{(3)}(\omega)$ Kerr nonlinearity. In typical structures this response is on the order of a few femtoseconds with a fixed frequency response. In contrast, the effective $\chi^{(3)}$ Kerr nonlinearity in photonic crystal waveguides (PhCWGs) is largely determined by the geometrical parameters of the structure and is consequently tunable over a wide range. Here we show self-steepening based on group-velocity (group-index) modulation for the first time, giving rise to a new physical mechanism for generating this effect. Further, we demonstrate that periodic media such as PhCWGS can exhibit self-steepening coefficients two orders of magnitude larger than typical systems. At these magnitudes the self-steepening strongly affects the nonlinear pulse dynamics even for picosecond pulses. Due to interaction with additional higher-order nonlinearities in the semiconductor materials under consideration, we employ a generalized nonlinear Schr\"{o}dinger equation numerical model to describe the impact of self-steepening on the temporal and spectral properties of the optical pulses in practical systems, including new figures of merit. These results provide a theoretical description for recent experimental results presented in [Scientific Reports 3, 1100 (2013) and Phys. Rev. A 87, 041802 (2013)]. More generally, these observations apply to all periodic media due to the rapid group-velocity variation characteristic of these structures.


title: Signal-Flow Based Runge-Kutta Methods for the Simulation of Complex Networks

abstract: Complex dynamical networks appear in a wide range of physical, biological, and engineering systems. The coupling of subsystems with varying time scales often results in multirate behavior. During the simulation of highly integrated circuits, for example, only a few elements underlie changing signals whereas the major part -- usually up to 80 or even 90 per cent -- remains latent. Standard integration schemes discretize the entire circuit with a single step size which is mainly limited by the accuracy requirements of the rapidly changing subcircuits. It is of a particular interest to speed up the simulation without a significant loss of accuracy. By exploiting the latency of the system, only a fraction of the equations has to be formulated and solved at a given time point. G\"unther and Rentrop suggest that multirate strategies must be based both on the numerical information of the integration scheme and on the topology of the circuit. In this paper, we will introduce a directed graph describing the interdependency of the underlying system and propose Runge--Kutta methods which utilize the signal flow of the system in order to identify and exploit inactive regions. Furthermore, we describe an extension of these methods to identify and exploit periodic subsystems.


title: Star junctions and watermelons of pure or random quantum Ising chains : finite-size properties of the energy gap at criticality

abstract: We consider $M \geq 2$ pure or random quantum Ising chains of $N$ spins when they are coupled via a single star junction at their origins or when they are coupled via two star junctions at the their two ends leading to the watermelon geometry. The energy gap is studied via a sequential self-dual real-space renormalization procedure that can be explicitly solved in terms of Kesten variables containing the initial couplings and and the initial transverse fields. In the pure case at criticality, the gap is found to decay as a power-law $\Delta_M \propto N^{-z(M)} $ with the dynamical exponent $z(M)=\frac{M}{2}$ for the single star junction (the case $M=2$ corresponds to $z=1$ for a single chain with free boundary conditions) and $z(M)=M-1$ for the watermelon (the case $M=2$ corresponds to $z=1$ for a single chain with periodic boundary conditions). In the random case at criticality, the gap follows the Infinite Disorder Fixed Point scaling $\ln \Delta_M = -N^{\psi} g$ with the same activated exponent $\psi=\frac{1}{2}$ as the single chain corresponding to $M=2$, and where $g$ is an $O(1)$ random positive variable, whose distribution depends upon the number $M$ of chains and upon the geometry (star or watermelon).


title: Covariant mutually unbiased bases

abstract: The connection between maximal sets of mutually unbiased bases (MUBs) in a prime-power dimensional Hilbert space and finite phase-space geometries is well known. In this article we classify MUBs according to their degree of covariance with respect to the natural symmetries of a finite phase-space, which are the group of its affine symplectic transformations. We prove that there exist maximal sets of MUBs that are covariant with respect to the full group only in odd prime-power dimensional spaces, and in this case their equivalence class is actually unique. Despite this limitation, we show that in even-prime power dimension covariance can still be achieved by restricting to proper subgroups of the symplectic group, that constitute the finite analogues of the oscillator group. For these subgroups, we explicitly construct the unitary operators yielding the covariance.


title: On the Minimum Distance of Array-Based Spatially-Coupled Low-Density Parity-Check Codes

abstract: An array low-density parity-check (LDPC) code is a quasi-cyclic LDPC code specified by two integers $q$ and $m$, where $q$ is an odd prime and $m \leq q$. The exact minimum distance, for small $q$ and $m$, has been calculated, and tight upper bounds on it for $m \leq 7$ have been derived. In this work, we study the minimum distance of the spatially-coupled version of these codes. In particular, several tight upper bounds on the optimal minimum distance for coupling length at least two and $m=3,4,5$, that are independent of $q$ and that are valid for all values of $q \geq q_0$ where $q_0$ depends on $m$, are presented. Furthermore, we show by exhaustive search that by carefully selecting the edge spreading or unwrapping procedure, the minimum distance (when $q$ is not very large) can be significantly increased, especially for $m=5$.


title: A NLTE analysis of the hot subdwarf O star BD+28 4211. II. The optical spectrum

abstract: We present the second part of our detailed analysis of the hot sdO and spectroscopic standard star BD+28 4211, in which we focus on the optical spectrum. In the first part of our study, we determined the abundances of some 11 metals detected in the atmosphere of BD+28 4211 using UV spectra of the star and corroborated the fundamental parameters estimated in past studies (Teff $\sim$ 82,000 K, log g $\sim$ 6.2, and solar N(He)/N(H)). In this work, we aim at rederiving these secured parameters on the sole basis of high-quality optical spectra. A first grid of non-LTE line-blanketed model atmospheres, including metals with the abundances derived from the UV spectrum, does not give satisfactory results when we apply a standard simultaneous fitting procedure to the observed H and He lines of our optical spectra. The line profiles are not finely reproduced and the resulting effective temperatures, in particular, are too low by $\sim$10,000 K. We next investigate the probable cause of this failure, that is, the importance of missing opacity sources on the atmospheric stratification. We compare line profiles computed from models with artificially boosted metallicities, from solar abundances to 15$\times$ these values. We find that the structural effects saturate for a metallicity of $\sim$10x solar, and use this to compute a second full grid of models and synthetic spectra. This metal-enriched grid allows us to achieve significantly improved spectral fits with models having the expected parameters. Our test case thus reveals that there is still a need for models with enhanced metallicity for better estimating the atmospheric parameters of objects such as hot subdwarfs and hot white dwarfs if only optical spectra are available.


title: Optimal convergence of adaptive FEM for eigenvalue clusters in mixed form

abstract: It is shown that the h-adaptive mixed finite element method for the discretization of eigenvalue clusters of the Laplace operator produces optimal convergence rates in terms of nonlinear approximation classes. The results are valid for the typical mixed spaces of Raviart-Thomas or Brezzi-Douglas-Marini type with arbitrary fixed polynomial degree in two and three space dimensions.


title: Highly Nonlinear Luminescence Induced by Gold Nanoparticles on Glass Surfaces with Continuous-Wave Laser Illumination

abstract: We report on highly nonlinear luminescence being observed from individual spherical gold nanoparticles immobilized on a glass surface and illuminated by continuous-wave (CW) lasers with relatively low power. The nonlinear luminescence shows optical super-resolution beyond the diffraction limit in three dimensions compared to the scatting of the excitation laser light. The luminescence intensity from most nanoparticles is proportional to the 5th--7th power of the excitation laser power and has wide excitation and emission spectra across the visible wavelength range. Strong nonlinear luminescence is only observed near the glass surface. High optical nonlinearity excited by low CW laser power is related to a long-lived dark state of the gold nanoparticles, where the excitation light is strongly absorbed. This phenomenon has potential biological applications in super-resolution and deep tissue imaging.


title: Interacting multiple zero mode formulation and its application to a system consisting of a dark soliton in a condensate

abstract: To formulate the zero modes in a finite-size system with spontaneous breakdown of symmetries in quantum field theory is not trivial, for in the naive Bogoliubov theory, one encounters difficulties such as phase diffusion, the absence of a definite criterion for determining the ground state, and infrared divergences. A new interacting zero mode formulation that has been proposed for systems with a single zero mode to avoid these difficulties is extended to general systems with multiple zero modes. It naturally and definitely gives the interactions among the quantized zero modes, the consequences of which can be observed experimentally. In this paper, as a typical example, we consider an atomic Bose-Einstein condensed system with a dark soliton that contains two zero modes corresponding to spontaneous breakdown of the U(1) gauge and translational symmetries. Then we evaluate the standard deviations of the zero mode operators and see how the mutual interaction between the two zero modes affects them.


title: Two-loop stability of singlet extensions of the SM with dark matter

abstract: We present a study of the high energy stability of a minimal complex singlet extension of the Standard Model with or without dark matter (CxSM). We start by obtaining the beta functions of the couplings of the theory from the effective potential and then perform the RGE evolution for the allowed parameter space of the model at the electroweak scale. We obtain the scale up to which the model survives and combine this information with all the LHC measurements as well as bounds from dark matter detection experiments as well as the relic density best measurement from Planck. We conclude that scenarios which solve both the dark matter and stability problems must contain a dark particle heavier than 50 GeV and a new visible state (mixing with the SM-like Higgs) heavier than 170 GeV.


title: Casimir-Polder effect with thermally excited surfaces

abstract: We take a closer look at the fundamental Casimir-Polder interaction between quantum particles and dispersive dielectric surfaces with surface polariton or plasmon resonances. Linear response theory shows that in the near field, van der Waals, regime the free energy shift of a particle contains a thermal component that depends exclusively on the population/excitation of the evanescent surface polariton/plasmon modes. Our work makes evident the link between particle surface interaction and near field thermal emission and demonstrates how this can be used to engineer Casimir-Polder forces. We also examine how the exotic effects of surface waves are washed out as the distance from the surface increases. In the case of molecules or excited state atoms, far field approximations result in a classical dipole-dipole interaction which depends on the surface reflectivity and the mean number of photons at the frequency of the atomic/molecular transition. Finally we present numerical results for the CP interaction between Cs atoms and various dielectric surfaces with a single polariton resonance and discuss the implications of temperature and retardation effects for specific spectroscopic experiments.


title: Information Gathering in Networks via Active Exploration

abstract: How should we gather information in a network, where each node's visibility is limited to its local neighborhood? This problem arises in numerous real-world applications, such as surveying and task routing in social networks, team formation in collaborative networks and experimental design with dependency constraints. Often the informativeness of a set of nodes can be quantified via a submodular utility function. Existing approaches for submodular optimization, however, require that the set of all nodes that can be selected is known ahead of time, which is often unrealistic. In contrast, we propose a novel model where we start our exploration from an initial node, and new nodes become visible and available for selection only once one of their neighbors has been chosen. We then present a general algorithm NetExp for this problem, and provide theoretical bounds on its performance dependent on structural properties of the underlying network. We evaluate our methodology on various simulated problem instances as well as on data collected from social question answering system deployed within a large enterprise.


title: Constraints on the Georgi-Machacek Model by Current LHC Data

abstract: In the Georgi-Machacek model with a custodial symmetry in the Higgs potential and vacuum alignment, the triplet vacuum expectation value is allowed to be of ${\cal O}(10)$ GeV, which leads to the possibility of significant modifications in the couplings of the SM-like Higgs bosons $h$ with other SM particles. In this talk given at the HPNP2015 conference held at Toyama University, we review constraints on the model based on the latest LHC data of the SM-like Higgs boson, like-sign $W$ boson events, and searches for additional neutral Higgs bosons. In particular, we concentrate on the parameter space for small mixing angle $\alpha$ between the two custodial singlets. It is pointed out that constraints from the non-SM custodial singlet are most constraining and those from the 5-plet are independent of $\alpha$. While currently there is no constraint from the 3-plet, we show that its $f \bar f$ and $\gamma\gamma$ channels through the gluon fusion production can be very promising for searches or constraining in the mass range between 160 GeV and 350 GeV because of its gauge-phobic property.


title: Symplectic 4-manifolds via Lorentzian geometry

abstract: We observe that, in dimension four, symplectic forms may be obtained via Lorentzian geometry; in particular, null vector fields can give rise to exact symplectic forms. That a null vector field is nowhere vanishing yet orthogonal to itself is essential to this construction. Specifically, we show that on a Lorentzian 4-manifold $(M,g)$, if ${\boldsymbol k}$ is a complete null vector field with geodesic flow along which $\text{Ric}({\boldsymbol k},{\boldsymbol k}) > 0$, and if $f$ is any smooth function on $M$ with ${\boldsymbol k}(f)$ nowhere vanishing, then $dg(e^f{\boldsymbol k},\cdot)$ is a symplectic form and ${\boldsymbol k}/{\boldsymbol k}(f)$ is a Liouville vector field; any null surface to which ${\boldsymbol k}$ is tangent is then a Lagrangian submanifold. Even if the Ricci curvature condition is not satisfied, one can still construct such symplectic forms with additional information from ${\boldsymbol k}$; we give an example of this, with ${\boldsymbol k}$ a complete Liouville vector field, on the maximally extended "rapidly rotating" Kerr spacetime.


title: From the granular Leidenfrost state to buoyancy-driven convection

abstract: Grains inside a vertically vibrated box undergo a transition from a density inverted and horizontally homogeneous state, referred to as the granular Leidenfrost state, to a buoyancy-driven convective state. We perform a simulational study of the precursors of such a transition, and quantify their dynamics as the bed of grains is progressively fluidized. The transition is preceded by transient convective states, which increase their correlation time as the transition point is approached. Increasingly correlated convective flows lead to density fluctuations, as quantified by the structure factor, that also shows critical behaviour near the transition point. The amplitude of the modulations in the vertical velocity field are seen to be best described by a quintic supercritical amplitude equation with an additive noise term. The validity of such an amplitude equation, and previously observed collective semi-periodic oscillations of the bed of grains, suggests a new interpretation of the transition analogous to a coupled chain of vertically vibrated damped oscillators. Increasing the size of the container shows metastability of convective states, as well as an overall invariant critical behaviour close to the transition.


title: Redshift and spatial distribution of the intermediate gamma-ray bursts

abstract: One of the most important task of the Gamma-Ray Burst field is the classification of the bursts. Many researches have proven the existence of the third kind (intermediate duration) of GRBs in the BATSE data. Recent works have analyzed BeppoSax and Swift observations and can also identify three types of GRBs in the data sets. However, the class memberships are probabilistic we have enough observed redshifts to calculate the redshift and spatial distribution of the intermediate GRBs. They are significantly farther than the short bursts and seems to be closer than the long ones.


title: Variability of young stellar objects: accretion, disks, outflows and magnetic activity

abstract: This article represents a short review of the variability characteristics of young stellar objects. Variability is a key property of young stars. Two major origins may be distinguished: a scaled-up version of the magnetic activity seen on main-sequence stars and various processes related to circumstellar disks, accretion and outflows.


title: Methyl Acetate and its singly deuterated isotopomers in the interstellar medium

abstract: Methyl acetate (CH_3COOCH_3) has been recently observed by IRAM 30 m radio telescope in Orion though the presence of its deuterated isotopomers is yet to be confirmed. We therefore study the properties of various forms of methyl acetate, namely, CH_3COOCH_3, CH_2DCOOCH_3 and CH_3COOCH_2D. Our simulation reveals that these species could be produced efficiently both in gas as well as in ice phases. Production of methyl acetate could follow radical-radical reaction between acetyl (CH_3CO) and methoxy (CH_3O) radicals. To predict abundances of CH_3COOCH_3 along with its two singly deuterated isotopomers and its two isomers (ethyl formate and hydroxyacetone), we prepare a gas-grain chemical network to study chemical evolution of these molecules. Since gas phase rate coefficients for methyl acetate and its related species were unknown, either we consider similar rate coefficients for similar types of reactions (by following existing data bases) or we carry out quantum chemical calculations to estimate the unknown rate coefficients. For the surface reactions, we use adsorption energies of reactants from some earlier studies. Moreover, we perform quantum chemical calculations to obtain spectral properties of methyl acetate in infrared and sub-millimeter regions. We prepare two catalog files for the rotational transitions of CH_2DCOOCH_3 and CH_3COOCH_2D in JPL format, which could be useful for their detection in regions of interstellar media where CH_3COOCH_3 has already been observed.


title: Entropy production and detailed balance for a class of quantum Markov semigroups

abstract: We give an explicit entropy production formula for a class of quantum Markov semigroups, arising in the weak coupling limit of a system coupled with reservoirs, whose generators $\mathcal{L}$ are sums of other generators $\mathcal{L}_\omega$ associated with positive Bohr frequencies $\omega$ of the system. As a consequence, we show that any such semigroup satisfies the quantum detailed balance condition with respect to an invariant state if and only if all semigroups generated by each $\mathcal{L}_\omega$ so do with respect to the same invariant state.


title: Generic Multiuser Coordinated Beamforming for Underlay Spectrum Sharing

abstract: The beamforming techniques have been recently studied as possible enablers for underlay spectrum sharing. The existing beamforming techniques have several common limitations: they are usually system model specific, cannot operate with arbitrary number of transmit/receive antennas, and cannot serve arbitrary number of users. Moreover, the beamforming techniques for underlay spectrum sharing do not consider the interference originating from the incumbent primary system. This work extends the common underlay sharing model by incorporating the interference originating from the incumbent system into generic combined beamforming design that can be applied on interference, broadcast or multiple access channels. The paper proposes two novel multiuser beamforming algorithms for user fairness and sum rate maximization, utilizing newly derived convex optimization problems for transmit and receive beamformers calculation in a recursive optimization. Both beamforming algorithms provide efficient operation for the interference, broadcast and multiple access channels, as well as for arbitrary number of antennas and secondary users in the system. Furthermore, the paper proposes a successive transmit/receive optimization approach that reduces the computational complexity of the proposed recursive algorithms. The results show that the proposed complexity reduction significantly improves the convergence rates and can facilitate their operation in scenarios which require agile beamformers computation.


title: Enhancement of the $H^\pm W^\mp Z$ vertex in the three scalar doublet model

abstract: We compute one-loop induced trilinear vertices with physical charged Higgs bosons $H^\pm$ and ordinary gauge bosons, i.e., $H^\pm W^\mp Z$ and $H^\pm W^\mp \gamma$, in the model with two active plus one inert scalar doublet fields under a $Z_2(\text{unbroken})\times \tilde{Z}_2(\text{softly-broken})$ symmetry. The $Z_2$ and $\tilde{Z}_2$ symmetries are introduced to guarantee the stability of a dark matter candidate and to forbid the flavour changing neutral current at the tree level, respectively. The dominant form factor $F_Z$ of the $H^\pm W^\mp Z$ vertex can be enhanced by non-decoupling effects of extra scalar boson loop contributions. We find that, in such a model, $|F_Z|^2$ can be one order of magnitude larger than that predicted in two Higgs doublet models under the constraints from vacuum stability, perturbative unitarity and the electroweak precision observables. In addition, the branching fraction of the $H^\pm \to W^\pm Z$ $(H^\pm \to W^\pm \gamma)$ mode can be of order 10~(1)\% level when the mass of $H^\pm$ is below the top quark mass. Such a light $H^\pm$ is allowed by the so-called Type-I and Type-X Yukawa interactions which appear under the classification of the $\tilde{Z}$ charge assignment of the quarks and leptons. We also calculate the cross sections for the processes $H^\pm \to W^\pm Z$ and $H^\pm \to W^\pm \gamma$ onset by the top quark decay $t\to H^\pm b$ and electroweak $H^\pm$ production at the LHC.


title: Processes iterated ad libitum

abstract: Consider the $n$th iterated Brownian motion $I^{(n)}=B_n \circ\cdots \circ B_1$. Curien and Konstantopoulos proved that for any distinct numbers $t_i\neq 0$, $(I^{(n)}(t_1),\dots,I^{(n)}(t_k))$ converges in distribution to a limit $I[k]$ independent of the $t_i$'s, exchangeable, and gave some elements on the limit occupation measure of $I^{(n)}$. Here, we prove under some conditions, finite dimensional distributions of $n$th iterated two-sided stable processes converge, and the same holds the reflected Brownian motions. We give a description of the law of $I[k]$, of the finite dimensional distributions of $I^{(n)}$, as well as those of the iterated reflected Brownian motion iterated ad libitum.


title: Situational Object Boundary Detection

abstract: Intuitively, the appearance of true object boundaries varies from image to image. Hence the usual monolithic approach of training a single boundary predictor and applying it to all images regardless of their content is bound to be suboptimal. In this paper we therefore propose situational object boundary detection: We first define a variety of situations and train a specialized object boundary detector for each of them using [Dollar and Zitnick 2013]. Then given a test image, we classify it into these situations using its context, which we model by global image appearance. We apply the corresponding situational object boundary detectors, and fuse them based on the classification probabilities. In experiments on ImageNet, Microsoft COCO, and Pascal VOC 2012 segmentation we show that our situational object boundary detection gives significant improvements over a monolithic approach. Additionally, our method substantially outperforms [Hariharan et al. 2011] on semantic contour detection on their SBD dataset.


title: The three different regimes in coulombic friction

abstract: de Gennes identified three regimes in the phenomenon of the Langevin equation wich includes Coulombic friction. Here we extend and precise this phenomenon to a constant external force.


title: Second Harmonic Coherent Driving of a Spin Qubit in a Si/SiGe Quantum Dot

abstract: We demonstrate coherent driving of a single electron spin using second harmonic excitation in a Si/SiGe quantum dot. Our estimates suggest that the anharmonic dot confining potential combined with a gradient in the transverse magnetic field dominates the second harmonic response. As expected, the Rabi frequency depends quadratically on the driving amplitude and the periodicity with respect to the phase of the drive is twice that of the fundamental harmonic. The maximum Rabi frequency observed for the second harmonic is just a factor of two lower than that achieved for the first harmonic when driving at the same power. Combined with the lower demands on microwave circuitry when operating at half the qubit frequency, these observations indicate that second harmonic driving can be a useful technique for future quantum computation architectures.


title: Analytic derivation of electrostrictive tensors and their application to optical force density calculations

abstract: Using multiple scattering theory, we derived for the first time analytical formulas for electrostrictive tensors for two dimensional metamaterial systems. The electrostrictive tensor terms are found to depend explicitly on the symmetry of the underlying lattice of the metamaterial and they also depend explicitly on the direction of a local effective wave vector. These analytical results enable us to calculate light induced body forces inside a composite system (metamaterial) using the Helmholtz stress tensor within the effective medium formalism in the sense that the fields used in the stress tensor are those obtained by solving the macroscopic Maxwell equation with the microstructure of the metamaterial replaced by an effective medium. Our results point to some fundamental questions of using an effective medium theory to determine optical force density. In particular, the fact that Helmholtz tensor carries electrostrictive terms that are explicitly symmetry dependent means that the standard effective medium parameters cannot give sufficient information to determine body force density, even though they can give the correct total force. A more challenging issue is that the electrostrictive terms are related to a local effective wave vector, and it is not always obtainable in systems with boundary reflections within the context of a standard effective medium approach.


title: Almost sure local wellposedness of energy critical fractional Schr\"odinger equations with hartree nonlinearity

abstract: We consider a Cauchy problem of energy-critical fractional Schr\"odinger equation with Hartree nonlinearity below the energy space. Using a method of randomization of functions on $\mathbb{R}^d$ associated with the Wiener decomposition, introduced by \'{A}. B\'{e}nyi, T. Oh, and O. Pocovnicu \cite{beohpo1,beohpo2}, we prove that the Cauchy problem is almost surely locally well-posed. Our result includes Hartree Schr\"odinger equation ($\alpha = 2$).


title: Stochastic differential equations with variable structure driven by multiplicative Gaussian noise and sliding mode dynamic

abstract: This work is concerned with existence of weak solutions to discon- tinuous stochastic differential equations driven by multiplicative Gaus- sian noise and sliding mode control dynamics generated by stochastic differential equations with variable structure, that is with jump nonlin- earity. The treatment covers the finite dimensional stochastic systems and the stochastic diffusion equation with multiplicative noise.


title: The Essentially Entangled Component of Multipartite Mixed Quantum States, its Properties and an Efficient Algorithm for its Extraction

abstract: We introduce with geometric means a density matrix decomposition of a multipartite quantum system of a finite dimension into two density matrices: a separable one, also known as the best separable approximation, and an essentially entangled one, which contains no product states components. We show that this convex decomposition can be achieved in practice with the help of an algorithm based on linear programming, which in the general case scales polynomially with the dimension of the multipartite system. Furthermore, we suggest methods for analyzing the multipartite entanglement content of the essentially entangled component and derive analytically an upper bound for its rank. We illustrate the algorithm at an example of a composed system of total dimension 12 undergoing loss of coherence due to classical noise and we trace the time evolution of its essentially entangled component. We suggest a "geometric" description of entanglement dynamics and show how it explains the well-known phenomena of sudden death and revival of multipartite entanglement.


title: Multilevel Monte Carlo simulation of a diffusion with non-smooth drift

abstract: We show that Lasso and Bayesian Lasso are very close when the sparsity is large and the noise is small. Then we propose to solve Bayesian Lasso using multivalued stochastic differential equation. We obtain three discretizations algorithms, and propose a method for calculating the cost of Monte-Carlo (MC), multilevel Monte Carlo (MLMC) and MCMC algorithms.


title: An Entropy Stable Central Solver for Euler Equations

abstract: An exact discontinuity capturing central solver developed recently, named MOVERS (Method of Optimal Viscosity for Enhanced Resolution of Shocks, J Computat Phys 2009;228:770-798), is analyzed and improved further to make it entropy stable. MOVERS, which is designed to capture steady shocks and contact discontinuities exactly by enforcing the Rankine-Hugoniot jump condition directly in the discretization process, is a low diffusive algorithm in a simple central discretization framework, free of complicated Riemann solvers and flux splittings. However, this algorithm needs an entropy fix to avoid nonsmoothness in the expansion regions. The entropy conservation equation is used as a guideline to introduce an optimal numerical diffusion in the smooth regions and a limiter based switchover is introduced for numerical diffusion based on jump conditions at the large gradients. The resulting new scheme is entropy stable, accurate and captures steady discontinuities exactly while avoiding an entropy fix.


title: Hybrid Numerical Solvers for Massively Parallel Eigenvalue Computation and Their Benchmark with Electronic Structure Calculations

abstract: Optimally hybrid numerical solvers were constructed for massively parallel generalized eigenvalue problem (GEP).The strong scaling benchmark was carried out on the K computer and other supercomputers for electronic structure calculation problems in the matrix sizes of M = 10^4-10^6 with upto 105 cores. The procedure of GEP is decomposed into the two subprocedures of the reducer to the standard eigenvalue problem (SEP) and the solver of SEP. A hybrid solver is constructed, when a routine is chosen for each subprocedure from the three parallel solver libraries of ScaLAPACK, ELPA and EigenExa. The hybrid solvers with the two newer libraries, ELPA and EigenExa, give better benchmark results than the conventional ScaLAPACK library. The detailed analysis on the results implies that the reducer can be a bottleneck in next-generation (exa-scale) supercomputers, which indicates the guidance for future research. The code was developed as a middleware and a mini-application and will appear online.

categories: physics.comp-ph cond-mat.mtrl-sci cs.NA math.NA


title: Search for Interstellar Adenine

abstract: It is long debated if pre-biotic molecules are indeed present in the interstellar medium. Despite substantial works pointing to their existence, pre-biotic molecules are yet to be discovered with a complete confidence. In this paper, our main aim is to study the chemical evolution of interstellar adenine under various circumstances. We prepare a large gas-grain chemical network by considering various pathways for the formation of adenine. Majumdar et al. (2012) proposed that in the absence of adenine detection, one could try to trace two precursors of adenine, namely, HCCN and NH_2CN. Recently Merz et al. (2014), proposed another route for the formation of adenine in interstellar condition. They proposed two more precursor molecules. But it was not verified by any accurate gas-grain chemical model. Neither was it known if the production rate would be high or low. Our paper fills this important gap. We include this new pathways to find that the contribution through this pathways for the formation of Adenine is the most dominant one in the context of interstellar medium. We propose that observers may look for the two precursors (C_3NH and HNCNH) in the interstellar media which are equally important for predicting abundances of adenine. We perform quantum chemical calculations to find out spectral properties of adenine and its two new precursor molecules in infrared, ultraviolet and sub-millimeter region. Our present study would be useful for predicting abundance of adenine.


title: Integer Quantum Hall Effect in Graphene

abstract: We study the quantum Hall effect in a monolayer graphene by using an approach based on thermodynamical properties. This can be done by considering a system of Dirac particles in an electromagnetic field and taking into account of the edges effect as a pseudo-potential varying continuously along the $x$ direction. At low temperature and in the weak electric field limit, we explicitly determine the thermodynamical potential. With this, we derive the particle numbers in terms of the quantized flux and therefore the Hall conductivity immediately follows.

categories: cond-mat.mes-hall hep-th math-ph math.MP quant-ph


title: Asymmetric truncated Toeplitz operators and Toeplitz operators with matrix symbol

abstract: Truncated Toeplitz operators and their asymmetric versions are studied in the context of the Hardy space $H^p$ of the half-plane for $1<p<\infty$. It is shown that they are equivalent after extension to $2 \times 2$ matricial Toeplitz operators, which allows one to deduce information about their invertibility properties. Shifted model spaces are presented in the context of invariant subspaces, allowing one to deduce new Beurling--Lax theorems.


title: Kinetics of self-assembly via facilitated diffusion: formation of the transcription complex

abstract: We present an analytically solvable model for self-assembly of a molecular complex on a filament. The process is driven by a seed molecule that undergoes facilitated diffusion, which is a search strategy that combines diffusion in three-dimensions and one-dimension. Our study is motivated by single molecule level observations revealing the dynamics of transcription factors that bind to the DNA at early stages of transcription. We calculate the probability that a complex made up of a given number of molecules is completely formed, as well as the distribution of completion times, upon the binding of a seed molecule at a target site on the filament (without explicitly modeling the three-dimensional diffusion that precedes binding). We compare two different mechanisms of assembly where molecules bind in sequential and random order. Our results indicate that while the probability of completion is greater for random binding, the completion time scales exponentially with the size of the complex; in contrast, it scales as a power-law or slower for sequential binding, asymptotically. Furthermore, we provide model predictions for the dissociation and residence times of the seed molecule, which are observables accessible in single molecule tracking experiments.


title: Tur\'an type inequalities for regular Coulomb wave functions

abstract: Tur\'an, Mitrinovi\'c-Adamovi\'c and Wilker type inequalities are deduced for regular Coulomb wave functions. The proofs are based on a Mittag-Leffler expansion for the regular Coulomb wave function, which may be of independent interest. Moreover, some complete monotonicity results concerning the Coulomb zeta functions and some interlacing properties of the zeros of Coulomb wave functions are given.


title: Correlations of mutual positions of charge density waves nodes in side-by-side placed InAs wires measured with scanning gate microscopy

abstract: We investigate the correlations of mutual positions of charge density waves nodes in side-by-side placed InAs nanowires in presence of a conductive atomic force microscope tip served as a mobile gate at helium temperatures. Scanning gate microscopy scans demonstrate mutual correlation of positions of charge density waves nodes of two wires. A general mutual shift of the nodes positions and "crystal lattice mismatch" defect were observed. These observations demonstrate the crucial role of Coulomb interaction in formation of charge density waves in InAs nanowires.


title: Isoptic curves of generalized conic sections in the hyperbolic plane

abstract: After having investigated the real conic sections and their isoptic curves in the hyperbolic plane $\bH^2$ we consider the problem of the isoptic curves of generalized conic sections in the extended hyperbolic plane. This topic is widely investigated in the Euclidean plane $\BE^2$ (see for example \cite{Lo}), but in the hyperbolic and elliptic planes there are few results (see \cite{CsSz1}, \cite{CsSz2} and \cite{CsSz3}). In this paper we recall the former results on isoptic curves in the hyperbolic plane geometry, and define the notion of the generalised hyperbolic angle between proper and non-proper straight lines, summarize the notions of generalized hyperbolic conic sections classified by K.~Fladt in \cite{KF1} and \cite{KF2} and gy E.~Moln\'ar in \cite{M81}. Furthermore, we determine and visualize the generalized isoptic curves to all hyperbolic conic sections. We use for the computations the classical models which are based on the projective interpretation of the hyperbolic geometry and in this manner the isoptic curves can be visualized on the Euclidean screen of computer.


title: A Framework for Managing Evolving Information Resources on the Data Web

abstract: The web of data has brought forth the need to preserve and sustain evolving information within linked datasets; however, a basic requirement of data preservation is the maintenance of the datasets' structural characteristics as well. As open data are often found using different and/or heterogeneous data models and schemata from one source to another, there is a need to reconcile these mismatches and provide common denominations of interpretation on a multitude of levels, in order to be able to preserve and manage the evolution of the generated resources. In this paper, we present a linked data approach for the preservation and archiving of open heterogeneous datasets that evolve through time, at both the structural and the semantic layer. We first propose a set of re-quirements for modelling evolving linked datasets. We then proceed on concep-tualizing a modelling framework for evolving entities and place these in a 2x2 model space that consists of the semantic and the temporal dimensions.


title: Correlation functions of the KdV hierarchy and applications to intersection numbers over $\overline{\mathcal M}_{g,n}$

abstract: We derive an explicit generating function of correlations functions of an arbitrary tau-function of the KdV hierarchy. In particular applications, our formulation gives closed formul\ae\ of a new type for the generating series of intersection numbers of $\psi$-classes as well as of mixed $\psi$- and $\kappa$-classes in full genera.


title: Monte Carlo simulation to investigate the formation of molecular hydrogen and its deuterated forms

abstract: $H_2$ is the most abundant interstellar species. Its deuterated forms ($HD$ and $D_2$) are also significantly abundant. Huge abundances of these molecules could be explained by considering the chemistry occurring on the interstellar dust. Because of its simplicity, Rate equation method is widely used to study the formation of grain-surface species. However, since recombination efficiency of formation of any surface species are heavily dependent on various physical and chemical parameters, Monte Carlo method would be best method suited to take care of randomness of the processes. We perform Monte Carlo simulation to study the formation of $H_2$, $HD$ and $D_2$ on interstellar ices. Adsorption energies of surface species are the key inputs for the formation of any species on interstellar dusts but binding energies of deuterated species are yet to known with certainty. A zero point energy correction exists between hydrogenated and deuterated species which should be considered while modeling the chemistry on the interstellar dusts. Following some earlier studies, we consider various sets of adsorption energies to study the formation of these species in diverse physical circumstances. As expected, noticeable difference in these two approaches (Rate equation method and Monte Carlo method) is observed for production of these simple molecules on interstellar ices. We introduce two factors, namely, $S_f$ and $\beta$ to explain these discrepancies: $S_f$ is a scaling factor, which could be used to correlate discrepancies between Rate equation and Monte Carlo methods. $\beta$ factor indicates the formation efficiency under various circumstances. Higher values of $\beta$ indicates a lower production efficiency. We found that $\beta$ increases with a decrease in rate of accretion from gas phase to grain phase.


title: On Pairwise Compatibility of Some Graph (Super)Classes

abstract: A graph G=(V,E) is a pairwise compatibility graph (PCG) if there exists an edge-weighted tree T and two non-negative real numbers `d' and `D' such that each leaf `u' of T is a node of V and the edge `(u,v) belongs to E' iff `d <= d_T(u, v) <= D' where d_T(u, v) is the sum of weights of the edges on the unique path from `u' to `v' in T. The main issue on these graphs consists in characterizing them. In this note we prove the inclusion in the PCG class of threshold tolerance graphs and the non-inclusion of a number of intersection graphs, such as disk and grid intersection graphs, circular arc and tolerance graphs. The non-inclusion of some superclasses (trapezoid, permutation and rectangle intersection graphs) follows.


title: Gaussian and non-Gaussian fluctuations for mesoscopic linear statistics in determinantal processes

abstract: We study mesoscopic linear statistics for a class of determinantal point processes which interpolates between Poisson and Gaussian Unitary Ensemble statistics. These processes are obtained by modifying the spectrum of the correlation kernel of the GUE eigenvalue process. An example of such a system comes from considering the distribution of non-colliding Brownian motions in a cylindrical geometry, or a grand canonical ensemble of free fermions in a quadratic well at positive temperature. When the scale of the modification of the spectrum of the GUE kernel, related to the size of the cylinder or the temperature, is different from the scale in the mesoscopic linear statistic, we get a central limit theorem of either Poisson or GUE type. On the other hand, in the critical regime where the scales are the same, we get a non-Gaussian process in the limit. Its distribution is characterized by explicit but complicated formulae for the cumulants of smooth linear statistics. These results rely on an asymptotic sine-kernel approximation of the GUE kernel which is valid at all mesoscopic scales, and a generalization of cumulant computations of Soshnikov for the sine process. Analogous determinantal processes on the circle are also considered with similar results.


title: Pointed Hopf Algebras with Triangular Decomposition -- A Characterization of Multiparameter Quantum Groups

abstract: In this paper, we present an approach to the definition of multiparameter quantum groups by studying Hopf algebras with triangular decomposition. Classifying all of these Hopf algebras which are of what we call weakly separable type over a group, we obtain a class of pointed Hopf algebras which can be viewed as natural generalizations of multiparameter deformations of universal enveloping algebras of Lie algebras. These Hopf algebras are instances of a new version of braided Drinfeld doubles, which we call asymmetric braided Drinfeld doubles. This is a generalization of an earlier result by Benkart and Witherspoon (2004) who showed that two-parameter quantum groups are Drinfeld doubles. It is possible to recover a Lie algebra from these doubles in the case where the group is free abelian and the parameters are generic. The Lie algebras arising are generated by Lie subalgebras isomorphic to sl2.


title: Realization independent single time-delay dynamical model interpolation and $\mathcal{H}_2$-optimal approximation

abstract: In this paper, the realization-free model approximation problem, as stated in \cite{mayo2007framework,beattie2012realization}, is revisited in the case where the interpolating model might be time-delay dependent. To this aim, the Loewner framework, initially settled for delay-free realization, is firstly generalized to the single delay case. Secondly, the (infinite) model approximation $\mathcal{H}_2$ optimality conditions are established through the use of the Lambert functions. Finally, a numerically effective iterative scheme, named \textbf{dTF-IRKA}, similar to the \textbf{TF-IRKA} \cite{beattie2012realization}, is proposed to reach a part of the aforementioned optimality conditions. The proposed method validity and interest are assessed on different numerical examples.


title: The two-star model: exact solution in the sparse regime and condensation transition

abstract: The $2$-star model is the simplest exponential random graph model that displays complex behavior, such as degeneracy and phase transition. Despite its importance, this model has been solved only in the regime of dense connectivity. In this work we solve the model in the finite connectivity regime, far more prevalent in real world networks. We show that the model undergoes a condensation transition from a liquid to a condensate phase along the critical line corresponding, in the ensemble parameters space, to the Erd\"os-R\'enyi graphs. In the fluid phase the model can produce graphs with a narrow degree statistics, ranging from regular to Erd\"os-R\'enyi graphs, while in the condensed phase, the "excess" degree heterogeneity condenses on a single site with degree $\sim\sqrt{N}$. This shows the unsuitability of the two-star model, in its standard definition, to produce arbitrary finitely connected graphs with degree heterogeneity higher than Erd\"os-R\'enyi graphs and suggests that non-pathological variants of this model may be attained by softening the global constraint on the two-stars, while keeping the number of links hardly constrained.


title: $k$-extendibility of high-dimensional bipartite quantum states

abstract: The idea of detecting the entanglement of a given bipartite state by searching for symmetric extensions of this state was first proposed by Doherty, Parrilo and Spedialeri. The complete family of separability tests it generates, often referred to as the hierarchy of $k$-extendibility tests, has already proved to be most promising. The goal of this paper is to try and quantify the efficiency of this separability criterion in typical scenarios. For that, we essentially take two approaches. First, we compute the average width of the set of $k$-extendible states, in order to see how it scales with the one of separable states. And second, we characterize when random-induced states are, depending on the ancilla dimension, with high probability violating or not the $k$-extendibility test, and compare the obtained result with the corresponding one for entanglement vs separability. The main results can be precisely phrased as follows: on $\mathbf{C}^d\otimes\mathbf{C}^d$, when $d$ grows, the average width of the set of $k$-extendible states is equivalent to $(2/\sqrt{k})/d$, while random states obtained as partial traces over an environment $\mathbf{C}^s$ of uniformly distributed pure states are violating the $k$-extendibility test with probability going to $1$ if $s<((k-1)^2/4k)d^2$. Both statements converge to the conclusion that, if $k$ is fixed, $k$-extendibility is asymptotically a weak approximation of separability, even though any of the other well-studied separability relaxations is outperformed by $k$-extendibility as soon as $k$ is above a certain (dimension independent) value.


title: Epistemic nature of quantum reasoning

abstract: Doubts are raised concerning the usual interpretation of the alleged failure, by quantum mechanics, of the distributive law of classical logic. The difficulty raised by incompatible sets of observables is overcome within an epistemic enrichment of classical logic that provides the means for distinguishing between the value of a variable and its observation while retaining the classical connectives.


title: 3-D Velocity Regulation for Nonholonomic Source Seeking Without Position Measurement

abstract: We consider a three-dimensional problem of steering a nonholonomic vehicle to seek an unknown source of a spatially distributed signal field without any position measurement. In the literature, there exists an extremum seeking-based strategy under a constant forward velocity and tunable pitch and yaw velocities. Obviously, the vehicle with a constant forward velocity may exhibit certain overshoots in the seeking process and can not slow down even it approaches the source. To resolve this undesired behavior, this paper proposes a regulation strategy for the forward velocity along with the pitch and yaw velocities. Under such a strategy, the vehicle slows down near the source and stays within a small area as if it comes to a full stop, and controllers for angular velocities become succinct. We prove the local exponential convergence via the averaging technique. Finally, the theoretical results are illustrated with simulations.


title: Symmetry reduction, integrability and reconstruction in k-symplectic field theory

abstract: We investigate the reduction process of a k-symplectic field theory whose Lagrangian is invariant under a symmetry group. We give explicit coordinate expressions of the resulting reduced partial differential equations, the so-called Lagrange-Poincare field equations. We discuss two issues about reconstructing a solution from a given solution of the reduced equations. The first one is an interpretation of the integrability conditions, in terms of the curvatures of some connections. The second includes the introduction of the concept of a k-connection to provide a reconstruction method. We show that an invariant Lagrangian, under suitable regularity conditions, defines a `mechanical' k-connection.


title: The dichotomy structure of Y chromosome Haplogroup N

abstract: Haplogroup N-M231 of human Y chromosome is a common clade from Eastern Asia to Northern Europe, being one of the most frequent haplogroups in Altaic and Uralic-speaking populations. Using newly discovered bi-allelic markers from high-throughput DNA sequencing, we largely improved the phylogeny of Haplogroup N, in which 16 subclades could be identified by 33 SNPs. More than 400 males belonging to Haplogroup N in 34 populations in China were successfully genotyped, and populations in Northern Asia and Eastern Europe were also compared together. We found that all the N samples were typed as inside either clade N1-F1206 (including former N1a-M128, N1b-P43 and N1c-M46 clades), most of which were found in Altaic, Uralic, Russian and Chinese-speaking populations, or N2-F2930, common in Tibeto-Burman and Chinese-speaking populations. Our detailed results suggest that Haplogroup N developed in the region of China since the final stage of late Paleolithic Era.


title: The role of the Wigner distribution function in iterative ptychography

abstract: Ptychography employs a set of diffraction patterns that capture redundant information about an illuminated specimen as a localized beam is moved over the specimen. The robustness of this method comes from the redundancy of the dataset that in turn depends on the amount of oversampling and the form of the illumination. Although the role of oversampling in ptychography is fairly well understood, the same cannot be said of the illumination structure. This paper provides a vector space model of ptychography that accounts for the illumination structure in a way that highlights the role of the Wigner distribution function in iterative ptychography.


title: Interpreting the X-ray state transitions of Cygnus X-1

abstract: We present a novel method for interpreting observations of high-mass X-ray binaries (HMXBs) based on a combination of spectroscopic data and numerical results from a radiation hydrodynamic model of stellar winds. We calculate synthetic Doppler tomograms of predicted emission in low/hard and high/soft X-ray states and compare them with Doppler tomograms produced using spectra of Cygnus X-1, a prototype of HMXBs. Emission from HMXBs is determined by local conditions within the circumstellar medium, namely density, temperature, and ionization state. These quantities depend strongly on the X-ray state of the systems. By increasing intensity of an X-ray emission produced by the compact companion in the HMXB model, we achieved a complete redistribution of the circumstellar medium in the vicinity of the modelled system. These changes (which simulate the transitions between two major spectral states) are also apparent in the synthetic Doppler tomograms which are in good agreement with the observations.


title: On the Dirichlet problem in billiard spaces

abstract: The constrained Dirichlet boundary value problem $\ddot x=f(t,x)$, $x(0)=x(T)$, is studied in billiard spaces, where impacts occur in boundary points. Therefore we develop the research on impulsive Dirichlet problems with state-dependent impulses. Inspiring simple examples lead to an approach enabling to obtain both the existence and multiplicity results in one dimensional billiards. Several observations concerning the multidimensional case are also given.


title: A Luna \'etale slice theorem for algebraic stacks

abstract: We prove that every algebraic stack, locally of finite type over an algebraically closed field with affine stabilizers, is \'etale-locally a quotient stack in a neighborhood of a point with a linearly reductive stabilizer group. The proof uses an equivariant version of Artin's algebraization theorem proved in the appendix. We provide numerous applications of the main theorems.


title: Comparison between RHD simulation of supercritical accretion flows and steady model with outflows

abstract: We apply our two-dimensional (2D), radially self-similar steady-state accretion flow model to the analysis of hydrodynamic simulation results of supercritical accretion flows. Self-similarity is checked and the input parameters for the model calculation, such as advective factor and heat capacity ratio, are obtained from time-averaged simulation data. Solutions of the model are then calculated and compared with the simulation results. We find that in the converged region of the simulation, excluding the part too close to the black hole, the radial distribution of azimuthal velocity $v_\phi$, density $\rho$ and pressure $p$ basically follows the self-similar assumptions, i.e. they are roughly proportional to $r^{-0.5}$, $r^{-n}$, and $r^{-(n+1)}$, respectively, where $n\sim0.85$ for the mass injection rate of $1000L_\mathrm{E}/c^2$, and $n\sim0.74$ for $3000L_\mathrm{E}/c^2$. The distribution of $v_r$ and $v_\theta$ agrees less with self-similarity, possibly due to convective motions in the $r\theta$ plane. The distribution of velocity, density and pressure in $\theta$ direction obtained by the steady model agrees well with the simulation results {within the calculation boundary of the steady model}. Outward mass flux in the simulations is overall directed toward polar angle of 0.8382 rad ($\sim 48.0^\circ$) for $1000L_\mathrm{E}/c^2$, and 0.7852 rad ($\sim 43.4^\circ$) for $3000L_\mathrm{E}/c^2$, and $\sim$94\% of the mass inflow are driven away as outflow, while outward momentum and energy fluxes are focused around the polar axis. Part of these fluxes lie in the region that are not calculated by the steady model, and special attention should be paid when the model is applied.


title: A compression algorithm for the combination of PDF sets

abstract: The current PDF4LHC recommendation to estimate uncertainties due to parton distribution functions (PDFs) in theoretical predictions for LHC processes involves the combination of separate predictions computed using PDF sets from different groups, each of which comprises a relatively large number of either Hessian eigenvectors or Monte Carlo (MC) replicas. While many fixed-order and parton shower programs allow the evaluation of PDF uncertainties for a single PDF set at no additional CPU cost, this feature is not universal, and moreover the a posteriori combination of the predictions using at least three different PDF sets is still required. In this work, we present a strategy for the statistical combination of individual PDF sets, based on the MC representation of Hessian sets, followed by a compression algorithm for the reduction of the number of MC replicas. We illustrate our strategy with the combination and compression of the recent NNPDF3.0, CT14 and MMHT14 NNLO PDF sets. The resulting Compressed Monte Carlo PDF (CMC-PDF) sets are validated at the level of parton luminosities and LHC inclusive cross-sections and differential distributions. We determine that around 100 replicas provide an adequate representation of the probability distribution for the original combined PDF set, suitable for general applications to LHC phenomenology.


title: Real-time imaging of density ducts between the plasmasphere and ionosphere

abstract: Ionization of the Earth's atmosphere by sunlight forms a complex, multi-layered plasma environment within the Earth's magnetosphere, the innermost layers being the ionosphere and plasmasphere. The plasmasphere is believed to be embedded with cylindrical density structures (ducts) aligned along the Earth's magnetic field, but direct evidence for these remains scarce. Here we report the first direct wide-angle observation of an extensive array of field-aligned ducts bridging the upper ionosphere and inner plasmasphere, using a novel ground-based imaging technique. We establish their heights and motions by feature-tracking and parallax analysis. The structures are strikingly organized, appearing as regularly-spaced, alternating tubes of overdensities and underdensities strongly aligned with the Earth's magnetic field. These findings represent the first direct visual evidence for the existence of such structures.


title: On the BLM optimal renormalization scale setting for semihard processes

abstract: The BFKL approach for the investigation of semihard processes is plagued by large next-to-leading corrections, both in the kernel of the universal BFKL Green's function and in the process-dependent impact factors, as well as by large uncertainties in the renormalization scale setting. All that calls for some optimization procedure of the perturbative series. In this respect, one of the most common methods is the Brodsky-Lepage-Mackenzie (BLM) one, that eliminates the renormalization scale ambiguity by absorbing the non-conformal $\beta_0$-terms into the running coupling. In this paper, we apply BLM scale setting procedure directly to the amplitudes (cross sections) of several semihard processes. We show that, due to the presence of $\beta_0$-terms in the next-to-leading expressions for the impact factors, the optimal renormalization scale is not universal, but depends both on the energy and on the type of process in question.


title: Signatures of a two million year old supernova in the spectra of cosmic ray protons, antiprotons and positrons

abstract: The locally observed cosmic ray spectrum has several puzzling features, such as the excess of positrons and antiprotons above $\sim 20$ GeV and the discrepancy in the slopes of the spectra of cosmic ray protons and heavier nuclei in the TeV-PeV energy range. We show that these features are consistently explained by a nearby source which was active $\sim 2$ Myr ago and has injected $(1-2)\times 10^{50}$ erg in cosmic rays. The transient nature of the source and its overall energy budget point to the supernova origin of this local cosmic ray source. The age of the supernova suggests that the local cosmic ray injection was produced by the same supernova that has deposited $^{60}$Fe isotopes in the deep ocean crust.


title: How are Forbush decreases related to interplanetary magnetic field enhancements ?

abstract: Aims. Forbush decrease (FD) is a transient decrease followed by a gradual recovery in the observed galactic cosmic ray intensity. We seek to understand the relationship between the FDs and near-Earth interplanetary magnetic field (IMF) enhancements associated with solar coronal mass ejections (CMEs). Methods. We use muon data at cutoff rigidities ranging from 14 to 24 GV from the GRAPES-3 tracking muon telescope to identify FD events. We select those FD events that have a reasonably clean profile, and magnitude > 0.25%. We use IMF data from ACE/WIND spacecrafts. We look for correlations between the FD profile and that of the one hour averaged IMF. We ask if the diffusion of high energy protons into the large scale magnetic field is the cause of the lag observed between the FD and the IMF. Results. The enhancement of the IMF associated with FDs occurs mainly in the shock-sheath region, and the turbulence level in the magnetic field is also enhanced in this region. The observed FD profiles look remarkably similar to the IMF enhancement profiles. The FDs typically lag the IMF enhancement by a few hours. The lag corresponds to the time taken by high energy protons to diffuse into the magnetic field enhancement via cross-field diffusion. Conclusions. Our findings show that high rigidity FDs associated with CMEs are caused primarily by the cumulative diffusion of protons across the magnetic field enhancement in the turbulent sheath region between the shock and the CME.


title: Speculative Segmented Sum for Sparse Matrix-Vector Multiplication on Heterogeneous Processors

abstract: Sparse matrix-vector multiplication (SpMV) is a central building block for scientific software and graph applications. Recently, heterogeneous processors composed of different types of cores attracted much attention because of their flexible core configuration and high energy efficiency. In this paper, we propose a compressed sparse row (CSR) format based SpMV algorithm utilizing both types of cores in a CPU-GPU heterogeneous processor. We first speculatively execute segmented sum operations on the GPU part of a heterogeneous processor and generate a possibly incorrect results. Then the CPU part of the same chip is triggered to re-arrange the predicted partial sums for a correct resulting vector. On three heterogeneous processors from Intel, AMD and nVidia, using 20 sparse matrices as a benchmark suite, the experimental results show that our method obtains significant performance improvement over the best existing CSR-based SpMV algorithms. The source code of this work is downloadable at https://github.com/bhSPARSE/Benchmark_SpMV_using_CSR


title: A Practical and Worst-Case Efficient Algorithm for Divisor Methods of Apportionment

abstract: Proportional apportionment is the problem of assigning seats to parties according to their relative share of votes. Divisor methods are the de-facto standard solution, used in many countries. In recent literature, there are two algorithms that implement divisor methods: one by Cheng and Eppstein (ISAAC, 2014) has worst-case optimal running time but is complex, while the other (Pukelsheim, 2014) is relatively simple and fast in practice but does not offer worst-case guarantees. We demonstrate that the former algorithm is much slower than the other in practice and propose a novel algorithm that avoids the shortcomings of both. We investigate the running-time behavior of the three contenders in order to determine which is most useful in practice.


title: On Petviashvili type methods for traveling wave computations: Acceleration techniques

abstract: In this paper a family of fixed point algorithms, generalizing the \PM method, is considered. A previous work studied the convergence of the methods. Presented here is a second part of the analysis, concerning the introduction of some acceleration techniques into the iterative procedures. The purpose of the research is two-fold: one is improving the performance of the methods in case of convergence and the second one is widening their application when generating traveling waves in nonlinear dispersive wave equations, transforming some divergent into convergent cases. Two families of acceleration techniques are considered: the vector extrapolation methods and the Anderson acceleration methods. A comparative study through several numerical experiments is carried out.


title: Linear interpolation method in ensemble Kohn-Sham and range-separated density-functional approximations for excited states

abstract: Gross-Oliveira-Kohn density functional theory (GOK-DFT) for ensembles is in principle very attractive, but has been hard to use in practice. A novel, practical model based on GOK-DFT for the calculation of electronic excitation energies is discussed. The new model relies on two modifications of GOK-DFT: use of range separation and use of the slope of the linearly-interpolated ensemble energy, rather than orbital energies. The range-separated approach is appealing as it enables the rigorous formulation of a multi-determinant state-averaged DFT method. In the exact theory, the short-range density functional, that complements the long-range wavefunction-based ensemble energy contribution, should vary with the ensemble weights even when the density is held fixed. This weight dependence ensures that the range-separated ensemble energy varies linearly with the ensemble weights. When the (weight-independent) ground-state short-range exchange-correlation functional is used in this context, curvature appears thus leading to an approximate weight-dependent excitation energy. In order to obtain unambiguous approximate excitation energies, we propose to interpolate linearly the ensemble energy between equiensembles. It is shown that such a linear interpolation method (LIM) can be rationalized and that it effectively introduces weight dependence effects. As proof of principle, LIM has been applied to He, Be, H$_2$ in both equilibrium and stretched geometries as well as the stretched HeH$^+$ molecule. Very promising results have been obtained for both single (including charge transfer) and double excitations with spin-independent short-range local and semi-local functionals. Even at the Kohn--Sham ensemble DFT level, that is recovered when the range-separation parameter is set to zero, LIM performs better than standard time-dependent DFT.


title: A test of hypotheses for random graph distributions built from EEG data

abstract: The theory of random graphs is being applied in recent years to model neural interactions in the brain. While the probabilistic properties of random graphs has been extensively studied in the literature, the development of statistical inference methods for this class of objects has received less attention. In this work we propose a non-parametric test of hypotheses to test if two samples of random graphs were originated from the same probability distribution. We show how to compute efficiently the test statistic and we study its performance on simulated data. We apply the test to compare graphs of brain functional network interactions built from electroencephalographic (EEG) data collected during the visualization of point light displays depicting human locomotion.


title: Effect of in-medium parameters of rho meson in its photoproduction reactions on nuclei

abstract: There exist model calculations showing the modification of the hadronic parameters of $\rho$ meson in the nuclear environment. From these parameters, we extract the $\rho$ meson nucleus optical potential and show the medium effect due to this potential on the $\rho$ meson mass distribution spectra in the photonuclear reaction. The calculated results reproduced reasonably the measured $e^+e^-$ invariant mass, i.e., $\rho$ meson mass, distribution spectra in the $( \gamma, \rho^0 \to e^+e^- )$ reaction on nuclei.


title: Wide enough Latin rectangles are perfects

abstract: Given two integers $m$ and $n$ with $m\leq n$, a Latin rectangle of size $m\times n$ is a bi-dimensional array with $m$ rows and $n$ columns filled with symbols from an alphabet with $n$ symbols, such that each row contains a permutation of the alphabet and each column contains no repeated symbols. Two rows $a$ and $b$ of a Latin rectangle $R$ define a permutation $R_{a,b}$ assigning the symbol $y$ to the symbol $x$ if they are in the same column, $x$ is in row $a$ and $y$ is in row $b$. A Latin rectangle $R$ is perfect is the permutation $R_{a,b}$ is cyclic, for each pair of rows $a$ and $b$. We prove that for each integer $m$ and each large enough odd integer $n$ there is a perfect Latin rectangle $R$ of size $m\times n$. It is a partial (asymptotic) answer to a well-known conjecture which says that the same property holds for each odd integer $m\leq n$.


title: Master equation approach to transient quantum transport incorporating with initial correlations

abstract: In this paper, the exact transient quantum transport of non-interacting nanostructures is investigated in the presence of initial system-lead correlations and initial lead-lead correlations for a device system coupled to general electronic leads. The exact master equation incorporating with initial correlations is derived through the extended quantum Langevin equation. The effects of the initial correlations are manifested through the time-dependent fluctuations contained explicitly in the exact master equation. The transient transport current incorporating with initial correlations is obtained from the exact master equation. The resulting transient transport current can be expressed in terms of the single-particle propagating and correlation Green functions of the device system. We show that the initial correlations can affect quantum transport not only in the transient regime, but also in the steady-state limit when system-lead couplings are strong enough so that electron localized bound states occur in the device system.