Liquid-Vapor Equilibrium of Multicomponent Cryogenic Systems

Liquid-vapor and solid-vapor equilibria at low to moderate pressures and low temperatures are important in many solar system environments, including the surface and clouds of Titan, the clouds of Uranus and Neptune, and the surfaces of Mars and Triton. The familiar cases of ideal behavior are limiting cases of a general thermodynamic representation for the vapor pressure of each component in a homogeneous multicomponent system. The fundamental connections of laboratory measurements to thermodynamic models are through the Gibbs-Duhem relation and the Gibbs-Helmholtz relation. Using laboratory measurements of the total pressure, temperature, and compositions of the liquid and vapor phases at equilibrium, the values of these parameters can be determined. The resulting model for vapor-liquid equilibrium can then conveniently and accurately be used to calculate pressures, compositions, condensation altitudes, and their dependencies on changing climatic conditions. A specific system being investigated is CH4-C2H6-N2, at conditions relevant to Titan's surface and atmosphere. Discussed are: the modeling of existing data on CH4-N2, with applications to the composition of Titan's condensate clouds; some new measurements on the CH4-C2H6 binary, using a high-precision static/volumetric system, and on the C2H6-N2 binary, using the volumetric system and a sensitive cryogenic flow calorimeter; and describe a new cryogenic phase-equilibrium vessel with which we are beginning a detailed, systematic study of the three constituent binaries and the ternary CH4-C2H6-N2 system at temperatures ranging from 80 to 105 K and pressures from 0.1 to 7 bar

Van der Waals loops and the melting transition in two dimensions

Evidence for the existence of van der Waals loops in pressure p versus volume v plots has for some time supported the belief that melting in two dimensions is a first order phase transition. We report rather accurate equilibrium p(v) curves for systems of hard disks obtained from long Monte Carlo simulations. These curves, obtained in the constant volume ensemble, using periodic boundary conditions, exhibit well defined van der Waals loops. We illustrate their existence for finite systems that are known to undergo a continuous transition in the thermodynamic limit. To this end, we obtain magnetization m versus applied field curves from Monte Carlo simulations of the 2D Ising model, in the constant m ensemble, at the critical point. Whether van der Waals loops for disk systems behave in the thermodynamic limit as they do for the 2D Ising model at the critical point cannot be ruled out. Thus, the often made claim that melting in 2D is a first order phase transition, based on the evidence that van der Waals loops exist, is not sound.Comment: 10 pages, 6 Postscript figures (submitted to Phys.Rev.E). For related work, see http://pipe.unizar.es/~jf

Liquid-Solid Phase Transition of the System with Two particles in a Rectangular Box

We study the statistical properties of two hard spheres in a two dimensional rectangular box. In this system, the relation like Van der Waals equation loop is obtained between the width of the box and the pressure working on side walls. The auto-correlation function of each particle's position is calculated numerically. By this calculation near the critical width, the time at which the correlation become zero gets longer according to the increase of the height of the box. Moreover, fast and slow relaxation processes like $\alpha$ and $\beta$ relaxations observed in supper cooled liquid are observed when the height of the box is sufficiently large. These relaxation processes are discussed with the probability distribution of relative position of two particles.Comment: 6 figure

THE OAK ORCHARD SOIL WATER ASSESSMENT TOOL A decision support system for watershed management Part 1: Calibration and Validation

A hydrologic model (SWAT) was developed and calibrated for the Oak Orchard watershed to evaluate sources and sinks of sediment and nutrients. The model included the most important anthropogenic features that impacted water flow and nonpoint source pollution in the watershed. These features included reservoirs at the Iroquois National Wildlife Refuge, Waterport and Medina; point sources such as the Erie Canal, US Gypsum, Allen Canning, wastewater treatment plants at Medina, Oakfield and Elba, and tiledrains at the mucklands, an intensely farmed area that was drained to combat malaria in the 19th century. The model included point sources for every subbasin so that the effects of future point sources can be evaluated. The model was calibrated for waterflow and sediment using observed loading data collected by Makarewicz and Lewis (2000, 2009). To achieve the proper water balance observed at the watershed, seasonal inputs of water had to be added from the Erie Canal and the Onondaga escarpment. This water came from outside of the watershed. The resulting calibration had a Nash-Sutcliffe (NS) prediction efficiency of 0.81 for the calibration period (1997-1999). The total cumulative sediment loading was within 2%, of observed and the monthly sediment loads fell within the uncertainty of the observed data (NS=0.31). Cumulative total phosphorous loads were within 2% of observed and the NS prediction efficiency was 0.91. The model validated very poorly in the 2008 time period primarily because of inaccurate precipitation data and incorrect groundwater fluxes from the escarpment. Further research needs to evaluate the timing and amount of groundwater flow from the escarpment because it has a significant impact on monthly flows in this watershed. It is likely that other watersheds that are nestled against the Onondaga escarpment are impacted by spring flows from this geologic feature

Elastic moduli, dislocation core energy and melting of hard disks in two dimensions

Elastic moduli and dislocation core energy of the triangular solid of hard disks of diameter $\sigma$ are obtained in the limit of vanishing dislocation- antidislocation pair density, from Monte Carlo simulations which incorporates a constraint, namely that all moves altering the local connectivity away from that of the ideal triangular lattice are rejected. In this limit, we show that the solid is stable against all other fluctuations at least upto densities as low as $\rho \sigma^2 = 0.88$. Our system does not show any phase transition so diverging correlation lengths leading to finite size effects and slow relaxations do not exist. The dislocation pair formation probability is estimated from the fraction of moves rejected due to the constraint which yields, in turn, the core energy E_c and the (bare) dislocation fugacity y. Using these quantities, we check the relative validity of first order and Kosterlitz-Thouless-Halperin-Nelson-Young (KTHNY) melting scenarios and obtain numerical estimates of the typical expected transition densities and pressures. We conclude that a KTHNY transition from the solid to a hexatic phase preempts the solid to liquid first order transition in this system albeit by a very small margin, easily masked by crossover effects in unconstrained ``brute-force'' simulations with small number of particles.Comment: 17 pages, 8 figure

Arecibo pulsar survey using ALPHA: Probing radio pulsar intermittency and transients

We present radio transient search algorithms, results, and statistics from the ongoing Arecibo Pulsar ALFA (PALFA) survey of the Galactic plane. We have discovered seven objects through a search for isolated dispersed pulses. All of these objects are Galactic and have measured periods between 0.4 and 4.7 s. One of the new discoveries has a duty cycle of 0.01%, smaller than that of any other radio pulsar. We discuss the impact of selection effects on the detectability and classification of intermittent sources, and compare the efficiencies of periodicity and single-pulse (SP) searches for various pulsar classes. For some cases we find that the apparent intermittency is likely to be caused by off-axis detection or a short time window that selects only a few bright pulses and favors detection with our SP algorithm. In other cases, the intermittency appears to be intrinsic to the source. No transients were found with DMs large enough to require that they originate from sources outside our Galaxy. Accounting for the on-axis gain of the ALFA system, as well as the low gain but large solid-angle coverage of far-out sidelobes, we use the results of the survey so far to place limits on the amplitudes and event rates of transients of arbitrary origin. © 2009. The American Astronomical Society. All rights reserved.