An assessment of key model parametric uncertainties in projections of Greenland Ice Sheet behavior

Lack of knowledge about the values of ice sheet model input parameters introduces substantial uncertainty into projections of Greenland Ice Sheet contributions to future sea level rise. Computer models of ice sheet behavior provide one of several means of estimating future sea level rise due to mass loss from ice sheets. Such models have many input parameters whose values are not well known. Recent studies have investigated the effects of these parameters on model output, but the range of potential future sea level increases due to model parametric uncertainty has not been characterized. Here, we demonstrate that this range is large, using a 100-member perturbed-physics ensemble with the SICOPOLIS ice sheet model. Each model run is spun up over 125 000 yr using geological forcings and subsequently driven into the future using an asymptotically increasing air temperature anomaly curve. All modeled ice sheets lose mass after 2005 AD. Parameters controlling surface melt dominate the model response to temperature change. After culling the ensemble to include only members that give reasonable ice volumes in 2005 AD, the range of projected sea level rise values in 2100 AD is ~40 % or more of the median. Data on past ice sheet behavior can help reduce this uncertainty, but none of our ensemble members produces a reasonable ice volume change during the mid-Holocene, relative to the present. This problem suggests that the model's exponential relation between temperature and precipitation does not hold during the Holocene, or that the central-Greenland temperature forcing curve used to drive the model is not representative of conditions around the ice margin at this time (among other possibilities). Our simulations also lack certain observed physical processes that may tend to enhance the real ice sheet's response. Regardless, this work has implications for other studies that use ice sheet models to project or hindcast the behavior of the Greenland Ice Sheet

Spin-waves in the J1a−J1b−J2J_{1a}-J_{1b}-J_{2} orthorombic square-lattice Heisenberg models: Application to the iron pnictide materials

Motivated by the observation of spatially anisotropic exchange constants in the iron pnictide materials, we study the spin-wave spectra of the $J_{1a}-J_{1b}-J_{2}$ Heisenberg models on a square-lattice with nearest neighbor exchange $J_{1a}$ along x and $J_{1b}$ along y axis and a second neighbor exchange $J_2$. We focus on the regime, where the spins order at ($\pi,0$), and compute the spectra by systematic expansions around the Ising limit. We study both spin-half and spin-one Heisenberg models as well as a range of parameters to cover various cases proposed for the iron pnictide materials. The low-energy spectra have anisotropic spin-wave velocities and are renormalized with respect to linear spin-wave theory by up to 20 percent, depending on parameters. Extreme anisotropy, consisting of a ferromagnetic $J_{1b}=- |J_F|$, is best distinguished from a weak anisotropy ($J_{1a}\approx J_{1b}=J_1$, $J_2>J_1/2$) by the nature of the spin-waves near the wavevectors ($0,\pi$) or ($\pi,\pi$). The reported spectra for the pnictide material CaFe$_2$As$_2$ clearly imply such an extreme anisotropy.Comment: 6 pages, 10 figure

A New Sub-Period-Minimum Cataclysmic Variable With Partial Hydrogen Depletion And Evidence Of Spiral Disk Structure

We present time-resolved spectroscopy and photometry of CSS 120422:111127+571239 (=SBS 1108+574), a recently discovered SU UMa-type dwarf nova whose 55 minute orbital period is well below the cataclysmic variable (CV) period minimum of similar to 78 minutes. In contrast with most other known CVs, its spectrum features He I emission of comparable strength to the Balmer lines, implying a hydrogen abundance less than 0.1 of long-period CVs-but still at least 10 times higher than that in AM CVn stars. Together, the short orbital period and remarkable helium-to-hydrogen ratio suggest that mass transfer in CSS 120422 began near the end of the donor star's main-sequence lifetime, meaning that this CV is a strong candidate progenitor of an AM CVn system as described by Podsiadlowski et al. Moreover, a Doppler tomogram of the Ha line reveals two distinct regions of enhanced emission. While one is the result of the stream-disk impact, the other is probably attributable to spiral disk structure generated when material in the outer disk achieves a 2:1 orbital resonance with respect to the donor.NSF AST-1211196, AST-9987045Department of Physics at the University of Notre DameNSF Telescope System Instrumentation Program (TSIP)Ohio Board of RegentsOhio State University Office of ResearchAstronom

Is the Yb2Ti2O7 pyrochlore a quantum spin ice?

We use numerical linked cluster (NLC) expansions to compute the specific heat, C(T), and entropy, S(T), of a quantum spin ice model of Yb2Ti2O7 using anisotropic exchange interactions recently determined from inelastic neutron scattering measurements and find good agreement with experimental calorimetric data. In the perturbative weak quantum regime, this model has a ferrimagnetic ordered ground state, with two peaks in C(T): a Schottky anomaly signalling the paramagnetic to spin ice crossover followed at lower temperature by a sharp peak accompanying a first order phase transition to the ferrimagnetic state. We suggest that the two C(T) features observed in Yb2Ti2O7 are associated with the same physics. Spin excitations in this regime consist of weakly confined spinon-antispinon pairs. We suggest that conventional ground state with exotic quantum dynamics will prove a prevalent characteristic of many real quantum spin ice materials.Comment: 8 pages (two-column), 9 figure

Mutant enrichment by filtration concentration: a variation for the selection of temperature-conditional heterocaryons

Mutant enrichment by filtration concentration: a variation for the selection of temperature-conditional heterocaryon

Collaborative Delivery with Energy-Constrained Mobile Robots

We consider the problem of collectively delivering some message from a specified source to a designated target location in a graph, using multiple mobile agents. Each agent has a limited energy which constrains the distance it can move. Hence multiple agents need to collaborate to move the message, each agent handing over the message to the next agent to carry it forward. Given the positions of the agents in the graph and their respective budgets, the problem of finding a feasible movement schedule for the agents can be challenging. We consider two variants of the problem: in non-returning delivery, the agents can stop anywhere; whereas in returning delivery, each agent needs to return to its starting location, a variant which has not been studied before. We first provide a polynomial-time algorithm for returning delivery on trees, which is in contrast to the known (weak) NP-hardness of the non-returning version. In addition, we give resource-augmented algorithms for returning delivery in general graphs. Finally, we give tight lower bounds on the required resource augmentation for both variants of the problem. In this sense, our results close the gap left by previous research.Comment: 19 pages. An extended abstract of this paper was published at the 23rd International Colloquium on Structural Information and Communication Complexity 2016, SIROCCO'1

Some Problems in Probabilistic Tomography

Given probability distributions F1 , F2 , . . ., Fk on R and distinct directions θ1, . . ., θk, one may ask whether there is a probability measure μ on R2 such that the marginal of μ in direction θj is Fj, j = 1, . . ., k. For example for k = 3 we ask what the marginal of μ at 45° can be if the x and y marginals are each say standard normal? In probabilistic language, if X and Y are each standard normal with an arbitrary joint distribution, what can the distribution of X + Y or X - Y be? This type of question is familiar to probabilists and is also familiar (except perhaps in that μ is positive) to tomographers, but is difficult to answer in special cases. The set of distributions for Z = X - Y is a convex and compact set, C, which contains the single point mass Z ≡ 0 since X ≡ Y, standard normal, is possible. We show that Z can be 3-valued, Z=0, ±a for any a, each with positive probability, but Z cannot have any (genuine) two-point distribution. Using numerical linear programming we present convincing evidence that Z can be uniform on the interval [-ε, ε] for ε small and give estimates for the largest such ε. The set of all extreme points of C seems impossible to determine explicitly. We also consider the more basic question of finding the extreme measures on the unit square with uniform marginals on both coordinates, and show that not every such measure has a support which has only one point on each horizontal or vertical line, which seems surprising

Large Scale Inhomogeneities from the QCD Phase Transition

We examine the first-order cosmological QCD phase transition for a large class of parameter values, previously considered unlikely. We find that the hadron bubbles can nucleate at very large distance scales, they can grow as detonations as well as deflagrations, and that the phase transition may be completed without reheating to the critical temperature. For a subset of the parameter values studied, the inhomogeneities generated at the QCD phase transition might have a noticeable effect on nucleosynthesis.Comment: 15 LaTeX pages + 6 PostScript figures appended at the end of the file, HU-TFT-94-1