Low-lying level structure of 56^{56}Cu and its implications on the rp process

The low-lying energy levels of proton-rich $^{56}$Cu have been extracted using in-beam $\gamma$-ray spectroscopy with the state-of-the-art $\gamma$-ray tracking array GRETINA in conjunction with the S800 spectrograph at the National Superconducting Cyclotron Laboratory at Michigan State University. Excited states in $^{56}$Cu serve as resonances in the $^{55}$Ni(p,$\gamma$)$^{56}$Cu reaction, which is a part of the rp-process in type I x-ray bursts. To resolve existing ambiguities in the reaction Q-value, a more localized IMME mass fit is used resulting in $Q=639\pm82$~keV. We derive the first experimentally-constrained thermonuclear reaction rate for $^{55}$Ni(p,$\gamma$)$^{56}$Cu. We find that, with this new rate, the rp-process may bypass the $^{56}$Ni waiting point via the $^{55}$Ni(p,$\gamma$) reaction for typical x-ray burst conditions with a branching of up to $\sim$40$\%$. We also identify additional nuclear physics uncertainties that need to be addressed before drawing final conclusions about the rp-process reaction flow in the $^{56}$Ni region.Comment: 8 pages, accepted for Phys. Rev.

Instantaneous Liquid Interfaces

We describe and illustrate a simple procedure for identifying a liquid interface from atomic coordinates. In particular, a coarse grained density field is constructed, and the interface is defined as a constant density surface for this coarse grained field. In applications to a molecular dynamics simulation of liquid water, it is shown that this procedure provides instructive and useful pictures of liquid-vapor interfaces and of liquid-protein interfaces.Comment: 15 pages, 4 figure

Gravitational Waves From Known Pulsars: Results From The Initial Detector Era

We present the results of searches for gravitational waves from a large selection of pulsars using data from the most recent science runs (S6, VSR2 and VSR4) of the initial generation of interferometric gravitational wave detectors LIGO (Laser Interferometric Gravitational-wave Observatory) and Virgo. We do not see evidence for gravitational wave emission from any of the targeted sources but produce upper limits on the emission amplitude. We highlight the results from seven young pulsars with large spin-down luminosities. We reach within a factor of five of the canonical spin-down limit for all seven of these, whilst for the Crab and Vela pulsars we further surpass their spin-down limits. We present new or updated limits for 172 other pulsars (including both young and millisecond pulsars). Now that the detectors are undergoing major upgrades, and, for completeness, we bring together all of the most up-to-date results from all pulsars searched for during the operations of the first-generation LIGO, Virgo and GEO600 detectors. This gives a total of 195 pulsars including the most recent results described in this paper.United States National Science FoundationScience and Technology Facilities Council of the United KingdomMax-Planck-SocietyState of Niedersachsen/GermanyAustralian Research CouncilInternational Science Linkages program of the Commonwealth of AustraliaCouncil of Scientific and Industrial Research of IndiaIstituto Nazionale di Fisica Nucleare of ItalySpanish Ministerio de Economia y CompetitividadConselleria d'Economia Hisenda i Innovacio of the Govern de les Illes BalearsNetherlands Organisation for Scientific ResearchPolish Ministry of Science and Higher EducationFOCUS Programme of Foundation for Polish ScienceRoyal SocietyScottish Funding CouncilScottish Universities Physics AllianceNational Aeronautics and Space AdministrationOTKA of HungaryLyon Institute of Origins (LIO)National Research Foundation of KoreaIndustry CanadaProvince of Ontario through the Ministry of Economic Development and InnovationNational Science and Engineering Research Council CanadaCarnegie TrustLeverhulme TrustDavid and Lucile Packard FoundationResearch CorporationAlfred P. Sloan FoundationAstronom

Quantum Monte Carlo simulations of solids

Published versio

Constraining the Neutron Star Compactness: Extraction of the 23Al(p,γ) Reaction Rate for the rp Process

The $^{23}$Al($p,\gamma$)$^{24}$Si reaction is among the most important reactions driving the energy generation in Type-I X-ray bursts. However, the present reaction-rate uncertainty limits constraints on neutron star properties that can be achieved with burst model-observation comparisons. Here, we present a novel technique for constraining this important reaction by combining the GRETINA array with the neutron detector LENDA coupled to the S800 spectrograph at the National Superconducting Cyclotron Laboratory. The $^{23}$Al($d,n$) reaction was used to populate the astrophysically important states in $^{24}$Si. This enables a measurement in complete kinematics for extracting all relevant inputs necessary to calculate the reaction rate. For the first time, a predicted close-lying doublet of a 2$_2^+$ and (4$_1^+$,0$_2^+$) state in $^{24}$Si was disentangled, finally resolving conflicting results from two previous measurements. Moreover, it was possible to extract spectroscopic factors using GRETINA and LENDA simultaneously. This new technique may be used to constrain other important reaction rates for various astrophysical scenarios

Determination of Alkali and Halide Monovalent Ion Parameters for Use in Explicitly Solvated Biomolecular Simulations

Alkali (Li+, Na+, K+, Rb+, and Cs+) and halide (F−, Cl−, Br−, and I−) ions play an important role in many biological phenomena, roles that range from stabilization of biomolecular structure, to influence on biomolecular dynamics, to key physiological influence on homeostasis and signaling. To properly model ionic interaction and stability in atomistic simulations of biomolecular structure, dynamics, folding, catalysis, and function, an accurate model or representation of the monovalent ions is critically necessary. A good model needs to simultaneously reproduce many properties of ions, including their structure, dynamics, solvation, and moreover both the interactions of these ions with each other in the crystal and in solution and the interactions of ions with other molecules. At present, the best force fields for biomolecules employ a simple additive, nonpolarizable, and pairwise potential for atomic interaction. In this work, we describe our efforts to build better models of the monovalent ions within the pairwise Coulombic and 6-12 Lennard-Jones framework, where the models are tuned to balance crystal and solution properties in Ewald simulations with specific choices of well-known water models. Although it has been clearly demonstrated that truly accurate treatments of ions will require inclusion of nonadditivity and polarizability (particularly with the anions) and ultimately even a quantum mechanical treatment, our goal was to simply push the limits of the additive treatments to see if a balanced model could be created. The applied methodology is general and can be extended to other ions and to polarizable force-field models. Our starting point centered on observations from long simulations of biomolecules in salt solution with the AMBER force fields where salt crystals formed well below their solubility limit. The likely cause of the artifact in the AMBER parameters relates to the naive mixing of the Smith and Dang chloride parameters with AMBER-adapted Åqvist cation parameters. To provide a more appropriate balance, we reoptimized the parameters of the Lennard-Jones potential for the ions and specific choices of water models. To validate and optimize the parameters, we calculated hydration free energies of the solvated ions and also lattice energies (LE) and lattice constants (LC) of alkali halide salt crystals. This is the first effort that systematically scans across the Lennard-Jones space (well depth and radius) while balancing ion properties like LE and LC across all pair combinations of the alkali ions and halide ions. The optimization across the entire monovalent series avoids systematic deviations. The ion parameters developed, optimized, and characterized were targeted for use with some of the most commonly used rigid and nonpolarizable water models, specifically TIP3P, TIP4PEW, and SPC/E. In addition to well reproducing the solution and crystal properties, the new ion parameters well reproduce binding energies of the ions to water and the radii of the first hydration shells

First searches for optical counterparts to gravitational-wave candidate events

During the Laser Interferometer Gravitational-wave Observatory and Virgo joint science runs in 2009-2010, gravitational wave (GW) data from three interferometer detectors were analyzed within minutes to select GW candidate events and infer their apparent sky positions. Target coordinates were transmitted to several telescopes for follow-up observations aimed at the detection of an associated optical transient. Images were obtained for eight such GW candidates. We present the methods used to analyze the image data as well as the transient search results. No optical transient was identified with a convincing association with any of these candidates, and none of the GW triggers showed strong evidence for being astrophysical in nature. We compare the sensitivities of these observations to several model light curves from possible sources of interest, and discuss prospects for future joint GW-optical observations of this type