title: Non-Axisymmetric Structure in the Satellite Dwarf Galaxy NGC2976: Implications for its Dark/Bright Mass Distribution and Evolution

abstract: We present the result of an extensive search for non-axisymmetric structures in the dwarf satellite galaxy of M81: NGC 2976, using multiwavelength archival observations. The galaxy is known to present kinematic evidence for a bysimmetric distortion, anyhow the stellar bar presence is controversial, due to the possible interpretation of NGC 2976 presenting an elliptical disk triggered by a prolate dark matter halo. We applied diagnostics in order to detect stellar bars or spiral arms. The m=2 fourier phase has a jump around 60 arcsecs consistent with a central bar and bisymmetric arms. The CO, 3.6 $\mu$ surface brightness and the dust lanes are consistent with a gas rich central bar and possibly with gaseous spiral arms. The kinematic jumps related with the dust lanes suggest that the bar perturbation in the disk kinematics is non-negligible and the reported non-circular motions, the central gas excess and the nuclear X-ray source (AGN/Starburst) might be produced by the central bar. SPH simulations of disks inside triaxial dark halos suggest that the two symmetric spots at 130 arcsecs and the narrow arms may be produced by gas at turning points in an elliptical disk, alternatively the potential ellipticity can be originated by tidally induced strong stellar bar/arms. We conclude that, the gas rich bar and spiral arms triggered by tidal stirring, and primordial halo triaxiality, can explain most of the non-circular motions, mass redistribution and nuclear activity. Distinguishing between them requires detailed modeling of environmental effects. An analysis similar to this may reveal such kind of structures in other nearby dwarf satellite galaxies, and if confirmed, the same evolutionary scenario will be applicable to them. This implies biases constraining their dark matter distribution and also making comparison against theoretical predictions for isolated galaxies.


title: The Peculiar Pulsar Population of the Central Parsec

abstract: Pulsars orbiting the Galactic center black hole, Sgr A*, would be potential probes of its mass, distance and spin, and may even be used to test general relativity. Despite predictions of large populations of both ordinary and millisecond pulsars in the Galactic center, none have been detected within 25 pc by deep radio surveys. One explanation has been that hyperstrong temporal scattering prevents pulsar detections, but the recent discovery of radio pulsations from a highly magnetized neutron star (magnetar) within 0.1 pc shows that the temporal scattering is much weaker than predicted. We argue that an intrinsic deficit in the ordinary pulsar population is the most likely reason for the lack of detections to date: a "missing pulsar problem" in the Galactic center. In contrast, we show that the discovery of a single magnetar implies efficient magnetar formation in the region. If the massive stars in the central parsec form magnetars rather than ordinary pulsars, their short lifetimes could explain the missing pulsars. Efficient magnetar formation could be caused by strongly magnetized progenitors, or could be further evidence of a top-heavy initial mass function. Furthermore, current high-frequency surveys should already be able to detect bright millisecond pulsars, given the measured degree of temporal scattering.


title: CMB Lensing Power Spectrum Biases from Galaxies and Clusters using High-angular Resolution Temperature Maps

abstract: The lensing power spectrum from cosmic microwave background (CMB) temperature maps will be measured with unprecedented precision with upcoming experiments, including upgrades to ACT and SPT. Achieving significant improvements in cosmological parameter constraints, such as percent level errors on sigma_8 and an uncertainty on the total neutrino mass of approximately 50 meV, requires percent level measurements of the CMB lensing power. This necessitates tight control of systematic biases. We study several types of biases to the temperature-based lensing reconstruction signal from foreground sources such as radio and infrared galaxies and the thermal Sunyaev-Zel'dovich effect from galaxy clusters. These foregrounds bias the CMB lensing signal due to their non-Gaussian nature. Using simulations as well as some analytical models we find that these sources can substantially impact the measured signal if left untreated. However, these biases can be brought to the percent level if one masks galaxies with fluxes at 150 GHz above 1 mJy and galaxy clusters with masses above M_vir = 10^14 M_sun. To achieve such percent level bias, we find that only modes up to a maximum multipole of l_max ~ 2500 should be included in the lensing reconstruction. We also discuss ways to minimize additional bias induced by such aggressive foreground masking by, for example, exploring a two-step masking and in-painting algorithm.


title: Anomaly inflow and thermal equilibrium

abstract: Using the anomaly inflow mechanism, we compute the flavor/Lorentz non-invariant contribution to the partition function in a background with a U(1) isometry. This contribution is a local functional of the background fields. By identifying the U(1) isometry with Euclidean time we obtain a contribution of the anomaly to the thermodynamic partition function from which hydrostatic correlators can be efficiently computed. Our result is in line with, and an extension of, previous studies on the role of anomalies in a hydrodynamic setting. Along the way we find simplified expressions for Bardeen-Zumino polynomials and various transgression formulae


title: TYC 3159-6-1: a runaway blue supergiant

abstract: We report the results of optical spectroscopy of a candidate evolved massive star in the Cygnus X region, TYC 3159-6-1, revealed via detection of its curious circumstellar nebula in archival data of the Spitzer Space Telecope. We classify TYC 3159-6-1 as an O9.5-O9.7 Ib star and derive its fundamental parameters by using the stellar atmosphere code FASTWIND. The He and CNO abundances in the photosphere of TYC 3159-6-1 are consistent with the solar abundances, suggesting that the star only recently evolved off the main sequence. Proper motion and radial velocity measurements for TYC 3159-6-1 show that it is a runaway star. We propose that Dolidze 7 is its parent cluster. We discuss the origin of the nebula around TYC 3159-6-1 and suggest that it might be produced in several successive episodes of enhanced mass-loss rate (outbursts) caused by rotation of the star near the critical, \Omega-limit.


title: Modeling Self-Subtraction in Angular Differential Imaging: Application to the HD 32297 Debris Disk

abstract: We present a new technique for forward-modeling self-subtraction of spatially extended emission in observations processed with angular differential imaging (ADI) algorithms. High-contrast direct imaging of circumstellar disks is limited by quasi-static speckle noise and ADI is commonly used to suppress those speckles. However, the application of ADI can result in self-subtraction of the disk signal due to the disk's finite spatial extent. This signal attenuation varies with radial separation and biases measurements of the disk's surface brightness, thereby compromising inferences regarding the physical processes responsible for the dust distribution. To compensate for this attenuation, we forward-model the disk structure and compute the form of the self-subtraction function at each separation. As a proof of concept, we apply our method to 1.6 and 2.2 micron Keck AO NIRC2 scattered-light observations of the HD 32297 debris disk reduced using a variant of the "locally optimized combination of images" (LOCI) algorithm. We are able to recover disk surface brightness that was otherwise lost to self-subtraction and produce simplified models of the brightness distribution as it appears with and without self-subtraction. From the latter models, we extract radial profiles for the disk's brightness, width, midplane position, and color that are unbiased by self-subtraction. Our analysis of these measurements indicates a break in the brightness profile power law at r~110 AU and a disk width that increases with separation from the star. We also verify disk curvature that displaces the midplane by up to 30 AU towards the northwest relative to a straight fiducial midplane.


title: From cosmic ray source to the Galactic pool

abstract: The Galactic cosmic ray spectrum is a remarkably straight power law. Our current understanding is that the dominant sources that accelerate cosmic rays up to the knee ($3 \times 10^{15}$ eV) or perhaps even the ankle ($3 \times 10^{18}$ eV), are young Galactic supernova remnants. In theory, however, there are various reasons why the spectrum may be different for different sources, and may not even be a power law if nonlinear shock acceleration applies during the most efficient stages of acceleration. We show how the spectrum at the accelerator translates to the spectrum that make up the escaping cosmic rays that replenish the Galactic pool of cosmic rays. We assume that cosmic ray confinement, and thus escape, is linked to the level of magnetic field amplification, and that the magnetic field is amplified by streaming cosmic rays according to the non-resonant hybrid or resonant instability. When a fixed fraction of the energy is transferred to cosmic rays, it turns out that a source spectrum that is flatter than $E^{-2}$ will result in a $E^{-2}$ escape spectrum, whereas a steeper source spectrum will result in an escape spectrum with equal steepening. This alleviates some of the concern that may arise from expected flat or concave cosmic ray spectra associated with nonlinear shock modification.


title: Multiplicativity of completely bounded $p$-norms implies a strong converse for entanglement-assisted capacity

abstract: The fully quantum reverse Shannon theorem establishes the optimal rate of noiseless classical communication required for simulating the action of many instances of a noisy quantum channel on an arbitrary input state, while also allowing for an arbitrary amount of shared entanglement of an arbitrary form. Turning this theorem around establishes a strong converse for the entanglement-assisted classical capacity of any quantum channel. This paper proves the strong converse for entanglement-assisted capacity by a completely different approach and identifies a bound on the strong converse exponent for this task. Namely, we exploit the recent entanglement-assisted "meta-converse" theorem of Matthews and Wehner, several properties of the recently established sandwiched Renyi relative entropy (also referred to as the quantum Renyi divergence), and the multiplicativity of completely bounded $p$-norms due to Devetak et al. The proof here demonstrates the extent to which the Arimoto approach can be helpful in proving strong converse theorems, it provides an operational relevance for the multiplicativity result of Devetak et al., and it adds to the growing body of evidence that the sandwiched Renyi relative entropy is the correct quantum generalization of the classical concept for all $\alpha>1$.


title: Charming the Higgs

abstract: We show that current Higgs data permit a significantly enhanced Higgs coupling to charm pairs, comparable to the Higgs to bottom pairs coupling in the Standard Model, without resorting to additional new physics sources in Higgs production. With a mild level of the latter current data even allow for the Higgs to charm pairs to be the dominant decay channel. An immediate consequence of such a large charm coupling is a significant reduction of the Higgs signal strengths into the known final states as in particular into bottom pairs. This might reduce the visible vector-boson associated Higgs production rate to a level that could compromise the prospects of ever observing it. We however demonstrate that a significant fraction of this reduced signal can be recovered by jet-flavor-tagging targeted towards charm-flavored jets. Finally we argue that an enhanced Higgs to charm pairs coupling can be obtained in various new physics scenarios in the presence of only a mild accidental cancellation between various contributions.


title: An advanced scattered moonlight model for Cerro Paranal

abstract: The largest natural source of light at night is the Moon, and it is the major contributor to the astronomical sky background. Being able to accurately predict the sky background, including scattered moonlight is important for scheduling astronomical observations. We have developed an improved scattered moonlight model, in which the components are computed with a better physical understanding as opposed to the simple empirical fit in the frequently used photometric model of Krisciunas & Schaefer (1991). Our spectroscopic model can better trace the spectral trends of scattered moonlight for any position of the Moon and target observation. This is the first scattered moonlight model that we know of which is this physical and versatile. We have incorporated an observed solar spectrum, accurate lunar albedo fit, and elaborate scattering and absorption calculations that include scattering off of molecules and aerosols. It was designed for Cerro Paranal, but can be modified for any location with known atmospheric properties. Throughout the optical range, the uncertainty is less than 20%. This advanced scattered moonlight model can predict the amount of scattered moonlight for any given geometry of the Moon and target, and lunar phase for the entire optical spectrum.


title: What Next-Generation 21 cm Power Spectrum Measurements Can Teach Us About the Epoch of Reionization

abstract: A number of experiments are currently working towards a measurement of the 21 cm signal from the Epoch of Reionization. Whether or not these experiments deliver a detection of cosmological emission, their limited sensitivity will prevent them from providing detailed information about the astrophysics of reionization. In this work, we consider what types of measurements will be enabled by a next-generation of larger 21 cm EoR telescopes. To calculate the type of constraints that will be possible with such arrays, we use simple models for the instrument, foreground emission, and the reionization history. We focus primarily on an instrument modeled after the $\sim 0.1~\rm{km}^2$ collecting area Hydrogen Epoch of Reionization Array (HERA) concept design, and parameterize the uncertainties with regard to foreground emission by considering different limits to the recently described "wedge" footprint in k-space. Uncertainties in the reionization history are accounted for using a series of simulations which vary the ionizing efficiency and minimum virial temperature of the galaxies responsible for reionization, as well as the mean free path of ionizing photons through the IGM. Given various combinations of models, we consider the significance of the possible power spectrum detections, the ability to trace the power spectrum evolution versus redshift, the detectability of salient power spectrum features, and the achievable level of quantitative constraints on astrophysical parameters. Ultimately, we find that $0.1~\rm{km}^2$ of collecting area is enough to ensure a very high significance ($\gtrsim30\sigma$) detection of the reionization power spectrum in even the most pessimistic scenarios. This sensitivity should allow for meaningful constraints on the reionization history and astrophysical parameters, especially if foreground subtraction techniques can be improved and successfully implemented.


title: V4641 Sgr: a candidate precessing microblazar

abstract: The X-ray spectrum of the Galactic X-ray binary V4641 Sgr in outburst has been found to exhibit a remarkably broad emission feature above 4 keV, with inferred equivalent widths up to 2 keV. Such a feature was first detected during the X-ray flaring activity associated with the giant outburst that the source experienced in 1999 September. The extraordinarily large equivalent width line was then ascribed to reflection/reprocessing of fluorescent Fe emission within an extended optically thick outflow enshrouding the binary system as a result of a short-lived, super-Eddington accretion episode. Making use of new and archival X-ray observations, we show here that a similar feature persists over four orders of magnitude in luminosity, down to Eddington ratios as low as log(L_Edd) = -4.5, where the existence of an optically thick envelope appears at odds with any viable accretion flow model. Possible interpretations for this highly unusual X-ray spectrum include a blend of Doppler shifted/boosted Fe lines from unresolved X-ray jets (a la SS433), or, the first Galactic analog of a blazar spectrum, where the >4 keV emission would correspond to the onset of the Inverse Compton hump. Either requires a low inclination angle of the jet with respect to the line of sight, in agreement with the estimates for the 1999 superluminal jet (i_jet<10 deg). The fast variability of the feature, combined with the high orbital axis inclination (60 deg< i_orb<71 deg), argue for a rapidly precessing accretion flow around V4641 Sgr, possibly leading to a transient microblazar behavior.


title: A feasible roadmap for unsupervised deconvolution of two-source mixed gene expressions

abstract: Tissue heterogeneity is a major confounding factor in studying individual populations that cannot be resolved directly by global profiling. Experimental solutions to mitigate tissue heterogeneity are expensive, time consuming, inapplicable to existing data, and may alter the original gene expression patterns. Here we ask whether it is possible to deconvolute two-source mixed expressions (estimating both proportions and cell-specific profiles) from two or more heterogeneous samples without requiring any prior knowledge. Supported by a well-grounded mathematical framework, we argue that both constituent proportions and cell-specific expressions can be estimated in a completely unsupervised mode when cell-specific marker genes exist, which do not have to be known a priori, for each of constituent cell types. We demonstrate the performance of unsupervised deconvolution on both simulation and real gene expression data, together with perspective discussions.


title: Calibration of semi-analytic models of galaxy formation using Particle Swarm Optimization

abstract: We present a fast and accurate method to select an optimal set of parameters in semi-analytic models of galaxy formation and evolution (SAMs). Our approach compares the results of a model against a set of observables applying a stochastic technique called Particle Swarm Optimization (PSO), a self-learning algorithm for localizing regions of maximum likelihood in multidimensional spaces that outperforms traditional sampling methods in terms of computational cost. We apply the PSO technique to the SAG semi-analytic model combined with merger trees extracted from a standard $\Lambda$CDM N-body simulation. The calibration is performed using a combination of observed galaxy properties as constraints, including the local stellar mass function and the black hole to bulge mass relation. We test the ability of the PSO algorithm to find the best set of free parameters of the model by comparing the results with those obtained using a MCMC exploration. Both methods find the same maximum likelihood region, however the PSO method requires one order of magnitude less evaluations. This new approach allows a fast estimation of the best-fitting parameter set in multidimensional spaces, providing a practical tool to test the consequences of including other astrophysical processes in SAMs.


title: Dusty supernovae running the thermodynamics of the matter reinserted within young and massive super stellar clusters

abstract: Following the observational and theoretical evidence that points at core collapse supernovae as major producers of dust, here we calculate the hydrodynamics of the matter reinserted within young and massive super stellar clusters under the assumption of gas and dust radiative cooling. The large supernova rate expected in massive clusters allows for a continuous replenishment of dust immersed in the high temperature thermalized reinserted matter and warrants a stationary presence of dust within the cluster volume during the type II supernova era. We first show that such a balance determines the range of dust to gas mass ratio and this the dust cooling law. We then search for the critical line that separates stationary cluster winds from the bimodal cases in the cluster mechanical luminosity (or cluster mass) vs cluster size parameter space. In the latter, strong radiative cooling reduces considerably the cluster wind mechanical energy output and affects particularly the cluster central regions, leading to frequent thermal instabilities that diminish the pressure and inhibit the exit of the reinserted matter. Instead matter accumulates there and is expected to eventually lead to gravitational instabilities and to further stellar formation with the matter reinserted by former massive stars. The main outcome of the calculations is that the critical line is almost two orders of magnitude or more, depending on the assumed value of the adiabatic wind terminal speed, lower than when only gas radiative cooling is applied. And thus, many massive clusters are predicted to enter the bimodal regime.


title: Characterization of superconducting RF breakdown by two-mode excitation

abstract: We show that thermal and magnetic contributions to the breakdown of superconductivity in radiofrequency (RF) fields can be separated by applying two RF modes simultaneously to a superconducting surface. We develop a simple model that illustrates how mode-mixing RF data can be related to properties of the superconductor. Within our model the data can be described by a single parameter, which can be derived either from RF or thermometry data. Our RF and thermometry data are in good agreement with the model. We propose to use mode-mixing technique to decouple thermal and magnetic effects on RF breakdown of superconductors.


title: Optical/Near-IR Polarization Survey of Sh 2-29: Magnetic Fields, Dense Cloud Fragmentations and Anomalous Dust Grain Sizes

abstract: Sh 2-29 is a conspicuous star-forming region marked by the presence of massive embedded stars as well as several notable interstellar structures. In this research, our goals were to determine the role of magnetic fields and to study the size distribution of interstellar dust particles within this turbulent environment. We have used a set of optical and near-infrared polarimetric data obtained at OPD/LNA (Brazil) and CTIO (Chile), correlated with extinction maps, 2MASS data and images from DSS and Spitzer. The region's most striking feature is a swept out interstellar cavity whose polarimetric maps indicate that magnetic field lines were dragged outwards, pilling up along its borders. This led to a higher magnetic strength value ($\approx400\,\mu$G) and an abrupt increase in polarization degree, probably due to an enhancement in alignment efficiency. Furthermore, dense cloud fragmentations with peak $A_{V}$ between 20 and 37 mag were probably triggered by its expansion. The presence of $24\,\mu$m point-like sources indicates possible newborn stars inside this dense environment. A statistical analysis of the angular dispersion function revealed areas where field lines are aligned in a well-ordered pattern, seemingly due to compression effects from the HII region expansion. Finally, Serkowski function fits were used to study the ratio of the total-to-selective extinction, reveling a dual population of anomalous grain particles sizes. This trend suggests that both effects of coagulation and fragmentation of interstellar grains are present in the region.


title: Evidence that All States Are Unitarily Equivalent to X States of the Same Entanglement

abstract: Strong numerical evidence is presented suggesting that all two-qubit mixed states are equivalent to X states by a single entanglement-preserving unitary (EPU) transformation, so that the concurrence of such an X state equals that of the original general state. An X-state parameterization of a general two-qubit state is given, allowing all states to have their concurrence parametrically specified. A new kind of entanglement measure is proposed, relating a general state's entanglement to that of a pure state in the same system. New states called "H States" are presented, having fully parametric concurrence and purity, with the intention of using them to construct entanglement-preserving depolarization channels, which may aid development of the new entanglement measure. A theory of "true-generalized" X states (TGX states) is proposed for the general case of $N$-partite systems. While such states do not generally have the literal "X" shape, evidence is shown that they are the true generalizations of X states in larger systems, since they appear to always be EPU-equivalent to general states of all ranks, whereas literal X states generally are not. An example of this is given for $2\times 3$, including the proposition of the $2\times 3$ maximally entangled mixed states (MEMS). If the claim that TGX states are universal is valid, then any entanglement measure may be computable in a simpler form by using the EPU-equivalence between general states and TGX states.


title: WIMP physics with ensembles of direct-detection experiments

abstract: The search for weakly-interacting massive particle (WIMP) dark matter is multi-pronged. Ultimately, the WIMP-dark-matter picture will only be confirmed if different classes of experiments see consistent signals and infer the same WIMP properties. In this work, we review the ideas, methods, and status of direct-detection searches. We focus in particular on extracting WIMP physics (WIMP interactions and phase-space distribution) from direct-detection data in the early discovery days when multiple experiments see of order dozens to hundreds of events. To demonstrate the essential complementarity of different direct-detection experiments in this context, we create mock data intended to represent the data from the near-future Generation 2 experiments. We consider both conventional supersymmetry-inspired benchmark points (with spin-independent and -dependent elastic cross sections just below current limits), as well as benchmark points for other classes of models (inelastic and effective-operator paradigms). We also investigate the effect on parameter estimation of loosening or dropping the assumptions about the local WIMP phase-space distribution. We arrive at two main conclusions. Firstly, teasing out WIMP physics with experiments depends critically on having a wide set of detector target materials, spanning a large range of target nuclear masses and spin-dependent sensitivity. It is also highly desirable to obtain data from low-threshold experiments. Secondly, a general reconstruction of the local WIMP velocity distribution, which will only be achieved if there are multiple experiments using different target materials, is critical to obtaining a robust and unbiased estimate of the WIMP mass.


title: Cosmic Frontier Indirect Dark Matter Detection Working Group Summary

abstract: As part of the Snowmass process, the Cosmic Frontier Indirect-Detection subgroup (CF2) has drawn on input from the Cosmic Frontier and the broader Particle Physics community to produce this document. The purposes of this report are to identify opportunities for dark matter science through indirect detection, to give an overview of the primary scientific drivers for indirect searches for dark matter, and to survey current and planned experiments that have, as a large part of their scientific program, the goal of searching for indirect (or astrophysical) signatures of dark matter. We primarily address existing experiments with a large U.S. role, or future experiments where a U.S. contribution is sought. We also address the limitations of this technique, and answer the tough questions relevant to this subgroup posed by the HEP community through the Snowmass process.


title: A complete classification of equational classes of threshold functions included in clones

abstract: The class of threshold functions is known to be characterizable by functional equations or, equivalently, by pairs of relations, which are called relational constraints. It was shown by Hellerstein that this class cannot be characterized by a finite number of such objects. In this paper, we investigate classes of threshold functions which arise as intersections of the class of all threshold functions with clones of Boolean functions, and provide a complete classification of such intersections in respect to whether they have finite characterizations. Moreover, we provide a characterizing set of relational constraints for each class of threshold functions arising in this way.


title: On Convexification of Range Measurement Based Sensor and Source Localization Problems

abstract: This manuscript is a preliminary pre-print version of a journal submission by the authors, revisiting the problem of range measurement based localization of a signal source or a sensor. The major geometric difficulty of the problem comes from the non-convex structure of optimization tasks associated with range measurements, noting that the set of source locations corresponding to a certain distance measurement by a fixed point sensor is non-convex both in two and three dimensions. Differently from various recent approaches to this localization problem, all starting with a non-convex geometric minimization problem and attempting to devise methods to compensate the non-convexity effects, we suggest a geometric strategy to compose a convex minimization problem first, that is equivalent to the initial non-convex problem, at least in noise-free measurement cases. Once the convex equivalent problem is formed, a wide variety of convex minimization algorithms can be applied. The paper also suggests a gradient based localization algorithm utilizing the introduced convex cost function for localization. Furthermore, the effects of measurement noises are briefly discussed. The design, analysis, and discussions are supported by a set of numerical simulations.


title: Nucleon form factors with light Wilson quarks

abstract: We present nucleon observables - primarily isovector vector form factors - from calculations using 2+1 flavors of Wilson quarks. One ensemble is used for a dedicated high-precision study of excited-state effects using five source-sink separations between 0.7 and 1.6 fm. We also present results from a larger set of calculations that include an ensemble with pion mass 149 MeV and box size 5.6 fm, which nearly eliminates the uncertainty associated with extrapolation to the physical pion mass. The results show agreement with experiment for the vector form factors, which occurs only when excited-state contributions are reduced. Finally, we show results from a subset of ensembles that have pion mass 254 MeV with varying temporal and spatial box sizes, which we use for a controlled study of finite-volume effects and a test of the "$m_\pi L=4$" rule of thumb.


title: Latest RHIC Results on Ultra-Peripheral Collisions

abstract: With RHIC running in its second phase at higher luminosities the new data sets collected so far by PHENIX and STAR are allowing improvements in the study of vector meson photo-production in Ultra peripheral Collision events in Au+Au at the highest energy. This is a brief summary of what has been accomplished so far by both collaborations.


title: Twisted relative Cuntz-Krieger algebras associated to finitely aligned higher-rank graphs

abstract: To each finitely aligned higher-rank graph $\Lambda$ and each $\mathbb{T}$-valued 2-cocycle on $\Lambda$, we associate a family of twisted relative Cuntz-Krieger algebras. We show that each of these algebras carries a gauge action, and prove a gauge-invariant uniqueness theorem. We describe an isomorphism between the fixed point algebras for the gauge actions on the twisted and untwisted relative Cuntz-Krieger algebras. We show that the quotient of a twisted relative Cuntz-Krieger algebra by a gauge-invariant ideal is canonically isomorphic to a twisted relative Cuntz-Krieger algebra associated to a subgraph. We use this to provide a complete graph-theoretic description of the gauge-invariant ideal structure of each twisted relative Cuntz-Krieger algebra.


title: Broadband and efficient plasmonic control in the near-infrared and visible via strong interference of surface plasmon polaritons

abstract: Broadband and tunable control of surface plasmon polaritons in the near-infrared and visible spectrum is demonstrated theoretically and numerically with a pair of phased nanoslits. We establish with simulations supported by a coupled wave model that by dividing the incident power equally between two input channels, the maximum plasmon intensity deliverable to either side of the nanoslit pair is twice that for an isolated slit. For a broadband source, a compact device with nanoslit separation of the order of a tenth of the wavelength is shown to steer nearly all the generated plasmons to one side for the same phase delay, thereby achieving a broadband unidirectional plasmon launcher. The reported effect can be applied to the design of ultra-broadband and efficient tunable plasmonic devices.


title: Summary of the EDS Blois 2013 Workshop

abstract: I give a personal overview of some highlights of the EDS Blois 2013 Conference on elastic and diffractive scattering in Saariselka, Finland.


title: Scaling SVM and Least Absolute Deviations via Exact Data Reduction

abstract: The support vector machine (SVM) is a widely used method for classification. Although many efforts have been devoted to develop efficient solvers, it remains challenging to apply SVM to large-scale problems. A nice property of SVM is that the non-support vectors have no effect on the resulting classifier. Motivated by this observation, we present fast and efficient screening rules to discard non-support vectors by analyzing the dual problem of SVM via variational inequalities (DVI). As a result, the number of data instances to be entered into the optimization can be substantially reduced. Some appealing features of our screening method are: (1) DVI is safe in the sense that the vectors discarded by DVI are guaranteed to be non-support vectors; (2) the data set needs to be scanned only once to run the screening, whose computational cost is negligible compared to that of solving the SVM problem; (3) DVI is independent of the solvers and can be integrated with any existing efficient solvers. We also show that the DVI technique can be extended to detect non-support vectors in the least absolute deviations regression (LAD). To the best of our knowledge, there are currently no screening methods for LAD. We have evaluated DVI on both synthetic and real data sets. Experiments indicate that DVI significantly outperforms the existing state-of-the-art screening rules for SVM, and is very effective in discarding non-support vectors for LAD. The speedup gained by DVI rules can be up to two orders of magnitude.


title: The Small Numbers of Large Kuiper Belt Objects

abstract: We explore the brightness distribution of the largest and brightest (m(R)<22) Kuiper belt objects (KBOs). We construct a luminosity function of the dynamically excited or hot Kuiper belt (orbits with inclinations > 5 degrees) from the very brightest to m(R)=23. We find for m(R)< 23, a single slope appears to describe the luminosity function. We estimate ~12 KBOs brighter than m(R)~19.5 are present in the Kuiper belt today. With 9 bodies already discovered this suggests that the inventory of bright KBOs is nearly complete.


title: Rank-finiteness for modular categories

abstract: We prove a rank-finiteness conjecture for modular categories: up to equivalence, there are only finitely many modular categories of any fixed rank. Our technical advance is a generalization of the Cauchy theorem in group theory to the context of spherical fusion categories. For a modular category $\mathcal{C}$ with $N=ord(T)$, the order of the modular $T$-matrix, the Cauchy theorem says that the set of primes dividing the global quantum dimension $D^2$ in the Dedekind domain $\mathbb{Z}[e^{\frac{2\pi i}{N}}]$ is identical to that of $N$.


title: Nonlinear Fourier analysis for discontinuous conductivities: computational results

abstract: Two reconstruction methods of Electrical Impedance Tomography (EIT) are numerically compared for nonsmooth conductivities in the plane based on the use of complex geometrical optics (CGO) solutions to D-bar equations involving the global uniqueness proofs for Calder\'on problem exposed in [Nachman; Annals of Mathematics 143, 1996] and [Astala and P\"aiv\"arinta; Annals of Mathematics 163, 2006]: the Astala-P\"aiv\"arinta theory-based "low-pass transport matrix method" implemented in [Astala et al.; Inverse Problems and Imaging 5, 2011] and the "shortcut method" which considers ingredients of both theories. The latter method is formally similar to the Nachman theory-based regularized EIT reconstruction algorithm studied in [Knudsen, Lassas, Mueller and Siltanen; Inverse Problems and Imaging 3, 2009] and several references from there. New numerical results are presented using parallel computation with size parameters larger than ever, leading mainly to two conclusions as follows. First, both methods can approximate piecewise constant conductivities better and better as the cutoff frequency increases, and there seems to be a Gibbs-like phenomenon producing ringing artifacts. Second, the transport matrix method loses accuracy away from a (freely chosen) pivot point located outside of the object to be studied, whereas the shortcut method produces reconstructions with more uniform quality.


title: B and L at the SUSY Scale, Dark Matter and R-parity Violation

abstract: We present a simple theory where baryon and lepton numbers are spontaneously broken at the supersymmetry scale. In this context R-parity must be spontaneously broken but the theory still contains a stable field which can play the role of the cold dark matter of the Universe. We discuss the spectrum of the theory, the properties of the dark matter candidate and the predictions for direct detection experiments. This theory provides a concrete example of exotic supersymmetric signatures associated with having the simultaneous presence of R-parity violating and missing energy signals at the Large Hadron Collider.


title: Multiplicative Arithmetic Functions of Several Variables: A Survey

abstract: We survey general properties of multiplicative arithmetic functions of several variables and related convolutions, including the Dirichlet convolution and the unitary convolution. We introduce and investigate a new convolution, called gcd convolution. We define and study the convolutes of arithmetic functions of several variables, according to the different types of convolutions. We discuss the multiple Dirichlet series and Bell series and present certain arithmetic and asymptotic results of some special multiplicative functions arising from problems in number theory, group theory and combinatorics. We give a new proof to obtain the asymptotic density of the set of ordered $r$-tuples of positive integers with pairwise relatively prime components and consider a similar question related to unitary divisors.


title: The NEXT experiment. Towards phase I

abstract: Phase I of the NEXT-100 $0\nu\beta\beta$ experiment (NEW) is scheduled for data taking in 2015 at Laboratorio Subterr\'aneo de Canfranc in the Spanish Pyrenees. Thanks to the light proportional technique, NEW anticipates an outstanding energy resolution nearing the Fano factor in Xenon (0.5-1%FWHM@$Q_{\beta\beta,^{136}Xe}$), with a TPC-design that allows tracking and identification of the double end-blob feature of the $0\nu\beta\beta$ decay. When properly mastered, the combination of these two assets can suppress the irreducible $2\nu\beta\beta$ and (single-blob) $\gamma$ backgrounds from natural radioactivity to minute levels, of the order of $5\times{10^{-4}}$ ckky. Given our knowledge of the available phase-space as obtained from neutrino oscillation experiments, this feat will expectedly allow for a sensitivity to the effective electron neutrino mass of $m_{\beta\beta}\simeq 30$ meV for exposures at the 20 ton $\times$ year scale. Hence, ultimately, a full survey of the inverse hierarchy of the neutrino mass ordering appears to be within reach for a ton-scale experiment based on this technology. NEW, with 10 kg of Xenon 90%-enriched in $^{136}$Xe, sets an unprecedented scale for gaseous Xenon TPCs and will be an important milestone for its anticipated upgrades (100 kg and 1 ton). I briefly summarize the status of the NEXT experiment, from the main results obtained with $\sim 1$ kg prototypes that substantiate the concept, to the ongoing works for deploying its first phase.


title: Asymptotic distribution for the birthday problem with multiple coincidences, via an embedding of the collision process

abstract: We study the random variable B(c,n), which counts the number of balls that must be thrown into n equally-sized bins in order to obtain c collisions. The asymptotic expected value of B(1,n) is the well-known $\sqrt{n\pi/2}$ appearing in the solution to the birthday problem; the limit distribution and asymptotic moments of B(1,n) are also well known. We calculate the distribution and moments of B(c,n) asymptotically as n goes to infinity and c = O(n). Our main tools are an embedding of the collision process, realizing the process as a deterministic function of the standard Poisson process, and a central limit result by Renyi.


title: A note on using Bayes priors for Weibull distribution

abstract: In this note, the practical use of priors for Bayes estimators, of the two parameters of the Weibull reliability model, is discussed in a technological context. The meaning of the priors as expression of virtual data samples is analyzed. The implications of physics of failures are also highlighted. The whole analysis shows a rational way to convert a really available technological knowledge into prior information, effectively and without a long elicitation process. A large Monte Carlo study, on both complete and censored samples of small size, shows the good properties of some estimators, which can exploit such a new approach.


title: Extremal Eigenvalues and Eigenvectors of Deformed Wigner Matrices

abstract: We consider random matrices of the form $H = W + \lambda V$, $\lambda\in\mathbb{R}^+$, where $W$ is a real symmetric or complex Hermitian Wigner matrix of size $N$ and $V$ is a real bounded diagonal random matrix of size $N$ with i.i.d.\ entries that are independent of $W$. We assume subexponential decay for the matrix entries of $W$ and we choose $\lambda \sim 1$, so that the eigenvalues of $W$ and $\lambda V$ are typically of the same order. Further, we assume that the density of the entries of $V$ is supported on a single interval and is convex near the edges of its support. In this paper we prove that there is $\lambda_+\in\mathbb{R}^+$ such that the largest eigenvalues of $H$ are in the limit of large $N$ determined by the order statistics of $V$ for $\lambda>\lambda_+$. In particular, the largest eigenvalue of $H$ has a Weibull distribution in the limit $N\to\infty$ if $\lambda>\lambda_+$. Moreover, for $N$ sufficiently large, we show that the eigenvectors associated to the largest eigenvalues are partially localized for $\lambda>\lambda_+$, while they are completely delocalized for $\lambda<\lambda_+$. Similar results hold for the lowest eigenvalues.


title: Lack of BV Bounds for Approximate Solutions to the p-system with Large Data

abstract: We consider front tracking approximate solutions to the p-system of isentropic gas dynamics. At interaction times, the outgoing wave fronts have the same strength as in the exact solution of the Riemann problem, but some error is allowed in their speed. For large BV initial data, we construct examples showing that the total variation of these approximate solutions can become arbitrarily large, or even blow up in finite time. This happens even if the density of the gas remains uniformly positive.


title: A Determinantal Formula for Catalan Tableaux and TASEP Probabilities

abstract: We present a determinantal formula for the steady state probability of each state of the TASEP (Totally Asymmetric Simple Exclusion Process) with open boundaries, a 1D particle model that has been studied extensively and displays rich combinatorial structure. These steady state probabilities are computed by the enumeration of Catalan tableaux, which are certain Young diagrams filled with $\alpha$'s and $\beta$'s that satisfy some conditions on the rows and columns. We construct a bijection from the Catalan tableaux to weighted lattice paths on a Young diagram, and from this we enumerate the paths with a determinantal formula, building upon a formula of Narayana that counts unweighted lattice paths on a Young diagram. Finally, we provide a formula for the enumeration of Catalan tableaux that satisfy a given condition on the rows, which corresponds to the steady state probability that in the TASEP on a lattice with $n$ sites, precisely $k$ of the sites are occupied by particles. This formula is an $\alpha\ /\ \beta$ generalization of the Narayana numbers.


title: On the foreshocks of strong earthquakes in the light of catastrophe theory

abstract: It is useful to consider the earthquakes in terms of catastrophe theory. In the paper, we illustrate this statement by analysis foreshocks preceding the strong earthquakes. We focused on the so-called catastrophe flags, and on the triggers that cause the critical transitions. The analysis shows a sharp increase of foreshock activity 3 h before the mainshock. This result is reminiscent of the well-known activation of the ULF magnetic precursors of earthquakes. Furthermore, we found that the characteristic frequency of foreshock sequence decreases a few hours before the mainshock, which is consistent with a prediction of the catastrophe theory. Finally, the analysis testifies that the strong foreshocks seems may be triggers of mainshocks. The idea is that the surface waves propagating outwards from the foreshock return back to the vicinity of the epicenter after having made a complete revolution around the Earth and induce there the mainshock.


title: GeV - PeV Neutrino Production and Oscillation in hidden jets from GRBs

abstract: Long gamma-ray bursts have been widely associated with collapsing massive stars in the framework of collapsar model. High-energy neutrinos and photons can be produced in the internal shocks of middle relativistic jets from core-collapse supernova. Although photons can hardly escape, high-energy neutrinos could be the only signature when the jets are hidden. We show that using suitable parameters, high-energy neutrinos in GeV - PeV range can be produced in the hidden jet inside the collapsar, thus demonstrating that these objects are candidates to produce neutrinos with energies between 1 - 10 PeV which were observed with IceCube. On the other hand, due to matter effects, high-energy neutrinos may oscillate resonantly from one flavor to another before leaving the star. Using two (solar, atmospheric and accelerator parameters) and three neutrino mixing, we study the possibility of resonant oscillation for these neutrinos created in internal shocks. Also we compute the probabilities of neutrino oscillations in the matter at different distances along the jet (before leaving the star) and after in vacuum, on their path to Earth. Finally, neutrino flavor ratios on Earth are estimated.


title: Real-Time Planning with Primitives for Dynamic Walking over Uneven Terrain

abstract: We present an algorithm for receding-horizon motion planning using a finite family of motion primitives for underactuated dynamic walking over uneven terrain. The motion primitives are defined as virtual holonomic constraints, and the special structure of underactuated mechanical systems operating subject to virtual constraints is used to construct closed-form solutions and a special binary search tree that dramatically speed up motion planning. We propose a greedy depth-first search and discuss improvement using energy-based heuristics. The resulting algorithm can plan several footsteps ahead in a fraction of a second for both the compass-gait walker and a planar 7-Degree-of-freedom/five-link walker.


title: On the Convergence of Decentralized Gradient Descent

abstract: Consider the consensus problem of minimizing $f(x)=\sum_{i=1}^n f_i(x)$ where each $f_i$ is only known to one individual agent $i$ out of a connected network of $n$ agents. All the agents shall collaboratively solve this problem and obtain the solution subject to data exchanges restricted to between neighboring agents. Such algorithms avoid the need of a fusion center, offer better network load balance, and improve data privacy. We study the decentralized gradient descent method in which each agent $i$ updates its variable $x_{(i)}$, which is a local approximate to the unknown variable $x$, by combining the average of its neighbors' with the negative gradient step $-\alpha \nabla f_i(x_{(i)})$. The iteration is $$x_{(i)}(k+1) \gets \sum_{\text{neighbor} j \text{of} i} w_{ij} x_{(j)}(k) - \alpha \nabla f_i(x_{(i)}(k)),\quad\text{for each agent} i,$$ where the averaging coefficients form a symmetric doubly stochastic matrix $W=[w_{ij}] \in \mathbb{R}^{n \times n}$. We analyze the convergence of this iteration and derive its converge rate, assuming that each $f_i$ is proper closed convex and lower bounded, $\nabla f_i$ is Lipschitz continuous with constant $L_{f_i}$, and stepsize $\alpha$ is fixed. Provided that $\alpha < O(1/L_h)$ where $L_h=\max_i\{L_{f_i}\}$, the objective error at the averaged solution, $f(\frac{1}{n}\sum_i x_{(i)}(k))-f^*$, reduces at a speed of $O(1/k)$ until it reaches $O(\alpha)$. If $f_i$ are further (restricted) strongly convex, then both $\frac{1}{n}\sum_i x_{(i)}(k)$ and each $x_{(i)}(k)$ converge to the global minimizer $x^*$ at a linear rate until reaching an $O(\alpha)$-neighborhood of $x^*$. We also develop an iteration for decentralized basis pursuit and establish its linear convergence to an $O(\alpha)$-neighborhood of the true unknown sparse signal.


title: Decaying Domain Walls in an Extended Gravity Model and Cosmology

abstract: We investigate cosmological consequences of an extended gravity model which belongs to the same class studied by Accetta and Steinhardt in an extended inflationary scenario. But we do not worry about inflation in our model; instead, we focus on a topological object formed during cosmological phase transitions. Although domain walls appear during first-order phase transitions such as QCD transition, they decay at the end of the phase transition. Therefore the "domain wall problem" does not exist in the suitable range of pameters and, on the contrary, the "fragments" of walls may become seeds of dark matter. A possible connection to "oscillating universe" model offered by Morikawa et al. is also discussed.


title: On the Contou-Carrere Symbol for Surfaces

abstract: This is a preliminary report on the Contou-Carrere symbol for surfaces. It consists of two parts. In the first part, we recall technical results needed to define the symbol. The second part is where we compute all components of the Coutou-Carrere symbol for surfaces, using iterated integrals over membranes.


title: The q=-1 phenomenon via homology concentration

abstract: We introduce a homological approach to exhibiting instances of Stembridge's q=-1 phenomenon. This approach is shown to explain two important instances of the phenomenon, namely that of partitions whose Ferrers diagrams fit in a rectangle of fixed size and that of plane partitions fitting in a box of fixed size. A more general framework of invariant and coinvariant complexes with coefficients taken mod 2 is developed, and as a part of this story an analogous homological result for necklaces is conjectured.


title: Repulsive interactions of a lipid membrane with graphene in composite materials

abstract: The van der Waals interaction between a lipid membrane and a substrate covered by a graphene sheet is investigated using the Lifshitz theory. The reflection coefficients are obtained for a layered planar system submerged in water. The dielectric response properties of the involved materials are also specified and discussed. Our calculations show that a graphene covered substrate can repel the biological membrane in water. This is attributed to the significant changes in the response properties of the system due to the monolayer graphene. It is also found that the van der Waals interaction is mostly dominated by the presence of graphene, while the role of the particular substrate is secondary.

categories: cond-mat.mes-hall cond-mat.mtrl-sci cond-mat.soft quant-ph


title: Scissors resonance in the quasi-continuum of Th, Pa and U isotopes

abstract: The gamma-ray strength function in the quasi-continuum has been measured for 231-233Th, 232,233Pa and 237-239U using the Oslo method. All eight nuclei show a pronounced increase in gamma strength at omega_SR approx 2.4 MeV, which is interpreted as the low-energy M1 scissors resonance (SR). The total strength is found to be B_SR = 9-11 mu_N^2 when integrated over the 1 - 4 MeV gamma-energy region. The SR displays a double-hump structure that is theoretically not understood. Our results are compared with data from (gamma, gamma') experiments and theoretical sum-rule estimates for a nuclear rigid-body moment of inertia.


title: Survey of Period Variations of Superhumps in SU UMa-Type Dwarf Novae. V: The Fifth Year (2012-2013)

abstract: Continuing the project described by Kato et al. (2009a, arXiv:0905.1757), we collected times of superhump maxima for SU UMa-type dwarf novae mainly observed during the 2012-2013 season. We found three objects (V444 Peg, CSS J203937 and MASTER J212624) having strongly positive period derivatives despite the long orbital period (Porb). By using the period of growing stage (stage A) superhumps, we obtained mass ratios for six objects. We characterized nine new WZ Sge-type dwarf novae. We made a pilot survey of the decline rate of slowly fading part of SU UMa-type and WZ Sge-type outbursts. The decline time scale was found to generally follow the expected Porb^(1/4) dependence and WZ Sge-type outbursts also generally follow this trend. There are some objects which show slower decline rates, and we consider these objects good candidates for period bouncers. We also studied unusual behavior in some objects, including BK Lyn which made a transition from an ER UMa-type state to the novalike (standstill) state in 2013 and unusually frequent occurrence of superoutbursts in NY Ser and CR Boo. We applied least absolute shrinkage and selection operator (Lasso) power spectral analysis, which has been proven to be very effective in analyzing the Kepler data, to ground-based photometry of BK Lyn and detected the dramatic disappearance of the signal of negative superhumps in 2013. We suggested that the mass-transfer rates did not vary strongly between the ER UMa-type state and novalike state in BK Lyn, and this transition was less likely caused by a systematic variation of the mass-transfer rate.


title: Evidence for dipole nature of the low-energy $\gamma$ enhancement in $^{56}$Fe

abstract: The $\gamma$-ray strength function of $^{56}$Fe has been measured from proton-$\gamma$ coincidences for excitation energies up to $\approx 11$ MeV. The low-energy enhancement in the $\gamma$-ray strength function, which was first discovered in the ($^3$He,$\alpha\gamma$)$^{56}$Fe reaction, is confirmed with the ($p,p^\prime\gamma$)$^{56}$Fe experiment reported here. Angular distributions of the $\gamma$ rays give for the first time evidence that the enhancement is dominated by dipole transitions.


title: Application of optimization method to the $x^4$ model in the Tsallis nonextensive statistics

abstract: We study the effects of the environment described by the Tsallis nonextensive statistics on physical quantities using an optimization method in the case of small deviation from the Boltzmann-Gibbs statistics. The $x^4$ model is used and the density operator is restricted to be a gaussian form. The variational parameter is the frequency $\Omega$ of a particle in the optimization method. We obtain an approximate expression of free energy and of the expectation value of $\beta m \Omega^2 x^2 /2$, where $\beta$ is the inverse of the temperature and $m$ is the mass of a particle. Numerically, the optimized frequency is estimated and the expectation value of $\beta m \Omega^2 x^2 /2$ is calculated. The effects of the Tsallis nonextensive statistic for small deviation from the Boltzmann-Gibbs statistics are found: 1) the frequency modulation of a particle and 2) the variation of the expectation value of $\beta m \Omega^2 x^2 /2$ at high temperature.


title: Transitional $\gamma$ strength in Cd isotopes

abstract: The level densities and $\gamma$-ray strength functions of $^{105,106,111,112}$Cd have been extracted from particle-$\gamma$ coincidence data using the Oslo method. The level densities are in very good agreement with known levels at low excitation energy. The $\gamma$-ray strength functions display no strong enhancement for low $\gamma$ energies. However, more low-energy strength is apparent for $^{105,106}$Cd than for $^{111,112}$Cd. For $\gamma$ energies above $\approx$ 4 MeV, there is evidence for some extra strength, similar to what has been previously observed for the Sn isotopes. The origin of this extra strength is unclear; it might be due to $E1$ and $M1$ transitions originating from neutron skin oscillations or the spin-flip resonance, respectively.


title: Isomorphism classes of finite dimensional connected Hopf algebras in positive characteristic

abstract: We classify all finite-dimensional connected Hopf algebras with large abelian primitive spaces. We show that they are Hopf algebra extensions of restricted enveloping algebras of certain restricted Lie algebras. For any abelian matched pair associated with these extensions, we construct a cohomology group, which classifies all the extensions up to equivalence. Moreover, we present a $1$-$1$ correspondence between the isomorphism classes and a group quotient of the cohomology group deleting some exceptional points, where the group respects the automorphisms of the abelian matched pair and the exceptional points represent those restricted Lie algebra extensions.


title: The $\tau_2$ model and parafermions

abstract: Paul Fendley has recently found a "parafermionic" way to diagonalise a simple solvable hamiltonian associated with the chiral Potts model. Here we indicate how this method generalizes to the $\tau_2$ model with open boundaries and make some comments.


title: Detection of ultra high energy neutrinos by IceCube: Sterile neutrino scenario

abstract: The short-baseline neutrino oscillation experiments, the excess of radiation from the measurement of the cosmic microwave background radiation, the necessity of the nonbaryonic dark matter candidate and the depletion of the neutrino flux in IceCube all seem to hint at new physics beyond the standard model. An economical way to address these issues is to invoke the existence of sterile neutrinos. We present simple extensions of the standard model with additional three sterile neutrinos and discuss the corresponding PMNS like neutrino flavor mixing matrix. The noteworthy features of the sterile neutrino scenario advocated here is that the lightest one is almost degenerate with one of the active neutrinos, the second sterile has mass of order eV and the heaviest one is in the keV range. In the present scenario, the short baseline anomaly is explained through $\Delta m^2\sim 1\, {\rm eV^2}$, the depletion of muon neutrino flux in IceCube is explained through $\Delta m^2\sim 4.0\times 10^{-16}\, {\rm eV^2}$ and the dark matter problem is addressed through $\Delta m^2\sim 1\, {\rm keV^2}$. Our proposed mixing matrix is also compatible with the observed neutrino oscillation data. We show that the high energy muon and the tau neutrino fluxes from Gamma Ray Bursts can be depleted in IceCube by as much as 38% and 43% respectively. These substantial depletion in both muon and tau neutrino fluxes is due to their small but sizable mixing with the sterile neutrinos.


title: Modulation of galactic cosmic rays during the unusual solar minimum between cycles 23 and 24

abstract: During the recent solar minimum between cycles 23 and 24 (solar minimum $P_{23/24}$) the intensity of Galactic Cosmic Rays (GCRs) measured at the Earth was the highest ever recorded since space age. It is the purpose of this paper to resolve the most plausible mechanism for this unusually high intensity. A GCR transport model in three-dimensional heliosphere based on a simulation of Markov stochastic process is used to find the relation of cosmic ray modulation to various transport parameters, including solar wind (SW) speed, distance of heliospheric boundary, magnitude of interplanetary magnetic field (IMF) at the Earth, tilt angle of heliospheric current sheet (HCS), values of parallel and perpendicular diffusion coefficients. We calculate GCR proton energy spectra at the Earth for the last three solar minima $P_{21/22}$, $P_{22/23}$, and $P_{23/24}$, with the transport parameters obtained from observations. Besides weak IMF magnitude and slow SW speed, we find that a possible low magnetic turbulence, which increases the parallel diffusion and reduces the perpendicular diffusion in the polar direction, might be an additional possible mechanism for the high GCR intensity in the solar minimum $P_{23/24}$.


title: Feasibility Study of Real-Time Planning for Stereotactic Radiosurgery

abstract: Purpose: 3D rotational setup errors in radiotherapy are often ignored by most clinics due to inability to correct or simulate them accurately and efficiently. There are two types of rotation-related problems in a clinical setting. One is to assess the affected dose distribution in real-time if correction is not applied and the other one is to correct the rotational setup errors prior to the initiation of the treatment. Here, we present our analytical solutions to both problems.


title: Triply-resonant micro-optical parametric oscillators based on Kerr nonlinearity: nonlinear loss, unequal resonance-port couplings, and coupled-cavity implementations

abstract: We develop a theoretical model of triply-resonant optical parametric oscillators (OPOs) based on degenerate four-wave mixing (FWM) that includes physics and degrees of freedom relevant to microphotonic (on-chip) device implementations, including nonlinear loss, a general resonant mode field structure, and mode-selective coupling to external ports. The coupled mode theory model addresses the effect of two-photon absorption and free-carrier absorption on parametric gain and oscillation thresholds, and ultimately on the optimum design for an OPO. The model goes beyond a typical free-space cavity configuration by incorporating a full modal analysis that admits distributed modes with non-uniform field distribution, relevant to photonic microcavity systems on chip. This leads to a generalization of the concept of nonlinear figure of merit (NFOM) to a vector of coefficients. In addition, by considering unconstrained signal, pump and idler resonance coupling strengths to excitation ports, not usually available in simple cavity geometry, we show that the efficiency-maximizing design will have unequal external Q for the three resonances. We arrive at generalized formulas for OPO oscillation threshold that include nonlinear absorption and free carrier lifetime, and provide a normalized solution to the design problem in the presence of nonlinear loss in terms of optimum choice of coupling. Based on the results, we suggest a family of coupled-cavity systems to implement optimum FWM, where control of resonant wavelengths can be separated from optimizing nonlinear conversion efficiency, and where furthermore pump, signal, and idler coupling to bus waveguides can be controlled independently, using interferometric cavity supermode coupling as an example. Using the generalized NFOM, we address the efficiency of single and multi-cavity geometry, as well as standing and traveling wave excitation.


title: Three independent one-dimensional margins for single-fraction frameless stereotactic radiosurgery brain cases using CBCT

abstract: Setting a proper margin is crucial for not only delivering the required radiation dose to a target volume, but also reducing the unnecessary radiation to the adjacent organs at risk. This study investigated the independent one-dimensional symmetric and asymmetric margins between the clinical target volume (CTV) and the planning target volume (PTV) for linac-based single-fraction frameless stereotactic radiosurgery (SRS).


title: Toda lattice hierarchy and Goldstein-Petrich flows for plane curves

abstract: A relation between the Goldstein-Petrich hierarchy for plane curves and the Toda lattice hierarchy is investigated. A representation formula for plane curves is given in terms of a special class of $\tau$-functions of the Toda lattice hierarchy. A representation formula for discretized plane curves is also discussed.


title: Compound kernel estimates for the transition probability density of a L\'evy process in $\rn$

abstract: We construct in the small-time setting the upper and lower estimates for the transition probability density of a L\'evy process in $\rn$. Our approach relies on the complex analysis technique and the asymptotic analysis of the inverse Fourier transform of the characteristic function of the respective process.


title: Concentration of small Willmore spheres in Riemannian 3-manifolds

abstract: Given a 3-dimensional Riemannian manifold $(M,g)$, we prove that if $(\Phi_k)$ is a sequence of Willmore spheres (or more generally area-constrained Willmore spheres), having Willmore energy bounded above uniformly strictly by $8 \pi$, and Hausdorff converging to a point $\bar{p}\in M$, then $Scal(\bar{p})=0$ and $\nabla Scal(\bar{p})=0$ (resp. $\nabla Scal(\bar{p})=0$). Moreover, a suitably rescaled sequence smoothly converges, up to subsequences and reparametrizations, to a round sphere in the euclidean 3-dimensional space. This generalizes previous results of Lamm and Metzger contained in \cite{LM1}-\cite{LM2}. An application to the Hawking mass is also established.


title: Cellular network entropy as the energy potential in Waddington's differentiation landscape

abstract: Differentiation is a key cellular process in normal tissue development that is significantly altered in cancer. Although molecular signatures characterising pluripotency and multipotency exist, there is, as yet, no single quantitative mark of a cellular sample's position in the global differentiation hierarchy. Here we adopt a systems view and consider the sample's network entropy, a measure of signaling pathway promiscuity, computable from a sample's genome-wide expression profile. We demonstrate that network entropy provides a quantitative, in-silico, readout of the average undifferentiated state of the profiled cells, recapitulating the known hierarchy of pluripotent, multipotent and differentiated cell types. Network entropy further exhibits dynamic changes in time course differentiation data, and in line with a sample's differentiation stage. In disease, network entropy predicts a higher level of cellular plasticity in cancer stem cell populations compared to ordinary cancer cells. Importantly, network entropy also allows identification of key differentiation pathways. Our results are consistent with the view that pluripotency is a statistical property defined at the cellular population level, correlating with intra-sample heterogeneity, and driven by the degree of signaling promiscuity in cells. In summary, network entropy provides a quantitative measure of a cell's undifferentiated state, defining its elevation in Waddington's landscape.


title: The reduction of Laplace equation in certain Riemannian spaces and the resulting Type II hidden symmetries

abstract: We prove a general theorem which allows the determination of Lie symmetries of Laplace equation in a general Riemannian space using the conformal group of the space. Algebraic computing is not necessary. We apply the theorem in the study of the reduction of Laplace equation in certain classes of Riemannian spaces which admit a gradient Killing vector, a gradient Homothetic vector and a special Conformal Killing vector. In each reduction we identify the source of Type II hidden symmetries. We find that in general the Type II hidden symmetries of Laplace equation are directly related to the transition of the CKVs from the space where the original equation is defined to the space where the reduced equation resides. In particular we consider the reduction of Laplace equation (i.e. the wave equation) in Minkowski space and obtain the results of all previous studies in a straightforward manner. We consider the reduction of Laplace equation in spaces which admit Lie point symmetries generated from a non gradient HV and a proper CKV and we show that the reduction with these vectors does not produce Type II hidden symmetries. We apply the results to General Relativity and consider the reduction of Laplace equation in locally rotational symmetric space times (LRS) and in algebraically special vacuum solutions of Einstein's equations which admit a homothetic algebra acting simply transitively. In each case we determine the Type II hidden symmetries.


title: Generalized wreath products of graphs and groups

abstract: Inspired by the definition of generalized wreath product of permutation groups, we define the generalized wreath product of graphs, containing the classical Cartesian and wreath product of graphs as particular cases. We prove that the generalized wreath product of Cayley graphs of finite groups is the Cayley graph of the generalized wreath product of the corresponding groups.


title: Spatiotemporal Binary Interaction and Designer quasi particle condensates

abstract: We introduce a new integrable model to investigate the dynamics of two component quasi particle condensates with spatio temporal interaction strengths. We derive the associated Lax-pair of the coupled GP equation and construct matter wave solitons. We show that the spatio temporal binary interaction strengths not only facilitate the stabilization of the condensates, but also enables one to fabricate condensates with desirable densities, geometries and properties leading to the so called "designer quasi particle condensates".


title: QCDGPU: open-source package for Monte Carlo lattice simulations on OpenCL-compatible multi-GPU systems

abstract: The multi-GPU open-source package QCDGPU for lattice Monte Carlo simulations of pure SU(N) gluodynamics in external magnetic field at finite temperature and O(N) model is developed. The code is implemented in OpenCL, tested on AMD and NVIDIA GPUs, AMD and Intel CPUs and may run on other OpenCL-compatible devices. The package contains minimal external library dependencies and is OS platform-independent. It is optimized for heterogeneous computing due to the possibility of dividing the lattice into non-equivalent parts to hide the difference in performances of the devices used. QCDGPU has client-server part for distributed simulations. The package is designed to produce lattice gauge configurations as well as to analyze previously generated ones. QCDGPU may be executed in fault-tolerant mode. Monte Carlo procedure core is based on PRNGCL library for pseudo-random numbers generation on OpenCL-compatible devices, which contains several most popular pseudo-random number generators.


title: Elliptic Curves over Real Quadratic Fields are Modular

abstract: We prove that all elliptic curves defined over real quadratic fields are modular.


title: Resonant waves in elastic structured media: dynamic homogenisation versus Green's functions

abstract: We address an important issue of a dynamic homogenisation in vector elasticity for a doubly periodic mass-spring elastic lattice. The notion of logarithmically growing resonant waves is used in a complete analysis of star-shaped wave forms induced by an oscillating point force. We note that the dispersion surfaces for Floquet-Bloch waves in an elastic lattice main contain critical points of the saddle type. Based on the local quadratic approximations of the frequency, as a function of wave vector components, we deduce properties of a transient asymptotic solution as the contribution of the point source to the wave form. In this way, we describe local Green's functions as localized wave forms corresponding to the resonant frequency. The peculiarity of the problem lies in the fact that, at the same resonant frequency, the Taylor quadratic approximations for different groups of the resonant points are different, and hence we deal with different local Green's functions. Thus, there is no uniformly defined homogenisation procedure for a given resonant frequency. Instead, the continuous approximation of the wave field can be obtained through the asymptotic analysis of the lattice Green's functions.


title: Addition Chains Meet Postage Stamps: Reducing the Number of Multiplications

abstract: We introduce stamp chains. A stamp chain is a finite integer sequence that is both an addition chain and an additive 2-basis, i.e., a solution to the postage stamp problem. We provide a simple method for converting known postage stamp solutions of length k into stamp chains of length k+1. Using stamp chains, we construct an algorithm that computes u(x^i) for i=1,...,n in less than n-1 multiplications, if u is a function that can be computed at zero cost, and if there exists another zero-cost function v such that v(a,b) = u(ab). This can substantially reduce the computational cost of repeated multiplication, as illustrated by application examples related to matrix multiplication and data clustering using subset convolution. In addition, we report the extremal postage stamp solutions of length k=24.


title: Quasi-elastic scattering, RPA, 2p2h and neutrino--energy reconstruction

abstract: We discuss some nuclear effects, RPA correlations and 2p2h (multinucleon) mechanisms, on charged-current neutrino-nucleus reactions that do not produce a pion in the final state. We study a wide range of neutrino energies, from few hundreds of MeV up to 10 GeV. We also examine the influence of 2p2h mechanisms on the neutrino energy reconstruction.


title: Multiple spiral patterns in the transitional disk of HD 100546

abstract: Protoplanetary disks around young stars harbor many structures related to planetary formation. Of particular interest, spiral patterns were discovered among several of these disks and are expected to be the sign of gravitational instabilities leading to giant planets formation or gravitational perturbations caused by already existing planets. In this context, the star HD100546 presents some specific characteristics with a complex gas and dusty disk including spirals as well as a possible planet in formation. The objective of this study is to analyze high contrast and high angular resolution images of this emblematic system to shed light on critical steps of the planet formation. We retrieved archival images obtained at Gemini in the near IR (Ks band) with the instrument NICI and processed the data using advanced high contrast imaging technique taking advantage of the angular differential imaging. These new images reveal the spiral pattern previously identified with HST with an unprecedented resolution, while the large-scale structure of the disk is mostly erased by the data processing. The single pattern at the southeast in HST images is now resolved into a multi-armed spiral pattern. Using two models of a gravitational perturber orbiting in a gaseous disk we attempted to bring constraints on the characteristics of this perturber assuming each spiral being independent and we derived qualitative conclusions. The non-detection of the northeast spiral pattern observed in HST allows to put a lower limit on the intensity ratio between the two sides of the disk, which if interpreted as forward scattering yields a larger anisotropic scattering than derived in the visible. Also, we found that the spirals are likely spatially resolved with a thickness of about 5-10AU. Finally, we did not detect the candidate forming planet recently discovered in the Lp band, with a mass upper limit of 16-18 MJ.


title: Nominal Regular Expressions for Languages over Infinite Alphabets. Extended Abstract

abstract: We propose regular expressions to abstractly model and study properties of resource-aware computations. Inspired by nominal techniques -- as those popular in process calculi -- we extend classical regular expressions with names (to model computational resources) and suitable operators (for allocation, deallocation, scoping of, and freshness conditions on resources). We discuss classes of such nominal regular expressions, show how such expressions have natural interpretations in terms of languages over infinite alphabets, and give Kleene theorems to characterise their formal languages in terms of nominal automata.


title: Mechanism of preferential adsorption of SO$_2$ into two microporous paddle wheel frameworks M(bdc)(ted)0.5

abstract: The selective adsorption of a corrosive gas, SO$_2$, into two microporous pillared paddle-wheel frameworks M(bdc)(ted)0.5 [M = Ni, Zn; bdc = 1,4-benzenedicarboxylate; ted = triethylenediamine] is studied by volumetric adsorption measurements and a combination of in-situ infrared spectroscopy and ab initio density functional theory (DFT) calculations. The uptake of SO$_2$ in M(bdc)(ted)0.5 at room temperature is quite significant, 9.97 mol/kg at 1.13 bar. The major adsorbed SO$_2$ molecules contributing to the isotherm measurements are characterized by stretching bands at 1326 and 1144 cm$^{-1}$. Theoretical calculations including van der Waals interactions (based on vdW-DF) suggest that two adsorption configurations are possible for these SO$_2$ molecules. One geometry involves an SO$_2$ molecule bonded through its sulfur atom to the oxygen atom of the paddle-wheel building unit and its two oxygen atoms to the C-H groups of the organic linkers by formation of hydrogen bonds. Such a configuration results in a distortion of the benzene rings, which is consistent with the experimentally observed shift of the ring deformation mode. In the other geometry, SO$_2$ establishes hydrogen bonding with -CH$_2$ group of the ted linker through its two oxygen atoms simultaneously. The vdW-DF-simulated frequency shifts of the SO$_2$ stretching bands in these two configurations are similar and in good agreement with spectroscopically measured values of physisorbed SO$_2$.In addition, the IR spectra reveal the presence of another minor species, characterized by stretching modes at 1242 and 1105 cm$^{-1}$ and causing significant perturbations of MOFs vibrational modes (CH$_x$ and carboxylate groups). This species is more strongly bound, requiring a higher temperature ($\sim$150 $^\circ$C) to remove it than for the main physisorbed species.


title: Parameter estimation of monomial-exponential sums

abstract: We propose a numerical method, based upon matrix-pencils, for the identification of parameters and coefficients of a monomial-exponential sum. We note that this method can be considered an extension of the numerical methods for the parameter estimation of exponential sums. The application of the method is applied to several examples, some already present in the literature and others, to our knowledge, never considered before.


title: A topological necessary condition for the existence of compact Clifford-Klein forms

abstract: We provide a necessary condition for the existence of a compact Clifford-Klein form of a given homogeneous space of reductive type. The key to the proof is to combine a result of Kobayashi-Ono with an elementary fact that certain two different Clifford-Klein forms have the same cohomology ring. We give some examples, SL(p+q, R)/SO(p, q) (p, q: odd) for instance, of homogeneous spaces which do not admit compact Clifford-Klein forms.


title: Conformal Transformations and Weak Field Limit of Scalar-Tensor Gravity

abstract: The weak field limit of scalar tensor theories of gravity is discussed in view of conformal transformations. Specifically, we consider how physical quantities, like gravitational potentials derived in the Newtonian approximation for the same scalar-tensor theory, behave in the Jordan and in the Einstein frame. The approach allows to discriminate features that are invariant under conformal transformations and gives contributions in the debate of selecting the true physical frame. As a particular example, the case of $f(R)$ gravity is considered.


title: Charged Higgs Pair Production in a General Two Higgs Doublet Model at $e^+e^-$ and $\mu^+\mu^-$ Linear Colliders

abstract: In this paper, charged Higgs pair production through $\ell^+ \ell^- \rightarrow H^+ H^-$ where $\ell = e$ or $\mu$, is studied within the framework of a general Two Higgs Doublet Model (2HDM). The analysis is relevant to a future $e^+e^-$ or $\mu^+\mu^-$ collider operating at center of mass energy of $\sqrt{s}=500$ GeV. Two different scenarios of small and large $\alpha$ values is studied. Here $\alpha$ is the parameter which diagonalizes the neutral CP-even Higgs boson mass matrix. Within the Minimal Supersymmetric Standard Model (MSSM), cross section of this process is almost the same at $e^+e^-$ and $\mu^+\mu^-$ colliders. It is shown that at $e^+e^-$ colliders within a general 2HDM, cross section is not sensitive to the mass of neutral Higgs bosons, however, it can acquire large values up to several picobarn at $\mu^+\mu^-$ colliders with the presence of heavy neutral Higgs bosons. A scan over Higgs boson mass parameter space is performed to analyze the effect of large masses of neutral Higgs bosons involved in the s-channel propagator and thus in the total cross section of this process.


title: NACRE II: an update of the NACRE compilation of charged-particle-induced thermonuclear reaction rates for nuclei with mass number $A < 16$

abstract: An update of the NACRE compilation [Angulo et al., Nucl. Phys. A 656 (1999) 3] is presented. This new compilation, referred to as NACRE II, reports thermonuclear reaction rates for 34 charged-particle induced, two-body exoergic reactions on nuclides with mass number $A<16$, of which fifteen are particle-transfer reactions and the rest radiative capture reactions. When compared with NACRE, NACRE II features in particular (1) the addition to the experimental data collected in NACRE of those reported later, preferentially in the major journals of the field by early 2013, and (2) the adoption of potential models as the primary tool for extrapolation to very low energies of astrophysical $S$-factors, with a systematic evaluation of uncertainties. As in NACRE, the rates are presented in tabular form for temperatures in the $10^{6}$ $\simeq\leq$ T $\leq$ $10^{10}$ K range. Along with the 'adopted' rates, their low and high limits are provided. The new rates are available in electronic form as part of the Brussels Library (BRUSLIB) of nuclear data. The NACRE II rates also supersede the previous NACRE rates in the Nuclear Network Generator (NETGEN) for astrophysics. [http://www.astro.ulb.ac.be/databases.html.]


title: Semi-Decentralized Approximation of Optimal Control of Distributed Systems Based on a Functional Calculus

abstract: This paper discusses a new approximation method for operators which are solution to an operational Riccati equation (ORE). The latter is derived from the theory of optimal control of linear problems posed in Hilbert spaces. The approximation is based on the functional calculus of self-adjoint operators and the Cauchy formula. Under a number of assumptions the approximation is suitable for implementation on a semi-decentralized computing architecture in view of real-time control. Our method is particularly applicable to problems in optimal control of systems governed by partial differential equations with distributed observation and control. Some relatively academic applications are presented for illustration. More realistic examples relating to microsystem arrays have already been published.


title: Possible solution to the riddle of HD 82943 multiplanet system: the three-planet resonance 1:2:5?

abstract: We carry out a new analysis of the published radial velocity data for the planet-hosting star HD82943. We include the recent Keck/HIRES measurements as well as the aged but much more numerous CORALIE data. We find that the CORALIE radial velocity measurements are polluted by a systematic annual variation which affected the robustness of many previous results. We show that after purging this variation, the residuals still contain a clear signature of an additional $\sim 1100$ days periodicity. The latter variation leaves significant hints in all three independent radial velocity subsets that we analysed: the CORALIE data, the Keck data acquired prior to a hardware upgrade, and the Keck data taken after the upgrade. We mainly treat this variation as a signature of a third planet in the system, although we cannot rule out other interpretations, such as long-term stellar activity. We find it easy to naturally obtain a stable three-planet radial-velocity fit close to the three-planet mean-motion resonance 1:2:5, with the two main planets (those in the 1:2 resonance) in an aligned apsidal corotation. The dynamical status of the third planet is still uncertain: it may reside in as well as slightly out of the 5:2 resonance. We obtain the value of $\sim 1075$ days for its orbital period and of $\sim 0.3 M_{\rm Jup}$ for its minimum mass, while the eccentric parameters are uncertain.


title: Blow-up of Smooth Solutions to the Euler-Poisson Equations

abstract: In this paper, the finite time blow-up of smooth solutions to the Cauchy problem for full Euler-Poisson equations and isentropic Euler-Poisson equations with repulsive forces or attractive forces in high dimensions $(n\geq3)$ is proved for a large class of initial data. It is not required that the initial data has compact support or contains vacuum in any finite regions.


title: Higher-order Changhee numbers and polynomials

abstract: In this paper, we consider the higher-order Changhee numbers and polynomials which are derived from the fermionic p-adic integral on Zp and give some relations between higher-order Changhee polynomials and special polynomials.


title: The use of gravitational lenses in the study of distant galaxy mergers

abstract: Gravlenses are efficiently explored for detecting the most distant galaxies (up to z=10 redshifts). As an example of the role played by gravlenses we refer to the observation of the galaxy merger at z=3 (Borys, et al; Berciano Alba, et al). We derived solutions for the Smoluchowski kinetic equation for the mass function of galaxies, which describes mergers in differential approximation (minor mergers). It is shown that the evolution of the slope of luminosity function observed in the Ultra Deep Hubble Field (Bouwence et al) can be described as a result of explosive evolution driven by galaxy mergers.


title: Scattering from dust in molecular clouds: Constraining the dust grain size distribution through near-infrared cloudshine and infrared coreshine

abstract: Context. The largest grains (0.5-1 micron) in the interstellar size distribution are efficient in scattering near- and mid-infrared radiation. These wavelengths are therefore particularly well suited to probe the still uncertain high-end of the size distribution. Aims. We investigate the change in appearance of a typical low-mass molecular core from the Ks (2.2 micron) band to the Spitzer IRAC 3.6 and 8 micron bands, and compare with model calculations, which include variations of the grain size distribution. Methods. We combine Spitzer IRAC and ground-based near-infrared observations to characterize the scattered light observed at the near- and mid-infrared wavelengths from the core L260. Using a spherical symmetric model core, we perform radiative transfer calculations to study the impact of various dust size distributions on the intensity profiles across the core. Results. The observed scattered light patterns in the Ks and 3.6 micron bands are found to be similar. By comparison with radiative transfer models the two profiles place constraints on the relative abundance of small and large (more than 0.25 micron) dust grains. The scattered light profiles are found to be inconsistent with an interstellar silicate grain distribution extending only to 0.25 micron and large grains are needed to reach the observed fluxes and the flux ratios. The shape of the Ks band surface brightness profile limits the largest grains to 1-1.5 micron. Conclusions. In addition to observing coreshine in the Spitzer IRAC channels, the combination with ground-based near-infrared observations are suited to constrain the properties of large grains in cores.


title: Evanescent Artificial Gauge Potentials for Neutral Atoms

abstract: We show that atoms interacting with evanescent light fields, generated at the interface of a dielectric with vacuum, experience artificial gauge potentials. These potentials depend crucially on the physical parameters which characterize the evanescent fields most notably the refractive index of the dielectric material and the angle of incidence of the laser beam totally internaly reflected at the interface. Gauge fields are derived for various evanescent light fields and for both two-level and three-level systems. The use of such artificial gauge potentials for the manipulation of atoms trapped at the interfaces is pointed out and discussed.


title: Multi-shell effective interactions

abstract: Background: Effective interactions, either derived from microscopic theories or based on fitting selected properties of nuclei in specific mass regions, are widely used inputs to shell-model studies of nuclei. Until recently, most shell-model calculations have been confined to a single oscillator shell. Recent interest in nuclei away from the stability line, requires however larger shell-model spaces. Since the derivation of microscopic effective interactions has been limited to degenerate model spaces, there are both conceptual and practical limits to present shell-model calculations that utilize such interactions. Purpose: The aim of this work is to present a novel microscopic method to calculate effective interactions for the nuclear shell model. Its main difference from existing theories is that it can be applied not only to degenerate model spaces but also to non-degenerate model spaces. Methods: The formalism is presented in the form of many-body perturbation theory based on the recently developed Extended Kuo-Krenciglowa method. Our method enables us to microscopically construct effective interactions not only in one oscillator shell but also for several oscillator shells. Results: We present numerical results using effective interactions within (i) a single oscillator shell like the sd-shell and the pf-shell, and (ii) two major shells like the sdf7p3-shell and the pfg9-shell. We also present energy levels of several nuclei which have two valence nucleons on top of a given closed-shell core. Conclusions: Our results show that the present method works excellently in shell-model spaces that comprise several oscillator shells, as well as in a single oscillator shell. We show in particular that the microscopic inter-shell interactions are much more attractive than has been expected by degenerate perturbation theory. The consequences for shell-model studies are discussed.


title: Finiteness of hitting times under taboo

abstract: We consider a continuous-time Markov chain with a finite or countable state space. For a site y and subset H of the state space, the hitting time of y under taboo H is defined to be infinite if the process trajectory hits H before y, and the first hitting time of y otherwise. We investigate the probability that such times are finite. In particular, if the taboo set is finite, an efficient iterative scheme reduces the study to the known case of a singleton taboo. A similar procedure applies in the case of finite complement of the taboo set. The study is motivated by classification of branching processes with finitely many catalysts. Keywords and phrases: Markov chain, hitting time, taboo probabilities, catalytic branching process.


title: Crossover from the weak to strong-field behavior of the longitudinal interlayer magnetoresistance in quasi-two-dimensional conductors

abstract: We investigate the monotonic growth of longitudinal interlayer magnetoresistance $\bar{R}_{zz}(B_z) $, analytically and numerically in the self-consistent Born approximation. We show that in a weak magnetic field the monotonic part of $\bar{R}_{zz}(B_z)$ is almost constant and starts to grow only above the crossover field $B_{c}$, when the Landau levels (LL) become isolated, i.e. when the LL separation becomes greater than the LL broadening. In higher field $B_{z}>>B_{c}$, $\bar{R}_{zz}(B_{z}) \propto B_{z}^{1/2}$ in agreement with previous works.


title: Nonlinear Bethe-Heitler Pair Creation in an Intense Two-Mode Laser Field

abstract: We investigate electron-positron pair creation in the interaction of a nuclear Coulomb field and a highly intense two-mode laser field. For bichromatic laser fields, we examine the differences arising for commensurable and incommensurable frequencies in a continuous variation of the laser frequency ratio and the quantum interference effects, which may occur in the commensurable case. We show that the interference manifests in the angular distributions and the total pair-production rates of the created particles. Additionally, by varying the amplitudes of the two modes we study pair creation in a monochromatic laser wave of arbitrarily elliptical polarization.


title: Conservation laws, exact travelling waves and modulation instability for an extended nonlinear Schr\"odinger equation

abstract: We study various properties of solutions of an extended nonlinear Schr\"{o}dinger (ENLS) equation, which arises in the context of geometric evolution problems -- including vortex filament dynamics -- and governs propagation of short pulses in optical fibers and nonlinear metamaterials. For the periodic initial-boundary value problem, we derive conservation laws satisfied by local in time, weak $H^2$ (distributional) solutions, and establish global existence of such weak solutions. The derivation is obtained by a regularization scheme under a balance condition on the coefficients of the linear and nonlinear terms -- namely, the Hirota limit of the considered ENLS model. Next, we investigate conditions for the existence of traveling wave solutions, focusing on the case of bright and dark solitons. The balance condition on the coefficients is found to be essential for the existence of exact analytical soliton solutions; furthermore, we obtain conditions which define parameter regimes for the existence of traveling solitons for various linear dispersion strengths. Finally, we study the modulational instability of plane waves of the ENLS equation, and identify important differences between the ENLS case and the corresponding NLS counterpart. The analytical results are corroborated by numerical simulations, which reveal notable differences between the bright and the dark soliton propagation dynamics, and are in excellent agreement with the analytical predictions of the modulation instability analysis.


title: Computation Over Gaussian Networks With Orthogonal Components

abstract: Function computation of arbitrarily correlated discrete sources over Gaussian networks with orthogonal components is studied. Two classes of functions are considered: the arithmetic sum function and the type function. The arithmetic sum function in this paper is defined as a set of multiple weighted arithmetic sums, which includes averaging of the sources and estimating each of the sources as special cases. The type or frequency histogram function counts the number of occurrences of each argument, which yields many important statistics such as mean, variance, maximum, minimum, median, and so on. The proposed computation coding first abstracts Gaussian networks into the corresponding modulo sum multiple-access channels via nested lattice codes and linear network coding and then computes the desired function by using linear Slepian-Wolf source coding. For orthogonal Gaussian networks (with no broadcast and multiple-access components), the computation capacity is characterized for a class of networks. For Gaussian networks with multiple-access components (but no broadcast), an approximate computation capacity is characterized for a class of networks.


title: Singleton sets random attractor for stochastic FitzHugh-Nagumo lattice equations driven by fractional Brownian motions

abstract: The paper is devoted to the study of the dynamical behavior of the solutions of stochastic FitzHugh-Nagumo lattice equations, driven by fractional Brownian motions, with Hurst parameter greater than $1/2$. Under some usual dissipativity conditions, the system considered here features different dynamics from the same one perturbed by Brownian motion. In our case, the random dynamical system has a unique random equilibrium, which constitutes a singleton sets random attractor.


title: Efficient Information Theoretic Clustering on Discrete Lattices

abstract: We consider the problem of clustering data that reside on discrete, low dimensional lattices. Canonical examples for this setting are found in image segmentation and key point extraction. Our solution is based on a recent approach to information theoretic clustering where clusters result from an iterative procedure that minimizes a divergence measure. We replace costly processing steps in the original algorithm by means of convolutions. These allow for highly efficient implementations and thus significantly reduce runtime. This paper therefore bridges a gap between machine learning and signal processing.