We present the first measurement of the planet frequency beyond the "snow line," for the planet-to-star mass-ratio interval −4.5 < log q < −2, corresponding to the range of ice giants to gas giants. We find d 2 N pl d log q d log s = (0.36 ± 0.15) dex −2 at the mean mass ratio q = 5 × 10 −4 with no discernible deviation from a flat (Öpik's law) distribution in logprojected separation s. The determination is based on a sample of six planets detected from intensive follow-up observations of high-magnification (A > 200) microlensing events during 2005-2008. The sampled host stars have a typical mass M host ∼ 0.5 M , and detection is sensitive to planets over a range of planet-star-projected separations (s −1 max R E , s max R E), where R E ∼ 3.5 AU (M host /M) 1/2 is the Einstein radius and s max ∼ (q/10 −4.3) 1/3. This corresponds to deprojected separations roughly three times the "snow line." We show that the observations of these events have the properties of a "controlled experiment," which is what permits measurement of absolute planet frequency. High-magnification events are rare, but the survey-plus-follow-up high-magnification channel is very efficient: half of all high-mag events were successfully monitored and half of these yielded planet detections. The extremely high sensitivity of high-mag events leads to a policy of monitoring them as intensively as possible, independent of whether they show evidence of planets. This is what allows us to construct an unbiased sample. The planet frequency derived from microlensing is a factor 8 larger than the one derived from Doppler studies at factor ∼25 smaller star-planet separations (i.e., periods 2-2000 days). However, this difference is basically consistent with the gradient derived from Doppler studies (when extrapolated well beyond the separations from which it is measured). This suggests a universal separation distribution across 2 dex in planet-star separation, 2 dex in mass ratio, and 0.3 dex in host mass. Finally, if all planetary systems were "analogs" of the solar system, our sample would have yielded 18.2 planets (11.4 "Jupiters," 6.4 "Saturns," 0.3 "Uranuses," 0.2 "Neptunes") including 6.1 systems with two or more planet detections. This compares to six planets including one twoplanet system in the actual sample, implying a first estimate of 1/6 for the frequency of solar-like systems.

Understanding the relation between underlying matter distribution and biased tracers such as galaxies or dark matter halos is essential to extract cosmological information from ongoing or future galaxy redshift surveys. At sufficiently large scales such as the Baryon Acoustic Oscillation (BAO) scale, a standard approach for the bias problem on the basis of the perturbation theory (PT) is to assume the 'local bias' model in which the density field of biased tracers is deterministically expanded in terms of matter density field at the same position. The higher-order bias parameters are then determined by combining the power spectrum with higher-order statistics such as the bispectrum.As is pointed out by recent studies, however, nonlinear gravitational evolution naturally induces nonlocal bias terms even if initially starting only with purely local bias. As a matter of fact, previous works showed that the second-order nonlocal bias term, which corresponds to the gravitational tidal field, is important to explain the characteristic scale-dependence of the bispectrum. In this paper we extend the nonlocal bias term up to third order, and investigate whether the PT-based model including nonlocal bias terms can simultaneously explain the power spectrum and the bispectrum of simulated halos in N -body simulations. The bias renormalization procedure ensures that only one additional term is necessary to be introduced to the power spectrum as a next-to-leading order correction, even if third-order nonlocal bias terms are taken into account. We show that the power spectrum, including density and momentum, and the bispectrum between halo and matter in Nbody simulations can be simultaneously well explained by the model including up to third-order nonlocal bias terms at k < ∼ 0.1h/Mpc. Also, the results are in a good agreement with theoretical predictions of a simple coevolution picture, although the agreement is not perfect. These trend can be found for a wide range of halo mass, 0.7 < ∼ M halo [10 13 M⊙/h] < ∼ 20 at various redshifts, 0 ≤ z ≤ 1. These demonstrations clearly show a failure of the local bias model even at such large scales, and we conclude that nonlocal bias terms should be consistently included in order to accurately model statistics of halos.

Searches for extrasolar planets have uncovered an astonishing diversity of planetary systems, yet the frequency of solar system analogs remains unknown. The gravitational microlensing planet search method is potentially sensitive to multiple-planet systems containing analogs of all the solar system planets except Mercury. We report the detection of a multiple-planet system with microlensing. We identify two planets with masses of approximately 0.71 and approximately 0.27 times the mass of Jupiter and orbital separations of approximately 2.3 and approximately 4.6 astronomical units orbiting a primary star of mass approximately 0.50 solar mass at a distance of approximately 1.5 kiloparsecs. This system resembles a scaled version of our solar system in that the mass ratio, separation ratio, and equilibrium temperatures of the planets are similar to those of Jupiter and Saturn. These planets could not have been detected with other techniques; their discovery from only six confirmed microlensing planet detections suggests that solar system analogs may be common.