title: The JCMT Gould Belt Survey: A First Look at Dense Cores in Orion B

abstract: We present a first look at the SCUBA-2 observations of three sub-regions of the Orion B molecular cloud: LDN 1622, NGC 2023/2024, and NGC 2068/2071, from the JCMT Gould Belt Legacy Survey. We identify 29, 564, and 322 dense cores in L1622, NGC 2023/2024, and NGC 2068/2071 respectively, using the SCUBA-2 850 micron map, and present their basic properties, including their peak fluxes, total fluxes, and sizes, and an estimate of the corresponding 450 micron peak fluxes and total fluxes, using the FellWalker source extraction algorithm. Assuming a constant temperature of 20 K, the starless dense cores have a mass function similar to that found in previous dense core analyses, with a Salpeter-like slope at the high-mass end. The majority of cores appear stable to gravitational collapse when considering only thermal pressure; indeed, most of the cores which have masses above the thermal Jeans mass are already associated with at least one protostar. At higher cloud column densities, above 1-2 x 10^23 cm^-2, most of the mass is found within dense cores, while at lower cloud column densities, below 1 x 10^23 cm^-2, this fraction drops to 10% or lower. Overall, the fraction of dense cores associated with a protostar is quite small (<8%), but becomes larger for the densest and most centrally concentrated cores. NGC 2023 / 2024 and NGC 2068/2071 appear to be on the path to forming a significant number of stars in the future, while L1622 has little additional mass in dense cores to form many new stars.


title: Untangling Performance from Success

abstract: Fame, popularity and celebrity status, frequently used tokens of success, are often loosely related to, or even divorced from professional performance. This dichotomy is partly rooted in the difficulty to distinguish performance, an individual measure that captures the actions of a performer, from success, a collective measure that captures a community's reactions to these actions. Yet, finding the relationship between the two measures is essential for all areas that aim to objectively reward excellence, from science to business. Here we quantify the relationship between performance and success by focusing on tennis, an individual sport where the two quantities can be independently measured. We show that a predictive model, relying only on a tennis player's performance in tournaments, can accurately predict an athlete's popularity, both during a player's active years and after retirement. Hence the model establishes a direct link between performance and momentary popularity. The agreement between the performance-driven and observed popularity suggests that in most areas of human achievement exceptional visibility may be rooted in detectable performance measures.


title: Non-Convex Distributed Optimization

abstract: We study distributed non-convex optimization on a time-varying multi-agent network. Each node has access to its own smooth local cost function, and the collective goal is to minimize the sum of these functions. We generalize the results obtained previously to the case of non-convex functions. Under some additional technical assumptions on the gradients we prove the convergence of the distributed push-sum algorithm to some critical point of the objective function. By utilizing perturbations on the update process, we show the almost sure convergence of the perturbed dynamics to a local minimum of the global objective function. Our analysis shows that this noised procedure converges at a rate of $O(1/t)$.


title: Sums of Quadratic residues and nonresidues

abstract: It is well known that when a prime $p$ is congruent to 1 modulo 4, the sum of the quadratic residues equals the sum of the quadratic nonresidues. In this note we give analogous results for the case where $p$ is congruent to 3 modulo 4.


title: Pattern completion in symmetric threshold-linear networks

abstract: Threshold-linear networks are a common class of firing rate models that describe recurrent interactions among neurons. Unlike their linear counterparts, these networks generically possess multiple stable fixed points (steady states), making them viable candidates for memory encoding and retrieval. In this work, we characterize stable fixed points of general threshold-linear networks with constant external drive, and discover constraints on the co-existence of fixed points involving different subsets of active neurons. In the case of symmetric networks, we prove the following antichain property: if a set of neurons $\tau$ is the support of a stable fixed point, then no proper subset or superset of $\tau$ can support a stable fixed point. Symmetric threshold-linear networks thus appear to be well suited for pattern completion, since the dynamics are guaranteed not to get "stuck" in a subset or superset of a stored pattern. We also show that for any graph G, we can construct a network whose stable fixed points correspond precisely to the maximal cliques of G. As an application, we design network decoders for place field codes, and demonstrate their efficacy for error correction and pattern completion. The proofs of our main results build on the theory of permitted sets in threshold-linear networks, including recently-developed connections to classical distance geometry.


title: Very weak solutions of the Stokes problem in a convex polygon

abstract: Motivated by the study of the corner singularities in the so-called cavity flow, we establish in this article, the existence and uniqueness of solutions in $L^2(\Omega)^2$ for the Stokes problem in a domain $\Omega,$ when $\Omega$ is a smooth domain or a convex polygon. We establish also a trace theorem and show that the trace of $u$ can be arbitrary in $L^2(\partial\Omega)^2.$ The results are also extended to the linear evolution Stokes problem.


title: Estimating Learning Effects: A Short-Time Fourier Transform Regression Model for MEG Source Localization

abstract: Magnetoencephalography (MEG) has a high temporal resolution well-suited for studying perceptual learning. However, to identify where learning happens in the brain, one needs to ap- ply source localization techniques to project MEG sensor data into brain space. Previous source localization methods, such as the short-time Fourier transform (STFT) method by Gramfort et al.([Gramfort et al., 2013]) produced intriguing results, but they were not designed to incor- porate trial-by-trial learning effects. Here we modify the approach in [Gramfort et al., 2013] to produce an STFT-based source localization method (STFT-R) that includes an additional regression of the STFT components on covariates such as the behavioral learning curve. We also exploit a hierarchical L 21 penalty to induce structured sparsity of STFT components and to emphasize signals from regions of interest (ROIs) that are selected according to prior knowl- edge. In reconstructing the ROI source signals from simulated data, STFT-R achieved smaller errors than a two-step method using the popular minimum-norm estimate (MNE), and in a real-world human learning experiment, STFT-R yielded more interpretable results about what time-frequency components of the ROI signals were correlated with learning.


title: Strongly interacting blow up bubbles for the mass critical NLS

abstract: We construct a new class of multi-solitary wave solutions for the mass critical two dimensional nonlinear Schrodinger equation (NLS). Given any integer K>1, there exists a global (for positive time) solution of (NLS) that decomposes asymptotically into a sum of solitary waves centered at the vertices of a K-sided regular polygon and concentrating at a logarithmic rate in large time. This solution blows up in infinite time with logarithmic rate. Using the pseudo-conformal transform, this yields the first example of solution blowing up in finite time with a rate strictly above the pseudo-conformal one. Such solution concentrates K bubbles at a point. These special behaviors are due to strong interactions between the waves, in contrast with previous works on multi-solitary waves of (NLS) where interactions do not affect the blow up rate.


title: Compressive hyperspectral imaging via adaptive sampling and dictionary learning

abstract: In this paper, we propose a new sampling strategy for hyperspectral signals that is based on dictionary learning and singular value decomposition (SVD). Specifically, we first learn a sparsifying dictionary from training spectral data using dictionary learning. We then perform an SVD on the dictionary and use the first few left singular vectors as the rows of the measurement matrix to obtain the compressive measurements for reconstruction. The proposed method provides significant improvement over the conventional compressive sensing approaches. The reconstruction performance is further improved by reconditioning the sensing matrix using matrix balancing. We also demonstrate that the combination of dictionary learning and SVD is robust by applying them to different datasets.


title: Muon g-2 through a flavor structure on soft SUSY terms

abstract: In this work we analyze the possibility to explain the muon anomalous magnetic moment discrepancy within theory and experiment through lepton flavor violation processes. We propose a flavor extended MSSM by considering a hierarchical family structure for the trilinear scalar Soft-Supersymmetric terms of the Lagranagian, present at the SUSY breaking scale. We obtain analytical results for the rotation mass matrix, with the consequence of having non-universal slepton masses and the possibility of leptonic flavour mixing. The one-loop supersymmetric contributions to the leptonic flavour violating process $\tau \to \mu\gamma$ are calculated in the physical basis, with slepton flavour mixed states, instead of using the well known Mass Insertion Method. We present the regions in parameter space where the muon g-2 problem is either entirely solved or partially reduced through the contribution of these flavor violating processes.


title: Existence of self-accelerating fronts for a non-local reaction-diffusion equations

abstract: We describe the accelerated propagation wave arising from a non-local reaction-diffusion equation. This equation originates from an ecological problem, where accelerated biological invasions have been documented. The analysis is based on the comparison of this model with a related local equation, and on the analysis of the dynamics of the solutions of this second model thanks to probabilistic methods.


title: Velocity resolved [CII] spectroscopy of the center and the BCLMP302 region of M33 (HerM33es)

abstract: We aim to understand the contribution of the ionized, atomic and molecular phases of the ISM to the [CII] emission from clouds near the dynamical center and the BCLMP302 HII region in the north of the nearby galaxy M33 at a spatial resolution of 50pc. We combine high resolution [CII] spectra taken with the HIFI spectrometer onboard the Herschel satellite with [CII] Herschel-PACS maps and ground-based observations of CO(2-1) and HI. All data are at a common spatial resolution of 50pc. Typically, the [CII] lines have widths intermediate between the narrower CO(2-1) and broader HI line profiles. We decomposed the [CII] spectra in terms of contribution from molecular and atomic gas detected in CO(2-1) and HI, respectively. We find that the relative contribution of molecular and atomic gas traced by CO(2-1) and HI varies depends mostly on the local physical conditions and geometry. We estimate that 11-60% and 5-34% of the [CII] intensities in the center and in BCLMP302, respectively, arise at velocities showing no CO(2-1) or HI emission and could arise in CO-dark molecular gas. The deduced strong variation in the [CII] emission not associated with CO and HI cannot be explained in terms of differences in A_v, far-ultraviolet radiation field, and metallicity between the two studied regions. Hence the relative amounts of diffuse (CO-dark) and dense molecular gas possibly vary on spatial scales smaller than 50pc. Based on the emission measure observed at radio wavelengths we estimate the contribution of ionized gas at a few positions to lie between 10-25%. The correlations between the intensities of tracers corresponding to the same velocity range as [CII], differ from the correlation derived from PACS data. The results in this paper emphasize the need for velocity-resolved observations to discern the contribution of different components of the ISM to [CII] emission. (abridged)


title: Shape-constrained uncertainty quantification in unfolding steeply falling elementary particle spectra

abstract: The high energy physics unfolding problem is an important statistical inverse problem in data analysis at the Large Hadron Collider (LHC) at CERN. The goal of unfolding is to make nonparametric inferences about a particle spectrum from measurements smeared by the finite resolution of the particle detectors. Previous unfolding methods use ad hoc discretization and regularization, resulting in confidence intervals that can have significantly lower coverage than their nominal level. Instead of regularizing using a roughness penalty or stopping iterative methods early, we impose physically motivated shape constraints: positivity, monotonicity, and convexity. We quantify the uncertainty by constructing a nonparametric confidence set for the true spectrum, consisting of all those spectra that satisfy the shape constraints and that predict the observations within an appropriately calibrated level of fit. Projecting that set produces simultaneous confidence intervals for all functionals of the spectrum, including averages within bins. The confidence intervals have guaranteed conservative frequentist finite-sample coverage in the important and challenging class of unfolding problems for steeply falling particle spectra. We demonstrate the method using simulations that mimic unfolding the inclusive jet transverse momentum spectrum at the LHC. The shape-constrained intervals provide usefully tight conservative inferences, while the conventional methods suffer from severe undercoverage.


title: Butcher series: A story of rooted trees and numerical methods for evolution equations

abstract: Butcher series appear when Runge-Kutta methods for ordinary differential equations are expanded in power series of the step size parameter. Each term in a Butcher series consists of a weighted elementary differential, and the set of all such differentials is isomorphic to the set of rooted trees, as noted by Cayley in the mid 19th century. A century later Butcher discovered that rooted trees can also be used to obtain the order conditions of Runge-Kutta methods, and he found a natural group structure, today known as the Butcher group. It is now known that many numerical methods also can be expanded in Butcher series; these are called B-series methods. A long-standing problem has been to characterize, in terms of qualitative features, all B-series methods. Here we tell the story of Butcher series, stretching from the early work of Cayley, to modern developments and connections to abstract algebra, and finally to the resolution of the characterization problem. This resolution introduces geometric tools and perspectives to an area traditionally explored using analysis and combinatorics.


title: Innovation Pursuit: A New Approach to Subspace Clustering

abstract: In subspace clustering, a group of data points belonging to a union of subspaces are assigned membership to their respective subspaces. This paper presents a new approach dubbed Innovation Pursuit (iPursuit) to the problem of subspace clustering using a new geometrical idea whereby subspaces are identified based on their relative novelties. We present two frameworks in which the idea of innovation pursuit is used to distinguish the subspaces. Underlying the first framework is an iterative method that finds the subspaces consecutively by solving a series of simple linear optimization problems, each searching for a direction of innovation in the span of the data potentially orthogonal to all subspaces except for the one to be identified in one step of the algorithm. A detailed mathematical analysis is provided establishing sufficient conditions for iPursuit to correctly cluster the data. The proposed approach can provably yield exact clustering even when the subspaces have significant intersections. It is shown that the complexity of the iterative approach scales only linearly in the number of data points and subspaces, and quadratically in the dimension of the subspaces. The second framework integrates iPursuit with spectral clustering to yield a new variant of spectral-clustering-based algorithms. The numerical simulations with both real and synthetic data demonstrate that iPursuit can often outperform the state-of-the-art subspace clustering algorithms, more so for subspaces with significant intersections, and that it significantly improves the state-of-the-art result for subspace-segmentation-based face clustering.


title: Lattice Surfaces and smallest triangles

abstract: We calculate the area of the smallest triangle and the area of the smallest virtual triangle for many known lattice surfaces. We show that our list of the lattice surfaces for which the area of the smallest virtual triangle greater than .05 is complete. In particular, this means that there are no new lattice surfaces for which the area of the smallest virtual triangle is greater than .05. Our method follows an algorithm described by Smillie and Weiss and improves on it in certain respects.


title: Existence of solution to a nonlinear first-order dynamic equation on time scales

abstract: We prove existence of solution to a nonlinear first-order nabla dynamic equation on an arbitrary bounded time scale with boundary conditions, where the right-hand side of the dynamic equation is a continuous function.


title: The Atacama Cosmology Telescope: Dynamical masses for 44 SZ-selected galaxy clusters over 755 square degrees

abstract: We present galaxy velocity dispersions and dynamical mass estimates for 44 galaxy clusters selected via the Sunyaev-Zel'dovich (SZ) effect by the Atacama Cosmology Telescope. Dynamical masses for 18 clusters are reported here for the first time. Using \Nbody\ simulations, we model the different observing strategies used to measure the velocity dispersions and account for systematic effects resulting from these strategies. We find that the galaxy velocity distributions may be treated as isotropic, and that an aperture correction of up to 7 per cent in the velocity dispersion is required if the spectroscopic galaxy sample is sufficiently concentrated towards the cluster centre. Accounting for the radial profile of the velocity dispersion in simulations enables consistent dynamical mass estimates regardless of the observing strategy. Cluster masses $M_{200}$ are in the range $(1-15)\times10^{14}M_\odot$. Comparing with masses estimated from the SZ distortion assuming a gas pressure profile derived from X-ray observations gives a mean SZ-to-dynamical mass ratio of $1.10\pm0.13$, but there is an additional 0.14 systematic uncertainty due to the unknown velocity bias; the statistical uncertainty is dominated by the scatter in the mass-velocity dispersion scaling relation. This ratio is consistent with previous determinations at these mass scales.


title: Introduction of the Residue Number Arithmetic Logic Unit With Brief Computational Complexity Analysis

abstract: Digital System Research has pioneered the mathematics and design for a new class of computing machine using residue numbers. Unlike prior art, the new breakthrough provides methods and apparatus for general purpose computation using several new residue based fractional representations. The result is that fractional arithmetic may be performed without carry. Additionally, fractional operations such as addition, subtraction and multiplication of a fraction by an integer occur in a single clock period, regardless of word size. Fractional multiplication is of the order O(p), where p equals the number of residues. More significantly, complex operations, such as sum of products, may be performed in an extended format, where fractional products are performed and summed using single clock instructions, regardless of word width, and where a normalization operation with an execution time of the order O(p) is performed as a final step.


title: Pure point diffraction and Poisson summation

abstract: We prove that the diffraction formula for regular model sets is equivalent to the Poisson Summation Formula for the underlying lattice. This is achieved using Fourier analysis of unbounded measures on locally compact abelian groups as developed by Argabright and de Lamadrid. We also obtain diffraction results for certain classes of non-regular so-called weak model sets.


title: Formulae for flavour neutrino masses and its application to texture two zeros

abstract: We demonstrate the usefulness of flavour neutrino masses expressed in terms of $M_{ee},M_{e\mu}$ and $M_{e\tau}$. The analytical expressions for the flavour neutrino masses, mass eigenstates and physical CP-violating Majorana phases for texture two zeros are obtained exactly.


title: Dengue disease: a multiobjective viewpoint

abstract: During the last decades, the global prevalence of dengue progressed dramatically. It is a disease that is now endemic in more than one hundred countries of Africa, America, Asia, and the Western Pacific. In this paper, we present a mathematical model for the dengue disease transmission described by a system of ordinary differential equations and propose a multiobjective approach to find the most effective ways of controlling the disease. We use evolutionary multiobjective optimization (EMO) algorithms to solve the resulting optimization problem, providing the performance comparison of different algorithms. The obtained results show that the multiobjective approach is an effective tool to solve the problem, giving higher quality and wider range of solutions compared to the traditional technique. The obtained trade-offs provide a valuable information about the dynamics of infection transmissions and can be used as an input in the process of planning the intervention measures by the health authorities. Additionally, a suggested hybrid EMO algorithm produces highly superior performance compared to five other state-of-the-art EMO algorithms, being indispensable to efficiently optimize the proposed model.


title: Regularity of Solutions to a Parabolic Free Boundary Problem with Variable Coefficients

abstract: The main result of this paper is to prove that viscosity solutions to a parabolic free boundary problem with variable coefficients are Lipschitz continuous under the assumptions that the solution has a Lipschitz free boundary and satisfies a non-degeneracy condition. This is the optimal regularity for this problem.


title: Study on Green's function on TI surface

abstract: In theory of superconducting junctions, Green's function has an important role to obtain Andreev bound states, local density of states and Josephson current in a systematic way. In this article, we show how to construct Green's function on the surface of topological insulator following McMillan's formalism where the energy spectrum of electrons obeys a linear dispersion. For a model of superconductor (S)/ferromagnet (F)/normal metal (N) junction, we show that the generation of Majorana Fermion gives rise to the enhanced local density of states and pair amplitude of odd-frequency pairing. We also derive an extended Furusaki-Tsukada's formula of d.c. Josephson current in S/F/S junctions. The obtained Josephson current depends on the direction and magnitude of the magnetization.


title: Regularity of the Free Boundary for a Bernoulli-Type Parabolic Problem with Variable Coefficients

abstract: This paper studies the regularity of the free boundary for viscosity solutions to a parabolic Bernoulli-type free boundary problem with variable coefficients. The main result is that Lipschitz free boundaries are $C^1$ with a normal vector that varies with a H\"older modulus of continuity in both space and in time. As a consequence, the viscosity solution satisfies the free boundary problem in a classical sense.


title: Shifted moments of L functions and moments of theta functions

abstract: Assuming the Riemann Hypothesis, Soundararajan showed that $\displaystyle{\int_{0}^{T} \vert \zeta(1/2 + it)\vert^{2k} \ll T(\log T)^{k^2 + \epsilon}}$ . His method was used by Chandee to obtain upper bounds for shifted moments of the Riemann Zeta function. Building on these ideas, we obtain, conditionally, upper bounds for shifted moments of Dirichlet $L$- functions which allow us to derive upper bounds for moments of theta functions.


title: Magnetospheric ULF waves driven by external sources

abstract: The multi-spacecraft missions (Cluster and THEMIS) observations allowed to collect large data base for Ultra Low Frequency (ULF) waves properties, their localization, and sources. Mainly here we focused on these recent results. Studying of the source and characteristics of ULF waves can help in the understanding of the interaction and energy transport from solar wind to the magnetosphere. Here we present peculiarities of ULF waves generated by different solar wind phenomenon: surface magnetopause instability, magnetosphere cavity modes and solar wind dynamic pressure sudden impulses (SI) penetration into the magnetosphere. Permanent observations of ULF waves involve existence of the permanent source and, as the previous studies showed, the contributions to Pc4-Pc5 ULF wave power from the external sources are larger than the contribution from internal magnetosphere sources. The Kelvin-Helmholtz instability (KHI) can generate on the magnetosphere flanks classical ULF resonant waves with spatially localized amplitude maximum. As observations show the constraint satisfaction of KHI development is quite rare. SI in solar wind dynamic pressure generate ULF waves with different polarization and frequency close to the frequency of the local field line resonance (FLR). Wide range of temporal and amplitude characteristics of solar wind dynamics can generate magnetosphere cavity modes and magnetosonic perturbations which penetrate through the magnetosphere and can couple with the local FLR modes. The observed dependence of ULF waves properties on their localization are well correspond to these sources and their occurrence.


title: From Interstellar Clouds to Stars

abstract: I review (1) Physics of Star Formation & Open Questions; (2) Structure & Dynamics of Star-Forming Clouds & Young Clusters; (3) Star Formation Rates: Observations & Theoretical Implications.


title: Viscous Cosmology

abstract: We discuss the possibility to implement a viscous cosmological model, attributing to the dark matter component a behaviour described by bulk viscosity. Since bulk viscosity implies negative pressure, this rises the possibility to unify the dark sector. At the same time, the presence of dissipative effects may alleviate the so called small scale problems in the $\Lambda$CDM model. While the unified viscous description for the dark sector does not lead to consistent results, the non-linear behaviour indeed improves the situation with respect to the standard cosmological model.


title: Rate equation model of phototransduction into the membranous disks of mouse rod cells

abstract: A theoretical model was developed to investigate the rod phototransduction process in the mouse. In particular, we explored the biochemical reactions of several chemical components that contribute to the signaling process into/around the membranous disks in the outer segments of the rod cells. We constructed a rate equation model incorporating the molecular crowding effects of rhodopsin according to experimental results, which may hinder the diffusion of molecules on the disk mem- brane. The present model could effectively reproduce and explain the mechanisms of the following phenomena observed in experiments. First, the activations and relaxation of the wild-type mouse rod cell progressed more slowly than those of mutant cells containing half the amount of rhodopsin on the disk membrane. Second, the strong photoactivated state of the cell was sustained for a longer period when the light stimuli were strong. Finally, the lifetime of photoactivation exhibited a logarithmic increase with increasing light strength given exposure to strong light stimuli.


title: Do the Fermi Gamma-Ray Burst Monitor and Swift Burst Alert Telescope see the Same Short Gamma-Ray Bursts?

abstract: Compact binary system mergers are expected to generate gravitational radiation detectable by ground-based interferometers. A subset of these, the merger of a neutron star with another neutron star or a black hole, are also the most popular model for the production of short gamma-ray bursts (GRBs). The Swift Burst Alert Telescope (BAT) and the Fermi Gamma-ray Burst Monitor (GBM) trigger on short GRBs (SGRBs) at rates that reflect their relative sky exposures, with the BAT detecting 10 per year compared to about 45 for GBM. We examine the SGRB populations detected by Swift BAT and Fermi GBM. We find that the Swift BAT triggers on weaker SGRBs than Fermi GBM, providing they occur close to the center of the BAT field-of-view, and that the Fermi GBM SGRB detection threshold remains flatter across its field-of-view. Overall, these effects combine to give the instruments the same average sensitivity, and account for the SGRBs that trigger one instrument but not the other. We do not find any evidence that the BAT and GBM are detecting significantly different populations of SGRBs. Both instruments can detect untriggered SGRBs using ground searches seeded with time and position. The detection of SGRBs below the on-board triggering sensitivities of Swift BAT and Fermi GBM increases the possibility of detecting and localizing the electromagnetic counterparts of gravitational wave events seen by the new generation of gravitational wave detectors.


title: An Explicit Upper Bound for Modulus of Divided Difference on a Jordan Arc in the Complex Plane

abstract: An explicit upper bound is derived for the modulus of divided difference for a smooth(not necessarily analytic) function defined on a smooth Jordan arc (or a smooth Jordan curve) in the complex plane. As an immediate application, an error estimate for complex polynomial interpolation on a Jordan arc (or a Jordan curve) is given, which extends the well-known error estimate for polynomial interpolation on the unit interval. Moreover, this upper bound is independent of the parametrization of the curve.


title: Time-resolved observation of coherent excitonic nonlinear response with a table-top narrowband THz pulse wave

abstract: By combining a tilted-pulse-intensity-front scheme using a LiNbO3 crystal and a chirped-pulse-beating method, we generated a narrowband intense terahertz (THz) pulse, which had a maximum electric field of more than 10 kV/cm at around 2 THz, a bandwidth of about 50 GHz,and frequency tunability from 0.5 to 2 THz. By performing THz-pump and near-infrared-probe experiments on GaAs quantum wells, we observed that the resonant excitation of the intraexcitonic 1s-2p transition induces a clear and large Autler-Townes splitting. Our time-resolved measurements show that the splitting energy observed in the rising edge region of electric field is larger than in the constant region. This result implies that the splitting energy depends on the time-averaged THz field over the excitonic dephasing time rather than that at the instant of the exciton creation by a probe pulse.


title: Vertical Distribution Relations For Special Cycles on Unitary Shimura Varieties

abstract: We consider cycles on a 3-dimensional Shimura varieties attached to a unitary group, defined over extensions of a CM field $E$, which appear in the context of the conjectures of Gan, Gross, and Prasad \cite{gan-gross-prasad}. We establish a vertical distribution relation for these cycles over an anticyclotomic extension of $E$, complementing the horizontal distribution relation of \cite{jetchev:unitary}, and use this to define a family of norm-compatible cycles over these fields, thus obtaining a universal norm construction similar to the Heegner $\Lambda$-module constructed from Heegner points.


title: Mean-Field Inference in Gaussian Restricted Boltzmann Machine

abstract: A Gaussian restricted Boltzmann machine (GRBM) is a Boltzmann machine defined on a bipartite graph and is an extension of usual restricted Boltzmann machines. A GRBM consists of two different layers: a visible layer composed of continuous visible variables and a hidden layer composed of discrete hidden variables. In this paper, we derive two different inference algorithms for GRBMs based on the naive mean-field approximation (NMFA). One is an inference algorithm for whole variables in a GRBM, and the other is an inference algorithm for partial variables in a GBRBM. We compare the two methods analytically and numerically and show that the latter method is better.


title: Next Generation Muon g-2 Experiments

abstract: I report on the progress of two new muon anomalous magnetic moment experiments, which are in advanced design and construction phases. The goal of Fermilab E989 is to reduce the experimental uncertainty of $a_\mu$ from Brookhaven E821 by a factor of 4; that is, $\delta a_\mu \sim 16 \times 10^{-11}$, a relative uncertainty of 140~ppb. The method follows the same magic-momentum storage ring concept used at BNL, and pioneered previously at CERN, but muon beam preparation, storage ring internal hardware, field measuring equipment, and detector and electronics systems are all new or upgraded significantly. In contrast, J-PARC E34 will employ a novel approach based on injection of an ultra-cold, low-energy, muon beam injected into a small, but highly uniform magnet. Only a small magnetic focusing field is needed to maintain storage, which distinguishes it from CERN, BNL and Fermilab. E34 aims to roughly match the previous BNL precision in their Phase~1 installation.


title: Excluded volume effect enhances the homology pairing of model chromosomes

abstract: To investigate the structural dynamics of the homology pairing of polymers, we mod- eled the scenario of homologous chromosome pairings during meiosis in Schizosaccharomyces pombe, one of the simplest model organisms of eukaryotes. We consider a simple model consist- ing of pairs of homologous polymers with the same structures that are confined in a cylindrical container, which represents the local parts of chromosomes contained in an elongated nucleus of S. pombe. Brownian dynamics simulations of this model showed that the excluded volume effects among non-homological chromosomes and the transitional dynamics of nuclear shape serve to enhance the pairing of homologous chromosomes.


title: The $\mathcal{L}$-invariant, the dual $\mathcal{L}$-invariant, and families

abstract: Given a rank two trianguline family of $(\varphi,\Gamma)$-modules having a noncrystalline semistable member, we compute the Fontaine--Mazur $\mathcal{L}$-invariant of that member in terms of the logarithmic derivative, with respect to the Sen weight, of the value at p of the trianguline parameter. This generalizes prior work, in the case of Galois representations, due to Greenberg--Stevens and Colmez.


title: Quantum cluster algorithm for frustrated Ising models in a transverse field

abstract: Working within the stochastic series expansion framework, we introduce and characterize a new quantum cluster algorithm for quantum Monte Carlo simulations of transverse field Ising models with frustrated Ising exchange interactions. As a demonstration of the capabilities of this new algorithm, we show that a relatively small, ferromagnetic next-nearest neighbour coupling drives the transverse field Ising antiferromagnet on the triangular lattice from an antiferromagnetic three-sublattice ordered state at low temperature to a ferrimagnetic three-sublattice ordered state.


title: The Indian Spontaneous Expression Database for Emotion Recognition

abstract: Automatic recognition of spontaneous facial expressions is a major challenge in the field of affective computing. Head rotation, face pose, illumination variation, occlusion etc. are the attributes that increase the complexity of recognition of spontaneous expressions in practical applications. Effective recognition of expressions depends significantly on the quality of the database used. Most well-known facial expression databases consist of posed expressions. However, currently there is a huge demand for spontaneous expression databases for the pragmatic implementation of the facial expression recognition algorithms. In this paper, we propose and establish a new facial expression database containing spontaneous expressions of both male and female participants of Indian origin. The database consists of 428 segmented video clips of the spontaneous facial expressions of 50 participants. In our experiment, emotions were induced among the participants by using emotional videos and simultaneously their self-ratings were collected for each experienced emotion. Facial expression clips were annotated carefully by four trained decoders, which were further validated by the nature of stimuli used and self-report of emotions. An extensive analysis was carried out on the database using several machine learning algorithms and the results are provided for future reference. Such a spontaneous database will help in the development and validation of algorithms for recognition of spontaneous expressions.


title: Probabilistic Integration: A Role in Statistical Computation?

abstract: A research frontier has emerged in scientific computation, wherein numerical error is regarded as a source of epistemic uncertainty that can be modelled. This raises several statistical challenges, including the design of statistical methods that enable the coherent propagation of probabilities through a (possibly deterministic) computational work-flow. This paper examines the case for probabilistic numerical methods in routine statistical computation. Our focus is on numerical integration, where a probabilistic integrator is equipped with a full distribution over its output that reflects the presence of an unknown numerical error. Our main technical contribution is to establish, for the first time, rates of posterior contraction for these methods. These show that probabilistic integrators can in principle enjoy the "best of both worlds", leveraging the sampling efficiency of Monte Carlo methods whilst providing a principled route to assess the impact of numerical error on scientific conclusions. Several substantial applications are provided for illustration and critical evaluation, including examples from statistical modelling, computer graphics and a computer model for an oil reservoir.

categories: stat.ML cs.NA math.NA math.ST stat.CO stat.TH


title: Role of anthropogenic direct heat emissions in global warming

abstract: The anthropogenic emissions of greenhouse gases (GHG) are widely realized as the predominant drivers of global warming, but the huge and increasing anthropogenic direct heat emissions (AHE) has not gained enough attention in terms of its role in the warming of the climate system. Based on two reasonable assumptions of (1) AHE eventually transfers to the Earth energy system and (2) the net warming is only driven by the net radioactive forcing (RF) from either GHG or other causes, we analyzed the role of AHE in global warming. The mean annual total AHE of the four main sources including energy consumption, residual heat of electricity generation, biomass decomposition by land use and cover change (LUCC) and food consumption was estimated to be 4.41*10^20 J in 1970-2010, accounting for 6.23% of the net annual heat increase of the Earth reported by IPCC AR5 for the period. The mean annual radioactive forcing (RF) by AHE was up to 29.94 mW m^(-2) globally in 1981-2010, less than the annual net increase of total GHG (34.80 mW m^(-2)) but more that of CO2 (24.43 mW m^(-2). The results suggested that AHE played a great role in global warming, urging us to pay more attention to DHE for a better evaluation of RF and more reasonable energy policy including non-GHG energy in the future (Keywords: climate change, anthropogenic direct heat emission, energy, radioactive forcing, heat balance, global warming).


title: Non-cyclic graphs of (non)orientable genus one

abstract: Let $G$ be a finite non-cyclic group. The non-cyclic graph $\Gamma_G$ of $G$ is the graph whose vertex set is $G\setminus Cyc(G)$, two distinct vertices being adjacent if they do not generate a cyclic subgroup, where $Cyc(G)=\{a\in G: \langle a,b\rangle\ \text{is cyclic for each } b\in G\}$. In this paper, we classify all finite non-cyclic groups $G$ such that $\Gamma_G$ has (non)orientable genus one.


title: Atom-Field Entanglement in Cavity QED: Nonlinearity and Saturation

abstract: We investigate the degree of entanglement between an atom and a driven cavity mode, in the presence of dissipation. Previous work has shown that in the limit of weak driving fields, the steady state entanglement is proportional to the square of the driving intensity. This quadratic dependence is due to the generation of entanglement by the creation of pairs of photons/excitations. In this work we investigate the entanglement between an atom and a cavity in the presence of multiple photons. Nonlinearity of the atomic response is needed to generate entanglement, but as that nonlinearity saturates the entanglement vanishes. We posit that this is due to spontaneous emission, which puts the atom in the ground state and the atom-field state into a direct product state. An intermediate value of the driving field, near the field that saturates the atomic response, optimizes the atom-field entanglement. In a parameter regime for which multiphoton resonances occur, we find that entanglement recurs at those resonances. \cite{Shamailov-10} In this regime, we find that the entanglement decreases with increaing photon number. We also investigate, in the bimodal regime of Alsing et. al.\cite{Alsing-91-DressedState} \cite{Alsing-92-DSS}\cite{Carm-15-X}, the entanglement as a function of atom and/or cavity detuning. Here we find that there is evidence of a phase transition in the entanglement, which occurs at $2\epsilon/g \geq 1$.


title: On the non-thermal kappa-distributed electrons in planetary nebulae and HII regions: the kappa index and its correlations with other nebular properties

abstract: Recently, a suspicion arose that the free electrons in planetary nebulae (PNe) and HII regions might have non-thermal energy distributions. In this scenario, a kappa index is introduced to characterize the electron energy distributions, with smaller kappa values indicating larger deviations from Maxwell-Boltzmann distributions. Assuming that this is the case, we determine the kappa values for a sample of PNe and HII regions by comparing the intensities of [OIII] collisionally excited lines and the hydrogen Balmer jump. We find the average kappa indices of PNe and HII regions to be 27 and 32, respectively. Correlations between the resultant kappa values and various physical properties of the nebulae are examined to explore the potential origin of non-thermal electrons in photoionized gaseous nebulae. However, no positive result is obtained. Thus the current analysis does not lend to support to the idea that kappa-distributed electrons are present in PNe and HII regions.


title: Entropy approximation versus uniqueness of equilibrium for a dense affine space of continuous functions

abstract: We show that for a $\mathbb{Z}^{l}$-action (or $(\N\cup\{0\})^l$-action) on a non-empty compact metrizable space $\Omega$, the existence of a affine space dense in the set of continuous functions on $\Omega$ constituted by elements admitting a unique equilibrium state implies that each invariant measure can be approximated weakly$^*$ and in entropy by a sequence of measures which are unique equilibrium states.


title: A Literature Survey of various Fingerprint De-noising Techniques to justify the need of a new De-noising model based upon Pixel Component Analysis

abstract: Image Preprocessing is a vital step in the field of image processing for biometric pattern recognition. This paper studies and reviews various classical and modern fingerprint image de-noising models. The various model used for de-noising ranges widely from transform matrix using frequency, histogram model de-noising, de-noising by introducing Gabor filter and its types to enhance fingerprint images. The output efficiency of various de-noising model proposed earlier is calculated on the basis of SNR (signal to noise ratio) and MSE (mean square error rate). Our simulated experimental results indicates that incorporating the de-noising model based on Gabor filter inside domain of wavelet ranges with composite method only betters MSE (Mean Square Error). Improved MSE without significant improvement in SNR improves the fingerprint images only by a little margin which is non-optimal in nature. Thus the objective of this research paper is to build an optimal de-noising model for fingerprint images so that its usage in biometric authentication can be more robust in nature.


title: An optical spectrum of a large isolated gas-phase PAH cation: C78H26+

abstract: A gas-phase optical spectrum of a large polycyclic aromatic hydrocarbon (PAH) cation - C78H26 +- in the 410-610 nm range is presented. This large all-benzenoid PAH should be large enough to be stable with respect to photodissociation in the harsh conditions prevailing in the interstellar medium (ISM). The spectrum is obtained via multi-photon dissociation (MPD) spectroscopy of cationic C78H26 stored in the Fourier Transform Ion Cyclotron Resonance (FT-ICR) cell using the radiation from a mid-band optical parametric oscillator (OPO) laser. The experimental spectrum shows two main absorption peaks at 431 nm and 516 nm, in good agreement with a theoretical spectrum computed via time-dependent density functional theory (TD-DFT). DFT calculations indicate that the equilibrium geometry, with the absolute minimum energy, is of lowered, nonplanar C2 symmetry instead of the more symmetric planar D2h symmetry that is usually the minimum for similar PAHs of smaller size. This kind of slightly broken symmetry could produce some of the fine structure observed in some diffuse interstellar bands (DIBs). It can also favor the folding of C78H26 + fragments and ultimately theformation of fullerenes. This study opens up the possibility to identify the most promising candidates for DIBs amongst large cationic PAHs.


title: Mott transitions in a three-component Falicov-Kimball model: A slave boson mean-field study

abstract: Metal-insulator transitions in a three-component Falicov-Kimball model are investigated within the Kotliar-Ruckenstein slave boson mean-field approach. The model describes a mixture of two interacting fermion atom species loaded into an optical lattice at ultralow temperature. One species is two-component atoms, which can hop in the optical lattice, and the other is single-component atoms, which are localized. Different correlation-driven metal-insulator transitions are observed depending on the atom filling conditions and local interactions. These metal-insulator transitions are classified by the band renormalization factors and the double occupancies of the atom species. The filling conditions and the critical value of the local interactions for these metal-insulator transitions are also analytically established. The obtained results not only are in good agreement with the dynamical mean-field theory for the three-component Falicov-Kimball model but also clarify the nature and properties of the metal-insulator transitions in a simple physics picture.


title: Derivation of the Hall and Extended Magnetohydrodynamics Brackets

abstract: There are several plasma models intermediate in complexity between ideal magnetohydrodynamics (MHD) and two-fluid theory, with Hall and Extended MHD being two important examples. In this paper we investigate several aspects of these theories, with the ultimate goal of deriving the noncanonical Poisson brackets used in their Hamiltonian formulations. We present fully Lagrangian actions for each, as opposed to the fully Eulerian, or mixed Eulerian-Lagrangian, actions that have appeared previously. As an important step in this process we exhibit each theory's two advected fluxes (in analogy to ideal MHD's advected magnetic flux), discovering also that with the correct choice of gauge they have corresponding Lie-dragged potentials resembling the electromagnetic vector potential, and associated conserved helicities. Finally, using the Euler-Lagrange map, we show how to derive the noncanonical Eulerian brackets from canonical Lagrangian ones.

categories: physics.plasm-ph astro-ph.IM math-ph math.MP


title: A Combinatorial Problem Related to Sparse Systems of Equations

abstract: Nowadays sparse systems of equations occur frequently in science and engineering. In this contribution we deal with sparse systems common in cryptanalysis. Given a cipher system, one converts it into a system of sparse equations, and then the system is solved to retrieve either a key or a plaintext. Raddum and Semaev proposed new methods for solving such sparse systems. It turns out that a combinatorial MaxMinMax problem provides bounds on the average computational complexity of sparse systems. In this paper we initiate a study of a linear algebra variation of this MaxMinMax problem.


title: Modified Non-Euclidian Transformation on the SO(2N+2)/U(N+1) Grassmannian and SO(2N+1) Random Phase Approximation for Unified Description of Bose and Fermi Type Collective Excitations

abstract: In a slight different way from the previous one, we propose a modified non-Euclidian transformation on the SO(2N+2)/U(N+1) Grassmannian which give the projected SO(2N+1) Tamm-Dancoff equation. We derive a classical time dependent (TD) SO(2N+1) Lagrangian which, through the Euler-Lagrange equation of motion for SO(2N+2)/U(N+1) coset variables, brings another form of the previous extended-TD Hartree-Bogoliubov (HB) equation. The SO(2N+1) random phase approximation (RPA) is derived using Dyson representation for paired and unpaired operators. In the SO(2N) HB case, one boson and two boson excited states are realized. We, however, stress non existence of a higher RPA vacuum. An integrable system is given by a geometrical concept of zero-curvature, i.e., integrability condition of connection on the corresponding Lie group. From the group theoretical viewpoint, we show the existence of a symplectic two-form omega.


title: Fundamental Limits of Exciton-Exciton Annihilation for Light Emission in Transition Metal Dichalcogenide Monolayers

abstract: We quantitatively illustrate the fundamental limit that exciton-exciton annihilation (EEA) may impose to the light emission of monolayer transition metal dichalcogenide (TMDC) materials. The EEA in TMDC monolayers shows dependence on the interaction with substrates as its rate increases from 0.1 cm2/s (0.05 cm2/s) to 0.3 cm2/s (0.1 cm2/s) with the substrates removed for WS2 (MoS2) monolayers. It turns to be the major pathway of exciton decay and dominates the luminescence efficiency when the exciton density is beyond 1010 cm-2 in suspended monolayers or 1011 cm-2 in supported monolayers. This sets an upper limit on the density of injected charges in light emission devices for the realization of optimal luminescence efficiency. The strong EEA rate also dictates the pumping threshold for population inversion in the monolayers to be 12-18 MW/cm2 (optically) or 2.5-4x105 A/cm2 (electrically).


title: Urban Skylines: building heights and shapes as measures of city size

abstract: The shape of buildings plays a critical role in the energy efficiency, lifestyles, land use and infrastructure systems of cities. Thus, as most of the world's cities continue to grow and develop, understanding the interplay between the characteristics of urban environments and the built form of cities is essential to achieve local and global sustainability goals. Here, we compile and analyze the most extensive data set of building shapes to date, covering more than 4.8 million individual buildings across several major cities in North America. We show that average building height increases systematically with city size and follows theoretical predictions derived from urban scaling theory. We also study the allometric relationship between surface area and volume of buildings in terms of characteristic shape parameters. This allows us to demonstrate that the reported trend towards higher (and more voluminous) buildings effectively decreases the average surface-to-volume ratio, suggesting potentially significant energy savings with growing city size. At the same time, however, the surface-to-volume ratio increases in the downtown cores of large cities, due to shape effects and specifically to the proliferation of tall, needlelike buildings. Thus, the issue of changes in building shapes with city size and associated energy management problem is highly heterogeneous. It requires a systematic approach that includes the factors that drive the form of built environments, entangling physical, infrastructural and socioeconomic aspects of cities.


title: A New Statistical Framework for Genetic Pleiotropic Analysis of High Dimensional Phenotype Data

abstract: The widely used genetic pleiotropic analysis of multiple phenotypes are often designed for examining the relationship between common variants and a few phenotypes. They are not suited for both high dimensional phenotypes and high dimensional genotype (next-generation sequencing) data. To overcome these limitations, we develop sparse structural equation models (SEMs) as a general framework for a new paradigm of genetic analysis of multiple phenotypes. To incorporate both common and rare variants into the analysis, we extend the traditional multivariate SEMs to sparse functional SEMs. To deal with high dimensional phenotype and genotype data, we employ functional data analysis and the alternative direction methods of multiplier (ADMM) techniques to reduce data dimension and improve computational efficiency. Using large scale simulations we showed that the proposed methods have higher power to detect true causal genetic pleiotropic structure than other existing methods. Simulations also demonstrate that the gene-based pleiotropic analysis has higher power than the single variant-based pleiotropic analysis. The proposed method is applied to exome sequence data from the NHLBI Exome Sequencing Project (ESP) with 11 phenotypes, which identifies a network with 137 genes connected to 11 phenotypes and 341 edges. Among them, 114 genes showed pleiotropic genetic effects and 45 genes were reported to be associated with phenotypes in the analysis or other cardiovascular disease (CVD) related phenotypes in the literature.


title: Besov spaces of self-affine lattice tilings and pointwise regularity

abstract: We investigate Besov spaces of self-affine tilings of ${\Bbb R}^{n}$ and discuss various characterizations of those Besov spaces. We see what is a finite set of functions which generates the Besov spaces from a view of multiresolution approximation on self-affine lattice tilings of ${\Bbb R}^{n}$. Using this result we give a generalization of already known characterizations of Besov spaces given by wavelet expansion and we apply to study the pointwise H${\ddot {\rm o}}$lder space. Furthermore we give descriptions of scaling exponents measured by Besov spaces, and estimations of a pointwise H${\ddot {\rm o}}$lder exponent to compute the pointwise scaling exponent of several oscillatory functions.


title: Multivariate Functional Regression Models for Epistasis Analysis

abstract: To date, most genetic analyses of phenotypes have focused on analyzing single traits or, analyzing each phenotype independently. However, joint epistasis analysis of multiple complementary traits will increase statistical power, and hold the key to understanding the complicated genetic structure of the complex diseases. Despite their importance in uncovering the genetic structure of complex traits, the statistical methods for identifying epistasis in multiple phenotypes remains fundamentally unexplored. To fill this gap, we formulate a test for interaction between two gens in multiple quantitative trait analysis as a multiple functional regression (MFRG) in which the genotype functions (genetic variant profiles) are defined as a function of the genomic position of the genetic variants. We use large scale simulations to calculate its type I error rates for testing interaction between two genes with multiple phenotypes and to compare its power with multivariate pair-wise interaction analysis and single trait interaction analysis by a single variate functional regression model. To further evaluate its performance, the MFRG for epistasis analysis is applied to five phenotypes and exome sequence data from the NHLBI Exome Sequencing Project (ESP) to detect pleiotropic epistasis. A total of 136 pairs of genes that formed a genetic interaction network showed significant evidence of epistasis influencing five traits. The results demonstrate that the joint interaction analysis of multiple phenotypes has much higher power to detect interaction than the interaction analysis of single trait and may open a new direction to fully uncovering the genetic structure of multiple phenotypes.


title: Contact geometric descriptions of vector fields on dually flat spaces and their applications in electric circuit models and nonequilibrium statistical mechanics

abstract: Contact geometry has been applied to various mathematical sciences, and it has been proposed that a contact manifold and a strictly convex function induce a dually flat space that is used in information geometry. Here, such a dually flat space is related to a Legendre submanifold in a contact manifold. In this paper contact geometric descriptions of vector fields on dually flat spaces are proposed on the basis of the theory of contact Hamiltonian vector fields. Based on these descriptions, two ways of lifting vector fields on Legendre submanifolds to contact manifolds are given. For some classes of these lifted vector fields, invariant measures in contact manifolds and stability analysis around Legendre submanifolds are explicitly given. Throughout this paper, Legendre duality is explicitly stated. In addition, to show how to apply these general methodologies to applied mathematical disciplines, electric circuit models and some examples taken from nonequilibrium statistical mechanics are analyzed.

categories: math-ph cond-mat.stat-mech math.MP nlin.CD


title: Real and virtual $N\bar{N}$ pair production near the threshold

abstract: Nucleon-antinucleon optical potential, which explains the experimental data for the processes $e^{+}e^{-}\rightarrow p\bar{p}$ and $e^{+}e^{-}\rightarrow\mbox{pions}$ near the threshold of $p\bar{p}$ pair production, is suggested. To obtain this potential we have used the available experimental data for $p\bar{p}$ scattering, $p\bar{p}$ pair production in $e^{+}e^{-}$ annihilation, and the ratio of electromagnetic form factors of a proton in the timelike region. It turns out that final-state interaction via the optical potential allows one to reproduce the available experimental data with good accuracy. Our results for the cross sections of $e^{+}e^{-}\to6\pi$ process near the threshold of $p\bar{p}$ pair production are in agreement with the recent experiments.


title: Pion Production via Proton Synchrotron Radiation in Strong Magnetic Fields in Relativistic Field Theory: Scaling Relations and Angular Distributions

abstract: We study pion production by proton synchrotron radiation in the presence of a strong magnetic field when the Landau numbers of the initial and final protons are $n_{i,f} \sim 10^4 - 10^5$. We find in our relativistic field theory calculations that the pion decay width depends only on the field strength parameter which previously was only conjectured based upon semi-classical arguments. Moreover, we also find new results that the decay width satisfies a robust scaling relation, and that the polar angular distribution of emitted pion momenta is very narrow and can be easily obtained. This scaling implies that one can infer the decay width in more realistic magnetic fields of $10^{15}$G, where $n_{i,f} \sim 10^{12} - 10^{13}$, from the results for $n_{i,f} \sim 10^4 - 10^5$. The resultant pion intensity and angular distributions for realistic magnetic field strengths are presented and their physical implications discussed.


title: Necessary optimality conditions for optimal control problems with nonsmooth mixed state and control constraints

abstract: In this paper we study an optimal control problem with nonsmooth mixed state and control constraints. In most of the existing results, the necessary optimality condition for optimal control problems with mixed state and control constraints are derived under the Mangasarian-Fromovitz condition and under the assumption that the state and control constraint functions are smooth. In this paper we derive necessary optimality conditions for problems with nonsmooth mixed state and control constraints under constraint qualifications based on pseudo-Lipschitz continuity and calmness of certain set-valued maps. The necessary conditions are stratified, in the sense that they are asserted on precisely the domain upon which the hypotheses (and the optimality) are assumed to hold. Moreover necessary optimality conditions with an Euler inclusion taking an explicit multiplier form are derived for certain cases.


title: Renormalization Scheme Dependence in b->u Semileptonic Decays

abstract: After reviewing how the renormalization group equation can be used to sum logarithmic corrections to the decay rate for the semi-leptonic process b->u when using minimal subtraction, we consider renormalization scheme dependence for this calculation when employing this renor- malization scheme. In this calculation, an ambiguity resides in the running strong coupling and in the running b quark mass. The ambiguity usually associated with the renormalization mass scale mu is shown to cancel. It is demonstrated how in one renormalization scheme, there are only leading-log contributions to this decay rate. Another choice of renormalization scheme results in the decay rate being expressed in terms of the two-loop contribution to the beta-function associated with the strong coupling and the one-loop contribution to the anomalous mass dimension as well as a set of renormalization scheme invariant parameters.


title: Nonresonance adiabatic photon confinement in spherical mirror system. Experimental study

abstract: New nonresonance approach of photon accumulation in two spherical mirrors has been experimentally demonstrated. In this work, we have received a high accumulation coefficient and shown good effectiveness of this technique for creating photo neutralizer of negative ion beams. This efficiency in such concept is generally determined by reflectance mirrors and is practically not dependent of input emission quality and does not require high precision adjusting the optical elements.


title: Bond Graph Modelling of Chemoelectrical Energy Transduction

abstract: Energy-based bond graph modelling of biomolecular systems is extended to include chemoelectrical trans- duction thus enabling integrated thermodynamically-compliant modelling of chemoelectrical systems in general and excitable membranes in particular. Our general approach is illustrated by recreating a well-known model of an excitable membrane. This model is used to investigate the energy consumed during a membrane action potential thus contributing to the current debate on the trade-off between the speed of an action potential event and energy consumption. The influx of Na+ is often taken as a proxy for energy consumption; in contrast, this paper presents an energy based model of action potentials. As the energy based approach avoids the assumptions underlying the proxy approach it can be directly used to compute energy consumption in both healthy and diseased neurons. These results are illustrated by comparing the energy consumption of healthy and degenerative retinal ganglion cells using both simulated and in vitro data.


title: SU(3) gauge theory with four degenerate fundamental fermions on the lattice

abstract: As a part of the project studying large $N_f$ QCD, the LatKMI Collaboration has been investigating the SU(3) gauge theory with four fundamental fermions (four-flavor QCD). The main purpose of studying four-flavor QCD is to provide a qualitative comparison to $N_f= 8$, $12$, $16$ QCD; however, a quantitative comparison to real-world QCD is also interesting. To make such comparisons more meaningful, it is desirable to use the same kind of lattice action consistently, so that qualitative difference of different theories are less affected by artifacts of lattice discretization. Here, we adopt the highly-improved staggered quark action with the tree-level Symanzik gauge action (HISQ/tree), which is exactly the same as the setup for our simulations for $SU(3)$ gauge theories with $N_f=8$, $12$ and $16$ fundamental fermions~\cite{Aoki:2013xza, Aoki:2012eq, Aoki:2014oma}. In the next section, we show the fermion mass dependence of $F_\pi$, $\langle\bar{\psi}\psi\rangle$, $M_\pi$, $M_\rho$, $M_N$ and their chiral extrapolations. In section 3, preliminary results of the measurement of the mass of the flavor-singlet scalar bound state will be reported.


title: Large deviations for near-extreme eigenvalues in the beta-ensembles

abstract: For beta ensembles with convex poynomial potentials, we prove a large deviation principle for the empirical spectral distribution seen from the rightmost particle. This modified spectral distribution was introduced by Perret and Schehr (J. Stat. Phys. 2014) to study the crowding near the maximal eigenvalue, in the case of the GUE. We prove also convergence of fluctuations.


title: Applicability of Taylor's hypothesis in thermally driven turbulence

abstract: In this paper, we show that in the presence of large-scale circulation (LSC), Taylor's hypothesis can be invoked to deduce the energy spectrum in thermal convection using real space probes, a popular experimental tool. We perform numerical simulation of turbulent convection in a cube and observe that the velocity field follows Kolmogorov's spectrum ($k^{-5/3}$). We also record the velocity time series using real space probes near the lateral walls. The corresponding frequency spectrum exhibits Kolmogorov's spectrum ($f^{-5/3}$), thus validating Taylor's hypothesis with the steady LSC playing the role of a mean velocity field. The aforementioned findings based on real space probes provide valuable inputs for experimental measurements used for studying the spectrum of convective turbulence.


title: Systematic study of complete fusion suppression in reactions involving weakly bound nuclei at energies above the Coulomb barrier

abstract: Complete fusion excitation functions of reactions involving breakup are studied by using the empirical coupled-channel (ECC) model with breakup effects considered. An exponential function with two parameters is adopted to describe the prompt-breakup probability in the ECC model. These two parameters are fixed by fitting the measured prompt-breakup probability or the complete fusion cross sections. The suppression of complete fusion at energies above the Coulomb barrier is studied by comparing the data with the predictions from the ECC model without the breakup channel considered. The results show that the suppression of complete fusion are roughly independent of the target for the reactions involving the same projectile.


title: Triplet Spike Time Dependent Plasticity: A floating-gate Implementation

abstract: Synapse plays an important role of learning in a neural network; the learning rules which modify the synaptic strength based on the timing difference between the pre- and post-synaptic spike occurrence is termed as Spike Time Dependent Plasticity (STDP). The most commonly used rule posits weight change based on time difference between one pre- and one post spike and is hence termed doublet STDP (DSTDP). However, D-STDP could not reproduce results of many biological experiments; a triplet STDP (T-STDP) that considers triplets of spikes as the fundamental unit has been proposed recently to explain these observations. This paper describes the compact implementation of a synapse using single floating-gate (FG) transistor that can store a weight in a nonvolatile manner and demonstrate the triplet STDP (T-STDP) learning rule by modifying drain voltages according to triplets of spikes. We describe a mathematical procedure to obtain control voltages for the FG device for T-STDP and also show measurement results from a FG synapse fabricated in TSMC 0.35um CMOS process to support the theory. Possible VLSI implementation of drain voltage waveform generator circuits are also presented with simulation results.


title: A New Infinite Family of Hemisystems of the Hermitian Surface

abstract: In this paper, we construct an infinite family of hemisystems of the Hermitian surface $\mathsf{H}(3,q^2)$. In particular, we show that for every odd prime power $q$ congruent to $3$ modulo $4$, there exists a hemisystem of $\mathsf{H}(3,q^2)$ admitting $C_{(q^3+1)/4} : C_3$.


title: Stability Boundaries and Sufficient Stability Conditions for Stably Stratified, Monotonic Shear Flows

abstract: Linear stability of inviscid, parallel, and stably stratified shear flow is studied under the assumption of smooth strictly monotonic profiles of shear flow and density, so that the local Richardson number is positive everywhere. The marginally unstable modes are systematically found by solving a one-parameter family of regular Sturm-Liouville problems, which can determine the stability boundaries more efficiently than solving the Taylor-Goldstein equation directly. By arguing for the non-existence of a marginally unstable mode, we derive new sufficient conditions for stability, which generalize the Rayleigh-Fj{\o}rtoft criterion for unstratified shear flows.


title: Modeling Human Understanding of Complex Intentional Action with a Bayesian Nonparametric Subgoal Model

abstract: Most human behaviors consist of multiple parts, steps, or subtasks. These structures guide our action planning and execution, but when we observe others, the latent structure of their actions is typically unobservable, and must be inferred in order to learn new skills by demonstration, or to assist others in completing their tasks. For example, an assistant who has learned the subgoal structure of a colleague's task can more rapidly recognize and support their actions as they unfold. Here we model how humans infer subgoals from observations of complex action sequences using a nonparametric Bayesian model, which assumes that observed actions are generated by approximately rational planning over unknown subgoal sequences. We test this model with a behavioral experiment in which humans observed different series of goal-directed actions, and inferred both the number and composition of the subgoal sequences associated with each goal. The Bayesian model predicts human subgoal inferences with high accuracy, and significantly better than several alternative models and straightforward heuristics. Motivated by this result, we simulate how learning and inference of subgoals can improve performance in an artificial user assistance task. The Bayesian model learns the correct subgoals from fewer observations, and better assists users by more rapidly and accurately inferring the goal of their actions than alternative approaches.


title: Neural Enquirer: Learning to Query Tables with Natural Language

abstract: We proposed Neural Enquirer as a neural network architecture to execute a natural language (NL) query on a knowledge-base (KB) for answers. Basically, Neural Enquirer finds the distributed representation of a query and then executes it on knowledge-base tables to obtain the answer as one of the values in the tables. Unlike similar efforts in end-to-end training of semantic parsers, Neural Enquirer is fully "neuralized": it not only gives distributional representation of the query and the knowledge-base, but also realizes the execution of compositional queries as a series of differentiable operations, with intermediate results (consisting of annotations of the tables at different levels) saved on multiple layers of memory. Neural Enquirer can be trained with gradient descent, with which not only the parameters of the controlling components and semantic parsing component, but also the embeddings of the tables and query words can be learned from scratch. The training can be done in an end-to-end fashion, but it can take stronger guidance, e.g., the step-by-step supervision for complicated queries, and benefit from it. Neural Enquirer is one step towards building neural network systems which seek to understand language by executing it on real-world. Our experiments show that Neural Enquirer can learn to execute fairly complicated NL queries on tables with rich structures.


title: Viscous singular shock profiles for the Keyfitz-Kranzer system

abstract: It was shown by Schecter (2004, J. Differential Equations, 205, 185-210), using the methods of Geometric Singular Perturbation Theory, that the Dafermos regularization $u_t+f(u)_x= \epsilon tu_{xx}$ for the Keyfitz-Kranzer system admits an unbounded family of solutions. Inspired by that work, in this paper we provide a more intuitive approach which leads to a stronger result. In addition to the existence of viscous profiles, we also prove the weak convergence and show that the maximum of the solution is of order $\epsilon^{-2}$. This asymptotic behavior is distinct from that obtained in the author's recent work (arXiv:1512.00394) on a system modeling two-phase fluid flow, for which the maximum of the viscous solution is of order $\exp(\epsilon^{-1})$.


title: Duality at Lattice Scale around the Takhtajan-Babujain Point in the Bilinear Biquadratic Spin-1 Chain

abstract: In the literature, it is known that near the Takhtajan-Babujain point of the bilinear-biquadratic spin-$1$ chain, the low energy theory is captured by three Majorana fermions with a mass sign change, or three Ising chains near the transition. In this paper, we study how the spin-$1$, the Ising spins, and the fermions are related at lattice scale. We construct the spin-$1$ operator out from the Ising spins and fermions, showing that when the Ising chains are in the ordered phase, the spin-$1$ chain is in the dimer phase, and when the Ising chains are in the paramagnetic phase, the spin-$1$ chain is in the Haldane phase. Through this mapping we can see how unit cell doubling in the dimer phase and the boundary spin-$\frac12$ zero mode in the Haldane phase take place, even though in the Ising transition neither phenomena exist.


title: Compact pluricanonical manifolds are Vaisman

abstract: A locally conformally Kahler manifold is a Hermitian manifold $(M,I,\omega)$ satisfying $d\omega=\theta\wedge \omega$, where $\theta$ is a closed 1-form, called the Lee form of $M$. It is called pluricanonical if $\nabla\theta$ is of Hodge type $(2,0)+(0,2)$, where $\nabla$ is the Levi-Civita connection, and Vaisman if $\nabla\theta=0$. We show that a compact LCK manifold is pluricanonical if and only if the Lee form has constant length and the Kahler form of its covering admits an automorphic potential. Using a degenerate Monge-Ampere equation and the classification of surfaces of Kahler rank one, due to Brunella, Chiose and Toma, we show that any pluricanonical metric on a compact manifold is Vaisman. Several errata to our previous work are given in the last Section.


title: Posterior Belief Assessment: Extracting Meaningful Subjective Judgements from Bayesian Analyses with Complex Statistical Models

abstract: In this paper, we are concerned with attributing meaning to the results of a Bayesian analysis for a problem which is sufficiently complex that we are unable to assert a precise correspondence between the expert probabilistic judgements of the analyst and the particular forms chosen for the prior specification and the likelihood for the analysis. In order to do this, we propose performing a finite collection of additional Bayesian analyses under alternative collections of prior and likelihood modelling judgements that we may also view as representative of our prior knowledge and the problem structure, and use these to compute posterior belief assessments for key quantities of interest. We show that these assessments are closer to our true underlying beliefs than the original Bayesian analysis and use the temporal sure preference principle to establish a probabilistic relationship between our true posterior judgements, our posterior belief assessment and our original Bayesian analysis to make this precise. We exploit second order exchangeability in order to generalise our approach to situations where there are infinitely many alternative Bayesian analyses we might consider as informative for our true judgements so that the method remains tractable even in these cases. We argue that posterior belief assessment is a tractable and powerful alternative to robust Bayesian analysis. We describe a methodology for computing posterior belief assessments in even the most complex of statistical models and illustrate with an example of calibrating an expensive ocean model in order to quantify uncertainty about global mean temperature in the real ocean.


title: A new enlightenment about gaps between primes

abstract: The idea of generating prime numbers through sequence of sets of co-primes was the starting point of this paper that ends up by proving two conjectures, the existence of infinitely many twin primes and the Goldbach conjecture. The main idea of our approach is summarized on the creation and on the analyzing sequence of sets of distinct co-primes with the first $n$ primes, $\left\{ p_i :\, i\leq n \right\}$, and the important properties of the modulus linear combination of the co-prime sets, $H=\left(1,p_{n+1},..., \Pi_{i=1}^n p_i-1\right) $, that gives sets of even numbers $\{0,2,4,..., \Pi_{i=1}^n p_i -2 \}$. Furthermore, by generalizing our approach, the Polignac conjecture "the existence of infinitely many cousin primes, $p_{n+1}-p_{n}=4$," and the statement that "every even integer can be expressed as a difference of two primes," are derived as well.


title: Contraction based stabilization of nonlinear singularly perturbed systems and application to high gain feedback

abstract: Recent development of contraction theory based analysis of singularly perturbed system has opened the door for inspecting differential behavior of multi time-scale systems. In this paper a contraction theory based framework is proposed for stabilization of singularly perturbed systems. The primary objective is to design a feedback controller to achieve bounded tracking error for both standard and non-standard singularly perturbed systems. This framework provides relaxation over traditional quadratic Lyapunov based method as there is no need to satisfy interconnection conditions during controller design algorithm. Moreover, the stability bound does not depend on smallness of singularly perturbed parameter. Combined with high gain scaling, the proposed technique is shown to assure contraction of approximate feedback linearizable systems. These findings extend the class of nonlinear systems which can be made contracting.


title: Transition metals doped CuAlSe2 for promising intermediate band materials

abstract: Introducing an isolated intermediate band (IB) into a wide band gap semiconductor can potentially improve the optical absorption of the material beyond the Shockley-Queisser limitation for solar cells. Here, we present a systematic study of the thermodynamic stability, electronic structures, and optical properties of transition metals (M = Ti, V, and Fe) doped CuAlSe2 for potential IB thin film solar cells, by adopting the first-principles calculation based on the hybrid functional method. We found from chemical potential analysis that for all dopants considered, the stable doped phase only exits when the Al atom is substituted. More importantly, with this substitution, the IB feature is determined by $3d$ electronic nature of M^{3+} ion, and the electronic configuration of 3d^1 can drive a optimum IB that possesses half-filled character and suitable subbandgap from valence band or conduction band. We further show that Ti-doped CuAlSe2 is the more promising candidate for IB materials since the resulted IB in it is half filled and extra absorption peaks occurs in the optical spectrum accompanied with a largely enhanced light absorption intensity. The result offers a understanding for IB induced by transition metals into CuAlSe2 and is significant to fabricate the related IB materials.


title: A note on $n$-centralizer finite rings

abstract: Let $R$ be a finite ring and let $\Cent(R)$ denote the set of all distinct centralizers of $R$. $R$ is called an $n$-centralizer ring if $|Cent(R)| = n$. In this paper, we characterize $n$-centralizer finite rings for $n \leq 7$.


title: Single crystal 27Al-NMR study of the cubic Gamma3 ground doublet system PrTi2Al20

abstract: We report results of 27Al-NMR measurements on a single crystal of PrTi2Al20, in which the ground state of Pr3+ ions in the crystalline electric field is the nonmagnetic Gamma3 doublet. From the analysis of NMR spectra with the magnetic field applied precisely along the <111> and <100> directions, we determined the electric field gradient tensors for all three inequivalent Al sites (Al(1) ~ Al(3) sites) and the anisotropic Knight shifts for the Al(3) sites. The hyperfine coupling tensor at the Al(3) sites is strongly anisotropic and much larger than the classical dipolar coupling, indication importance of the anisotropic hybridization between the conduction and f electron states.


title: Interplay of fission modes in mass distribution of light actinide nuclei 225,227Pa

abstract: Fission-fragment mass distributions were measured for 225,227Pa nuclei formed in fusion reactions of 19F + 206, 208Pb around fusion barrier energies. Mass-angle correlations do not indicate any quasi-fission like events in this bombarding energy range. Mass distributions were fitted by Gaussian distribution and mass variance extracted. At below-barrier energies, the mass variance was found to increase with decrease in energy for both nuclei. Results from present work were compared with existing data for induced fission of 224, 226Th and 228U around barrier energies. Enhancement in mass variance of 225, 227Pa nuclei at below-barrier energies shows evidence for presence of asymmetric fission events mixed with symmetric fission events. This is in agreement with the results of mass distributions of nearby nuclei 224, 226Th and 228U where two-mode fission process was observed. Two-mode feature of fission arises due to the shell effects changing the landscape of the potential energy surfaces at low excitation energies. The excitation-energy dependence of the mass variance gives strong evidence for survival of microscopic shell effects in fission of light actinide nuclei 225, 227Pa with initial excitation energy ~30 - 50 MeV.


title: Sequential Bayesian Model Selection of Regular Vine Copulas

abstract: Regular vine copulas can describe a wider array of dependency patterns than the multivariate Gaussian copula or the multivariate Student's t copula. This paper presents two contributions related to model selection of regular vine copulas. First, our pair copula family selection procedure extends existing Bayesian family selection methods by allowing pair families to be chosen from an arbitrary set of candidate families. Second, our method represents the first Bayesian model selection approach to include the regular vine density construction in its scope of inference. The merits of our approach are established in a simulation study that benchmarks against methods suggested in current literature. A real data example about forecasting of portfolio asset returns for risk measurement and investment allocation illustrates the viability and relevance of the proposed scheme.


title: A Study on Artificial Intelligence IQ and Standard Intelligent Model

abstract: Currently, potential threats of artificial intelligence (AI) to human have triggered a large controversy in society, behind which, the nature of the issue is whether the artificial intelligence (AI) system can be evaluated quantitatively. This article analyzes and evaluates the challenges that the AI development level is facing, and proposes that the evaluation methods for the human intelligence test and the AI system are not uniform; and the key reason for which is that none of the models can uniformly describe the AI system and the beings like human. Aiming at this problem, a standard intelligent system model is established in this study to describe the AI system and the beings like human uniformly. Based on the model, the article makes an abstract mathematical description, and builds the standard intelligent machine mathematical model; expands the Von Neumann architecture and proposes the Liufeng - Shiyong architecture; gives the definition of the artificial intelligence IQ, and establishes the artificial intelligence scale and the evaluation method; conduct the test on 50 search engines and three human subjects at different ages across the world, and finally obtains the ranking of the absolute IQ and deviation IQ ranking for artificial intelligence IQ 2014.


title: Coqblin-Schrieffer Model for an Ultra-cold Gas of Ytterbium atoms with Metastable States

abstract: Motivated by the impressive recent advance in manipulating cold ytterbium atoms we explore and substantiate the feasibility of realizing the Coqblin-Schrieffer model in a gas of cold fermionic $^{173}$Yb atoms. Making use of different AC polarizabillity of the electronic ground state (electronic configuration $^1S_0$) and the long lived metastable state (electronic configuration $^3P_0$), it is substantiated that the latter can be localized and serve as a magnetic impurity while the former remains itinerant. The exchange mechanism between the itinerant $^1S_0$ and the localized $^3P_0$ atoms is analyzed and shown to be antiferromagnetic. The ensuing SU(6) symmetric Coqblin-Schrieffer Hamiltonian is constructed, and, using the calculated exchange constant $J$, perturbative RG analysis yield the Kondo temperature $T_K$ that is experimentally accessible. A number of thermodynamic measurable observables are calculated in the weak coupling regime $T>T_K$ (using perturbative RG analysis) and in the strong coupling regime $T<T_K$ (employing known Bethe ansatz techniques).


title: A bidirectional dual-comb ring laser for simple and robust dual-comb spectroscopy

abstract: Fourier-transform spectroscopy is an indispensable tool for analyzing chemical samples in scientific research as well as chemical and pharmaceutical industries. Recently, its measurement speed, sensitivity, and precision have been shown to be significantly enhanced by using dual frequency combs. However, wide acceptance of this technique is hindered by its requirement for two frequency combs and active stabilization of the combs. Here we overcome this predicament with a Kerr-lens mode-locked bidirectional ring laser that generates two frequency combs with slightly different pulse repetition rates and a tunable yet highly stable rate difference. This peculiar lasing principle builds on a slight difference in optical cavity length between two counter-propagating lasing modes due to Kerr lensing. Since these combs are produced by the one and same laser cavity, their relative coherence stays passively stable without the need for active stabilization. To show its utility, we demonstrate broadband dual-comb spectroscopy with the single laser.


title: Near-Infrared Polarization Source Catalog of the Northeastern Regions of the Large Magellanic Cloud

abstract: We present a near-infrared band-merged photometric and polarimetric catalog for the 39$\arcmin$ $\times$ 69$\arcmin$ fields on the northeastern part of the Large Magellanic Cloud (LMC), which were observed using SIRPOL, an imaging polarimeter of the InfraRed Survey Facility (IRSF). This catalog lists 1,858 sources brighter than 14 mag at $H$ band with polarization signal-to-noise ratio greater than three in the $J$, $H$, or $K_s$ bands. Based on the relationship between the extinction and the polarization degree, we argue that the polarization mostly arises from dichroic extinctions caused by local interstellar dust in the LMC. This catalog allows us to map polarization structures to examine the global geometry of the local magnetic field, and to show a statistical analysis of polarization of each field to understand its polarization properties. At the selected fields with coherent polarization position angles, we estimate magnetic field strengths in the range of 3$-$25 $\mu$G using the Chandrasekhar-Fermi method. This implies the presence of large-scale magnetic fields on a scale of around one hundred parsecs. When comparing mid and far-infrared dust emission maps, we confirmed that the polarization patterns are well aligned with molecular clouds around the star-forming regions.


title: Pion-induced production of the $Z_c(3900)$ off a nuclear target

abstract: We investigate the possibility to study the charmoniumlike state $Z_c(3900)$ through the pion-induced production off a nuclear target. By using a high-energy pion beam, the $Z_c(3900)$ can be produced off a proton or nucleus though the Primakoff effect. The production amplitude is calculated in an effective Lagrangian approach combined with the vector dominance model. The total cross sections of the $p(\pi^-, Z^-_c(3900))$ and $p(\pi^-, Z^-_c(3900)\to J/\psi\pi^-)$ reactions are calculated, and their order of magnitude is about 0.1 and 0.01nb, respectively, with an assumption of branch ratio 10\% for the $Z_c(3900)$ decay in $J/\psi\pi$ channel. If the proton target is replaced by a nuclear target, the production of the $Z_c(3900)$ enhances obviously. The predicted total cross sections for the $A(\pi^-, Z^-_c(3900))$ and $A(\pi^-, Z^-_c(3900)\to J/\psi\pi^-)$ reactions with $A=^{12}$C or $^{208}$Pb are on the order of magnitude of 100 and 10 nb, respectively, which is about one thousand times larger than the cross sections off a proton target. Based on these the results, we suggest the experimental study of the $Z_c(3900)$ by using high-energy pion beams with a nuclear target at facilities such as COMPASS and J-PARC.


title: Bayesian non-parametric inference for $\Lambda$-coalescents: consistency and a parametric method

abstract: We investigate Bayesian non-parametric inference of the $\Lambda$-measure of $\Lambda$-coalescent processes with recurrent mutation, parametrised by probability measures on the unit interval. We give verifiable criteria on the prior for posterior consistency when observations form a time series, and prove that any non-trivial prior is inconsistent when all observations are contemporaneous. We then show that the likelihood given a data set of size $n \in \mathbb{N}$ is constant across $\Lambda$-measures whose leading $n - 2$ moments agree, and focus on inferring truncated sequences of moments. We provide a large class of functionals which can be extremised using finite computation given a credible region of posterior truncated moment sequences, and a pseudo-marginal Metropolis-Hastings algorithm for sampling the posterior. Finally, we compare the efficiency of the exact and noisy pseudo-marginal algorithms with and without delayed acceptance acceleration using a simulation study.

categories: stat.ME math.PR math.ST q-bio.PE stat.CO stat.TH


title: Cavity quantum electrodynamics with ferromagnetic magnons in a small yttrium-iron-garnet sphere

abstract: Hybridizing collective spin excitations and a cavity with high cooperativity provides a new research subject in the field of cavity quantum electrodynamics and can also have potential applications to quantum information. Here we report an experimental study of cavity quantum electrodynamics with ferromagnetic magnons in a small yttrium-iron-garnet (YIG) sphere at both cryogenic and room temperatures. We observe for the first time a strong coupling of the same cavity mode to both a ferromagnetic-resonance (FMR) mode and a magnetostatic (MS) mode near FMR in the quantum limit. This is achieved at a temperature ~ 22 mK, where the average microwave photon number in the cavity is less than one. At room temperature, we also observe strong coupling of the cavity mode to the FMR mode in the same YIG sphere and find a slight increase of the damping rate of the FMR mode. These observations reveal the extraordinary robustness of the FMR mode against temperature. However, the MS mode becomes unobservable at room temperature in the measured transmission spectrum of the microwave cavity containing the YIG sphere. Our numerical simulations show that this is due to a drastic increase of the damping rate of the MS mode.


title: Fast Low-Rank Matrix Learning with Nonconvex Regularization

abstract: Low-rank modeling has a lot of important applications in machine learning, computer vision and social network analysis. While the matrix rank is often approximated by the convex nuclear norm, the use of nonconvex low-rank regularizers has demonstrated better recovery performance. However, the resultant optimization problem is much more challenging. A very recent state-of-the-art is based on the proximal gradient algorithm. However, it requires an expensive full SVD in each proximal step. In this paper, we show that for many commonly-used nonconvex low-rank regularizers, a cutoff can be derived to automatically threshold the singular values obtained from the proximal operator. This allows the use of power method to approximate the SVD efficiently. Besides, the proximal operator can be reduced to that of a much smaller matrix projected onto this leading subspace. Convergence, with a rate of O(1/T) where T is the number of iterations, can be guaranteed. Extensive experiments are performed on matrix completion and robust principal component analysis. The proposed method achieves significant speedup over the state-of-the-art. Moreover, the matrix solution obtained is more accurate and has a lower rank than that of the traditional nuclear norm regularizer.


title: Nanostrucure of Hybrid Plasmonic-Potonic Crystal Formed on Gel-Immobilized Colloidal Crystal Observer by AFM after Drying

abstract: Aiming at fabrication of hybrid plasmonic-photonic crystals, gel-immobilized colloidal crystals made of a polystyrene colloidal suspension and an N-(Hydroxy methyl)acrylamid-based gel were immersed into an aqueous dispersion of gold nanoparticles. Atomic force microscope (AFM) observations have been performed for the gel-immobilized colloidal crystals with gold nanoparticles deposited on their surfaces. In the present study, the diameter of a colloidal sphere was c.a. 190 nm. The diameter of a gold nanoparticle was the same as in a preliminary study, c.a. 40 nm. Various immersion times up to two hours were tested. Surface of a sample of 2 hr immersion has been observed. Prior to the AFM observation, the sample was dried in a desiccator for 18 hrs. We have identified a face-centered cubic {111} structure of a colloidal crystal of nearly close packing. Nanoparticles isolated with one another have been observed on the surface of gel-immobilized crystal, which can be regarded as gold nanoparticles from their sizes.


title: Flux jumps at pulsed field magnetization

abstract: Flux jumps lead to 10-20 times reduction of the shielding currents at a pulsed field magnetization of a melt-grown HTS annuli. The induced circular currents and the related trapped field are small after a field pulse since the duration of the pulse is shorter than the temperature relaxation time (<<1 s). Moreover, the residual trapped field in the hole of the annuli has the opposite direction in respect to the applied field. It is shown that the small value of the induced circular currents at the relatively high average critical current density is related to a narrow angular region of the flux motion (approximately 7 deg).


title: Star formation towards the southern Cometary HII region IRAS 17256-3631

abstract: IRAS 17256-3631 is a southern Galactic massive star forming region located at a distance of 2 kpc. In this paper, we present a multiwavelength investigation of the embedded cluster, the HII region, as well as the parent cloud. Radio images at 325, 610 and 1372 MHz were obtained using GMRT, India while the near-infrared imaging and spectroscopy were carried out using UKIRT and Mt. Abu Infrared Telescope, India. The near-infrared K-band image reveals the presence of a partially embedded infrared cluster. The spectral features of the brightest star in the cluster, IRS-1, spectroscopically agrees with a late O or early B star and could be the driving source of this region. Filamentary H_2 emission detected towards the outer envelope indicates presence of highly excited gas. The parent cloud is investigated at far-infrared to millimeter wavelengths and eighteen dust clumps have been identified. The spectral energy distributions (SEDs) of these clumps have been fitted as modified blackbodies and the best-fit peak temperatures are found to range from 14-33 K, while the column densities vary from 0.7-8.5x10^22 cm^-2. The radio maps show a cometary morphology for the distribution of ionized gas that is density bounded towards the north-west and ionization bounded towards the south-east. This morphology is better explained with the champagne flow model as compared to the bow shock model. Using observations at near, mid and far-infrared, submillimeter and radio wavelengths, we examine the evolutionary stages of various clumps.