title: Quasars can be used to verify the parallax zero-point of the Tycho-Gaia Astrometric Solution

abstract: Context. The Gaia project will determine positions, proper motions, and parallaxes for more than one billion stars in our Galaxy. It is known that Gaia's two telescopes are affected by a small but significant variation of the basic angle between them. Unless this variation is taken into account during data processing, e.g. using on-board metrology, it causes systematic errors in the astrometric parameters, in particular a shift of the parallax zero-point. Previously, we suggested an early reduction of Gaia data for the subset of Tycho-2 stars (Tycho-Gaia Astrometric Solution; TGAS). Aims. We aim to investigate whether quasars can be used to independently verify the parallax zero-point already in early data reductions. This is not trivially possible as the observation interval is too short to disentangle parallax and proper motion for the quasar subset. Methods. We repeat TGAS simulations but additionally include simulated Gaia observations of quasars from ground-based surveys. All observations are simulated with basic angle variations. To obtain a full astrometric solution for the quasars in TGAS we explore the use of prior information for their proper motions. Results. It is possible to determine the parallax zero-point for the quasars with a few {\mu}as uncertainty, and it agrees to a similar precision with the zero-point for the Tycho-2 stars. The proposed strategy is robust even for quasars exhibiting significant fictitious proper motion due to a variable source structure, or when the quasar subset is contaminated with stars misidentified as quasars. Conclusions. Using prior information about quasar proper motions we could provide an independent verification of the parallax zero-point in early solutions based on less than one year of Gaia data.


title: Turbulent Cosmic Ray Reacceleration at Radio Relics and Halos in Clusters of Galaxies

abstract: Radio relics are synchrotron emission found on the periphery of galaxy clusters. From the position and the morphology, it is often believed that the relics are generated by cosmic-ray (CR) electrons accelerated at shocks through a diffusive shock acceleration (DSA) mechanism. However, some radio relics have harder spectra than the prediction of the standard DSA model. One example is observed in the cluster 1RXS J0603.3+4214, which is often called the ``Toothbrush Cluster''. Interestingly, the position of the relic is shifted from that of a possible shock. In this study, we show that these discrepancies in the spectrum and the position can be solved if turbulent (re)acceleration is very effective behind the shock. This means that for some relics turbulent reacceleration may be the main mechanism to produce high-energy electrons, contrary to the common belief that it is the DSA. Moreover, we show that for efficient reacceleration, the effective mean free path of the electrons has to be much smaller than their Coulomb mean free path. We also study the merging cluster 1E 0657-56 or the ``Bullet Cluster'', in which a radio relic has not been found at the position of the prominent shock ahead of the bullet. We indicate that a possible relic at the shock is obscured by the observed large radio halo that is generated by strong turbulence behind the shock. We propose a simple explanation if the morphological differences of radio emission among the Toothbrush, the Bullet, and the Sausage (CIZA J2242.8+5301) Clusters.


title: The physical properties of z>2 Lyman limit systems: new constraints for feedback and accretion models

abstract: We study the physical properties of a homogeneous sample of 157 optically-thick absorption line systems at redshifts ~1.8-4.4, selected from a high-dispersion spectroscopic survey of Lyman limit systems (LLSs). By means of multiple ionisation models and Bayesian techniques, we derive the posterior probability distribution functions for the density, metallicity, temperature, and dust content of the absorbing gas. We find that z>2 LLSs are highly ionised with ionisation parameters between -3<log U<-2, depending on the HI column density. LLSs are characterised by low temperatures (T<5x10^4 K) and reside in dust-poor environments. Between z~2.5-3.5, ~80% of the LLSs have physical densities between n(H)~10^-3.5-10^-2 cm^-3 for the assumed UV background, but we caution that a degeneracy between the ionisation parameter and the intensity of the radiation field prevents robust inference on the density and sizes of LLSs. Conversely, metallicity estimates are less sensitive to the assumptions behind ionisation corrections. LLSs at z>2 are characterised by a broad unimodal distribution over >4 orders of magnitude, with a peak at log Z/Zsun~-2. LLSs are metal poor, significantly less enriched than DLAs, with ~70% of the metallicity PDF below log Z/Zsun<-1.5. The median metallicity of super LLSs with log N(HI)>19 rapidly evolves with redshift, with a ten-fold increase between z~2.1-3.6 (~1.5 Gyr). Based on this sample, we find that LLSs at z=2.5-3.5 account for ~15% of all the metals produced by UV-selected galaxies. The implications for theories of cold gas accretion and metal ejection from galaxies are also discussed.


title: Non-linearity and environmental dependence of the star forming galaxies Main Sequence

abstract: Using data from four deep fields (COSMOS, AEGIS, ECDFS, and CDFN), we study the correlation between the position of galaxies in the star formation rate (SFR) versus stellar mass plane and local environment at $z<1.1$. To accurately estimate the galaxy SFR, we use the deepest available Spitzer/MIPS 24 and Herschel/PACS datasets. We distinguish group environments ( $M_{halo}\sim$10$^{12.5-14.2}$$M_{\odot}$) based on the available deep X-ray data and lower halo mass environments based on the local galaxy density. We confirm that the Main Sequence (MS) of star forming galaxies is not a linear relation and there is a flattening towards higher stellar masses ( $M_*>10^{10.4-10.6}$ $M_{\odot}$), across all environments. At high redshift ( $0.5<z<1.1$ ), the MS varies little with environment. At low redshift ( $0.15<z<0.5$ ), group galaxies tend to deviate from the mean MS towards the region of quiescence with respect to isolated galaxies and less-dense environments. We find that the flattening of the MS toward low SFR is due to an increased fraction of bulge dominated galaxies at high masses. Instead, the deviation of group galaxies from the MS at low redshift is caused by a large fraction of red disk dominated galaxies which are not present in the lower density environments. Our results suggest that above a mass threshold ( $\sim10^{10.4}-10^{10.6}$$M_{\odot}$ ) stellar mass, morphology and environment act together in driving the evolution of the SF activity towards lower level. The presence of a dominating bulge and the associated quenching processes are already in place beyond $z\sim$1. The environmental effects appear, instead, at lower redshifts and have a long time-scale.


title: One Leptoquark to Rule Them All: A Minimal Explanation for $R_{D^{(*)}}$, $R_K$ and $(g-2)_\mu$

abstract: We show that by adding a single new scalar particle to the Standard Model, a TeV-scale leptoquark with the quantum numbers of a right-handed down quark, one can explain in a natural way three of the most striking anomalies of particle physics: the violation of lepton universality in $\bar B\to\bar K\ell^+\ell^-$ decays, the enhanced $\bar B\to D^{(*)}\tau\bar\nu$ decay rates, and the anomalous magnetic moment of the muon. Constraints from other precision measurements in the flavor sector can be satisfied without fine-tuning. Our model predicts enhanced $\bar B\to\bar K^{(*)}\nu\bar\nu$ decay rates and a new-physics contribution to $B_s\!-\!\bar B_s$ mixing close to the current central fit value.


title: Shapiro Delays at the Quadrupole Order for Tests of the No-Hair Theorem Using Pulsars around Spinning Black Holes

abstract: One avenue for testing the no-hair theorem is obtained through timing a pulsar orbiting close to a black hole and fitting for quadrupolar effects on the time-of-arrival of pulses. If deviations from the Kerr quadrupole are measured, then the no-hair theorem is invalidated. To this end, we derive an expression for the light travel time delay for a pulsar orbiting in a black-hole spacetime described by the Butterworth-Ipser metric, which has an arbitrary spin and quadrupole moment. We consider terms up to the quadrupole order in the black-hole metric and derive the time-delay expression in a closed analytic form. This allows for fast computations that are useful in fitting time-of-arrival observations of pulsars orbiting close to astrophysical black holes.


title: Disorder Operators in Chern-Simons-Fermion Theories

abstract: Building on the recent progress in solving Chern-Simons-matter theories in the planar limit, we compute the scaling dimensions of a large class of disorder ("monopole") operators in $U(N)_k$ Chern-Simons-fermion theories at all 't Hooft couplings $\lambda = N/k$. We find that the lowest-dimension operator of this sort has dimension $\frac23 k^{3/2}$. We comment on the implications of these results to analyzing maps of fermionic disorder operators under 3D bosonization.


title: Mining Planet Search Data for Binary Stars: The $\psi^1$ Draconis system

abstract: Several planet-search groups have acquired a great deal of data in the form of time-series spectra of several hundred nearby stars with time baselines of over a decade. While binary star detections are generally not the goal of these long-term monitoring efforts, the binary stars hiding in existing planet search data are precisely the type that are too close to the primary star to detect with imaging or interferometry techniques. We use a cross-correlation analysis to detect the spectral lines of a new low-mass companion to $\psi^1$ Draconis A, which has a known roughly equal-mass companion at ${\sim}680$ AU. We measure the mass of $\psi^1$ Draconis C as $M_2 = 0.70 \pm 0.07 M_{\odot}$, with an orbital period of ${\sim}20$ years. This technique could be used to characterize binary companions to many stars that show large-amplitude modulation or linear trends in radial velocity data.


title: The Higgs Mass in the MSSM at two-loop order beyond minimal flavour violation

abstract: Soft supersymmetry-breaking terms provide a wealth of new potential sources of flavour violation, which are tightly constrained by precision experiments. This has posed a challenge to construct flavour models which both explain the structure of the Standard Model Yukawa couplings and also predict soft-breaking patterns that are compatible with these constraints. While such models have been studied in great detail, the impact of flavour violating soft terms on the Higgs mass at the two-loop level has been assumed to be small or negligible. In this letter, we show that large flavour violation in the up-squark sector can give a positive or negative mass shift to the SM-like Higgs of several GeV, without being in conflict with other observations. We investigate in which regions of the parameter space these effects can be expected.


title: Possible resolution of the domain wall problem in the NMSSM

abstract: We discuss a possibility that the domain wall problem in the next-to-minimal supersymmetric standard model is alleviated without introducing a small explicit $Z_3$ breaking term by analyzing the evolution of the singlet scalar field within an inflationary paradigm. The singlet scalar field which explains the $\mu$-term tracks a time-varying minimum of the effective potential after inflation and slowly rolls down to its global minimum if there exist sufficiently large negative Hubble-induced corrections on the effective potential for the singlet field, which arise through supergravity. As a consequence, the whole Universe is confined within a single domain during and after inflation, which prevents the formation of domain walls. This will further constrain the history of the early Universe along with the Higgs-singlet coupling.


title: Quantum annealing speedup over simulated annealing on random Ising chains

abstract: We show clear evidence of a quadratic speedup of a quantum annealing (QA) Schroedinger dynamics over a Glauber master-equation simulated annealing (SA) for a random Ising model in one dimension, via an equal-footing exact deterministic dynamics of the Jordan-Wigner fermionized problems. This is remarkable, in view of the arguments of Katzgraber et al., PRX 4, 021008 (2014), since SA does not encounter any phase transition, while QA does. We also find a second remarkable result: that a "quantum-inspired" imaginary-time Schroedinger QA provides a further exponential speedup, i.e., an asymptotic residual error decreasing as a power-law $\tau^{-\mu}$ of the annealing time $\tau$.


title: Metamorphosis of SN 2014C: Delayed Interaction Between a Hydrogen Poor Core-collapse Supernova and a Nearby Circumstellar Shell

abstract: We present optical observations of supernova SN 2014C, which underwent an unprecedented slow metamorphosis from H-poor type Ib to H-rich type IIn over the course of one year. The observed spectroscopic evolution is consistent with the supernova having exploded in a cavity before encountering a massive shell of the progenitor star's stripped hydrogen envelope. Possible origins for the circumstellar shell include a brief Wolf-Rayet fast wind phase that overtook a slower red supergiant wind, eruptive ejection, or confinement of circumstellar material by external influences of neighboring stars. An extended high velocity Halpha absorption feature seen in near-maximum light spectra implies that the progenitor star was not completely stripped of hydrogen at the time of core collapse. Archival pre-explosion Subaru Telescope Suprime-Cam and Hubble Space Telescope Wide Field Planetary Camera 2 images of the region obtained in 2009 show a coincident source that is most likely a compact massive star cluster in NGC 7331 that hosted the progenitor system. By comparing the emission properties of the source with stellar population models that incorporate interacting binary stars we estimate the age of the host cluster to be 30 - 300 Myr, and favor ages closer to 30 Myr in light of relatively strong Halpha emission. SN 2014C is the best-observed member of a class of core-collapse supernovae that fill the gap between events that interact strongly with dense, nearby environments immediately after explosion and those that never show signs of interaction. Better understanding of the frequency and nature of this intermediate population can contribute valuable information about the poorly understood final stages of stellar evolution.


title: Towards a Theory of the QCD String

abstract: We construct a new model of four-dimensional relativistic strings with integrable dynamics on the worldsheet. In addition to translational modes this model contains a single massless pseudoscalar worldsheet field - the worldsheet axion. The axion couples to a topological density which counts the self-intersection number of a string. The corresponding coupling is fixed by integrability to $Q=\sqrt{7\over 16\pi}\approx 0.37$. We argue that this model is a member of a larger family of relativistic non-critical integrable string models. This family includes and extends conventional non-critical strings described by the linear dilaton CFT. Intriguingly, recent lattice data in $SU(3)$ and $SU(5)$ gluodynamics reveals the presence of a massive pseudoscalar axion on the worldsheet of confining flux tubes. The value of the corresponding coupling, as determined from the lattice data, is equal to $Q_L\approx0.38\pm0.04$.


title: Effective Window Function for Lagrangian Halos

abstract: The window function for protohalos in Lagrangian space is often assumed to be a tophat in real space. We measure this profile directly and find that it is more extended than a tophat but less extended than a Gaussian; its shape is well-described by rounding the edges of the tophat by convolution with a Gaussian that has a scale length about 5 times smaller. This effective window $W_{\rm eff}$ is particularly simple in Fourier space, and has an analytic form in real space. Together with the excursion set bias parameters, $W_{\rm eff}$ describes the scale-dependence of the Lagrangian halo-matter cross correlation up to $kR_{\rm Lag} \sim 10 $, where $R_{\rm Lag}$ is the Lagrangian size of the protohalo. Moreover, with this $W_{\rm eff}$, all the spectral moments of the power spectrum are finite, allowing a straightforward estimate of the excursion set peak mass function. This estimate requires a prescription of the critical overdensity enclosed within a protohalo if it is to collapse, which we calibrate from simulations. We find that the resulting estimate of halo abundances is only accurate to about 20%, and we discuss why: A tophat in `infall time' towards the protohalo center need not correspond to a tophat in the initial spatial distribution, so models in which infall rather than smoothed overdensity is the relevant variable may be more accurate.


title: Radiatively Induced Fermi Scale in Grand Unified Theories

abstract: We consider Grand Unified Theories in which the hierarchy between the unification and the Fermi scale emerges radiatively. Within the Pati-Salam framework, we show that it is possible to construct a viable model where the Higgs is an elementary pseudo-Goldstone boson, and the correct hierarchy is generated.


title: The evolution of the equivalent width of the Ha emission line and specific star-formation rate in star-forming galaxies at 1<z<5

abstract: We present the results of a study which uses spectral energy distribution (SED) fitting to investigate the evolution of the equivalent width (EW) of the Halpha emission line in star-forming galaxies over the redshift interval 1<z<5. After first demonstrating the ability of our SED-fitting technique to recover EW(Ha) using a sample of galaxies at z~1.3 with EW(Ha) measurements from 3D-HST grism spectroscopy, we proceed to apply our technique to samples of spectroscopically confirmed and photometric-redshift selected star-forming galaxies at z>=1 in the CANDELS UDS and GOODS-S fields. Confining our analysis to a constant stellar mass range (9.5<log(M/Msun)<10.5), we find that the median EW(Ha) evolves only modestly with redshift, reaching a rest-frame value of EW(Ha)=301+/-30 Angs by redshift z~4.5. Furthermore, using estimates of star-formation rate (SFR) based on both UV luminosity and Ha line flux, we use our galaxy samples to compare the evolution of EW(Ha) and specific star-formation rate (sSFR). Our results indicate that over the redshift range 1<z<5, the evolution displayed by EW(Ha) and sSFR is consistent, and can be adequately parameterized as: propto (1+z)^(1.0+/-0.2). As a consequence, over this redshift range we find that the sSFR and rest-frame EW(Ha) of star-forming galaxies with stellar masses M~10^(10) Msun are related by: EW(Ha)/Ang=(63+/-7)sSFR/Gyr^(-1). Given the current uncertainties in measuring the SFRs of high-redshift galaxies, we conclude that EW(Ha) provides a useful independent tracer of sSFR for star-forming galaxies out to redshifts of z=5.


title: Main-belt comets: sublimation-driven activity in the asteroid belt

abstract: Our knowledge of main-belt comets (MBCs), which exhibit comet-like activity likely due to the sublimation of volatile ices, yet orbit in the main asteroid belt, has increased greatly since the discovery of the first known MBC, 133P/Elst-Pizarro, in 1996, and their recognition as a new class of solar system objects after the discovery of two more MBCs in 2005. I review work that has been done over the last 10 years to improve our understanding of these enigmatic objects, including the development of systematic discovery methods and diagnostics for distinguishing MBCs from disrupted asteroids (which exhibit comet-like activity due to physical disruptions such as impacts or rotational destabilization). I also discuss efforts to understand the dynamical and thermal properties of these objects.


title: Charged black strings in a five-dimensional Kasner universe

abstract: We construct time-dependent charged black string solutions in five-dimensional Einstein-Maxwell theory. In the far region, the spacetime approaches a five-dimensional Kasner universe with a expanding three-dimensional space and a shrinking extra dimension. Near the event horizon, the spacetime is approximately static and has a smooth event horizon. We also study the motion of test particles around the black string and show the existence of quasi-circular orbits. Finally, we briefly discuss the stability of this spacetime.


title: A nebular analysis of the central Orion Nebula with MUSE

abstract: A nebular analysis of the central Orion Nebula and its main structures is presented. We exploit MUSE integral field observations in the wavelength range 4595-9366 \r{A} to produce the first O, S and N ionic and total abundance maps of a region spanning 6' x 5' with a spatial resolution of 0.2". We use the S$_{23}$ ( = ([SII]$\lambda$6717,31+[SIII]$\lambda$9068)/H$\beta$) parameter, together with [OII]/[OIII] as an indicator of the degree of ionisation, to distinguish between the various small-scale structures. The only Orion Bullet covered by MUSE is HH 201, which shows a double component in the [FeII]$\lambda$8617 line throughout indicating an expansion, and we discuss a scenario in which this object is undergoing a disruptive event. We separate the proplyds located south of the Bright Bar into four categories depending on their S$_{23}$ values, propose the utility of the S$_{23}$ parameter as an indicator of the shock-contribution to the excitation of line-emitting atoms, and show that the MUSE data is able to identify the proplyds associated with disks and microjets. We compute the second order structure function for the H$\alpha$, [OIII]$\lambda$5007, [SII]$\lambda$6731 and [OI]$\lambda$6300 emission lines to analyse the turbulent velocity field of the region covered with MUSE. We find that the spectral and spatial resolution of MUSE is not able to faithfully reproduce the structure functions of previous works.


title: NuSTAR and XMM-Newton Observations of the Hard X-Ray Spectrum of Centaurus A

abstract: We present simultaneous XMM-Newton and NuSTAR observations spanning 3-78 keV of the nearest radio galaxy, Centaurus A (Cen A). The accretion geometry around the central engine in Cen A is still debated, and we investigate possible configurations using detailed X-ray spectral modeling. NuSTAR imaged the central region of Cen A with sub-arcminute resolution at X-ray energies above 10 keV for the first time, but finds no evidence for an extended source or other off-nuclear point-sources. The XMM-Newton and NuSTAR spectra agree well and can be described with an absorbed power-law with a photon index {\Gamma} = 1.815 +/- 0.005 and a fluorescent Fe K{\alpha} line in good agreement with literature values. The spectrum does not require a high-energy exponential rollover, with a constraint of E_fold > 1 MeV. A thermal Comptonization continuum describes the data well, with parameters that agree with values measured by INTEGRAL, in particular an electron temperature kTe between ~100-300 keV, seed photon input temperatures between 5-50 eV. We do not find evidence for reflection or a broad iron line and put stringent upper limits of R < 0.01 on the reflection fraction and accretion disk illumination. We use archival Chandra data to estimate the contribution from diffuse emission, extra-nuclear point-sources, and the outer X-ray jet to the observed NuSTAR and XMM-Newton X-ray spectra and find the contribution to be negligible. We discuss different scenarios for the physical origin of the observed hard X-ray spectrum, and conclude that the inner disk is replaced by an advection-dominated accretion flow or that the X-rays are dominated by synchrotron self-Compton emission from the inner regions of the radio jet or a combination thereof.


title: Graphs that are simultaneously efficient open domination and efficient closed domination graphs

abstract: A graph is an efficient open (resp.\ closed) domination graph if there exists a subset of vertices whose open (resp.\ closed) neighborhoods partition its vertex set. Graphs that are efficient open as well as efficient closed (shortly EOCD graphs) are investigated. The structure of EOCD graphs with respect to their efficient open and efficient closed dominating sets is explained. It is shown that the decision problem regarding whether a graph is an EOCD graph is an NP-complete problem. A recursive description that constructs all EOCD trees is given and EOCD graphs are characterized among the Sierpi\'nski graphs.


title: Active Asteroids: Main-Belt Comets and Disrupted Asteroids

abstract: The study of active asteroids has attracted a great deal of interest in recent years since the recognition of main-belt comets (which orbit in the main asteroid belt, but exhibit comet-like activity due to the sublimation of volatile ices) as a new class of comets in 2006, and the discovery of the first disrupted asteroids (which, unlike MBCs, exhibit comet-like activity due to a physical disruption such as an impact or rotational destabilization, not sublimation) in 2010. In this paper, I will briefly discuss key areas of interest in the study of active asteroids.


title: Sequential Electronic and Structural Transitions in VO2 Observed Using X-ray Absorption Spectromicroscopy

abstract: The popular dual electronic (Mott) and structural (Peierls) transitions in VO2 are explored using x-ray absorption spectromicroscopy with high spatial and spectral resolutions. It is found that during both heating and cooling, the electronic transition always precedes the structural Peierls transition. Between the two transitions, there are intermediate states that are spectrally isolated here.


title: LoCuSS: Testing hydrostatic equilibrium in galaxy clusters

abstract: We test the assumption of hydrostatic equilibrium in an X-ray luminosity selected sample of 50 galaxy clusters at $0.15<z<0.3$ from the Local Cluster Substructure Survey (LoCuSS). Our weak-lensing measurements of $M_{500}$ control systematic biases to sub-4 per cent, and our hydrostatic measurements of the same achieve excellent agreement between XMM-Newton and Chandra. The mean ratio of X-ray to lensing mass for these 50 clusters is $\beta_{\rm X}=0.95\pm0.05$, and for the 44 clusters also detected by Planck, the mean ratio of Planck mass estimate to LoCuSS lensing mass is $\beta_{\rm P}=0.95\pm0.04$. Based on a careful like-for-like analysis, we find that LoCuSS, the Canadian Cluster Comparison Project (CCCP), and Weighing the Giants (WtG) agree on $\beta_{\rm P}\simeq0.9-0.95$ at $0.15<z<0.3$. This small level of hydrostatic bias disagrees at $\sim5\sigma$ with the level required to reconcile Planck cosmology results from the cosmic microwave background and galaxy cluster counts.


title: Accelerator physics and technology research toward future multi-MW proton accelerators

abstract: Recent P5 report indicated the accelerator-based neutrino and rare decay physics research as a centrepiece of the US domestic HEP program. Operation, upgrade and development of the accelerators for the near-term and longer-term particle physics program at the Intensity Frontier face formidable challenges. Here we discuss accelerator physics and technology research toward future multi-MW proton accelerators.


title: Characterization of electronic structure of periodically strained graphene

abstract: We induced periodic biaxial tensile strain in polycrystalline graphene by wrapping it over a substrate with repeating pillar-like structures with a periodicity of 600 nm. Using Raman spectroscopy, we determined to have introduced biaxial strains in graphene in the range of 0.4% to 0.7%. Its band structure was characterized using photoemission from valance bands, shifts in the secondary electron emission, and x-ray absorption from the carbon 1s levels to the unoccupied graphene conduction bands. It was observed that relative to unstrained graphene, strained graphene had a higher work function and higher density of states in the valence and conduction bands. We measured the conductivity of the strained and unstrained graphene in response to a gate voltage and correlated the changes in their behavior to the changes in the electronic structure. From these sets of data, we propose a simple band diagram representing graphene with periodic biaxial strain.


title: A Multi-stage Stochastic Programming Approach for Network Capacity Expansion with Multiple Sources of Capacity

abstract: In networks, there are often more than one source of capacity. The capacities can be permanently or temporarily owned by the decision maker. Depending on the nature of sources, we identify the permanent capacity, spot market capacity and contract capacity. We use a scenario tree to model the uncertainty, and build a multi-stage stochastic integer program that can incorporate multiple sources and multiple types of capacities in a general network. We propose two solution methodologies for the problem. Firstly, we design an asymptotically convergent approximation algorithm. Secondly, we design a cutting plane algorithm based on Benders decomposition to find tight bounds for the problem. The numerical experiments show superb performance of the proposed algorithms compared with commercial software.


title: UV Spectroscopy of Star-Grazing Comets within the 49 Ceti Debris Disk

abstract: We present analysis of time-variable, shifted absorption features in far-UV spectra of the unusual 49 Ceti debris disk. This nearly edge-on disk is one of the brightest known, and is one of the very few containing detectable amounts of circumstellar gas as well as dust. In our two visits of Hubble Space Telescope STIS spectra, variable absorption features are seen on the wings of lines arising from C II and C IV, but not for any of the other circumstellar absorption lines. Similar variable features have long been seen in spectra of the well-studied $\beta$ Pictoris debris disk and attributed to the transits of star-grazing comets. We calculate the velocity ranges and apparent column densities of the 49 Cet variable gas, which appears to be moving at velocities of tens to hundreds of km s$^{-1}$ relative to the central star. The velocities of the gas in the redshifted variable event in Visit 2 show that the maximum distances of the infalling gas at the time of transit are about 0.05 to 0.2 AU from the central star. A preliminary attempt at a composition analysis of the redshifted event suggests that the C/O ratio in the infalling gas may be super-solar, as it is in the bulk of the stable disk gas.


title: An Overview of the New Test Stand for H- Ion Sources at FNAL

abstract: A new test stand at FNAL is being constructed to carry out experiments to develop and upgrade the present magnetron-type sources of H- ions of up to 80 mA at 35 keV in the context of the Proton Improvement Plan (PIP). The aim of this plan is to provide high-power proton beams for the experiments at FNAL. The technical details of the construction and layout of this test stand are presented, along with a prospective set of diagnostics to monitor the sources.


title: Characteristics of undulatory locomotion in granular media

abstract: Undulatory locomotion is ubiquitous in nature and observed in different media, from the swimming of flagellated microorganisms in biological fluids, to the slithering of snakes on land, or the locomotion of sandfish lizards in sand. Despite the similarity in the undulating pattern, the swimming characteristics depend on the rheological properties of different media. Analysis of locomotion in granular materials is relatively less developed compared with fluids partially due to a lack of validated force models but recently a resistive force theory in granular media has been proposed and shown useful in studying the locomotion of a sand-swimming lizard. Here we employ the proposed model to investigate the swimming characteristics of a slender filament, of both finite and infinite length, undulating in a granular medium and compare the results with swimming in viscous fluids. In particular, we characterize the effects of drifting and pitching in terms of propulsion speed and efficiency for a finite sinusoidal swimmer. We also find that, similar to Lighthill's results using resistive force theory in viscous fluids, the sawtooth swimmer is the optimal waveform for propulsion speed at a given power consumption in granular media. The results complement our understanding of undulatory locomotion and provide insights into the effective design of locomotive systems in granular media.


title: X-ray long-term variations in the low-luminosity AGN NGC835 and its circumnuclear emission

abstract: Obscured active galactic nuclei (AGNs) are thought to be very common in the Universe. Observations and surveys have shown that the number of sources increases for near galaxies and at the low-luminosity regime (the so-called LLAGNs). Furthermore, many AGNs show changes in their obscuration properties at X-rays that may suggest a configuration of clouds very close to the accretion disk. However, these variations could also be due to changes in the intrinsic continuum of the source. It is therefore important to study nearby AGN to better understand the locus and distribution of clouds in the neighbourhood of the nucleus. We aim to study the nuclear obscuration of LLAGN NGC835 and its extended emission using mid-infrared observations. We present mid-infrared 11.5 microns imaging of the LLAGN galaxy NGC835 obtained with the instrument CanariCam in the Gran Telescopio CANARIAS (GTC), archival Spitzer/IRS spectroscopy, and archival Chandra data observed in 2000, 2008, and 2013. The GTC/CanariCam 11.5 microns image reveals faint extended emission out to ~6 arcsec. We obtained a nuclear flux of F(11.5 microns) ~18 mJy, whereas the extended emission accounts for 90% of the total flux within the 6 arcsec. This means that the low angular resolution (~4 arcsec) IRS spectrum is dominated by this extended emission and not by the AGN, clearly seen in the Spitzer/IRS spectrum. Although the extended soft X-ray emission shows some resemblance with that of the mid-infrared, the knots seen at X-rays are mostly located in the inner side of this mid-infrared emission. The nuclear X-ray spectrum of the source has undergone a spectral change between 2000/2008 and 2013. We argue that this is most probably due to changes in the hydrogen column density from ~ 8x10E+23 cm-2 to ~ 3x10E+23 cm-2. NGC835 therefore is one of the few LLAGN, together with NGC1052, in which changes in the absorber can be claimed.


title: Texture and Wettability of Metallic Lotus Leaves

abstract: Superhydrophobic surfaces with the self-cleaning behavior of lotus leaves are sought for drag reduction and phase change heat transfer applications. These superrepellent surfaces have traditionally been fabricated by random or deterministic texturing of a hydrophobic material. Recently, superrepellent surfaces have also been made from hydrophilic materials, by deterministic texturing using photolithography, without low-surface energy coating. Here, we show that hydrophilic materials can also be made superrepellent to water by chemical texturing, a stochastic rather than deterministic process. These metallic surfaces are the first analog of lotus leaves, in terms of wettability, texture and repellency. A mechanistic model is also proposed to describe the influence of multiple tiers of roughness on wettability and repellency. This demonstrated ability to make hydrophilic materials superrepellent without deterministic structuring or additional coatings opens the way to large scale and robust manufacturing of superrepellent surfaces.


title: Quasi-Hamiltonian description of classical spin

abstract: A family of Lagrange functions is considered, each producing the classical relativistic free spinning particle equation of motion of the third order. On this grounds a generalized Hamilton-Ostrohrads'kyj description of the free relativistic spherical top is proposed, which comply with the Pirani supplementary conditions.


title: Active compensation of aperture discontinuities for WFIRST-AFTA: analytical and numerical comparison of propagation methods and preliminary results with a WFIRST-AFTA-like pupil

abstract: The new frontier in the quest for the highest contrast levels in the focal plane of a coronagraph is now the correction of the large diffractive artifacts effects introduced at the science camera by apertures of increasing complexity. The coronagraph for the WFIRST/AFTA mission will be the first of such instruments in space with a two Deformable Mirrors wavefront control system. Regardless of the control algorithm for these multi Deformable Mirrors, they will have to rely on quick and accurate simulation of the propagation effects introduced by the out-of-pupil surface. In the first part of this paper, we present the analytical description of the different approximations to simulate these propagation effects. In Annex A, we prove analytically that, in the special case of surfaces inducing a converging beam, the Fresnel method yields high fidelity for simulations of these effects. We provide numerical simulations showing this effect. In the second part, we use these tools in the framework of the Active Compensation of Aperture Discontinuities technique (ACAD) applied to pupil geometries similar to WFIRST-AFTA. We present these simulations in the context of the optical layout of the High-contrast imager for Complex Aperture Telescopes, which will test ACAD on a optical bench. The results of this analysis show that using the ACAD method, an apodized pupil lyot coronagraph and the performance of our current deformable mirrors, we are able to obtain, in numerically simulations, a dark hole with an AFTA-like pupil. Our numerical simulation shows that we can obtain contrast better than $2.10^{-9}$ in monochromatic light and better than 3.e-8 with 10% bandwidth between 5 and 14 lambda/D.


title: On the Matsumoto-Yor property in free probability

abstract: We study the Matsumoto-Yor property in free probability. We prove that the limiting empirical eigenvalue distribution of the GIG matrices and the Marchenko-Pastur distribution have the free Matsumoto-Yor property. Finally we characterize these distributions by a regression properties in free probability.


title: Fermilab Muon Campus g-2 Cryogenic Distribution Remote Control System

abstract: The Muon Campus (MC) is able to measure Muon g-2 with high precision and comparing its value to the theoretical prediction. The MC has four 300 KW screw compressors and four liquid helium refrigerators. The centerpiece of the Muon g-2 experiment at Fermilab is a large, 50-foot-diameter superconducting muon storage ring. This one-of-a-kind ring, made of steel, aluminum and superconducting wire, was built for the previous g-2 experiment at Brookhaven. Due to each subsystem has to be far away from each other and be placed in the distant location, therefore, Siemens Process Control System PCS7-400, Automation Direct DL205 & DL05 PLC, Synoptic and Fermilab ACNET HMI are the ideal choices as the MC g-2 cryogenic distribution real-time and on-Line remote control system. This paper presents a method which has been successfully used by many Fermilab distribution cryogenic real-time and On-Line remote control systems.


title: On the use of Perfectly Matched Layers at corners for scattering problems with sign-changing coefficients

abstract: We investigate in a $2$D setting the scattering of time-harmonic electromagnetic waves by a plasmonic device, represented as a non dissipative bounded and penetrable obstacle with a negative permittivity. Using the $\textrm{T}$-coercivity approach, we first prove that the problem is well-posed in the classical framework $H^1_{\text{loc}} $ if the negative permittivity does not lie in some critical interval whose definition depends on the shape of the device. When the latter has corners, for values inside the critical interval, unusual strong singularities for the electromagnetic field can appear. In that case, well-posedness is obtained by imposing a radiation condition at the corners to select the outgoing black-hole plasmonic wave, that is the one which carries energy towards the corners. A simple and systematic criterion is given to define what is the outgoing solution. Finally, we propose an original numerical method based on the use of Perfectly Matched Layers at the corners. We emphasize that it is necessary to design an $\textit{ad hoc}$ technique because the field is too singular to be captured with standard finite element methods.


title: Interface control of electronic transport across the magnetic phase transition in SrRuO_3/SrTiO_3 heterointerface

abstract: The emerging material class of complex-oxides, where manipulation of physical properties lead to new functionalities at their heterointerfaces, is expected to open new frontiers in Spintronics. For example, SrRuO_3 is a promising material where external stimuli like strain, temperature and structural distortions control the stability of electronic and magnetic states, across its magnetic phase transition, useful for Spintronics. Despite this, not much has been studied to understand such correlations in SrRuO_3. Here we explore the influence of electron-lattice correlation to electron-transport, at interfaces between SrRuO_3 and Nb:SrTiO_3 across its ferromagnetic transition, using a nanoscale transport probe and first-principles calculations. We find that the geometrical reconstructions at the interface and hence modifications in electronic structures dominate the transmission across its ferromagnetic transition, eventually flipping the charge-transport length-scale in SrRuO_3. This approach can be easily extended to other devices where competing ground states can lead to different functional properties across their heterointerfaces.


title: Crystal Ball: On the Future High Energy Colliders

abstract: High energy particle colliders have been in the forefront of particle physics for more than three decades. At present the near term US, European and international strategies of the particle physics community are centered on full exploitation of the physics potential of the Large Hadron Collider (LHC) through its high-luminosity upgrade (HL-LHC). A number of next generation collider facilities have been proposed and are currently under consideration for the medium- and far-future of the accelerator-based high energy physics. In this paper we offer a uniform approach to evaluation of various accelerators based on the feasibility of their energy reach, performance reach and cost range. We briefly review such post-LHC options as linear e+e- colliders in Japan (ILC) or at CERN (CLIC), muon collider, and circular lepton or hadron colliders in China (CepC/SppC) and Europe (FCC). We conclude with a look into ultimate energy reach accelerators based on plasmas and crystals, and some perspectives for the far future of accelerator-based particle physics.


title: A Multiresolution Ensemble Kalman Filter using Wavelet Decomposition

abstract: We present a method of using classical wavelet based multiresolution analysis to separate scales in model and observations during data assimilation with the ensemble Kalman filter. In many applications, the underlying physics of a phenomena involve the interaction of features at multiple scales. Blending of observational and model error across scales can result in large forecast inaccuracies since large errors at one scale are interpreted as inexact data at all scales. Our method uses a transformation of the observation operator in order to separate the information from different scales of the observations. This naturally induces a transformation of the observation covariance and we put forward several algorithms to efficiently compute the transformed covariance. Another advantage of our multiresolution ensemble Kalman filter is that scales can be weighted independently to adjust each scale's effect on the forecast. To demonstrate feasibility we present applications to a one dimensional Kuramoto-Sivashinsky (K-S) model with scale dependent observation noise and an application involving the forecasting of solar photospheric flux. The latter example demonstrates the multiresolution ensemble Kalman filter's ability to account for scale dependent model error. Modeling of photospheric magnetic flux transport is accomplished by the Air Force Data Assimilative Photospheric Transport (ADAPT) model.


title: A new approach to barrier-top fission dynamics

abstract: We proposed a calculational framework for describing induced fission that avoids the Bohr-Wheeler assumption of well-defined fission channels. The building blocks of our approach are configurations that form a discrete, orthogonal basis and can be characterized by both energy and shape. The dynamics is to be determined by interaction matrix elements between the states rather than by a Hill-Wheeler construction of a collective coordinate. Within our approach, several simple limits can be seen: diffusion; quantized conductance; and ordinary decay through channels. The specific proposal for the discrete basis is to use the $K^\pi$ quantum numbers of the axially symmetric Hartree-Fock approximation to generate the configurations. Fission paths would be determined by hopping from configuration to configuration via the residual interaction. We show as an example the configurations needed to describe a fictitious fission decay $^{32}{\rm S} \rightarrow ^{16}{\rm O} + ^{16}{\rm O}$. We also examine the geometry of the path for fission of $^{236}$U, measuring distances by the number of jumps needed to go to a new $K^\pi$ partition.


title: Separations in query complexity using cheat sheets

abstract: We show a power 2.5 separation between bounded-error randomized and quantum query complexity for a total Boolean function, refuting the widely believed conjecture that the best such separation could only be quadratic (from Grover's algorithm). We also present a total function with a power 4 separation between quantum query complexity and approximate polynomial degree, showing severe limitations on the power of the polynomial method. Finally, we exhibit a total function with a quadratic gap between quantum query complexity and certificate complexity, which is optimal (up to log factors). These separations are shown using a new, general technique that we call the cheat sheet technique. The technique is based on a generic transformation that converts any (possibly partial) function into a new total function with desirable properties for showing separations. The framework also allows many known separations, including some recent breakthrough results of Ambainis et al., to be shown in a unified manner.


title: The mathematics of superoscillations

abstract: In the past 50 years, quantum physicists have discovered, and experimentally demonstrated, a phenomenon which they termed superoscillations. Aharonov and his collaborators showed that superoscillations naturally arise when dealing with weak values, a notion that provides a fundamentally different way to regard measurements in quantum physics. From a mathematical point of view, superoscillating functions are a superposition of small Fourier components with a bounded Fourier spectrum, which result, when appropriately summed, in a shift that can be arbitrarily large, and well outside the spectrum. Purpose of this work is twofold: on one hand we provide a self-contained survey of the existing literature, in order to offer a systematic mathematical approach to superoscillations; on the other hand, we obtain some new and unexpected results, by showing that superoscillating sequences can be seen of as solutions to a large class of convolution equations and can therefore be treated within the theory of Analytically Uniform spaces. In particular, we will also discuss the persistence of the superoscillatory behavior when superoscillating sequences are taken as initial values of the Schr\"odinger equation and other equations.


title: Noninertial effects on the quantum dynamics of scalar bosons

abstract: The noninertial effect of rotating frames on the quantum dynamics of scalar bosons embedded in the background of a cosmic string is considered. In this work, scalar bosons are described by the Duffin--Kemmer--Petiau (DKP) formalism. Considering the DKP oscillator in this background the combined effects of a rotating frames and cosmic string on the equation of motion, energy spectrum and DKP spinor are analyzed and discussed in detail. Additionally, the effect of rotating frames on the scalar bosons localization is studied.


title: Large-scale CO (J=4-3) Mapping toward the Orion-A Giant Molecular Cloud

abstract: We have mapped the Orion-A Giant Molecular Cloud in the CO (J=4-3) line with the Tsukuba 30-cm submillimeter telescope.The map covered a 7.125 deg^2 area with a 9' resolution, including main components of the cloud such as Orion Nebula, OMC-2/3, and L1641-N. The most intense emission was detected toward the Orion KL region. The integrated intensity ratio between CO (J=4-3) and CO (J=1-0) was derived using data from the Columbia-Univ. de Chile CO survey, which was carried out with a comparable angular resolution. The ratio was r_{4-3/1-0} ~ 0.2 in the southern region of the cloud and 0.4-0.8 at star forming regions. We found a trend that the ratio shows higher value at edges of the cloud. In particular the ratio at the north-eastern edge of the cloud at (l, b) = (208.375 deg, -19.0 deg) shows the specific highest value of 1.1. The physical condition of the molecular gas in the cloud was estimated by non-LTE calculation. The result indicates that the kinetic temperature has a gradient from north (Tkin=80 K) to south (20 K). The estimation shows that the gas associated with the edge of the cloud is warm (Tkin~60 K), dense (n_{H_2}~10^4 cm^{-3}), and optically thin, which may be explained by heating and sweeping of interstellar materials from OB clusters.


title: Prediction of Solar Flare Size and Time-to-Flare Using Support Vector Machine Regression

abstract: We study the prediction of solar flare size and time-to-flare using 38 features describing magnetic complexity of the photospheric magnetic field. This work uses support vector regression to formulate a mapping from the 38-dimensional feature space to a continuous-valued label vector representing flare size or time-to-flare. When we consider flaring regions only, we find an average error in estimating flare size of approximately half a \emph{geostationary operational environmental satellite} (\emph{GOES}) class. When we additionally consider non-flaring regions, we find an increased average error of approximately 3/4 a \emph{GOES} class. We also consider thresholding the regressed flare size for the experiment containing both flaring and non-flaring regions and find a true positive rate of 0.69 and a true negative rate of 0.86 for flare prediction. The results for both of these size regression experiments are consistent across a wide range of predictive time windows, indicating that the magnetic complexity features may be persistent in appearance long before flare activity. This is supported by our larger error rates of some 40 hr in the time-to-flare regression problem. The 38 magnetic complexity features considered here appear to have discriminative potential for flare size, but their persistence in time makes them less discriminative for the time-to-flare problem.


title: Stop Wasting My Gradients: Practical SVRG

abstract: We present and analyze several strategies for improving the performance of stochastic variance-reduced gradient (SVRG) methods. We first show that the convergence rate of these methods can be preserved under a decreasing sequence of errors in the control variate, and use this to derive variants of SVRG that use growing-batch strategies to reduce the number of gradient calculations required in the early iterations. We further (i) show how to exploit support vectors to reduce the number of gradient computations in the later iterations, (ii) prove that the commonly-used regularized SVRG iteration is justified and improves the convergence rate, (iii) consider alternate mini-batch selection strategies, and (iv) consider the generalization error of the method.


title: Radiative transfer with long-range interactions: regularity and asymptotics

abstract: This work is devoted to radiative transfer equations with long-range interactions. Such equations arise in the modeling of high frequency wave propagation in random media with long-range dependence. In the regime we consider, the singular collision operator modeling the interaction between the wave and the medium is conservative, and as a consequence wavenumbers take values on the unit sphere. Our goals are to investigate the regularizing effects of grazing collisions, the diffusion limit, and the peaked forward limit. As in the case where wavenumbers take values in R^{d+1}, we show that the transport operator is hypoelliptic, so that the solutions are infinitely differerentiable in all variables. Using probabilistic techniques, we show as well that the diffusion limit can be carried on as in the case of a regular collision operator, and as a consequence that the diffusion coefficient is non-zero and finite. We finally consider the regime where grazing collisions are dominant.


title: A closer view of the IGR J11014-6103 outflows

abstract: The complex X-ray system IGR J11014-6103 (a.k.a. the Lighthouse nebula) is composed of a bow-shock pulsar wind nebula (PWN) as well as large-scale jet-like features, all launched by IGR J11014-6103 which is moving supersonically in the interstellar medium. Previous observations suggested that the jet features stem from a ballistic jet of relativistic particles. In order to confirm the nature of the jet and the marginally detected counter-jet, we obtained a new deep 250 ks Chandra observation of the Lighthouse nebula. We performed detailed spatial and spectral analysis of all X-ray components of the system. The X-ray PWN is now better resolved and shows a clear bi-modal morphology. The overall helical pattern of the main jet is confirmed. However, there are large deviations from a simple helical model at small and large scales. Significant extended emission is now detected, encompassing the main jet all along its length. The brightness dip of the main jet at ~50" distance from the pulsar is confirmed, the extended emission however prevents conclusions about the coherence at this position of the jet. The counter-jet is now detected at high statistical significance. In addition, we found two small-scale "arcs" departing from the pulsar towards the jets. We also looked for possible bow-shock emission in front of the pulsar, with a short VLT/FORS2 H-alpha observation. No clear emission is found, most likely because of contamination from a diffuse surrounding nebulosity. The results of our X-ray analysis show that current expectations from a ballistic nature of the jets can explain satisfactorily some of the observational evidences but cannot fully reproduce the observations. The alternative scenario (diffusion of particles along preexisting magnetic field lines in the surrounding medium) however also continues to suffer from conflicts with the observations.


title: Regularity theory and extension problem for fractional nonlocal parabolic equations and the master equation

abstract: We develop the regularity theory for solutions to space-time nonlocal equations driven by fractional powers of the heat operator $$(\partial_t-\Delta)^su(t,x)=f(t,x),\quad\hbox{for}~0<s<1.$$ This nonlocal equation of order $s$ in time and $2s$ in space arises in Nonlinear Elasticity, Semipermeable Membranes, Continuous Time Random Walks and Mathematical Biology. It plays for space-time nonlocal equations like the generalized master equation the same role as the fractional Laplacian for nonlocal in space equations. We obtain a pointwise integro-differential formula for $(\partial_t-\Delta)^su(t,x)$ and parabolic maximum principles. A novel extension problem to characterize this nonlocal equation with a local degenerate parabolic equation is proved. We show parabolic interior and boundary Harnack inequalities, and an Almgrem-type monotonicity formula. H\"older and Schauder estimates for the space-time Poisson problem are deduced using a new characterization of parabolic H\"older spaces. Our methods involve the \textit{parabolic language of semigroups} and the Cauchy Integral Theorem, which are original to define the fractional powers of $\partial_t-\Delta$. Though we mainly focus in the equation $(\partial_t-\Delta)^su=f$, applications of our ideas to variable coefficients, discrete Laplacians and Riemannian manifolds are stressed out.


title: SecureD: A Secure Dual Core Embedded Processor

abstract: Security of embedded computing systems is becoming of paramount concern as these devices become more ubiquitous, contain personal information and are increasingly used for financial transactions. Security attacks targeting embedded systems illegally gain access to the information in these devices or destroy information. The two most common types of attacks embedded systems encounter are code-injection and power analysis attacks. In the past, a number of countermeasures, both hardware- and software-based, were proposed individually against these two types of attacks. However, no single system exists to counter both of these two prominent attacks in a processor based embedded system. Therefore, this paper, for the first time, proposes a hardware/software based countermeasure against both code-injection attacks and power analysis based side-channel attacks in a dual core embedded system. The proposed processor, named SecureD, has an area overhead of just 3.80% and an average runtime increase of 20.0% when compared to a standard dual processing system. The overhead were measured using a set of industry standard application benchmarks, with two encryption and five other programs.


title: The pro-nilpotent group topology on a free group

abstract: In this paper, we work on the pro-nilpotent group topology of a free group. First we investigate the closure of the product of finitely many subgroups of a free group in the pro-nilpotent group topology. We present an algorithm for the calculation of the closure in the pro-nilpotent group topology of the product of finitely many finitely generated subgroups of a free group. Then we deduce that the nil-closure of a rational subset is computable. We also prove that the pseudovarieties V malcev G_nil, for every decidable pseudovariety of monoids V, and J * G_nil are decidable.


title: Self-Trapping of Diskoseismic Corrugation Modes in Neutron Star Spacetimes

abstract: We examine the effects of higher-order multipole contributions of rotating neutron star (NS) spacetimes on the propagation of corrugation (c-)modes within a thin accretion disk. We find that the Lense-Thirring precession frequency, which determines the propagation region of the low-frequency fundamental corrugation modes, can experience a turnover allowing for c-modes to become self-trapped for sufficiently high dimensionless spin $j$ and quadrupole rotational deformability $\alpha$. If such self-trapping c-modes can be detected, e.g. through phase-resolved spectroscopy of the iron line for a high-spin low-mass accreting neutron star, this could potentially constrain the spin-induced NS quadrupole and the NS equation of state.


title: Black holes in an expanding universe and supersymmetry

abstract: This paper analyzes the supersymmetric solutions to five and six-dimensional minimal (un)gauged supergravities for which the bilinear Killing vector constructed from the Killing spinor is null. We focus on the spacetimes which admit an additional ${\rm SO}(1,1)$ boost symmetry. Upon the toroidal dimensional reduction along the Killing vector corresponding to the boost, we show that the solution in the ungauged case describes a charged, nonextremal black hole in a Friedmann-Lema\^itre-Robertson-Walker (FLRW) universe with an expansion driven by a massless scalar field. For the gauged case, the solution corresponds to a charged, nonextremal black hole embedded conformally into a Kantowski-Sachs universe. It turns out that these dimensional reductions break supersymmetry since the bilinear Killing vector and the Killing vector corresponding to the boost fail to commute. This represents a new mechanism of supersymmetry breaking that has not been considered in the literature before.


title: The Low Energy Constants of $SU(2)$ Partially Quenched Chiral Perturbation Theory from $N_{f}=2+1$ Domain Wall QCD

abstract: We have performed fits of the pseudoscalar masses and decay constants, from a variety of RBC-UKQCD domain wall fermion ensembles, to $SU(2)$ partially quenched chiral perturbation theory at next-to leading order (NLO) and next-to-next-to leading order (NNLO). We report values for 9 NLO and 8 linearly independent combinations of NNLO partially quenched low energy constants, which we compare to other lattice and phenomenological determinations. We discuss the size of successive terms in the chiral expansion and use our large set of low energy constants to make predictions for mass splittings due to QCD isospin breaking effects and the S-wave $\pi \pi$ scattering lengths. We conclude that, for the range of pseudoscalar masses explored in this work, $115~\mathrm{MeV} \lesssim m_{\rm PS} \lesssim 430~\mathrm{MeV}$, the NNLO $SU(2)$ expansion is quite robust and can fit lattice data with percent-scale accuracy.


title: The NuSTAR view of reflection and absorption in NGC 7582

abstract: NGC 7582 is a well-studied X-ray bright Seyfert 2 with moderately heavy ($N_{\text{H}}\sim10^{23}-10^{24}$~cm$^{-2}$), highly variable absorption and strong reflection spectral features. The spectral shape changed around the year 2000, dropping in observed flux and becoming much more highly absorbed. Two scenarios have been put forth to explain this spectral change: 1) the central X-ray source partially ``shut off'' around this time, decreasing in intrinsic luminosity, with a delayed decrease in reflection features due to the light-crossing time of the Compton-thick material or 2) the source became more heavily obscured, with only a portion of the power law continuum leaking through. NuSTAR observed NGC~7582 twice in 2012, two weeks apart, in order to quantify the reflection using high-quality data above 10 keV. We find that the most plausible scenario is that NGC 7582 has recently become more heavily absorbed by a patchy torus with a covering fraction of $\sim\,80-90\%$ and an equatorial column density of $\sim 3 \times10^{24}$ cm$^{-2}$. We find the need for an additional highly variable full-covering absorber with $N_{\text{H}}= 4-6 \times10^{23}$ cm$^{-2}$ in the line of sight, possibly associated with a hidden broad line region.


title: Nucleon Polarisabilities at and Beyond Physical Pion Masses

abstract: We examine the results of Chiral Effective Field Theory ($\chi$EFT) for the scalar- and spin-dipole polarisabilities of the proton and neutron, both for the physical pion mass and as a function of $m_\pi$. This provides chiral extrapolations for lattice-QCD polarisability computations. We include both the leading and sub-leading effects of the nucleon's pion cloud, as well as the leading ones of the $\Delta(1232)$ resonance and its pion cloud. The analytic results are complete at N$^2$LO in the $\delta$-counting for pion masses close to the physical value, and at leading order for pion masses similar to the Delta-nucleon mass splitting. In order to quantify the truncation error of our predictions and fits as $68$\% degree-of-belief intervals, we use a Bayesian procedure recently adapted to EFT expansions. At the physical point, our predictions for the spin polarisabilities are, within respective errors, in good agreement with alternative extractions using experiments and dispersion-relation theory. At larger pion masses we find that the chiral expansion of all polarisabilities becomes intrinsically unreliable as $m_\pi$ approaches about $300\;$MeV---as has already been seen in other observables. $\chi$EFT also predicts a substantial isospin splitting above the physical point for both the electric and magnetic scalar polarisabilities; and we speculate on the impact this has on the stability of nucleons. Our results agree very well with emerging lattice computations in the realm where $\chi$EFT converges. Curiously, for the central values of some of our predictions, this agreement persists to much higher pion masses. We speculate on whether this might be more than a fortuitous coincidence.


title: Weighted Sum-Throughput Maximization for MIMO Broadcast Channel: Energy Harvesting Under System Imperfection

abstract: In this work, a MIMO broadcast channel under the energy harvesting (EH) constraint and the peak power constraint is investigated. The transmitter is equipped with a hybrid energy storage system consisting of a perfect super capacitor (SC) and an inefficient battery, where both elements have limited energy storage capacities. In addition, the effect of data processing circuit power consumption is also addressed. To be specific, two extreme cases are studied here, where the first assumes ideal/zero circuit power consumption and the second considers a positive constant circuit power consumption where the circuit is always operating at its highest power level. The performance of these two extreme cases hence serve as the upper bound and the lower bound of the system performance in practice, respectively. In this setting, the offline scheduling with ideal and maximum circuit power consumptions are investigated. The associated optimization problems are formulated and solved in terms of weighted throughput optimization. Further, we extend to a general circuit power consumption model. To complement this work, some intuitive online policies are presented for all cases. Interestingly, for the case with maximum circuit power consumption, a close-to-optimal online policy is presented and its performance is shown to be comparable to its offline counterpart in the numerical results.


title: Recovering hard-to-find object instances by sampling context-based object proposals

abstract: In this paper we focus on improving object detection performance in terms of recall. We propose a post-detection stage during which we explore the image with the objective of recovering missed detections. This exploration is performed by sampling object proposals in the image. We analyze four different strategies to perform this sampling, giving special attention to strategies that exploit spatial relations between objects. In addition, we propose a novel method to discover higher-order relations between groups of objects. Experiments on the challenging KITTI dataset show that our proposed relations-based proposal generation strategies can help improving recall at the cost of a relatively low amount of object proposals.


title: One Cyclic Codes over $\mathbb{F}_{p^k} + v\mathbb{F}_{p^k} + v^2\mathbb{F}_{p^k} + ... + v^r\mathbb{F}_{p^k}$

abstract: In this paper, we investigate cyclic code over the ring $\mathbb{F}_{p^k} + v\mathbb{F}_{p^k} + v^2\mathbb{F}_{p^k} + ... + v^r\mathbb{F}_{p^k}$, where $v^{r+1}=v$, $p$ a prime number, $r>1$ and $\gcd(r,p)=1$, we prove as generalisation of P. Sol\'e et al. in 2015 that these codes are principally generated, give generator polynomial and idempotent depending on idempotents over this ring as response to an open problem related by J. QIAN et al. in 2005. we also give a gray map and proprieties of the related dual code.


title: Statistically-Consistent k-mer Methods for Phylogenetic Tree Reconstruction

abstract: Frequencies of $k$-mers in sequences are sometimes used as a basis for inferring phylogenetic trees without first obtaining a multiple sequence alignment. We show that a standard approach of using the squared-Euclidean distance between $k$-mer vectors to approximate a tree metric can be statistically inconsistent. To remedy this, we derive model-based distance corrections for orthologous sequences without gaps, which lead to consistent tree inference. The identifiability of model parameters from $k$-mer frequencies is also studied. Finally, we report simulations showing the corrected distance out-performs many other $k$-mer methods, even when sequences are generated with an insertion and deletion process. These results have implications for multiple sequence alignment as well, since $k$-mer methods are usually the first step in constructing a guide tree for such algorithms.


title: False Discoveries Occur Early on the Lasso Path

abstract: In regression settings where explanatory variables have very low correlations and there are relatively few effects, each of large magnitude, we expect the Lasso to find the important variables with few errors, if any. This paper shows that in a regime of linear sparsity---meaning that the fraction of variables with a non-vanishing effect tends to a constant, however small---this cannot really be the case, even when the design variables are stochastically independent. We demonstrate that true features and null features are always interspersed on the Lasso path, and that this phenomenon occurs no matter how strong the effect sizes are. We derive a sharp asymptotic trade-off between false and true positive rates or, equivalently, between measures of type I and type II errors along the Lasso path. This trade-off states that if we ever want to achieve a type II error (false negative rate) under a critical value, then anywhere on the Lasso path the type I error (false positive rate) will need to exceed a given threshold so that we can never have both errors at a low level at the same time. Our analysis uses tools from approximate message passing (AMP) theory as well as novel elements to deal with a possibly adaptive selection of the Lasso regularizing parameter.


title: Local equilibrium in bird flocks

abstract: The correlated motion of flocks is an instance of global order emerging from local interactions. An essential difference with analogous ferromagnetic systems is that flocks are active: animals move relative to each other, dynamically rearranging their interaction network. The effect of this off-equilibrium element is well studied theoretically, but its impact on actual biological groups deserves more experimental attention. Here, we introduce a novel dynamical inference technique, based on the principle of maximum entropy, which accodomates network rearrangements and overcomes the problem of slow experimental sampling rates. We use this method to infer the strength and range of alignment forces from data of starling flocks. We find that local bird alignment happens on a much faster timescale than neighbour rearrangement. Accordingly, equilibrium inference, which assumes a fixed interaction network, gives results consistent with dynamical inference. We conclude that bird orientations are in a state of local quasi-equilibrium over the interaction length scale, providing firm ground for the applicability of statistical physics in certain active systems.

categories: q-bio.PE cond-mat.stat-mech physics.bio-ph q-bio.QM


title: Density and bond-orientational relaxations in supercooled water

abstract: Recent computational studies have reported evidence of a metastable liquid-liquid phase transition (LLPT) in molecular models of water under deeply supercooled conditions. A competing hypothesis suggests, however, that non-equilibrium artifacts associated with coarsening of the stable crystal phase have been mistaken for an LLPT in these models. Such artifacts are posited to arise due to a separation of time scales in which density fluctuations in the supercooled liquid relax orders of magnitude faster than those associated with bond-orientational order. Here, we use molecular simulation to investigate the relaxation of density and bond-orientational fluctuations in three molecular models of water (ST2, TIP5P and TIP4P/2005) in the vicinity of their reported LLPT. For each model, we find that density is the slowly relaxing variable under such conditions. We also observe similar behavior in the coarse-grained mW model of water. Our findings therefore challenge the key physical assumption underlying the competing hypothesis.e find that density relaxes significantly faster than bond-orientational order, as incorrectly predicted by this competing hypothesis.


title: An Action Language for Multi-Agent Domains: Foundations

abstract: In multi-agent domains (MADs), an agent's action may not just change the world and the agent's knowledge and beliefs about the world, but also may change other agents' knowledge and beliefs about the world and their knowledge and beliefs about other agents' knowledge and beliefs about the world. The goals of an agent in a multi-agent world may involve manipulating the knowledge and beliefs of other agents' and again, not just their knowledge/belief about the world, but also their knowledge about other agents' knowledge about the world. Our goal is to present an action language (mA+) that has the necessary features to address the above aspects in representing and RAC in MADs. mA+ allows the representation of and reasoning about different types of actions that an agent can perform in a domain where many other agents might be present -- such as world-altering actions, sensing actions, and announcement/communication actions. It also allows the specification of agents' dynamic awareness of action occurrences which has future implications on what agents' know about the world and other agents' knowledge about the world. mA+ considers three different types of awareness: full-, partial- awareness, and complete oblivion of an action occurrence and its effects. This keeps the language simple, yet powerful enough to address a large variety of knowledge manipulation scenarios in MADs. The semantics of mA+ relies on the notion of state, which is described by a pointed Kripke model and is used to encode the agent's knowledge and the real state of the world. It is defined by a transition function that maps pairs of actions and states into sets of states. We illustrate properties of the action theories, including properties that guarantee finiteness of the set of initial states and their practical implementability. Finally, we relate mA+ to other related formalisms that contribute to RAC in MADs.


title: Topology of two-row Springer fibers for the even orthogonal and symplectic group

abstract: We construct an explicit topological model (similar to the topological Springer fibers appearing in work of Khovanov and Russell) for every two-row Springer fiber associated with the even orthogonal group and prove that the respective topological model is homeomorphic to its corresponding Springer fiber. This confirms a conjecture by Ehrig and Stroppel concerning the topology of the equal-row Springer fiber for the even orthogonal group. Moreover, we show that every two-row Springer fiber for the symplectic group is homeomorphic (even isomorphic as an algebraic variety) to a connected component of a certain two-row Springer fiber for the even orthogonal group.


title: Quantum Walks on Generalized Quadrangles

abstract: We study the transition matrix of a quantum walk on strongly regular graphs. It is proposed by Emms, Hancock, Severini and Wilson in 2006, that the spectrum of $S^+(U^3)$, a matrix based on the amplitudes of walks in the quantum walk, distinguishes strongly regular graphs. We probabilistically compute the spectrum of the line intersection graphs of two non-isomorphic generalized quadrangles of order $(5^2,5)$ under this matrix and thus provide strongly regular counter-examples to the conjecture.


title: Two Polarization Entangled Sources from the Same Semiconductor Chip

abstract: Generating nonclassical states of photons such as polarization entangled states on a monolithic chip is a crucial step towards practical applications of optical quantum information processing such as quantum computing and quantum key distribution. Here we demonstrate two polarization entangled photon sources in a single monolithic semiconductor waveguide. The first source is achieved through the concurrent utilization of two spontaneous parametric down-conversion (SPDC) processes, namely type-0 and type-I SPDC processes. The chip can also generate polarization entangled photons via the type-II SPDC process, enabling the generation of both co-polarized and cross-polarized entangled photons in the same device. In both cases, polarization entanglement is generated directly on the chip without the use of any off-chip compensation, interferometry or bandpass filtering. This enables direct, chip-based generation of both Bell states $(|H,H\rangle+|V,V\rangle)/\sqrt{2}$ and $(|H,V\rangle+|V,H\rangle)/\sqrt{2}$ simultaneously utilizing the same pump source. In addition, based on compound semiconductors, this chip can be directly integrated with it own pump laser. This technique ushers an era of self-contained, integrated, electrically pumped, room-temperature polarization entangled photon sources.


title: Network comparison using directed graphlets

abstract: With recent advances in high-throughput cell biology the amount of cellular biological data has grown drastically. Such data is often modeled as graphs (also called networks) and studying them can lead to new insights into molecule-level organization. A possible way to understand their structure is by analysing the smaller components that constitute them, namely network motifs and graphlets. Graphlets are particularly well suited to compare networks and to assess their level of similarity but are almost always used as small undirected graphs of up to five nodes, thus limiting their applicability in directed networks. However, a large set of interesting biological networks such as metabolic, cell signaling or transcriptional regulatory networks are intrinsically directional, and using metrics that ignore edge direction may gravely hinder information extraction. The applicability of graphlets is extended to directed networks by considering the edge direction of the graphlets. We tested our approach on a set of directed biological networks and verified that they were correctly grouped by type using directed graphlets. However, enumerating all graphlets in a large network is a computationally demanding task. Our implementation addresses this concern by using a state-of-the-art data structure, the g-trie, which is able to greatly reduce the necessary computation. We compared our tool, gtrieScanner, to other state-of-the art methods and verified that it is the fastest general tool for graphlet counting.


title: Sequential Detection of Market shocks using Risk-averse Agent Based Models

abstract: This paper considers a statistical signal processing problem involving agent based models of financial markets which at a micro-level are driven by socially aware and risk- averse trading agents. These agents trade (buy or sell) stocks by exploiting information about the decisions of previous agents (social learning) via an order book in addition to a private (noisy) signal they receive on the value of the stock. We are interested in the following: (1) Modelling the dynamics of these risk averse agents, (2) Sequential detection of a market shock based on the behaviour of these agents. Structural results which characterize social learning under a risk measure, CVaR (Conditional Value-at-risk), are presented and formulation of the Bayesian change point detection problem is provided. The structural results exhibit two interesting prop- erties: (i) Risk averse agents herd more often than risk neutral agents (ii) The stopping set in the sequential detection problem is non-convex. The framework is validated on data from the Yahoo! Tech Buzz game dataset.


title: Enhanced Low-Rank Matrix Approximation

abstract: This letter proposes to estimate low-rank matrices by formulating a convex optimization problem with non-convex regularization. We employ parameterized non-convex penalty functions to estimate the non-zero singular values more accurately than the nuclear norm. A closed-form solution for the global optimum of the proposed objective function (sum of data fidelity and the non-convex regularizer) is also derived. The solution reduces to singular value thresholding method as a special case. The proposed method is demonstrated for image denoising.


title: Diagonalization of the Finite Hilbert Transform on two adjacent intervals

abstract: We study the interior problem of tomography. The starting point is the Gelfand-Graev formula, which converts the tomographic data into the finite Hilbert transform (FHT) of an unknown function $f$ along a collection of lines. Pick one such line, call it the $x$-axis, and assume that the function to be reconstructed depends on a one-dimensional argument by restricting $f$ to the line. Let $\Omega_1$ be the interval where $f$ is supported, and $\Omega_2$ be the interval where the Hilbert transform of $f$ can be computed using the Gelfand-Graev formula. The equation we study is $H_1 f=g|_{\Omega_2}$, where $H_1$ is the FHT that integrates over $\Omega_1$ and gives the result on $\Omega_2$, i.e. $H_1: L^2(\Omega_1)\to L^2(\Omega_2)$. In the case of the interior problem the tomographic data are truncated, and $\Omega_1$ is no longer a subset of $\Omega_2$. In this paper we consider the case when the intervals $\Omega_1=(a_1,0)$ and $\Omega_2=(0,a_2)$ are adjacent. Here $a_1 < 0 < a_2$. First we find a differential operator $L$ that commutes with $H_1$. Using the Titchmarsh-Weyl theory, we show that $L$ has only continuous spectrum and obtain two isometric transformations $U_1$, $U_2$, such that $U_2 H_1 U_1^*$ is the multiplication operator with the function $\sigma(\lambda)$, $\lambda\geq(a_1^2+a_2^2)/8$. Here $\lambda$ is the spectral parameter. Then we show that $\sigma(\lambda)\to0$ as $\lambda\to\infty$ exponentially fast. We also obtain the leading asymptotic behavior of the kernels involved in the integral operators $U_1$, $U_2$ as $\lambda\to\infty$. When the intervals are symmetric, i.e. $-a_1=a_2$, the operators $U_1$, $U_2$ are obtained explicitly in terms of hypergeometric functions.


title: Lattice study for conformal windows of SU(2) and SU(3) gauge theories with fundamental fermions

abstract: We present our investigation of SU(2) gauge theory with 8 flavours, and SU(3) gauge theory with 12 flavours. For the SU(2) case, at strong bare coupling, $\beta \lesssim 1.45$, the distribution of the lowest eigenvalue of the Dirac operator can be described by chiral random matrix theory for the Gaussian symplectic ensemble. Our preliminary result indicates that the chiral phase transition in this theory is of bulk nature. For the SU(3) theory, we use high-precision lattice data to perform the step-scaling study of the coupling, $g_{{\rm GF}}$, in the Gradient Flow scheme. We carefully examine the reliability of the continuum extrapolation in the analysis, and conclude that the scaling behaviour of this SU(3) theory is not governed by possible infrared conformality at $g_{{\rm GF}}^{2} \lesssim 6$.


title: Energy Efficient Resource Allocation for Control Data Separation Architecture based H-CRAN with Heterogeneous Fronthaul

abstract: Control data separation architecture (CDSA) is a more efficient architecture to overcome the overhead issue than the conventional cellular networks, especially for the huge bursty traffic like Internet of Things, and over-the-top (OTT) content service. In this paper, we study the optimization issue of network energy efficiency of the CDSA-based heterogeneous cloud radio access networks (H-CRAN) networks, which has heterogeneous fronthaul between control base station (CBS) and data base stations (DBSs). We first present a modified power consumption model for the CDSA-based H-CRAN, and then formulate the optimization problem with constraint of overall capacity of wireless fronthaul. We work out the resource assignment and power allocation by the convex relaxation approach Using fractional programming method and Lagrangian dual decomposition method, we derive the close-form optimal solution and verify it by comprehensive system-level simulation. The simulation results show that our proposed algorithm has 8% EE gain compared to the static algorithm, and the CDSA-based H-CRAN networks can achieve up to 16% EE gain compared to the conventional network even under strict fronthaul capacity limit.


title: On a divisibility relation for Lucas sequences

abstract: In this note, we study the divisibility relation $U_m\mid U_{n+k}^s-U_n^s$, where ${\bf U}:=\{U_n\}_{n\ge 0}$ is the Lucas sequence of characteristic polynomial $x^2-ax\pm 1$ and $k,m,n,s$ are positive integers.


title: Observation of time-invariant coherence in a room temperature quantum simulator

abstract: The ability to live in coherent superpositions is a signature trait of quantum systems and constitutes an irreplaceable resource for quantum-enhanced technologies. However, decoherence effects usually destroy quantum superpositions. It has been recently predicted that, in a composite quantum system exposed to dephasing noise, quantum coherence in a transversal reference basis can stay protected for indefinite time. This can occur for a class of quantum states independently of the measure used to quantify coherence, and requires no control on the system during the dynamics. Here, such an invariant coherence phenomenon is observed experimentally in two different setups based on nuclear magnetic resonance at room temperature, realising an effective quantum simulator of two- and four-qubit spin systems. Our study further reveals a novel interplay between coherence and various forms of correlations, and highlights the natural resilience of quantum effects in complex systems.

categories: quant-ph cond-mat.mes-hall nucl-ex physics.data-an


title: Edge states and phase diagram for graphene under polarized light

abstract: We investigate the topological phase transitions in graphene under the modulation of circularly polarized light, by analyzing the changes of edge states and its topological structures. A full phase diagram, with up to ten different topological phases, is presented in the parameter space spanned by the driving frequency and light strength. We find that the high-Chern number behavior is very common in the driven system. While the one-photon resonance can create the chiral edge states in the $\pi-$gap, the two-photon resonance will induce the counter-propagating edge modes in the zero-energy gap. When the driving light strength is strong, the number and even the chirality of the edge states may change in the $\pi-$gap. The robustness of the edge states to disorder potential are also examined. We close by discussing the feasibility of experimental proposals.


title: Improving Covariate Balance in 2^K Factorial Designs via Rerandomization

abstract: Factorial designs are widely used in agriculture, engineering, and the social sciences to study the causal effects of several factors simultaneously on a response. The objective of such a design is to estimate all factorial effects of interest, which typically include main effects and interactions among factors. To estimate factorial effects with high precision when a large number of pre-treatment covariates are present, balance among covariates across treatment groups should be ensured. We propose utilizing rerandomization to ensure covariate balance in factorial designs. Although both factorial designs and rerandomization have been discussed before, the combination has not. Here, theoretical properties of rerandomization for factorial designs are established, and empirical results are explored using an application from the New York Department of Education.


title: Multi-lingual Geoparsing based on Machine Translation

abstract: Our method for multi-lingual geoparsing uses monolingual tools and resources along with machine translation and alignment to return location words in many languages. Not only does our method save the time and cost of developing geoparsers for each language separately, but also it allows the possibility of a wide range of language capabilities within a single interface. We evaluated our method in our LanguageBridge prototype on location named entities using newswire, broadcast news and telephone conversations in English, Arabic and Chinese data from the Linguistic Data Consortium (LDC). Our results for geoparsing Chinese and Arabic text using our multi-lingual geoparsing method are comparable to our results for geoparsing English text with our English tools. Furthermore, experiments using our machine translation approach results in accuracy comparable to results from the same data that was translated manually.


title: Persistence of centrality in random growing trees

abstract: We investigate properties of node centrality in random growing tree models. We focus on a measure of centrality that computes the maximum subtree size of the tree rooted at each node, with the most central node being the tree centroid. For random trees grown according to a preferential attachment model, a uniform attachment model, or a diffusion processes over a regular tree, we prove that a single node persists as the tree centroid after a finite number of steps, with probability 1. Furthermore, this persistence property generalizes to the top $K \ge 1$ nodes with respect to the same centrality measure. We also establish necessary and sufficient conditions for the size of an initial seed graph required to ensure persistence of a particular node with probability $1-\epsilon$, as a function of $\epsilon$: In the case of preferential and uniform attachment models, we derive bounds for the size of an initial hub constructed around the special node. In the case of a diffusion process over a regular tree, we derive bounds for the radius of an initial ball centered around the special node. Our necessary and sufficient conditions match up to constant factors for preferential attachment and diffusion tree models.


title: Distributed User Association in Energy Harvesting Small Cell Networks: A Competitive Market Model with Uncertainty

abstract: We consider a distributed user association problem in the downlink of a small cell network, where small cells obtain the required energy for providing wireless services to users through ambient energy harvesting. Since energy harvesting is opportunistic in nature, the amount of harvested energy is a random variable, without a priori known characteristics. We model the network as a competitive market with uncertainty, where self-interested small cells, modeled as consumers, are willing to maximize their utility scores by selecting users, represented by commodities. The utility scores of small cells depend on the amount of harvested energy, formulated as natures' state. Under this model, the problem is to assign users to small cells, so that the aggregate network utility is maximized. The solution is the general equilibrium under uncertainty, also called Arrow-Debreu equilibrium. We show that in our setting, such equilibrium not only exists, but also is unique and is Pareto optimal in the sense of expected aggregate network utility. We use the Walras' tatonnement process with some modifications in order to implement the equilibrium efficiently.


title: $h^0(125GeV) \to c \bar{c}$ as a test case for quark flavor violation in the MSSM

abstract: We calculate the decay width of $h^0(125GeV) \to c \bar{c}$ in the Minimal Supersymmetric Standard Model (MSSM) with non-minimal quark flavor violation (QFV) at full one-loop level. We adopt the $\overline{\rm DR}$ renormalization scheme. We study the effects of the mixing of the second and third squark generations (i.e. scharm-stop mixing) on the decay width, respecting the experimental constraints from B-meson data, the Higgs mass measurement and supersymmetric (SUSY) particle searches. We show that the decay width $\Gamma (h^0 \to c \bar{c})$ at the full one-loop level is very sensitive to the SUSY QFV parameters. In a scenario with large scharm-stop mixing, the decay width can differ up to $\sim \pm 35\%$ from its SM prediction. After taking into account the experimental and theoretical uncertainties of the decay width, we conclude that these QFV SUSY effects can be observed at a future $e^+ e^-$ collider such as ILC (International Linear Collider).


title: Characterizing star cluster formation with WISE: 652 newly found star clusters and candidates

abstract: We report the discovery of 652 star clusters, stellar groups and candidates in the Milky Way with WISE. Most of the objects are projected close to Galactic Plane and are embedded clusters. The present sample complements a similar study (Paper I) which provided 437 star clusters and alike. We find evidence that star formation processes span a wide range of sizes, from populous dense clusters to small compact embedded ones, sparse stellar groups or in relative isolation. The present list indicates multiple stellar generations during the embedded phase, with giant molecular clouds collapsing into several clumps composing an embedded cluster aggregate. We investigate the field star decontaminated Colour Magnitude Diagrams and Radial Density Profiles of 9 cluster candidates in the list, and derive their parameters, confirming them as embedded clusters.


title: An efficient time-stepping scheme for ab initio molecular dynamics simulations

abstract: In ab initio molecular dynamics simulations of real-world problems, the simple Verlet method is still widely used for integrating the equations of motion, while more efficient algorithms are routinely used in classical molecular dynamics. We show that if the Verlet method is used in conjunction with pre- and postprocessing, the accuracy of the time integration is significantly improved with only a small computational overhead. The validity of the processed Verlet method is demonstrated in several examples including ab initio molecular dynamics simulations of liquid water. The structural properties obtained from the processed Verlet method are found to be sufficiently accurate even for large time steps close to the stability limit. This approach results in a 2x performance gain over the standard Verlet method for a given accuracy.


title: Experimental Demonstration of Single Electron Transistors Featuring SiO2 PEALD in Ni-SiO2-Ni Tunnel Junctions

abstract: We report the use of plasma-enhanced atomic layer deposition (PEALD) to fabricate single-electron transistors (SETs) featuring ultra-thin (~1 nm) tunnel-transparent SiO2 in Ni-SiO2-Ni tunnel junctions. We show that as a result of the O2 plasma steps in PEALD of SiO2, the top surface of the underlying Ni electrode is oxidized. Additionally, the bottom surface of the upper Ni layer is also oxidized where it is in contact with the deposited SiO2, most likely as a result of oxygen-containing species on the surface of the SiO2. Due to the presence of these surface parasitic layers of NiO, which exhibit features typical of thermally activated transport, the resistance of Ni-SiO2-Ni tunnel junctions is drastically increased. Moreover, the transport mechanism is changed from quantum tunneling through the dielectric barrier to one consistent with thermally activated resistors in series with tunnel junctions. The reduction of NiO to Ni is therefore required to restore the metal-insulator-metal (MIM) structure of the junctions. Rapid thermal annealing in a forming gas ambient at elevated temperatures is presented as a technique to reduce both parasitic oxide layers. This method is of great interest for devices that rely on MIM tunnel junctions with ultra-thin barriers. Using this technique, we successfully fabricated MIM SETs with minimal trace of parasitic NiO component. We demonstrate that the properties of the tunnel barrier in nanoscale tunnel junctions can be evaluated by electrical characterization of SETs.


title: Statistical properties of $^{243}$Pu, and $^{242}$Pu(n,$\gamma$) cross section calculation

abstract: The level density and gamma-ray strength function (gammaSF) of 243Pu have been measured in the quasi-continuum using the Oslo method. Excited states in 243Pu were populated using the 242Pu(d,p) reaction. The level density closely follows the constant-temperature level density formula for excitation energies above the pairing gap. The gammaSF displays a double-humped resonance at low energy as also seen in previous investigations of actinide isotopes. The structure is interpreted as the scissors resonance and has a centroid of omega_{SR}=2.42(5)MeV and a total strength of B_{SR}=10.1(15)mu_N^2, which is in excellent agreement with sum-rule estimates. The measured level density and gammaSF were used to calculate the 242Pu(n,gamma) cross section in a neutron energy range for which there were previously no measured data.


title: Dynamic Response in a Finite Size Composite Multiferroic Thin Film

abstract: Composite multiferroics, heterostructures of ferromagnetic (FM) and ferroelectric (FE) materials, are characterized by a remarkable magnetoelectric effect at the interface. Previous work has supported the ferromagnetic structure with magnetic spins and the ferroelectric with pseudospins which act as electric dipoles in a microscopic model, coupled with a magnetoelectric interaction [J. Appl. Phys. 118, 124109 (2015)]. In this work, by solving the stochastic Landau-Lifshitz-Gilbert equation, the electric-field-induced magnetization switching in a twisted boundary condition has been studied, and a behavior of domain wall in the ferromagnetic structure is discussed.