title: Generating Merger Trees for Dark Matter Haloes: A Comparison of Methods

abstract: Halo merger trees describe the hierarchical mass assembly of dark matter haloes, and are the backbone for modeling galaxy formation and evolution. Merger trees constructed using Monte Carlo algorithms based on the extended Press-Schechter (EPS) formalism are complementary to those extracted from N-body simulations, and have the advantage that they are not trammeled by limited numerical resolution and uncertainties in identifying (sub)haloes and linking them between snapshots. This paper compares multiple EPS-based merger tree algorithms to simulation results using four diagnostics: progenitor mass function (PMF), mass assembly history (MAH), merger rate per descendant halo, and the unevolved subhalo mass function (USMF). In general, algorithms based on spherical collapse yield major-merger rates that are too high by a factor of two, resulting in MAHs that are systematically offset. Assuming ellipsoidal collapse solves most of these issues, but the particular algorithm investigated here that incorporates ellipsoidal collapse dramatically overpredicts the minor-merger rate for massive haloes. The only algorithm in our comparison that yields MAHs, merger rates, and USMFs in good agreement with simulations, is that by Parkinson et al. (2008). However, this is not a true EPS-based algorithm as it draws its progenitor masses from a PMF calibrated against simulations, rather than `predicted' by EPS. Finally we emphasize that the benchmarks used to test the EPS algorithms are obtained from simulations and are hampered by significant uncertainties themselves. In particular, MAHs and halo merger rates obtained from simulations by different authors reveal discrepancies that easily exceed 50 percent, even when based on the same simulation. Given this status quo, merger trees constructed using the Parkinson et al. algorithm are as accurate as those extracted from N-body simulations.


title: Cosmological Constraints from the Anisotropic Clustering Analysis using BOSS DR9

abstract: Our observations of the Universe are fundamentally anisotropic, with data from galaxies separated transverse to the line of sight coming from the same epoch while that from galaxies separated parallel to the line of sight coming from different times. Moreover, galaxy velocities along the line of sight change their redshift, giving redshift space distortions. We perform a full two-dimensional anisotropy analysis of galaxy clustering data, fitting in a substantially model independent manner the angular diameter distance D_A, Hubble parameter H, and growth rate ddelta/dln a without assuming a dark energy model. The results demonstrate consistency with LCDM expansion and growth, hence also testing general relativity. We also point out the interpretation dependence of the effective redshift z_eff, and its cosmological impact for next generation surveys.


title: ADHM Construction of Noncommutative Instantons

abstract: We discuss the Atiyah-Drinfeld-Hitchin-Manin (ADHM) construction of U(N) instantons in noncommutative (NC) space and prove the one-to-one correspondence between moduli spaces of the noncommutative instantons and the ADHM data, together with an origin of the instanton number for U(1). We also give a derivation of the ADHM construction from the viewpoint of the Nahm transformation of instantons on four-tori. This article is a composite version of [23] and [24].


title: Hidden Starbursts and Active Galactic Nuclei at 0<z<4 from the Herschel-VVDS-CFHTLS-D1 Field: Inferences on Coevolution and Feedback

abstract: We investigate of the properties of $\sim$2000 Herschel/SPIRE-selected galaxies from $0<z<4$ using a combination of extensive spectroscopy, deep imaging from CFHT, VLA, Spitzer, XMM-Newton, and Herschel, and well-calibrated SED fitting. Herschel galaxies are observed to span a range of stellar masses, colors, and absolute magnitudes equivalent to galaxies undetected in SPIRE. Though many Herschel galaxies appear to be in transition, such galaxies are largely consistent with normal star-forming galaxies when rest-frame colors are utilized. The nature of the star-forming "main sequence" is studied and we warn against adopting this framework unless the main sequence is determined precisely. Herschel galaxies at different total infrared luminosities ($L_{TIR}$) are compared. Bluer colors, larger nebular extinctions, and larger contributions from younger stellar populations are observed for galaxies with larger $L_{TIR}$, suggesting that low-$L_{TIR}$ galaxies are undergoing rejuvenated starbursts while galaxies with higher $L_{TIR}$ are forming a larger percentage of their stellar mass. A variety of methods are used to select powerful active galactic nuclei (AGN). Galaxies hosting AGN are observed to be undergoing starbursts more commonly and vigorously than a matched sample of galaxies without powerful AGN and, additionally, the fraction of galaxies with an AGN increases with increasing star formation rate at all redshifts. At all redshifts ($0<z<4$) the most prodigious star-forming galaxies are found to contain the highest fraction of powerful AGN. For redshift bins that allow a comparison ($z>0.5$), the highest $L_{TIR}$ galaxies in a given redshift bin are unobserved by SPIRE at subsequently lower redshifts, a trend linked to downsizing. In conjunction with other results, this evidence is used to argue for prevalent AGN-driven quenching in starburst galaxies across cosmic time.


title: RUN DMC: An efficient, parallel code for analyzing Radial Velocity Observations using N-body Integrations and Differential Evolution Markov chain Monte Carlo

abstract: In the 20+ years of Doppler observations of stars, scientists have uncovered a diverse population of extrasolar multi-planet systems. A common technique for characterizing the orbital elements of these planets is Markov chain Monte Carlo (MCMC), using a Keplerian model with random walk proposals and paired with the Metropolis-Hastings algorithm. For approximately a couple of dozen planetary systems with Doppler observations, there are strong planet-planet interactions due to the system being in or near a mean-motion resonance (MMR). An N-body model is often required to accurately describe these systems. Further computational difficulties arise from exploring a high-dimensional parameter space ($\sim$7 x number of planets) that can have complex parameter correlations. To surmount these challenges, we introduce a differential evolution MCMC (DEMCMC) applied to radial velocity data while incorporating self-consistent N-body integrations. Our Radial velocity Using N-body DEMCMC (RUN DMC) algorithm improves upon the random walk proposal distribution of the traditional MCMC by using an ensemble of Markov chains to adaptively improve the proposal distribution. We describe the methodology behind the algorithm, along with results of tests for accuracy and performance. We find that most algorithm parameters have a modest effect on the rate of convergence. However, the size of the ensemble can have a strong effect on performance. We show that the optimal choice depends on the number of planets in a system, as well as the computer architecture used and the resulting extent of parallelization. While the exact choices of optimal algorithm parameters will inevitably vary due to the details of individual planetary systems, we offer recommendations for choosing the DEMCMC algorithm's algorithmic parameters that result in excellent performance for a wide variety of planetary systems.


title: The stellar mass function and efficiency of galaxy formation with a varying initial mass function

abstract: Several recent observational studies have concluded that the initial mass function (IMF) of stars varies systematically with galaxy properties such as velocity dispersion. In this paper, we investigate the effect of linking the circular velocity of galaxies, as determined from the Fundamental Plane and Tully-Fisher relations, to the slope of the IMF with parameterizations guided by several of these studies. For each empirical relation, we generate stellar masses of ~600,000 SDSS galaxies at z ~ 0.1, by fitting the optical photometry to large suites of synthetic stellar populations that sample the full range of galaxy parameters. We generate stellar mass functions and examine the stellar-to-halo mass relations using sub-halo abundance matching. At the massive end, the stellar mass functions become a power law, instead of the familiar exponential decline. As a result, it is a generic feature of these models that the central galaxy stellar-to-halo mass relation is significantly flatter at high masses (slope ~ -0.3 to -0.4) than in the case of a universal IMF (slope ~ -0.6). We find that regardless of whether the IMF varies systematically in all galaxies or just early types, there is still a well-defined peak in the central stellar-to-halo mass ratio at halo masses of ~ 10E12 solar masses. In general, the IMF variations explored here lead to significantly higher integrated stellar densities if the assumed dependence on circular velocity applies to all galaxies, including late-types; in fact the more extreme cases can be ruled out, as they imply an unphysical situation in which the stellar fraction exceeds the universal baryon fraction.


title: Topological gaps without masses in driven graphene-like systems

abstract: We illustrate the possibility of realizing band gaps in graphene-like systems that fall outside the existing classification of gapped Dirac Hamiltonians in terms of masses. As our primary example we consider a band gap arising due to time-dependent distortions of the honeycomb lattice. By means of an exact, invertible, and transport-preserving mapping to a time-independent Hamiltonian, we show that the system exhibits Chern-insulating phases with quantized Hall conductivities $\pm e^2/h$. The chirality of the corresponding gapless edge modes is controllable by both the frequency of the driving and the manner in which sublattice symmetry is broken by the dynamical lattice modulations. Finally, we discuss a promising possible realization of this physics in photonic lattices.


title: Star Formation in Orion's L1630 Cloud: an Infrared and Multi-epoch X-ray Study

abstract: X-ray emission is characteristic of young stellar objects (YSOs) and is known to be highly variable. We investigate, via an infrared and multi-epoch X-ray study of the L1630 dark cloud, whether and how X-ray variability in young stellar objects is related to protostellar evolutionary state. We have analyzed 11 Chandra X-ray Observatory observations, obtained over the course of four years and totaling ~240 ks exposure time, targeting the eruptive Class I YSO V1647 Ori in L1630. We used 2MASS and Spitzer data to identify and classify IR counterparts to L1630 X-ray sources and identified a total of 52 X-ray emitting YSOs with IR counterparts, including 4 Class I sources and 1 Class 0/I source. We have detected cool (< 3 MK) plasma, possibly indicative of accretion shocks, in three classical T Tauri stars. A subsample of 27 X-ray-emitting YSOs were covered by 9 of the 11 Chandra observations targeting V1647 Ori and vicinity. For these 27 YSOs, we have constructed X-ray light curves spanning approximately four years. These light curves highlight the variable nature of pre-main sequence X-ray emitting young stars; many of the L1630 YSOs vary by orders of magnitude in count rate between observations. We discuss possible scenarios to explain apparent trends between various X-ray spectral properties, X-ray variance and YSO classification.


title: Topological Criticality in the Chiral-Symmetric AIII Class at Strong Disorder

abstract: The chiral AIII symmetry class in the periodic table of topological insulators contains topological phases classified by a winding number $\nu$ for each odd space-dimension. An open problem for this class is the characterization of the phases and phase-boundaries in the presence of strong disorder. In this work, we derive a covariant real-space formula for $\nu$ and, using an explicit 1-dimensional disordered topological model, we show that $\nu$ remains quantized and non-fluctuating when disorder is turned on, even though the bulk energy-spectrum is completely localized. Furthermore, $\nu$ remains robust even after the insulating gap is filled with localized states, but when the disorder is increased even further, an abrupt change of $\nu$ to a trivial value is observed. Using exact analytic calculations, we show that this marks a critical point where the localization length diverges. As such, in the presence of disorder, the AIII class displays a markedly different physics from everything known to date, with robust invariants being carried entirely by localized states and bulk extended states emerging from an absolutely localized spectrum. Detailed maps and a clear physical description of the phases and phase boundaries are presented based on numerical and exact analytic calculations.


title: Jet Veto Clustering Logarithms Beyond Leading Order

abstract: Many experimental analyses separate events into exclusive jet bins, using a jet algorithm to cluster the final state and then veto on jets. Jet clustering induces logarithmic dependence on the jet radius R in the cross section for exclusive jet bins, a dependence that is poorly controlled due to the non-global nature of the clustering. At jet radii of experimental interest, the leading order (LO) clustering effects are numerically significant, but the higher order effects are currently unknown. We rectify this situation by calculating the most important part of the next-to-leading order (NLO) clustering logarithms of R for any 0-jet process, which enter as $O(\alpha_s^3)$ corrections to the cross section. The calculation blends subtraction methods for NLO calculations with factorization properties of QCD and soft-collinear effective theory (SCET). We compare the size of the known LO and new NLO clustering logarithms and find that the impact of the NLO terms on the 0-jet cross section in Higgs production is small. This brings clustering effects under better control and may be used to improve uncertainty estimates on cross sections with a jet veto.


title: Zero-Field muSR Search for a Time-Reversal-Symmetry-Breaking Mixed Pairing State in Superconducting Ba1-xKxFe2As2

abstract: We report the results of a zero-field muon spin relaxation (ZF-muSR) study of superconducting Ba1-xKxFe2As2 (0.5 < x < 0.9) in search of weak spontaneous internal magnetic fields associated with proposed time-reversal-symmetry breaking mixed pairing states. The measurements were performed on polycrystalline samples, which do not exhibit the mesoscopic phase separation previously observed in single crystals of Ba1-xKxFe2As2. No evidence of spontaneous internal magnetic fields is found in any of the samples at temperatures down to T ~ 0.02 K.


title: Timing variations in the secondary eclipse of NN Ser

abstract: The eclipsing white dwarf plus main-sequence binary NN Serpentis provides one of the most convincing cases for the existence of circumbinary planets around evolved binaries. The exquisite timing precision provided by the deep eclipse of the white dwarf has revealed complex variations in the eclipse arrival times over the last few decades. These variations have been interpreted as the influence of two planets in orbit around the binary. Recent studies have proved that such a system is dynamically stable over the current lifetime of the binary. However, the existence of such planets is by no means proven and several alternative mechanisms have been proposed that could drive similar variations. One of these is apsidal precession, which causes the eclipse times of eccentric binaries to vary sinusoidally on many year timescales. In this paper we present timing data for the secondary eclipse of NN Ser and show that they follow the same trend seen in the primary eclipse times, ruling out apsidal precession as a possible cause for the variations. This result leaves no alternatives to the planetary interpretation for the observed period variations, although we still do not consider their existence as proven. Our data limits the eccentricity of NN Ser to e<0.001. We also detect a 3.3+/-1.0 second delay in the arrival times of the secondary eclipses relative to the best planetary model. This delay is consistent with the expected 2.84+/-0.04 second Romer delay of the binary, and is the first time this effect has been detected in a white dwarf plus M dwarf system.


title: On The Formation Of Multiple Absorption Troughs In Broad Absorption Line QSOs

abstract: We present theoretical CIV 1548,1550 absorption line profiles for QSOs calculated assuming the accretion disk wind (ADW) scenario. The results suggest that the multiple absorption troughs seen in many QSOs may be due to the discontinuities in the ion balance of the wind (caused by X-rays), rather than discontinuities in the density/velocity structure. The profiles are calculated from a 2.5D time-dependent hydrodynamic simulation of a line-driven disk wind for a typical QSO black hole mass, a typical QSO luminosity, and for a standard Shakura-Sunyaev disk. We include the effects of ionizing X-rays originating from within the inner disk radius by assuming that the wind is shielded from the X-rays from a certain viewing angle up to 90o ("edge on"). In the shielded region we assume constant ionization equilibrium, and thus constant line-force parameters. In the non-shielded region we assume that both the line-force and the C+3 populations are nonexistent. The model, at viewing angles close to the angle that separates the shielded region and the non-shielded region, produces absorption lines with multiple troughs. The steady nature of accretion disk winds, in turn, may account for the steady nature of the absorption structure observed in multiple-trough broad absorption line QSOs. The model parameters are Mbh = 10^9 Msun and Ldisk = 10^47 erg s-1.


title: Evidence for High-Energy Extraterrestrial Neutrinos at the IceCube Detector

abstract: We report on results of an all-sky search for high-energy neutrino events interacting within the IceCube neutrino detector conducted between May 2010 and May 2012. The search follows up on the previous detection of two PeV neutrino events, with improved sensitivity and extended energy coverage down to approximately 30 TeV. Twenty-six additional events were observed, substantially more than expected from atmospheric backgrounds. Combined, both searches reject a purely atmospheric origin for the twenty-eight events at the $4\sigma$ level. These twenty-eight events, which include the highest energy neutrinos ever observed, have flavors, directions, and energies inconsistent with those expected from the atmospheric muon and neutrino backgrounds. These properties are, however, consistent with generic predictions for an additional component of extraterrestrial origin.


title: General aspects of the vapor growth of semiconductor crystals - a study based on DFT simulations of the NH3/NH2 covered GaN(0001) surface in hydrogen ambient

abstract: Vapor growth of semiconductors is analyzed using recently obtained dependence of the adsorption energy on the electron charge transfer between the surface adsorbed species and the bulk [Krukowski et al. J. Appl. Phys. 114 (2013) 063507, Kempisty et al. ArXiv 1307.5778 (2013)]. Ab initio calculations were performed to study the physical properties of GaN(0001) surface in ammonia-rich conditions, i.e. covered by mixture of NH3 molecules and NH2 radicals. The Fermi level is pinned at valence band maximum (VBM) and conduction band minimum (CBM) for full coverage by NH3 molecules and NH2 radicals, respectively. For the crossover content of ammonia of about 25% monolayer (ML), the Fermi level is unpinned. It was shown that hydrogen adsorption energy depends on the doping in the bulk for the unpinned Fermi level, i.e. for this coverage. Surface structure thermodynamic and mechanical stability criteria are defined and compared. Mechanical stability of the coverage of such surfaces was checked by determination of the desorption energy of hydrogen molecules. Thermodynamic stability analysis indicates that initally equilibrium hydrogen vapor partial pressure steeply increases with NH3 content to attain the crossover NH3/NH2 coverage, i.e. the unpinned Fermi level condition. For such condition the entire range of experimentally accessible pressures belongs showing that vapor growth of semiconductor crystals occurs predominantly for unpinned Fermi level at the surface, i.e. for flat bands. Accordingly, adsorption energy of most species depends on the doping in the bulk that is basis of the possible molecular scenario explaining dependence of the growth and the doping of semiconductor crystals on the doping in the bulk


title: PENLAB: A MATLAB solver for nonlinear semidefinite optimization

abstract: PENLAB is an open source software package for nonlinear optimization, linear and nonlinear semidefinite optimization and any combination of these. It is written entirely in MATLAB. PENLAB is a young brother of our code PENNON \cite{pennon} and of a new implementation from NAG \cite{naglib}: it can solve the same classes of problems and uses the same algorithm. Unlike PENNON, PENLAB is open source and allows the user not only to solve problems but to modify various parts of the algorithm. As such, PENLAB is particularly suitable for teaching and research purposes and for testing new algorithmic ideas. In this article, after a brief presentation of the underlying algorithm, we focus on practical use of the solver, both for general problem classes and for specific practical problems.


title: Evolution of linear perturbations in spherically symmetric dust spacetimes

abstract: We present results from a numerical code implementing a new method to solve the master equations describing the evolution of linear perturbations in a spherically symmetric but inhomogeneous background. This method can be used to simulate several configurations of physical interest, such as relativistic corrections to structure formation, the lensing of gravitational waves and the evolution of perturbations in a cosmological void model. This paper focuses on the latter problem, i.e. structure formation in a Hubble scale void in the linear regime. This is considerably more complicated than linear perturbations of a homogeneous and isotropic background because the inhomogeneous background leads to coupling between density perturbations and rotational modes of the spacetime geometry, as well as gravitational waves. Previous analyses of this problem ignored this coupling in the hope that the approximation does not affect the overall dynamics of structure formation in such models. We show that for a giga-parsec void, the evolution of the density contrast is well approximated by the previously studied decoupled evolution only for very large-scale modes. However, the evolution of the gravitational potentials within the void is inaccurate at more than the $10\%$ level, and is even worse on small scales.


title: An Improved Variable Step-size Affine Projection Sign Algorithm for Echo Cancellation

abstract: This paper proposes an improved variable step-size (VSS) algorithm for the recently introduced affine projection sign algorithm (APSA) based on the recovery of the near-end signal energy in the error signal. Simulation results demonstrate that, compared to the previous VSS for APSA, the proposed approach provides both more robustness to impulse interference and better tracking ability of echo path change.


title: Solar Eruptive Events (SEE) 2020 Mission Concept

abstract: Major solar eruptive events (SEEs), consisting of both a large flare and a near simultaneous large fast coronal mass ejection (CME), are the most powerful explosions and also the most powerful and energetic particle accelerators in the solar system, producing solar energetic particles (SEPs) up to tens of GeV for ions and hundreds of MeV for electrons. The intense fluxes of escaping SEPs are a major hazard for humans in space and for spacecraft. Furthermore, the solar plasma ejected at high speed in the fast CME completely restructures the interplanetary medium (IPM) - major SEEs therefore produce the most extreme space weather in geospace, the interplanetary medium, and at other planets. Thus, understanding the flare/CME energy release process(es) and the related particle acceleration processes are major goals in Heliophysics. To make the next major breakthroughs, we propose a new mission concept, SEE 2020, a single spacecraft with a complement of advanced new instruments that focus directly on the coronal energy release and particle acceleration sites, and provide the detailed diagnostics of the magnetic fields, plasmas, mass motions, and energetic particles required to understand the fundamental physical processes involved.


title: A Generalized Markov-Chain Modelling Approach to $(1,\lambda )$-ES Linear Optimization

abstract: The manuscript generalizes several recent results of the 2nd author concerning Markov-Chain Modelling of $(1,\lambda )$-ES Linear Optimization.


title: NuSTAR Observations of GRB130427A establish a single component synchrotron afterglow origin for the late optical to multi-GeV emission

abstract: GRB 130427A occurred in a relatively nearby galaxy; its prompt emission had the largest GRB fluence ever recorded. The afterglow of GRB 130427A was bright enough for the Nuclear Spectroscopic Telescope ARray (NuSTAR) to observe it in the 3-79 keV energy range long after its prompt emission (~1.5 and 5 days). This range, where afterglow observations were previously not possible, bridges an important spectral gap. Combined with Swift, Fermi and ground-based optical data, NuSTAR observations unambiguously establish a single afterglow spectral component from optical to multi-GeV energies a day after the event, which is almost certainly synchrotron radiation. Such an origin of the late-time Fermi/LAT > 10 GeV photons requires revisions in our understanding of collisionless relativistic shock physics.


title: 2MASS Photometric Redshift catalog: a comprehensive three-dimensional census of the whole sky

abstract: Key cosmological applications require the three-dimensional galaxy distribution on the entire celestial sphere. These include measuring the gravitational pull on the Local Group, estimating the large-scale bulk flow and testing the Copernican principle. However, the largest all-sky redshift surveys -- the 2MRS and IRAS PSCz -- have median redshifts of only z=0.03 and sample the very local Universe. There exist all-sky galaxy catalogs reaching much deeper -- SuperCOSMOS in the optical, 2MASS in the near-IR and WISE in the mid-IR -- but these lack complete redshift information. At present, the only rapid way towards larger 3D catalogs covering the whole sky is through photometric redshift techniques. In this paper we present the 2MASS Photometric Redshift catalog (2MPZ) containing 1 million galaxies, constructed by cross-matching 2MASS XSC, WISE and SuperCOSMOS all-sky samples and employing the artificial neural network approach (the ANNz algorithm), trained on such redshift surveys as SDSS, 6dFGS and 2dFGRS. The derived photometric redshifts have errors nearly independent of distance, with an all-sky accuracy of sigma_z=0.015 and a very small percentage of outliers. In this way, we obtain redshift estimates with a typical precision of 12% for all the 2MASS XSC galaxies that lack spectroscopy. In addition, we have made an early effort towards probing the entire 3D sky beyond 2MASS, by pairing up WISE with SuperCOSMOS and training the ANNz on GAMA redshift data reaching currently to z_med~0.2. This has yielded photo-z accuracies comparable to those in the 2MPZ. These all-sky photo-z catalogs, with a median z~0.1 for the 2MPZ, and significantly deeper for future WISE-based samples, will be the largest and most complete of their kind for the foreseeable future.


title: Strong-Field Many-Body Physics and the Giant Enhancement in the High-Harmonic Spectrum of Xenon

abstract: We resolve an open question about the origin of the giant enhancement in the high-harmonic generation (HHG) spectrum of atomic xenon around 100 eV. By solving the many-body time-dependent Schr\"odinger equation with all orbitals in the 4d, 5s, and 5p shells active, we demonstrate the enhancement results truly from collective many-body excitation induced by the returning photoelectron via two-body interchannel interactions. Without the many-body interactions, which promote a 4d electron into the 5p vacancy created by strong-field ionization, no collective excitation and no enhancement in the HHG spectrum exist.


title: New detections of radio minihalos in cool cores of galaxy clusters

abstract: Cool cores of some galaxy clusters exhibit faint radio "minihalos." Their origin is unclear; their study has been limited by their small number. We undertook a systematic search for minihalos in a large sample of X-ray luminous clusters with high-quality radio data. In this paper, we report four new minihalos (A 478, ZwCl 3146, RXJ 1532.9+3021 and A 2204), and five candidates, found in the reanalyzed archival Very Large Array observations. The radio luminosities of our minihalos and candidates are in the range $10^{23-25}$ W Hz$^{-1}$ at 1.4 GHz, consistent with this type of radio sources. Their sizes (40-160 kpc in radius) are somewhat smaller than those of the previously known minihalos. We combine our new detections with previously known minihalos, obtaining a total sample of 21 objects, and briefly compare the cluster radio properties to the average X-ray temperature and the total masses estimated from Planck. We find that nearly all clusters hosting minihalos are hot and massive. Beyond that, there is no clear correlation between the minihalo radio power and cluster temperature or mass (in contrast with the giant radio halos found in cluster mergers, whose radio luminosity correlates with the cluster mass). Chandra X-ray images indicate gas sloshing in the cool cores of most of our clusters, with minihalos contained within the sloshing regions in many of them. This supports the hypothesis that radio-emitting electrons are reaccelerated by sloshing. Advection of relativistic electrons by the sloshing gas may also play a role in the formation of the less-extended minihalos.


title: Efficiency Gains from Using Auxiliary Variables in Imputation

abstract: Imputation models sometimes use auxiliary variables that, though not part of the planned analysis, can improve the accuracy of imputed values and the efficiency of point estimates. A recent article, using evidence from simulations, argued that the use of auxiliary variables in imputation did not improve efficiency. We review the simulation results and find that the use of auxiliary variables did improve efficiency; under some conditions the efficiency gain was equivalent to increasing the sample size by a quarter. We give an example from our own research where the efficiency gained from auxiliary variables was equivalent to increasing the sample size by three quarters, and pushed some estimates from statistical insignificance to significance. For auxiliary variables to make a difference, there must be a lot of missing data, some estimates must be near the border of significance, and the auxiliary variables must be excellent predictors of the missing values.


title: Attractor Control Using Machine Learning

abstract: We propose a general strategy for feedback control design of complex dynamical systems exploiting the nonlinear mechanisms in a systematic unsupervised manner. These dynamical systems can have a state space of arbitrary dimension with finite number of actuators (multiple inputs) and sensors (multiple outputs). The control law maps outputs into inputs and is optimized with respect to a cost function, containing physics via the dynamical or statistical properties of the attractor to be controlled. Thus, we are capable of exploiting nonlinear mechanisms, e.g. chaos or frequency cross-talk, serving the control objective. This optimization is based on genetic programming, a branch of machine learning. This machine learning control is successfully applied to the stabilization of nonlinearly coupled oscillators and maximization of Lyapunov exponent of a forced Lorenz system. We foresee potential applications to most nonlinear multiple inputs/multiple outputs control problems, particulary in experiments.


title: Site-Selective Electronic Correlation in alpha-Plutonium Metal

abstract: An understanding of the phase diagram of elemental plutonium (Pu) must include both the effects of the strong directional bonding and the high density of states of the Pu 5f electrons, as well as how that bonding weakens under the influence of strong electronic correlations. We present for the first time electronic-structure calculations of the full 16-atom per unit cell alpha-phase structure within the framework of density functional theory (DFT) together with dynamical mean-field theory (DMFT). Our calculations demonstrate that Pu atoms sitting on different sites within the alpha-Pu crystal structure have a strongly varying site dependence of the localization-delocalization correlation effects of their 5f electrons and a corresponding effect on the bonding and electronic properties of this complicated metal. In short, alpha-Pu has the capacity to simultaneously have multiple degrees of electron localization/delocalization of Pu 5f electrons within a pure single-element material.


title: On the $p$-integrality of $A$-hypergeometric series

abstract: Let $A$ be a set of $N$ vectors in ${\mathbb Z}^n$ and let $v$ be a vector in ${\mathbb C}^N$ that has minimal negative support for $A$. Such a vector $v$ gives rise to a formal series solution of the $A$-hypergeometric system with parameter $\beta = Av$. If $v$ lies in ${\mathbb Q}^n$, then this series has rational coefficients. Let $p$ be a prime number. We characterize those $v$ whose coordinates are rational, $p$-integral, and lie in the closed interval $[-1,0]$ for which the corresponding normalized series solution has $p$-integral coefficients.


title: An Introduction to QBism with an Application to the Locality of Quantum Mechanics

abstract: We give an introduction to the QBist interpretation of quantum mechanics. We note that it removes the paradoxes, conundra, and pseudo-problems that have plagued quantum foundations for the past nine decades. As an example, we show in detail how it eliminates quantum "nonlocality".


title: GRB 130427A: a Nearby Ordinary Monster

abstract: Long-duration Gamma-Ray Bursts (GRBs) are an extremely rare outcome of the collapse of massive stars, and are typically found in the distant Universe. Because of its intrinsic luminosity ($L\sim 3 \times 10^{53}$ erg s$^{-1}$) and its relative proximity ($z=0.34$), GRB 130427A was a unique event that reached the highest fluence observed in the gamma-ray band. Here we present a comprehensive multiwavelength view of GRB 130427A with Swift, the 2-m Liverpool and Faulkes telescopes and by other ground-based facilities, highlighting the evolution of the burst emission from the prompt to the afterglow phase. The properties of GRB 130427A are similar to those of the most luminous, high-redshift GRBs, suggesting that a common central engine is responsible for producing GRBs in both the contemporary and the early Universe and over the full range of GRB isotropic energies.


title: The kinetic Sunyaev-Zel'dovich effect from the diffuse gas in the Local Group

abstract: Since the Local Group (LG) of galaxies moves with a bulk velocity with respect to the cosmic microwave background radiation (CMB), free electrons in its gaseous halo should imprint large-scale non-primordial temperature shifts in the CMB via the kinetic Sunyaev-Zel'dovich (kSZ) effect. By modeling the distribution of gas in the LG halo and using its inferred velocity with respect to the CMB, we calculate the resulting kSZ signal from the diffuse LG medium. We find that it is dominated by a hot spot $\sim 10^\circ$ in size in the direction of M31, where the optical depth of free electrons is the greatest. By performing a correlation analysis, we find no statistical evidence that the kSZ signal from model of the LG halo is embedded in the CMB temperature map measured by the Planck satellite. We constrain the amount of mass in the LG halo by limiting the kSZ temperature shift around the hot spot to be smaller than the observed temperature shift in the Planck map. We find the tightest constraints for models where the halo mass is highly concentrated, with the mass limited to roughly $2.5-5\times 10^{12}$M$_\odot$, but note that halos with such high concentrations are rare.


title: Positive isotropic curvature and self-duality in dimension 4

abstract: We study a positivity condition for the curvature of oriented Riemannian 4-manifolds: The half-$PIC$ condition. It is a slight weakening of the positive isotropic curvature ($PIC$) condition introduced by M. Micallef and J. Moore. We observe that the half-$PIC$ condition is preserved by the Ricci flow and satisfies a maximality property among all Ricci flow invariant positivity conditions on the curvature of oriented 4-manifolds. We also study some geometric and topological aspects of half-$PIC$ manifolds.


title: Cylinders in singular del Pezzo surfaces

abstract: For each del Pezzo surface $S$ with du Val singularities, we determine whether it admits a $(-K_S)$-polar cylinder or not. If it allows one, then we present an effective $\mathbb{Q}$-divisor $D$ that is $\mathbb{Q}$-linearly equivalent to $-K_S$ and such that the open set $S\setminus\mathrm{Supp}(D)$ is a cylinder. As a corollary, we classify all the del Pezzo surfaces with du Val singularities that admit nontrivial $\mathbb{G}_a$-actions on their affine cones defined by their anticanonical divisors.


title: The Nuclear Cluster of the Milky Way: Total Mass and Luminosity

abstract: Here we present the fundamental properties of the nuclear cluster of the Milky Way. First, we derive its structural properties by constructing a density map of the central 1000'' using extinction-corrected star counts. We can describe the data with a two-component model built from Sersic profiles. The inner nearly spherical component is the nuclear cluster. The outer, strongly flattened component can be identified with the stellar component of the circumnuclear zone. Second, we enlarge the radius inside which detailed dynamics are available from $1\,$pc to $4\,$pc. We use more than 10000 individual proper motions and more than 2700 radial velocities. We determine the cluster mass by means of isotropic spherical Jeans modeling. We get a nuclear cluster mass within 100'' of $M_{100''}=(6.11 \pm 0.52|_{\mathrm{fix} R_0}\pm 0.97|_{R_0}) \times 10^6 $ M$_{\odot}$, which corresponds to a total cluster mass of M$_{\mathrm{NC}}=(13.08 \pm 2.51|_{\mathrm{fix} R_0}\pm 2.08|_{R_0}) \times 10^6 $ M$_{\odot}$. By combination of our mass with the flux we calculate $M/L=0.50 \pm 0.12 M_{\odot}/L_{\odot,\mathrm{Ks}}$ for the central 100''. That is broadly consistent with a Chabrier IMF. With its mass and a luminosity of M$_{\mathrm{Ks}}=-15.30\pm0.26$ the nuclear cluster is a bright and massive specimen with a typical size.


title: Multi-faced Black Janus and Entanglement

abstract: To the second order of the deformation parameter, we construct the black brane solutions, which are dual to the multiple interface deformations of conformal field theories. We compute their thermodynamic quantities from the gravity and the field theory sides, and find a precise agreement, confirming the validity of the AdS/CFT correspondence. The correlation of two separate interfaces induces a Casimir energy and an interesting form of correlated entropy contribution. We comment on the properties of the interface lattice system too.


title: Del Pezzo surfaces and local inequalities

abstract: I prove new local inequality for divisors on smooth surfaces, describe its applications, and compare it to a similar local inequality that is already known by experts.


title: Cosmological solutions in F(T) gravity with the presence of spinor fields

abstract: The presence of spinor fields is considered in the framework of some extensions of teleparallel gravity, where the Weitzenb\"ock connection is assumed. Some well known models as the Chaplygin gas and its generalizations are reconstructed in terms of a spinor field in the framework of teleparallel gravity. In addition, the {\Lambda}CDM model is also realized with the presence of a spinor field where a simple self-intercating term is considered and the corresponding action is reconstructed. Other cosmological solutions and the reconstruction of the gravitational action in terms of the scalar torsion is studied.


title: From a melt of rings to chromosome territories: The role of topological constraints in genome folding

abstract: We review pro and contra of the hypothesis that generic polymer properties of topological constraints are behind many aspects of chromatin folding in eukaryotic cells. For that purpose, we review, first, recent theoretical and computational findings in polymer physics related to concentrated, topologically-simple (unknotted and unlinked) chains or a system of chains. Second, we review recent experimental discoveries related to genome folding. Understanding in these fields is far from complete, but we show how looking at them in parallel sheds new light on both.


title: Flat-relative optimal extraction. A quick and efficient algorithm for stabilised spectrographs

abstract: Optimal extraction is a key step in processing the raw images of spectra as registered by two-dimensional detector arrays to a one-dimensional format. Previously reported algorithms reconstruct models for a mean one-dimensional spatial profile to assist a properly weighted extraction. We outline a simple optimal extraction algorithm including error propagation, which is very suitable for stabilised, fibre-fed spectrographs and does not model the spatial profile shape. A high signal-to-noise, master-flat image serves as reference image and is directly used as an extraction profile mask. Each extracted spectral value is the scaling factor relative to the cross-section of the unnormalised master-flat which contains all information about the spatial profile as well as pixel-to-pixel variations, fringing, and blaze. The extracted spectrum is measured relative to the flat spectrum. Using echelle spectra of the HARPS spectrograph we demonstrate a competitive extraction performance in terms of signal-to-noise and show that extracted spectra can be used for high precision radial velocity measurement. Pre- or post-flat-fielding of the data is not necessary, since all spectrograph inefficiencies inherent to the extraction mask are automatically accounted for. Also the reconstruction of the mean spatial profile by models is not needed, thereby reducing the number of operations to extract spectra. Flat-relative optimal extraction is a simple, efficient, and robust method that can be applied easily to stabilised, fibre-fed spectrographs.


title: The Range of Variation of the Mass of the Most Massive Star in Stellar Clusters Derived from 35 Million Monte Carlo Simulations

abstract: A growing fraction of Simple Stellar Population (SSP) models, in an aim to create more realistic simulations capable of including stochastic variation in their outputs, begin their simulations with a distribution of discrete stars following a power-law function of masses. Careful attention is needed to create a correctly sampled Initial Mass Function (IMF) and in this contribution we provide a solid mathematical method called MASSCLEAN IMF Sampling for doing so. We then use our method to perform $10$ $million$ MASSCLEAN Monte Carlo stellar cluster simulations to determine the most massive star in a mass distribution as a function of the total mass of the cluster. We find a maximum mass range is predicted, not a single maximum mass. This maximum mass range is (a) dependent on the total mass of the cluster and (b) independent of an upper stellar mass limit, $M_{limit}$, for $unsaturated$ clusters and comes out naturally using our IMF sampling method. We then turn our analysis around, now starting with our new $25$ $million$ simulated cluster database, to constrain the highest mass star from the observed integrated colors of a sample of 40 low-mass LMC stellar clusters of known age and mass. Finally, we present an analytical description of the maximum mass range of the most massive star as a function of the cluster's total mass, and present a new $M_{max}-M_{cluster}$ relation.


title: Atomic representations of Cuntz algebras

abstract: To a representation of $\O_N$ (the Cuntz algebra with $N$ generators) we associate a projection valued measure and we study the case when this measure has atoms. The main technical tool are the spaces invariant for all the operators $S_i^*$. We classify the purely atomic representations and find when such representations are permutative. Applications include: wavelet representations, representations generated by finitely correlated states, representation associated to Hadamard triples (or fractal spectral measures) and representations associated to generalized Walsh bases.


title: HERWIRI2: Exponentiated Electroweak Corrections in a Hadronic Event Generator

abstract: Reaching the desired precision level for W and Z processes at the LHC will require a mixture of higher-order QCD and electroweak corrections. HERWIRI2 is a step in implementing QED x QCD exponentiation in a hadronic event generator. This program implements leading electroweak corrections and coherent exclusive exponentiation in a HERWIG environment. We discuss the status of the program, recent tests, and future developments.


title: Instability of magnetic fields in electroweak plasma driven by neutrino asymmetries

abstract: The magnetohydrodynamics (MHD) is modified to incorporate the parity violation in the Standard Model leading to a new instability of magnetic fields in the electroweak plasma in the presence of nonzero neutrino asymmetries. The main ingredient for such a modified MHD is the antisymmetric part of the photon polarization tensor in plasma, where the parity violating neutrino interaction with charged leptons is present. We calculate this contribution to the polarization tensor connected with the Chern-Simons term in effective Lagrangian of the electromagnetic field. The general expression for such a contribution which depends on the temperature and the chemical potential of plasma as well as on the photon's momentum is derived. The instability of a magnetic field driven by the electron neutrino asymmetry for the $\nu$-burst during the first second of a supernova explosion can amplify a seed magnetic field of a protostar, and, perhaps, can explain the generation of strongest magnetic fields in magnetars. The growth of a cosmological magnetic field driven by the neutrino asymmetry density $\Delta n_{\nu}=n_{\nu} - n_{\bar{\nu}}\neq 0$ is provided by a lower bound on $|\xi_{\nu_e}|=|\mu_{\nu_e}|/T$ which is consistent with the well-known Big Bang nucleosynthesis (upper) bound on neutrino asymmetries in a hot universe plasma.


title: Study of Possible Ultraviolet Zero of the Beta Function in Gauge Theories with Many Fermions

abstract: We study the possibility of an ultraviolet (UV) zero in the $n$-loop beta function of U(1) and non-Abelian gauge theories with $N_f$ fermions for large $N_f$. The effect of scheme transformations on the coefficients of different powers of $N_f$ in the $n$-loop term in the beta function is calculated. A general scheme-independent criterion is given for determining whether or not the $n$-loop beta function has a UV zero for large $N_f$. We compare the results with exact integral representations of the leading terms in the beta functions for the respective Abelian and non-Abelian theories in the limit $N_f \to \infty$ limit with $N_f \alpha$ finite. As part of this study, new analytic and numerical results are presented for certain coefficients, denoted $b_{n,n-1}$, that control the large-$N_f$ behavior at $n$-loop order in the beta function. We also investigate various test functions incorporating a power-law and essential UV zero in the beta function and determine their manifestations in series expansions in powers of coupling and in powers of $1/N_f$.


title: Quantum Hamiltonian Learning Using Imperfect Quantum Resources

abstract: Identifying an accurate model for the dynamics of a quantum system is a vexing problem that underlies a range of problems in experimental physics and quantum information theory. Recently, a method called quantum Hamiltonian learning has been proposed by the present authors that uses quantum simulation as a resource for modeling an unknown quantum system. This approach can, under certain circumstances, allow such models to be efficiently identified. A major caveat of that work is the assumption of that all elements of the protocol are noise-free. Here, we show that quantum Hamiltonian learning can tolerate substantial amounts of depolarizing noise and show numerical evidence that it can tolerate noise drawn from other realistic models. We further provide evidence that the learning algorithm will find a model that is maximally close to the true model in cases where the hypothetical model lacks terms present in the true model. Finally, we also provide numerical evidence that the algorithm works for non-commuting models. This work illustrates that quantum Hamiltonian learning can be performed using realistic resources and suggests that even imperfect quantum resources may be valuable for characterizing quantum systems.


title: Automated Feature Identification in Web Applications

abstract: Market-driven software intensive product development companies have been more and more experiencing the problem of feature expansion over time. Product managers face the challenge of identifying and locating the high value features in an application and weeding out the ones of low value from the next releases. Currently, there are few methods and tools that deal with feature identification and they address the problem only partially. Therefore, there is an urgent need of methods and tools that would enable systematic feature reduction to resolve issues resulting from feature creep. This paper presents an approach and an associated tool to automate feature identification for web applications. For empirical validation, a multiple case study was conducted using three well known web applications: Youtube, Google and BBC. The results indicate that there is a good potential for automating feature identification in web applications.


title: The SOPHIE search for northern extrasolar planets: VI. Three new hot Jupiters in multi-planet extrasolar systems

abstract: We present high-precision radial-velocity measurements of three solar-type stars: HD 13908, HD 159243, and HIP 91258. The observations were made with the SOPHIE spectrograph at the 1.93-m telescope of Observatoire de Haute-Provence (France). They show that these three bright stars host exoplanetary systems composed of at least two companions. HD 13908 b is a planet with a minimum mass of 0.865+-0.035 Mjup, on a circular orbit with a period of 19.382+-0.006 days. There is an outer massive companion in the system with a period of 931+-17 days, e = 0.12+-0.02, and a minimum mass of 5.13+-0.25 Mjup. The star HD 159243, also has two detected companions with respective masses, periods, and eccentricities of Mp = 1.13+-0.05 and 1.9+-0.13 Mjup, $P$ = 12.620+-0.004 and 248.4+-4.9 days, and e = 0.02+-0.02 and 0.075+-0.05. Finally, the star HIP 91258 has a planetary companion with a minimum mass of 1.068+-0.038 Mjup, an orbital period of 5.0505+-0.0015 days, and a quadratic trend indicating an outer planetary or stellar companion that is as yet uncharacterized. The planet-hosting stars HD 13908, HD 159243, and HIP 91258 are main-sequence stars of spectral types F8V, G0V, and G5V, respectively, with moderate activity levels. HIP 91258 is slightly over-metallic, while the two other stars have solar-like metallicity. The three systems are discussed in the frame of formation and dynamical evolution models of systems composed of several giant planets.


title: Student's Video Production as Formative Assessment

abstract: Learning assessments are subject of discussions that envolve theoretical and practical approaches. To measure learning in physics by high school students, either qualitatively or quantitatively, is a process in which it should be possible to identify not only the concepts and contents students failed to achieve but also the reasons of the failure. We propose that students' video production offers a very effective formative assessment to teachers: as a formative assessment, it produces information that allows the understanding of where and when the learning process succeeded or failed, of identifying, as a subject or as a group, the defficiencies or misunderstandings related to the theme under analysis and their interpretation by students, and it provides also a different kind of assessment, related to some other life skills, like the ability to carry a project to its conclusion and to work cooperatively. In this paper, we describe the use of videos produced by high school students as an assessment resource. The students were asked to prepare a short video, which was then presented to the whole group and discussed. The videos reveal aspects of students' difficulties that usually do not appear in formal assessments such as tests and questionnaires. After the use of the videos as an component of classroom assessments and the use of the discussions to rethink learning activities in the group, the videos were analysed and classified in various categories. This analysis showed a strong correlation between the technical quality of the video and the content quality of the students argumentation. Also, it was shown that the students don't prepare their video based on quick and easy production, they usually choose forms of video production that require careful planning and implementation, and this reflects directly on the overall quality of the video and of the learning process.


title: Topology and purity for torsors

abstract: We study the homotopy theory of the classifying space of the complex projective linear groups to prove that purity fails for $PGL_p$-torsors on regular noetherian schemes when $p$ is a prime. Extending our previous work when $p=2$, we obtain a negative answer to a question of Colliot-Th\'el\`ene and Sansuc, for all $PGL_p$. We also give a new example of the failure of purity for the cohomological filtration on the Witt group, which is the first example of this kind of a variety over an algebraically closed field.


title: A Study of Error Variance Estimation in Lasso Regression

abstract: Variance estimation in the linear model when $p > n$ is a difficult problem. Standard least squares estimation techniques do not apply. Several variance estimators have been proposed in the literature, all with accompanying asymptotic results proving consistency and asymptotic normality under a variety of assumptions. It is found, however, that most of these estimators suffer large biases in finite samples when true underlying signals become less sparse with larger per element signal strength. One estimator seems to be largely neglected in the literature: a residual sum of squares based estimator using Lasso coefficients with regularisation parameter selected adaptively (via cross-validation). In this paper, we review several variance estimators and perform a reasonably extensive simulation study in an attempt to compare their finite sample performance. It would seem from the results that variance estimators with adaptively chosen regularisation parameters perform admirably over a broad range of sparsity and signal strength settings. Finally, some intial theoretical analyses pertaining to these types of estimators are proposed and developed.


title: The index of an Eisenstein ideal and multiplicity one

abstract: Mazur's fundamental work on Eisenstein ideals of prime level has a variety of arithmetic applications. In this article, we generalize some of his work to square-free level. More specifically, we attempt to compute the index of an Eisenstein ideal and the dimension of the m-torsion of the modular Jacobian variety, where m is an Eisenstein maximal ideal. In many cases, the dimension of the m-torsion is 2, in other words, a multiplicity one theorem holds.


title: Characterizing derivations for any nest algebras on Banach spaces by their behaviors at an injective operator

abstract: Let ${\mathcal N}$ be a nest on a complex Banach space $X$ and let $\mbox{ Alg}{\mathcal N}$ be the associated nest algebra. We say that an operator $Z\in \mbox{ Alg}{\mathcal N}$ is an all-derivable point of $\mbox{ Alg}{\mathcal N}$ if every linear map $\delta$ from $\mbox{ Alg}{\mathcal N}$ into itself derivable at $Z$ (i.e. $\delta$ satisfies $\delta(A)B+A\delta(B)=\delta(Z)$ for any $A,B \in \mbox{ Alg}{\mathcal N}$ with $AB=Z$) is a derivation. In this paper, it is shown that every injective operator and every operator with dense range in $\mbox{Alg}{\mathcal N}$ are all-derivable points of $\mbox{Alg}{\mathcal N}$ without any additional assumption on the nest.


title: Free probability, planar algebras, and the multicolored Guionnet-Jones-Shlyakhtenko construction

abstract: This is the graduate thesis of Michael Hartglass at UC Berkeley. A multicolored variant of the Guionnet-Jones-Shlyakhtenko construction is explored.


title: Charged Leptons

abstract: This is the report of the Intensity Frontier Charged Lepton Working Group of the 2013 Community Summer Study "Snowmass on the Mississippi", summarizing the current status and future experimental opportunities in muon and tau lepton studies and their sensitivity to new physics. These include searches for charged lepton flavor violation, measurements of magnetic and electric dipole moments, and precision measurements of the decay spectrum and parity-violating asymmetries.


title: Nonlinear travelling waves on non-Euclidean spaces

abstract: We study travelling wave solutions, that is, solutions of the form $v(t, x) = e^{i\lambda t}u(g(t)x)$, to nonlinear Schr\"odinger and Klein-Gordon equations on Riemannian manifolds, both compact and non-compact ones, with emphasis on the NLKG. Here $g(t)$ represents a one-parameter family of isometries generated by a Killing field $X$ and a case of particular interest is when $X$ has length $\leq 1$, which leads in certain settings to hypoelliptic operators with loss of at least one derivative. In the compact case, we establish existence of travelling wave solutions via "energy" minimization methods and prove that at least compact isotropic manifolds have \emph{genuinely} travelling waves. We establish certain sharp regularity estimates on low dimensional spheres that improve results in ~\cite{T1} and carry out the subelliptic analysis for NLKG on spheres of higher dimensions utilizing their homogeneous coset space properties. Such subelliptic phenomenon have no parallel in the setting of flat spaces. We will also study related phenomenon on complete noncompact manifolds with certain symmetry assumptions using concentration-compactness type arguments.


title: Topological Equivalence and Curvature Convergence: Bezier Surface Approximation

abstract: A set of control points can determine a Bezier surface and a triangulated surface simultaneously. We prove that the triangulated surface becomes homeomorphic and ambient isotopic to the Bezier surface via subdivision. We also show that the total Gaussian curvature of the triangulated surface converges to the total Gaussian curvature of the Bezier surface.


title: Uniqueness of diffusion operators and capacity estimates

abstract: Let $\Omega$ be a connected open subset of $\Ri^d$. We analyze $L_1$-uniqueness of real second-order partial differential operators $H=-\sum^d_{k,l=1}\partial_k\,c_{kl}\,\partial_l$ and $K=H+\sum^d_{k=1}c_k\,\partial_k+c_0$ on $\Omega$ where $c_{kl}=c_{lk}\in W^{1,\infty}_{\rm loc}( \Omega), c_k\in L_{\infty,{\rm loc}}(\Omega)$, $c_0\in L_{2,{\rm loc}}(\Omega)$ and $C(x)=(c_{kl}(x))>0$ for all $x\in\Omega$. Boundedness properties of the coefficients are expressed indirectly in terms of the balls $B(r)$ associated with the Riemannian metric $C^{-1}$ and their Lebesgue measure $|B(r)|$. \noindent\hspace{10mm}First we establish that if the balls $B(r)$ are bounded, the T\"acklind condition $\int^\infty_Rdr\,r(\log|B(r)|)^{-1}=\infty$ is satisfied for all large $R$ and $H$ is Markov unique then $H$ is $L_1$-unique. If, in addition, $C(x)\geq \kappa\, (c^{T}\!\otimes\, c)(x)$ for some $\kappa>0$ and almost all $x\in\Omega$, $\divv c\in L_{\infty,{\rm loc}}(\Omega)$ is upper semi-bounded and $c_0$ is lower semi-bounded then $K$ is also $L_1$-unique. \noindent\hspace{10mm}Secondly, if the $c_{kl}$ extend continuously to functions which are locally bounded on $\partial\Omega$ and if the balls $B(r)$ are bounded we characterize Markov uniqueness of $H$ in terms of local capacity estimates and boundary capacity estimates. For example, $H$ is Markov unique if and only if for each bounded subset $A$ of $\overline\Omega$ there exist $\eta_n \in C_c^\infty(\Omega)$ satisfying $\lim_{n\to\infty} \|\one_A\Gamma(\eta_n)\|_1 = 0$, where $\Gamma(\eta_n)=\sum^d_{k,l=1}c_{kl}\,(\partial_k\eta_n)\,(\partial_l\eta_n)$, and $\lim_{n\to\infty}\|\one_A (\one_\Omega-\eta_n )\, \varphi\|_2 = 0$ for each $\varphi \in L_2(\Omega)$ or if and only if $\capp(\partial\Omega)=0$.


title: RESOLVE: A new algorithm for aperture synthesis imaging of extended emission in radio astronomy

abstract: We present RESOLVE, a new algorithm for radio aperture synthesis imaging of extended and diffuse emission in total intensity. The algorithm is derived using Bayesian statistical inference techniques, estimating the surface brightness in the sky assuming a priori log-normal statistics. RESOLVE not only estimates the measured sky brightness in total intensity, but also its spatial correlation structure, which is used to guide the algorithm to an optimal reconstruction of extended and diffuse sources. For a radio interferometer, it succeeds in deconvolving the effects of the instrumental point spread function during this process. Additionally, RESOLVE provides a map with an uncertainty estimate of the reconstructed surface brightness. Furthermore, with RESOLVE we introduce a new, optimal visibility weighting scheme that can be viewed as an extension to robust weighting. In tests using simulated observations, the algorithm shows improved performance against two standard imaging approaches for extended sources, Multiscale-CLEAN and the Maximum Entropy Method.


title: Two holes in the t-J model form a bound state for any nonzero J/t

abstract: Determination of the parameter regime in which two holes in the t-J model form a bound state represents a long standing open problem in the field of strongly correlated systems. By applying and systematically improving the exact diagonalization method defined over a limited functional space (EDLFS), we show that the average distance between two holes scales as $\langle d \rangle \sim 2 (J/t)^{-1/4}$ for J/t < 0.15, therefore providing strong evidence that two holes in the t-J model form the bound state for any nonzero J/t. However, the symmetry of such bound pair in the ground state is p-wave. This state is consistent with phase separation at finite hole filling, as observed in a recent study [Maska et al, Phys. Rev. B 85, 245113 (2012)].


title: Radiative neutron capture on 16O at astrophysical energies

abstract: The total cross sections of the 16O(n,gamma)17O radiative capture process is described in the framework of the modified potential cluster model with the classification of orbital states according to Young tableaux. The possibility to describe the experimental data for the total cross sections of the neutron radiative capture on 16O at the energies from 0.025 eV to 1.5 MeV is shown on the basis of E1 and M1 processes. The transitions to the bound states of 17O in the n16O channel from 2S1/2, 2P3/2 and 2D3/2 waves of the n16O scattering have been considered. These transitions have been considered in the resonance energy region at 433 and 1000 keV.


title: Baryon Number Violation

abstract: This report, prepared for the Community Planning Study - Snowmass 2013 - summarizes the theoretical motivations and the experimental efforts to search for baryon number violation, focussing on nucleon decay and neutron-antineutron oscillations. Present and future nucleon decay search experiments using large underground detectors, as well as planned neutron-antineutron oscillation search experiments with free neutron beams are highlighted.


title: Matrix Convex Hulls of Free Semialgebraic Sets

abstract: This article resides in the realm of the noncommutative (free) analog of real algebraic geometry - the study of polynomial inequalities and equations over the real numbers - with a focus on matrix convex sets $C$ and their projections $\hat C$. A free semialgebraic set which is convex as well as bounded and open can be represented as the solution set of a Linear Matrix Inequality (LMI), a result which suggests that convex free semialgebraic sets are rare. Further, Tarski's transfer principle fails in the free setting: The projection of a free convex semialgebraic set need not be free semialgebraic. Both of these results, and the importance of convex approximations in the optimization community, provide impetus and motivation for the study of the free (matrix) convex hull of free semialgebraic sets. This article presents the construction of a sequence $C^{(d)}$ of LMI domains in increasingly many variables whose projections $\hat C^{(d)}$ are successively finer outer approximations of the matrix convex hull of a free semialgebraic set $D_p=\{X: p(X)\succeq0\}$. It is based on free analogs of moments and Hankel matrices. Such an approximation scheme is possibly the best that can be done in general. Indeed, natural noncommutative transcriptions of formulas for certain well known classical (commutative) convex hulls does not produce the convex hulls in the free case. This failure is illustrated on one of the simplest free nonconvex $D_p$. A basic question is which free sets $\hat S$ are the projection of a free semialgebraic set $S$? Techniques and results of this paper bear upon this question which is open even for convex sets.


title: T2HK: J-PARC upgrade plan for future and beyond T2K

abstract: Upgrades to the J-PARC accelerators and the neutrino experimental facility are of vital importance to the Tokai to Hyper-Kamiokande project (T2HK), which aims to explore CP asymmetry by the $\nu_e$ appearance. In this talk I present overview of the T2HK project, current status of the beam operation at J-PARC, and prospect to realize the rated 750 kW operation in coming years.


title: Non-naturally reductive Einstein metrics on the compact simple Lie group $F_4$

abstract: Based on the representation theory and the study on the involutions of compact simple Lie groups, we show that $F_4$ admits non-naturally reductive Einstein metrics.


title: Conformality at large number of fermion flavors and composite Higgs

abstract: The many open questions concerning the existence of IR and UV fixed points in gauge theories as a function of the number of fermion flavors and bare coupling are briefly reviewed and discussed. It is pointed out that only a small subset of potential IR-conformal gauge theories, i.e. theories whose IR behavior is determined by an IR fixed point, has so far been examined. It is suggested that the naturally light scalar composites that seem to appear generically as lowest states in the non-QCD-like spectrum of such theories provide a natural basis for composite Higgs models, where the composite Higgs is not a NG boson and some of the usual fine-tuning problems are evaded.


title: Texture descriptor combining fractal dimension and artificial crawlers

abstract: Texture is an important visual attribute used to describe images. There are many methods available for texture analysis. However, they do not capture the details richness of the image surface. In this paper, we propose a new method to describe textures using the artificial crawler model. This model assumes that each agent can interact with the environment and each other. Since this swarm system alone does not achieve a good discrimination, we developed a new method to increase the discriminatory power of artificial crawlers, together with the fractal dimension theory. Here, we estimated the fractal dimension by the Bouligand-Minkowski method due to its precision in quantifying structural properties of images. We validate our method on two texture datasets and the experimental results reveal that our method leads to highly discriminative textural features. The results indicate that our method can be used in different texture applications.


title: Non-Archimedean meromorphic solutions of functional equations

abstract: In this paper, we discuss meromorphic solutions of functional equations over non-Archimedean fields, and prove analogues of the Clunie lemma, Malmquist-type theorem and Mokhon'ko theorem.


title: Low-light-level four-wave mixing by quantum interference

abstract: We observed electromagnetically-induced-transparency-based four-wave mixing (FWM) in the pulsed regime at low light levels. The FWM conversion efficiency of 3.8(9)% was observed in a four-level system of cold 87Rb atoms using a driving laser pulse with a peak intensity of approximately 80 {\mu}W/cm^2, corresponding to an energy of approximately 60 photons per atomic cross section. Comparison between the experimental data and the theoretical predictions proposed by Harris and Hau [Phys. Rev. Lett. 82, 4611 (1999)] showed good agreement. Additionally, a high conversion efficiency of 46(2)% was demonstrated when applying this scheme using a driving laser intensity of approximately 1.8 mW/cm^2. According to our theoretical predictions, this FWM scheme can achieve a conversion efficiency of nearly 100% when using a dense medium with an optical depth of 500.


title: Bruckner--Garg-type results with respect to Haar null sets in $C[0,1]$

abstract: A set $\mathcal{A}\subset C[0,1]$ is \emph{shy} or \emph{Haar null } (in the sense of Christensen) if there exists a Borel set $\mathcal{B}\subset C[0,1]$ and a Borel probability measure $\mu$ on $C[0,1]$ such that $\mathcal{A}\subset \mathcal{B}$ and $\mu\left(\mathcal{B}+f\right) = 0$ for all $f \in C[0,1]$. The complement of a shy set is called a \emph{prevalent} set. We say that a set is \emph{Haar ambivalent} if it is neither shy nor prevalent. The main goal of the paper is to answer the following question: What can we say about the topological properties of the level sets of the prevalent/non-shy many $f\in C[0,1]$? The classical Bruckner--Garg Theorem characterizes the level sets of the generic (in the sense of Baire category) $f\in C[0,1]$ from the topological point of view. We prove that the functions $f\in C[0,1]$ for which the same characterization holds form a Haar ambivalent set. In an earlier paper we proved that the functions $f\in C[0,1]$ for which positively many level sets with respect to the Lebesgue measure $\lambda$ are singletons form a non-shy set in $C[0,1]$. The above result yields that this set is actually Haar ambivalent. Now we prove that the functions $f\in C[0,1]$ for which positively many level sets with respect to the occupation measure $\lambda\circ f^{-1}$ are not perfect form a Haar ambivalent set in $C[0,1]$. We show that for the prevalent $f\in C[0,1]$ for the generic $y\in f([0,1])$ the level set $f^{-1}(y)$ is perfect. Finally, we answer a question of Darji and White by showing that the set of functions $f \in C[0,1]$ for which there exists a perfect $P_f\subset [0,1]$ such that $f'(x) = \infty$ for all $x \in P_f$ is Haar ambivalent.


title: Atomic ionization of germanium due to neutrino magnetic moments

abstract: An ab initio calculation of atomic ionization of germanium (Ge) by neutrinos was carried out in the framework of multiconfiguration relativistic random phase approximation. The main goal is to provide a more accurate cross section formula than the conventional one, which is based on the free electron approximation, for searches of neutrino magnetic moments with Ge detectors whose threshold is reaching down to the sub-keV regime. Limits derived with both methods are compared, using reactor neutrino data taken with low threshold germanium detectors.


title: New result for the neutron 9be capture at astrophysical energies

abstract: The total cross sections of the n9Be --> gamma10Be radiative capture are described in the framework of the modified potential cluster model (MPCM) with the classification of orbital states according to Young tableaux at thermal and astrophysical energies.


title: The Functional Determinant and the Partition Function in Geometric Flows

abstract: We propose the use of the functional determinant of geometric operators in constructing an entropy functional associated to geometric flows. Our approach is based on the direct computation of the partition function, with a well-defined set of microstates and macrostates in the canonical ensemble. The approach is motivated by a fundamental enigma in Perelman's derivation of his famous $\mathcal{W}$-entropy. The defining feature of our entropy is that the energy of each microstate in the partition function is invariant along the associated geometric flow - a clue that could be inferred from Perelman's work. Moreover, the monotonicity of our entropy along the associated geometric flow is then a natural result in the statistical mechanics framework. While we will not argue in a completely rigorous manner, we will use the formalism to derive an explicit formula for an entropy associated to conformal flows on a closed surface based on the Polyakov formula for the determinant of the Laplacian. We also discuss possible extensions of our results to more general operators and manifolds.


title: Constraining Non-thermal and Thermal properties of Dark Matter

abstract: We describe the evolution of Dark Matter (DM) abundance from the very onset of its creation from inflaton decay under the assumption of an instantaneous reheating. Based on the initial conditions such as the inflaton mass and its decay branching ratio to DM, reheating temperature, and the DM mass and interaction rate with the thermal bath, the DM particles can either thermalize (fully/partially) with the primordial bath or remain non-thermal throughout their evolution history. In the thermal case, the final abundance is set by the standard freeze-out mechanism for large annihilation rates, irrespective of the initial conditions. For smaller annihilation rates, it can be set by the freeze-in mechanism, also independent of the initial abundance, provided it is small to begin with. For even smaller interaction rates, the DM decouples while being non-thermal, and the relic abundance will be essentially set by the initial conditions. We put model-independent constraints on the DM mass and annihilation rate from over-abundance by exactly solving the relevant Boltzmann equations, and identify the thermal freeze-out, freeze-in and non-thermal regions of the allowed parameter space. We highlight a generic fact that inflaton decay to DM inevitably leads to an overclosure of the Universe for a large range of DM parameter space, and thus poses a stringent constraint that must be taken into account while constructing models of DM. For the thermal DM region, we also show the complementary constraints from indirect DM search experiments, Big Bang Nucleosynthesis, Cosmic Microwave Background, Planck measurements, and theoretical limits due to the unitarity of S-matrix. For the non-thermal DM scenario, we show the allowed parameter space in terms of the inflaton and DM masses for a given reheating temperature, and compute the comoving free-streaming length to identify the hot, warm and cold DM regimes.


title: Stationary resonances of rapidly-rotating Kerr black holes

abstract: The Klein-Gordon equation for a massive scalar field in the background of a rapidly-rotating Kerr black hole is studied analytically. In particular, we derive a simple formula for the stationary (marginally-stable) resonances of the field in the black-hole spacetime. The analytically derived formula is shown to agree with direct numerical computations of the resonances. Our results provide an upper bound on the instability regime of rapidly-rotating Kerr black holes to massive scalar perturbations.


title: Simple, locally finite dimensional Lie algebras in positive characteristic

abstract: We prove two structure theorems for simple, locally finite dimensional Lie algebras over an algebraically closed field of characteristic $p$ which give sufficient conditions for the algebras to be of the form $[R^{(-)}, R^{(-)}] / (Z(R) \cap [R^{(-)}, R^{(-)}])$ or $[K(R, *), K(R, *)]$ for a simple, locally finite dimensional associative algebra $R$ with involution $*$. The first proves that a condition we introduce, known as locally nondegenerate, along with the existence of an ad-nilpotent element suffice. The second proves that a uniformly ad-integrable Lie algebra is of this type if the characteristic of the ground field is sufficiently large. Lastly we construct a simple, locally finite dimensional associative algebra $R$ with involution $*$ such that $K(R, *) \ne [K(R, *), K(R, *)]$ to demonstrate the necessity of considering the commutator in the first two theorems.


title: A Solution of the Hubbard Model

abstract: We report a ground-state solution for the two-dimensional fermionic Hubbard model, which is obtained via a numerical variational method. The two ingredients in this approach are tensor network states and the time-evolving block decimation. We easily handle the horizontal hopping in the Hamiltonian, and we proceed further to observe the fermion-exchange effect caused by the vertical hopping. By requiring no divergence and no convergence to zero for the ground state, we successively determine the ground-state energy per site as a function of the chemical potential and the lattice length. In addition, we observe saturation in the behavior of the ground-state energy as the lattice length increases.


title: Robust Estimation under Heavy Contamination using Enlarged Models

abstract: In data analysis, contamination caused by outliers is inevitable, and robust statistical methods are strongly demanded. In this paper, our concern is to develop a new approach for robust data analysis based on scoring rules. The scoring rule is a discrepancy measure to assess the quality of probabilistic forecasts. We propose a simple way of estimating not only the parameter in the statistical model but also the contamination ratio of outliers. Estimating the contamination ratio is important, since one can detect outliers out of the training samples based on the estimated contamination ratio. For this purpose, we use scoring rules with an extended statistical models, that is called the enlarged models. Also, the regression problems are considered. We study a complex heterogeneous contamination, in which the contamination ratio of outliers in the dependent variable may depend on the independent variable. We propose a simple method to obtain a robust regression estimator under heterogeneous contamination. In addition, we show that our method provides also an estimator of the expected contamination ratio that is available to detect the outliers out of training samples. Numerical experiments demonstrate the effectiveness of our methods compared to the conventional estimators.


title: Formulation for zero mode of Bose-Einstein condensate beyond Bogoliubov approximation

abstract: It is shown for the Bose-Einstein condensate of cold atomic system that the new unperturbed Hamiltonian, which includes not only the first and second powers of the zero mode operators but also the higher ones, determines a unique and stationary vacuum at zero temperature. From the standpoint of quantum field theory, it is done in a consistent manner that the canonical commutation relation of the field operator is kept. In this formulation, the condensate phase does not diffuse and is robust against the quantum fluctuation of the zero mode. The standard deviation for the phase operator depends on the condensed atom number with the exponent of $-1/3$, which is universal for both homogeneous and inhomogeneous systems.


title: Reducing beam hardening effects and metal artefacts using Medipix3RX: With applications from biomaterial science

abstract: This paper discusses methods for reducing beam hardening effects using spectral data for biomaterial applications. A small-animal spectral scanner operating in the diagnostic energy range was used. We investigate the use of photon-processing features of the Medipix3RX ASIC in reducing beam hardening and associated artefacts. A fully operational charge summing mode was used during the imaging routine. We present spectral data collected for metal alloy samples, its analysis using algebraic 3D reconstruction software and volume visualisation using a custom volume rendering software. Narrow high energy acquisition using the photon-processing detector revealed substantial reduction in beam hardening effects and metal artefacts.


title: Dynamic Partitioning-based JPEG Decompression on Heterogeneous Multicore Architectures

abstract: With the emergence of social networks and improvements in computational photography, billions of JPEG images are shared and viewed on a daily basis. Desktops, tablets and smartphones constitute the vast majority of hardware platforms used for displaying JPEG images. Despite the fact that these platforms are heterogeneous multicores, no approach exists yet that is capable of joining forces of a system's CPU and GPU for JPEG decoding. In this paper we introduce a novel JPEG decoding scheme for heterogeneous architectures consisting of a CPU and an OpenCL-programmable GPU. We employ an offline profiling step to determine the performance of a system's CPU and GPU with respect to JPEG decoding. For a given JPEG image, our performance model uses (1) the CPU and GPU performance characteristics, (2) the image entropy and (3) the width and height of the image to balance the JPEG decoding workload on the underlying hardware. Our run-time partitioning and scheduling scheme exploits task, data and pipeline parallelism by scheduling the non-parallelizable entropy decoding task on the CPU, whereas inverse cosine transformations (IDCTs), color conversions and upsampling are conducted on both the CPU and the GPU. Our kernels have been optimized for GPU memory hierarchies. We have implemented the proposed method in the context of the libjpeg-turbo library, which is an industrial-strength JPEG encoding and decoding engine. Libjpeg-turbo's hand-optimized SIMD routines for ARM and x86 constitute a competitive yardstick for the comparison to the proposed approach. Retro-fitting our method with libjpeg-turbo provided insights on the software-engineering aspects of re-engineering legacy code for heterogeneous multicores.


title: Statistics and Classification of the Microwave Zebra Patterns Associated with Solar Flares

abstract: The microwave zebra pattern (ZP) is the most interesting, intriguing, and complex spectral structure frequently observed in solar flares. A comprehensive statistical study will certainly help us to understand the formation mechanism, which is not exactly clear now. This work presents a comprehensive statistical analysis on a big sample with 202 ZP events collected from observations at the Chinese Solar Broadband Radio Spectrometer at Huairou and the Ondrejov Radiospectrograph in Czech Republic at frequencies of 1.00 - 7.60 GHz during 2000 - 2013. After investigating the parameter properties of ZPs, such as the occurrence in flare phase, frequency range, polarization degree, duration, etc., we find that the variation of zebra stripe frequency separation with respect to frequency is the best indicator for a physical classification of ZPs. Microwave ZPs can be classified into 3 types: equidistant ZP, variable-distant ZP, and growing-distant ZP, possibly corresponding to mechanisms of Bernstein wave model, whistler wave model, and double plasma resonance model, respectively. This statistical classification may help us to clarify the controversies between the existing various theoretical models, and understand the physical processes in the source regions.


title: Quadratic Twists of Elliptic Curves with 3-Selmer Rank 1

abstract: A weaker form of a 1979 conjecture of Goldfeld states that for every elliptic curve $E/\mathbb{Q}$, a positive proportion of its quadratic twists $E^{(d)}$ have rank 1. Using tools from Galois cohomology, we give criteria on E and d which force a positive proportion of the quadratic twists of E to have 3-Selmer rank 1 and global root number -1. We then give four nonisomorphic infinite families of elliptic curves $E_{m, n}$ which satisfy these criteria. Conditional on the rank part of the Birch and Swinnerton-Dyer conjecture, this verifies the aforementioned conjecture for infinitely many elliptic curves. Our elliptic curves are easy to give explicitly and we state precisely which quadratic twists d to use. Furthermore, our methods have the potential of being generalized to elliptic curves over other number fields.


title: The evolution of the power law k-essence cosmology

abstract: We investigate the evolution of the power law k-essence field in FRWL spacetime. The autonomous dynamical system and critical points are obtained. The corresponding cosmological parameters, such as $\Omega _{\phi }$ and $w_{\phi }$, are calculated at these critical points. We find it is possible to achieve an equation of state crossing through $-1$ for k-essence field. The results we obtained indicate that the power law k-essence dark energy model can be compatible with observations.


title: The busy function: a new analytic function for describing the integrated 21-cm spectral profile of galaxies

abstract: Accurate parametrization of galaxies detected in the 21-cm HI emission is of fundamental importance to the measurement of commonly used indicators of galaxy evolution, including the Tully-Fisher relation and the HI mass function. Here, we propose a new analytic function, named the 'busy function', that can be used to accurately describe the characteristic double-horn HI profile of many galaxies. The busy function is a continuous, differentiable function that consists of only two basic functions, the error function, erf(x), and a polynomial, |x|^n, of degree n >= 2. We present the basic properties of the busy function and illustrate its great flexibility in fitting a wide range of HI profiles from the Gaussian profiles of dwarf galaxies to the broad, asymmetric double-horn profiles of spiral galaxies. Applications of the busy function include the accurate and efficient parametrization of observed HI spectra of galaxies and the construction of spectral templates for simulations and matched filtering algorithms. We demonstrate the busy function's power by automatically fitting it to the HI spectra of 1000 galaxies from the HIPASS Bright Galaxy Catalog, using our own C/C++ implementation, and comparing the resulting parameters with the catalogued ones. We also present two methods for determining the uncertainties of observational parameters derived from the fit.


title: Schmidt Games and Nondense forward Orbits of certain Partially Hyperbolic Systems

abstract: Let $f: M \to M$ be a partially hyperbolic diffeomorphism with conformality on unstable manifolds. Consider a set of points with nondense forward orbit: $E(f, y) := \{ z\in M: y\notin \overline{\{f^k(z), k \in \mathbb{N}\}}\}$ for some $y \in M$. Define $E_{x}(f, y) := E(f, y) \cap W^u(x)$ for any $x\in M$. Following a method of Broderick-Fishman-Kleinbock, we show that $E_x(f,y)$ is a winning set of Schmidt games played on $W^u(x)$ which implies that $E_x(f,y)$ has full Hausdorff dimension equal to $\dim W^u(x)$. Furthermore we show that for any nonempty open set $V \subset M$, $E(f, y) \cap V$ has full Hausdorff dimension equal to $\dim M$, by constructing measures supported on $E(f, y)\cap V$ with lower pointwise dimension converging to $\dim M$ and with conditional measures supported on $E_x(f,y)\cap V$. The results can be extended to the set of points with forward orbit staying away from a countable subset of $M$.


title: Arbitrarily high order Convected Scheme solution of the Vlasov-Poisson system

abstract: The Convected Scheme (CS) is a `forward-trajectory' semi-Lagrangian method for solution of transport equations, which has been most often applied to the kinetic description of plasmas and rarefied neutral gases. In its simplest form, the CS propagates the solution by advecting the `moving cells' along their characteristic trajectories, and by remapping them on the mesh at the end of the time step. The CS is conservative, positivity preserving, simple to implement, and not subject to time step restriction to maintain stability. Recently [Y. G\"u\c{c}l\"u and W.N.G. Hitchon, 2012] a new methodology was introduced for reducing numerical diffusion, based on a modified equation analysis: the remapping error was compensated by applying small corrections to the final position of the moving cells prior to remapping. While the spatial accuracy was increased from 2nd to 4th order, the new scheme retained the important properties of the original method, and was shown to be simple and efficient for constant advection problems. Here the CS is applied to the solution of the Vlasov-Poisson system: the Vlasov equation is split into two constant advection equations, one in configuration space and one in velocity space, and high order time accuracy is achieved by proper composition of the operators. The splitting procedure enables us to use the constant advection solver, which we extend to arbitrarily high order of accuracy: a new improved procedure is given, which makes the calculation of the corrections straightforward. Focusing on periodic domains, we describe a spectrally accurate scheme based on the fast Fourier transform; the proposed implementation is strictly conservative and positivity preserving. The ability to correctly reproduce the system dynamics, as well as resolving small-scale features in the solution, is shown in classical 1D-1V test cases, both in the linear and the non-linear regimes.


title: On_the_implementation_of_the_conditions_of_Inertial_Confinement_ Fusion by bombarding the target a macro particle

abstract: The acceleration of lithium tube segments with the length one centimeter, diameter sixteen microns wall thickness one nanometer is considered. These segments are electrically charged by proton beams produced by an electron beam source. Then, they are accelerated by the traveling wave field in a spiral waveguide. The segments are next sent to a target where they are compressed by three hundred times in the longitudinal direction and compressing target radially, so the conditions for thermonuclear fusion are realized.


title: Mapping Topographic Structure in White Matter Pathways with Level Set Trees

abstract: Fiber tractography on diffusion imaging data offers rich potential for describing white matter pathways in the human brain, but characterizing the spatial organization in these large and complex data sets remains a challenge. We show that level set trees---which provide a concise representation of the hierarchical mode structure of probability density functions---offer a statistically-principled framework for visualizing and analyzing topography in fiber streamlines. Using diffusion spectrum imaging data collected on neurologically healthy controls (N=30), we mapped white matter pathways from the cortex into the striatum using a deterministic tractography algorithm that estimates fiber bundles as dimensionless streamlines. Level set trees were used for interactive exploration of patterns in the endpoint distributions of the mapped fiber tracks and an efficient segmentation of the tracks that has empirical accuracy comparable to standard nonparametric clustering methods. We show that level set trees can also be generalized to model pseudo-density functions in order to analyze a broader array of data types, including entire fiber streamlines. Finally, resampling methods show the reliability of the level set tree as a descriptive measure of topographic structure, illustrating its potential as a statistical descriptor in brain imaging analysis. These results highlight the broad applicability of level set trees for visualizing and analyzing high-dimensional data like fiber tractography output.


title: Interstellar chemistry of nitrogen hydrides in dark clouds

abstract: The aim of the present work is to perform a comprehensive analysis of the interstellar chemistry of nitrogen, focussing on the gas-phase formation of the smallest polyatomic species and in particular nitrogen hydrides. We present a new chemical network in which the kinetic rates of critical reactions have been updated based on recent experimental and theoretical studies, including nuclear spin branching ratios. Our network thus treats the different spin symmetries of the nitrogen hydrides self-consistently together with the ortho and para forms of molecular hydrogen. This new network is used to model the time evolution of the chemical abundances in dark cloud conditions. The steady-state results are analysed, with special emphasis on the influence of the overall amounts of carbon, oxygen, and sulphur. Our calculations are also compared with Herschel/HIFI observations of NH, NH$_2$, and NH$_3$ detected towards the external envelope of the protostar IRAS 16293-2422. The observed abundances and abundance ratios are reproduced for a C/O gas-phase elemental abundance ratio of $\sim0.8$, provided that the sulphur abundance is depleted by a factor larger than 2. The ortho-to-para ratio of H$_2$ in these models is $\sim10^{-3}$. Our models also provide predictions for the ortho-to-para ratios of NH$_2$ and NH$_3$ of $\sim2.3$ and $\sim0.7$ respectively. We conclude that the abundances of nitrogen hydrides in dark cloud conditions are consistent with the gas-phase synthesis predicted with our new chemical network.


title: Id\`elic class field theory for 3-manifolds

abstract: Following the analogies between 3-dimensional topology and number theory, we study an id\`elic form of class field theory for 3-manifolds. For a certain set $\mathcal{K}$ of knots in a 3-manifold $M$, we first present a local theory for each knot in $\mathcal{K}$, which is analogous to local class field theory, and then, getteing together over all knots in $\mathcal{K}$, we give an analogue of id\`elic global class field theory for an integral homology sphere $M$.


title: About_the_electrodynamic_acceleration_of_cylinder-shaped_particles

abstract: A possibility of electrodynamic acceleration of particles from the initial zero velocity to the final velocity ten kilometeres per second the acceleration length five meters is considered. After the electrostatic preacceleration particles are accelerated at the trailing edge of the voltage pulse six megavolts, which runs along the spiral turns. Accelerated particles have the diameter six microns, and length one centimeter. Because of a pointed cone at the head the particles can move in the air almost without loss of velocity, penetrating into aluminum and water as deep as ten centimeter and water one meter respectively.


title: Randall-Sundrum scenario with a small curvature and Drell-Yan process at the LHC

abstract: The Randall-Sundrum-like scenario with the small curvature kappa (RSSC model) is studied in detail in comparison with the original RS1 model. In the framework of the RSSC model, the p_t distributions for the dilepton production at the LHC are calculated. Both dielectron and dimuon events are taken into account. The important feature of calculations is the account of the widths of massive graviton excitations. For the summary statistics taken at 7 TeV (L = 5 fb^(-1)) and 8 TeV (L = 20 fb^(-1)), the exclusion limit on the 5-dimensional gravity scale M_5 is set to be 6.84 TeV at 95% C.L. For sqrt(s) = 13 TeV and integrated luminosity 30 fb^(-1), the LHC search limit is found to be 10.16 TeV. These bounds on M_5 are independent of kappa (up to powerlike corrections), provided kappa << M_5.


title: Distributed Connectivity Decomposition

abstract: We present time-efficient distributed algorithms for decomposing graphs with large edge or vertex connectivity into multiple spanning or dominating trees, respectively. As their primary applications, these decompositions allow us to achieve information flow with size close to the connectivity by parallelizing it along the trees. More specifically, our distributed decomposition algorithms are as follows: (I) A decomposition of each undirected graph with vertex-connectivity $k$ into (fractionally) vertex-disjoint weighted dominating trees with total weight $\Omega(\frac{k}{\log n})$, in $\widetilde{O}(D+\sqrt{n})$ rounds. (II) A decomposition of each undirected graph with edge-connectivity $\lambda$ into (fractionally) edge-disjoint weighted spanning trees with total weight $\lceil\frac{\lambda-1}{2}\rceil(1-\varepsilon)$, in $\widetilde{O}(D+\sqrt{n\lambda})$ rounds. We also show round complexity lower bounds of $\tilde{\Omega}(D+\sqrt{\frac{n}{k}})$ and $\tilde{\Omega}(D+\sqrt{\frac{n}{\lambda}})$ for the above two decompositions, using techniques of [Das Sarma et al., STOC'11]. Moreover, our vertex-connectivity decomposition extends to centralized algorithms and improves the time complexity of [Censor-Hillel et al., SODA'14] from $O(n^3)$ to near-optimal $\tilde{O}(m)$. As corollaries, we also get distributed oblivious routing broadcast with $O(1)$-competitive edge-congestion and $O(\log n)$-competitive vertex-congestion. Furthermore, the vertex connectivity decomposition leads to near-time-optimal $O(\log n)$-approximation of vertex connectivity: centralized $\widetilde{O}(m)$ and distributed $\tilde{O}(D+\sqrt{n})$. The former moves toward the 1974 conjecture of Aho, Hopcroft, and Ullman postulating an $O(m)$ centralized exact algorithm while the latter is the first distributed vertex connectivity approximation.