title: Parity-Violating Hydrodynamics in 2+1 Dimensions

abstract: We study relativistic hydrodynamics of normal fluids in two spatial dimensions. When the microscopic theory breaks parity, extra transport coefficients appear in the hydrodynamic regime, including the Hall viscosity, and the anomalous Hall conductivity. In this work we classify all the transport coefficients in first order hydrodynamics. We then use properties of response functions and the positivity of entropy production to restrict the possible coefficients in the constitutive relations. All the parity-breaking transport coefficients are dissipationless, and some of them are related to the thermodynamic response to an external magnetic field and to vorticity. In addition, we give a holographic example of a strongly interacting relativistic fluid where the parity-violating transport coefficients are computable.


title: A Criticism on "A Mathematician's Apology" by G. H. Hardy

abstract: This paper is a criticism on "A Mathematician's Apology" by G. H. Hardy.


title: Unstable electromagnetic modes in strongly magnetized plasmas

abstract: The electromagnetic modes possibly unstable in strongly magnetized plasmas are identified. The regime where this instability might stand out compared to the incoherent electron-cyclotron radiation is explored. These modes are relevant to the inertial confinement fusion and the gamma ray burst.


title: NLO corrections for the dipole factorization of DIS structure functions at low x

abstract: The NLO generalization of the dipole factorization formula for the structure functions F_2 and F_L at low x is calculated using light front perturbation theory. That result gives some interesting insight into the kinematics of initial state parton showers in mixed space.


title: Bilocal versus non-bilocal correlations in entanglement swapping experiments

abstract: Entanglement swapping is a process by which two initially independent quantum systems can become entangled and generate nonlocal correlations. To characterize such correlations, we compare them to those predicted by bilocal models, where systems that are initially independent are described by uncorrelated states. We extend in this paper the analysis of bilocal correlations initiated in [Phys. Rev. Lett. 104, 170401 (2010)]. In particular, we derive new Bell-type inequalities based on the bilocality assumption in different scenarios, we study their possible quantum violations, and analyze their resistance to experimental imperfections. The bilocality assumption, being stronger than Bell's standard local causality assumption, lowers the requirements for the demonstration of quantumness in entanglement swapping experiments.


title: Exploiting boundary states of imperfect spin chains for high-fidelity state transfer

abstract: We study transfer of a quantum state through XX spin chains with static imperfections. We combine the two standard approaches for state transfer based on (i) modulated couplings between neighboring spins throughout the spin chain and (ii) weak coupling of the outermost spins to an unmodulated spin chain. The combined approach allows us to design spin chains with modulated couplings and localized boundary states, permitting high-fidelity state transfer in the presence of random static imperfections of the couplings. The modulated couplings are explicitly obtained from an exact algorithm using the close relation between tridiagonal matrices and orthogonal polynomials [Linear Algebr. Appl. 21, 245 (1978)]. The implemented algorithm and a graphical user interface for constructing spin chains with boundary states (spinGUIn) are provided as Supplemental Material.


title: Stability analysis of collisionless plasmas with specularly reflecting boundary

abstract: In this paper we provide sharp criteria for linear stability or instability of equilibria of collisionless plasmas in the presence of boundaries. Specifically, we consider the relativistic Vlasov-Maxwell system with specular reflection at the boundary for the particles and with the perfectly conducting boundary condition for the electromagnetic field. Here we initiate our investigation in the simple geometry of radial and longitudinal symmetry.


title: Practical experimental certification of computational quantum gates via twirling

abstract: Due to the technical difficulty of building large quantum computers, it is important to be able to estimate how faithful a given implementation is to an ideal quantum computer. The common approach of completely characterizing the computation process via quantum process tomography requires an exponential amount of resources, and thus is not practical even for relatively small devices. We solve this problem by demonstrating that twirling experiments previously used to characterize the average fidelity of quantum memories efficiently can be easily adapted to estimate the average fidelity of the experimental implementation of important quantum computation processes, such as unitaries in the Clifford group, in a practical and efficient manner with applicability in current quantum devices. Using this procedure, we demonstrate state-of-the-art coherent control of an ensemble of magnetic moments of nuclear spins in a single crystal solid by implementing the encoding operation for a 3 qubit code with only a 1% degradation in average fidelity discounting preparation and measurement errors. We also highlight one of the advances that was instrumental in achieving such high fidelity control.


title: The Kepler light curves of V1504 Cygni and V344 Lyrae: A study of the Outburst Properties

abstract: We examine the Kepler light curves of V1504 Cyg and V344 Lyr, encompassing ~736 d at 1 min cadence. During this span each system exhibited ~64-65 outbursts, including six superoutbursts. We find that, in both systems, the normal outbursts between two superoutbursts increase in duration over time by a factor ~1.2-1.9, and then reset to a small value after the following superoutburst. In both systems the trend of quiescent intervals between normal outbursts is to increase to a local maximum about half way through the supercycle - the interval from one superoutburst to the next - and then to decrease back to a small value by the time of the next superoutburst. This is inconsistent with Osaki's thermal-tidal model, which predicts a monotonic increase in the quiescent intervals between normal outbursts during a supercycle. Also, most of the normal outbursts have an asymmetric, fast-rise/slower-decline shape, consistent with outbursts triggered at large radii. The exponential rate of decay of the plateau phase of the superoutbursts is 8 d/mag for V1504 Cyg and 12 d/mag for V344 Lyr. This time scale gives a direct measure of the viscous time scale in the outer accretion disk given the expectation that the entire disk is in the hot, viscous state during superoutburst. The resulting constraint on the Shakura-Sunyaev parameter, alpha_{hot} ~ 0.1, is consistent with the value inferred from the fast dwarf nova decays. By looking at the slow decay rate for superoutbursts, which occur in systems below the period gap, in combination with the slow decay rate in one long outburst above the period gap (in U Gem), we infer a steep dependence of the decay rate on orbital period for long outbursts. This implies a steep dependence of alpha_{cold} on orbital period, consistent with tidal torquing as being the dominant angular momentum transport mechanism in quiescent disks in interacting binary systems.


title: Slime mould imitation of Belgian transport networks: redundancy, bio-essential motorways, and dissolution

abstract: Belgium is amongst few artificial countries, established on purpose, when Dutch and French speaking parts were joined in a single unit. This makes Belgium a particularly interesting testbed for studying bio-inspired techniques for simulation and analysis of vehicular transport networks. We imitate growth and formation of a transport network between major urban areas in Belgium using the acellular slime mould Physarum polycephalum. We represent the urban areas with the sources of nutrients. The slime mould spans the sources of nutrients with a network of protoplasmic tubes. The protoplasmic tubes represent the motorways. In an experimental laboratory analysis we compare the motorway network approximated by P. polycephalum and the man-made motorway network of Belgium. We evaluate the efficiency of the slime mould network and the motorway network using proximity graphs.


title: The Friedmann-Lemaitre-Robertson-Walker Big Bang singularities are well behaved

abstract: We show that the Big Bang singularity of the Friedmann-Lemaitre-Robertson-Walker model does not raise major problems to General Relativity. We prove a theorem showing that the Einstein equation can be written in a non-singular form, which allows the extension of the spacetime before the Big Bang. The physical interpretation of the fields used is discussed. These results follow from our research on singular semi-Riemannian geometry and singular General Relativity.

categories: gr-qc astro-ph.CO math-ph math.DG math.MP


title: Pole prescription of higher order induced vertices in Lipatov's QCD effective action

abstract: We investigate Lipatov's QCD effective action for the QCD high energy limit and propose a pole prescription for higher order induced vertices. The latter can be used in the evaluation of loop corrections to high energy factorized matrix elements within the effective action approach. The proposed prescription respects the symmetry properties of the unregulated vertices. Explicit results are presented up to third order in the gauge coupling, while an iterative procedure for higher orders is proposed.


title: Detecting Triaxiality in the Galactic Dark Matter halo through Stellar Kinematics

abstract: Assuming the dark matter halo of the Milky Way as a non-spherical potential (i.e. triaxial, prolate, oblate), we show how the assembling process of the Milky Way halo, may have left long lasting stellar halo kinematic fossils only due to the shape of the dark matter halo. In contrast with tidal streams, associated with recent satellite accretion events, these stellar kinematic groups will typically show inhomogeneous chemical and stellar population properties. However, they may be dominated by a single accretion event for certain mass assembling histories. If the detection of these peculiar kinematic stellar groups is confirmed, they would be the smoking gun for the predicted triaxiality of dark halos in cosmological galaxy formation scenarios.


title: Polynomial Bridgeland Stable Objects and Reflexive Sheaves

abstract: On a smooth projective threefold, we show that there are only two isomorphism types for the moduli of stable objects with respect to Bayer's standard polynomial Bridgeland stability - the moduli of Gieseker-stable sheaves and the moduli of PT-stable objects - under the following assumptions: no two of the stability vectors are collinear, and the degree and rank of the objects are relatively prime. We also describe a close relation between the intersection of the moduli spaces of PT-stable and dual-PT-stable objects, and the moduli of reflexive sheaves.


title: On Fano foliations

abstract: In this paper we address Fano foliations on complex projective varieties. These are foliations whose anti-canonical class is ample. We focus our attention on a special class of Fano foliations, namely del Pezzo foliations on complex projective manifolds. We show that these foliations are algebraically integrable, with one exceptional case when the ambient space is a projective space. We also provide a classification of del Pezzo foliations with mild singularities.


title: Coarse-Grained Modeling of a Deformable Nematic Vesicle

abstract: We develop a coarse-grained particle-based model to simulate membranes with nematic liquid-crystal order. The coarse-grained particles form vesicles which, at low temperature, have orientational order in the local tangent plane. As the strength of coupling between the nematic director and the vesicle curvature increases, the vesicles show a morphology transition from spherical to prolate and finally to a tube. We also observe the shape and defect arrangement around the tips of the prolate vesicle.


title: Kepler-20: A Sun-like Star with Three Sub-Neptune Exoplanets and Two Earth-size Candidates

abstract: We present the discovery of the Kepler-20 planetary system, which we initially identified through the detection of five distinct periodic transit signals in the Kepler light curve of the host star 2MASSJ19104752+4220194. We find a stellar effective temperature Teff=5455+-100K, a metallicity of [Fe/H]=0.01+-0.04, and a surface gravity of log(g)=4.4+-0.1. Combined with an estimate of the stellar density from the transit light curves we deduce a stellar mass of Mstar=0.912+-0.034 Msun and a stellar radius of Rstar=0.944^{+0.060}_{-0.095} Rsun. For three of the transit signals, our results strongly disfavor the possibility that these result from astrophysical false positives. We conclude that the planetary scenario is more likely than that of an astrophysical false positive by a factor of 2e5 (Kepler-20b), 1e5 (Kepler-20c), and 1.1e3 (Kepler-20d), sufficient to validate these objects as planetary companions. For Kepler-20c and Kepler-20d, the blend scenario is independently disfavored by the achromaticity of the transit: From Spitzer data gathered at 4.5um, we infer a ratio of the planetary to stellar radii of 0.075+-0.015 (Kepler-20c) and 0.065+-0.011 (Kepler-20d), consistent with each of the depths measured in the Kepler optical bandpass. We determine the orbital periods and physical radii of the three confirmed planets to be 3.70d and 1.91^{+0.12}_{-0.21} Rearth for Kepler-20b, 10.85 d and 3.07^{+0.20}_{-0.31} Rearth for Kepelr-20c, and 77.61 d and 2.75^{+0.17}_{-0.30} Rearth for Kepler-20d. From multi-epoch radial velocities, we determine the masses of Kepler-20b and Kepler-20c to be 8.7\+-2.2 Mearth and 16.1+-3.5 Mearth, respectively, and we place an upper limit on the mass of Kepler-20d of 20.1 Mearth (2 sigma).


title: Unsolved problems in particle physics

abstract: I consider selected (most important according to my own choice) unsolved problems in particle theory, both those related to extensions of the Standard Model (neutrino oscillations, which probably do not fit the usual three-generation scheme; indications in favour of new physics from astrophysical observations; electroweak symmetry breaking and hierarchy of parameters) and those which appear in the Standard Model (description of strong interactions at low and intermediate energies).


title: Creation and vanishing of superposition states

abstract: The evolution of a measured system and an experimental apparatus is presented in an unified form. Conditions under which the state of such a total system forms, evaluates and declines from a superposition of states are defined. The problem of definite outcomes in the quantum-mechanical measurement theory is considered.


title: Cosmological Radio Emission induced by WIMP Dark Matter

abstract: We present a detailed analysis of the radio synchrotron emission induced by WIMP dark matter annihilations and decays in extragalactic halos. We compute intensity, angular correlation, and source counts and discuss the impact on the expected signals of dark matter clustering, as well as of other astrophysical uncertainties as magnetic fields and spatial diffusion. Bounds on dark matter microscopic properties are then derived, and, depending on the specific set of assumptions, they are competitive with constraints from other indirect dark matter searches. At GHz frequencies, dark matter sources can become a significant fraction of the total number of sources with brightness below the microJansky level. We show that, at this level of fluxes (which are within the reach of the next-generation radio surveys), properties of the faint edge of differential source counts, as well as angular correlation data, can become an important probe for WIMPs.


title: VERITAS Observations of day-scale flaring of M87 in 2010 April

abstract: VERITAS has been monitoring the very-high-energy (VHE; >100GeV) gamma-ray activity of the radio galaxy M87 since 2007. During 2008, flaring activity on a timescale of a few days was observed with a peak flux of (0.70 +- 0.16) X 10^{-11} cm^{-2} s^{-1} at energies above 350GeV. In 2010 April, VERITAS detected a flare from M87 with peak flux of (2.71 +- 0.68) X 10^{-11} cm^{-2} s^{-1} for E>350GeV. The source was observed for six consecutive nights during the flare, resulting in a total of 21 hr of good quality data. The most rapid flux variation occurred on the trailing edge of the flare with an exponential flux decay time of 0.90^{+0.22}_{-0.15} days. The shortest detected exponential rise time is three times as long, at 2.87^{+1.65}_{-0.99} days. The quality of the data sample is such that spectral analysis can be performed for three periods: rising flux, peak flux, and falling flux. The spectra obtained are consistent with power-law forms. The spectral index at the peak of the flare is equal to 2.19 +- 0.07. There is some indication that the spectrum is softer in the falling phase of the flare than the peak phase, with a confidence level corresponding to 3.6 standard deviations. We discuss the implications of these results for the acceleration and cooling rates of VHE electrons in M87 and the constraints they provide on the physical size of the emitting region.


title: A comprehensive error rate for multiple testing

abstract: The higher criticism of a family of tests starts with the individual uncorrected p-values of each test. It then requires a procedure for deciding whether the collection of p-values indicates the presence of a real effect and if possible selects the ones that deserve closer scrutiny. This paper investigates procedures in which the ordered p-values are compared to an arbitrary positive and non-decreasing threshold sequence.


title: Vertical Graphene Base Transistor

abstract: We present a novel, graphene-based device concept for high-frequency operation: a hot electron graphene base transistor (GBT). Simulations show that GBTs have high current on/off ratios and high current gain. Simulations and small-signal models indicate that it potentially allows THz operation. Based on energy band considerations we propose a specific materials solution that is compatible with SiGe process lines.


title: Fermat-type equations of signature (13,13,p) via Hilbert cuspforms

abstract: In this paper we prove that equations of the form $x^{13} + y^{13} = Cz^{p}$ have no non-trivial primitive solutions (a,b,c) such that $13 \nmid c$ if $p > 4992539$ for an infinite family of values for $C$. Our method consists in relating a solution (a,b,c) to the previous equation to a solution (a,b,c_1) of another Diophantine equation with coefficients in $\Q(\sqrt{13})$. We then construct Frey-curves associated with (a,b,c_1) and we prove modularity of them in order to apply the modular approach via Hilbert cusp forms over $\Q(\sqrt{13})$. We also prove a modularity result for elliptic curves over totally real cyclic number fields of interest by itself.


title: A Very Common Fallacy in Quantum Mechanics: Superposition, Delayed Choice, Quantum Erasers, Retrocausality, and All That

abstract: There is a very common fallacy, here called the separation fallacy, that is involved in the interpretation of quantum experiments involving a certain type of separation such as the: double-slit experiments, which-way interferometer experiments, polarization analyzer experiments, Stern-Gerlach experiments, and quantum eraser experiments. It is the separation fallacy that leads not only to flawed textbook accounts of these experiments but to flawed inferences about retrocausality in the context of "delayed choice" versions of separation experiments.


title: Complexity and Algorithms for Euler Characteristic of Simplicial Complexes

abstract: We consider the problem of computing the Euler characteristic of an abstract simplicial complex given by its vertices and facets. We show that this problem is #P-complete and present two new practical algorithms for computing Euler characteristic. The two new algorithms are derived using combinatorial commutative algebra and we also give a second description of them that requires no algebra. We present experiments showing that the two new algorithms can be implemented to be faster than previous Euler characteristic implementations by a large margin.

categories: cs.CG cs.CC cs.DS cs.MS cs.SC math.AC math.CO


title: The three-loop form factor in N=4 super Yang-Mills

abstract: In this paper we study the Sudakov form factor in N=4 super Yang-Mills theory to the three-loop order. The latter is expressed in terms of planar and non-planar loop integrals. We show that it is possible to choose a representation in which each loop integral has uniform transcendentality. We verify analytically the expected exponentiation of the infrared divergences with the correct values of the three-loop cusp and collinear anomalous dimensions in dimensional regularisation. We find that the form factor in N=4 super Yang-Mills can be related to the leading transcendentality part of the quark and gluon form factors in QCD. We also study the ultraviolet properties of the form factor in D>4 dimensions, and find unexpected cancellations, resulting in an improved ultraviolet behaviour.


title: Unstable manifolds of Euler equations

abstract: We consider a steady state $v_{0}$ of the Euler equation in a fixed bounded domain in $\mathbf{R}^{n}$. Suppose the linearized Euler equation has an exponential dichotomy of unstable and center-stable subspaces. By rewriting the Euler equation as an ODE on an infinite dimensional manifold of volume preserving maps in $W^{k, q}$, $(k>1+\frac{n}{q})$, the unstable (and stable) manifolds of $v_{0}$ are constructed under certain spectral gap condition which is verified for both 2D and 3D examples. In particular, when the unstable subspace is finite dimensional, this implies the nonlinear instability of $v_{0}$ in the sense that arbitrarily small $W^{k, q}$ perturbations can lead to $L^{2}$ growth of the nonlinear solutions.


title: Mysteries around the graph Laplacian eigenvalue 4

abstract: We describe our current understanding on the phase transition phenomenon of the graph Laplacian eigenvectors constructed on a certain type of unweighted trees, which we previously observed through our numerical experiments. The eigenvalue distribution for such a tree is a smooth bell-shaped curve starting from the eigenvalue 0 up to 4. Then, at the eigenvalue 4, there is a sudden jump. Interestingly, the eigenvectors corresponding to the eigenvalues below 4 are semi-global oscillations (like Fourier modes) over the entire tree or one of the branches; on the other hand, those corresponding to the eigenvalues above 4 are much more localized and concentrated (like wavelets) around junctions/branching vertices. For a special class of trees called starlike trees, we obtain a complete understanding of such phase transition phenomenon. For a general graph, we prove the number of the eigenvalues larger than 4 is bounded from above by the number of vertices whose degrees is strictly higher than 2. Moreover, we also prove that if a graph contains a branching path, then the magnitudes of the components of any eigenvector corresponding to the eigenvalue greater than 4 decay exponentially from the branching vertex toward the leaf of that branch.


title: Coincidences between OVI and OVII Lines: Insights from High Resolution Simulations of the Warm-Hot Intergalactic Medium

abstract: With high resolution (0.46kpc/h), adaptive mesh-refinement Eulerian cosmological hydrodynamic simulations we compute properties of O VI and O VII absorbers from the warm-hot intergalactic medium (WHIM). Our new simulations are in broad agreement with previous simulations, with ~40% of the intergalactic medium being in the WHIM at z=0. It is found (1) The amount of gas in the WHIM at temperature below and above 10^6K is about equal within uncertainties. (1) Our simulations are in excellent agreement with observed properties of O VI absorbers, with respect to the line incidence rate and Doppler width-column density relation. (2) Velocity structures within absorbing regions are a significant, and for large Doppler width clouds, a dominant contributor to the Doppler widths of both O VI and O VII absorbers. A non-negligible fraction (in number and mass) of O VI and O VII clouds can arise from gas of temperature lower than 10^5, until the Doppler width is well in excess of 100km/s. (3) Strong O VI absorbers are predominantly collisionally ionized. About (61%, 57%, 39%) of O VI absorbers in the column density ranges of log N(OVI) cm^2=(12.5-13,13-14,>14) have temperature lower than 10^5K. (4) Quantitative prediction is made for the presence of broad and shallow O VI lines, which current observations may have largely missed. Upcoming observations by COS may be able to provide a test. (5) The reported 3 sigma upper limit on the mean column density of coincidental O VII lines at the location of detected O VI lines by Yao et al is above the predicted value by a factor of 2.5-4. (6) The claimed observational detection of O VII lines by Nicastro et al, if true, is 2 sigma above what our simulations predict.


title: Leptogenesis with heavy neutrino flavours: from density matrix to Boltzmann equations

abstract: Leptogenesis with heavy neutrino flavours is discussed within a density matrix formalism. We write the density matrix equation, describing the generation of the matter-antimatter asymmetry, for an arbitrary choice of the right-handed (RH) neutrino masses. For hierarchical RH neutrino masses lying in the fully flavoured regimes, this reduces to multiple-stage Boltzmann equations. In this case we recover and extend results previously derived within a quantum state collapse description. We confirm the generic existence of phantom terms. However, taking into account the effect of gauge interactions, we show that they are washed out at the production with a wash-out rate that is halved compared to that one acting on the total asymmetry. In the N_1-dominated scenario they cancel without contributing to the final baryon asymmetry. In other scenarios they do not in general and they have to be taken into account. We also confirm that there is a (orthogonal) component in the asymmetry produced by the heavier RH neutrinos which completely escapes the washout from the lighter RH neutrinos and show that phantom terms additionally contribute to it. The other (parallel) component is washed out with the usual exponential factor, even for weak washout. Finally, as an illustration, we study the two RH neutrino model in the light of the above findings, showing that phantom terms can contribute to the final asymmetry also in this case.


title: Evidence for 6{\Lambda}H

abstract: Evidence for the neutron-rich hypernucleus 6{\Lambda}H is presented from the FINUDA experiment at DA{\Phi}NE, Frascati, studying ({\pi}+, {\pi}-) pairs in coincidence from the K- +6Li \rightarrow 6 H+{\pi}+ production reaction followed by 6{\Lambda}H \rightarrow 6He + {\pi}- weak decay. The production rate of 6{\Lambda}H undergoing this two-body {\pi}- decay is determined to be (2.9\pm2.0)\cdot10-6/K-. Its binding energy, evaluated jointly from production and decay, is B{\Lambda}(6{\Lambda}H) = (4.0\pm1.1) MeV with respect to 5H+{\Lambda}. A systematic difference of (0.98 \pm 0.74) MeV between B{\Lambda} values derived separately from decay and from production is tentatively assigned to the 6{\Lambda}H 0+g.s. \rightarrow 1+ excitation.


title: Contrasting Probabilistic Scoring Rules

abstract: There are several scoring rules that one can choose from in order to score probabilistic forecasting models or estimate model parameters. Whilst it is generally agreed that proper scoring rules are preferable, there is no clear criterion for preferring one proper scoring rule above another. This manuscript compares and contrasts some commonly used proper scoring rules and provides guidance on scoring rule selection. In particular, it is shown that the logarithmic scoring rule prefers erring with more uncertainty, the spherical scoring rule prefers erring with lower uncertainty, whereas the other scoring rules are indifferent to either option.


title: BFKL Pomeron calculus: nucleus-nucleus scattering

abstract: In this paper the action of the BFKL Pomeron calculus is re-written in momentum representation, and the equations of motion for nucleus-nucleus collisions are derived, in this representation. We found the semi-classical solutions to these equations, outside of the saturation domain. Inside this domain these equations reduce to the set of delay differential equations, and their asymptotic solutions are derived.


title: Three-dimensional shapelets and an automated classification scheme for dark matter haloes

abstract: We extend the two-dimensional Cartesian shapelet formalism to d-dimensions. Concentrating on the three-dimensional case, we derive shapelet-based equations for the mass, centroid, root-mean-square radius, and components of the quadrupole moment and moment of inertia tensors. Using cosmological N-body simulations as an application domain, we show that three-dimensional shapelets can be used to replicate the complex sub-structure of dark matter halos and demonstrate the basis of an automated classification scheme for halo shapes. We investigate the shapelet decomposition process from an algorithmic viewpoint, and consider opportunities for accelerating the computation of shapelet-based representations using graphics processing units (GPUs).


title: Generating SU(Nc) pure gauge lattice QCD configurations on GPUs with CUDA

abstract: The starting point of any lattice QCD computation is the generation of a Markov chain of gauge field configurations. Due to the large number of lattice links and due to the matrix multiplications, generating SU(Nc) lattice QCD configurations is a highly demanding computational task, requiring advanced computer parallel architectures such as clusters of several Central Processing Units (CPUs) or Graphics Processing Units (GPUs). In this paper we present and explore the performance of CUDA codes for NVIDIA GPUs to generate SU(Nc) lattice QCD pure gauge configurations. Our implementation in one GPU uses CUDA and in multiple GPUs uses OpenMP and CUDA. We present optimized CUDA codes SU(2), SU(3) and SU(4). We also show a generic SU(Nc) code for Nc$\,\geq 4$ and compare it with the optimized version of SU(4). Our codes are publicly available for free use by the lattice QCD community.


title: An application of the method of moments to volatility estimation using daily high, low, opening and closing prices

abstract: We use the expectation of the range of an arithmetic Brownian motion and the method of moments on the daily high, low, opening and closing prices to estimate the volatility of the stock price. The daily price jump at the opening is considered to be the result of the unobserved evolution of an after-hours virtual trading day.The annualized volatility is used to calculate Black-Scholes prices for European options, and a trading strategy is devised to profit when these prices differ flagrantly from the market prices.


title: Continuants and some decompositions into squares

abstract: In 1855 H. J. S. Smith proved Fermat's two-square using the notion of palindromic continuants. In his paper, Smith constructed a proper representation of a prime number $p$ as a sum of two squares, given a solution of $z^2+1\equiv0\pmod{p}$, and vice versa. In this paper, we extend the use of continuants to proper representations by sums of two squares in rings of polynomials on fields of characteristic different from 2. New deterministic algorithms for finding the corresponding proper representations are presented. Our approach will provide a new constructive proof of the four-square theorem and new proofs for other representations of integers by quaternary quadratic forms.


title: Intractability of the Minimum-Flip Supertree problem and its variants

abstract: Computing supertrees is a central problem in phylogenetics. The supertree method that is by far the most widely used today was introduced in 1992 and is called Matrix Representation with Parsimony analysis (MRP). Matrix Representation using Flipping (MRF)}, which was introduced in 2002, is an interesting variant of MRP: MRF is arguably more relevant that MRP and various efficient implementations of MRF have been presented. From a theoretical point of view, implementing MRF or MRP is solving NP-hard optimization problems. The aim of this paper is to study the approximability and the fixed-parameter tractability of the optimization problem corresponding to MRF, namely Minimum-Flip Supertree. We prove strongly negative results.


title: Extended Josephson Relation and Abrikosov lattice deformation

abstract: From the point of view of time-dependent Ginzburg Landau (TDGL) theory, a Josephson-like relation is derived for an Abrikosov vortex lattice accelerated and deformed by applied fields. Beginning with a review of the Josephson Relation derived from the two ingredients of a lattice-kinematics assumption in TDGL theory and gauge invariance, we extend the construction to accommodate a time-dependent applied magnetic field, a floating-kernel formulation of normal current, and finally lattice deformation due to the electric field and inertial effects of vortex-lattice motion. The resulting Josephson-like relation, which we call an Extended Josephson Relation, applies to a much wider set of experimental conditions than the original Josephson Relation, and is explicitly compatible with the considerations of TDGL theory.


title: The Constraints of Conformal Symmetry on RG Flows

abstract: If the coupling constants in QFT are promoted to functions of space-time, the dependence of the path integral on these couplings is highly constrained by conformal symmetry. We begin the present note by showing that this idea leads to a new proof of Zamolodchikov's theorem. We then review how this simple observation also leads to a derivation of the a-theorem. We exemplify the general procedure in some interacting theories in four space-time dimensions. We concentrate on Banks-Zaks and weakly relevant flows, which can be controlled by ordinary and conformal perturbation theories, respectively. We compute explicitly the dependence of the path integral on the coupling constants and extract the change in the a-anomaly (this agrees with more conventional computations of the same quantity). We also discuss some general properties of the sum rule found in arXiv:1107.3987 and study it in several examples.


title: Implementation of a Parallel Tree Method on a GPU

abstract: The kd-tree is a fundamental tool in computer science. Among other applications, the application of kd-tree search (by the tree method) to the fast evaluation of particle interactions and neighbor search is highly important, since the computational complexity of these problems is reduced from O(N^2) for a brute force method to O(N log N) for the tree method, where N is the number of particles. In this paper, we present a parallel implementation of the tree method running on a graphics processing unit (GPU). We present a detailed description of how we have implemented the tree method on a Cypress GPU. An optimization that we found important is localized particle ordering to effectively utilize cache memory. We present a number of test results and performance measurements. Our results show that the execution of the tree traversal in a force calculation on a GPU is practical and efficient.


title: Investigation of the Formation and Separation of An EUV Wave from the Expansion of A Coronal Mass Ejection

abstract: We address the nature of EUV waves through direct observations of the formation of a diffuse wave driven by the expansion of a coronal mass ejection (CME) and its subsequent separation from the CME front. The wave and the CME on 2011 June 7 were well observed by Atmospheric Imaging Assembly onboard Solar Dynamic Observatory. Following the solar eruption onset, marked by the beginning of the rapid increasing of the CME velocity and the X-ray flux of accompanying flare, the CME exhibits a strong lateral expansion. During this impulsive expansion phase, the expansion speed of the CME bubble increases from 100 km s$^{-1}$ to 450 km s$^{-1}$ in only six minutes. An important finding is that a diffuse wave front starts to separate from the front of the expanding bubble shortly after the lateral expansion slows down. Also a type-II burst is formed near the time of the separation. After the separation, two distinct fronts propagate with different kinematic properties. The diffuse front travels across the entire solar disk; while the sharp front rises up, forming the CME ejecta with the diffuse front ahead of it. These observations suggest that the previously termed EUV wave is a composite phenomenon and driven by the CME expansion. While the CME expansion is accelerating, the wave front is cospatial with the CME front, thus the two fronts are indiscernible. Following the end of the acceleration phase, the wave moves away from the CME front with gradually an increasing distance between them.


title: Dimension and measure for typical random fractals

abstract: We consider the dimension and measure of typical attractors of random iterated function systems (RIFSs). We define a RIFS to be a finite set of (deterministic) iterated function systems (IFSs) acting on the same metric space and, for a given RIFS, we define a continuum of random attractors corresponding to each sequence of deterministic IFSs. Much work has been done on computing the 'almost sure' dimensions of these random attractors. Here we compute the typical dimensions (in the sense of Baire) and observe that our results are in stark contrast to those obtained using the probabilistic approach. Furthermore, we examine the typical Hausdorff and packing measures of the random attractors and give a number of examples to illustrate some of the strange phenomena that can occur. The only restriction we impose on the maps is that they are bi-Lipschitz and we obtain our dimension results without assuming any separation conditions.


title: Source monitoring for continuous-variable quantum key distribution

abstract: The noise in optical source needs to be characterized for the security of continuous-variable quantum key distribution (CVQKD). Two feasible schemes, based on either active optical switch or passive beamsplitter are proposed to monitor the variance of source noise, through which, Eve's knowledge can be properly estimated. We derive the security bounds for both schemes against collective attacks in the asymptotic case, and find that the passive scheme performs better.


title: Functional Integration and High Energy Scattering of Particles with Anomalous Magnetic Moments in Quantum Field Theory

abstract: The functional integration method is used for studying the scattering of a scalar pion on nucleon with the anomalous magnetic moment in the framework of nonrenomalizable quantum field theory. In the asymptotic region s {\to} {\infty}, |t| << s the representation of eikonal type for the amplitude of pion-nucleon scattering is obtained. The anomalous magnetic moment leads to additional terms in the amplitude which describe the spin flips in the scattering process. It is shown that the renormalization problem does not arise in the asymptotic s {\to} {\infty} since the unrenomalized divergences disappear in this approximation. Coulomb interference is considered as an application.


title: Orthogonality Catastrophe in Quantum Sticking

abstract: The probability that a particle will stick to a surface is fundamental to a variety of processes in surface science, including catalysis, epitaxial growth, and corrosion. At ultralow energies, how particles scatter or stick to a surface affects the performance of atomic clocks, matter-wave interferometers, atom chips and other quantum information processing devices. In this energy regime, the sticking probability is influenced by a distinctly quantum mechanical effect: quantum reflection, a result of matter wave coherence, suppresses the probability of finding the particle near the surface and reduces the sticking probability. We find that another quantum effect can occur, further shaping the sticking probability: the orthogonality catastrophe, a result of the change in the quantum ground state of the surface in the presence of a particle, can dramatically alter the probability for quantum sticking and create a superreflective surface at low energies.


title: Synchronization of Huygens' clocks and the Poincare method

abstract: We study two models of connected pendulum clocks synchronizing their oscillations, a phenomenon originally observed by Huygens. The oscillation angles are assumed to be small so that the pendulums are modeled by harmonic oscillators, clock escapements are modeled by the van der Pol terms. The mass ratio of the pendulum bobs to their casings is taken as a small parameter. Analytic conditions for existence and stability of synchronization regimes, and analytic expressions for their stable amplitudes and period corrections are derived using the Poincare theorem on existence of periodic solutions in autonomous quasi-linear systems. The anti-phase regime always exists and is stable under variation of the system parameters. The in-phase regime may exist and be stable, exist and be unstable, or not exist at all depending on parameter values. As the damping in the frame connecting the clocks is increased the in-phase stable amplitude and period are decreasing until the regime first destabilizes and then disappears. The results are most complete for the traditional three degrees of freedom model, where the clock casings and the frame are consolidated into a single mass.


title: Computational studies for reduced graphene oxide in hydrogen-rich environment

abstract: We employ molecular dynamic simulations to study the reduction process of graphene-oxide (GO) in a chemically active environment enriched with hydrogen. We examine the concentration and pressure of hydrogen gas as a function of temperature in which abstraction of oxygen is possible with minimum damage to C-sp$^2$ bonds hence preserving the integrity of the graphene sheet. Through these studies we find chemical pathways that demonstrate beneficiary mechanisms for the quality of graphene including formation of water as well as suppression of carbonyl pair holes in favor of hydroxyl and epoxy formation facilitated by hydrogen gas in the environment.

categories: cond-mat.mtrl-sci physics.atm-clus physics.chem-ph physics.comp-ph


title: Handling a large bound for a problem on the generalized Pillai equation $\pm r a^x \pm s b^y = c$

abstract: We consider $N$, the number of solutions $(x,y,u,v)$ to the equation $ (-1)^u r a^x + (-1)^v s b^y = c $ in nonnegative integers $x, y$ and integers $u, v \in \{0,1\}$, for given integers $a>1$, $b>1$, $c>0$, $r>0$ and $s>0$. Previous work showed that there are nine essentially distinct $(a,b,c,r,s)$ for which $N \ge 4$, except possibly for cases in which the solutions have $r$, $a$, $x$, $s$, $b$, and $y$ each bounded by $8 \cdot 10^{14}$ or $2 \cdot 10^{15}$. In this paper we show that there are no further cases with $N \ge 4$ within these bounds. We note that $N = 3$ for an infinite number of $(a,b,c,r,s)$, even if we eliminate from consideration cases which are directly derived from other cases in one of several completely designated ways. Our work differs from previous work in that we allow $x$ and $y$ to be zero and also allow choices of $(u,v)$ other than $(0,1)$.


title: The equation $|p^x \pm q^y| = c$ in nonnegative $x$, $y$

abstract: We improve earlier work on the title equation (where $p$ and $q$ are primes and $c$ is a positive integer) by allowing $x$ and $y$ to be zero as well as positive. Earlier work on the title equation showed that, with listed exceptions, there are at most two solutions in positive integers $x$ and $y$, using elementary methods. Here we show that, with listed exceptions, there are at most two solutions in nonnegative integers $x$ and $y$, but the proofs are dependent on nonelementary work of Mignotte, Bennett, Luca, and Szalay. In order to provide some of our results with purely elementary proofs, we give short elementary proofs of the results of Luca, made possible by an elementary lemma which also has an application to the familiar equation $x^2 + C = y^n$. We also give shorter simpler proofs of Szalay's results. A summary of results on the number of solutions to the generalized Pillai equation $(-1)^u r a^x + (-1)^v s b^y = c$ is also given.


title: When is the Cuntz-Krieger algebra of a higher-rank graph approximately finite-dimensional?

abstract: We investigate the question: when is a higher-rank graph C*-algebra approximately finite dimensional? We prove that the absence of an appropriate higher-rank analogue of a cycle is necessary. We show that it is not in general sufficient, but that it is sufficient for higher-rank graphs with finitely many vertices. We give a detailed description of the structure of the C*-algebra of a row-finite locally convex higher-rank graph with finitely many vertices. Our results are also sufficient to establish that if the C*-algebra of a higher-rank graph is AF, then its every ideal must be gauge-invariant. We prove that for a higher-rank graph C*-algebra to be AF it is necessary and sufficient for all the corners determined by vertex projections to be AF. We close with a number of examples which illustrate why our question is so much more difficult for higher-rank graphs than for ordinary graphs.


title: Two Earth-sized planets orbiting Kepler-20

abstract: Since the discovery of the first extrasolar giant planets around Sun-like stars, evolving observational capabilities have brought us closer to the detection of true Earth analogues. The size of an exoplanet can be determined when it periodically passes in front of (transits) its parent star, causing a decrease in starlight proportional to its radius. The smallest exoplanet hitherto discovered has a radius 1.42 times that of the Earth's radius (R Earth), and hence has 2.9 times its volume. Here we report the discovery of two planets, one Earth-sized (1.03R Earth) and the other smaller than the Earth (0.87R Earth), orbiting the star Kepler-20, which is already known to host three other, larger, transiting planets. The gravitational pull of the new planets on the parent star is too small to measure with current instrumentation. We apply a statistical method to show that the likelihood of the planetary interpretation of the transit signals is more than three orders of magnitude larger than that of the alternative hypothesis that the signals result from an eclipsing binary star. Theoretical considerations imply that these planets are rocky, with a composition of iron and silicate. The outer planet could have developed a thick water vapour atmosphere.


title: Compact source of narrow-band counterpropagating polarization-entangled photon pairs using a single dual-periodically poled crystal

abstract: We propose a scheme for the generation of counterpropagating polarization-entangled photon pairs from a dual-periodically poled crystal. Compared with the usual forward-wave type source, this source, in the backward-wave way, has a much narrower bandwidth. With a 2-cm-long bulk crystal, the bandwidths of the example sources are estimated to be 3.6 GHz, and the spectral brightnesses are more than 100 pairs/(s GHz mW). Two concurrent quasi-phase-matched spontaneous parametric down-conversion processes in a single crystal enable our source to be compact and stable. This scheme does not rely on any state projection and applies to both degenerate and non-degenerate cases, facilitating applications of the entangled photons.


title: Single Spin Asymmetry through QCD Instantons

abstract: We revisit the effects of QCD instantons in semi-inclusive deep inelastic scattering (SIDIS). We show that large single spin asymmetry (SSA) effects can be induced in longitudinally and transversely polarized proton targets. The results are in agreement with most of the reported data for pion and kaon production. The same effects are found to be important in polarized proton on proton scattering for both charged and chargeless pion productions. The results agree with the reported data in a wide range of \sqrt{s} = 19.4-200 GeV. We predict the SSA for \pi^\pm production in p_\uparrow p in the collider range of \sqrt{s} = 62.4-500 GeV. The backward \pi^{\pm} and \pi^0 productions for the SSA in p_\uparrow p collisions are predicted to coincide at large \sqrt{s}.


title: Cooperative Algorithms for MIMO Amplify-and-Forward Relay Networks

abstract: Interference alignment is a signaling technique that provides high multiplexing gain in the interference channel. It can be extended to multi-hop interference channels, where relays aid transmission between sources and destinations. In addition to coverage extension and capacity enhancement, relays increase the multiplexing gain in the interference channel. In this paper, three cooperative algorithms are proposed for a multiple-antenna amplify-and-forward (AF) relay interference channel. The algorithms design the transmitters and relays so that interference at the receivers can be aligned and canceled. The first algorithm minimizes the sum power of enhanced noise from the relays and interference at the receivers. The second and third algorithms rely on a connection between mean square error and mutual information to solve the end-to-end sum-rate maximization problem with either equality or inequality power constraints via matrix-weighted sum mean square error minimization. The resulting iterative algorithms converge to stationary points of the corresponding optimization problems. Simulations show that the proposed algorithms achieve higher end-to-end sum-rates and multiplexing gains that existing strategies for AF relays, decode-and-forward relays, and direct transmission. The first algorithm outperforms the other algorithms at high signal-to-noise ratio (SNR) but performs worse than them at low SNR. Thanks to power control, the third algorithm outperforms the second algorithm at the cost of overhead.


title: Binomial ARMA count series from renewal processes

abstract: This paper describes a new method for generating stationary integer-valued time series from renewal processes. We prove that if the lifetime distribution of renewal processes is nonlattice and the probability generating function is rational, then the generated time series satisfy causal and invertible ARMA type stochastic difference equations. The result provides an easy method for generating integer-valued time series with ARMA type autocovariance functions. Examples of generating binomial ARMA(p,p-1) series from lifetime distributions with constant hazard rates after lag p are given as an illustration. An estimation method is developed for the AR(p) cases.


title: Characteristic Polynomials and Fixed Spaces of Semisimple Elements

abstract: Answering a question of Frank Calegari, we extend some of our earlier results on dimension of fixed point spaces of elements in irreducible linear groups. We consider characteristic polynomials rather than just fixed spaces.


title: Effect of dynamical traps on chaotic transport in a meandering jet flow

abstract: We continue our study of chaotic mixing and transport of passive particles in a simple model of a meandering jet flow [Prants, et al, Chaos {\bf 16}, 033117 (2006)]. In the present paper we study and explain phenomenologically a connection between dynamical, topological, and statistical properties of chaotic mixing and transport in the model flow in terms of dynamical traps, singular zones in the phase space where particles may spend arbitrary long but finite time [Zaslavsky, Phys. D {\bf 168--169}, 292 (2002)]. The transport of passive particles is described in terms of lengths and durations of zonal flights which are events between two successive changes of sign of zonal velocity. Some peculiarities of the respective probability density functions for short flights are proven to be caused by the so-called rotational-islands traps connected with the boundaries of resonant islands (including the vortex cores) filled with the particles moving in the same frame and the saddle traps connected with periodic saddle trajectories. Whereas, the statistics of long flights can be explained by the influence of the so-called ballistic-islands traps filled with the particles moving from a frame to frame.


title: Gamma-Dirichlet Structure and Two Classes of Measure-valued Processes

abstract: The Gamma-Dirichlet structure corresponds to the decomposition of the gamma process into the independent product of a gamma random variable and a Dirichlet process. This structure allows us to study the properties of the Dirichlet process through the gamma process and vice versa. In this article, we begin with a brief review of existing results concerning the Gamma-Dirichlet structure. New results are obtained for the large deviations of the jump sizes of the gamma process and the quasi-invariance of the two-parameter Poisson-Dirichlet distribution. The laws of the gamma process and the Dirichlet process are the respective reversible measures of the measure-valued branching diffusion with immigration and the Fleming-Viot process with parent independent mutation. We view the relation between these two classes of measure-valued processes as the dynamical Gamma-Dirichlet structure. Other results of this article include the derivation of the transition function of the Fleming-Viot process with parent independent mutation from the transition function of the measure-valued branching diffusion with immigration, and the establishment of the reversibility of the latter. One of these is related to an open problem by Ethier and Griffiths and the other leads to an alternative proof of the reversibility of the Fleming-Viot process.


title: Some Ramsey-type results on intrinsic linking of n-complexes

abstract: Define the complete n-complex on N vertices to be the n-skeleton of an (N-1)-simplex. We show that embeddings of sufficiently large complete n-complexes in R^{2n+1} necessarily exhibit complicated linking behaviour, thereby extending known results on embeddings of large complete graphs in R^3 (the case n=1) to higher dimensions. In particular, we prove the existence of links of the following types: r-component links, with the linking pattern of a chain, necklace or keyring; 2-component links with linking number at least lambda in absolute value; and 2-component links with linking number a non-zero multiple of a given integer q. For fixed n the number of vertices required for each of our results grows at most polynomially with respect to the parameter r, lambda or q.


title: Lifting in Frattini Covers and A Characterization of Finite Solvable Groups

abstract: We prove a lifting theorem for odd Frattini covers of finite groups. Using this, we characterize solvable groups and more generally p-solvable groups in terms of containing a triple of elements of distinct prime power orders with product 1. This is also related to various questions about Fried's modular tower program and properties of Hurwitz spaces of covers of curves.


title: Experimental observation of the spin Hall effect of light on a nano-metal film via weak measurements

abstract: We theorize the spin Hall effect of light (SHEL) on a nano-metal film and demonstrate it experimentally via weak measurements. A general propagation model to describe the relationship between the spin-orbit coupling and the thickness of the metal film is established. It is revealed that the spin-orbit coupling in the SHEL can be effectively modulated by adjusting the thickness of the metal film, and the transverse displacement is sensitive to the thickness of metal film in certain range for horizontal polarization light. Importantly, a large negative transverse shift can be observed as a consequence of the combined contribution of the ratio and the phase difference of Fresnel coefficients.


title: Adequate Subgroups II

abstract: The notion of adequate subgroups was introduced by Jack Thorne. It is a weakening of the notion of big subgroup used by Wiles and Taylor in proving automorphy lifting theorems for certain Galois representations. Using this idea, Thorne was able to prove some new lifting theorems. In an appendix to Thorne's article, it was shown in that certain groups were adequate. One of the key aspects was the question of whether the span of the semsimple elements in the group is the full endomorphism ring of an absolutely irreducible module. We show that this is the case in prime characteristic p for p-solvable groups as long the dimension is not divisible by p. We also observe that the condition holds for certain infinite groups. Finally, we present the first examples showing that this condition need not hold and give a negative answer to a question of Richard Taylor.


title: Time-reversal symmetry breaking phase in the Hubbard model: a VCA study

abstract: We study the stability of the time-reversal symmetry breaking staggered flux phase of a single band Hubbard model, within the Variational Cluster Approach (VCA). For intermediate and small values of the interaction $U$, we find metastable solutions for the staggered flux phase, with a maximum current per bond at $U\approx 3.2$. However, allowing for antiferromagnetic and superconducting long-range order it turns out that in the region at and close to half filling the antiferromagnetic phase is the most favorable energetically. The effect of nearest-neighbour interaction is also considered. Our results show that a negative nearest-neighbour interaction and finite doping favors the stability of the staggered-flux phase. We also present preliminary results for the three-band Hubbard model obtained with a restricted set of variational parameters. For this case, no spontaneous time-reversal symmetry breaking phase is found in our calculations.


title: Effect of both Z and Z'-mediated flavor-changing neutral currents on the baryonic rare decay $lambda_b$ decays into $lambda l^+l^- $

abstract: We study the effect of both Z and Z'-mediated flavor-changing neutral currents (FCNCs) on the $ lambda_b arrow lambda l^+l^- (l = mu, tau) $ rare decay. We find the branching ratio is reasonably enhanced from its standard model value due to the effect of both Z and Z'-mediated FCNCs, and gives the possibility of new physics beyond the standard model. The contribution of Z'-boson depends upon the precise value of the mass of Z' boson.


title: Gauge theory webs and surfaces

abstract: We analyze the perturbative cusp and closed polygons of Wilson lines for massless gauge theories in coordinate space, and express them as exponentials of two-dimensional integrals. These integrals have geometric interpretations, which link renormalization scales with invariant distances.


title: Minimax bounds for estimation of normal mixtures

abstract: This paper deals with minimax rates of convergence for estimation of density functions on the real line. The densities are assumed to be location mixtures of normals, a global regularity requirement that creates subtle difficulties for the application of standard minimax lower bound methods. Using novel Fourier and Hermite polynomial techniques, we determine the minimax optimal rate - slightly larger than the parametric rate - under squared error loss. For Hellinger loss, we provide a minimax lower bound using ideas modified from the squared error loss case.


title: Existence of Smooth Solutions to Coupled Chemotaxis-Fluid Equations

abstract: We consider a system coupling the parabolic-parabolic Keller-Segel equations to the in- compressible Navier-Stokes equations in spatial dimensions two and three. We establish the local existence of regular solutions and present some blow-up criteria. For two dimensional Navier-Stokes-Keller-Segel equations, regular solutions constructed locally in time are, in reality, extended globally under some assumptions pertinent to experimental observation in [20] on the consumption rate and chemotactic sensitivity. We also show the existence of global weak solutions in spatially three dimensions with rather restrictive consumption rate and chemotactic sensitivity.


title: Metallicity and Temperature Indicators in M dwarf K band Spectra: Testing New & Updated Calibrations With Observations of 133 Solar Neighborhood M dwarfs

abstract: We present K band spectra for 133 nearby (d < 33 parsecs) M dwarfs, including 18 M dwarfs with reliable metallicity estimates (as inferred from an FGK type companion), 11 M dwarf planet hosts, more than 2/3 of the M dwarfs in the Northern 8 pc sample, and several M dwarfs from the LSPM catalog. From these spectra, we measure equivalent widths of the Ca and Na lines, and a spectral index quantifying the absorption due to H2O opacity (the H2O-K2 index). Using empirical spectral types standards and synthetic models, we calibrate the H2O-K2 index as an indicator of an M dwarf's spectral type and effective temperature. We also present a revised relationship that estimates the [Fe/H] and [M/H] metallicities of M dwarfs from their Na I, Ca I, and H2O-K2 measurements. Comparisons to model atmosphere provide a qualitative validation of our approach, but also reveal an overall offset between the atomic line strengths predicted by models as compared to actual observations. Our metallicity estimates also reproduce expected correlations with Galactic space motions and H alpha emission line strengths, and return statistically identical metallicities for M dwarfs within a common multiple system. Finally, we find systematic residuals between our H2O-based spectral types and those derived from optical spectral features with previously known sensitivity to stellar metallicity, such as TiO, and identify the CaH1 index as a promising optical index for diagnosing the metallicities of near-solar M dwarfs.


title: Search for Radiative Decays of Cosmic Background Neutrino using Cosmic Infrared Background Energy Spectrum

abstract: We propose to search for the neutrino radiative decay by fitting a photon energy spectrum of the cosmic infrared background to a sum of the photon energy spectrum from the neutrino radiative decay and a continuum. By comparing the present cosmic infrared background energy spectrum observed by AKARI and Spitzer to the photon energy spectrum expected from neutrino radiative decay with a maximum likelihood method, we obatined a lifetime lower limit of $3.1 \times 10^{12}$ to $3.8 \times 10^{12}$ years at 95% confidence level for the third generation neutrino $\nu_3$ in the $\nu_3$ mass range between 50 \mmev and 150 \mmev under the present constraints by the neutrino oscillation measurements. In the left-right symmetric model, the minimum lifetime of $\nu_3$ is predicted to be $1.5 \times 10^{17}$ years for $m_3$ of 50 \mmev. We studied the feasibility of the observation of the neutrino radiative decay with a lifetime of $1.5 \times 10^{17}$ years, by measuring a continuous energy spectrum of the cosmic infrared background.


title: A note on the extended dToda hierarchy

abstract: We give a derivation of dispersionless Hirota equations for the extended dispersionless Toda hierarchy. We show that the dispersionless Hirota equations are nothing but a direct consequence of the genus-zero topological recursion relation for the topological $CP^1$ model. Using the dispersionless Hirota equations we compute the two point functions and express the result in terms of Catalan number.


title: Correcting the apparent mutation rate acceleration at shorter time scales under a Jukes-Cantor model

abstract: At macroevolutionary time scales, and for a constant mutation rate, there is an expected linear relationship between time and the number of inferred neutral mutations (the "molecular clock"). However, at shorter time scales a number of recent studies have observed an apparent acceleration in the rate of molecular evolution. We study this apparent acceleration under a Jukes-Cantor model applied to a randomly mating population, and show that, under the model, it arises as a consequence of ignoring short term effects due to existing diversity within the population. The acceleration can be accounted for by adding the correction term h_0e^{-4mu*t/3} to the usual Jukes-Cantor formula p(t)=(3/4)(1-e^{-4mu*t/3}), where h_0 is the expected heterozygosity in the population at time t=0. The true mutation rate mu may then be recovered, even if h_0 is not known, by estimating mu and h_0 simultaneously using least squares. Rate estimates made without the correction term (that is, incorrectly assuming the population to be homogeneous) will result in a divergent rate curve of the form mu_{div}=mu+C/t, so that the mutation rate appears to approach infinity as the time scale approaches zero. While our quantitative results apply only to the Jukes-Cantor model, it is reasonable to suppose that the qualitative picture that emerges also applies to more complex models. Our study therefore demonstrates the importance of properly accounting for any ancestral diversity, as it may otherwise play a dominant role in rate overestimation.


title: Lower Bounds for Laplacian and Fractional Laplacian Eigenvalues

abstract: In this paper, we investigate eigenvalues of Laplacian on a bounded domain in an $n$-dimensional Euclidean space and obtain a sharper lower bound for the sum of its eigenvalues, which gives an improvement of results due to A. D. Melas [15]. On the other hand, for the case of fractional Laplacian $(-\Delta)^{\alpha/2}|_{D}$, where $\alpha\in(0,2]$, we obtain a sharper lower bound for the sum of its eigenvalues, which gives an improvement of results due to S.Y. Yolcu and T. Yolcu [23].


title: A Constructive Approach to Reduced-Form Auctions with Applications to Multi-Item Mechanism Design

abstract: We provide a constructive proof of Border's theorem [Bor91, HR15a] and its generalization to reduced-form auctions with asymmetric bidders [Bor07, MV10, CKM13]. Given a reduced form, we identify a subset of Border constraints that are necessary and sufficient to determine its feasibility. Importantly, the number of these constraints is linear in the total number of bidder types. In addition, we provide a characterization result showing that every feasible reduced form can be induced by an ex-post allocation rule that is a distribution over ironings of the same total ordering of the union of all bidders' types. We show how to leverage our results for single-item reduced forms to design auctions with heterogeneous items and asymmetric bidders with valuations that are additive over items. Appealing to our constructive Border's theorem, we obtain polynomial-time algorithms for computing the revenue-optimal auction. Appealing to our characterization of feasible reduced forms, we characterize feasible multi-item allocation rules.


title: On the Boundedness of the Carleson Operator near $L^1$

abstract: Based on the tile discretization elaborated by the author in "The Polynomial Carleson Operator", we develop a Calderon-Zygmund type decomposition of the Carleson operator. As a consequence, through a unitary method that makes no use of extrapolation techniques, we recover the previously known results regarding the largest rearrangement invariant space of functions with almost everywhere convergent Fourier series.


title: Notes on the Schreier graphs of the Grigorchuk group

abstract: The paper is concerned with the space of the marked Schreier graphs of the Grigorchuk group and the action of the group on this space. In particular, we describe an invariant set of the Schreier graphs corresponding to the action on the boundary of the binary rooted tree and dynamics of the group action restricted to this invariant set.


title: Pseudodifferential operators on manifolds with fibred corners

abstract: One way to geometrically encode the singularities of a stratified pseudomanifold is to endow its interior with an iterated fibred cusp metric. For such a metric, we develop and study a pseudodifferential calculus generalizing the \Phi-calculus of Mazzeo and Melrose. Our starting point is the observation, going back to Melrose, that a stratified pseudomanifold can be `resolved' into a manifold with fibred corners. This allows us to define pseudodifferential operators as conormal distributions on a suitably blown-up double space. Various symbol maps are introduced, leading to the notion of full ellipticity. This is used to construct refined parametrices and to provide criteria for the mapping properties of operators such as Fredholmness or compactness. We also introduce a semiclassical version of the calculus and use it to establish a Poincar\'e duality between the K-homology of the stratified pseudomanifold and the K-group of fully elliptic operators.


title: Discrete-time Ruijsenaars-Schneider system and Lagrangian 1-form structure

abstract: We study the Lagrange formalism of the (rational) Ruijsenaars-Schneider (RS) system, both in discrete time as well as in continuous time, as a further example of a Lagrange 1-form structure in the sense of the recent paper [24]. The discrete-time model of the RS system was established some time ago arising via an Ansatz of a Lax pair, and was shown to lead to an exactly integrable correspondence (multivalued map)[15]. In this paper we consider an extended system representing a family of commuting flows of this type, and establish a connection with the lattice KP system. In the Lagrangian 1-form structure of this extended model, the closure relation is verified making use of the equations of motion. Performing successive continuum limits on the RS system, we establish the Lagrange 1-form structure for the corresponding continuum case of the RS model.


title: Totally Bipartite/ABipartite Leonard pairs and Leonard triples of Bannai/Ito type

abstract: This paper is about three classes of objects: Leonard pairs, Leonard triples, and the finite-dimensional irreducible modules for an algebra $\mathcal{A}$. Let $\K$ denote an algebraically closed field of characteristic zero. Let $V$ denote a vector space over $\K$ with finite positive dimension. A Leonard pair on $V$ is an ordered pair of linear transformations in End$(V)$ such that for each of these transformations there exists a basis for $V$ with respect to which the matrix representing that transformation is diagonal and the matrix representing the other transformation is irreducible tridiagonal. Whenever the tridiagonal matrices are bipartite, the Leonard pair is said to be totally bipartite. A mild weakening yields a type of Leonard pair said to be totally almost bipartite. A Leonard pair is said to be totally B/AB whenever it is totally bipartite or totally almost bipartite. The notion of a Leonard triple and the corresponding notion of totally B/AB are similarly defined. There are families of Leonard pairs and Leonard triples said to have Bannai/Ito type. The Leonard pairs and Leonard triples of interest to us are totally B/AB and of Bannai/Ito type. Let $\mathcal{A}$ denote the unital associative $\K$-algebra defined by generators $x,y,z$ and relations\[xy+yx=2z,\qquad yz+zy=2x,\qquad zx+xz=2y.\]The algebra $\mathcal{A}$ has a presentation involving generators $x,y$ and relations\[x^{2}y+2xyx+yx^{2}=4y,\qquad y^{2}x+2yxy+xy^{2}=4x.\] In this paper we obtain the following results. We classify up to isomorphism the totally B/AB Leonard pairs of Bannai/Ito type. We classify up to isomorphism the totally B/AB Leonard triples of Bannai/Ito type. We classify up to isomorphism the finite-dimensional irreducible $\mathcal{A}$-modules. We show that these three classes of objects are essentially in one-to-one correspondence, and describe these correspondences in detail.


title: Self-Index based on LZ77 (thesis)

abstract: Domains like bioinformatics, version control systems, collaborative editing systems (wiki), and others, are producing huge data collections that are very repetitive. That is, there are few differences between the elements of the collection. This fact makes the compressibility of the collection extremely high. For example, a collection with all different versions of a Wikipedia article can be compressed up to the 0.1% of its original space, using the Lempel-Ziv 1977 (LZ77) compression scheme. Many of these repetitive collections handle huge amounts of text data. For that reason, we require a method to store them efficiently, while providing the ability to operate on them. The most common operations are the extraction of random portions of the collection and the search for all the occurrences of a given pattern inside the whole collection. A self-index is a data structure that stores a text in compressed form and allows to find the occurrences of a pattern efficiently. On the other hand, self-indexes can extract any substring of the collection, hence they are able to replace the original text. One of the main goals when using these indexes is to store them within main memory. In this thesis we present a scheme for random text extraction from text compressed with a Lempel-Ziv parsing. Additionally, we present a variant of LZ77, called LZ-End, that efficiently extracts text using space close to that of LZ77. The main contribution of this thesis is the first self-index based on LZ77/LZ-End and oriented to repetitive texts, which outperforms the state of the art (the RLCSA self-index) in many aspects. Finally, we present a corpus of repetitive texts, coming from several application domains. We aim at providing a standard set of texts for research and experimentation, hence this corpus is publicly available.


title: Limit Theorems for Quantum Walks on the Union of Planes

abstract: We extend the construction given by [Chisaki et.al, arXiv:1009.1306v1] from lines to planes, and obtain the associated limit theorems for quantum walks on such a graph.


title: Uncertainty relations for cosmological particle creation and existence of large fluctuations in reheating

abstract: We derive an uncertainty relation for the energy density and pressure of a quantum scalar field in a time-dependent, homogeneous and isotropic, classical background, which implies the existence of large fluctuations comparable to their vacuum expectation values. A similar uncertainty relation is known to hold for the field square since the field can be viewed as a Gaussian random variable. We discuss possible implications of these results for the reheating process in scalar field driven inflationary models, where reheating is achieved by the decay of the coherently oscillating inflaton field. Specifically we argue that the evolution after backreaction can seriously be altered by the existence of these fluctuations. For example, in one model the coherence of the inflaton oscillations is found to be completely lost in a very short time after backreaction starts. Therefore we argue that entering a smooth phase in thermal equilibrium is questionable in such models and reheating might destroy the smoothness attained by inflation.


title: Limits on the Gas Disk Content of Two "Evolved" T Tauri Stars

abstract: We derived upper limits of the circumstellar gas disk masses around the T Tauri stars St 34 and RX J0432.8+1735 in order to place constraints on theories of planet formation and to explore the evolution of the gas-to-dust ratio during the epoch of disk dissipation around young sun-like stars. Since sub-millimeter lines of ^{12}CO trace of the cold, outer regions of circumstellar disks, we observed ^{12}CO J=2-1 emission with the 10 m Sub-Millimeter Telescope (SMT) for two carefully chosen targets. St 34 is a rare classical T Tauri star with an age of 8\pm3 Myr, and RX J0432.8+1735 is a rare weak-emission T Tauri star with far-infrared excess. Both exhibit radial space motion enabling us to distinguish disk emission from ambient cloud material. Assuming a ^{12}CO excitation temperature of 20 K, a ^{12}CO line-width of 5 km s^{-1}, and optically-thin emission, we derive 3{\sigma} upper limits on the H_{2} circumstellar disk mass for St 34 and RX J0432.8+1735 to be <4.20 M\odot for both disks. Placing these results in the context of other studies, we discuss their implications on planet formation models.


title: Phase transitions for random states and a semi-circle law for the partial transpose

abstract: For a system of N identical particles in a random pure state, there is a threshold k_0 = k_0(N) ~ N/5 such that two subsystems of k particles each typically share entanglement if k > k_0, and typically do not share entanglement if k < k_0. By "random" we mean here "uniformly distributed on the sphere of the corresponding Hilbert space." The analogous phase transition for the positive partial transpose (PPT) property can be described even more precisely. For example, for N qubits the two subsystems of size k are typically in a PPT state if k < k_1 := N/4 - 1/2, and typically in a non-PPT state if k > k_1. Since, for a given state of the entire system, the induced state of a subsystem is given by the partial trace, the above facts can be rephrased as properties of random induced states. An important step in the analysis depends on identifying the asymptotic spectral density of the partial transposes of such random induced states, a result which is interesting in its own right.


title: Multi-Wavelength Photometric and Polarimetric Observations of the Outburst of 3C 454.3 in Dec. 2009

abstract: In December 2009, the bright blazar, 3C 454.3 exhibited a strong outburst in the optical, X-ray and gamma-ray regions. We performed photometric and polarimetric monitoring of this outburst in the optical and near-infrared bands with TRISPEC and HOWPol attached to the Kanata telescope. We also observed this outburst in the infrared band with AKARI, and the radio band with the 32-m radio telescope of Yamaguchi University. The object was in an active state from JD 2455055 to 2455159. It was 1.3 mag brighter than its quiescent state before JD 2455055 in the optical band. After the end of the active state in JD 2455159, a prominent outburst was observed in all wavelengths. The outburst continued for two months. Our optical and nearinfrared polarimetric observations revealed that the position angle of the polarization (PA) apparently rotated clockwise by 240 degrees within 11 d in the active state (JD 2455063-2455074), and after this rotation, PA remained almost constant during our monitoring. In the outburst state, PA smoothly rotated counterclockwise by 350 degrees within 35 d (JD 2455157-2455192). Thus, we detected two distinct rotation events of polarization vector in opposite directions. We discuss these two events compared with the past rotation events observed in 2005, 2007 and 2008.


title: Equivariant semiprojectivity

abstract: We define equivariant semiprojectivity for C*-algebras equipped with actions of compact groups. We prove that the following examples are equivariantly semiprojective: arbitrary finite dimensional C*-algebras with arbitrary actions of compact groups; the Cuntz algebras ${\mathcal{O}}_d$ and extended Cuntz algebras $E_d,$ for finite $d,$ with quasifree actions of compact groups; the Cuntz algebra ${\mathcal{O}}_{\infty}$ with any quasifree action of a finite group. For actions of finite groups, we prove that equivariant semiprojectivity is equivalent to a form of equivariant stability of generators and relations. We also prove that if $G$ is finite, then $C^* (G)$ is graded semiprojective.


title: The dark side of benzene: interference vs. interaction

abstract: We present the study of the linear conductance vs.\ applied gate voltage for an interacting six site ring structure, which is threaded by a flux of $\pi$ and coupled to a left and a right lead. This ring structure is designed to have a vanishing conductance for all gate voltages and temperatures provided interactions are ignored. Therefore this system is an ideal testbed to study the interplay of interaction and interference. First we find a Kondo type resonance for rather large hopping parameter. Second, we find additional resonance peaks which can be explained by a population blocking mechanism. To this end we have to extend the Kubo approach within the Density Matrix Renormalization Group method to handle degenerate states.


title: Shear viscosity and the nucleation of antikaon condensed matter in protoneutron stars

abstract: We study shear viscosities of different species in hot and neutrino-trapped dense matter relevant to protoneutron stars. It is found that the shear viscosities of neutrons, protons and electrons in neutrino-trapped matter are of the same orders of magnitude as the corresponding shear viscosities in neutrino-free matter. Above all, the shear viscosity due to neutrinos is higher by several orders of magnitude than that of other species in neutrino-trapped matter. Next we investigate the effect of shear viscosity in particular, neutrino shear viscosity on the thermal nucleation rate of droplets of antikaon condensed matter in protoneutron stars. The first-order phase transition from hadronic matter to antikaon condensed matter is driven by the thermal nucleation process. We compute the equation of state used for the calculation of shear viscosity and thermal nucleation time within the relativistic mean field model. Neutrino shear viscosity enhances the prefactor in the nucleation rate by several orders of magnitude compared with the $T^4$ approximation of earlier calculations. Consequently the thermal nucleation time in the $T^4$ approximation overestimates our result. Furthermore, the thermal nucleation of an antikaon droplet might be possible in neutrino-trapped matter before neutrino diffusion takes place.


title: Third order operator with periodic coefficients on the real line

abstract: We consider the third order operator with periodic coefficients on the real line. This operator is used in the integration of the non-linear evolution Boussinesq equation. For the minimal smoothness of the coefficients we prove that: 1) the operator is self-adjoint and it is decomposable into the direct integral, 2) the spectrum is absolutely continuous, fills the whole real axis, and has multiplicity one or three, 3) the Lyapunov function, analytic on a three-sheeted Riemann surface, is constructed and researched, 4) the spectrum of multiplicity three is bounded and it is described in terms of some entire function (the discriminant).


title: Sub-cycle waveform synthesis for effective control of electron localization

abstract: We have proposed and demonstrated an effective scheme to steer both the electronic and molecular motions of hydrogen molecule and its isotopes. Our scheme lies in the idea of sub-cycle waveform shaping (SWS) with synthesizing coherent IR pulses of different wavelengths. The quantum model simulation indicates a remarkable enhancement of the electron localization asymmetry compared with the few-cycle pulse. 68% and 89% electrons of the dissociating fragmentation can be localized to one nucleus for H_2^+ and H_2, respectively.


title: An Information Theoretical Analysis of Kinase Activated Phosphorylation Dephosphorylation Cycle

abstract: Signal transduction, the information processing mechanism in biological cells, is carried out by a network of biochemical reactions. The dynamics of driven biochemical reactions can be studied in terms of nonequilibrium statistical physics. Such systems may also be studied in terms of Shannon's information theory. We combine these two perspectives in this study of the basic units (modules) of cellular signaling: the phosphorylation dephosphorylation cycle (PdPC) and the guanosine triphosphatase (GTPase). We show that the channel capacity is zero if and only if the free energy expenditure of biochemical system is zero. In fact, a positive correlation between the channel capacity and free energy expenditure is observed. In terms of the information theory, a linear signaling cascade consisting of multiple steps of PdPC can function as a distributed "multistage code". With increasing number of steps in the cascade, the system trades channel capacity with the code complexity. Our analysis shows that while a static code can be molecular structural based; a biochemical communication channel has to have energy expenditure.


title: Generalized repeated interaction model and transfer functions

abstract: Using a scheme involving a lifting of a row contraction we introduce a toy model of repeated interactions between quantum systems. In this model there is an outgoing Cuntz scattering system involving two wandering subspaces. We associate to this model an input/output linear system which leads to a transfer function. This transfer function is a multi-analytic operator, and we show that it is inner if we assume that the system is observable. Finally it is established that transfer functions coincide with characteristic functions of associated liftings.


title: Parametrization of coarse grained force fields for dynamic property of ethylene glycol oligomers/water binary mixtures

abstract: To evaluate shear viscosity of ethylene glycol oligomers (EGO)/water binary mixture by means of coarse-grained molecular dynamics (CG-MD) simulations, we proposed the self-diffusion-coefficient-based parameterization of non-bonded interactions among CG particles. Our parameterization procedure consists of three steps: 1) determination of bonded potentials, 2) scaling for time and solvent diffusivity, and 3) optimization of Lennard-Jones parameters to reproduce experimental self-diffusion coefficient and density data. With the determined parameters and the scaling relations, we evaluated shear viscosities of aqueous solutions of EGOs with various molecular weights and concentrations. Our simulation result are in close agreement with the experimental data. The largest simulation in this article corresponds to a 1.2 microseconds atomistic simulation for 100,000 atoms. Our CG model with the parameterization scheme for CG particles may be useful to study the dynamic properties of a liquid which contains relatively low molecular weight polymers or oligomers.

categories: math-ph cond-mat.soft math.MP physics.comp-ph


title: Measurement of charged particle multiplicities in $pp$ collisions at ${\sqrt{s} =7}$TeV in the forward region

abstract: The charged particle production in proton-proton collisions is studied with the LHCb detector at a centre-of-mass energy of ${\sqrt{s} =7}$TeV in different intervals of pseudorapidity $\eta$. The charged particles are reconstructed close to the interaction region in the vertex detector, which provides high reconstruction efficiency in the $\eta$ ranges $-2.5<\eta<-2.0$ and $2.0<\eta<4.5$. The data were taken with a minimum bias trigger, only requiring one or more reconstructed tracks in the vertex detector. By selecting an event sample with at least one track with a transverse momentum greater than 1 GeV/c a hard QCD subsample is investigated. Several event generators are compared with the data; none are able to describe fully the multiplicity distributions or the charged particle density distribution as a function of $\eta$. In general, the models underestimate the charged particle production.


title: Fluid dynamics of dilatant fluid

abstract: Dense mixture of granules and liquid often shows a sever shear thickening and is called a dilatant fluid. We construct a fluid dynamics model for the dilatant fluid by introducing a phenomenological state variable for a local state of dispersed particles. With simple assumptions for an equation of the state variable, we demonstrate that the model can describe basic features of the dilatant fluid such as the stress-shear rate curve that represents discontinuous severe shear thickening, hysteresis upon changing shear rate, instantaneous hardening upon external impact. Analysis of the model reveals that the shear thickening fluid shows an instability in a shear flow for some regime and exhibits {\it the shear thickening oscillation}, i.e. the oscillatory shear flow alternating between the thickened and the relaxed states. Results of numerical simulations are presented for one and two-dimensional systems.


title: Atomic homodyne detection of continuous variable entangled twin-atom states

abstract: Historically, the completeness of quantum theory has been questioned using the concept of bipartite continuous variable entanglement. The non-classical correlations (entanglement) between the two subsystems imply that the observables of one subsystem are determined by the measurement choice on the other, regardless of their distance. Nowadays, continuous variable entanglement is regarded as an essential resource allowing for quantum enhanced measurement resolution, the realization of quantum teleportation and quantum memories, or the demonstration of the Einstein-Podolsky-Rosen paradox. These applications rely on techniques to manipulate and detect coherences of quantum fields, the quadratures. While in optics coherent homodyne detection of quadratures is a standard technique, for massive particles a corresponding method was missing. Here we report on the realization of an atomic analog to homodyne detection for the measurement of matter-wave quadratures. The application of this technique to a quantum state produced by spin-changing collisions in a Bose-Einstein condensate reveals continuous variable entanglement, as well as the twin-atom character of the state. With that we present a new system in which continuous variable entanglement of massive particles is demonstrated. The direct detection of atomic quadratures has applications not only in experimental quantum atom optics but also for the measurement of fields in many-body systems of massive particles.


title: Extension of the Schiff theorem to ions and molecules

abstract: According to the Schiff theorem the nuclear electric dipole moment (EDM) is screened in neutral atoms. In ions this screening is incomplete. We extend a derivation of the Schiff theorem to ions and molecules. The finite nuclear size effects are considered including Z^2 alpha^2 corrections to the nuclear Schiff moment which are significant in all atoms and molecules of experimental interest. We show that in majority of ionized atoms the nuclear EDM contribution to the atomic EDM dominates while in molecules the contribution of the Schiff moment dominates. We also consider the screening of electron EDM in ions.


title: Discovery of the Youngest Molecular Outflow associated with an Intermediate-mass protostellar Core, MMS-6/OMC-3

abstract: We present sub-arcsecond resolution HCN (4-3) and CO (3-2) observations made with the Submillimeter Array (SMA), toward an extremely young intermediate-mass protostellar core, MMS 6-main, located in the Orion Molecular Cloud 3 region (OMC-3). We have successfully imaged a compact molecular outflow lobe (~1500 AU) associated with MMS6-main, which is also the smallest molecular outflow ever found in the intermediate-mass protostellar cores. The dynamical time scale of this outflow is estimated to be <100 yr. The line width dramatically increases downstream at the end of the molecular outflow ({\Delta}v~25 km s^{-1}), and clearly shows the bow-shock type velocity structure. The estimated outflow mass (~10^{-4} M_{sun}) and outflow size are approximately 2-4 orders and 1-3 orders of magnitude smaller, while the outflow force (~10^{-4} M_{sun} km s^{-1} yr^{-1}) is similar, as compared to the other molecular outflows studied in OMC-2/3. These results show that MMS 6-main is a protostellar core at the earliest evolutionary stage, most likely shortly after the 2nd core formation.


title: Evaluating Network Models: A Likelihood Analysis

abstract: Many models are put forward to mimic the evolution of real networked systems. A well-accepted way to judge the validity is to compare the modeling results with real networks subject to several structural features. Even for a specific real network, we cannot fairly evaluate the goodness of different models since there are too many structural features while there is no criterion to select and assign weights on them. Motivated by the studies on link prediction algorithms, we propose a unified method to evaluate the network models via the comparison of the likelihoods of the currently observed network driven by different models, with an assumption that the higher the likelihood is, the better the model is. We test our method on the real Internet at the Autonomous System (AS) level, and the results suggest that the Generalized Linear Preferential (GLP) model outperforms the Tel Aviv Network Generator (Tang), while both two models are better than the Barab\'asi-Albert (BA) and Erd\"os-R\'enyi (ER) models. Our method can be further applied in determining the optimal values of parameters that correspond to the maximal likelihood. Experiment indicates that the parameters obtained by our method can better capture the characters of newly-added nodes and links in the AS-level Internet than the original methods in the literature.


title: The Formation of Supermassive Black Holes from Low-Mass Pop III Seeds

abstract: The existence of 10$^9$ M$_{\odot}$ black holes (BH) in massive galaxies by $z \ sim$ 7 is one of the great unsolved mysteries in cosmological structure formation. One theory argues that they originate from the black holes of Pop III stars at $z \sim$ 20 and then accrete at the Eddington limit down to the epoch of reionization, which requires that they have constant access to rich supplies of fuel. Because early numerical simulations suggested that Pop III stars were $\gtrsim$ 100 M$_{\odot}$, the supermassive black hole seeds considered up to now were 100 - 300 M$_{\odot}$. However, there is a growing numerical and observational consensus that some Pop III stars were tens of solar masses, not hundreds, and that 20 - 40 M$_{\odot}$ black holes may have been much more plentiful at high redshift. However, we find that natal kicks imparted to 20 - 40 M$_{\odot}$ Pop III BHs during formation eject them from their halos and hence their fuel supply, precluding them from Eddington-limit growth. Consequently, supermassive black holes are far less likely to form from low-mass Pop III stars than from very massive ones.


title: Strings from geometric tachyon in Rindler space and black hole thermodynamics

abstract: The dynamics of a probe particle or wrapped brane moving in the two-dimensional Rindler space can be described by a time-dependent tachyon field theory. Using knowledge of tachyon condensation, we learn that the infalling brane gets thermalised and produces open string pairs at the Hagedorn temperature when entering into the near-horizon Rindler wedge. It is shown that the Hagedorn temperature of the infalling brane is equal to the Hawking temperature of the host black hole detected in the same time coordinate. The infalling brane will decay completely into closed strings, mainly massive modes, when it reaches the horizon in infinitely long time as observed by observers at spatial infinity. Preliminary estimates indicate that the degeneracy of states of the closed strings emitted from the infalling brane should be responsible for the increased entropy in the host black hole due to absorption of the brane.


title: Mass spectrum of the minimal SUSY B-L model

abstract: The origin of R-parity in supersymmetric models can be explained if \BL is part of the gauge group. We discuss the mass spectrum of the minimal $U(1)_Y \times U(1)_{B-L}$ model based on a GUT implementation using CMSSM-like boundary conditions. Here we focus in particular on the Higgs and neutralino sectors in this class of models. While the neutralinos can have masses as low as 100 GeV, we show that the requirement of being consistent with existing bounds on the $Z'$ implies that in general the sfermions have masses in the multi-TeV range. In the extended Higgs sector we show the existence of a second light state which, however, will be difficult to observe, while having at the same time a SM-like Higgs in a mass range of 123-126 GeV.. Moreover, we propose a set of benchmark scenarios for phenomenological studies. On the technical side we demonstrate that gauge kinetic mixing effects can be quite important, affecting in particular the Higgs and the neutralino sectors. Not only can they shift the mass of the lightest neutralino by about 10 per-cent but also they can change the nature of neutralinos and Higgs bosons in a significant way.


title: Formal pseudodifferential operators and Witten's r-spin numbers

abstract: We derive an effective recursion for Witten's r-spin intersection numbers, using Witten's conjecture relating r-spin numbers to the Gel'fand-Dikii hierarchy (Theorem 4.1). Consequences include closed-form descriptions of the intersection numbers (for example, in terms of gamma functions: Propositions 5.2 and 5.4, Corollary 5.5). We use these closed-form descriptions to prove Harer-Zagier's formula for the Euler characteristic of M_{g,1}. Finally in Section 6, we extend Witten's series expansion formula for the Landau-Ginzburg potential to study r-spin numbers in the small phase space in genus zero. Our key tool is the calculus of formal pseudodifferential operators, and is partially motivated by work of Brezin and Hikami.


title: A time-variable, phase-dependent emission line in the X-ray spectrum of the isolated neutron star RX J0822-4300

abstract: RX J0822-4300 is the Central Compact Object associated with the Puppis A supernova remnant. Previous X-ray observations suggested RX J0822-4300 to be a young neutron star with a weak dipole field and a peculiar surface temperature distribution dominated by two antipodal spots with different temperatures and sizes. An emission line at 0.8 keV was also detected. We performed a very deep (130 ks) observation with XMM-Newton, which allowed us to study in detail the phase-resolved properties of RX J0822-4300. Our new data confirm the existence of a narrow spectral feature, best modelled as an emission line, only seen in the `Soft' phase interval - when the cooler region is best aligned to the line of sight. Surprisingly, comparison of our recent observations to the older ones yields evidence for a variation in the emission line component, which can be modelled as a decrease in the central energy from ~0.80 keV in 2001 to ~0.73 keV in 2009--2010. The line could be generated via cyclotron scattering of thermal photons in an optically thin layer of gas, or - alternatively - it could originate in low-rate accretion by a debris disk. In any case, a variation in energy, pointing to a variation of the magnetic field in the line emitting region, cannot be easily accounted for.