58 research outputs found
Ion Pair Potentials-of-Mean-Force in Water
Recent molecular simulation and integral equation results alkali-halide ion
pair potentials-of-mean-force in water are discussed. Dielectric model
calculations are implemented to check that these models produce that
characteristic structure of contact and solvent-separated minima for oppositely
charged ions in water under physiological thermodynamic conditions. Comparison
of the dielectric model results with the most current molecular level
information indicates that the dielectric model does not, however, provide an
accurate description of these potentials-of-mean-force. We note that linear
dielectric models correspond to modelistic implementations of second-order
thermodynamic perturbation theory for the excess chemical potential of a
distinguished solute molecule. Therefore, the molecular theory corresponding to
the dielectric models is second-order thermodynamic perturbation theory for
that excess chemical potential. The second-order, or fluctuation, term raises a
technical computational issue of treatment of long-ranged interactions similar
to the one which arises in calculation of the dielectric constant of the
solvent. It is contended that the most important step for further development
of dielectric models would be a separate assessment of the first-order
perturbative term (equivalently the {\it potential at zero charge} ) which
vanishes in the dielectric models but is generally nonzero. Parameterization of
radii and molecular volumes should then be based of the second-order
perturbative term alone. Illustrative initial calculations are presented and
discussed.Comment: 37 pages and 8 figures. LA-UR-93-420
Ion Sizes and Finite-Size Corrections for Ionic-Solvation Free Energies
Free energies of ionic solvation calculated from computer simulations exhibit
a strong system size dependence. We perform a finite-size analysis based on a
dielectric-continuum model with periodic boundary conditions. That analysis
results in an estimate of the Born ion size. Remarkably, the finite-size
correction applies to systems with only eight water molecules hydrating a
sodium ion and results in an estimate of the Born radius of sodium that agrees
with the experimental value.Comment: 2 EPS figure
Waveform-based simulated annealing of crosshole transmission data: a semi-global method for estimating seismic anisotropy
We successfully apply the semi-global inverse method of simulated annealing to determine the best-fitting 1-D anisotropy model for use in acoustic frequency domain waveform tomography. Our forward problem is based on a numerical solution of the frequency domain acoustic wave equation, and we minimize wavefield phase residuals through random perturbations to a 1-D vertically varying anisotropy profile. Both real and synthetic examples are presented in order to demonstrate and validate the approach. For the real data example, we processed and inverted a cross-borehole data set acquired by Vale Technology Development (Canada) Ltd. in the Eastern Deeps deposit, located in Voisey's Bay, Labrador, Canada. The inversion workflow comprises the full suite of acquisition, data processing, starting model building through traveltime tomography, simulated annealing and finally waveform tomography. Waveform tomography is a high resolution method that requires an accurate starting model. A cycle-skipping issue observed in our initial starting model was hypothesized to be due to an erroneous anisotropy model from traveltime tomography. This motivated the use of simulated annealing as a semi-global method for anisotropy estimation. We initially tested the simulated annealing approach on a synthetic data set based on the Voisey's Bay environment; these tests were successful and led to the application of the simulated annealing approach to the real data set. Similar behaviour was observed in the anisotropy models obtained through traveltime tomography in both the real and synthetic data sets, where simulated annealing produced an anisotropy model which solved the cycle-skipping issue. In the real data example, simulated annealing led to a final model that compares well with the velocities independently estimated from borehole logs. By comparing the calculated ray paths and wave paths, we attributed the failure of anisotropic traveltime tomography to the breakdown of the ray-theoretical approximation in the vicinity of strong velocity discontinuitie
Climatic Cycles Recorded in Glacially Influenced Rhythmites of the Gowganda Formation, Huronian Supergroup
The Gowganda Formation of the 2.45–2.2 Ga Huronian Supergroup contains glacially-induced, varve-like rhythmites that potentially preserve a detailed record of climatic conditions during the Paleoproterozoic Era. Four rhythmic couplet thickness records were measured at two outcrops near Wharncliffe, Ontario for the purpose of time-series analysis. The couplets, which range from 1 to 32 mm thick, are composed of alternating layers of siltstone and claystone. Time-series analysis of the couplet thickness records using the MTM Toolkit of Mann and Lees (1996) consistently revealed periodicities in the range of 2.2–2.9 couplets per cycle, which is consistent with climatic cycles such as the quasi-biennial oscillation (QBO) and the El Niño Southern Oscillation (ENSO) observed in modern times. This periodicity suggests that the rhythmic couplets represent annual deposits (i.e. varves). Evidence for the presence of cycles at 3.0– 4.9 couplets, 6.6–6.9 couplets, 8.8–9.2 couplets, 22.8 couplets, and 30.1–31.0 couplets were also observed in some couplet thickness records; however, the presence of these longer term cycles was inconsistent from site to site
Origin of entropy convergence in hydrophobic hydration and protein folding
An information theory model is used to construct a molecular explanation why
hydrophobic solvation entropies measured in calorimetry of protein unfolding
converge at a common temperature. The entropy convergence follows from the weak
temperature dependence of occupancy fluctuations for molecular-scale volumes in
water. The macroscopic expression of the contrasting entropic behavior between
water and common organic solvents is the relative temperature insensitivity of
the water isothermal compressibility. The information theory model provides a
quantitative description of small molecule hydration and predicts a negative
entropy at convergence. Interpretations of entropic contributions to protein
folding should account for this result.Comment: Phys. Rev. Letts. (in press 1996), 3 pages, 3 figure
Upper limits on the strength of periodic gravitational waves from PSR J1939+2134
The first science run of the LIGO and GEO gravitational wave detectors
presented the opportunity to test methods of searching for gravitational waves
from known pulsars. Here we present new direct upper limits on the strength of
waves from the pulsar PSR J1939+2134 using two independent analysis methods,
one in the frequency domain using frequentist statistics and one in the time
domain using Bayesian inference. Both methods show that the strain amplitude at
Earth from this pulsar is less than a few times .Comment: 7 pages, 1 figure, to appear in the Proceedings of the 5th Edoardo
Amaldi Conference on Gravitational Waves, Tirrenia, Pisa, Italy, 6-11 July
200
Improving the sensitivity to gravitational-wave sources by modifying the input-output optics of advanced interferometers
We study frequency dependent (FD) input-output schemes for signal-recycling
interferometers, the baseline design of Advanced LIGO and the current
configuration of GEO 600. Complementary to a recent proposal by Harms et al. to
use FD input squeezing and ordinary homodyne detection, we explore a scheme
which uses ordinary squeezed vacuum, but FD readout. Both schemes, which are
sub-optimal among all possible input-output schemes, provide a global noise
suppression by the power squeeze factor, while being realizable by using
detuned Fabry-Perot cavities as input/output filters. At high frequencies, the
two schemes are shown to be equivalent, while at low frequencies our scheme
gives better performance than that of Harms et al., and is nearly fully
optimal. We then study the sensitivity improvement achievable by these schemes
in Advanced LIGO era (with 30-m filter cavities and current estimates of
filter-mirror losses and thermal noise), for neutron star binary inspirals, and
for narrowband GW sources such as low-mass X-ray binaries and known radio
pulsars. Optical losses are shown to be a major obstacle for the actual
implementation of these techniques in Advanced LIGO. On time scales of
third-generation interferometers, like EURO/LIGO-III (~2012), with
kilometer-scale filter cavities, a signal-recycling interferometer with the FD
readout scheme explored in this paper can have performances comparable to
existing proposals. [abridged]Comment: Figs. 9 and 12 corrected; Appendix added for narrowband data analysi
- …