title: Effect of Electrical Properties on Gd modified BiFeO3-PbZrO3

abstract: The 0.5(BiGdxFe1-xO3)-0.5(PbZrO3) composite was synthesized using a high temperature solid-state reaction technique. Preliminary X-ray structural analysis confirms the formation of the composite. The dielectric constant and loss tangent have been studied. The hysteresis loop suggest that the material is lossy. The impedance parameters were studied using an impedance analyzer in a wide range of frequency (102-106 Hz) at different temperatures for all samples. The Nyquist plot suggests the contribution of bulk effect as well as grain boundary effect and the bulk resistance deceases with rise in temperature for all samples. The electrical transport confirms the presence of hopping mechanism in the material. The dc conductivity increases with rise in temperature. The frequency variation of ac conductivity shows that the compound obeys Jonschers universal power law and confirms the Small Polaron (SP) tunneling effect due to low activation energy for all samples. Temperature dependence of dc and ac conductivity indicates that electrical conduction in the materials are thermally activated process.


title: Deriving accurate peculiar velocities (even at high redshift)

abstract: The way that peculiar velocities are often inferred from measurements of distances and redshifts makes an approximation, v_p = cz-H_0 D, that gives significant errors even at relatively low redshifts (overestimates peculiar velocity by ~ 100 km/s at z~0.04). Here we demonstrate where the approximation breaks down, the systematic offset it introduces, and how the exact calculation should be implemented.


title: Strong Stability of Cotangent Bundles of Cyclic Covers

abstract: Let $X$ be a smooth projective variety over an algebraically closed field $k$ of characteristic $p>0$ of $\dim X\geq 4$ and Picard number $\rho(X)=1$. Suppose that $X$ satisfies $H^i(X,F^{m*}_X(\Omg^j_X)\otimes\Ls^{-1})=0$ for any ample line bundle $\Ls$ on $X$, and any nonnegative integers $m,i,j$ with $0\leq i+j<\dim X$, where $F_X:X\rightarrow X$ is the absolute Frobenius morphism. We prove that by procedures combining taking smooth hypersurfaces of dimension $\geq 3$ and cyclic covers along smooth divisors, if the resulting smooth projective variety $Y$ has ample (resp. nef) canonical bundle $\omega_Y$, then $\Omg_Y$ is strongly stable $($resp. strongly semistable$)$ with respect to any polarization.


title: Independence and Matchings in $\sigma$-hypergraphs

abstract: Let $\sigma$ be a partition of the positive integer $r$. A $\sigma$-hypergraph $H=H(n,r,q|\sigma)$ is an $r$-uniform hypergraph on $nq$ vertices which are partitioned into $n$ classes $V_1, V_2, \ldots, V_n$ each containing $q$ vertices. An $r$-subset $K$ of vertices is an edge of the hypergraph if the partition of $r$ formed by the non-zero cardinalities $|K\cap V_i|, 1\leq i \leq n,$ is $\sigma$. In earlier works we have considered colourings of the vertices of $H$ which are constrained such that any edge has at least $\alpha$ and at most $\beta$ vertices of the same colour, and we have shown that interesting results can be obtained by varying $\alpha, \beta$ and the parameters of $H$ appropriately. In this paper we continue to investigate the versatility of $\sigma$-hypergraphs by considering two classical problems: independence and matchings. We first demonstrate an interesting link between the constrained colourings described above and the $k$-independence number of a hypergraph, that is, the largest cardinality of a subset of vertices of a hypergraph not containing $k+1$ vertices in the same edge. We also give an exact computation of the $k$-independence number of the $\sigma$-hypergraph $H$. We then present results on maximum, and sometimes perfect, matchings in $H$. These results often depend on divisibility relations between the parameters of $H$ and on the highest common factor of the parts of $\sigma$.


title: Computing Strong Nash Equilibria for Multiplayer Games

abstract: An heuristic approach to compute strong Nash (Aumann) equilibria is presented. The method is based on differential evolution and three variants of a generative relation for strong Nash equilibria characterization. Numerical experiments performed on the minimum effort game for up to 150 players illustrate the efficiency of the approach. The advantages and disadvantages of each variant is discussed in terms of precision and running time.


title: Estimation of Optimized Energy and Latency Constraints for Task Allocation in 3d Network on Chip

abstract: In Network on Chip (NoC) rooted system, energy consumption is affected by task scheduling and allocation schemes which affect the performance of the system. In this paper we test the pre-existing proposed algorithms and introduced a new energy skilled algorithm for 3D NoC architecture. An efficient dynamic and cluster approaches are proposed along with the optimization using bio-inspired algorithm. The proposed algorithm has been implemented and evaluated on randomly generated benchmark and real life application such as MMS, Telecom and VOPD. The algorithm has also been tested with the E3S benchmark and has been compared with the existing mapping algorithm spiral and crinkle and has shown better reduction in the communication energy consumption and shows improvement in the performance of the system. On performing experimental analysis of proposed algorithm results shows that average reduction in energy consumption is 49%, reduction in communication cost is 48% and average latency is 34%. Cluster based approach is mapped onto NoC using Dynamic Diagonal Mapping (DDMap), Crinkle and Spiral algorithms and found DDmap provides improved result. On analysis and comparison of mapping of cluster using DDmap approach the average energy reduction is 14% and 9% with crinkle and spiral.


title: A Structural Approach to Coordinate-Free Statistics

abstract: We consider the question of learning in general topological vector spaces. By exploiting known (or parametrized) covariance structures, our Main Theorem demonstrates that any continuous linear map corresponds to a certain isomorphism of embedded Hilbert spaces. By inverting this isomorphism and extending continuously, we construct a version of the Ordinary Least Squares estimator in absolute generality. Our Gauss-Markov theorem demonstrates that OLS is a "best linear unbiased estimator", extending the classical result. We construct a stochastic version of the OLS estimator, which is a continuous disintegration exactly for the class of "uncorrelated implies independent" (UII) measures. As a consequence, Gaussian measures always exhibit continuous disintegrations through continuous linear maps, extending a theorem of the first author. Applying this framework to some problems in machine learning, we prove a useful representation theorem for covariance tensors, and show that OLS defines a good kriging predictor for vector-valued arrays on general index spaces. We also construct a support-vector machine classifier in this setting. We hope that our article shines light on some deeper connections between probability theory, statistics and machine learning, and may serve as a point of intersection for these three communities.


title: The wavelet transforms in Gelfand-Shilov spaces

abstract: We describe local and global properties of wavelet transforms of ultradifferentiable functions. The results are given in the form of continuity properties of the wavelet transform on Gelfand-Shilov type spaces and their duals. In particular, we introduce a new family of highly time-scale localized spaces on the upper half-space. We study the wavelet synthesis operator (the left-inverse of the wavelet transform) and obtain the resolution of identity (Calder\'{o}n reproducing formula) in the context of ultradistributions.


title: Upper bounds on the charge susceptibility of many-electron systems coupled to the quantized radiation field

abstract: We extend the Kubo-Kishi theorem concerning the charge susceptibility of the Hubbard model in the following way: (i) The electron-photon interaction is taken into account. (ii) Not only on-site but also general Coulomb repulsion is considered.


title: Sandpile groups of generalized de Bruijn and Kautz graphs and circulant matrices over finite fields

abstract: A maximal minor $M$ of the Laplacian of an $n$-vertex Eulerian digraph $\Gamma$ gives rise to a finite group $\mathbb{Z}^{n-1}/\mathbb{Z}^{n-1}M$ known as the sandpile (or critical) group $S(\Gamma)$ of $\Gamma$. We determine $S(\Gamma)$ of the generalized de Bruijn graphs $\Gamma=\mathrm{DB}(n,d)$ with vertices $0,\dots,n-1$ and arcs $(i,di+k)$ for $0\leq i\leq n-1$ and $0\leq k\leq d-1$, and closely related generalized Kautz graphs, extending and completing earlier results for the classical de Bruijn and Kautz graphs. Moreover, for a prime $p$ and an $n$-cycle permutation matrix $X\in\mathrm{GL}_n(p)$ we show that $S(\mathrm{DB}(n,p))$ is isomorphic to the quotient by $\langle X\rangle$ of the centralizer of $X$ in $\mathrm{PGL}_n(p)$. This offers an explanation for the coincidence of numerical data in sequences A027362 and A003473 of the OEIS, and allows one to speculate upon a possibility to construct normal bases in the finite field $\mathbb{F}_{p^n}$ from spanning trees in $\mathrm{DB}(n,p)$.


title: Hierarchical clustering of DNA k-mer counts in RNA-seq fastq files reveals batch effects

abstract: Batch effects, artificial sources of variation due to experimental design, are a widespread phenomenon in high throughput data. Therefore, mechanisms for detection of batch effects are needed requiring comparison of multiple samples. We apply hierarchical clustering (HC) on DNA k-mer counts of multiple RNA-seq derived Fastq files. Ideally, HC generated trees reflect experimental treatment groups and thus may indicate experimental effects, but clustering of preparation groups indicates the presence of batch effects. In order to provide a simple applicable tool we implemented sequential analysis of Fastq reads with low memory usage in an R package (seqTools) available on Bioconductor. DNA k-mer counts were analysed on 61 Fastq files containing RNA-seq data from two cell types (dermal fibroblasts and Jurkat cells) sequenced on 8 different Illumina Flowcells. Results: Pairwise comparison of all Flowcells with hierarchical clustering revealed strong Flowcell based tree separation in 6 (21 %) and detectable Flowcell based clustering in 17 (60.7 %) of 28 Flowcell comparisons. In our samples, batch effects were also present in reads mapped to the human genome. Filtering reads for high quality (Phred >30) did not remove the batch effects. Conclusions: Hierarchical clustering of DNA k-mer counts provides a quality criterion and an unspecific diagnostic tool for RNA-seq experiments.


title: Kernels in tropical geometry and a Jordan-H\"older Theorem

abstract: A correspondence exists between affine tropical varieties and algebraic objects, following the classical Zariski correspondence between irreducible affine varieties and the prime spectrum of the coordinate algebra in affine algebraic geometry. Although in this context the natural analog of the polynomial ring over a field is the polynomial semiring over a semifield (without a zero element), one obtains homomorphic images of coordinate algebras via congruences rather than ideals, which complicates the algebraic theory considerably. In this paper, we pass to the semifield $F(\lambda_1, \dots, \lambda_n)$ of fractions of the polynomial semiring, for which there already exists a well developed theory of kernels, which are normal convex subgroups; this approach enables us to switch the structural roles of addition and multiplication and makes available much of the extensive theory of chains of homomorphisms of groups, including the Jordan-Holder theory. The parallel of the zero set now is the 1-set. These notions are refined in the language of supertropical algebra to $\nu$-kernels and $1^\nu$-sets, lending more precision to the theory. In analogy to Hilbert's celebrated Nullstellensatz which provides a correspondence between radical ideals and zero sets, we develop a correspondence between $1^\nu$-sets and a well-studied class of $\nu$-kernels of the rational semifield called polars, originating from the theory of lattice-ordered groups. This correspondence becomes simpler and more applicable when restricted to a special kind of kernel, called principal, intersected with the kernel generated by $F$. We utilize this theory to study tropical roots in tropical geometry. As an application, we develop composition series and convexity degree, leading to a tropical version of the Jordan-H\"{o}lder theorem.


title: The Allen-Cahn equation with dynamic boundary conditions and mass constraints

abstract: The Allen-Cahn equation, coupled with dynamic boundary conditions, has recently received a good deal of attention. The new issue of this paper is the setting of a rather general mass constraint which may involve either the solution inside the domain or its trace on the boundary. The system of nonlinear partial differential equations can be formulated as variational inequality. The presence of the constraint in the evolution process leads to additional terms in the equation and the boundary condition containing a suitable Lagrange multiplier. A well-posedness result is proved for the related initial value problem.


title: The Deep Full-Stokes Radio Sky

abstract: The new broad-band capabilities of large radio interferometers such as the GMRT and JVLA allow for long-integration mosaic imaging observations to create ultra-deep full-polarization images of the sky over wide frequency ranges. Achieving rms sensitivities of order 1 $\mu$Jy, these observations explore the radio source population at flux densities well below the regime dominated by classical radio galaxies and Active Galactic Nuclei. We present initial results from radio sources revealed with deep mosaicking observations with the GMRT and JVLA at respectively 0.6 and 5 GHz, and evidence that the $\mu$Jy sensitivity level marks the transition to detection of polarized emission from a population of sources dominated by emission from magnetic fields in the disks of starburst and normal galaxies.


title: Classical reactive scattering in a quantum spirit: Improving the shape of rotational state distributions in the quantum regime

abstract: For triatomic chemical reactions under single-collision conditions, we propose a new quasi-classical trajectory (QCT) approach to rotational-state distributions of particular interest in the quantum regime where only a few rotational states are available to the products. Our method is directly inspired from the amendments to be introduced in classical phase space theory (PST) in order to make it in exact agreement with quantum PST. The method is applied to the D$^+$ + H$_2$ and H$^+$ + D$_2$ reactions and the population of the rotational ground state is found to be in much closer agreement with the exact quantum one than the same population obtained by means of standard QCT calculations. The impact on the whole distribution is all the stronger as the number of available states is small. Last but not least, the shape of the distribution appears to be controlled to a large extent by three factors, respectively called parity, edge and rotational shift factors.


title: Time-resolved Spectral Study Of Fermi GRBs Having Single Pulses

abstract: We analyze gamma-ray bursts (GRBs) detected by \textit{Fermi}/Gamma-Ray Burst Monitor (GBM) and having single pulse. We fit the light curves with a model having exponential rise and decay parts. We perform a detailed time-resolved spectroscopy using four models: Band, blackbody with a power-law (BBPL), multicolour blackbody with a power-law (mBBPL) and two blackbodies with a power-law (2BBPL). We find that models other than the BBPL give better $\chi_{red}^2$ for the "hard-to-soft" (HTS) pulses, while for the "intensity tracking" (IT) pulses, the BBPL model is statistically as good as the other models. Interestingly, the energy at the peak of the spectrum resulting from the BBPL model ($\sim3kT$), is always lower than that of the $\nu F_{\nu}$ spectrum of the Band function. The values of the low energy photon index ($\alpha$) of the Band function are often higher than the fundamental single particle synchrotron limit, especially for the HTS pulses. Specifically we find two extreme cases --- for GRB~110817A (HTS GRB) $\alpha$ is always higher, while for GRB~100528A (IT GRB) $\alpha$ is always within the synchrotron regime. The PL component of the BBPL model always starts with a delay compared to the BB component, and it lingers at the later part of the prompt emission phase. For three HTS GRBs, namely, GRB~081224, GRB~100707A and GRB~110721A this behaviour is particularly significant and interestingly there are reported LAT detections for them. Finally, we argue that various evidences hint that neither BBPL nor Band model is acceptable, while 2BBPL and mBBPL are the most acceptable models for the set of GRBs we have analyzed.


title: Global existence for semilinear wave equations with the critical blow-up term in high dimensions

abstract: We are interested in almost global existence cases in the general theory for nonlinear wave equations, which are caused by critical exponents of nonlinear terms. Such situations can be found in only three cases in the theory, cubic terms in two space dimensions, quadratic terms in three space dimesions and quadratic terms including a square of unknown functions itself in four space dimensions. Except for the last case, criterions to classify nonlinear terms into the almost global, or global existence case, are well-studied and known to be so-called null condition and non-positive condition. Our motivation of this work is to find such a kind of the criterion in four space dimensions. In our previous paper, an example of the non-single term for the almost global existence case is introduced. In this paper, we show an example of the global existence case. These two examples have nonlinear integral terms which are closely related to derivative loss due to high dimensions. But it may help us to describe the final form of the criterion.


title: Postulation of disjoint unions of lines and a multiple point, II

abstract: We study the postulation of a general union $X\subset \mathbb {P}^3$ of one m-point $mP$ and $t$ disjoint lines. We prove that it has the expected Hilbert function, proving a conjecture by E. Carlini, M. V. Catalisano and A. V. Geramita.


title: Newton's problem of minimal resistance under the single-impact assumption

abstract: A parallel flow of non-interacting point particles is incident on a body at rest. When hitting the body's surface, the particles are reflected elastically. Assume that each particle hits the body at most once (SIC condition); then the force of resistance of the body along the flow direction can be written down in a simple analytical form. The problem of minimal resistance within this model was first considered by Newton (1687) in the class of bodies with a fixed length M along the flow direction and with a fixed maximum orthogonal cross section, under the additional conditions that the body is convex and rotationally symmetric. Here we solve the problem (first stated by Ferone, Buttazzo, and Kawohl in 1995) for the wider class of bodies satisfying SIC and with the additional conditions removed. The scheme of solution is inspired by Besicovitch's method of solving the Kakeya problem. If the maximum cross section is a disc, the decrease of resistance as compared with the original Newton problem is more than twofold; the ratio tends to 2 as M goes to 0 and to 81/4 as M goes to infinity. We also prove that the infimum of resistance is 0 for a wider class of bodies with both single and double impacts allowed.


title: On the notion of 'retractable modules' in the context of algebras

abstract: This is a survey on the usage of the module theoretic notion of a "retractable module" in the study of algebras with actions. We explain how classical results can be interpreted using module theory and end the paper with some open questions.


title: Exploring the phase diagram of fully turbulent Taylor-Couette flow

abstract: Direct numerical simulations of Taylor-Couette flow (TC). Shear Reynolds numbers of up to $3\cdot10^5$, corresponding to Taylor numbers of $Ta=4.6\cdot10^{10}$, were reached. Effective scaling laws for the torque are presented. The transition to the ultimate regime, in which asymptotic scaling laws (with logarithmic corrections) for the torque are expected to hold up to arbitrarily high driving, is analysed for different radius ratios $\eta$, different aspect ratios $\Gamma$ and different rotation ratios $Ro$. It is shown that the transition is approximately independent of $Ro$ and $\Gamma$, but depends significantly on $\eta$. We furthermore calculate the local angular velocity profiles and visualize different flow regimes that depend both on the shearing of the flow, and the Coriolis force originating from the outer cylinder rotation. Two main regimes are distinguished, based on the magnitude of the Coriolis force, namely the co-rotating and weakly counter-rotating regime dominated by Rayleigh-unstable regions, and the strongly counter-rotating regime where a mixture of stable and unstable regions exist. Furthermore, an analogy between $\eta$ and outer-cylinder rotation is revealed, namely that smaller gaps behave like a wider gap with co-rotating cylinders, and that wider gaps behave like smaller gaps with weakly counter-rotating cylinders. Finally, the effect of $\Gamma$ on the effective torque versus $Ta$ scaling is analysed and it is shown that different branches of the torque-versus-$Ta$ relationships associated to different aspect ratios are found to cross within $15%$ of the $Re$ associated to the transition to the ultimate regime. The paper culminates in phase diagram in the inner vs outer $Re$ number parameter space and in the $Ta$ vs $Ro$ parameter space, which can be seen as the extension of the Andereck \emph{et al.} phase diagram towards the ultimate regime.


title: Implicit yield function formulation for granular and rock-like materials

abstract: The constitutive modelling of granular, porous and quasi-brittle materials is based on yield (or damage) functions, which may exhibit features (for instance, lack of convexity, or branches where the values go to infinity, or false elastic domains) preventing the use of efficient return-mapping integration schemes. This problem is solved by proposing a general construction strategy to define an implicitly defined convex yield function starting from any convex yield surface. Based on this implicit definition of the yield function, a return-mapping integration scheme is implemented and tested for elastic-plastic (or -damaging) rate equations. The scheme is general and, although it introduces a numerical cost when compared to situations where the scheme is not needed, is demonstrated to perform correctly and accurately.


title: Is Consciousness Computable? Quantifying Integrated Information Using Algorithmic Information Theory

abstract: In this article we review Tononi's (2008) theory of consciousness as integrated information. We argue that previous formalizations of integrated information (e.g. Griffith, 2014) depend on information loss. Since lossy integration would necessitate continuous damage to existing memories, we propose it is more natural to frame consciousness as a lossless integrative process and provide a formalization of this idea using algorithmic information theory. We prove that complete lossless integration requires noncomputable functions. This result implies that if unitary consciousness exists, it cannot be modelled computationally.


title: On the minimization of Dirichlet eigenvalues

abstract: Results are obtained for two minimization problems: $$I_k(c)=\inf \{\lambda_k(\Omega): \Omega\ \textup{open, convex in}\ \mathbb{R}^m,\ \mathcal{T}(\Omega)= c \},$$ and $$J_k(c)=\inf\{\lambda_k(\Omega): \Omega\ \textup{quasi-open in}\ \mathbb{R}^m, |\Omega|\le 1, \mathcal {P}(\Omega)\le c \},$$ where $c>0$, $\lambda_k(\Omega)$ is the $k$'th eigenvalue of the Dirichlet Laplacian acting in $L^2(\Omega)$, $|\Omega|$ denotes the Lebesgue measure of $\Omega$, $\mathcal{P}(\Omega)$ denotes the perimeter of $\Omega$, and where $\mathcal{T}$ is in a suitable class set functions. The latter include for example the perimeter of $\Omega$, and the moment of inertia of $\Omega$ with respect to its center of mass.


title: On the proper use of the Schwarzschild and Ledoux criteria in stellar evolution computations

abstract: The era of detailed asteroseismic analyses opened by space missions such as CoRoT and $\textit{Kepler}$ has highlighted the need for stellar models devoid of numerical inaccuracies, in order to be able to diagnose which physical aspects are being ignored or poorly treated in standard stellar modeling. We tackle here the important problem of fixing convective zones boundaries in the frame of the local mixing length theory. First we show that the only correct way to locate a convective zone boundary is to find, at each iteration step, through interpolations or extrapolations from points $\textit{within the convective zone}$, the mass where the radiative luminosity is equal to the total one. We then discuss two misuses of the boundary condition and the way they affect stellar modeling and stellar evolution. The first one consists in applying the neutrality condition for convective instability on the $\textit{radiative}$ side of the convective boundary. The second way of misusing the boundary condition comes from the process of fixing the convective boundary through the search for a change of sign of a possibly \textit{discontinuous} function. We show that these misuses can lead to completely wrong estimates of convective core sizes with important consequences for the following evolutionary phases. We point out the advantages of using a double mesh point at each convective zone boundaries. The specific problem of a convective shell is discussed and some remarks concerning overshooting are given.


title: Large-area Epitaxial Monolayer MoS2

abstract: Two-dimensional semiconductors such as MoS2 are an emerging material family with wide-ranging potential applications in electronics, optoelectronics and energy harvesting. Large-area growth methods are needed to open the way to the applications. While significant progress to this goal was made, control over lattice orientation during growth still remains a challenge. This is needed in order to minimize or even avoid the formation of grain boundaries which can be detrimental to electrical, optical and mechanical properties of MoS2 and other 2D semiconductors. Here, we report on the uniform growth of high-quality centimeter-scale continuous monolayer MoS2 with control over lattice orientation. Using transmission electron microscopy we show that the monolayer film is composed of coalescing single islands that share a predominant lattice orientation due to an epitaxial growth mechanism. Raman and photoluminescence spectra confirm the high quality of the grown material. Optical absorbance spectra acquired over large areas show new features in the high-energy part of the spectrum, indicating that MoS2 could also be interesting for harvesting this region of the solar spectrum and fabrication of UV-sensitive photodetectors. Even though the interaction between the growth substrate and MoS2 is strong enough to induce lattice alignment, we can easily transfer the grown material and fabricate field-effect transistors on SiO2 substrates showing mobility superior to the exfoliated material.


title: Energy Loss of Solar $p$ Modes due to the excitation of Magnetic Sausage Tube Waves: Importance of Coupling the Upper Atmosphere

abstract: We consider damping and absorption of solar $p$ modes due to their energy loss to magnetic tube waves that can freely carry energy out of the acoustic cavity. The coupling of $p$ modes and sausage tube waves is studied in a model atmosphere composed of a polytropic interior above which lies an isothermal upper atmosphere. The sausage tube waves, excited by $p$ modes, propagate along a magnetic fibril which is assumed to be a vertically aligned, stratified, thin magnetic flux-tube. The deficit of $p$-mode energy is quantified through the damping rate, $\Gamma$ and absorption coefficient, $\alpha$. The variation of $\Gamma$ and $\alpha$ as a function of frequency and the tube's plasma properties is studied in detail. Previous similar studies have considered only a subphotospheric layer, modelled as a polytrope that has been truncated at the photosphere (Bogdan et al. (1996), Hindman & Jain 2008, Gascoyne et al. (2011)). Such studies have found that the resulting energy loss by the $p$ modes is very sensitive to the upper boundary condition, which because of the lack of a upper atmosphere have been imposed in a somewhat ad hoc manner. The model presented here avoids such problems by using an isothermal layer to model the overlying atmosphere (chromosphere), and consequently, allows us to analyse the propagation of $p$-mode driven sausage waves above the photosphere. In this paper we restrict our attention to frequencies below the acoustic cut-off frequency. We demonstrate the importance of coupling all waves (acoustic, magnetic) in the subsurface solar atmosphere with the overlying atmosphere in order to accurately model the interaction of solar $f$ and $p$ modes with sausage tube waves.


title: Two Procedures for Robust Monitoring of Probability Distributions of Economic Data Streams induced by Depth Functions

abstract: Data streams (streaming data) consist of transiently observed, evolving in time, multidimensional data sequences that challenge our computational and/or inferential capabilities. In this paper we propose user friendly approaches for robust monitoring of selected properties of unconditional and conditional distribution of the stream basing on depth functions. Our proposals are robust to a small fraction of outliers and/or inliers but sensitive to a regime change of the stream at the same time. Their implementations are available in our free R package DepthProc.


title: Slow-light enhanced gain in active photonic crystal waveguides

abstract: Slow light is a fascinating physical effect, raising fundamental questions related to our understanding of light-matter interactions as well as offering new possibilities for photonic devices. From the first demonstrations of slow light propagation in ultra-cold atomic gasses, solid-state Ruby and photonic crystal structures, focus has shifted to applications, with slow light offering the ability to enhance and control light-matter interactions. The demonstration of tuneable delay lines, enhanced nonlinearities and spontaneous emission, enlarged spectral sensitivity and increased phase shifts illustrate the possibilities enabled by slow light propagation, with microwave photonics emerging as one of the promising applications. Here, we demonstrate that slow light can be used to control and increase the gain coefficient of an active semiconductor waveguide. The effect was theoretically predicted but not yet experimentally demonstrated. These results show a route towards realizing ultra-compact optical amplifiers for linear and nonlinear applications in integrated photonics and prompts further research into the rich physics of such structures.


title: Geodesic Distance Function Learning via Heat Flow on Vector Fields

abstract: Learning a distance function or metric on a given data manifold is of great importance in machine learning and pattern recognition. Many of the previous works first embed the manifold to Euclidean space and then learn the distance function. However, such a scheme might not faithfully preserve the distance function if the original manifold is not Euclidean. Note that the distance function on a manifold can always be well-defined. In this paper, we propose to learn the distance function directly on the manifold without embedding. We first provide a theoretical characterization of the distance function by its gradient field. Based on our theoretical analysis, we propose to first learn the gradient field of the distance function and then learn the distance function itself. Specifically, we set the gradient field of a local distance function as an initial vector field. Then we transport it to the whole manifold via heat flow on vector fields. Finally, the geodesic distance function can be obtained by requiring its gradient field to be close to the normalized vector field. Experimental results on both synthetic and real data demonstrate the effectiveness of our proposed algorithm.


title: Input-to-State Stability, integral Input-to-State Stability, and $\mathcal{L}_2$-Gain Properties: Qualitative Equivalences and Interconnected Systems

abstract: Input-to-state stability (ISS) and $\mathcal{L}_2$-gain are well-known robust stability properties that continue to find wide application in the analysis and control of nonlinear dynamical systems and their interconnections. We investigate the relationship between ISS-type and $\mathcal{L}_2$-gain properties, demonstrating several qualitative equivalences between these two approaches. We subsequently present several new sufficient conditions for the stability of interconnected systems derived by exploiting these qualitative equivalences.


title: Separated design of encoder and controller for networked linear quadratic optimal control

abstract: For a networked control system, we consider the problem of encoder and controller design. We study a discrete-time linear plant with a finite horizon performance cost, comprising of a quadratic function of the states and controls, and an additive communication cost. We study separation in design of the encoder and controller, along with related closed-loop properties such as the dual effect and certainty equivalence. We consider three basic formats for encoder outputs: quantized samples, real-valued samples at event-triggered times, and real-valued samples over additive noise channels. If the controller and encoder are dynamic, then we show that the performance cost is minimized by a separated design: the controls are updated at each time instant as per a certainty equivalence law, and the encoder is chosen to minimize an aggregate quadratic distortion of the estimation error. This separation is shown to hold even though a dual effect is present in the closed-loop system. We also show that this separated design need not be optimal when the controller or encoder are to be chosen from within restricted classes.


title: Value sharing by an entire function with its derivatives

abstract: We prove a uniqueness theorem for an entire function, which shares certain values with its higher order derivatives.


title: On the informational completeness of local observables

abstract: For a general multipartite quantum state, we formulate a locally checkable condition, under which the expectation values of certain nonlocal observables are completely determined by the expectation values of some local observables. The condition is satisfied by ground states of gapped quantum many-body systems in two spatial dimensions, assuming a widely conjectured form of area law is correct. Its implications on the studies of gapped quantum many-body systems, quantum state tomography, and quantum state verification are discussed. These results are based on a partial characterization of states with small yet nonzero conditional mutual information, which may be of independent interest.

categories: quant-ph cond-mat.str-el math-ph math.MP


title: Corrected phase-type approximations of heavy-tailed queueing models in a Markovian environment

abstract: Significant correlations between arrivals of load-generating events make the numerical evaluation of the workload of a system a challenging problem. In this paper, we construct highly accurate approximations of the workload distribution of the MAP/G/1 queue that capture the tail behavior of the exact workload distribution and provide a bounded relative error. Motivated by statistical analysis, we consider the service times as a mixture of a phase-type and a heavy-tailed distribution. With the aid of perturbation analysis, we derive our approximations as a sum of the workload distribution of the MAP/PH/1 queue and a heavy-tailed component that depends on the perturbation parameter. We refer to our approximations as corrected phase-type approximations, and we exhibit their performance with a numerical study.


title: Cosmological bounces in spatially flat FRW spacetimes in metric f(R) gravity

abstract: The present work analyzes the various conditions in which there can be a bouncing universe solution in f(R) gravity. In the article an interesting method, to analyze the bouncing FRW solutions in a spatially flat universe using f(R) gravity models using an effective Einstein frame description of the process, is presented. The analysis shows that a cosmological bounce in the f(R) theory need not be described by an equivalent bounce in the Einstein frame description of the process where actually there may be no bounce at all. Nevertheless the Einstein frame description of the bouncing phenomena turns out to be immensely important as the dynamics of the bounce becomes amenable to logic based on general relativistic intuition. The theory of scalar cosmological perturbations in the bouncing universe models in $f(R)$ theories has also been worked out in the Einstein frame.


title: Understanding the $\mu$SR spectra of MnSi without magnetic polarons

abstract: Transverse-field muon-spin rotation ($\mu$SR) experiments were performed on a single crystal sample of the non-centrosymmetric system MnSi. The observed angular dependence of the muon precession frequencies matches perfectly the one of the Mn-dipolar fields acting on the muons stopping at a 4a position of the crystallographic structure. The data provide a precise determination of the magnetic dipolar tensor. In addition, we have calculated the shape of the field distribution expected below the magnetic transition temperature $T_C$ at the 4a muon-site when no external magnetic field is applied. We show that this field distribution is consistent with the one reported by zero-field $\mu$SR studies. Finally, we present ab initio calculations based on the density-functional theory which confirm the position of the muon stopping site inferred from transverse-field $\mu$SR. In view of the presented evidence we conclude that the $\mu$SR response of MnSi can be perfectly and fully understood without invoking a hypothetical magnetic polaron state.


title: Corrected phase-type approximations for the workload of the MAP/G/1 queue with heavy-tailed service times

abstract: In many applications, significant correlations between arrivals of load-generating events make the numerical evaluation of the load of a system a challenging problem. Here, we construct very accurate approximations of the workload distribution of the MAP/G/1 queue that capture the tail behavior of the exact workload distribution and provide a small relative error. Motivated by statistical analysis, we assume that the service times are a mixture of a phase-type and a heavy-tailed distribution. With the aid of perturbation analysis, we derive our approximations as a sum of the workload distribution of the MAP/PH/1 queue and a heavy-tailed component that depends on the perturbation parameter. We refer to our approximations as corrected phase-type approximations, and we exhibit their performance with a numerical study.


title: Asymptotic behavior of a relativistic diffusion in Robertson-Walker space-times

abstract: We determine the long-time asymptotic behavior of a relativistic diffusion taking values in the unitary tangent bundle of a Robertson-Walker space-time. We prove in particular that when approaching the explosion time of the diffusion, its projection on the base manifold almost surely converges to a random point of the causal boundary and we also describe the behavior of the tangent vector in the neighborhood of this limiting point.


title: An infinite family of prime knots with a certain property for the clasp number

abstract: The clasp number $c(K)$ of a knot $K$ is the minimum number of clasp singularities among all clasp disks bounded by $K$. It is known that the genus $g(K)$ and the unknotting number $u(K)$ are lower bounds of the clasp number, that is, $\max\{g(K),u(K)\} \leq c(K)$. Then it is natural to ask whether there exists a knot $K$ such that $\max\{g(K),u(K)\}<c(K)$. In this paper, we prove that there exists an infinite family of prime knots such that the question above is affirmative.


title: Discrete-Time Fractional-Order PID Controller: Definition, Tuning, Digital Realization and Experimental Results

abstract: In some of the complicated control problems we have to use the controllers that apply nonlocal operators to the error signal to generate the control. Currently, the most famous controller with nonlocal operators is the fractional-order PID (FOPID). Commonly, after tuning the parameters of FOPID controller, its transfer function is discretized (for realization purposes) using the so-called generating function. This discretization is the origin of some errors and unexpected results in feedback systems. It may even happen that the controller obtained by discretizing a FOPID controller works worse than a directly-tuned discrete-time classical PID controller. Moreover, FOPID controllers cannot directly be applied to the processes modeled by, e.g., the ARMA or ARMAX model. The aim of this paper is to propose a discrete-time version of the FOPID controller and discuss on its properties and applications. Similar to the FOPID controller, the proposed structure applies nonlocal operators (with adjustable memory length) to the error signal. Two methods for tuning the parameters of the proposed controller are developed and it is shown that the proposed controller has the capacity of solving complicated control problems with a high performance.


title: Contextual Semantic Parsing using Crowdsourced Spatial Descriptions

abstract: We describe a contextual parser for the Robot Commands Treebank, a new crowdsourced resource. In contrast to previous semantic parsers that select the most-probable parse, we consider the different problem of parsing using additional situational context to disambiguate between different readings of a sentence. We show that multiple semantic analyses can be searched using dynamic programming via interaction with a spatial planner, to guide the parsing process. We are able to parse sentences in near linear-time by ruling out analyses early on that are incompatible with spatial context. We report a 34% upper bound on accuracy, as our planner correctly processes spatial context for 3,394 out of 10,000 sentences. However, our parser achieves a 96.53% exact-match score for parsing within the subset of sentences recognized by the planner, compared to 82.14% for a non-contextual parser.


title: Moment Asymptotic Expansions of the Wavelet Transforms

abstract: Using distribution theory we present the moment asymptotic expansion of continuous wavelet transform in different distributional spaces for large and small values of dilation parameter $a$. We also obtain asymptotic expansions for certain wavelet transform.


title: Receivers for Diffusion-Based Molecular Communication: Exploiting Memory and Sampling Rate

abstract: In this paper, a diffusion-based molecular communication channel between two nano-machines is considered. The effect of the amount of memory on performance is characterized, and a simple memory-limited decoder is proposed and its performance is shown to be close to that of the best possible imaginable decoder (without any restriction on the computational complexity or its functional form), using Genie-aided upper bounds. This effect is specialized for the case of Molecular Concentration Shift Keying; it is shown that a four-bits memory achieved nearly the same performance as infinite memory. Then a general class of threshold decoders is considered and shown not to be optimal for Poisson channel with memory, unless SNR is higher than a value specified in the paper. Another contribution is to show that receiver sampling at a rate higher than the transmission rate, i.e., a multi-read system, can significantly improve the performance. The associated decision rule for this system is shown to be a weighted sum of the samples during each symbol interval. The performance of the system is analyzed using the saddle point approximation. The best performance gains are achieved for an oversampling factor of three.


title: Poisson boundary of a relativistic diffusion in curved space-times: an example

abstract: We study in details the long-time asymptotic behavior of a relativistic diffusion taking values in the unitary tangent bundle of a curved Lorentzian manifold, namely a spatially flat and fast expanding Robertson-Walker space-time. We prove in particular that the Poisson boundary of the diffusion can be identified with the causal boundary of the underlying manifold.


title: Coding Theoretic Construction of Quantum Ramp Secret Sharing

abstract: We show a construction of a quantum ramp secret sharing scheme from a nested pair of linear codes. Necessary and sufficient conditions for qualified sets and forbidden sets are given in terms of combinatorial properties of nested linear codes. An algebraic geometric construction for quantum secret sharing is also given.


title: On complexity of multiplication in finite soluble groups

abstract: We determine a reasonable upper bound for the complexity of collection from the left to multiply two elements of a finite soluble, or polycyclic, group by restricting attention to certain polycyclic presentations of the group.


title: On the collapse of trial solutions for a damped-driven non-linear Schr\"odinger equation

abstract: We consider the focusing 2D non-linear Schr\"odinger equation, perturbed by a damping term, and driven by multiplicative noise. We show that a physically motivated trial solution does not collapse for any admissible initial condition although the exponent of the non-linearity is critical. Our method is based on the construction of a global solution to a singular stochastic Hamiltonian system used to connect trial solution and Schr\"odinger equation.


title: Scaling and Inverse Scaling in Anisotropic Bootstrap percolation

abstract: In bootstrap percolation it is known that the critical percolation threshold tends to converge slowly to zero with increasing system size, or, inversely, the critical size diverges fast when the percolation probability goes to zero. To obtain higher-order terms (that is, sharp and sharper thresholds) for the percolation threshold in general is a hard question. In the case of two-dimensional anisotropic models, sometimes correction terms can be obtained from inversion in a relatively simple manner.

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


title: Test of the $h_1(1830)$ made of $K^{*}\bar{K}^*$ with the $\eta_c\rightarrow \phi K^{*}\bar{K}^{*}$ decay

abstract: We present a new reaction, complementary to the $J/\psi\rightarrow \eta K^{*0}\bar{K}^{*0}$ from which an $h_1$ resonance with mass around $1830$ MeV was reported from a BESIII experiment. The new reaction is $\eta_c\rightarrow \phi K^*\bar{K}^*$, or $\eta_c(2S)\rightarrow \phi K^*\bar{K}^*$. Using the information from the analysis of $J/\psi\rightarrow \eta K^{*0}\bar{K}^{*0}$, we find that the $K^*\bar{K}^*$ invariant mass distribution for those two $\eta_c$ decays exhibits a clear peak around $1830$ MeV perfectly distinguishable from what one obtains with pure phase space. We suggest the implementation of these reactions to assert the existence of this elusive resonance which, by its nature as a vector-vector molecule with $0^-(1^{+-})$ quantum numbers, only couples to the $K^*\bar{K}^*$ channel.


title: Sine-Gordon equation in higher dimensions: A fresh look at integrability

abstract: The Sine-Gordon equation is integrable in (1+1)-dimensional Minkowski and in 2-dimensional Euclidean spaces. In each case, it has a Lax pair, and a Hirota algorithm generates its N soliton solutions for all N greater than or equal to 1. The (1+2)-dimensional equation does not pass known integrability tests and does not have a Lax pair. Still, the Hirota algorithm generates N soliton solutions of that equation for all N greater than or equal to 1. Each multi-soliton solution propagates rigidly at a constant velocity, v. The solutions are divided into two unconnected subspaces: Solutions with v greater than or equal to c =1, and v smaller than c. Each subspace is connected by an invertible transformation (rotation plus dilation) to the space of soliton solutions of an integrable Sine-Gordon equation in two dimensions. The faster-than-light solutions are connected to the solutions in (1+1)-dimensional Minkowski space. The slower-than-light solutions are connected to the solutions in 2-dimensional Euclidean space by this transformation and also by Lorentz transformations. The Sine-Gordon equation in (1+3)-dimensional Minkowski space has a richer variety of solutions. Its slower-than-light solutions are connected to the solutions of the integrable equation in 2-dimensional Euclidean space. However, only a subset of its faster-than-light solutions is connected to the solutions of the integrable equation in (1+1)-dimensional Minkowski space.


title: Asteroseismic measurement of surface-to-core rotation in a main sequence A star, KIC 11145123

abstract: We have discovered rotationally split core g-mode triplets and surface p-mode triplets and quintuplets in a terminal age main sequence A star, KIC 11145123, that shows both $\delta$ Sct p-mode pulsations and $\gamma$ Dor g-mode pulsations. This gives the first robust determination of the rotation of the deep core and surface of a main sequence star, essentially model-independently. We find its rotation to be nearly uniform with a period near 100 d, but we show with high confidence that the surface rotates slightly faster than the core. A strong angular momentum transfer mechanism must be operating to produce the nearly rigid rotation, and a mechanism other than viscosity must be operating to produce a more rapidly rotating surface than core. Our asteroseismic result, along with previous asteroseismic constraints on internal rotation in some B stars, and measurements of internal rotation in some subgiant, giant and white dwarf stars, has made angular momentum transport in stars throughout their lifetimes an observational science.


title: Density scaling and quasiuniversality of flow-event statistics for athermal plastic flows

abstract: Athermal plastic flows were simulated for the Kob-Andersen binary Lennard-Jones system and its repulsive version in which the sign of the attractive terms is changed to a plus. Properties evaluated from simulations at different densities include the distributions of energy drops, stress drops, and strain intervals between the flow events. By reference to hidden scale invariance we show that simulations at a single density in conjunction with an equilibrium-liquid simulation at the same density allows one to predict the plastic flow-event properties at other densities. We furthermore demonstrate quasiuniversality of the flow-event statistics.


title: Dimensionality of social networks using motifs and eigenvalues

abstract: We consider the dimensionality of social networks, and develop experiments aimed at predicting that dimension. We find that a social network model with nodes and links sampled from an $m$-dimensional metric space with power-law distributed influence regions best fits samples from real-world networks when $m$ scales logarithmically with the number of nodes of the network. This supports a logarithmic dimension hypothesis, and we provide evidence with two different social networks, Facebook and LinkedIn. Further, we employ two different methods for confirming the hypothesis: the first uses the distribution of motif counts, and the second exploits the eigenvalue distribution.


title: On CR Paneitz operators and CR pluriharmonic functions

abstract: Let $(X,T^{1,0}X)$ be a compact orientable embeddable three dimensional strongly pseudoconvex CR manifold and let ${\rm P\,}$ be the associated CR Paneitz operator. In this paper, we show that (I) ${\rm P\,}$ is self-adjoint and ${\rm P\,}$ has $L^2$ closed range. Let $N$ and $\Pi$ be the associated partial inverse and the orthogonal projection onto ${\rm Ker\,}{\rm P\,}$ respectively, then $N$ and $\Pi$ enjoy some regularity properties. (II) Let $\hat{\mathcal{P}}$ and $\hat{\mathcal{P}_0}$ be the space of $L^2$ CR pluriharmonic functions and the space of real part of $L^2$ global CR functions respectively. Let $S$ be the associated Szeg\"o projection and let $\tau$, $\tau_0$ be the orthogonal projections onto $\hat{\mathcal{P}}$ and $\hat{\mathcal{P}_0}$ respectively. Then, $\Pi=S+\ol S+F_0$, $\tau=S+\ol S+F_1$, $\tau_0=S+\ol S+F_2$, where $F_0, F_1, F_2$ are smoothing operators on $X$. In particular, $\Pi$, $\tau$ and $\tau_0$ are Fourier integral operators with complex phases and $\hat{\mathcal{P}}^\perp\bigcap{\rm Ker\,}{\rm P\,}$, $\hat{\mathcal{P}_0}^\perp\bigcap\hat{\mathcal{P}}$, $\hat{\mathcal{P}_0}^\perp\bigcap{\rm Ker\,}{\rm P\,}$ are all finite dimensional subspaces of $C^\infty(X)$ (it is well-known that $\hat{\mathcal{P}_0}\subset\hat{\mathcal{P}}\subset{\rm Ker\,}{\rm P\,}$). (III) ${\rm Spec\,}{\rm P\,}$ is a discrete subset of $\Real$ and for every $\lambda\in{\rm Spec\,}{\rm P\,}$, $\lambda\neq0$, $\lambda$ is an eigenvalue of ${\rm P\,}$ and the associated eigenspace $H_\lambda({\rm P\,})$ is a finite dimensional subspace of $C^\infty(X)$.


title: Scale-dependent hemispherical asymmetry from general initial state during inflation

abstract: We consider a general initial state for inflation as the mechanism for generating scale-dependent hemispherical asymmetry. An observable scale-dependent non-Gaussianity is generated that leads to observable hemispherical asymmetry from the super-horizon long mode modulation. We show that the amplitude of dipole asymmetry falls off exponentially on small angular scales which can address the absence of dipole asymmetry at these scales.In addition, depending on the nature of non-vaccum initial state, the amplitude of the dipole asymmetry has oscillatory features which can be detected in a careful CMB map analysis. Furthermore, we show that the non-vacuum initial state provides a natural mechanism for enhancing the super horizon long mode perturbation as required to generate the dipole asymmetry.


title: BiSON data preparation: A correction for differential extinction and the weighted averaging of contemporaneous data

abstract: The Birmingham Solar Oscillations Network (BiSON) has provided high-quality high-cadence observations from as far back in time as 1978. These data must be calibrated from the raw observations into radial velocity and the quality of the calibration has a large impact on the signal-to-noise ratio of the final time series. The aim of this work is to maximise the potential science that can be performed with the BiSON data set by optimising the calibration procedure. To achieve better levels of signal-to-noise ratio we perform two key steps in the calibration process: we attempt a correction for terrestrial atmospheric differential extinction; and the resulting improvement in the calibration allows us to perform weighted averaging of contemporaneous data from different BiSON stations. The improvements listed produce significant improvement in the signal-to-noise ratio of the BiSON frequency-power spectrum across all frequency ranges. The reduction of noise in the power spectrum will allow future work to provide greater constraint on changes in the oscillation spectrum with solar activity. In addition, the analysis of the low-frequency region suggests we have achieved a noise level that may allow us to improve estimates of the upper limit of g-mode amplitudes.


title: Topics In Primitive Roots

abstract: This monograph considers a few topics in the theory of primitive roots g(p) modulo a prime p>=2. A few estimates of the least primitive roots g(p) and the least prime primitive roots g^*(p) modulo p, a large prime, are determined. One of the estimate here seems to sharpen the Burgess estimate g(p) << p^(1/4+e) for arbitrarily small number 3 > 0, to the smaller estimate g(p) <= p^(5/loglog p) uniformly for all large primes p => 2. The expected order of magnitude is g(p) <<(log p)^c, c>1 constant. The corresponding estimates for least prime primitive roots g^*(p) are slightly higher. Anotrher topic deals with an effective lower bound #{p <= x : ord(g)= p-1} >> x/log x for the number of primes p <= x with a fixed primitive root g != -1, b^2 for all large number x >1. The current results in the literature claim the lower bound #{p <= x : ord(g) = p-1} >> x/(log x)^2, and have restrictions on the minimal number of fixed integers to three or more.


title: A note on semilinear elliptic equation with biharmonic operator and multiple critical nonlinearities

abstract: We study the existence and non-existence of nontrivial weak solution of $$ {\Delta^2u-\mu\frac{u}{|x|^{4}} = \frac{|u|^{q_{\beta}-2}u}{|x|^{\beta}}+|u|^{q-2}u\quad\textrm{in ${\mathbb R}^N$,}} $$ where $N\geq 5$, $q_{\beta}=\frac{2(N-\beta)}{N-4}$, $0<\beta<4$, $1<q\leq 2^{**}$ and $\mu<\mu_1:=\big(\frac{N(N-4)}{4}\big)^2$. Using Pohozaev type of identity, we prove the non-existence result when $1<q< 2^{**}$. On the other hand when the equation has multiple critical nonlinearities i.e. $q=2^{**}$ and $-(N-2)^2\leq\mu<\mu_1$, we establish the existence of nontrivial solution using the Mountain-Pass theorem by Ambrosetti and Rabinowitz and the variational methods.


title: On plane waves in diluted relativistic cold plasmas

abstract: We briefly report on some exact results [G. Fiore, arXiv:1312.4665 preprint, to appear in J. Phys. A] regarding plane waves in a relativistic cold plasma. If the plasma, initially at rest, is reached by a transverse plane electromagnetic travelling-wave, then its motion has a very simple dependence on this wave in the limit of zero density, otherwise can be determined by an iterative procedure whose accuracy decreases with time or the plasma density. Thus one can describe in particular the impact of a very intense and short laser pulse onto a plasma and determine conditions for the "slingshot effect" [G. Fiore, R. Fedele, U. De Angelis, arXiv:1309.1400 preprint] to occur. The motion in vacuum of a charged test particle subject to a wave of the same kind is also determined, for any initial velocity.


title: Reverses and variations of Heinz inequality

abstract: Let $A, B$ be positive definite $n\times n$ matrices. We present several reverse Heinz type inequalities, in particular \begin{align*} \|AX+XB\|_2^2+ 2(\nu-1) \|AX-XB\|_2^2\leq \|A^{\nu}XB^{1-\nu}+A^{1-\nu}XB^{\nu}\|_2^2, \end{align*} where $X$ is an arbitrary $n \times n$ matrix, $\|\cdot\|_2$ is Hilbert-Schmidt norm and $\nu>1$. We also establish a Heinz type inequality involving the Hadamard product of the form \begin{align*} 2|||A^{1\over2}\circ B^{1\over2}|||\leq|||A^{s}\circ B^{1-t}+A^{1-s}\circ B^{t}||| \leq\max\{|||(A+B)\circ I|||,|||(A\circ B)+I|||\}, \end{align*} in which $s, t\in [0,1]$ and $|||\cdot|||$ is a unitarily invariant norm.


title: Surface vs bulk electronic structures of a moderately correlated topological insulator YbB6 revealed by ARPES

abstract: Systematic angle-resolved photoemission spectroscopy (ARPES) experiments have been carried out to investigate the bulk and (100) surface electronic structures of a topological mixed-valence insulator candidate, YbB6. The bulk states of YbB6 were probed with bulk-sensitive soft X-ray ARPES, which show strong three-dimensionality as required by cubic symmetry. Surprisingly the measured Yb 4f states are located around 1 eV and 2.3 eV below the Fermi level (EF), instead of being near EF as indicated by first principle calculations. The dispersive bands near EF are B 2p states, which hybridize with the 4f states. Using surface-sensitive vacuum-ultraviolet ARPES, we revealed two-dimensional surface states which form three electron-like Fermi surfaces (FSs) with Dirac-cone-like dispersions. The odd number of surface FSs gives the first indication that YbB6 is a moderately correlated topological insulator. The spin-resolved ARPES measurements provide further evidence that these surface states are spin polarized with spin locked to the crystal momentum. We have observed a second set of surface states with different topology (hole-like pockets). Clear folding of the bands suggest their origin from a 1X2 reconstructed surface. The topological property of the reconstructed surface states needs further studies.


title: Proceedings of the 6th European Conference on Python in Science (EuroSciPy 2013)

abstract: These are the proceedings of the 6th European Conference on Python in Science, EuroSciPy 2013, that was held in Brussels (21-25 August 2013).


title: Asymptotics of sharp constants of Markov-Bernstein inequalities in integral norm with Jacobi weight

abstract: The classical A. Markov inequality establishes a relation between the maximum modulus or the $L^{\infty}\left([-1,1]\right)$ norm of a polynomial $Q_{n}$ and of its derivative: $\|Q'_{n}\|\leqslant M_{n} n^{2}\|Q_{n}\|$, where the constant $M_{n}=1$ is sharp. The limiting behavior of the sharp constants $M_{n}$ for this inequality, considered in the space $L^{2}\left([-1,1], w^{(\alpha,\beta)}\right)$ with respect to the classical Jacobi weight $w^{(\alpha,\beta)}(x):=(1-x)^{\alpha}(x+1)^{\beta}$, is studied. We prove that, under the condition $|\alpha - \beta| < 4 $, the limit is $\lim_{n \to \infty} M_{n} = 1/(2 j_{\nu})$ where $j_{\nu}$ is the smallest zero of the Bessel function $J_{\nu}(x)$ and $2 \nu= \mbox{min}(\alpha, \beta) - 1$.


title: Adding SALT to Coupled Microcavities: the making of active photonic molecule lasers

abstract: A large body of work has accumulated over the years in the study of the optical properties of single and coupled microcavities for a variety of applications, ranging from filters to sensors and lasers. The focus has been mostly on the geometry of individual resonators and/or on their combination in arrangements often referred to as photonic molecules (PMs). Our primary concern will be the lasing properties of PMs as ideal candidates for the fabrication of integrated microlasers, photonic molecule lasers. Whereas most calculations on PM lasers have been based on cold-cavity (passive) modes, i.e. quasi-bound states, a recently formulated steady-state ab initio laser theory (SALT) offers the possibility to take into account the spectral properties of the underlying gain transition, its position and linewidth, as well as incorporating an arbitrary pump profile. We will combine two theoretical approaches to characterize the lasing properties of PM lasers: for two-dimensional systems, the generalized Lorenz-Mie theory will obtain the resonant modes of the coupled molecules in an active medium described by SALT. Not only is then the theoretical description more complete, the use of an active medium provides additional parameters to control, engineer and harness the lasing properties of PM lasers for ultra-low threshold and directional single-mode emission.


title: The Laplacian spectral excess theorem for distance-regular graphs

abstract: The spectral excess theorem states that, in a regular graph G, the average excess, which is the mean of the numbers of vertices at maximum distance from a vertex, is bounded above by the spectral excess (a number that is computed by using the adjacency spectrum of G), and G is distance-regular if and only if equality holds. In this note we prove the corresponding result by using the Laplacian spectrum without requiring regularity of G.


title: Un algorithme de test pour la connexit\'e temporelle des graphes dynamiques de faible densit\'e

abstract: We address the problem of testing whether a dynamic graph is temporally connected, i.e. a temporal path ({\em journey}) exists between all pairs of vertices. We consider a discrete version of the problem, where the topology is given as an evolving graph $\G=\{G_1,G_2,...,G_{k}\}$ in which only the set of (directed) edges varies. Two cases are studied, depending on whether a single edge or an unlimited number of edges can be crossed in a same $G_i$ (strict journeys {\it vs} non-strict journeys). For strict journeys, two existing algorithms designed for other problems can be adapted. However, we show that a dedicated approach achieves a better time complexity than one of these two algorithms in all cases, and than the other one for those graphs whose density is low at any time (though arbitrary over time). The time complexity of our algorithm is $O(k\mu n)$, where $k=|\G|$ is the number of time steps and $\mu=max(|E_i|)$ is the maximum {\em instant} density, to be contrasted with $m=|\cup E_i|$, the {\em cumulated} density. Indeed, it is not uncommon for a mobility scenario to satisfy, for instance, both $\mu=o(n)$ and $m=\Theta(n^2)$. We characterize the key values of $k, \mu$ and $m$ for which our algorithm should be used. For non-strict journeys, for which no algorithm is known, we show that a similar strategy can be used to answer the question, still in $O(k\mu n)$ time.


title: On the rate of convergence to equilibrium for the homogeneous Landau equation with soft potentials

abstract: We investigate in this work the rate of convergence to equilibrium of solutions to the spatially homogeneous Landau equation with soft potentials. Firstly, we prove a polynomial in time convergence using an entropy method with some new a priori estimates. Finally we prove an exponential in time convergence towards the equilibrium with the optimal rate given by the spectral gap of the associated linearized operator, combining new decay estimates for the semigroup generated by the linearized Landau operator in weighted $L^p$-spaces togheter with the polynomial decay described above.


title: Quenching the Anisotropic Heisenberg Chain: Exact Solution and Generalized Gibbs Ensemble Predictions

abstract: We study quenches in integrable spin-1/2 chains in which we evolve the ground state of the antiferromagnetic Ising model with the anisotropic Heisenberg Hamiltonian. For this nontrivially interacting situation, an application of the first-principles-based quench action method allows us to give an exact description of the postquench steady state in the thermodynamic limit. We show that a generalized Gibbs ensemble, implemented using all known local conserved charges, fails to reproduce the exact quench action steady state and to correctly predict postquench equilibrium expectation values of physical observables. This is supported by numerical linked-cluster calculations within the diagonal ensemble in the thermodynamic limit.


title: What sodium absorption lines tell us about type Ia supernovae

abstract: We propose that the sodium responsible for the variable NaI D absorption lines in some type Ia supernovae (SN Ia) originate mainly from dust residing at ~1 pc from the supernovae. In this Na-from-dust absorption (NaDA) model the process by which the SN Ia peak luminosity releases sodium from dust at 1 pc from the SN is similar to the processes by which solar radiation releases sodium from cometary dust when comets approach a distance of 1 AU from the Sun. The dust grains are not sublimated but rather stay intact, and release sodium by photon-stimulated desorption (PSD; or photo-sputtering). Some of the Na might start in the gas phase before the explosion. Weakening in absorption strength is caused by Na-ionizing radiation of the SN. We apply the NaDA model to SN 2006X and SN 2007le, and find it to comply better with the observed time variability of the NaI D absorption lines than the Na recombination model. The mass in the dusty shell of the NaDA model is much too high to be accounted for in the single-degenerate scenario for SN Ia. Therefore, the presence of variable NaI D lines in some SN Ia further weakens the already very problematic single-degenerate scenario for SN Ia.


title: VSCAN: An Enhanced Video Summarization using Density-based Spatial Clustering

abstract: In this paper, we present VSCAN, a novel approach for generating static video summaries. This approach is based on a modified DBSCAN clustering algorithm to summarize the video content utilizing both color and texture features of the video frames. The paper also introduces an enhanced evaluation method that depends on color and texture features. Video Summaries generated by VSCAN are compared with summaries generated by other approaches found in the literature and those created by users. Experimental results indicate that the video summaries generated by VSCAN have a higher quality than those generated by other approaches.


title: Super persistent transient in a master-slave configuration with Colpitts oscillators

abstract: Master-slave systems have been intensively investigated to model the application of chaos to communications. We considered Colpitts oscillators coupled according to a master-slave configuration to study chaos synchronisation. We revealed the existence of super persistent transients in this coupled system. Moreover, we showed that an additive noise added to the slave system may suppress chaos synchronisation. When synchronisation is not suppressed, the noise induces longer transients.


title: Cohomology of presheaves with oriented weak transfers

abstract: Over a field of characteristic zero, we establish the homotopy invariance of the Nisnevich cohomology of homotopy invariant presheaves with oriented weak transfers, and the agreement of Zariski and Nisnevich cohomology for such presheaves. This generalizes a foundational result in Voevodsky's theory of motives. The main idea is to find explicit smooth representatives of the correspondences which provide the input for Voevodsky's cohomological architecture.


title: Direct Bounds on Electroweak Scale Pseudo-Dirac Neutrinos from $\sqrt s=8$ TeV LHC Data

abstract: Seesaw models with a small lepton number breaking can naturally accommodate electroweak-scale pseudo-Dirac neutrinos with a sizable mixing with the active neutrinos, while satisfying the light neutrino oscillation data. Due to the smallness of the lepton number breaking parameter, the 'smoking gun' collider signature of same-sign dileptons is suppressed, and the heavy neutrinos in these models would manifest at the LHC dominantly through lepton number conserving trilepton final states. Using the recent CMS results for anomalous production of multilepton events at $\sqrt{s}$=8 TeV LHC with an integrated luminosity of $19.5~{\rm fb}^{-1}$, we derive direct upper bounds on the light-heavy neutrino mixing parameter as a function of the heavy Dirac neutrino mass. These limits extend the collider sensitivity to higher heavy neutrino masses up to about 500 GeV. In the lower mass range, our limits are comparable to the existing indirect constraints derived from Higgs and electroweak precision data.


title: Jet quenching effects on the anisotropic flow at RHIC

abstract: In this article we investigate how the energy and momentum deposited by partonic dijets in the quark-gluon plasma may perturb the geometry-induced hydrodynamic expansion of the bulk nuclear matter created in heavy ion collisions at RHIC. The coupling between the jets and the medium is done through a source term in the energy-momentum conservation equations for ideal hydrodynamics. We concentrate our attention at mid-rapidity and solve the equations event-by-event imposing boost-invariance. For $p_T \gtrsim 1$ GeV the anisotropic flow is found to be considerably enhanced, if the dijets deposit on average more than 12 GeV in the medium (or equivalently 6 GeV for each jet of the pair), which corresponds, in our model, to an average suppression greater than $65\%$ of the initial jet transverse energy.


title: New rigorous perturbation bounds for the LU and QR factorizations

abstract: Combining the modified matrix-vector equation approach with the technique of Lyapunov majorant function and the Banach fixed point principle, we obtain new rigorous perturbation bounds for the LU and QR factorizations with normwise or componentwise perturbations in the given matrix, where the componentwise perturbations have the form of backward errors resulting from the standard factorization algorithms. Each of the new rigorous perturbation bounds is a rigorous version of the first-order perturbation bound derived by the matrix-vector equation approach in the literature, and we present their explicit expressions. These bounds improve the results given by Chang and Stehl\'{e} [SIAM Journal on Matrix Analysis and Applications 2010; 31:2841--2859]. Moreover, we derive new tighter first-order perturbation bounds including two optimal ones for the LU factorization, and provide the explicit expressions of the optimal first-order perturbation bounds for the LU and QR factorizations.


title: Quasi-Periodic Variations in X-ray Emission and Long-Term Radio Observations: Evidence for a Two-Component Jet in Sw J1644+57

abstract: The continued observations of Sw J1644+57 in X-ray and radio bands accumulated a rich data set to study the relativistic jet launched in this tidal disruption event. The X-ray light curve of Sw J1644+57 from 5-30 days presents two kinds of quasi-periodic variations: a 200 second quasi-periodic oscillation (QPO) and a 2.7-day quasi-periodic variation. The latter has been interpreted by a precessing jet launched near the Bardeen-Petterson radius of a warped disk. Here we suggest that the $\sim$ 200s QPO could be associated with a second, narrower jet sweeping the observer line-of-sight periodically, which is launched from a spinning black hole in the misaligned direction with respect to the black hole's angular momentum. In addition, we show that this two-component jet model can interpret the radio light curve of the event, especially the re-brightening feature starting $\sim 100$ days after the trigger. From the data we infer that inner jet may have a Lorentz factor of $\Gamma_{\rm j} \sim 5.5$ and a kinetic energy of $E_{\rm k,iso} \sim 3.0 \times 10^{52} {\rm erg}$, while the outer jet may have a Lorentz factor of $\Gamma_{\rm j} \sim 2.5$ and a kinetic energy of $E_{\rm k,iso} \sim 3.0 \times 10^{53} {\rm erg}$.


title: BSM Primary Effects

abstract: Using the predictive power of the effective field theory approach, we present a physical parametrization of the leading effects beyond the SM (BSM), that give us at present the best way to constrain heavy new-physics at low-energies. We show that other BSM effects are not independent from these ones, and we provide the explicit correlations. This information is useful to know where to primarily look for new physics in future experiments, and to know how this new physics is related to previous measurements, most importantly in electroweak-symmetry breaking processes or Higgs physics.


title: Ergodicity of Approximate MCMC Chains with Applications to Large Data Sets

abstract: In many modern applications, difficulty in evaluating the posterior density makes performing even a single MCMC step slow. This difficulty can be caused by intractable likelihood functions, but also appears for routine problems with large data sets. Many researchers have responded by running approximate versions of MCMC algorithms. In this note, we develop quantitative bounds for showing the ergodicity of these approximate samplers. We then use these bounds to study the bias-variance trade-off of approximate MCMC algorithms. We apply our results to simple versions of recently proposed algorithms, including a variant of the "austerity" framework of Korratikara et al.


title: Bulldozing of granular material

abstract: We investigate the bulldozing motion of a granular sandpile driven forwards by a vertical plate. The problem is set up in the laboratory by emplacing the pile on a table rotating underneath a stationary plate; the continual circulation of the bulldozed material allows the dynamics to be explored over relatively long times, and the variation of the velocity with radius permits one to explore the dependence on bulldozing speed within a single experiment. We measure the time-dependent surface shape of the dune for a range of rotation rates, initial volumes and radial positions, for four granular materials, ranging from glass spheres to irregularly shaped sand. The evolution of the dune can be separated into two phases: a rapid initial adjustment to a state of quasi-steady avalanching perpendicular to the blade, followed by a much slower phase of lateral spreading and radial migration. The quasi-steady avalanching sets up a well-defined perpendicular profile with a nearly constant slope. This profile can be scaled by the depth against the bulldozer to collapse data from different times, radial positions and experiments onto common master curves that are characteristic of the granular material and depend on the local Froude number. The lateral profile of the dune along the face of the bulldozer varies more gradually with radial position, and evolves by slow lateral spreading. The spreading is asymmetrical, with the inward progress of the dune eventually arrested and its bulk migrating to larger radii. A one-dimensional depth-averaged model recovers the nearly linear perpendicular profile of the dune, but does not capture the finer nonlinear details of the master curves. A two-dimensional version of the model leads to an advection-diffusion equation that reproduces the lateral spreading and radial migration.


title: Energy levels, radiative rates and electron impact excitation rates for transitions in Si II

abstract: Energies for the lowest 56 levels, belonging to the 3s$^2$3p, 3s3p$^2$, 3p$^3$, 3s$^2$3d, 3s3p3d, 3s$^2$4$\ell$ and 3s$^2$5$\ell$ configurations of Si II, are calculated using the {\sc grasp} (General-purpose Relativistic Atomic Structure Package) code. Analogous calculations have also been performed (for up to 175 levels) using the Flexible Atomic Code ({\sc fac}). Furthermore, radiative rates are calculated for all E1, E2, M1 and M2 transitions. Extensive comparisons are made with available theoretical and experimental energy levels, and the accuracy of the present results is assessed to be better than 0.1 Ryd. Similarly, the accuracy for radiative rates (and subsequently lifetimes) is estimated to be better than 20% for most of the (strong) transitions. Electron impact excitation collision strengths are also calculated, with the Dirac Atomic R-matrix Code ({\sc darc}), over a wide energy range up to 13 Ryd. Finally, to determine effective collision strengths, resonances are resolved in a fine energy mesh in the thresholds region. These collision strengths are averaged over a Maxwellian velocity distribution and results listed over a wide range of temperatures, up to 10$^{5.5}$ K. Our data are compared with earlier $R$-matrix calculations and differences noted, up to a factor of two, for several transitions. Although scope remains for improvement, the accuracy for our results of collision strengths and effective collision strengths is assessed to be about 20% for a majority of transitions.


title: Additive Schwarz preconditioner for the general finite volume element discretization of symmetric elliptic problems

abstract: A symmetric and a nonsymmetric variant of the additive Schwarz preconditioner are proposed for the solution of a nonsymmetric system of algebraic equations arising from a general finite volume element discretization of symmetric elliptic problems with large jumps in the entries of the coefficient matrices across subdomains. It is shown in the analysis, that the convergence of the preconditioned GMRES iteration with the proposed preconditioners, depends polylogarithmically on the mesh parameters, in other words, the convergence is only weakly dependent on the mesh parameters, and it is robust with respect to the jumps in the coefficients.


title: Characterizations of sets of finite perimeter using heat kernels in metric spaces

abstract: The overarching goal of this paper is to link the notion of sets of finite perimeter (a concept associated with $N^{1,1}$-spaces) and the theory of heat semigroups (a concept related to $N^{1,2}$-spaces) in the setting of metric measure spaces whose measure is doubling and supports a $1$-Poincar\'e inequality. We prove a characterization of sets of finite perimeter in terms of a short time behavior of the heat semigroup in such metric spaces. We also give a new characterization of ${\rm BV}$ functions in terms of a near-diagonal energy in this general setting.


title: Split-Quaternions and the Dirac Equation

abstract: We show that Dirac 4-spinors admit an entirely equivalent formulation in terms of 2-spinors defined over the split-quaternions. In this formalism, a Lorentz transformation is represented as a $2 \times 2$ unitary matrix over the split-quaternions. The corresponding Dirac equation is then derived in terms of these 2-spinors. In this framework the $SO(3,2; {\bf R})$ symmetry of the Lorentz invariant scalar $\overline{\psi}\psi$ is manifest.


title: Kummer Subspaces of Tensor Products of Cyclic Algebras

abstract: We discuss the Kummer subspaces of tensor products of cyclic algebras, focusing mainly on the case of cyclic algebras of degree 3. We present a family of maximal spaces in the general case, classify all the monomial spaces in the case of tensor products of cyclic algebras of degree 3 using graph theory, and provide an upper bound for the dimension in the generic tensor product of cyclic algebras of degree 3.


title: On Hardness of Jumbled Indexing

abstract: Jumbled indexing is the problem of indexing a text $T$ for queries that ask whether there is a substring of $T$ matching a pattern represented as a Parikh vector, i.e., the vector of frequency counts for each character. Jumbled indexing has garnered a lot of interest in the last four years. There is a naive algorithm that preprocesses all answers in $O(n^2|\Sigma|)$ time allowing quick queries afterwards, and there is another naive algorithm that requires no preprocessing but has $O(n\log|\Sigma|)$ query time. Despite a tremendous amount of effort there has been little improvement over these running times. In this paper we provide good reason for this. We show that, under a 3SUM-hardness assumption, jumbled indexing for alphabets of size $\omega(1)$ requires $\Omega(n^{2-\epsilon})$ preprocessing time or $\Omega(n^{1-\delta})$ query time for any $\epsilon,\delta>0$. In fact, under a stronger 3SUM-hardness assumption, for any constant alphabet size $r\ge 3$ there exist describable fixed constant $\epsilon_r$ and $\delta_r$ such that jumbled indexing requires $\Omega(n^{2-\epsilon_r})$ preprocessing time or $\Omega(n^{1-\delta_r})$ query time.


title: Towards a Modular Recommender System for Research Papers written in Albanian

abstract: In the recent years there has been an increase in scientific papers publications in Albania and its neighboring countries that have large communities of Albanian speaking researchers. Many of these papers are written in Albanian. It is a very time consuming task to find papers related to the researchers' work, because there is no concrete system that facilitates this process. In this paper we present the design of a modular intelligent search system for articles written in Albanian. The main part of it is the recommender module that facilitates searching by providing relevant articles to the users (in comparison with a given one). We used a cosine similarity based heuristics that differentiates the importance of term frequencies based on their location in the article. We did not notice big differences on the recommendation results when using different combinations of the importance factors of the keywords, title, abstract and body. We got similar results when using only the title and abstract in comparison with the other combinations. Because we got fairly good results in this initial approach, we believe that similar recommender systems for documents written in Albanian can be build also in contexts not related to scientific publishing.


title: Distinguishing perfect set properties in separable metrizable spaces

abstract: All spaces are assumed to be separable and metrizable. Our main result is that the statement "For every space $X$, every closed subset of $X$ has the perfect set property if and only if every analytic subset of $X$ has the perfect set property" is equivalent to $\mathfrak{b}>\omega_1$ (hence, in particular, it is independent of $\mathsf{ZFC}$). This, together with a theorem of Solecki and an example of Miller, will allow us to determine the status of the statement "For every space $X$, if every $\mathbf{\Gamma}$ subset of $X$ has the perfect set property then every $\mathbf{\Gamma}'$ subset of $X$ has the perfect set property" as $\mathbf{\Gamma},\mathbf{\Gamma}'$ range over all pointclasses of complexity at most analytic or coanalytic. Along the way, we define and investigate a property of independent interest. We will say that a subset $W$ of $2^\omega$ has the Grinzing property if it is uncountable and for every uncountable $Y\subseteq W$ there exists an uncountable collection consisting of uncountable subsets of $Y$ with pairwise disjoint closures in $2^\omega$. The following theorems hold. (1) There exists a subset of $2^\omega$ with the Grinzing property. (2) Assume $\mathsf{MA}+\neg\mathsf{CH}$. Then $2^\omega$ has the Grinzing property. (3) Assume $\mathsf{CH}$. Then $2^\omega$ does not have the Grinzing property. The first result was obtained by Miller using a theorem of Todor\v{c}evi\'c, and is needed in the proof of our main result.


title: Geo-neutrinos and Earth Models

abstract: We present the current status of geo-neutrino measurements and their implications for radiogenic heating in the mantle. Earth models predict different levels of radiogenic heating and, therefore, different geo-neutrino fluxes from the mantle. Seismic tomography reveals features in the deep mantle possibly correlated with radiogenic heating and causing spatial variations in the mantle geo-neutrino flux at the Earth surface. An ocean-based observatory offers the greatest sensitivity to the mantle flux and potential for resolving Earth models and mantle features. Refinements to estimates of the geo-neutrino flux from continental crust reduce uncertainty in measurements of the mantle flux, especially measurements from land-based observatories. These refinements enable the resolution of Earth models using the combined measurements from multiple continental observatories.


abstract: The complex planetary synchronization structure of the solar system, which since Pythagoras of Samos (ca. 570-495 BC) is known as the music of the spheres, is briefly reviewed from the Renaissance up to contemporary research. Copernicus' heliocentric model from 1543 suggested that the planets of our solar system form a kind of mutually ordered and quasi-synchronized system. From 1596 to 1619 Kepler formulated preliminary mathematical relations of approximate commensurabilities among the planets, which were later reformulated in the Titius-Bode rule (1766-1772) that successfully predicted the orbital position of Ceres and Uranus. Following the discovery of the ~11 yr sunspot cycle, in 1859 Wolf suggested that the observed solar variability could be approximately synchronized with the orbital movements of Venus, Earth, Jupiter and Saturn. Modern research have further confirmed that: (1) the planetary orbital periods can be approximately deduced from a simple system of resonant frequencies; (2) the solar system oscillates with a specific set of gravitational frequencies, and many of them (e.g. within the range between 3 yr and 100 yr) can be approximately constructed as harmonics of a base period of ~178.38 yr; (3) solar and climate records are also characterized by planetary harmonics from the monthly to the millennia time scales. This short review concludes with an emphasis on the contribution of the author's research on the empirical evidences and physical modeling of both solar and climate variability based on astronomical harmonics. The general conclusion is that the solar system works as a resonator characterized by a specific harmonic planetary structure that synchronizes also the Sun's activity and the Earth's climate.


title: SIESTA-PEXSI: Massively parallel method for efficient and accurate \textit{ab initio} materials simulation without matrix diagonalization

abstract: We describe a scheme for efficient large-scale electronic-structure calculations based on the combination of the pole expansion and selected inversion (PEXSI) technique with the SIESTA method, which uses numerical atomic orbitals within the Kohn-Sham density functional theory (KSDFT) framework. The PEXSI technique can efficiently utilize the sparsity pattern of the Hamiltonian and overlap matrices generated in SIESTA, and for large systems has a much lower computational complexity than that associated with the matrix diagonalization procedure. The PEXSI technique can be used to evaluate the electron density, free energy, atomic forces, density of states and local density of states without computing any eigenvalue or eigenvector of the Kohn-Sham Hamiltonian. It can achieve accuracy fully comparable to that obtained from a matrix diagonalization procedure for general systems, including metallic systems at low temperature. The PEXSI method is also highly scalable. With the recently developed massively parallel PEXSI technique, we can make efficient use of more than $10,000$ processors on high performance machines. We demonstrate the performance and accuracy of the SIESTA-PEXSI method using several examples of large scale electronic structure calculations, including 1D, 2D and bulk problems with insulating, semi-metallic, and metallic character.


title: Analytical approach to the dynamics of facilitated spin models on random networks

abstract: Facilitated spin models were introduced some decades ago to mimic systems characterized by a glass transition. Recent developments have shown that a class of facilitated spin models is also able to reproduce characteristic signatures of the structural relaxation properties of glass-forming liquids. While the equilibrium phase diagram of these models can be calculated analytically, the dynamics are usually investigated numerically. Here we propose a new network-based approach, called approximate master equation (AME), to the dynamics of the Fredrickson-Andersen model. The approach correctly predicts the critical temperature at which the glass transition occurs. We also find excellent agreement between the theory and the numerical simulations for the transient regime, except in close proximity of the liquid-glass transition. Finally, we analytically characterize the critical clusters of the model and show that the departures between our AME approach and the Monte Carlo can be related to the large interface between frozen and unfrozen spins at temperatures close to the glass transition.


title: Determination of flow stress and the critical strain for the onset of dynamic recrystallization using a sine function

abstract: A new model has been developed to estimate the flow stress under hot deformation conditions up to the peak of the stress-strain curves. The model is constructed on the basis of the general form of sine functions by introducing an additional exponent. Besides, an equation is derived from the model, which could be used to predict the critical strain for the onset of dynamic recrystallization. The model has been examined for a 304 austenitic stainless steel. The estimated flow stress had an average percentage error of %1.66.


title: Chaotic Inflation in No-Scale Supergravity with String Inspired Moduli Stabilization

abstract: The simple chaotic inflation is highly consistent with the BICEP2 experiment, and no-scale supergravity can be realized naturally in various string compactifications. Thus, we construct a chaotic inflation model in no-scale supergravity inspired from Type IIB string compactification with an anomalous U(1)_X gauged symmetry. We introduce two moduli T_1 and T_2 which transform non-trivially under U(1)_X, and some pairs of fundamental quarks charged under the SU(N) x U(1)_X gauge group. The non-trivial transformations of moduli under U(1)_X lead to a moduli-dependent Fayet-Iliopoulos (FI) term. The modulus T_2 and the real component of T_1 are stabilized by the non-perturbative effect from quark condensation and the U(1)_X D-term. In particular, the stabilization from the anomalous U(1)_X D-term with moduli-dependent FI term is crucial for inflation since it gives heavy mass to the real component of the modulus T_1 while keeping its axionic part light. Choosing the proper parameters, we obtain a global Minkowski vacuum where the imaginary part of T_1 has a quadratic potential for chaotic inflation.


title: No Superluminal Signaling Implies Unconditionally Secure Bit Commitment

abstract: Bit commitment (BC) is an important cryptographic primitive for an agent to convince a mutually mistrustful party that she has already made a binding choice of 0 or 1 but only to reveal her choice at a later time. Ideally, a BC protocol should be simple, reliable, easy to implement using existing technologies, and most importantly unconditionally secure in the sense that its security is based on an information-theoretic proof rather than computational complexity assumption or the existence of a trustworthy arbitrator. Here we report such a provably secure scheme involving only one-way classical communications whose unconditional security is based on no superluminal signaling (NSS). Our scheme is inspired by the earlier works by Kent, who proposed two impractical relativistic protocols whose unconditional securities are yet to be established as well as several provably unconditionally secure protocols which rely on both quantum mechanics and NSS. Our scheme is conceptually simple and shows for the first time that quantum communication is not needed to achieve unconditional security for BC. Moreover, with purely classical communications, our scheme is practical and easy to implement with existing telecom technologies. This completes the cycle of study of unconditionally secure bit commitment based on known physical laws.


title: Liquid Phase 3D Printing for Quickly Manufacturing Metal Objects with Low Melting Point Alloy Ink

abstract: Conventional 3D printings are generally time-consuming and printable metal inks are rather limited. From an alternative way, we proposed a liquid phase 3D printing for quickly making metal objects. Through introducing metal alloys whose melting point is slightly above room temperature as printing inks, several representative structures spanning from one, two and three dimension to more complex patterns were demonstrated to be quickly fabricated. Compared with the air cooling in a conventional 3D printing, the liquid-phase-manufacturing offers a much higher cooling rate and thus significantly improves the speed in fabricating metal objects. This unique strategy also efficiently prevents the liquid metal inks from air oxidation which is hard to avoid otherwise in an ordinary 3D printing. Several key physical factors (like properties of the cooling fluid, injection speed and needle diameter, types and properties of the printing ink, etc.) were disclosed which would evidently affect the printing quality. In addition, a basic route to make future liquid phase 3D printer incorporated with both syringe pump and needle arrays was also suggested. The liquid phase 3D printing method, which owns potential values not available in a conventional modality, opens an efficient way for quickly making metal objects in the coming time.


title: Black hole radiation of spin-1 particles in (1+2) dimensions

abstract: The radiation of vector particles by black holes in (1+2) dimensions is investigated within the WKB approximation. We consider the process of quantum tunnelling of bosons through an event horizon of the black hole. The emission temperature for the Schwarzschild background geometry coincides with the Hawking temperature and for the Rindler spacetime the temperature is the Unruh temperature. We also obtain the radiation temperatures for the de Sitter spacetime.


title: Analysis of Digital Knapsack Based Sealed Bid Auction

abstract: The need of totally secure online auction has led to the invention of many auction protocols. But as new attacks are developed, auction protocols also require corresponding strengthening. We analyze the auction protocol based on the well-known mathematical public-key knapsack problem for the design of asymmetric public-key knapsack trapdoor cryptosystem. Even though the knapsack system is not cryptographically secure, it can be used in certain auction situations. We describe the limitations of the protocol like detecting and solving the tie between bidders, malicious behavior of participants and also selection of price set by the seller and offer solutions.