title: Solar Flares and Coronal Mass Ejections: A Statistically Determined Flare Flux-CME Mass Correlation

abstract: In an effort to examine the relationship between flare flux and corresponding CME mass, we temporally and spatially correlate all X-ray flares and CMEs in the LASCO and GOES archives from 1996 to 2006. We cross-reference 6,733 CMEs having well-measured masses against 12,050 X-ray flares having position information as determined from their optical counterparts. For a given flare, we search in time for CMEs which occur 10-80 minutes afterward, and we further require the flare and CME to occur within +/-45 degrees in position angle on the solar disk. There are 826 CME/flare pairs which fit these criteria. Comparing the flare fluxes with CME masses of these paired events, we find CME mass increases with flare flux, following an approximately log-linear, broken relationship: in the limit of lower flare fluxes, log(CME mass)~0.68*log(flare flux), and in the limit of higher flare fluxes, log(CME mass)~0.33*log(flare flux). We show that this broken power-law, and in particular the flatter slope at higher flare fluxes, may be due to an observational bias against CMEs associated with the most energetic flares: halo CMEs. Correcting for this bias yields a single power-law relationship of the form log(CME mass)~0.70*log(flare flux). This function describes the relationship between CME mass and flare flux over at least 3 dex in flare flux, from ~10^-7 to 10^-4 W m^-2.


title: The Banana Project. III. Spin-orbit Alignment in the Long-period Eclipsing Binary NY Cephei

abstract: Binaries are not always neatly aligned. Previous observations of the DI Her system showed that the spin axes of both stars are highly inclined with respect to one another and the orbital axis. Here we report on a measurement of the spin-axis orientation of the primary star of the NY Cep system, which is similar to DI Her in many respects: it features two young early-type stars (~6 Myr, B0.5V+B2V), in an eccentric and relatively long-period orbit (e=0.48, P=15.d3). The sky projections of the rotation vector and the spin vector are well-aligned (beta_p = 2 +- 4 degrees), in strong contrast to DI Her. Although no convincing explanation has yet been given for the misalignment of DI Her, our results show that the phenomenon is not universal, and that a successful theory will need to account for the different outcome in the case of NY Cep.


title: HETDEX pilot survey for emission-line galaxies - I. Survey design, performance, and catalog

abstract: We present a catalog of emission-line galaxies selected solely by their emission-line fluxes using a wide-field integral field spectrograph. This work is partially motivated as a pilot survey for the upcoming Hobby-Eberly Telescope Dark Energy Experiment (HETDEX). We describe the observations, reductions, detections, redshift classifications, line fluxes, and counterpart information for 397 emission-line galaxies detected over 169 sq.arcmin with a 3500-5800 Ang. bandpass under 5 Ang. full-width-half-maximum (FWHM) spectral resolution. The survey's best sensitivity for unresolved objects under photometric conditions is between 4-20 E-17 erg/s/sq.cm depending on the wavelength, and Ly-alpha luminosities between 3-6 E42 erg/s are detectable. This survey method complements narrowband and color-selection techniques in the search for high redshift galaxies with its different selection properties and large volume probed. The four survey fields within the COSMOS, GOODS-N, MUNICS, and XMM-LSS areas are rich with existing, complementary data. We find 104 galaxies via their high redshift Ly-alpha emission at 1.9<z<3.8, and the majority of the remainder objects are low redshift [OII]3727 emitters at z<0.56. The classification between low and high redshift objects depends on rest frame equivalent width, as well as other indicators, where available. Based on matches to X-ray catalogs, the active galactic nuclei (AGN) fraction amongst the Ly-alpha emitters (LAEs) is 6%. We also analyze the survey's completeness and contamination properties through simulations. We find five high-z, highly-significant, resolved objects with full-width-half-maximum sizes >44 sq.arcsec which appear to be extended Ly-alpha nebulae. We also find three high-z objects with rest frame Ly-alpha equivalent widths above the level believed to be achievable with normal star formation, EW(rest)>240 Ang.


title: The Sudden Death of the Nearest Quasar

abstract: Galaxy formation is significantly modulated by energy output from supermassive black holes at the centers of galaxies which grow in highly efficient luminous quasar phases. The timescale on which black holes transition into and out of such phases is, however, unknown. We present the first measurement of the shutdown timescale for an individual quasar using X-ray observations of the nearby galaxy IC 2497, which hosted a luminous quasar no more than 70,000 years ago that is still seen as a light echo in `Hanny's Voorwerp', but whose present-day radiative output is lower by at least 2 and more likely by over 4 orders of magnitude. This extremely rapid shutdown provides new insights into the physics of accretion in supermassive black holes, and may signal a transition of the accretion disk to a radiatively inefficient state.


title: A Model for Thermal Phase Variations of Circular and Eccentric Exoplanets

abstract: We present a semi-analytic model atmosphere for close-in exoplanets that captures the essential physics of phase curves: orbital and viewing geometry, advection, and re-radiation. We calibrate the model with the well-characterized transiting planet, HD 189733b, then compute light curves for seven of the most eccentric transiting planets. We present phase variations for a variety of different radiative times and wind speeds. In the limit of instant re-radiation, the light curve morphology is entirely dictated by the planet's eccentricity and argument of pericenter: the light curve maximum leads or trails the eclipse depending on whether the planet is receding from or approaching the star at superior conjunction, respectively. For a planet with non-zero radiative timescales, the phase peak occurs early for super- rotating winds, and late for sub-rotating winds. We find that for a circular orbit, the timing of the phase variation maximum with respect to superior conjunction indicates the direction of the dominant winds, but cannot break the degeneracy between wind speed and radiative time. For circular planets the phase minimum occurs half an orbit away from the phase maximum -despite the fact that the coolest longitudes are always near the dawn terminator- and therefore does not convey any additional information. In general, increasing the advective frequency or the radiative time has the effect of reducing the peak-to-trough amplitude of phase variations, but there are interesting exceptions to these trends. Lastly, eccentric planets with orbital periods significantly longer than their radiative time exhibit "ringing" whereby the hot spot generated at periastron rotates in and out of view. The existence of ringing makes it possible to directly measure the wind speed (the frequency of the ringing) and the radiative time constant (the damping of the ringing).


title: Hot-Dust-Poor Quasars in Mid-Infrared and Optically Selected Samples

abstract: We show that the Hot-Dust-Poor (HDP) quasars, originally found in the X-ray selected XMM-COSMOS type 1 AGN sample, are just as common in two samples selected at optical/infrared wavelengths: the Richards et al. Spitzer/SDSS sample ($8.7%\pm2.2%$), and the PG-quasar dominated sample of Elvis et al. ($9.5%\pm 5.0%$). The properties of the HDP quasars in these two samples are consistent with the XMM-COSMOS sample, except that, at the $99% (\sim 2.5\sigma)$ significance, a larger proportion of the HDP quasars in the Spitzer/SDSS sample have weak host galaxy contributions, probably due to the selection criteria used. Either the host-dust is destroyed (dynamically or by radiation), or is offset from the central black hole due to recoiling. Alternatively, the universality of HDP quasars in samples with different selection methods and the continuous distribution of dust covering factor in type 1 AGNs, suggest that the range of SEDs could be related to the range of tilts in warped fueling disks, as in the model of Lawrence and Elvis (2010), with HDP quasars having relatively small warps.


title: The HETDEX Pilot Survey. II. The Evolution of the Ly-alpha Escape Fraction from the UV Slope and Luminosity Function of 1.9<z<3.8 LAEs

abstract: We study the escape of Ly-alpha photons from Ly-alpha emitting galaxies (LAEs) and the overall galaxy population using a sample of 99 LAEs at 1.9<z<3.8 detected through integral-field spectroscopy of blank fields by the HETDEX Pilot Survey. For 89 LAEs showing counterparts in deep broad-band images we measure the rest-frame UV luminosity and the UV slope, which we use to estimate E(B-V) under the assumption of a constant intrinsic UV slope for LAEs. These two quantities are used to measure the dust-corrected star formation rate (SFR). A comparison between the observed Ly-alpha luminosity and that predicted by the dust-corrected SFR yields the Ly-alpha escape fraction. We also measure the Ly-alpha luminosity function. Integration of the luminosity function provides a measurement of the Ly-alpha luminosity density across our redshift range. We combine our data with that from other surveys at 0.3<z<7.7 to trace the evolution of the Ly-alpha luminosity density. We then compare it to that expected from the star-formation history of the universe in order to characterize the evolution of the Ly-alpha escape fraction of the overall galaxy population [abriged]


title: The HETDEX Pilot Survey III: The Low Metallicities of High-Redshift Lyman Alpha Galaxies

abstract: We present Keck/NIRSPEC spectroscopic observations of three Lyman alpha emitting galaxies (LAEs) at z ~ 2.3 discovered with the HETDEX pilot survey. We detect Halpha, [OIII], and Hbeta emission from two galaxies at z = 2.29 and 2.49, designated HPS194 and HPS256, respectively, representing the first detection of multiple rest-frame optical emission lines in galaxies at high-redshift selected on the basis of their Lyman alpha emission. The redshifts of the Lyman alpha emission from these galaxies are offset redward of the systemic redshifts by Delta_v = 162 +/- 37 (photometric) +/- 42 (systematic) km/s for HPS194, and Delta_v = 36 +/- 35 +/- 18 km/s for HPS256. An interpretation for HPS194 is that a large-scale outflow may be occurring in its interstellar medium. The emission line ratios imply that neither LAE hosts an active galactic nucleus. Using the upper limits on the [NII] emission we place meaningful constraints on the gas-phase metallicities in these two LAEs of Z < 0.17 and < 0.28 Zsol (1 sigma). Measuring the stellar masses of these objects via spectral energy distribution (SED) fitting (~ 10^10 and 6 x 10^8 Msol, respectively), we study the nature of LAEs in a mass-metallicity plane. At least one of these two LAEs appears to be more metal poor than continuum-selected star-forming galaxies at the same redshift and stellar mass, implying that objects exhibiting Lyman alpha emission may be systematically less chemically enriched than the general galaxy population. We use the SEDs of these two galaxies to show that neglecting the emission lines when fitting stellar population models to the observed photometry can result in overestimates of the population age by orders of magnitude, and the stellar mass by a factor of ~ 2. This effect is particularly important at z > 7, where similarly strong emission lines may masquerade in the photometry as a 4000 A break (abridged).


title: Interpreting the Evolution of the Size - Luminosity Relation for Disk Galaxies from Redshift 1 to the Present

abstract: A sample of very high resolution cosmological disk galaxy simulations is used to investigate the evolution of galaxy disk sizes back to redshift 1 within the Lambda CDM cosmology. Artificial images in the rest frame B band are generated, allowing for a measurement of disk scale lengths using surface brightness profiles as observations would, and avoiding any assumption that light must follow mass as previous models have assumed. We demonstrate that these simulated disks are an excellent match to the observed magnitude - size relation for both local disks, and for disks at z=1 in the magnitude/mass range of overlap. We disentangle the evolution seen in the population as a whole from the evolution of individual disk galaxies. In agreement with observations, our simulated disks undergo roughly 1.5 magnitudes/arcsec^2 of surface brightness dimming since z=1. We find evidence that evolution in the magnitude - size plane varies by mass, such that galaxies with M* > 10^9 M_sun undergo more evolution in size than luminosity, while dwarf galaxies tend to evolve potentially more in luminosity. The disks grow in such a way as to stay on roughly the same stellar mass - size relation with time. Finally, due to an evolving stellar mass - SFR relation, a galaxy at a given stellar mass (or size) at z=1 will reside in a more massive halo and have a higher SFR, and thus a higher luminosity, than a counterpart of the same stellar mass at z=0.


title: Short-lived star-forming giant clumps in cosmological simulations of z~2 disks

abstract: Many observed massive star-forming z\approx2 galaxies are large disks that exhibit irregular morphologies, with \sim1kpc, \sim10^(8-10)Msun clumps. We present the largest sample to date of high-resolution cosmological SPH simulations that zoom-in on the formation of individual M*\sim10^(10.5)Msun galaxies in \sim10^(12)Msun halos at z\approx2. Our code includes strong stellar feedback parameterized as momentum-driven galactic winds. This model reproduces many characteristic features of this observed class of galaxies, such as their clumpy morphologies, smooth and monotonic velocity gradients, high gas fractions (f_g\sim50%) and high specific star-formation rates (\gtrsim1Gyr^(-1)). In accord with recent models, giant clumps (Mclump\sim(5x10^8-10^9)Msun) form in-situ via gravitational instabilities. However, the galactic winds are critical for their subsequent evolution. The giant clumps we obtain are short-lived and are disrupted by wind-driven mass loss. They do not virialise or migrate to the galaxy centers as suggested in recent work neglecting strong winds. By phenomenologically implementing the winds that are observed from high-redshift galaxies and in particular from individual clumps, our simulations reproduce well new observational constraints on clump kinematics and clump ages. In particular, the observation that older clumps appear closer to their galaxy centers is reproduced in our simulations, as a result of inside-out formation of the disks rather than inward clump migration.


title: Optimizing future experimental probes of inflation

abstract: The discovery of many novel realizations of the inflationary universe paradigm has led to a degeneracy problem: many different inflationary Lagrangians generate the same perturbation spectra. Resolving this problem requires the future discovery of additional observables, beyond the scalar adiabatic and tensor two-point functions on CMB scales. One important source of degeneracy arises in models where the density perturbation is generated by a non-inflationary degree of freedom, for example, through curvatons or modulated reheating. We consider the curvaton scenario as representative of this class, and analyze the degeneracy with single field, canonical inflation that results if the curvaton goes undetected by future observations. We perform Monte Carlo potential reconstructions in the absence of distinguishing observables, such as non-Gaussiantities or isocurvature modes. The resulting degeneracy is considerable and the improved measurements of spectral parameters from future probes like CMBPol, offer little to better the situation. Given a degeneracy-breaking observation, the observables must still be inverted to obtain the inflationary potential, with different observations resulting in reconstructions of varying quality. We find that a future detection of isocurvature modes or a precision measurement of the tensor spectral index will enable the most successful reconstructions in the presence of curvatons.


title: Asteroseismology of the Transiting Exoplanet Host HD 17156 with HST FGS

abstract: Observations conducted with the Fine Guidance Sensor on Hubble Space Telescope (HST) providing high cadence and precision time-series photometry were obtained over 10 consecutive days in December 2008 on the host star of the transiting exoplanet HD 17156b. During this time 10^12 photons (corrected for detector deadtime) were collected in which a noise level of 163 parts per million per 30 second sum resulted, thus providing excellent sensitivity to detection of the analog of the solar 5-minute p-mode oscillations. For HD 17156 robust detection of p-modes supports determination of the stellar mean density of 0.5301 +/- 0.0044 g/cm^3 from a detailed fit to the observed frequencies of modes of degree l = 0, 1, and 2. This is the first star for which direct determination of the mean stellar density has been possible using both asteroseismology and detailed analysis of a transiting planet light curve. Using the density constraint from asteroseismology, and stellar evolution modeling results in M_star = 1.285 +/- 0.026 solar, R_star = 1.507 +/- 0.012 solar, and a stellar age of 3.2 +/- 0.3 Gyr.


title: The Hubble Space Telescope Andromeda Treasury Survey. I. A public Star Catalog and Atlas based on the ACS/WFC data

abstract: As a legacy of the Hubble Space Telescope, a multi-cycle public survey has recently been approved and began to scan an area of approximately 0.5 degrees x 1.5 degrees in one quadrant of M31, using two cameras and several filters. The purpose of this work is to immediately release to the community a catalog of the sources imaged by ACS/WFC. The list will be updated on an almost real-time basis, as new data arrive, and the process will continue for the next two years, until the survey is completed. Each update will contain a chunk of 1/46 of the entire ACS/WFC survey, by means of astrometrized stacked images, and a catalog providing: positions, magnitudes, and several quality parameters.


title: The Gamma-Ray Background Constrains the Origins of the Radio and X-Ray Backgrounds

abstract: Cosmic ray protons generate gamma-rays, neutrinos, and secondary electrons and positrons (e+/-) through pion-producing collisions with gas atoms. Any synchrotron or Inverse Compton (IC) radiation from secondary e+/- is therefore accompanied by pionic gamma-rays. Using the extragalactic gamma-ray background, we constrain the contribution of secondary e+/- to the cosmic radio, X-ray, and soft gamma-ray backgrounds. These bounds depend on the unknown hadronic contribution to the gamma-ray background and the backgrounds' source redshifts. With our assumptions, we find that IC-upscattered light from secondaries is <~ 1/4 of the MeV - GeV gamma-ray background and <~ 10% of the 0.5 keV - 1 MeV background (for sources at a redshift z <~ 10). The low intensity of the observed gamma-ray background is marginally inconsistent with a secondary e+/- origin for the radio background reported by ARCADE at ~3 GHz, unless the magnetic field strength in their sources is milliGauss or greater. These limits on the magnetic field strength are sensitive to uncertainties. However, any contribution to the gamma-ray background from sources not responsible for the ARCADE excess increases the inconsistency.


title: The First Galaxies: Assembly of Disks and Prospects for Direct Detection

abstract: The James Webb Space Telescope (JWST) will enable observations of galaxies at redshifts z > 10 and hence allow to test our current understanding of structure formation at very early times. Previous work has shown that the very first galaxies inside halos with virial temperatures T < 10^4 K and masses M < 10^8 M_sun at z > 10 are probably too faint, by at least one order of magnitude, to be detected even in deep exposures with JWST. The light collected with JWST may therefore be dominated by radiation from galaxies inside ten times more massive halos. We use cosmological zoomed smoothed particle hydrodynamics simulations to investigate the assembly of such galaxies and assess their observability with JWST. We compare two simulations that are identical except for the inclusion of non-equilibrium H/D chemistry and radiative cooling by molecular hydrogen. In both simulations a large fraction of the halo gas settles in two nested, extended gas disks which surround a compact massive gas core. The presence of molecular hydrogen allows the disk gas to reach low temperatures and to develop marked spiral structure but does not qualitatively change its stability against fragmentation. We post-process the simulated galaxies by combining idealized models for star formation with stellar population synthesis models to estimate the luminosities in nebular recombination lines as well as in the ultraviolet continuum. We demonstrate that JWST will be able to constrain the nature of the stellar populations in galaxies such as simulated here based on the detection of the He1640 recombination line. Extrapolation of our results to halos with masses both lower and higher than those simulated shows that JWST may find up to a thousand star-bursting galaxies in future deep exposures of the z > 10 universe.


title: Improving Transit Predictions of Known Exoplanets with TERMS

abstract: Transiting planet discoveries have largely been restricted to the short-period or low-periastron distance regimes due to the bias inherent in the geometric transit probability. Through the refinement of planetary orbital parameters, and hence reducing the size of transit windows, long-period planets become feasible targets for photometric follow-up. Here we describe the TERMS project that is monitoring these host stars at predicted transit times.


title: Precise Estimates of the Physical Parameters for the Exoplanet System HD-17156 Enabled by HST FGS Transit and Asteroseismic Observations

abstract: We present observations of three distinct transits of HD 17156b obtained with the Fine Guidance Sensors (FGS) on board the Hubble Space Telescope} (HST). We analyzed both the transit photometry and previously published radial velocities to find the planet-star radius ratio R_p/R_s = 0.07454 +/- 0.00035, inclination i=86.49 +0.24/-0.20 deg, and scaled semi-major axis a/R = 23.19 +0.32/-0.27. This last value translates directly to a mean stellar density determination of 0.522 +0.021/-0.018 g cm^-3. Analysis of asteroseismology observations by the companion paper of Gilliland et al. (2009) provides a consistent but significantly refined measurement of the stellar mean density. We compare stellar isochrones to this density estimate and find M_s = 1.275 +/- 0.018 M_sun and a stellar age of $3.37 +0.20/-0.47 Gyr. Using this estimate of M_s and incorporating the density constraint from asteroseismology, we model both the photometry and published radial velocities to estimate the planet radius R_p= 1.0870 +/- 0.0066 Jupiter radii and the stellar radius R_s = 1.5007 +/- 0.0076 R_sun. The planet radius is larger than that found in previous studies and consistent with theoretical models of a solar-composition gas giant of the same mass and equilibrium temperature. For the three transits, we determine the times of mid-transit to a precision of 6.2 s, 7.6 s, and 6.9 s, and the transit times for HD 17156 do not show any significant departures from a constant period. The joint analysis of transit photometry and asteroseismology presages similar studies that will be enabled by the NASA Kepler Mission.


title: Type Ia Supernovae as Stellar Endpoints and Cosmological Tools

abstract: Empirically, Type Ia supernovae are the most useful, precise, and mature tools for determining astronomical distances. Acting as calibrated candles they revealed the presence of dark energy and are being used to measure its properties. However, the nature of the SN Ia explosion, and the progenitors involved, have remained elusive, even after seven decades of research. But now new large surveys are bringing about a paradigm shift --- we can finally compare samples of hundreds of supernovae to isolate critical variables. As a result of this, and advances in modeling, breakthroughs in understanding all aspects of SNe Ia are finally starting to happen.

categories: astro-ph.CO astro-ph.GA astro-ph.HE astro-ph.SR


title: Travelling-waves consistent with turbulence-driven secondary flow in a square duct

abstract: We present numerically determined travelling-wave solutions for pressure-driven flow through a straight duct with a square cross-section. This family of solutions represents typical coherent structures (a staggered array of counter-rotating streamwise vortices and an associated low-speed streak) on each wall. Their streamwise average flow in the cross-sectional plane corresponds to an eight vortex pattern much alike the secondary flow found in the turbulent regime.


title: Stellar X-ray sources in the Chandra COSMOS survey

abstract: We present an analysis of the X-ray properties of a sample of solar- and late-type field stars identified in the Chandra Cosmic Evolution Survey (COSMOS), a deep (160ks) and wide (0.9 deg2) extragalactic survey. The sample of 60 sources was identified using both morphological and photometric star/galaxy separation methods. We determine X-ray count rates, extract spectra and light curves and perform spectral fits to determine fluxes and plasma temperatures. Complementary optical and near-IR photometry is also presented and combined with spectroscopy for 48 of the sources to determine spectral types and distances for the sample. We find distances ranging from 30pc to ~12kpc, including a number of the most distant and highly active stellar X-ray sources ever detected. This stellar sample extends the known coverage of the L_X-distance plane to greater distances and higher luminosities, but we do not detect as many intrinsically faint X-ray sources compared to previous surveys. Overall the sample is typically more luminous than the active Sun, representing the high-luminosity end of the disk and halo X-ray luminosity functions. The halo population appears to include both low-activity spectrally hard sources that may be emitting through thermal bremsstrahlung, as well as a number of highly active sources in close binaries.


title: Graphene, Nobel Prize and All that Jazz

abstract: Graphene, a single atomic layer of graphite, first isolated in 2004, has made a quantum leap in the exploration of the physics of two-dimensional electron systems. Since the initial report of its discovery, many thousands of papers have been published, attempting to explain every aspect of the exotic electronic properties of this system. The graphene euphoria has culminated with the 2010 Nobel Prize in physics being awarded jointly to Andre Geim and Konstantin Novoselov of the University of Manchester, UK, "for groundbreaking experiments regarding the two-dimensional material graphene". But, what are the properties of graphene, and how was it made? Why it is so exciting for so many researchers, and why the Nobel Prize?


title: Global-scale wreath-building dynamos in stellar convection zones

abstract: When stars like our Sun are young they rotate rapidly and are very magnetically active. We explore dynamo action in rapidly rotating suns with the 3-D MHD anelastic spherical harmonic (ASH) code. The magnetic fields built in these dynamos are organized on global-scales into wreath-like structures that span the convection zone. Wreath-building dynamos can undergo quasi-cyclic reversals of polarity and such behavior is common in the parameter space we have been able to explore. These dynamos do not appear to require tachoclines to achieve their spatial or temporal organization. Wreath-building dynamos are present to some degree at all rotation rates, but are most evident in the more rapidly rotating simulations.


title: Chemistry of a protoplanetary disk with grain settling and Lyman alpha radiation

abstract: We present results from a model of the chemical evolution of protoplanetary disks. In our models we directly calculate the changing propagation and penetration of a high energy radiation field with Lyman alpha radiation included. We also explore the effect on our models of including dust grain settling. We find that, in agreement with earlier studies, the evolution of dust grains plays a large role in determining how deep the UV radiation penetrates into the disk. Significant grain settling at the midplane leads to much smaller freeze-out regions and a correspondingly larger molecular layer, which leads to an increase in column density for molecular species such as CO, CN and SO. The inclusion of Lyman alpha radiation impacts the disk chemistry through specific species that have large photodissociation cross sections at 1216 A. These include HCN, NH3 and CH4, for which the column densities are decreased by an order of magnitude or more due to the presence of Lyman alpha radiation in the UV spectrum. A few species, such as CO2 and SO, are enhanced by the presence of Lyman alpha radiation, but rarely by more than a factor of a few.


title: Finite Model Finding for Parameterized Verification

abstract: In this paper we investigate to which extent a very simple and natural "reachability as deducibility" approach, originated in the research in formal methods in security, is applicable to the automated verification of large classes of infinite state and parameterized systems. The approach is based on modeling the reachability between (parameterized) states as deducibility between suitable encodings of states by formulas of first-order predicate logic. The verification of a safety property is reduced to a pure logical problem of finding a countermodel for a first-order formula. The later task is delegated then to the generic automated finite model building procedures. In this paper we first establish the relative completeness of the finite countermodel finding method (FCM) for a class of parameterized linear arrays of finite automata. The method is shown to be at least as powerful as known methods based on monotonic abstraction and symbolic backward reachability. Further, we extend the relative completeness of the approach and show that it can solve all safety verification problems which can be solved by the traditional regular model checking.


title: Corrections to scaling in entanglement entropy from boundary perturbations

abstract: We investigate the corrections to scaling of the Renyi entropies of a region of size l at the end of a semi-infinite one-dimensional system described by a conformal field theory when the corrections come from irrelevant boundary operators. The corrections from irrelevant bulk operators with scaling dimension x have been studied by Cardy and Calabrese (2010), and they found not only the expected corrections of the form l^(4-2x) but also unusual corrections that could not have been anticipated by finite-size scaling arguments alone. However, for the case of perturbations from irrelevant boundary operators we find that the only corrections that can occur to leading order are of the form l^(2-2x_b) for boundary operators with scaling dimension x_b < 3/2, and l^(-1) when x_b > 3/2. When x_b=3/2 they are of the form l^(-1)log(l). A marginally irrelevant boundary perturbation will give leading corrections going as log(l)^(-3). No unusual corrections occur when perturbing with a boundary operator.


title: Identificaci\'on de nuevos medicamentos a trav\'es de m\'etodos computacionales

abstract: Resumen: El desarrollo de nuevos medicamentos es un problema complejo que carece de una soluci\'on \'unica y autom\'atica desde un punto de vista computacional, debido a la carencia de programas que permitan manejar grandes vol\'umenes de informaci\'on que est\'an distribuidos a lo largo de todo el mundo entre m\'ultiples bases de datos. Por ello se describe una metodolog\'ia que permita realizar experimentos in silico para la identificaci\'on actual de nuevos medicamentos. Abstract: The development of new drugs is a problem that nowadays has no solution in terms of computational power due to the lack of software for handling the big volume of available information; besides, these data are stored in multiple formats and are distributed all around the world. To resolve that, a development of an in silico drug design methodology.


title: From Sparse Signals to Sparse Residuals for Robust Sensing

abstract: One of the key challenges in sensor networks is the extraction of information by fusing data from a multitude of distinct, but possibly unreliable sensors. Recovering information from the maximum number of dependable sensors while specifying the unreliable ones is critical for robust sensing. This sensing task is formulated here as that of finding the maximum number of feasible subsystems of linear equations, and proved to be NP-hard. Useful links are established with compressive sampling, which aims at recovering vectors that are sparse. In contrast, the signals here are not sparse, but give rise to sparse residuals. Capitalizing on this form of sparsity, four sensing schemes with complementary strengths are developed. The first scheme is a convex relaxation of the original problem expressed as a second-order cone program (SOCP). It is shown that when the involved sensing matrices are Gaussian and the reliable measurements are sufficiently many, the SOCP can recover the optimal solution with overwhelming probability. The second scheme is obtained by replacing the initial objective function with a concave one. The third and fourth schemes are tailored for noisy sensor data. The noisy case is cast as a combinatorial problem that is subsequently surrogated by a (weighted) SOCP. Interestingly, the derived cost functions fall into the framework of robust multivariate linear regression, while an efficient block-coordinate descent algorithm is developed for their minimization. The robust sensing capabilities of all schemes are verified by simulated tests.


title: Equivalent Versions of "Khajuraho" and "Lo-Shu" Magic Squares and the day 1st October 2010 (01.10.2010)

abstract: In this short note we shall give connection between the most perfect "Khajuraho" magic square of order 4x4 discovered in 10th century and the "Lo-Shu" magic square of order 3x3 with the day October 1, 2010, i.e., 01.10.2010. The day has only three digits 0, 1 and 2. Here we have given an equivalent version of Khajuraho magic square using only three digits 0, 1 and 2. If we write the above date in two parts, 0110 2010, interestingly, the sum of new magic square is the first part, i.e., 0110, and the numbers appearing in the magic square are from the second part. An equivalent version of "Lo-Shu" magic square of order 3x3 is also given


title: Holographic Non-Gaussianity

abstract: We investigate the non-Gaussianity of primordial cosmological perturbations within our recently proposed holographic description of inflationary universes. We derive a holographic formula that determines the bispectrum of cosmological curvature perturbations in terms of correlation functions of a holographically dual three-dimensional non-gravitational quantum field theory (QFT). This allows us to compute the primordial bispectrum for a universe which started in a non-geometric holographic phase, using perturbative QFT calculations. Strikingly, for a class of models specified by a three-dimensional super-renormalisable QFT, the primordial bispectrum is of exactly the factorisable equilateral form with f_nl^eq=5/36, irrespective of the details of the dual QFT. A by-product of this investigation is a holographic formula for the three-point function of the trace of the stress-energy tensor along general holographic RG flows, which should have applications outside the remit of this work.


title: Ray Transforms on a Conformal Class of Curves

abstract: We introduce a technique for recovering a sufficiently smooth function from its ray transform over a wide class of curves in a general region of Euclidean space. The method is based on a complexification of the underlying vector fields defining the initial transport and recasting the problem in terms of complex-analytic function theory. Explicit inversion formulae are then given in a unified form. The method is then used to give inversion formulae for the attenuated ray transform.


title: Weitzenb\"ock derivations of nilpotency 3

abstract: We consider a Weitzenb\"ock derivation $\Delta$ acting on a polynomial ring $R=K[\xi_1,\xi_2,...,\xi_m]$ over a field $K$ of characteristic 0. The $K$-algebra $R^\Delta = \{h \in R \mid \Delta(h) = 0\}$ is called the algebra of constants. Nowicki considered the case where the Jordan matrix for $\Delta$ acting on $R_1$, the degree 1 component of $R$, has only Jordan blocks of size 2. He conjectured (\cite{N}) that a certain set generates $R^{\Delta}$ in that case. Recently Koury (\cite{Kh}), Drensky and Makar-Limanov (\cite{DM}) and Kuroda (\cite{K}) have given proofs of Nowicki's conjecture. Here we consider the case where the Jordan matrix for $\Delta$ acting on $R_{1}$ has only Jordan blocks of size at most 3. Here we use combinatorial methods to give a minimal set of generators $\mathcal G$ for the algebra of constants $R^{\Delta}$. Moreover, we show how our proof yields an algorithm to express any $h \in R^\Delta$ as a polynomial in the elements of $\mathcal G$. In particular, our solution shows how the classical techniques of polarization and restitution may be used to augment the techniques of SAGBI bases to construct generating sets for subalgebras.


title: Laser phase noise effects on the dynamics of optomechanical resonators

abstract: We investigate theoretically the influence of laser phase noise on the cooling and heating of a generic cavity optomechanical system. We derive the back-action damping and heating rates and the mechanical frequency shift of the radiation pressure-driven oscillating mirror, and derive the minimum phonon occupation number for small laser linewidths. We find that in practice laser phase noise does not pose serious limitations to ground state cooling. We then consider the effects of laser phase noise in a parametric cavity driving scheme that minimizes the back-action heating of one of the quadratures of the mechanical oscillator motion. Laser linewidths narrow compared to the decay rate of the cavity field will not pose any problems in an experimental setting, but broader linewidths limit the practicality of this back-action evasion method.


title: Lectures on Scattering Amplitudes in String Theory

abstract: In these lecture notes, we take a closer look at the calculation of scattering amplitudes for the bosonic string. It is believed that string theories form the UV completions of (super)gravity theories. Support for this claim can be found in the (on-shell) scattering amplitudes of strings. On the other hand, studying these string scattering amplitudes opens a window on the UV behavior of the string theories themselves. In these short set of lectures, we discuss the two-dimensional Polyakov path integral for the string, and its gauge symmetries, the connection to Riemann surfaces and how to obtain some of the simplest string scattering amplitudes. We end with some comments on more advanced topics. For simplicity we limit ourselves to bosonic open string theory in 26 dimensions.


title: The WFC3 Galactic Bulge Treasury Program: Metallicity Estimates for the Stellar Population and Exoplanet Hosts

abstract: We present new UV-to-IR stellar photometry of four low-extinction windows in the Galactic bulge, obtained with the Wide Field Camera 3 on the Hubble Space Telescope (HST). Using our five bandpasses, we have defined reddening-free photometric indices sensitive to stellar effective temperature and metallicity. We find that the bulge populations resemble those formed via classical dissipative collapse: each field is dominated by an old (~10 Gyr) population exhibiting a wide metallicity range (-1.5 < [Fe/H] < 0.5). We detect a metallicity gradient in the bulge population, with the fraction of stars at super-solar metallicities dropping from 41% to 35% over distances from the Galactic center ranging from 0.3 to 1.2 kpc. One field includes candidate exoplanet hosts discovered in the SWEEPS HST transit survey. Our measurements for 11 of these hosts demonstrate that exoplanets in the distinct bulge environment are preferentially found around high-metallicity stars, as in the solar neighborhood, supporting the view that planets form more readily in metal-rich environments.


title: Leverage Bubble

abstract: Leverage is strongly related to liquidity in a market and lack of liquidity is considered a cause and/or consequence of the recent financial crisis. A repurchase agreement is a financial instrument where a security is sold simultaneously with an agreement to buy it back at a later date. Repurchase agreements (repos) market size is a very important element in calculating the overall leverage in a financial market. Therefore, studying the behavior of repos market size can help to understand a process that can contribute to the birth of a financial crisis. We hypothesize that herding behavior among large investors led to massive over-leveraging through the use of repos, resulting in a bubble (built up over the previous years) and subsequent crash in this market in early 2008. We use the Johansen-Ledoit-Sornette (JLS) model of rational expectation bubbles and behavioral finance to study the dynamics of the repo market that led to the crash. The JLS model qualifies a bubble by the presence of characteristic patterns in the price dynamics, called log-periodic power law (LPPL) behavior. We show that there was significant LPPL behavior in the market before that crash and that the predicted range of times predicted by the model for the end of the bubble is consistent with the observations.


title: Non-thermal radiation associated with astrophysical shock waves

abstract: The main goal of this thesis is to study the physical processes that can produce non-thermal emission at high energies in astrophysical objects capable to accelerate particles up to relativistic velocities. In particular, we have studied the gamma-ray emission produced in cosmic sources with different spatial scales, from young stellar objects to clusters of galaxies, going through microquasars and active galactic nuclei. In the former cases, we have modeled the gamma-ray emission using the radio data from the sources IRAS 16547-4247 and Abell 3376. In the latter, we have developed a specific radiation model based on the interaction of the inhomogeneities of the external medium with the jets generated by the compact object. Specifically, we have considered clumps of the massive stellar wind in microquasars, and clouds of the broad line region in active galactic nuclei, interacting with the jets of the sources. In all cases, the developed models allow us to make predictions testables with the new generation of instruments operating at high energies, such as the satellites Fermi and AGILE, and the Cherenkov telescopes HESS, MAGIC, and the forthcoming CTA.


title: Ultraproducts of Tannakian Categories and Generic Representation Theory of Unipotent Algebraic Groups

abstract: The principle of tannakian duality states that any neutral tannakian category is tensorially equivalent to the category Rep_k G of finite dimensional representations of some affine group scheme G and field k, and conversely. Originally motivated by an attempt to find a first-order explanation for generic cohomology of algebraic groups, we study neutral tannakian categories as abstract first-order structures and, in particular, ultraproducts of them. One of the main theorems of this dissertation is that certain naturally definable subcategories of these ultraproducts are themselves neutral tannakian categories, hence tensorially equivalent to Comod_A for some Hopf algebra A over a field k. We are able to give a fairly tidy description of the representing Hopf algebras of these categories, and explicitly compute them in several examples. For the second half of this dissertation we turn our attention to the representation theories of certain unipotent algebraic groups, namely the additive group G_a and the Heisenberg group H_1. The results we obtain for these groups in characteristic zero are not at all new or surprising, but in positive characteristic they perhaps are. In both cases we obtain that, for a given dimension n, if p is large enough with respect to n, all n-dimensional modules for these groups in characteristic p are given by commuting products of representations, with the constituent factors resembling representations of the same group in characteristic zero. We later use these results to extrapolate some generic cohomology results for these particular unipotent groups.


title: Triangle groups, automorphic forms, and torus knots

abstract: This paper's theme is the relation between several classical and well-known objects: triangle Fuchsian groups, quasi-homogeneous singularities of plane curves, torus knot complements in the 3-sphere. Torus knots are the only nontrivial knots whose complements admit transitive Lie group actions. In fact S^3\K_{p,q} is diffeomorphic to a coset space of the universal covering group of PSL_2(R) with respect to a discrete subgroup G contained in the preimage of a (p,q,\infty)-triangle Fuchsian group. The existence of such a diffeomorphism between is known from a general topological classification of Seifert fibred 3-manifolds. Our goal is to construct an explicit diffeomorphism using automorphic forms. Such a construction is previously known for the trefoil knot K_{2,3} and in fact S^3\K_{2,3} = SL_2(R)/SL_2(Z). The connection between the two sides of the diffeomorphism comes via singularities of plane curves.


title: Poisson smooth structures on stratified symplectic spaces

abstract: In this paper we introduce the notion of a smooth structure on a stratified space, the notion of a Poisson smooth structure and the notion of a weakly symplectic smooth structure on a stratified symplectic space, refining the concept of a stratified symplectic Poisson algebra introduced by Sjamaar and Lerman. We show that these smooth spaces possess several important properties, e.g. the existence of smooth partitions of unity. Furthermore, under mild conditions many properties of a symplectic manifold can be extended to a symplectic stratified space provided with a smooth Poisson structure, e.g. the existence and uniqueness of a Hamiltonian flow, the isomorphism between the Brylinski-Poisson homology and the de Rham homology, the existence of a Leftschetz decomposition on a symplectic stratified space. We give many examples of stratified symplectic spaces possessing a Poisson smooth structure which is also weakly symplectic.


title: Inertial forces and dissipation on accelerated boundaries

abstract: We study dissipative effects due to inertial forces acting on matter fields confined to accelerated boundaries in $1+1$, $2+1$, and $3+1$ dimensions. These matter fields describe the internal degrees of freedom of `mirrors' and impose, on the surfaces where they are defined, boundary conditions on a fluctuating `vacuum' field. We construct different models, involving either scalar or Dirac matter fields coupled to a vacuum scalar field, and use effective action techniques to calculate the strength of dissipation. In the case of massless Dirac fields, the results could be used to describe the inertial forces on an accelerated graphene sheet.


title: Computing local constants for CM elliptic curves

abstract: For E/k an elliptic curve with CM by O, we determine a formula for (a generalization of) the arithmetic local constant of [4] at almost all primes of good reduction. We apply this formula to the CM curves defined over Q and are able to describe extensions F/Q over which the O-rank of E grows.


title: The interactions of winds from massive young stellar objects

abstract: The supersonic stellar and disk winds possessed by massive young stellar objects will produce shocks when they collide against the interior of a pre-existing bipolar cavity (resulting from an earlier phase of jet activity). The shock heated gas emits thermal X-rays which may be observable by spaceborne observa- tories such as the Chandra X-ray Observatory. Hydrodynamical models are used to explore the wind-cavity interaction. Radiative transfer calculations are performed on the simulation output to produce synthetic X-ray observations, allowing constraints to be placed on model parameters through comparisons with observations. The model reveals an intricate interplay between the inflowing and outflowing material and is successful in reproducing the observed X-ray count rates from massive young stellar objects.


title: Fluctuation-driven Turing patterns

abstract: Models of diffusion driven pattern formation that rely on the Turing mechanism are utilized in many areas of science. However, many such models suffer from the defect of requiring fine tuning of parameters or an unrealistic separation of scales in the diffusivities of the constituents of the system in order to predict the formation of spatial patterns. In the context of a very generic model of ecological pattern formation, we show that the inclusion of intrinsic noise in Turing models leads to the formation of "quasi-patterns" that form in generic regions of parameter space and are experimentally distinguishable from standard Turing patterns. The existence of quasi-patterns removes the need for unphysical fine tuning or separation of scales in the application of Turing models to real systems.


title: Population imbalance and pairing in the BCS-BEC crossover of three-component ultracold fermions

abstract: We investigate the phase diagram and the BCS-BEC crossover of a homogeneous three-component ultracold Fermi gas with a U(3) invariant attractive interaction. We show that the system at sufficiently low temperatures exhibits population imbalance, as well as fermionic pairing. We describe the crossover in this system, connecting the weakly interacting BCS regime of the partially population-imbalanced fermion pairing state and the BEC limit with three weakly interacting species of molecules, including pairing fluctuations within a t-matrix calculation of the particle self-energies.


title: Efficient Triangle Counting in Large Graphs via Degree-based Vertex Partitioning

abstract: The number of triangles is a computationally expensive graph statistic which is frequently used in complex network analysis (e.g., transitivity ratio), in various random graph models (e.g., exponential random graph model) and in important real world applications such as spam detection, uncovering of the hidden thematic structure of the Web and link recommendation. Counting triangles in graphs with millions and billions of edges requires algorithms which run fast, use small amount of space, provide accurate estimates of the number of triangles and preferably are parallelizable. In this paper we present an efficient triangle counting algorithm which can be adapted to the semistreaming model. The key idea of our algorithm is to combine the sampling algorithm of Tsourakakis et al. and the partitioning of the set of vertices into a high degree and a low degree subset respectively as in the Alon, Yuster and Zwick work treating each set appropriately. We obtain a running time $O \left(m + \frac{m^{3/2} \Delta \log{n}}{t \epsilon^2} \right)$ and an $\epsilon$ approximation (multiplicative error), where $n$ is the number of vertices, $m$ the number of edges and $\Delta$ the maximum number of triangles an edge is contained. Furthermore, we show how this algorithm can be adapted to the semistreaming model with space usage $O\left(m^{1/2}\log{n} + \frac{m^{3/2} \Delta \log{n}}{t \epsilon^2} \right)$ and a constant number of passes (three) over the graph stream. We apply our methods in various networks with several millions of edges and we obtain excellent results. Finally, we propose a random projection based method for triangle counting and provide a sufficient condition to obtain an estimate with low variance.


title: Classical Topological Order in Kagome Ice

abstract: We examine the onset of classical topological order in a nearest-neighbor kagome ice model. Using Monte Carlo simulations, we characterize the topological sectors of the groundstate using a non-local cut measure which circumscribes the toroidal geometry of the simulation cell. We demonstrate that simulations which employ global loop updates that are allowed to wind around the periodic boundaries cause the topological sector to fluctuate, while restricted local loop updates freeze the simulation into one topological sector. The freezing into one topological sector can also be observed in the susceptibility of the real magnetic spin vectors projected onto the kagome plane. The ability of the susceptibility to distinguish between fluctuating and non-fluctuating topological sectors should motivate its use as a local probe of topological order in a variety of related model and experimental systems.


title: Astronomy with Radioactivities: Chapter 9, Nuclear Reactions

abstract: Nuclear reaction rates determine the abundances of isotopes in stellar burning processes. A multitude of reactions determine the reaction flow pattern which is described in terms of reaction network simulations. The reaction rates are determined by laboratory experiments supplemented by nuclear reaction and structure theory. We will discuss the experimental approach as well as the theoretical tools for obtaining the stellar reaction rates. A detailed analysis of a reaction is only possible for a few selected cases which will be highlighted in this section. The bulk of nuclear reaction processes is however described in terms of a statistical model approach, which relies on global nuclear structure and reaction parameters such as level density and mass and barrier penetration, respectively. We will discuss a variety of experimental facilities and techniques used in the field, this includes low energy stable beam experiments, measurements at radioactive beam accelerators, and neutron beam facilities.

categories: astro-ph.IM astro-ph.HE astro-ph.SR nucl-ex nucl-th


title: A Note on an R^2 Measure for Fixed Effects in the Generalized Linear Mixed Model

abstract: Using the LRT statistic, a model R^2 is proposed for the generalized linear mixed model for assessing the association between the correlated outcomes and fixed effects. The R^2 compares the full model to a null model with all fixed effects deleted.


title: Regularized Risk Minimization by Nesterov's Accelerated Gradient Methods: Algorithmic Extensions and Empirical Studies

abstract: Nesterov's accelerated gradient methods (AGM) have been successfully applied in many machine learning areas. However, their empirical performance on training max-margin models has been inferior to existing specialized solvers. In this paper, we first extend AGM to strongly convex and composite objective functions with Bregman style prox-functions. Our unifying framework covers both the $\infty$-memory and 1-memory styles of AGM, tunes the Lipschiz constant adaptively, and bounds the duality gap. Then we demonstrate various ways to apply this framework of methods to a wide range of machine learning problems. Emphasis will be given on their rate of convergence and how to efficiently compute the gradient and optimize the models. The experimental results show that with our extensions AGM outperforms state-of-the-art solvers on max-margin models.


title: Coupled quantized mechanical oscillators

abstract: The harmonic oscillator is one of the simplest physical systems but also one of the most fundamental. It is ubiquitous in nature, often serving as an approximation for a more complicated system or as a building block in larger models. Realizations of harmonic oscillators in the quantum regime include electromagnetic fields in a cavity [1-3] and the mechanical modes of a trapped atom [4] or macroscopic solid [5]. Quantized interaction between two motional modes of an individual trapped ion has been achieved by coupling through optical fields [6], and entangled motion of two ions in separate locations has been accomplished indirectly through their internal states [7]. However, direct controllable coupling between quantized mechanical oscillators held in separate locations has not been realized previously. Here we implement such coupling through the mutual Coulomb interaction of two ions held in trapping potentials separated by 40 um (similar work is reported in a related paper [8]). By tuning the confining wells into resonance, energy is exchanged between the ions at the quantum level, establishing that direct coherent motional coupling is possible for separately trapped ions. The system demonstrates a building block for quantum information processing and quantum simulation. More broadly, this work is a natural precursor to experiments in hybrid quantum systems, such as coupling a trapped ion to a quantized macroscopic mechanical or electrical oscillator [9-13].


title: Construction of New Delay-Tolerant Space-Time Codes

abstract: Perfect Space-Time Codes (STC) are optimal codes in their original construction for Multiple Input Multiple Output (MIMO) systems. Based on Cyclic Division Algebras (CDA), they are full-rate, full-diversity codes, have Non-Vanishing Determinants (NVD) and hence achieve Diversity-Multiplexing Tradeoff (DMT). In addition, these codes have led to optimal distributed space-time codes when applied in cooperative networks under the assumption of perfect synchronization between relays. However, they loose their diversity when delays are introduced and thus are not delay-tolerant. In this paper, using the cyclic division algebras of perfect codes, we construct new codes that maintain the same properties as perfect codes in the synchronous case. Moreover, these codes preserve their full-diversity in asynchronous transmission.


title: The magnetic field in the NGC 2024 FIR 5 dense core

abstract: We used the Submillimeter Array (SMA) to observe the thermal polarized dust emission from the protostellar source NGC 2024 FIR 5. The polarized emission outlines a partial hourglass morphology for the plane-of-sky component of the core magnetic field. Our data are consistent with previous BIMA maps, and the overall magnetic field geometries obtained with both instruments are similar. We resolve the main core into two components, FIR 5A and FIR 5B. A possible explanation for the asymmetrical field lies in depolarization effects due to the lack of internal heating from FIR 5B source, which may be in a prestellar evolutionary state. The field strength was estimated to be 2.2 mG, in agreement with previous BIMA data. We discuss the influence of a nearby H{\sc ii} region over the field lines at scales of $\sim 0.01$ pc. Although the hot component is probably compressing the molecular gas where the dust core is embedded, it is unlikely that the radiation pressure exceeds the magnetic tension. Finally, a complex outflow morphology is observed in CO (3 $\rightarrow$ 2) maps. Unlike previous maps, several features associated with dust condensations other than FIR 5 are detected.


title: NICMOS transmission spectroscopy of HD 189733b: controversy becomes confirmation

abstract: Spectral features corresponding to methane and water opacity were reported based on spectroscopic observations of HD 189733b with Hubble/NICMOS. Recently, these data, and other NICMOS exoplanet spectroscopy measurements, have been reexamined in Gibson et al. 2010, who claim that the features in the transmission spectra are due to uncorrected systematic errors and not molecular opacities. We examine the methods used by the Gibson team and show that, contrary to their claim, their results for the transmission spectrum of HD 189733b are in fact in agreement with the original results. In the case of HD 189733b, the most significant problem with the Gibson approach is a poorly determined instrument model, which causes (1) an increase in the formal uncertainty and (2) instability in the minimization process; although Gibson et al. do recover the correct spectrum, they cannot identify it due to the problems caused by a poorly determined instrument model. In the case of XO-1b, the Gibson method is fundamentally flawed because they omit the most important parameters from the instrument model. For HD 189733b, the Gibson team did not omit these parameters, which explains why they are able to reproduce previous results in this case, although with poor SNR.


title: Ross Ice Shelf in situ radio-frequency ice attenuation

abstract: We have measured the in situ average electric field attenuation length for radio-frequency signals broadcast vertically through the Ross Ice Shelf. We chose a location, Moore Embayment, south of Minna Bluff, known for its high reflectivity at the ice-sea interface. We confirmed specular reflection and used the return pulses to measure the average attenuation length from 75-1250 MHz over the round-trip distance of 1155 m. We find the average electric field attenuation length to vary from 500 m at 75 MHz to 300 m at 1250 MHz, with an experimental uncertainty of 55 to 15 m. We discuss the implications for neutrino telescopes that use the radio technique and include the Ross Ice Shelf as part of their sensitive volume.


title: Preserving multiple first integrals by discrete gradients

abstract: We consider systems of ordinary differential equations with known first integrals. The notion of a discrete tangent space is introduced as the orthogonal complement of an arbitrary set of discrete gradients. Integrators which exactly conserve all the first integrals simultaneously are then defined. In both cases we start from an arbitrary method of a prescribed order (say, a Runge-Kutta scheme) and modify it using two approaches: one based on projection and one based one local coordinates. The methods are tested on the Kepler problem.


title: Kerr-de Sitter Universe

abstract: It is now widely accepted that the universe as we understand it is accelerating in expansion and fits the de Sitter model rather well. As such, a realistic assumption of black holes must place them on a de Sitter background and not Minkowski as is typically done in General Relativity. The most astrophysically relevant black hole is the uncharged, rotating Kerr solution, a member of the more general Kerr-Newman metrics. A generalization of the rotating Kerr black hole to a solution of the Einstein's equation with a cosmological constant $\Lambda$ was discovered by Carter \cite{DWDW}. It is typically referred to as the Kerr-de Sitter spacetime. Here, we discuss the horizon structure of this spacetime and its dependence on $\Lambda$. We recall that in a $\La>0$ universe, the term `extremal black hole' refers to a black hole with angular momentum $J > M^2 $. We obtain explicit numerical results for the black hole's maximal spin value and get a distribution of admissible Kerr holes in the ($\Lambda$, spin) parameter space. We look at the conformal structure of the extended spacetime and the embedding of the 3-geometry of the spatial hypersurfaces. In analogy with Reissner-Nordstr\"{o}m -de Sitter spacetime, in particular by considering the Kerr-de Sitter causal structure as a distortion of the Reissner-Nordstr\"{o}m-de Sitter one, we show that spatial sections of the extended spacetime are 3-spheres containing 2-dimensional topologically spherical sections of the horizons of Kerr holes at the poles. Depending on how a $t=$ constant 3-space is defined these holes may be seen as black or white holes (four possible combinations).


title: Excited States of the Nucleon in 2+1 Flavor QCD

abstract: Recent developments on the determination of the spin-1/2 spectrum of the nucleon in full QCD are presented. Our focus is on the PACS-CS 2+1 flavor configurations made available through the ILDG. Using correlation matrix techniques, in which a wide variety of gauge-invariant Gaussian-smeared fermion-propagator sources and sinks are considered, excited states are determined. We consider several correlation matrices of various sizes, each constructed with a different set of basis interpolators, in order to demonstrate the invariance of the eigenstates on the basis choice. Of particular interest is the approach to the elusive Roper resonance and we report preliminary results in full QCD.


title: A novel method for evaluating correlation functions in lattice hadron spectroscopy

abstract: We describe a new approach for evaluating hadronic correlation functions which combines Laplacian-Heaviside quark smearing with a stochastic estimator of quark propagators. This method utilizes noise dilution in a new way to reduce the variance in correlators. The efficacy of the new algorithm is demonstrated on a number of systems, including disconnected diagrams and multi-hadron correlators, on a small lattice where comparisons with the results obtained with exactly determined quark propagators are possible. On larger lattice volumes, the use of exact propagators becomes prohibitively expensive, while the stochastic method is still computationally feasible.


title: Quantum Incompressibility of a Freely Falling Hydrogen Atom in a Circular Rydberg State, and a Gravitationally-Induced Charge Separation Effect in Superconducting Systems

abstract: Freely falling point-like objects converge toward the center of the Earth. Hence the gravitational field of the Earth is inhomogeneous, and possesses a tidal component. The free fall of an extended quantum mechanical object such as a hydrogen atom prepared in a high principal-quantum-number state, i.e. a circular Rydberg atom, is predicted to fall more slowly than a classical point-like object, when both objects are dropped from the same height above the Earth's surface. This indicates that, apart from transitions between quantum states, the atom exhibits a kind of quantum mechanical incompressibility during free fall in inhomogeneous, tidal gravitational fields like those of the Earth. A superconducting ring-like system with a persistent current circulating around it behaves like the circular Rydberg atom during free fall. Like the electronic wavefunction of the freely falling atom, the Cooper-pair wavefunction is quantum mechanically incompressible. The ions in the lattice of the superconductor, however, are not incompressible, since they do not possess a globally coherent quantum phase. The resulting difference during free fall in the response of the nonlocalizable Cooper pairs of electrons and the localizable ions to inhomogeneous gravitational fields is predicted to lead to a charge separation effect, which in turn leads to a large Coulomb force that opposes the convergence caused by the tidal gravitational force on the superconducting system. A "Cavendish-like" experiment is proposed for observing the charge separation effect induced by inhomogeneous gravitational fields in a superconducting circuit. The charge separation effect is determined to be limited by a pair-breaking process that occurs when low frequency gravitational perturbations are present.


title: Unquenching the Landau Gauge Lattice Propagators and the Gribov Problem

abstract: The gluon and ghost propagators are computed using both quenched and dynamical configurations for the same lattice spacings. The Wilson fermions are simulated at several quark masses. Furthermore, the effect of the Gribov copies is evaluated for all sets of configurations.


title: Effect of oxygen plasma etching on graphene studied with Raman spectroscopy and electronic transport

abstract: We report a study of graphene and graphene field effect devices after exposure to a series of short pulses of oxygen plasma. We present data from Raman spectroscopy, back-gated field-effect and magneto-transport measurements. The intensity ratio between Raman "D" and "G" peaks, I(D)/I(G) (commonly used to characterize disorder in graphene) is observed to increase approximately linearly with the number (N(e)) of plasma etching pulses initially, but then decreases at higher Ne. We also discuss implications of our data for extracting graphene crystalline domain sizes from I(D)/I(G). At the highest Ne measured, the "2D" peak is found to be nearly suppressed while the "D" peak is still prominent. Electronic transport measurements in plasma-etched graphene show an up-shifting of the Dirac point, indicating hole doping. We also characterize mobility, quantum Hall states, weak localization and various scattering lengths in a moderately etched sample. Our findings are valuable for understanding the effects of plasma etching on graphene and the physics of disordered graphene through artificially generated defects.


title: Variable Order Mixed H-Finite Element Method for Linear Elasticity with Weakly Imposed Symmetry. Ii. Affine and Curvilinear Elements in 2D

abstract: We continue our study on variable order Arnold-Falk-Winther elements for 2D elasticity in context of both affine and parametric curvilinear elements. We present an $h$-stability result for affine elements, and an asymptotic stability result for curvilinear elements. Both theoretical results are confirmed with numerical experiments.


title: Hadron plus photon production in polarized hadronic collisions at next-to-leading order accuracy

abstract: We compute the next-to-leading order QCD corrections to the polarized (and unpolarized) cross sections for the production of a hadron accompanied by an opposite-side prompt photon. This process, being studied at RHIC, permits us to reconstruct partonic kinematics using experimentally measurable variables. We study the correlation between the reconstructed momentum fractions and the true partonic ones, which in the polarized case might allow us to reveal the spin-dependent gluon distribution with a higher precision.


title: Stochastic Simulation of Process Calculi for Biology

abstract: Biological systems typically involve large numbers of components with complex, highly parallel interactions and intrinsic stochasticity. To model this complexity, numerous programming languages based on process calculi have been developed, many of which are expressive enough to generate unbounded numbers of molecular species and reactions. As a result of this expressiveness, such calculi cannot rely on standard reaction-based simulation methods, which require fixed numbers of species and reactions. Rather than implementing custom stochastic simulation algorithms for each process calculus, we propose to use a generic abstract machine that can be instantiated to a range of process calculi and a range of reaction-based simulation algorithms. The abstract machine functions as a just-in-time compiler, which dynamically updates the set of possible reactions and chooses the next reaction in an iterative cycle. In this short paper we give a brief summary of the generic abstract machine, and show how it can be instantiated with the stochastic simulation algorithm known as Gillespie's Direct Method. We also discuss the wider implications of such an abstract machine, and outline how it can be used to simulate multiple calculi simultaneously within a common framework.


title: Measurable Stochastics for Brane Calculus

abstract: We give a stochastic extension of the Brane Calculus, along the lines of recent work by Cardelli and Mardare. In this presentation, the semantics of a Brane process is a measure of the stochastic distribution of possible derivations. To this end, we first introduce a labelled transition system for Brane Calculus, proving its adequacy w.r.t. the usual reduction semantics. Then, brane systems are presented as Markov processes over the measurable space generated by terms up-to syntactic congruence, and where the measures are indexed by the actions of this new LTS. Finally, we provide a SOS presentation of this stochastic semantics, which is compositional and syntax-driven.


title: An Abstraction Theory for Qualitative Models of Biological Systems

abstract: Multi-valued network models are an important qualitative modelling approach used widely by the biological community. In this paper we consider developing an abstraction theory for multi-valued network models that allows the state space of a model to be reduced while preserving key properties of the model. This is important as it aids the analysis and comparison of multi-valued networks and in particular, helps address the well-known problem of state space explosion associated with such analysis. We also consider developing techniques for efficiently identifying abstractions and so provide a basis for the automation of this task. We illustrate the theory and techniques developed by investigating the identification of abstractions for two published MVN models of the lysis-lysogeny switch in the bacteriophage lambda.


title: Computational Modeling for the Activation Cycle of G-proteins by G-protein-coupled Receptors

abstract: In this paper, we survey five different computational modeling methods. For comparison, we use the activation cycle of G-proteins that regulate cellular signaling events downstream of G-protein-coupled receptors (GPCRs) as a driving example. Starting from an existing Ordinary Differential Equations (ODEs) model, we implement the G-protein cycle in the stochastic Pi-calculus using SPiM, as Petri-nets using Cell Illustrator, in the Kappa Language using Cellucidate, and in Bio-PEPA using the Bio-PEPA eclipse plug in. We also provide a high-level notation to abstract away from communication primitives that may be unfamiliar to the average biologist, and we show how to translate high-level programs into stochastic Pi-calculus processes and chemical reactions.


title: Aspects of multiscale modelling in a process algebra for biological systems

abstract: We propose a variant of the CCS process algebra with new features aiming at allowing multiscale modelling of biological systems. In the usual semantics of process algebras for modelling biological systems actions are instantaneous. When different scale levels of biological systems are considered in a single model, one should take into account that actions at a level may take much more time than actions at a lower level. Moreover, it might happen that while a component is involved in one long lasting high level action, it is involved also in several faster lower level actions. Hence, we propose a process algebra with operations and with a semantics aimed at dealing with these aspects of multiscale modelling. We study behavioural equivalences for such an algebra and give some examples.


title: Multiscale Bone Remodelling with Spatial P Systems

abstract: Many biological phenomena are inherently multiscale, i.e. they are characterized by interactions involving different spatial and temporal scales simultaneously. Though several approaches have been proposed to provide "multilayer" models, only Complex Automata, derived from Cellular Automata, naturally embed spatial information and realize multiscaling with well-established inter-scale integration schemas. Spatial P systems, a variant of P systems in which a more geometric concept of space has been added, have several characteristics in common with Cellular Automata. We propose such a formalism as a basis to rephrase the Complex Automata multiscaling approach and, in this perspective, provide a 2-scale Spatial P system describing bone remodelling. The proposed model not only results to be highly faithful and expressive in a multiscale scenario, but also highlights the need of a deep and formal expressiveness study involving Complex Automata, Spatial P systems and other promising multiscale approaches, such as our shape-based one already resulted to be highly faithful.


title: Modeling biological systems with delays in Bio-PEPA

abstract: Delays in biological systems may be used to model events for which the underlying dynamics cannot be precisely observed, or to provide abstraction of some behavior of the system resulting more compact models. In this paper we enrich the stochastic process algebra Bio-PEPA, with the possibility of assigning delays to actions, yielding a new non-Markovian process algebra: Bio-PEPAd. This is a conservative extension meaning that the original syntax of Bio-PEPA is retained and the delay specification which can now be associated with actions may be added to existing Bio-PEPA models. The semantics of the firing of the actions with delays is the delay-as-duration approach, earlier presented in papers on the stochastic simulation of biological systems with delays. These semantics of the algebra are given in the Starting-Terminating style, meaning that the state and the completion of an action are observed as two separate events, as required by delays. Furthermore we outline how to perform stochastic simulation of Bio-PEPAd systems and how to automatically translate a Bio-PEPAd system into a set of Delay Differential Equations, the deterministic framework for modeling of biological systems with delays. We end the paper with two example models of biological systems with delays to illustrate the approach.


title: Hybrid Calculus of Wrapped Compartments

abstract: The modelling and analysis of biological systems has deep roots in Mathematics, specifically in the field of ordinary differential equations (ODEs). Alternative approaches based on formal calculi, often derived from process algebras or term rewriting systems, provide a quite complementary way to analyze the behaviour of biological systems. These calculi allow to cope in a natural way with notions like compartments and membranes, which are not easy (sometimes impossible) to handle with purely numerical approaches, and are often based on stochastic simulation methods. Recently, it has also become evident that stochastic effects in regulatory networks play a crucial role in the analysis of such systems. Actually, in many situations it is necessary to use stochastic models. For example when the system to be described is based on the interaction of few molecules, when we are at the presence of a chemical instability, or when we want to simulate the functioning of a pool of entities whose compartmentalised structure evolves dynamically. In contrast, stable metabolic networks, involving a large number of reagents, for which the computational cost of a stochastic simulation becomes an insurmountable obstacle, are efficiently modelled with ODEs. In this paper we define a hybrid simulation method, combining the stochastic approach with ODEs, for systems described in CWC, a calculus on which we can express the compartmentalisation of a biological system whose evolution is defined by a set of rewrite rules.


title: Edge- and Node-Disjoint Paths in P Systems

abstract: In this paper, we continue our development of algorithms used for topological network discovery. We present native P system versions of two fundamental problems in graph theory: finding the maximum number of edge- and node-disjoint paths between a source node and target node. We start from the standard depth-first-search maximum flow algorithms, but our approach is totally distributed, when initially no structural information is available and each P system cell has to even learn its immediate neighbors. For the node-disjoint version, our P system rules are designed to enforce node weight capacities (of one), in addition to edge capacities (of one), which are not readily available in the standard network flow algorithms.


title: Lumpability Abstractions of Rule-based Systems

abstract: The induction of a signaling pathway is characterized by transient complex formation and mutual posttranslational modification of proteins. To faithfully capture this combinatorial process in a mathematical model is an important challenge in systems biology. Exploiting the limited context on which most binding and modification events are conditioned, attempts have been made to reduce the combinatorial complexity by quotienting the reachable set of molecular species, into species aggregates while preserving the deterministic semantics of the thermodynamic limit. Recently we proposed a quotienting that also preserves the stochastic semantics and that is complete in the sense that the semantics of individual species can be recovered from the aggregate semantics. In this paper we prove that this quotienting yields a sufficient condition for weak lumpability and that it gives rise to a backward Markov bisimulation between the original and aggregated transition system. We illustrate the framework on a case study of the EGF/insulin receptor crosstalk.


title: The baroclinic instability in the context of layered accretion. Self-sustained vortices and their magnetic stability in local compressible unstratified models of protoplanetary disks

abstract: Turbulence and angular momentum transport in accretion disks remains a topic of debate. With the realization that dead zones are robust features of protoplanetary disks, the search for hydrodynamical sources of turbulence continues. A possible source is the baroclinic instability (BI), which has been shown to exist in unmagnetized non-barotropic disks. We present shearing box simulations of baroclinicly unstable, magnetized, 3D disks, in order to assess the interplay between the BI and other instabilities, namely the magneto-rotational instability (MRI) and the magneto-elliptical instability. We find that the vortices generated and sustained by the baroclinic instability in the purely hydrodynamical regime do not survive when magnetic fields are included. The MRI by far supersedes the BI in growth rate and strength at saturation. The resulting turbulence is virtually identical to an MRI-only scenario. We measured the intrinsic vorticity profile of the vortex, finding little radial variation in the vortex core. Nevertheless, the core is disrupted by an MHD instability, which we identify with the magneto-elliptic instability. This instability has nearly the same range of unstable wavelengths as the MRI, but has higher growth rates. In fact, we identify the MRI as a limiting case of the magneto-elliptic instability, when the vortex aspect ratio tends to infinity (pure shear flow). We conclude that vortex excitation and self-sustenance by the baroclinic instability in protoplanetary disks is viable only in low ionization, i.e., the dead zone. Our results are thus in accordance with the layered accretion paradigm. A baroclinicly unstable dead zone should be characterized by the presence of large-scale vortices whose cores are elliptically unstable, yet sustained by the baroclinic feedback. As magnetic fields destroy the vortices and the MRI outweighs the BI, the active layers are unmodified.


title: Qualitative modelling and analysis of regulations in multi-cellular systems using Petri nets and topological collections

abstract: In this paper, we aim at modelling and analyzing the regulation processes in multi-cellular biological systems, in particular tissues. The modelling framework is based on interconnected logical regulatory networks a la Rene Thomas equipped with information about their spatial relationships. The semantics of such models is expressed through colored Petri nets to implement regulation rules, combined with topological collections to implement the spatial information. Some constraints are put on the the representation of spatial information in order to preserve the possibility of an enumerative and exhaustive state space exploration. This paper presents the modelling framework, its semantics, as well as a prototype implementation that allowed preliminary experimentation on some applications.


title: Bordered Floer homology and the spectral sequence of a branched double cover I

abstract: Given a link in the three-sphere, Z. Szab\'o and the second author constructed a spectral sequence starting at the Khovanov homology of the link and converging to the Heegaard Floer homology of its branched double-cover. The aim of this paper and its sequel is to explicitly calculate this spectral sequence, using bordered Floer homology. There are two primary ingredients in this computation: an explicit calculation of filtered bimodules associated to Dehn twists and a pairing theorem for polygons. In this paper we give the first ingredient, and so obtain a combinatorial spectral sequence from Khovanov homology to Heegaard Floer homology; in the sequel we show that this spectral sequence agrees with the previously known one.


title: A novel method to extract dark matter parameters from neutrino telescope data

abstract: Recently it has been shown that when the Dark Matter (DM) particles captured in the Sun directly annihilate into neutrino pairs, the oscillatory terms in the oscillation probability do not average to zero and can lead to a seasonal variation as the distance between the Sun and Earth changes in time. In this paper, we explore this feature as a novel method to extract information on the properties of dark matter. We show that by studying the variation of the flux over a few months, it would in principle be possible to derive the DM mass as well as new information on the flavor structure of the DM annihilation modes. In addition to analytic analysis, we present the results of our numerical calculations that take into account scattering and regeneration of neutrinos traversing the Sun.


title: Geometric phase for nonlinear coherent and squeezed state

abstract: The geometric phases for standard coherent states which are widely used in quantum optics have attracted a large amount of attention. Nevertheless, few physicists consider about the counterparts of non-linear coherent states, which are useful in the description of the motion of a trapped ion. In this paper, the non-unitary and non-cyclic geometric phases for two nonlinear coherent and one squeezed states are formulated respectively. Moreover, some of their common properties are discussed respectively, such as gauge invariance, non-locality and non-linear effects. The non-linear functions have dramatic impacts on the evolution of the corresponding geometric phases. They speed the evolution up or down. So this property may have application in controlling or measuring geometric phase. For the squeezed case, when the squeezed parameter r -> \infinity, the limiting value of the geometric phase is also determined by non-linear function at a given time and angular velocity. In addition, the geometric phases for standard coherent and squeezed states are obtained under a particular condition. When the time evolution undergoes a period, their corresponding cyclic geometric phases are achieved as well. And the distinction between the geometric phases of the two coherent states maybe regarded as a geometric criterion.


title: A New Email Retrieval Ranking Approach

abstract: Email Retrieval task has recently taken much attention to help the user retrieve the email(s) related to the submitted query. Up to our knowledge, existing email retrieval ranking approaches sort the retrieved emails based on some heuristic rules, which are either search clues or some predefined user criteria rooted in email fields. Unfortunately, the user usually does not know the effective rule that acquires best ranking related to his query. This paper presents a new email retrieval ranking approach to tackle this problem. It ranks the retrieved emails based on a scoring function that depends on crucial email fields, namely subject, content, and sender. The paper also proposes an architecture to allow every user in a network/group of users to be able, if permissible, to know the most important network senders who are interested in his submitted query words. The experimental evaluation on Enron corpus prove that our approach outperforms known email retrieval ranking approaches.


title: Measurements of vector magnetic field using multiple electromagnetically induced transparency resonances in Rb vapor

abstract: We study intensity distribution of transmission resonances under the conditions of electromagnetically induced transparency (EIT) with lin||lin configuration in 87Rb vapor. We demonstrate that for various orientations of the external magnetic field all seven resolved EIT resonances exhibit maxima or minima at certain directions of the linear polarization of the laser field. This effect can be used for development of a high-precision vector EIT magnetometer.


title: On the weighted forward reduced Entropy of Ricci flow

abstract: In this paper, we first introduce the weighted forward reduced volume of Ricci flow. The weighted forward reduced volume, which related to expanders of Ricci flow, is well-defined on noncompact manifolds and monotone non-increasing under Ricci flow. Moreover, we show that, just the same as the Perelman's reduced volume, the weighted reduced volume entropy has the value $(4\pi)^{\frac{n}{2}}$ if and only if the Ricci flow is the trivial flow on flat Euclidean space.


title: Spontaneously axisymmetry breaking phase in a binary mixture of spinor Bose-Einstein condensates

abstract: We study the ground state phases for a mixture of two atomic spin-1 Bose-Einstein condensates (BECs) in the presence of a weak magnetic (B-) field. The ground state is found to contain a broken-axisymmetry (BA) phase due to competitions among intra- and inter-species spin exchange interactions and the linear Zeeman shifts. This is in contrast to the case of a single species spin- 1 condensate, where the axisymmetry breaking results from competitions among the linear and quadratic Zeeman shifts and the intra-species ferromagnetic interaction. All other remaining ground state phases for the mixture are found to preserve axisymmetry. We further elaborate on the ground state phase diagram and calculate their Bogoliubov excitation spectra. For the BA phase, there exist three Goldstone modes attempting to restore the broken U(1) and SO(2) symmetries.


title: A Very Fast Algorithm for Matrix Factorization

abstract: We present a very fast algorithm for general matrix factorization of a data matrix for use in the statistical analysis of high-dimensional data via latent factors. Such data are prevalent across many application areas and generate an ever-increasing demand for methods of dimension reduction in order to undertake the statistical analysis of interest. Our algorithm uses a gradient-based approach which can be used with an arbitrary loss function provided the latter is differentiable. The speed and effectiveness of our algorithm for dimension reduction is demonstrated in the context of supervised classification of some real high-dimensional data sets from the bioinformatics literature.


title: Level Shifter Design for Low Power Applications

abstract: With scaling of Vt sub-threshold leakage power is increasing and expected to become significant part of total power consumption In present work three new configurations of level shifters for low power application in 0.35{\mu}m technology have been presented. The proposed circuits utilize the merits of stacking technique with smaller leakage current and reduction in leakage power. Conventional level shifter has been improved by addition of three NMOS transistors, which shows total power consumption of 402.2264pW as compared to 0.49833nW with existing circuit. Single supply level shifter has been modified with addition of two NMOS transistors that gives total power consumption of 108.641pW as compared to 31.06nW. Another circuit, contention mitigated level shifter (CMLS) with three additional transistors shows total power consumption of 396.75pW as compared to 0.4937354nW. Three proposed circuit's shows better performance in terms of power consumption with a little conciliation in delay. Output level of 3.3V has been obtained with input pulse of 1.6V for all proposed circuits.


title: Equation of state for dark energy in $f(T)$ gravity

abstract: We study the cosmological evolutions of the equation of state for dark energy $w_{\mathrm{DE}}$ in the exponential and logarithmic as well as their combination $f(T)$ theories. We show that the crossing of the phantom divide line of $w_{\mathrm{DE}} = -1$ can be realized in the combined $f(T)$ theory even though it cannot be in the exponential or logarithmic $f(T)$ theory. In particular, the crossing is from $w_{\mathrm{DE}} > -1$ to $w_{\mathrm{DE}} < -1$, in the opposite manner from $f(R)$ gravity models. We also demonstrate that this feature is favored by the recent observational data.


title: Generic Representation Theory of the Additive and Heisenberg Groups

abstract: In this paper we give an intimate connection between the characteristic zero representation theories of the Additive and Heisenberg groups, and their characteristic p >0 theories when p is much larger than the dimension a representation. In particular, if p >> dimension, then all characteristic p representations for these groups can be factored into commuting products of representations, with each factor arising from a representation of the Lie algebra of the group, one for each of the representation's Frobenius layers. In this sense, for a fixed dimension and large enough p, all representations for these groups look generically like representations for direct powers of themselves over a field of characteristic zero.


title: Theory of cooperation in a micro-organismal snow-drift game

abstract: We present a mean field model for the phase diagram of a community of micro-organisms, interacting through their metabolism so that they are, in effect, engaging in a cooperative social game. We show that as a function of the concentration of the nutrients glucose and histidine, the community undergoes a phase transition separating a state in which one strain is dominant to a state which is characterized by coexisting populations. Our results are in good agreement with recent experimental results, correctly predicting quantitative trends and the phase diagram.


title: Supernova dust for the extinction law in a young infrared galaxy at z = 1

abstract: We apply the supernova(SN) extinction curves to reproduce the observed properties of SST J1604+4304 which is a young infrared (IR) galaxy at z = 1. The SN extinction curves used in this work were obtained from models of unmixed ejecta of type II supernovae(SNe II) for the Salpeter initial mass function (IMF) with a mass range from 8 to 30 M_sun or 8 to 40 M_sun. The effect of dust distributions on the attenuation of starlight is investigated by performing the chi-square fitting method against various dust distributions. These are the commonly used uniform dust screen, the clumpy dust screen, and the internal dust geometry. We add to these geometries three scattering properties, namely, no-scattering, isotropic scattering, and forward-only scattering. Judging from the chi-square values, we find that the uniform screen models with any scattering property provide good approximations to the real dust geometry. Internal dust is inefficient to attenuate starlight and thus cannot be the dominant source of the extinction. We show that the SN extinction curves reproduce the data of SST J1604+4304 comparable to or better than the Calzetti extinction curve. The Milky Way extinction curve is not in satisfactory agreement with the data unless several dusty clumps are in the line of sight. This trend may be explained by the abundance of SN-origin dust in these galaxies; SN dust is the most abundant in the young IR galaxy at z = 1, abundant in local starbursts, and less abundant in the Galaxy. If dust in SST J1604+4304 is dominated by SN dust, the dust production rate is about 0.1 M_sun per SN.


title: Vector Magnetic Fields and Current Helicities in Coronal Holes and Quiet Regions

abstract: In the solar photosphere, many properties of coronal holes (CHs) are not known, especially vector magnetic fields. Using observations from \emph{Hinode}, we investigate vector magnetic fields, current densities and current helicities in two CHs and compare them with two normal quiet regions (QRs) for the first time. We find that, in the CHs and QRs, the areas where large current helicities are located are mainly co-spatial with strong vertical and horizontal field elements both in shape and location. In the CHs, horizontal magnetic fields, inclination angles, current densities and current helicities are larger than those in the QRs. The mean vertical current density and current helicity, averaged over all the observed areas including the CHs and QRs, are approximately 0.008 A m$^{-2}$ and 0.005 G$^{2}$ m$^{-1}$, respectively. The mean current density in magnetic flux concentrations where the vertical fields are stronger than 100 G is as large as 0.012 $\pm$ 0.001 A m$^{-2}$, consistent with that in the flare productive active regions. Our results imply that the magnetic fields, especially the strong fields, both in the CHs and QRs are nonpotential.


title: Convergence of vector bundles with metrics of Sasaki-type

abstract: If a sequence of Riemannian manifolds, $X_i$, converges in the pointed Gromov-Hausdorff sense to a limit space, $X_\infty$, and if $E_i$ are vector bundles over $X_i$ endowed with metrics of Sasaki-type with a uniform upper bound on rank, then a subsequence of the $E_i$ converges in the pointed Gromov-Hausdorff sense to a metric space, $E_\infty$. The projection maps $\pi_i$ converge to a limit submetry $\pi_\infty$ and the fibers converge to its fibers; the latter may no longer be vector spaces but are homeomorphic to $\R^k/G$, where $G$ is a closed subgroup of $O(k)$ ---called the {\em wane group}--- that depends on the basepoint and that is defined using the holonomy groups on the vector bundles. The norms $\mu_i=\|\cdot\|_i$ converges to a map $\mu_{\infty}$ compatible with the re-scaling in $\R^k/G$ and the $\R$-action on $E_i$ converges to an $\R-$action on $E_{\infty}$ compatible with the limiting norm. In the special case when the sequence of vector bundles has a uniform lower bound on holonomy radius (as in a sequence of collapsing flat tori to a circle), the limit fibers are vector spaces. Under the opposite extreme, e.g. when a single compact $n$-dimensional manifold is re-scaled to a point, the limit fiber is $\R^n/H$ where $H$ is the closure of the holonomy group of the compact manifold considered. An appropriate notion of parallelism is given to the limiting spaces by considering curves whose length is unchanged under the projection. The class of such curves is invariant under the $\R$-action and each such curve preserves norms. The existence of parallel translation along rectifiable curves with arbitrary initial conditions is also exhibited. Uniqueness is not true in general, but a necessary condition is given in terms of the aforementioned wane groups $G$.