title: The fraction of muon tracks in cosmic neutrinos

abstract: The study of the distintive signatures of the ultra high energy events recently seen by IceCube can allow to single the neutrino origin out. The detection of tau neutrinos would be a clear way to prove that they come from cosmic distances, but at the highest energies currently seen, about 1 PeV, an experimental characterization of tau events is difficult. The study of the fraction of the muon tracks seems more promising. In fact, for any initial composition, because of the occurrence of flavor oscillations and despite their uncertainties, the fraction of muon tracks in the cosmic neutrinos is smaller than the one of atmospheric neutrinos, even hypothesizing an arbitrarily large contribution from charmed mesons. A good understanding of the detection efficiencies and the optimization of the analysis cuts, along with a reasonable increase in the statistics, should provide us with a significant test of the cosmic origin of these events.


title: On existence and properties of strong solutions of one-dimensional stochastic equations with an additive noise

abstract: One-dimensional stochastic differential equations with additive L\'evy noise are considered. Conditions for existence and uniqueness of a strong solution are obtained. In particular, if the noise is a L\'evy symmetric stable process with $\alpha\in(1;2)$, then the measurability and boundedness of a drift term is sufficient for the existence of a strong solution. We also study continuous dependence of the strong solution on the initial value and the drift.


title: Derivation of the coefficient squared probability law in quantum mechanics

abstract: If one assumes there is probability of perception in quantum mechanics, then unitarity dictates that it must have the coefficient squared form, in agreement with experiment.


title: Equation of State and Viscosities from a Gravity Dual of the Gluon Plasma

abstract: Employing new precision data of the equation of state of the SU(3) Yang-Mills theory (gluon plasma) the dilaton potential of a gravity-dual model is adjusted in the temperature range $(1 - 10) T_c$ within a bottom-up approach. The ratio of bulk viscosity to shear viscosity follows then as $\zeta/\eta \approx \pi \Delta v_s^2$ for $\Delta v_s^2 < 0.2$ and achieves a maximum value of $0.94$ at $\Delta v_s^2 \approx 0.3$, where $\Delta v_s^2 \equiv 1/3 - v_s^2$ is the non-conformality measure and $v_s^2$ is the velocity of sound squared, while the ratio of shear viscosity to entropy density is known as $(4 \pi)^{-1}$ for the considered set-up with Hilbert action on the gravity side.


title: Cross-border Portfolio Investment Networks and Indicators for Financial Crises

abstract: Cross-border equity and long-term debt securities portfolio investment networks are analysed from 2002 to 2012, covering the 2008 global financial crisis. They serve as network-proxies for measuring the robustness of the global financial system and the interdependence of financial markets, respectively. Two early-warning indicators for financial crises are identified: First, the algebraic connectivity of the equity securities network, as a measure for structural robustness, drops close to zero already in 2005, while there is an over-representation of high-degree off-shore financial centres among the countries most-related to this observation, suggesting an investigation of such nodes with respect to the structural stability of the global financial system. Second, using a phenomenological model, the edge density of the debt securities network is found to describe, and even forecast, the proliferation of several over-the-counter-traded financial derivatives, most prominently credit default swaps, enabling one to detect potentially dangerous levels of market interdependence and systemic risk.


title: Circular unitary ensemble with highly oscillatory potential

abstract: We study the effect of highly oscillatory potentials to the eigenvalues of a random matrix. Consider the circular unitary ensembles with an external potential which is periodic with the period comparable to the average spacing of the eigenvalues. We show that in this case the density of states is periodic and does not converge in the large matrix limit, but the local correlation functions converge to some simple combinations of the sine kernel and the potential. We evaluate the correlation functions exactly and also asymptotically.


title: Eigendecomposition of Block Tridiagonal Matrices

abstract: Block tridiagonal matrices arise in applied mathematics, physics, and signal processing. Many applications require knowledge of eigenvalues and eigenvectors of block tridiagonal matrices, which can be prohibitively expensive for large matrix sizes. In this paper, we address the problem of the eigendecomposition of block tridiagonal matrices by studying a connection between their eigenvalues and zeros of appropriate matrix polynomials. We use this connection with matrix polynomials to derive a closed-form expression for the eigenvectors of block tridiagonal matrices, which eliminates the need for their direct calculation and can lead to a faster calculation of eigenvalues. We also demonstrate with an example that our work can lead to fast algorithms for the eigenvector expansion for block tridiagonal matrices.


title: The Martingale Property in the Context of Stochastic Differential Equations

abstract: This note studies the martingale property of a nonnegative, continuous local martingale Z, given as a nonanticipative functional of a solution to a stochastic differential equation. The condition states that Z is a (uniformly integrable) martingale if and only if an integral test of a related functional holds.


title: Constant mean curvature $k$-noids in homogeneous manifolds

abstract: For each $k\geq2$, we construct two families of surfaces with constant mean curvature $H$ for $H\in[0,1/2]$ in $\Sigma(\kappa)\times\R$ where $\kappa+4H^2\leq0$. The surfaces are invariant under $2\pi/k$-rotations about a vertical fiber of $\Sigma(\kappa)\times\R$, have genus zero, and a finite number of ends. The first family generalizes the notion of $k$-noids: It has $k$ ends, one horizontal and $k$ vertical symmetry planes. The second family is less symmetric and has two types of ends. Each surface arises as the conjugate (sister) surface of a minimal graph in a homogeneous 3-manifold. The domain of the graph is non-convex in the second family. For $\kappa=-1$ the surfaces with constant mean curvature $H$ arise from a minimal surface in $\widetilde{\PSL}_2(\R)$ for $H\in(0,1/2)$ and in $\Nil$ for H=1/2. For H=0, the conjugate surfaces are both minimal in a product space.


title: Stringy Membranes in AdS/CFT

abstract: We study membrane configurations in AdS_{7/4}xS^{4/7}. The membranes are wrapped around the compact manifold S^{4/7} and are dynamically equivalent to bosonic strings in AdS_5. We thus conveniently identify them as "Stringy Membranes". For the case of AdS_7xS^4, their construction is carried out by embedding the Polyakov action for classical bosonic strings in AdS_5, into the corresponding membrane action. Therefore, every string configuration in AdS_5 can be realized by an appropriately chosen stringy membrane in AdS_7xS^4. We discuss the possibility of this being also the case for stringy membranes in AdS_4xS^7/Z^k (k > 1 or k = 1). By performing a stability analysis to the constructed solutions, we find that the (membrane) fluctuations along their transverse directions are organized in multiple Lam\'{e} stability bands and gaps in the space of parameters of the configurations. In this membrane picture, strings exhibit a single band/gap structure.


title: Survey on QoE\QoS Correlation Models For Multimedia Services

abstract: This paper presents a brief review of some existing correlation models which attempt to map Quality of Service (QoS) to Quality of Experience (QoE) for multimedia services. The term QoS refers to deterministic network behaviour, so that data can be transported with a minimum of packet loss, delay and maximum bandwidth. QoE is a subjective measure that involves human dimensions; it ties together user perception, expectations, and experience of the application and network performance. The Holy Grail of subjective measurement is to predict it from the objective measurements; in other words predict QoE from a given set of QoS parameters or vice versa. Whilst there are many quality models for multimedia, most of them are only partial solutions to predicting QoE from a given QoS. This contribution analyses a number of previous attempts and optimisation techniquesthat can reliably compute the weighting coefficients for the QoS/QoE mapping.


title: Hybrid Optical and Electrical Network Flows Scheduling in Cloud Data Centres

abstract: Hybrid intra-data centre networks, with optical and electrical capabilities, are attracting research interest in recent years. This is attributed to the emergence of new bandwidth greedy applications and novel computing paradigms. A key decision to make in networks of this type is the selection and placement of suitable flows for switching in circuit network. Here, we propose an efficient strategy for flow selection and placement suitable for hybrid Intra-cloud data centre networks. We further present techniques for investigating bottlenecks in a packet networks and for the selection of flows to switch in circuit network. The bottleneck technique is verified on a Software Defined Network (SDN) testbed. We also implemented the techniques presented here in a scalable simulation experiment to investigate the impact of flow selection on network performance. Results obtained from scalable simulation experiment indicate a considerable improvement on average throughput, lower configuration delay, and stability of offloaded flows


title: Competing 4f-electron dynamics in Ce(Ru1-xFex)2Al10 (x=0 to 1): magnetic ordering emerging from the Kondo semiconducting state

abstract: We have carried out muon spin relaxation (muSR), neutron diffraction and inelastic neutron scattering (INS) investigations on polycrystalline samples of Ce(Ru1-xFex)2Al10 (x=0, 0.3, 0.5, 0.8 and 1) to investigate the nature of the ground state (magnetic ordered versus paramagnetic) and the origin of the spin gap formation as evident from the bulk measurements in the end members. Our zero-field muSR spectra clearly reveal coherent two-frequency oscillations at low temperature in x=0, 0.3 and 0.5 samples, which confirms the long-range magnetic ordering of the Ce-moment with TN=27, 26 and 21 K respectively. On the other hand the muSR spectra of x=0.8 and x=1 down to 1.4 K and 0.045 K, respectively exhibit a temperature independent Kubo-Toyabe term confirming a paramagnetic ground state. The long-range magnetic ordering in x=0.5 below 21 K has been confirmed through the neutron diffraction study. INS measurements of x=0 clearly reveal the presence of a sharp inelastic excitation near 8 meV between 5 K and 26 K, due to an opening of a gap in the spin excitation spectrum, which transforms into a broad response at and above 30 K. Interestingly, at 4.5 K the spin gap excitation broadens in x=0.3 and exhibits two clear peaks at 8.4(3) and 12.0(5) meV in x=0.5. In the x=0.8 sample, which remains paramagnetic down to 1.2 K, there is a clear signature of a spin gap of 10-12 meV at 7 K, with a strong Q-dependent intensity. Evidence of a spin gap of 12.5(5) meV has also been found in x=1. The observation of a spin gap in the paramagnetic samples (x=0.8 and 1) is an interesting finding in this study and it challenges our understanding of the origin of the semiconducting gap in CeT2Al10 (T=Ru and Os) compounds in terms of hybridization gap opening only a small part of the Fermi surface, gapped spin waves, or a spin-dimer gap.


title: q-symmetric functions and q-quasisymmetric functions

abstract: In this paper, we construct the q-analogue of Poirier-Reutenauer algebras, related deeply with other q-combinatorial Hopf algebras. As an application, we use them to realize the odd Schur functions defined in \cite{EK}, and naturally obtain the odd Littlewood-Richardson rule concerned in \cite{Ell}. Moreover, we construct the refinement of the odd Schur functions, called odd quasisymmetric Schur functions, parallel to the consideration in \cite{HLMW1}. All the q-Hopf algebras we discuss here provide the corresponding q-dual graded graphs in the sense of \cite{BLL}.


title: Convergence Analysis and Parallel Computing Implementation for the Multiagent Coordination Optimization Algorithm

abstract: In this report, a novel variation of Particle Swarm Optimization (PSO) algorithm, called Multiagent Coordination Optimization (MCO), is implemented in a parallel computing way for practical use by introducing MATLAB built-in function "parfor" into MCO. Then we rigorously analyze the global convergence of MCO by means of semistability theory. Besides sharing global optimal solutions with the PSO algorithm, the MCO algorithm integrates cooperative swarm behavior of multiple agents into the update formula by sharing velocity and position information between neighbors to improve its performance. Numerical evaluation of the parallel MCO algorithm is provided in the report by running the proposed algorithm on supercomputers in the High Performance Computing Center at Texas Tech University. In particular, the optimal value and consuming time are compared with PSO and serial MCO by solving several benchmark functions in the literature, respectively. Based on the simulation results, the performance of the parallel MCO is not only superb compared with PSO for solving many nonlinear, noncovex optimization problems, but also is of high efficiency by saving the computational time.


title: The 617 MHz - $\lambda$ 850 $\mu$m Correlation (Cosmic Rays and Cold Dust) in NGC 3044 and NGC 4157

abstract: We present the first maps of NGC 3044 and NGC 4157 at $\lambda\,450 \mu$m and $\lambda\,850 \mu$m from the JCMT as well as the first maps at 617 MHz from the GMRT. High latitude emission has been detected in both the radio continuum and sub-mm for NGC 3044 and in the radio continuum for NGC 4157, including several new features. The dust spectrum at long wavelengths required fitting with a two-temperature model for both galaxies, implying the presence of cold dust. Dust masses are $M_d\,=\,1.6\,\times\,10^8 M_\odot$ and $M_d\,=\,2.1\,\times\,10^7 M_\odot$ for NGC 3044 and NGC 4157, respectively, and are dominated by the cold component. There is a clear correlation between the 617 MHz and $\lambda\,850 \mu$m emission in the two galaxies. In the case of NGC 3044, this implies a relation between the non-thermal synchrotron emission and cold dust. The 617 MHz component represents an integration of massive star formation over the past $10^{7-8}$ yr and the $\lambda \,850 \mu$m emission represents heating from the diffuse interstellar radiation field (ISRF). The 617 MHz -- $\lambda\,850 \mu$m correlation improves when a smoothing kernel is applied to the $\lambda\,850 \mu$m data to account for differences between the CR electron diffusion scale and the mean free path of an ISRF photon to dust. The best-fit relation is $L_{617_{\rm MHz}}\,\propto\,{L_{850 \mu{\rm m}}}^{2.1\,\pm\,0.2}$ for NGC 3044. If variations in the cold dust emissivity are dominated by variations in dust density, and the synchrotron emission depends on magnetic field strength (a function of gas density) as well as CR electron generation (a function of massive star formation rate and therefore density via the Schmidt law) then the expected correlation for NGC 3044 is $L_{617_{\rm MHz}}\,\propto\,{L_{850 \mu{\rm m}}}^{2.2}$, in agreement with the observed correlation.


title: High Order Asymptotic Preserving DG-IMEX Schemes for Discrete-Velocity Kinetic Equations in a Diffusive Scaling

abstract: In this paper, we develop a family of high order asymptotic preserving schemes for some discrete-velocity kinetic equations under a diffusive scaling, that in the asymptotic limit lead to macroscopic models such as the heat equation, the porous media equation, the advection-diffusion equation, and the viscous Burgers equation. Our approach is based on the micro-macro reformulation of the kinetic equation which involves a natural decomposition of the equation to the equilibrium and non-equilibrium parts. To achieve high order accuracy and uniform stability as well as to capture the correct asymptotic limit, two new ingredients are employed in the proposed methods: discontinuous Galerkin spatial discretization of arbitrary order of accuracy with suitable numerical fluxes; high order globally stiffly accurate implicit-explicit Runge-Kutta scheme in time equipped with a properly chosen implicit-explicit strategy. Formal asymptotic analysis shows that the proposed scheme in the limit of epsilon -> 0 is an explicit, consistent and high order discretization for the limiting equation. Numerical results are presented to demonstrate the stability and high order accuracy of the proposed schemes together with their performance in the limit.


title: Quantum discord for two-qubit X states: Analytical formula with very small worst-case error

abstract: Quantum discord is a measure of quantum correlation beyond entanglement. Computing quantum discord for simple quantum states is a basic problem. An analytical formula of quantum discord for two-qubit X states is first claimed in [Ali, Rau, and Alber, Phys. Rev. A 81, 042105 (2010)], but later found to be not always correct. I observe numerically that the formula is valid with worst-case absolute error 0.0021. For symmetric two-qubit X states, I give a counterexample to the analytical formula derived in [F. F. Fanchini et al., Phys. Rev. A 81, 052107 (2010)], but observe that the formula is valid with worst-case absolute error 0.0006. The formula has been used in many research papers. The results in all these works are approximately correct, even if they may not be exactly correct.


title: Characterization of the degree of Musical non-Markovianity

abstract: Recently, as an aid for musical analysis, in computational musicology mathematical and informatics tools have been developed to quantitatively characterize some aspects of musical compositions. To a musical composition can be attributed by ear a certain amount of memory. This results associated to repetitions and similarities of the patterns in musical scores. To higher variations, a lower amount of memory is perceived. However musical memory of a score has never been quantitatively defined. Here we aim to give such a measure following an approach similar to that used in Physics to quantify the memory (non-Markovianity) of open quantum systems. We have applied this measure to some existing musical compositions, showing that the results obtained via this quantifier agree with what one expects by ear. The Musical non-Markovianity quantifier can thus be used as a new tool that can aid quantitative musical analysis.


title: Spin relaxation mechanism in graphene: resonant scattering by magnetic impurities

abstract: It is proposed that the observed small (100 ps) spin relaxation time in graphene is due to resonant scattering by local magnetic moments. At resonances, magnetic moments behave as spin hot spots: the spin-flip scattering rates are as large as the spin-conserving ones, as long as the exchange interaction is greater than the resonance width. Smearing of the resonance peaks by the presence of electron-hole puddles gives quantitative agreement with experiment, for about 1 ppm of local moments. While the local moments can come from a variety of sources, we specifically focus on hydrogen adatoms. We perform first-principles supercell calculations and introduce an effective Hamiltonian to obtain realistic input parameters for our mechanism.


title: An improved algorithm and a Fortran 90 module for computing the conical function $P^m_{-1/2+i\tau}(x)$

abstract: In this paper we describe an algorithm and a Fortran 90 module ({\bf Conical}) for the computation of the conical function $P^m_{-\tfrac12+i\tau}(x)$ for $x>-1$, $m \ge 0$, $\tau >0$. These functions appear in the solution of Dirichlet problems for domains bounded by cones; because of this, they are involved in a large number of applications in Engineering and Physics. In the Fortran 90 module, the admissible parameter ranges for computing the conical functions in standard IEEE double precision arithmetic are restricted to $(x,m,\tau) \in (-1,1) \times [0,\,40] \times [0,\,100]$ and $(x,m,\tau) \in (1,100) \times [0,\,100] \times [0,\,100]$. Based on tests of the three-term recurrence relation satisfied by these functions and direct comparison with Maple, we claim a relative accuracy close to $10^{-12}$ in the full parameter range, although a mild loss of accuracy can be found at some points of the oscillatory region of the conical functions. The relative accuracy increases to $10^{-13}\,-\,10^{-14}$ in the region of the monotonic regime of the functions where integral representations are computed ($-1<x<0$).


title: Fixed points for nilpotent actions on the plane and the Cartwright-Littlewood theorem

abstract: The goal of this paper is proving the existence and then localizing global fixed points for nilpotent groups generated by homeomorphisms of the plane satisfying a certain Lipschitz condition. The condition is inspired in a classical result of Bonatti for commuting diffeomorphisms of the 2-sphere and in particular it is satisfied by diffeomorphisms, not necessarily of class $C^{1}$, whose linear part at every point is uniformly close to the identity. In this same setting we prove a version of the Cartwright-Littlewood theorem, obtaining fixed points in any continuum preserved by a nilpotent action.


title: Scale-Free Networks with the Same Degree Distribution: Different Structural Properties

abstract: We have analysed some structural properties of scale-free networks with the same degree distribution. Departing from a degree distribution obtained from the Barab\'asi-Albert (BA) algorithm, networks were generated using four additional different algorithms a (Molloy-Reed, Kalisky, and two new models named A and B) besides the BA algorithm itself. For each network, we have calculated the following structural measures: average degree of the nearest neighbours, central point dominance, clustering coefficient, the Pearson correlation coefficient, and global efficiency. We found that different networks with the same degree distribution may have distinct structural properties. In particular, model B generates decentralized networks with a larger number of components, a smaller giant component size, and a low global efficiency when compared to the other algorithms, especially compared to the centralized BA networks that have all vertices in a single component, with a medium to high global efficiency. The other three models generate networks with intermediate characteristics between B and BA models. A consequence of this finding is that the dynamics of different phenomena on these networks may differ considerably.

categories: cs.SI cond-mat.dis-nn cond-mat.stat-mech math.DS physics.soc-ph


title: Topologie des vari\'et\'es alg\'ebriques r\'eelles de dimension 3

abstract: We review the state of the art of the classification of real uniruled threefolds


title: Sur la mod\'elisation de l'interaction entre polarons et cristaux quantiques

abstract: This text is a summary (in french) of recent works in collaboration with Mathieu Lewin, on the modelling of polarons in quantum crystals.


title: Isotropic realizability of electric fields around critical points

abstract: In this paper we study the isotropic realizability of a given regular gradient field $\nabla u$ as an electric field, namely when $\nabla u$ is solution of the equation $\div\left(\si\nabla u\right)=0$ for some isotropic conductivity $\si>0$. The case of a function $u$ without critical point was investigated in \cite{BMT} thanks to a gradient flow approach. The presence of a critical point needs a specific treatment according to the behavior of the dynamical system around the point. The case of a saddle point is the most favorable and leads us to a characterization of the local isotropic realizability through some boundedness condition involving the laplacian of $u$ along the gradient flow. The case of a sink or a source implies a strong maximum principle under the same boundedness condition. However, when the critical point is not hyperbolic the isotropic realizability is not generally satisfied even piecewisely in the neighborhood of the point. The isotropic realizability in the torus for periodic gradient fields is also discussed in particular when the trajectories of the gradient system are bounded.


title: Guided Random Forest in the RRF Package

abstract: Random Forest (RF) is a powerful supervised learner and has been popularly used in many applications such as bioinformatics. In this work we propose the guided random forest (GRF) for feature selection. Similar to a feature selection method called guided regularized random forest (GRRF), GRF is built using the importance scores from an ordinary RF. However, the trees in GRRF are built sequentially, are highly correlated and do not allow for parallel computing, while the trees in GRF are built independently and can be implemented in parallel. Experiments on 10 high-dimensional gene data sets show that, with a fixed parameter value (without tuning the parameter), RF applied to features selected by GRF outperforms RF applied to all features on 9 data sets and 7 of them have significant differences at the 0.05 level. Therefore, both accuracy and interpretability are significantly improved. GRF selects more features than GRRF, however, leads to better classification accuracy. Note in this work the guided random forest is guided by the importance scores from an ordinary random forest, however, it can also be guided by other methods such as human insights (by specifying $\lambda_i$). GRF can be used in "RRF" v1.4 (and later versions), a package that also includes the regularized random forest methods.


title: Constraints on the general solutions of Einstein cosmological equations by Hubble parameter times cosmic age: a historical perspective

abstract: In a historical perspective, compact solutions of Einstein's equations, including the cosmological constant and the curvature terms, are obtained, starting from two recent observational estimates of the Hubble's parameter (H0) and the "age" of the universe (t0). Cosmological implications for {\Lambda}CDM ({\Lambda} Cold Dark Matter), KOFL (k Open Friedman-Lemaitre), plus two mixed solutions are investigated, under the constraints imposed by the relatively narrow current uncertainties. Quantitative results obtained for the KOFL case seem to be compatible with matter density and the highest observed red-shifts from distant galaxies, while those obtained for the {\Lambda}CDM may be more difficult to reconcile.


title: Deep Learning using Linear Support Vector Machines

abstract: Recently, fully-connected and convolutional neural networks have been trained to achieve state-of-the-art performance on a wide variety of tasks such as speech recognition, image classification, natural language processing, and bioinformatics. For classification tasks, most of these "deep learning" models employ the softmax activation function for prediction and minimize cross-entropy loss. In this paper, we demonstrate a small but consistent advantage of replacing the softmax layer with a linear support vector machine. Learning minimizes a margin-based loss instead of the cross-entropy loss. While there have been various combinations of neural nets and SVMs in prior art, our results using L2-SVMs show that by simply replacing softmax with linear SVMs gives significant gains on popular deep learning datasets MNIST, CIFAR-10, and the ICML 2013 Representation Learning Workshop's face expression recognition challenge.


title: Deformations of period lattices of flexible polyhedral surfaces

abstract: In the end of the 19th century Bricard discovered a phenomenon of flexible polyhedra, that is, polyhedra with rigid faces and hinges at edges that admit non-trivial flexes. One of the most important results in this field is a theorem of Sabitov asserting that the volume of a flexible polyhedron is constant during the flexion. In this paper we study flexible polyhedral surfaces in the 3-space two-periodic with respect to translations by two non-colinear vectors that can vary continuously during the flexion. The main result is that the period lattice of a flexible two-periodic surface homeomorphic to a plane cannot have two degrees of freedom.


title: Characterization of a cryogenic beam source for atoms and molecules

abstract: We present a combined experimental and theoretical study of beam formation from a cryogenic buffer gas cell. Atoms and molecules are loaded into the cell by laser ablation of a target, and are cooled and swept out of the cell by a flow of cold helium. We study the thermalization and flow dynamics inside the cell and measure how the speed, temperature, divergence and extraction efficiency of the beam are influenced by the helium flow. We use a finite element model to simulate the flow dynamics and use the predictions of this model to interpret our experimental results.


title: On lattices, distinct distances, and the Elekes-Sharir framework

abstract: In this note we consider distinct distances determined by points in an integer lattice. We first consider Erdos's lower bound for the square lattice, recast in the setup of the so-called Elekes-Sharir framework \cite{ES11,GK11}, and show that, without a major change, this framework \emph{cannot} lead to Erdos's conjectured lower bound. This shows that the upper bound of Guth and Katz \cite{GK11} for the related 3-dimensional line-intersection problem is tight for this instance. The gap between this bound and the actual bound of Erdos arises from an application of the Cauchy-Schwarz inequality (which is an integral part of the Elekes-Sharir framework). Our analysis relies on two number-theoretic results by Ramanujan. We also consider distinct distances in rectangular lattices of the form $\{(i,j) \mid 0\le i\le n^{1-\alpha},\ 0\le j\le n^{\alpha}\}$, for some $0<\alpha<1/2$, and show that the number of distinct distances in such a lattice is $\Theta(n)$. In a sense, our proof "bypasses" a deep conjecture in number theory, posed by Cilleruelo and Granville \cite{CG07}. A positive resolution of this conjecture would also have implied our bound.


title: Testing Convergence for Global Accretion Disks

abstract: Global disk simulations provide a powerful tool for investigating accretion and the underlying magnetohydrodynamic turbulence driven by the magneto-rotational instability (MRI). Using them to predict accurately quantities such as stress, accretion rate, and surface brightness profile requires that purely numerical effects, arising from both resolution and algorithm, be understood and controlled. We use the flux-conservative Athena code to conduct a series of experiments on disks having a variety of magnetic topologies to determine what constitutes adequate resolution. We develop and apply several resolution metrics: Qz and Qphi, the ratio of the grid zone size to the characteristic MRI wavelength, alpha_mag, the ratio of the Maxwell stress to the magnetic pressure, and the ratio of radial to toroidal magnetic field energy. For the initial conditions considered here, adequate resolution is characterized by Qz > 15, Qphi > 20, alpha_mag = 0.45, and a field energy ratio of 0.2. These values are associated with > 35 zones per scaleheight, a result consistent with shearing box simulations. Numerical algorithm is also important. Use of the HLLE flux solver or second-order interpolation can significantly degrade the effective resolution compared to the HLLD flux solver and third-order interpolation. Resolution at this standard can be achieved only with large numbers of grid zones, arranged in a fashion that matches the symmetries of the problem and the scientific goals of the simulation.


title: Projective geometries in exponentially dense matroids. II

abstract: We show for each positive integer $a$ that, if $\mathcal{M}$ is a minor-closed class of matroids not containing all rank-$(a+1)$ uniform matroids, then there exists an integer $c$ such that either every rank-$r$ matroid in $\mathcal{M}$ can be covered by at most $r^c$ rank-$a$ sets, or $\mathcal{M}$ contains the GF$(q)$-representable matroids for some prime power $q$ and every rank-$r$ matroid in $\mathcal{M}$ can be covered by at most $cq^r$ rank-$a$ sets. In the latter case, this determines the maximum density of matroids in $\mathcal{M}$ up to a constant factor.


title: Interband spin-orbit coupling between anti-parallel spin states in Pb quantum well states

abstract: Using spin and angle-resolved photoemission spectroscopy we investigate a momentum region in Pb quantum well states on Si(111) where hybridization between Rashba-split bands alters the band structure significantly. Starting from the Rashba regime where the dispersion of the quasi-free two-dimensional electron gas is well described by two spin-polarized parabolas, we find a breakdown of the Rashba behavior which manifests itself (i) in a spin splitting that is no longer proportional to the in-plane momentum and (ii) in a reversal of the sign of the momentum splitting. Our experimental findings are well explained by including interband spin-orbit coupling that mixes Rashba-split states with anti-parallel rather than parallel spins. Similar results for Pb/Cu(111) reveal that the proposed hybridization scenario is independent on the supporting substrate.


title: Geometry of the $L_q$-centroid bodies of an isotropic log-concave measure

abstract: We study some geometric properties of the $L_q$-centroid bodies $Z_q(\mu)$ of an isotropic log-concave measure $\mu $ on ${\mathbb R}^n$. For any $2\ls q\ls\sqrt{n}$ and for $\varepsilon \in (\varepsilon_0(q,n),1)$ we determine the inradius of a random $(1-\varepsilon)n$-dimensional projection of $Z_q(\mu)$ up to a constant depending polynomially on $\varepsilon $. Using this fact we obtain estimates for the covering numbers $N(\sqrt{\smash[b]{q}}B_2^n,tZ_q(\mu))$, $t\gr 1$, thus showing that $Z_q(\mu)$ is a $\beta $-regular convex body. As a consequence, we also get an upper bound for $M(Z_q(\mu))$.


title: Differential function spectra, the differential Becker-Gottlieb transfer, and applications to differential algebraic K-theory

abstract: We develop differential algebraic K-theory for rings of integers in number fields and we construct a cycle map from geometrized bundles of modules over such a ring to the differential algebraic K-theory. We also treat some of the foundational aspects of differential cohomology, including differential function spectra and the differential Becker-Gottlieb transfer. We then state a transfer index conjecture about the equality of the Becker-Gottlieb transfer and the analytic transfer defined by Lott. In support of this conjecture, we derive some non-trivial consequences which are provable by independent means.


title: Lower bound for the maximum of some derivative of Hardy's function

abstract: Under the Riemann hypothesis, we use the distribution of zeros of the zeta function to get a lower bound for the maximum of some derivative of Hardy's function.


title: Performance of the Muon Identification at LHCb

abstract: The performance of the muon identification in LHCb is extracted from data using muons and hadrons produced in J/\psi->\mu\mu, \Lambda->p\pi and D^{\star}->\pi D0(K\pi) decays. The muon identification procedure is based on the pattern of hits in the muon chambers. A momentum dependent binary requirement is used to reduce the probability of hadrons to be misidentified as muons to the level of 1%, keeping the muon efficiency in the range of 95-98%. As further refinement, a likelihood is built for the muon and non-muon hypotheses. Adding a requirement on this likelihood that provides a total muon efficiency at the level of 93%, the hadron misidentification rates are below 0.6%.


title: Application of the Space-Time Method to Stimulated Raman Adiabatic Passage on the Simple Harmonic Oscillator

abstract: The space-time method is applied to a model system-the Simple Harmonic Oscillator in a laser field to simulate the Stimulated Raman Adiabatic Passage (STIRAP) process. The Space-Time method is a computational theory first introduced by Weatherford et. al. to solve Time-Dependent Systems with one boundary value and applied to electron spin system with invariant Hamiltonian [Journal of Molecular Structure {\bf 592} 47]. The implementation in the present work provides an efficient and general way to solve the Time-Dependent Schr{\"o}dinger Equation and can be applied to multi-state systems. The algorithm for simulating the Simple Harmonic Oscillator STIRAP can be applied to solve STIRAP problems for complex systems.


title: Innerness of continuous derivations on algebras of measurable operators affiliated with finite von Neumann algebras

abstract: This paper is devoted to derivations on the algebra $S(M)$ of all measurable operators affiliated with a finite von Neumann algebra $M.$ We prove that if $M$ is a finite von Neumann algebra with a faithful normal semi-finite trace $\tau$, equipped with the locally measure topology $t,$ then every $t$-continuous derivation $D:S(M)\rightarrow S(M)$ is inner. A similar result is valid for derivation on the algebra $S(M,\tau)$ of $\tau$-measurable operators equipped with the measure topology $t_{\tau}$.


title: Q-balls and the chiral vortical effect

abstract: The paper considers the possibility of the existence of the axial current, that is directed along the axis of rotation and arises due to the presence of the Q-ball.


title: Game of Life on Phyllosilicates: Gliders, Oscillators and Still Life

abstract: A phyllosilicate is a sheet of silicate tetrahedra bound by basal oxygens. A phyllosilicate automaton is a regular network of finite state machines --- silicon nodes and oxygen nodes --- which mimics structure of the phyllosilicate. A node takes states 0 and 1. Each node updates its state in discrete time depending on a sum of states of its three (silicon) or six (oxygen) neighbours. Phyllosilicate automata exhibit localizations attributed to Conway's Game of Life: gliders, oscillators, still lifes, and a glider gun. Configurations and behaviour of typical localizations, and interactions between the localizations are illustrated.


title: Central Limit Theorems for Classical Likelihood Ratio Tests for High-Dimensional Normal Distributions

abstract: For random samples of size n obtained from p-variate normal distributions, we consider the classical likelihood ratio tests (LRT) for their means and covariance matrices in the high-dimensional setting. These test statistics have been extensively studied in multivariate analysis and their limiting distributions under the null hypothesis were proved to be chi-square distributions as n goes to infinity and p remains fixed. In this paper, we consider the high-dimensional case where both p and n go to infinity with p=n/y in (0, 1]. We prove that the likelihood ratio test statistics under this assumption will converge in distribution to normal distributions with explicit means and variances. We perform the simulation study to show that the likelihood ratio tests using our central limit theorems outperform those using the traditional chi-square approximations for analyzing high-dimensional data.


title: Generalised Particle Filters with Gaussian Mixtures

abstract: Stochastic filtering is defined as the estimation of a partially observed dynamical system. A massive scientific and computational effort is dedicated to the development of numerical methods for approximating the solution of the filtering problem. Approximating the solution of the filtering problem with Gaussian mixtures has been a very popular method since the 1970s (see [1],[2],[46],[49]). Despite nearly fifty years of development, the existing work is based on the success of the numerical implementation and is not theoretically justified. This paper fills this gap and contains a rigorous analysis of a new Gaussian mixture approximation to the solution of the filtering problem. We deduce the L^2-convergence rate for the approximating system and show some numerical example to test the new algorithm.


title: Distributions of Angles in Random Packing on Spheres

abstract: This paper studies the asymptotic behaviors of the pairwise angles among n randomly and uniformly distributed unit vectors in R^p as the number of points n -> infinity, while the dimension p is either fixed or growing with n. For both settings, we derive the limiting empirical distribution of the random angles and the limiting distributions of the extreme angles. The results reveal interesting differences in the two settings and provide a precise characterization of the folklore that "all high-dimensional random vectors are almost always nearly orthogonal to each other". Applications to statistics and machine learning and connections with some open problems in physics and mathematics are also discussed.


title: Spatially distributed social complex networks

abstract: We propose a bare-bones stochastic model that takes into account both the geographical distribution of people within a country and their complex network of connections. The model, which is designed to give rise to a scale-free network of social connections and to visually resemble the geographical spread seen in satellite pictures of the Earth at night, gives rise to a power-law distribution for the ranking of cities by population size (but for the largest cities) and reflects the notion that highly connected individuals tend to live in highly populated areas. It also yields some interesting insights regarding Gibrat's law for the rates of city growth (by population size), in partial support of the findings in a recent analysis of real data [Rozenfeld et al., Proc. Natl. Acad. Sci. U.S.A. 105, 18702 (2008)]. The model produces a nontrivial relation between city population and city population density and a superlinear relationship between social connectivity and city population, both of which seem quite in line with real data.

categories: physics.soc-ph cond-mat.stat-mech cs.SI nlin.AO


title: Slime mould tactile sensor

abstract: Slime mould P. polycephalum is a single cells visible by unaided eye. The cells shows a wide spectrum of intelligent behaviour. By interpreting the behaviour in terms of computation one can make a slime mould based computing device. The Physarum computers are capable to solve a range of tasks of computational geometry, optimisation and logic. Physarum computers designed so far lack of localised inputs. Commonly used inputs --- illumination and chemo-attractants and -repellents --- usually act on extended domains of the slime mould's body. Aiming to design massive-parallel tactile inputs for slime mould computers we analyse a temporal dynamic of P. polycephalum's electrical response to tactile stimulation. In experimental laboratory studies we discover how the Physarum responds to application and removal of a local mechanical pressure with electrical potential impulses and changes in its electrical potential oscillation patterns.


title: On The Coordinated g-convex Dominated Functions

abstract: In this study, we define g-convex dominated functions on the co-ordinates and prove some Hadamard-type, Fejer-type inequalities for this new class of functions. We also give some results related to the functional H.


title: A Distributed Algorithm for Solving Positive Definite Linear Equations over Networks with Membership Dynamics

abstract: This paper considers the problem of solving a symmetric positive definite system of linear equations over a network of agents with arbitrary asynchronous interactions and membership dynamics. The latter implies that each agent is allowed to join and leave the network at any time, for infinitely many times, and lose all its memory upon leaving. We develop Subset Equalizing (SE), a distributed asynchronous algorithm for solving such a problem. To design and analyze SE, we introduce a novel time-varying Lyapunov-like function, defined on a state space with changing dimension, and a generalized concept of network connectivity, capable of handling such interactions and membership dynamics. Based on them, we establish the boundedness, asymptotic convergence, and exponential convergence of SE, along with a bound on its convergence rate. Finally, through extensive simulation, we show that SE is effective in a volatile agent network and that a special case of SE, termed Groupwise Equalizing, is significantly more bandwidth/energy efficient than two existing algorithms in multi-hop wireless networks.


title: Propagators in two-dimensional lattices

abstract: This paper is devoted to the computation of discrete propagators in two-dimensional crystals and their application to a number of time dependent problems. The methods to compute such kernels are provided by a tight-binding representation of Dirac matrices and the generalizations of Bessel functions. Diffusive effects of point-like distributions on crystalline sheets are studied in a second quantization scheme. In the last part, a compendium of propagators is presented. The cases of square, triangular and hexagonal arrays are covered.


title: Gauge invariant electromagnetic properties of fermions induced by CPT violation in the Standard Model Extension

abstract: Low-energy Lorentz-invariant quantities could receive contributions from a fundamental theory producing small Lorentz-violating effects. Within the Lorentz-violating extension of quantum electrodynamics, we investigate, perturbatively, the contributions to the one-loop $ff\gamma$ vertex from the $CPT$-violating axial coupling of a vector background field to fermions. We find that the resulting vertex function has a larger set of Lorentz structures than the one characterizing the usual, Lorentz invariant, parametrization of the $ff\gamma$ vertex. We prove gauge invariance of the resulting one-loop expression through a set of gauge invariant nonrenormalizable operators introducing new-physics effects at the first and second orders in Lorentz violation, and which generate tree-level contributions to the $ff\gamma$ vertex. Whereas loop contributions involving parameters that violate Lorentz invariance at the first order are $CPT$-odd, those arising at the second order are $CPT$-even, so that contributions to low-energy physics are restricted to emerge for the first time at the second order. In this context, we derive a contribution to anomalous magnetic moment of fermions, which we use to set a bound on Lorentz violation.


title: Interacting Particle Systems in Time-Dependent Geometries

abstract: Many complex structures and stochastic patterns emerge from simple kinetic rules and local interactions, and are governed by scale invariance properties in combination with effects of the global geometry. We consider systems that can be described effectively by space-time trajectories of interacting particles, such as domain boundaries in two-dimensional growth or river networks. We study trajectories embedded in time-dependent geometries, and the main focus is on uniformly expanding or decreasing domains for which we obtain an exact mapping to simple fixed domain systems while preserving the local scale invariance properties. This approach was recently introduced in [A. Ali et al., Phys. Rev. E 87, 020102(R) (2013)] and here we provide a detailed discussion on its applicability for self-affince Markovian models, and how it can be adapted to self-affine models with memory or explicit time dependence. The mapping corresponds to a non-linear time transformation which convergences to a finite value for a large class of trajectories, enabling an exact analysis of asymptotic properties in expanding domains. We further provide a detailed discussion of different particle interactions and generalized geometries. All our findings are based on exact computations and are illustrated numerically for various examples, including L\'evy processes and fractional Brownian motion.


title: Epidemic-like Proximity-based Traffic Offloading

abstract: Cellular networks are overloaded due to the mobile traffic surge, and mobile social network (MSNets) carrying information flow can help reduce cellular traffic load. If geographically-nearby users directly adopt WiFi or Bluetooth technology (i.e., leveraging proximity-based communication) for information spreading in MSNets, a portion of mobile traffic can be offloaded from cellular networks. For many delay-tolerant applications, it is beneficial for traffic offloading to pick some seed users as information sources, which help further spread the information to others in an epidemic-like manner using proximity-based communication. In this paper, we develop a theoretical framework to study the issue of choosing only k seed users so as to maximize the mobile traffic offloaded from cellular networks via proximity-based communication. We introduce a gossip-style social cascade (GSC) model to model the information diffusion process, which captures the epidemic-like nature of proximity-based communication and characterizes users' social participation as well. For static networks as a special-case study and mobile networks, we establish an equivalent view and a temporal mapping of the information diffusion process, respectively, leveraging virtual coupon collectors. We further prove the submodularity in the information diffusion and propose a greedy algorithm to choose the seed users for proximity-based traffic offloading, yielding a solution within about 63% of the optimal value to the traffic offloading maximization (TOM) problem. Experiments are carried out to study the offloading performance of our approach, illustrating that proximity-based communication can offload cellular traffic by over 60% with a small number of seed users and the greedy algorithm significantly outperforms the heuristic and random algorithms.


title: Multi-frequency subspace migration for imaging of perfectly conducting, arc-like cracks

abstract: Multi-frequency subspace migration imaging technique are usually adopted for the non-iterative imaging of unknown electromagnetic targets such as cracks in the concrete walls or bridges, anti-personnel mines in the ground, etc. in the inverse scattering problems. It is confirmed that this technique is very fast, effective, robust, and can be applied not only full- but also limited-view inverse problems if suitable number of incident and corresponding scattered field are applied and collected. But in many works, the application of such technique is somehow heuristic. Under the motivation of such heuristic application, this contribution analyzes the structure of imaging functional employed in the subspace migration imaging technique in two-dimensional inverse scattering when the unknown target is arbitrary shaped, arc-like perfectly conducting cracks located in the homogeneous two-dimensional space. Opposite to the Statistical approach based on the Statistical Hypothesis Testing, our approach is based on the fact that subspace migration imaging functional can be expressed by a linear combination of Bessel functions of integer order of the first kind. This is based on the structure of the Multi-Static Response (MSR) matrix collected in the far-field at nonzero frequency in either Transverse Magnetic (TM) mode or Transverse Electric (TE) mode. Explored expression of imaging functionals gives us certain properties of subspace migration and an answer of why multi-frequency enhances imaging resolution. Particularly, we carefully analyze the subspace migration and confirm some properties of imaging when a small number of incident field is applied. Consequently, we simply introduce a weighted multi-frequency imaging functional and confirm that which is an improved version of subspace migration in TM mode.


title: Predicted alternative structure for tantalum metal under high pressure and high temperature

abstract: First-principles simulations have been performed to investigate the phase stability of tantalum metal under high pressure and high temperature (HPHT). We searched its low-energy structures globally using our developed multi-algorithm collaborative (MAC) crystal structure prediction technique. The body-centred cubic (bcc) was found to be stable at pressure up to 300 GPa. The previously reported $\omega$ and \textit{A}15 structures were also reproduced successfully. More interestingly, we observed another phase (space group: \textit{Pnma}, 62) that is more stable than $\omega$ and \textit{A}15. Its stability is confirmed by its phonon spectra and elastic constants. For $\omega$-Ta, the calculated elastic constants and high-temperature phonon spectra both imply that it is neither mechanically nor dynamically stable. Thus, $\omega$ is not the structure to which bcc-Ta transits before melting. On the contrary, the good agreement of \textit{Pnma}-Ta shear sound velocities with experiment suggests \textit{Pnma} is the new structure of Ta implied by the discontinuation of shear sound velocities in recent shock experiment [J. Appl. Phys. \textbf{111}, 033511 (2012)].


title: Locality statistics for anomaly detection in time series of graphs

abstract: The ability to detect change-points in a dynamic network or a time series of graphs is an increasingly important task in many applications of the emerging discipline of graph signal processing. This paper formulates change-point detection as a hypothesis testing problem in terms of a generative latent position model, focusing on the special case of the Stochastic Block Model time series. We analyze two classes of scan statistics, based on distinct underlying locality statistics presented in the literature. Our main contribution is the derivation of the limiting distributions and power characteristics of the competing scan statistics. Performance is compared theoretically, on synthetic data, and on the Enron email corpus. We demonstrate that both statistics are admissible in one simple setting, while one of the statistics is inadmissible a second setting.


title: Intermediate-Age Stellar Populations in Classical QSO Host Galaxies

abstract: Although mergers and starbursts are often invoked in the discussion of QSO activity in the context of galaxy evolution, several studies have questioned their importance or even their presence in QSO host galaxies. Accordingly, we are conducting a study of z~0.2 QSO host galaxies previously classified as passively evolving elliptical galaxies. We present deep Keck LRIS spectroscopy of a sample of 15 hosts and model their stellar absorption spectra using stellar synthesis models. The high S/N of our spectra allow us to break various degeneracies that arise from different combinations of models, varying metallicities, and contamination from QSO light. We find that none of the host spectra can be modeled by purely old stellar populations and that the majority of the hosts (14/15) have a substantial contribution from intermediate-age populations with ages ranging from 0.7 to 2.4 Gyr. An average host spectrum is strikingly well fit by a combination of an old population and a 2.1 (+0.5, -0.7) Gyr population. The morphologies of the host galaxies suggest that these aging starbursts were induced during the early stages of the mergers that resulted in the elliptical-shaped galaxies that we observe. The current AGN activity likely corresponds to the late episodes of accretion predicted by numerical simulations, which occur near the end of the mergers, whereas earlier episodes may be more difficult to observe due to obscuration. Our off-axis observations prevent us from detecting any current star formation or young stellar populations that may be present in the central few kiloparsecs.


title: Experimental measurements of stress redistribution in flowing emulsions

abstract: We study how local rearrangements alter droplet stresses within flowing dense quasi-two-dimensional emulsions at area fractions $\phi \geq 0.88$. Using microscopy, we measure droplet positions while simultaneously using their deformed shape to measure droplet stresses. We find that rearrangements alter nearby stresses in a quadrupolar pattern: stresses on neighboring droplets tend to either decrease or increase depending on location. The stress redistribution is more anisotropic with increasing $\phi$. The spatial character of the stress redistribution influences where subsequent rearrangements occur. Our results provide direct quantitative support for rheological theories of dense amorphous materials that connect local rearrangements to changes in nearby stress.


title: Study of spatial resolution in a single GEM simulated by Monte-Carlo method

abstract: Spatial resolution is a significant factor in the GEM performance in view of X-rays radiography and UV, visible light imaging. Monte-Carlo method is used to investigate the spatial resolution determined by the transverse diffusion in the device. The simulation results indicate that the electrical parameters, such as the GEM voltages and the electric field at the drift and induction regions, only have minor effects on the spatial resolution. The geometrical parameters and the working gases chosen, on the other hand, are the main parameters that determine the spatial resolution. The spatial resolution is determined more on the drift and diffusion processes than on the avalanche process. Especially for the different working gases, the square root function of the ratio of the electron diffusion coefficient and the mobility has a significant effect on the spatial resolution.


title: KERT: Automatic Extraction and Ranking of Topical Keyphrases from Content-Representative Document Titles

abstract: We introduce KERT (Keyphrase Extraction and Ranking by Topic), a framework for topical keyphrase generation and ranking. By shifting from the unigram-centric traditional methods of unsupervised keyphrase extraction to a phrase-centric approach, we are able to directly compare and rank phrases of different lengths. We construct a topical keyphrase ranking function which implements the four criteria that represent high quality topical keyphrases (coverage, purity, phraseness, and completeness). The effectiveness of our approach is demonstrated on two collections of content-representative titles in the domains of Computer Science and Physics.


title: One can hear the area and curvature of boundary of a domain by hearing the Steklov eigenvalues

abstract: For a given bounded domain $\Omega$ with smooth boundary in a smooth Riemannian manifold $(\mathcal{M},g)$, we show that the Poisson type upper-estimate of the heat kernel associated to the Dirichlet-to-Neumann operator, the Sobolev trace inequality, the Log-Sobolev trace inequality, the Nash trace inequality, and the Rozenblum-Lieb-Cwikel type inequality are all equivalent. Upon decomposing the Dirichlet-to-Neumann operator into a sum of the square root of the Laplacian and a pseudodifferntial operator and by applying Grubb's method of symbolic calculus for the corresponding pseudodifferential heat kernel operators, we establish a procedure to calculate all the coefficients of the asymptotic expansion of the trace of the heat kernel associated to Dirichlet-to-Neumann operator as $t\to 0^+$. In particular, we explicitly give the first four coefficients of this asymptotic expansion. These coefficients give precise information regarding the area and curvatures of the boundary of the domain in terms of the spectrum of the Steklov problem.


title: The spectral variability of quasar SDSS J030639.57+000343.1

abstract: We compiled a sample of 60 quasars with spectroscopy on at least six epochs from the Sloan Digital Sky Survey (SDSS) to study the variabilities of the spectral shape, the continuum and the emission lines luminosity. In this paper, we present the results of SDSS J030639.57+000343.1. We found a strong anti-correlation between the continuum luminosity at $5100 \AA$ and the spectral index, implying a bluer-when-brighter trend. The luminosity of the broad $\rm H_\alpha$ line is proportion to the continuum luminosity at $5100 \AA$. Correspondingly, we did not find strong correlation between the equivalent width of broad $\rm H\alpha$ and the continuum luminosity, i.e. no baldwin effect of broad $\rm H\alpha$ in this source.


title: A long-term VLBA monitoring campaign of the v=1, J=1-0 SiO masers toward TX Cam - I. Morphology and Shock Waves

abstract: We present the latest and final version of the movie of the SiO masers toward the Mira variable TX Cam. The new version consists of 112 frames (78 successfully reduced epochs) with data covering almost three complete stellar cycles between 24 May 1997 to 25 January 2002, observed with the VLBA. In this paper we examine the global morphology, kinematics and variability of the masering zone. The morphology of the emission is confined in a structure that usually resembles a ring or an ellipse, with occasional deviations due to localised phenomena. The ring appears to be contracting and expanding, although for the first cycle contraction is not observed. The width and outer boundary of the masering zone follow the stellar pulsation. Our data seem to be consistent with a shock being created once per stellar cycle at maximum that propagates with a velocity of ~7 km/s. The difference in velocities along different axes strongly suggests that the outflow in TX Cam is bipolar. The contribution of projection is examined and our results are compared with the latest theoretical model.


title: GASPS - a Herschel survey of gas and dust in Protoplanetary Disks: Summary and Initial Statistics

abstract: GASPS is a far-infrared line and continuum survey of protoplanetary and young debris disks using PACS on the Herschel Space Observatory. The survey includes [OI] at 63 microns, as well as 70, 100 and 160um continuum, with the brightest objects also studied in [OI]145um, [CII]157um, H2O and CO. Targets included T Tauri stars and debris disks in 7 nearby young associations, and a sample of isolated Herbig AeBe stars. The aim was to study the global gas and dust content in a wide disk sample, systemically comparing the results with models. In this paper we review the main aims, target selection and observing strategy. We show initial results, including line identifications, sources detected, and a first statistical study. [OI]63um was the brightest line in most objects, by a factor of ~10. Detection rates were 49%, including 100% of HAeBe stars and 43% of T Tauri stars. Comparison with published dust masses show a dust threshold for [OI]63um detection of ~1e-5 M_solar. Normalising to 140pc distance, 32% with mass 1e-6 - 1e-5 M_solar, and a small number with lower mass were also detected. This is consistent with moderate UV excess and disk flaring. In most cases, continuum and line emission is spatially and spectrally unresolved, suggesting disk emission. ~10 objects were resolved, likely from outflows. Detection rates in [OI]145um, [CII]157um and CO J=18-17 were 20-40%, but [CII] was not correlated with disk mass, suggesting it arises instead from a compact envelope. [OI] detection rates in T Tauri associations of ages 0.3-4Myr were ~50%. ~2 stars were detectable in associations of 5-20Myr, with no detections in associations of age >20Myr. Comparing with the total number of young stars, and assuming a ISM-like gas/dust ratio, this indicates that ~18% of stars retain a gas-rich disk of total mass >1M_Jupiter for 1-4Myr, 1-7% keep such disks for 5-10Myr, and none remain beyond 10-20Myr.


title: Super sub-wavelength patterns in photon coincidence detection

abstract: High-precision measurements implemented by means of light is desired in all fields of science. However, light is a wave and Rayleigh criterion gives us a diffraction limitation in classical optics which restricts to get arbitrary high resolution. Sub-wavelength interference has a potential application in lithography to beat the classical Rayleigh limit of resolution. We carefully study the second-order correlation theory to get the physics behind sub-wavelength interference in photon coincidence detection. A Young's double-slit experiment with pseudo-thermal light is carried out to test the second-order correlation pattern. The result shows that when different scanning ways of two point detectors are chosen, one can get super sub-wavelength interference patterns. We then give a theoretical explanation to this surprising result, and find this explanation is also suitable for the result by using entangled light. Furthermore, we discuss the limitation of this kind of super sub-wavelength interference patterns in quantum lithography.


title: Hardening transition in a one-dimensional model for ferrogels

abstract: We introduce and investigate a coarse-grained model for quasi one-dimensional ferrogels. In our description the magnetic particles are represented by hard spheres with a magnetic dipole moment in their centers. Harmonic springs connecting these spheres mimic the presence of a cross-linked polymer matrix. A special emphasis is put on the coupling of the dipolar orientations to the elastic deformations of the matrix, where a memory effect of the orientations is included. Although the particles are displaced along one spatial direction only, the system already shows rich behavior: as a function of the magnetic dipole moment, we find a phase transition between "soft-elastic" states with finite interparticle separation and finite compressive elastic modulus on the one hand, and "hardened" states with touching particles and therefore diverging compressive elastic modulus on the other hand. Corresponding phase diagrams are derived neglecting thermal fluctuations of the magnetic particles. In addition, we consider a situation in which a spatially homogeneous magnetization is initially imprinted into the material. Depending on the strength of the magneto-mechanical coupling between the dipole orientations and the elastic deformations, the system then relaxes to a uniaxially ferromagnetic, an antiferromagnetic, or a spiral state of magnetization to minimize its energy. One purpose of our work is to provide a largely analytically solvable approach that can provide a benchmark to test future descriptions of higher complexity. From an applied point of view, our results could be exploited, for example, for the construction of novel damping devices of tunable shock absorbance.


title: Aboriginal Oral Traditions of Australian Impact Craters

abstract: We explore Aboriginal oral traditions that relate to Australian meteorite craters. Using the literature, first-hand ethnographic records, and fieldtrip data, we identify oral traditions and artworks associated with four impact sites: Gosses Bluff, Henbury, Liverpool, and Wolfe Creek. Oral traditions describe impact origins for Gosses Bluff, Henbury, and Wolfe Creek craters and non-impact origins of Liverpool crater, with Wolfe Creek and Henbury having both impact and non-impact origins in oral tradition. Three impact sites that are believed to have formed during human habitation of Australia - Dalgaranga, Veevers, and Boxhole - do not have associated oral traditions that are reported in the literature.


title: Pair production of a 125 GeV Higgs boson in MSSM and NMSSM at the ILC

abstract: In this work we investigate the Higgs pair production in the MSSM and NMSSM at the photon-photon collision of the ILC. We consider various experimental constraints and scan over the parameter space of the MSSM and NMSSM. Then we calculate the cross section of Higgs pair production in the allowed parameter space and compare the results with the predictions in the SM. We find that the large enhancement of the cross section in the MSSM is mainly due to the contributions from the loops mediated by the stau, while in the NMSSM it is mainly due to the contributions from the top-squark loops. For light $m_{\tilde{\tau}_1}$ and large $\mu\tan\beta$, the production rate can be enhanced by a factor of 18 in the MSSM (relative to the SM prediction). And for a large trilinear soft breaking parameter $A_t$ and a moderately light top-squark, it can also be enhanced by a factor of 2 in the NMSSM. Moreover, we also calculate the $\chi^2$ values with the LHC Higgs data and display the results for the parameter space with $\chi^2$ better than the SM value.


title: On sequences without geometric progressions

abstract: An improved upper bound is obtained for the density of sequences of positive integers that contain no k-term geometric progression.


title: Classical Hardness of Learning with Errors

abstract: We show that the Learning with Errors (LWE) problem is classically at least as hard as standard worst-case lattice problems, even with polynomial modulus. Previously this was only known under quantum reductions. Our techniques capture the tradeoff between the dimension and the modulus of LWE instances, leading to a much better understanding of the landscape of the problem. The proof is inspired by techniques from several recent cryptographic constructions, most notably fully homomorphic encryption schemes.


title: An efficient method for evaluating BEM singular integrals on curved elements with application in acoustic analysis

abstract: The polar coordinate transformation (PCT) method has been extensively used to treat various singular integrals in the boundary element method (BEM). However, the resultant integrands of the PCT tend to become nearly singular when (1) the aspect ratio of the element is large or (2) the field point is closed to the element boundary; thus a large number of quadrature points are needed to achieve a relatively high accuracy. In this paper, the first problem is circumvented by using a conformal transformation so that the geometry of the curved physical element is preserved in the transformed domain. The second problem is alleviated by using a sigmoidal transformation, which makes the quadrature points more concentrated around the near singularity. By combining the proposed two transformations with the Guiggiani's method in [M. Guiggiani, \emph{et al}. A general algorithm for the numerical solution of hypersingular boundary integral equations. \emph{ASME Journal of Applied Mechanics}, 59(1992), 604-614], one obtains an efficient and robust numerical method for computing the weakly-, strongly- and hyper-singular integrals in high-order BEM with curved elements. Numerical integration results show that, compared with the original PCT, the present method can reduce the number of quadrature points considerably, for given accuracy. For further verification, the method is incorporated into a 2-order Nystr\"om BEM code for solving acoustic Burton-Miller boundary integral equation. It is shown that the method can retain the convergence rate of the BEM with much less quadrature points than the existing PCT. The method is implemented in C language and freely available.


title: Belief propagation with multipoint correlations and its application to Inverse problem

abstract: We give explicit formulas of the Bethe approximation with multipoint correlations for systems with magnetic field. The obtained formulas include the closed form of the magnetization and the correlations between adjacent degrees of freedom. On the basis of our results, we propose a new iterative algorithm of the improved Bethe approximation. We confirm that the proposed technique is available for the random spin systems and indeed gives more accurate locations of the critical points. We discuss the possibility of the application of our method to the Inverse Ising model by use of these equations.


title: Some Topics in Quantum Games

abstract: This work, based on the author's MA thesis, concentrates on simultaneous move quantum games of two players. A numerical algorithm based on the method of best response functions, designed to search for pure strategy Nash equilibrium in quantum games is constructed. It is then used show that there is a critical degree of entanglement above which a pure strategy Nash Equilibrium does not exist. The behavior of the two payoffs approach the cooperative classical ones. Bayesian quantum games are defined, and it is shown that under certain conditions, there {\it is} a pure strategy Nash equilibrium in such games even when entanglement is maximal. The basic ingredients of a quantum game based on a two-players {\it three} decisions classical games are presented. This requires the definition of trits (instead of bits) and quantum trits (instead of quantum bits). It is shown that in this quantum game, there is no classical commensurability in the sense that the classical strategies are not obtained as a special case of the quantum strategies. An entanglement operator transforming a pure two qutrit state into a maximally entangled one is explicitly constructed.


title: Calibrating the fundamental plane with SDSS DR8 data

abstract: We present a calibration of the fundamental plane using SDSS Data Release 8. We analysed about 93000 elliptical galaxies up to $z<0.2$, the largest sample used for the calibration of the fundamental plane so far. We incorporated up-to-date K-corrections and used GalaxyZoo data to classify the galaxies in our sample. We derived independent fundamental plane fits in all five Sloan filters u, g, r, i and z. A direct fit using a volume-weighted least-squares method was applied to obtain the coefficients of the fundamental plane, which implicitly corrects for the Malmquist bias. We achieved an accuracy of 15% for the fundamental plane as a distance indicator. We provide a detailed discussion on the calibrations and their influence on the resulting fits. These re-calibrated fundamental plane relations form a well-suited anchor for large-scale peculiar-velocity studies in the nearby universe. In addition to the fundamental plane, we discuss the redshift distribution of the elliptical galaxies and their global parameters.


title: Inelastic transport dynamics through attractive impurity in charge Kondo regime

abstract: Within the frame of quantum dissipation theory, we develop a new hierarchical equations of motion theory, combined with the small polaron transformation. We fully investigate the electron transport of a single attractive impurity system with the strong electron-phonon coupling in charge- Kondo regime. Numerical results demonstrate the following facts: (i) The density of states curve shows that the attraction mechanism results in not only the charge-Kondo resonance (i.e. elastic pair-transition resonance), but also the inelastic pair-transition resonances and inelastic cotunneling resonances. These signals are separated discernibly in asymmetric levels about Fermi level; (ii) The differential conductance spectrum shows the distinct peaks or wiggles and the abnormal split of pair-transition peaks under the asymmetric bias; (iii) The improvement of bath-temperature can enhance the phonon-emission (absorption)-assisted sequential tunneling and also strengthen the signals of inelastic pair-transition under the asymmetric bias; (iv) The inelastic dynamics driven by the ramp-up time-dependent voltage presents clear steps which can be tailored by the duration-time and bath-temperature; (v) The linear-response spectrum obtained from the linear-response theory in the Liouville space reveals the excitation signals of electrons' dynamical transition. Briefly, the physics caused by attraction mechanism stems from the double-occupancy or vacant-occupation of impurity system


title: On the general homogenization of von Karman plate equations from 3D nonlinear elasticity

abstract: Starting from 3D elasticity equations we derive the model of the homogenized von K\'arm\'an plate by means of $\Gamma$-convergence. This generalizes the recent results, where the material oscillations were assumed to be periodic.


title: Nitrogen and argon doping of niobium for superconducting radio frequency cavities: a pathway to highly efficient accelerating structures

abstract: We report a surface treatment that systematically improves the quality factor of niobium radio frequency cavities beyond the expected limit for niobium. A combination of annealing in a partial pressure of nitrogen or argon gas and subsequent electropolishing of the niobium cavity surface leads to unprecedented low values of the microwave surface resistance, and an improvement in the efficiency of the accelerating structures up to a factor of 3, reducing the cryogenic load of superconducting cavities for both pulsed and continuous duty cycles. The field dependence of the surface resistance is reversed compared to standardly treated niobium.


title: Analytical Formulae for the Loss Factors and Wakefields of a Rectangular Accelerating Structure

abstract: Wakefields in a rectangular accelerating structure can be calculated in time domain by directly solving Maxwell's equations by a 3D code. In this paper, we will give analytical formulae to calculate the synchronous modes' loss factors. From these analytical formulae on can get the delta function wakefields. The relations between the loss factors (wakefields) and the structure geometrical dimensions are well established. These analytical expressions of loss factors can be used also in a single rectangular resonant cavity. It is shown that the potential application of a rectangular accelerating structure is to accelerate a flat beam in a linear collider.


title: Probability density function of the Cartesian x-coordinate of the random point inside the hypersphere

abstract: Consider randomly picked points inside the n-dimensional unit hypersphere centered at the origin of the Cartesian coordinate system. The Cartesian coordinates of the points are random variables, which form an n-dimensional vector for each point. Observing only the x-coordinate I obtained its probability density function (PDF). I show that it is related to the Gaussian distribution: in limit its companion PDF?? converges to the PDF of the standard normal distribution.


title: Revisiting Circular-Based Random Node Simulation

abstract: In literature, a stochastic model for spreading nodes in a cellular cell is available. Despite its existence, the current method does not offer any versatility in dealing with sectored layers. Of course, this needed adaptability could be created synthetically through heuristic means. However, due to selective sampling, such practice dissolves the true randomness sought. Hence, in this paper, a universal exact scattering model is derived. Also, as an alternative to exhaustive simulation, a generic close-form path-loss predictor between a node and a BS is obtained. Further, using these results, an algorithm based on the superposition principle is proposed. This will ensure greater emulation flexibility, and attain a heterogeneous spatial density.


title: Electrical switching dynamics and broadband microwave characteristics of VO2 RF devices

abstract: Vanadium dioxide is a correlated electron system that features a metal-insulator phase transition (MIT) above room temperature and is of interest in high speed switching devices. Here, we integrate VO2 into two-terminal coplanar waveguides and demonstrate a large resistance modulation of the same magnitude (>10^3) in both electrically (i.e. by bias voltage, referred to as E-MIT) and thermally (T-MIT) driven transitions. We examine transient switching characteristics of the E-MIT and observe two distinguishable time scales for switching. We find an abrupt jump in conductivity with a rise time of the order of 10 ns followed by an oscillatory damping to steady state on the order of several {\mu}s. We characterize the RF power response in the On state and find that high RF input power drives VO2 further into the metallic phase, indicating that electromagnetic radiation-switching of the phase transition may be possible. We measure S-parameter RF properties up to 13.5 GHz. Insertion loss is markedly flat at 2.95 dB across the frequency range in the On state and sufficient isolation of over 25 dB is observed in the Off state. We are able to simulate the RF response accurately using both lumped element and 3D electromagnetic models. Extrapolation of our results suggests that optimizing device geometry can reduce insertion loss further and maintain broadband flatness up to 40 GHz.


title: A realization for a $\mathbb{Q}$-Hermitian variation of Hodge structure of Calabi-Yau type with real multiplication

abstract: We show that the $\mathbb{Q}$-descents of the canonical $\mathbb{R}$-variation of Hodge structure of Calabi-Yau type over a tube domain of type $A$ can be realized as sub-variations of Hodge structure of certain $\mathbb{Q}$-variations of Hodge structure which are naturally associated to abelian varieties of (generalized) Weil type.


title: Chip-firing game and partial Tutte polynomial for Eulerian digraphs

abstract: The Chip-firing game is a discrete dynamical system played on a graph, in which chips move along edges according to a simple local rule. Properties of the underlying graph are of course useful to the understanding of the game, but since a conjecture of Biggs that was proved by Merino L\'opez, we also know that the study of the Chip-firing game can give insights on the graph. In particular, a strong relation between the partial Tutte polynomial $T_G(1,y)$ and the set of recurrent configurations of a Chip-firing game (with a distinguished sink vertex) has been established for undirected graphs. A direct consequence is that the generating function of the set of recurrent configurations is independent of the choice of the sink for the game, as it characterizes the underlying graph itself. In this paper we prove that this property also holds for Eulerian directed graphs (digraphs), a class on the way from undirected graphs to general digraphs. It turns out from this property that the generating function of the set of recurrent configurations of an Eulerian digraph is a natural and convincing candidate for generalizing the partial Tutte polynomial $T_G(1,y)$ to this class. Our work also gives some promising directions of looking for a generalization of the Tutte polynomial to general digraphs.


title: Phase structure of hot dense QCD by a histogram method

abstract: We study the phase structure of QCD at high temperature and density by lattice QCD simulations adopting a histogram method. The quark mass dependence and the chemical potential dependence of the nature of phase transition are investigated focusing on the probability distribution function (histogram). The shape of the distribution function changes with the quark mass and chemical potential. Through the shape of the distribution, the critical surface which separates the first order transition and crossover regions in the heavy quark region is determined for the (2+1)-flavor case. Moreover, we determined the critical point at finite density for two-flavor QCD with an intermediate quark mass, using a Gaussian approximation of the complex phase distribution of the quark determinant. The chemical potential dependence of the critical quark mass is also evaluated in the situation where two light quarks and many massive quarks exist. We find that the first order transition region becomes wider with the chemical potential in the many-flavor QCD.


title: Negative K-groups of abelian categories

abstract: We prove that negative K-groups of small abelian categories are trivial.


title: Generation of Multiple Dirac Cones in Graphene under Double-periodic and Quasiperiodic Potentials

abstract: We investigate generation of new Dirac cones in graphene under double-periodic and quasiperiodic superlattice potentials. We first show that double-periodic potentials generate the Dirac cones sporadically, following the Diophantine equation, in spite of the fact that double-periodic potentials are also periodic ones, for which previous studies predict consecutive appearance of the cones. The sporadic appearance is due to the fact that the dispersion relation of graphene is linear only up to an energy cutoff. We then show that quasiperiodic potentials generate the new Dirac cones densely with its density depending on the energy. We also extend the above predictions to other materials of Dirac electrons with different energy cutoffs of the linear dispersion.


title: Higher-order discrete variational problems with constraints

abstract: An interesting family of geometric integrators for Lagrangian systems can be defined using discretizations of the Hamilton's principle of critical action. This family of geometric integrators is called variational integrators. In this paper, we derive new variational integrators for higher-order lagrangian mechanical system subjected to higher-order constraints. From the discretization of the variational principles, we show that our methods are automatically symplectic and, in consequence, with a very good energy behavior. Additionally, the symmetries of the discrete Lagrangian imply that momenta is conserved by the integrator. Moreover, we extend our construction to variational integrators where the lagrangian is explicitly time-dependent. Finally, some motivating applications of higher-order problems are considered; in particular, optimal control problems for explicitly time-dependent underactuated systems and an interpolation problem on Riemannian manifolds.


title: Non-abelian Hodge theory for algebraic curves in characteristic p

abstract: Let G be a reductive group over an algebraically closed field of positive characteristic. Let C be a smooth projective curve over k. We give a description of the moduli space of flat G-bundles in terms of the moduli space of G-Higgs bundles over the Frobenius twist C' of C. This description can be regarded as the non-abelian Hodge theory for curves in positive characteristic.


title: Emergence of Reflectionless Scattering from Linearizations of Integrable PDEs around Solitons

abstract: We present four examples of integrable partial differential equations (PDEs) of mathematical physics that---when linearized around a stationary soliton---exhibit scattering without reflection at {\it all} energies. Starting from the most well-known and the most empirically relevant phenomenon of the transparency of one-dimensional bright bosonic solitons to Bogoliubov excitations, we proceed to the sine-Gordon, Korteweg-de Vries, and Liouville's equation whose stationary solitons also support our assertion. The proposed connection between integrability and reflectionless scattering seems to span at least two distinct paradigms of integrability: S-integrability in the first three cases, and C-integrability in the last one. We argue that the transparency of linearized integrable PDEs is necessary to ensure that they can support the transparency of stationary solitons in the original integrable PDEs. As contrasting cases, the analysis is further extended to cover two non-integrable systems: a sawtooth-Gordon and a $\phi^4$ model.

categories: cond-mat.quant-gas math-ph math.MP nlin.PS nlin.SI


title: Effect of substrate temperature on the spin transport property in C60-based spin valve devices

abstract: We report the effect of the substrate temperature on the magnetoresistance (MR) of the C60-based spin valve (SV) devices with the sandwich configuration of La0.67Sr0.33MnO3 (LSMO)/C60/cobalt (Co). The C60 interlayer deposited at different substrate temperatures resulted in four types of devices. We observed all types of devices showed a monotonic increase in MR ratio with the substrate temperature. And more interesting, an especially large MR (28.5%) were obtained in the device fabricated at higher substrate temperature, while for the other types of devices, the MR magnitudes were just about a few percent. Based on the I-V measurements, as well as SEM and AFM characteristics, we have obtained that the higher substrate temperature can cause many pits and hollows in the organic film, and these pits will increase the tunneling probability of spin-polarized carriers from one ferromagnetic electrode to the other.


title: The nucleation parameter for heavy-ion induced bubble nucleation in superheated emulsion detector

abstract: The values of the nucleation parameter, k, for bubble nucleation induced by high energy heavy ions 12C (180 MeV/u), 20Ne (400 MeV/u) and 28Si (350 MeV/u) in superheated emulsion detector are determined by comparing the experimentally obtained normalized count rates with those obtained from simulations done using the GEANT3.21 simulation code. The results show that the nucleation parameter depends on the mass number of the incident heavy ions, and decreases with increasing mass number.


title: The expected degree of minimal spanning forests

abstract: We give a lower bound on the expected degree of the free minimal spanning forest of a vertex transitive graph in terms of its spectral radius. This result answers a question of Lyons-Peres-Schramm and simplifies the Gaboriau-Lyons proof of the measurable-group-theoretic solution to von Neumann's problem. In the second part we study a relative version of the free minimal spanning forest. As a consequence of this study we can show that non-torsion unitarizable groups have fixed price one.


title: Kite Pseudo Effect Algebras

abstract: We define a new class of pseudo effect algebras, called kite pseudo effect algebras, which is connected with partially ordered groups not necessarily with strong unit. In such a case, starting even with an Abelian po-group, we can obtain a noncommutative pseudo effect algebra. We show how such kite pseudo effect algebras are tied with different types of the Riesz Decomposition Properties. Kites are so-called perfect pseudo effect algebras, and we define conditions when kite pseudo effect algebras have the least non-trivial normal ideal.


title: Remarks on the asymptotically discretely self-similar solutions of the Navier-Stokes and the Euler equations

abstract: We study scenarios of self-similar type blow-up for the incompressible Navier-Stokes and the Euler equations. The previous notions of the discretely (backward) self-similar solution and the asymptotically self-similar solution are generalized to the locally asymptotically discretely self-similar solution. We prove that there exists no such locally asymptotically discretely self-similar blow-up for the 3D Navier-Stokes equations if the blow-up profile is a time periodic function belonging to $C^1(\Bbb R ; L^3(\Bbb R^3)\cap C^2 (\Bbb R^3))$. For the 3D Euler equations we show that the scenario of asymptotically discretely self-similar blow-up is excluded if the blow-up profile satisfies suitable integrability conditions.


title: Boost Invariant Surfaces with Pointwise 1-Type Gauss Map in Minkowski 4-Space E4-1

abstract: In this paper, we study spacelike rotational surfaces which are called boost invariant surfaces in Minkowski 4-space E41. We give necessary and sufficient condition for flat spacelike rotational surface to have pointwise 1-type Gauss map. Also, we obtain a characterization for boost invariant marginally trapped surface with pointwise 1-type Gauss map.


title: Nonequilibrium dynamical mean-field theory based on weak-coupling perturbation expansions: Application to dynamical symmetry breaking in the Hubbard model

abstract: We discuss the general formalism and validity of weak-coupling perturbation theory as an impurity solver for nonequilibrium dynamical mean-field theory. The method is implemented and tested in the Hubbard model, using expansions up to fourth order for the paramagnetic phase at half filling and third order for the antiferromagnetic and paramagnetic phase away from half filling. We explore various types of weak-coupling expansions and examine the accuracy and applicability of the methods for equilibrium and nonequilibrium problems. We find that in most cases an expansion of local self-energy diagrams including all the tadpole diagrams with respect to the Weiss Green's function (bare-diagram expansion) gives more accurate results than other schemes such as self-consistent perturbation theory using the fully interacting Green's function (bold-diagram expansion). In the paramagnetic phase at half filling, the fourth-order bare expansion improves the result of the second-order expansion in the weak-coupling regime, while both expansions suddenly fail at some intermediate interaction strength. The higher-order bare perturbation is especially advantageous in the antiferromagnetic phase near half filling. We use the third-order bare perturbation expansion within the nonequilibrium dynamical mean-field theory to study dynamical symmetry breaking from the paramagnetic to the antiferromagnetic phase induced by an interaction ramp in the Hubbard model. The results show that the order parameter, after an initial exponential growth, exhibits an amplitude oscillation around a nonthermal value followed by a slow drift toward the thermal value. The transient dynamics seems to be governed by a nonthermal critical point, associated with a nonthermal universality class, which is distinct from the conventional Ginzburg-Landau theory.


title: Probabilistic Solutions to Differential Equations and their Application to Riemannian Statistics

abstract: We study a probabilistic numerical method for the solution of both boundary and initial value problems that returns a joint Gaussian process posterior over the solution. Such methods have concrete value in the statistics on Riemannian manifolds, where non-analytic ordinary differential equations are involved in virtually all computations. The probabilistic formulation permits marginalising the uncertainty of the numerical solution such that statistics are less sensitive to inaccuracies. This leads to new Riemannian algorithms for mean value computations and principal geodesic analysis. Marginalisation also means results can be less precise than point estimates, enabling a noticeable speed-up over the state of the art. Our approach is an argument for a wider point that uncertainty caused by numerical calculations should be tracked throughout the pipeline of machine learning algorithms.


title: Sejong Open Cluster Survey (SOS). 0. Target Selection and Data Analysis

abstract: Star clusters are superb astrophysical laboratories containing cospatial and coeval samples of stars with similar chemical composition. We have initiated the Sejong Open cluster Survey (SOS) - a project dedicated to providing homogeneous photometry of a large number of open clusters in the SAAO Johnson-Cousins' $UBVI$ system. To achieve our main goal, we have paid much attention to the observation of standard stars in order to reproduce the SAAO standard system. Many of our targets are relatively small, sparse clusters that escaped previous observations. As clusters are considered building blocks of the Galactic disk, their physical properties such as the initial mass function, the pattern of mass segregation, etc. give valuable information on the formation and evolution of the Galactic disk. The spatial distribution of young open clusters will be used to revise the local spiral arm structure of the Galaxy. In addition, the homogeneous data can also be used to test stellar evolutionary theory, especially concerning rare massive stars. In this paper we present the target selection criteria, the observational strategy for accurate photometry, and the adopted calibrations for data analysis such as color-color relations, zero-age main sequence relations, Sp - Mv relations, Sp - Teff relations, Sp - color relations, and Teff - BC relations. Finally we provide some data analysis such as the determination of the reddening law, the membership selection criteria, and distance determination.


title: Pre-saddle neutron multiplicity for fission reactions induced by heavy ions and light particles

abstract: Pre-saddle neutron multiplicity has been calculated for several fission reactions induced by heavy ions and light particles. Experimentally, it is impossible to determine the contribution of neutrons being emitted before the saddle point and those emitted between the saddle and the scission points. Determination of the pre-saddle neutron multiplicity in our research is based on the comparison between the experimental anisotropies and those predicted by the standard saddle-point statistical model. Analysis of the results shows that the pre-saddle neutron multiplicity depends on the fission barrier height and stability of the compound nucleus. In heavy ion induced fission, the number of pre-saddle neutrons decreases with increasing the excitation energy of the compound nucleus. A main cause of this behavior is due to a reduction in the ground state-to-saddle point transition time with increasing the excitation energy of the compound nucleus. Whereas in induced fission by light particles, the number of pre-saddle neutrons increases with increasing the excitation energy of the compound nucleus.


title: Bak-Tang-Wiesenfeld Model in the Finite Range Random Link Lattice

abstract: We consider the BTW model in random link lattices with finite range interaction (RLFRI). The degree distribution for nodes is considered to be uniform in the interval $(0,n_0)$. We numerically calculate the exponents of the distribution functions in terms of $(n_0,R)$ in which $R$ is the range of interactions. Dijkstra radius is utilized to calculate the fractal dimension of the avalanches. Our analysis shows that there is, at least one length scale ($r_0(n_0,R)$) in which the fractal dimension is changed. We find that for the scales smaller than $r_0(n_0,R)$, which is typically one decade, the fractal dimension is nearly independent of $n_0$ and $R$ and is equal to $1.4$, i.e. close to that of the BTW in the regular lattice ($1.25$). Using this fact and other analysis, we conclude that the BTW-type behaviors are dominant for small values of $n_0$ and $R$, whereas for large values of these parameters a new regime is seen in which the exponent of distribution function of avalanche masses is nearly $1.4$. We also numerically calculate the explicit form of the \textit{number of unstable nodes} (NUN) as a time dependent process and show that for regular lattice it is (up to a normalization) proportional to a one dimensional Weiner process and for RLFRI it acquires a drift term. Using this dynamical variable it is numerically shown that we can not continuously approach the regular lattice limit by decreasing $R$.