Monte Carlo Procedure for Protein Design

A new method for sequence optimization in protein models is presented. The approach, which has inherited its basic philosophy from recent work by Deutsch and Kurosky [Phys. Rev. Lett. 76, 323 (1996)] by maximizing conditional probabilities rather than minimizing energy functions, is based upon a novel and very efficient multisequence Monte Carlo scheme. By construction, the method ensures that the designed sequences represent good folders thermodynamically. A bootstrap procedure for the sequence space search is devised making very large chains feasible. The algorithm is successfully explored on the two-dimensional HP model with chain lengths N=16, 18 and 32.Comment: 7 pages LaTeX, 4 Postscript figures; minor change

The Lick AGN Monitoring Project 2011: Dynamical Modeling of the Broad Line Region in Mrk 50

We present dynamical modeling of the broad line region (BLR) in the Seyfert 1 galaxy Mrk 50 using reverberation mapping data taken as part of the Lick AGN Monitoring Project (LAMP) 2011. We model the reverberation mapping data directly, constraining the geometry and kinematics of the BLR, as well as deriving a black hole mass estimate that does not depend on a normalizing factor or virial coefficient. We find that the geometry of the BLR in Mrk 50 is a nearly face-on thick disk, with a mean radius of 9.6(+1.2,-0.9) light days, a width of the BLR of 6.9(+1.2,-1.1) light days, and a disk opening angle of 25\pm10 degrees above the plane. We also constrain the inclination angle to be 9(+7,-5) degrees, close to face-on. Finally, the black hole mass of Mrk 50 is inferred to be log10(M(BH)/Msun) = 7.57(+0.44,-0.27). By comparison to the virial black hole mass estimate from traditional reverberation mapping analysis, we find the normalizing constant (virial coefficient) to be log10(f) = 0.78(+0.44,-0.27), consistent with the commonly adopted mean value of 0.74 based on aligning the M(BH)-{\sigma}* relation for AGN and quiescent galaxies. While our dynamical model includes the possibility of a net inflow or outflow in the BLR, we cannot distinguish between these two scenarios.Comment: Accepted for publication in ApJ. 8 pages, 6 figure

The DEHVILS Survey Overview and Initial Data Release: High-Quality Near-Infrared Type Ia Supernova Light Curves at Low Redshift

While the sample of optical Type Ia Supernova (SN Ia) light curves (LCs) usable for cosmological parameter measurements surpasses 2000, the sample of published, cosmologically viable near-infrared (NIR) SN Ia LCs, which have been shown to be good "standard candles," is still $\lesssim$ 200. Here, we present high-quality NIR LCs for 83 SNe Ia ranging from $0.002 < z < 0.09$ as a part of the Dark Energy, H$_0$, and peculiar Velocities using Infrared Light from Supernovae (DEHVILS) survey. Observations are taken using UKIRT's WFCAM, where the median depth of the images is 20.7, 20.1, and 19.3 mag (Vega) for $Y$, $J$, and $H$-bands, respectively. The median number of epochs per SN Ia is 18 for all three bands ($YJH$) combined and 6 for each band individually. We fit 47 SN Ia LCs that pass strict quality cuts using three LC models, SALT3, SNooPy, and BayeSN and find scatter on the Hubble diagram to be comparable to or better than scatter from optical-only fits in the literature. Fitting NIR-only LCs, we obtain standard deviations ranging from 0.128-0.135 mag. Additionally, we present a refined calibration method for transforming 2MASS magnitudes to WFCAM magnitudes using HST CALSPEC stars that results in a 0.03 mag shift in the WFCAM $Y$-band magnitudes.Comment: 24 pages, 9 figures. Accepted by MNRA

A measurement of the millimetre emission and the Sunyaev-Zel'dovich effect associated with low-frequency radio sources

We present a statistical analysis of the millimetre-wavelength properties of 1.4GHz-selected sources and a detection of the Sunyaev–Zel’dovich (SZ) effect associated with the haloes that host them. We stack data at 148, 218 and 277GHz from the Atacama Cosmology Telescope at the positions of a large sample of radio AGN selected at 1.4GHz. The thermal SZ effect associated with the haloes that host the AGN is detected at the 5σ level through its spectral signature, representing a statistical detection of the SZ effect in some of the lowest mass haloes (average M 200 ≈ 10 13 M. h −1 70 ) studied to date. The relation between the SZ effect and mass (based on weak lensing measurements of radio galaxies) is consistent with that measured by Planck for local bright galaxies. In the context of galaxy evolution models, this study confirms that galaxies with radio AGN also typically support hot gaseous haloes. Adding Herschel observations allows us to show that the SZ signal is not significantly contaminated by dust emission. Finally, we analyse the contribution of radio sources to the angular power spectrum of the cosmic microwave background

A 2 R_⊕ Planet Orbiting the Bright Nearby K Dwarf Wolf 503

Since its launch in 2009, the Kepler telescope has found thousands of planets with radii between that of Earth and Neptune. Recent studies of the distribution of these planets have revealed a gap in the population near 1.5–2.0 R⊕, informally dividing these planets into "super-Earths" and "sub-Neptunes." The origin of this division is difficult to investigate directly because the majority of planets found by Kepler orbit distant, dim stars and are not amenable to radial velocity follow-up or transit spectroscopy, making bulk density and atmospheric measurements difficult. Here, we present the discovery and validation of a newly found 2.03^(+0.08)_(-0.07) R⊕ planet in direct proximity to the radius gap, orbiting the bright (J = 8.32 mag), nearby (D = 44.5 pc) high proper motion K3.5V star Wolf 503 (EPIC 212779563). We determine the possibility of a companion star and false positive detection to be extremely low using both archival images and high-contrast adaptive optics images from the Palomar observatory. The brightness of the host star makes Wolf 503b a prime target for prompt radial velocity follow-up, and with the small stellar radius (0.690 ± 0.025R⊙), it is also an excellent target for HST transit spectroscopy and detailed atmospheric characterization with JWST. With its measured radius near the gap in the planet radius and occurrence rate distribution, Wolf 503b offers a key opportunity to better understand the origin of this radius gap as well as the nature of the intriguing populations of "super-Earths" and "sub-Neptunes" as a whole

A 2 Earth Radius Planet Orbiting the Bright Nearby K-Dwarf Wolf 503

Since its launch in 2009, the Kepler telescope has found thousands of planets with radii between that of Earth and Neptune. Recent studies of the distribution of these planets have revealed a rift in the population near 1.5-2.0$R_{\bigoplus}$, informally dividing these planets into "super-Earths" and "sub-Neptunes". The origin of this division is not well understood, largely because the majority of planets found by Kepler orbit distant, dim stars and are not amenable to radial velocity follow-up or transit spectroscopy, making bulk density and atmospheric measurements difficult. Here, we present the discovery and validation of a newly found $2.03^{+0.08}_{-0.07}~R_{\bigoplus}$ planet in direct proximity to the radius gap, orbiting the bright ($J=8.32$~mag), nearby ($D=44.5$~pc) high proper motion star Wolf 503 (EPIC 212779563). We classify Wolf 503 as a K3.5V star and member of the thick disc population. We determine the possibility of a companion star and false positive detection to be extremely low using both archival images and high-contrast adaptive optics images from the Palomar observatory. The brightness of the host star makes Wolf 503b a prime target for prompt radial velocity follow-up, HST transit spectroscopy, as well as detailed atmospheric characterization with JWST. With its measured radius near the gap in the planet radius and occurrence rate distribution, Wolf 503b offers a key opportunity to better understand the origin of this radius gap as well as the nature of the intriguing populations of "super-Earths" and "sub-Neptunes" as a whole

WASP-107b’s Density Is Even Lower: A Case Study for the Physics of Planetary Gas Envelope Accretion and Orbital Migration

With a mass in the Neptune regime and a radius of Jupiter, WASP-107b presents a challenge to planet formation theories. Meanwhile, the planet's low surface gravity and the star's brightness also make it one of the most favorable targets for atmospheric characterization. Here, we present the results of an extensive 4 yr Keck/HIRES radial-velocity (RV) follow-up program of the WASP-107 system and provide a detailed study of the physics governing the accretion of the gas envelope of WASP-107b. We reveal that WASP-107b's mass is only 1.8 Neptune masses (M_b = 30.5 ± 1.7 M_⊕). The resulting extraordinarily low density suggests that WASP-107b has a H/He envelope mass fraction of >85% unless it is substantially inflated. The corresponding core mass of <4.6 M_⊕ at 3σ is significantly lower than what is traditionally assumed to be necessary to trigger massive gas envelope accretion. We demonstrate that this large gas-to-core mass ratio most plausibly results from the onset of accretion at gsim1 au onto a low-opacity, dust-free atmosphere and subsequent migration to the present-day a_b = 0.0566 ± 0.0017 au. Beyond WASP-107b, we also detect a second, more massive planet (M_c sin i = 0.36 ± 0.04MJ ) on a wide eccentric orbit (e _c = 0.28 ± 0.07) that may have influenced the orbital migration and spin–orbit misalignment of WASP-107b. Overall, our new RV observations and envelope accretion modeling provide crucial insights into the intriguing nature of WASP-107b and the system's formation history. Looking ahead, WASP-107b will be a keystone planet to understand the physics of gas envelope accretion

A consensus protocol for the recovery of mercury methylation genes from metagenomes

Mercury (Hg) methylation genes (hgcAB) mediate the formation of the toxic methylmercury and have been identified from diverse environments, including freshwater and marine ecosystems, Arctic permafrost, forest and paddy soils, coal-ash amended sediments, chlor-alkali plants discharges and geothermal springs. Here we present the first attempt at a standardized protocol for the detection, identification and quantification of hgc genes from metagenomes. Our Hg-cycling microorganisms in aquatic and terrestrial ecosystems (Hg-MATE) database, a catalogue of hgc genes, provides the most accurate information to date on the taxonomic identity and functional/metabolic attributes of microorganisms responsible for Hg methylation in the environment. Furthermore, we introduce "marky-coco", a ready-to-use bioinformatic pipeline based on de novo single-metagenome assembly, for easy and accurate characterization of hgc genes from environmental samples. We compared the recovery of hgc genes from environmental metagenomes using the marky-coco pipeline with an approach based on coassembly of multiple metagenomes. Our data show similar efficiency in both approaches for most environments except those with high diversity (i.e., paddy soils) for which a coassembly approach was preferred. Finally, we discuss the definition of true hgc genes and methods to normalize hgc gene counts from metagenomes