Long term variation of the solar diurnal anisotropy of galactic cosmic rays observed with the Nagoya multi-directional muon detector

We analyze the three dimensional anisotropy of the galactic cosmic ray (GCR) intensities observed independently with a muon detector at Nagoya in Japan and neutron monitors over four solar activity cycles. We clearly see the phase of the free-space diurnal anisotropy shifting toward earlier hours around solar activity minima in A>0 epochs, due to the reduced anisotropy component parallel to the mean magnetic field. The average parallel component is consistent with a rigidity independent spectrum, while the perpendicular component increases with GCR rigidity. We suggest that this harder spectrum of the perpendicular component is due to contribution from the drift streaming. We find that the bidirectional latitudinal density gradient is positive in A>0 epoch, while it is negative in A<0 epoch, in accord with the drift model prediction. The radial density gradient, on the other hand, varies with ~11-year cycle with maxima (minima) in solar maximum (minimum) periods, but there is no significant difference seen between average radial gradients in A>0 and A<0 epochs. The average parallel mean free path is larger in A0. We also find, however, that parallel mean free path (radial gradient) appears to persistently increase (decreasing) in the last three cycles of weakening solar activity. We suggest that simple differences between these parameters in A>0 and A<0 epochs are seriously biased by these long-term trends.Comment: accepted for the publication in the Astrophysical Journa

Tidal torques. A critical review of some techniques

We point out that the MacDonald formula for body-tide torques is valid only in the zeroth order of e/Q, while its time-average is valid in the first order. So the formula cannot be used for analysis in higher orders of e/Q. This necessitates corrections in the theory of tidal despinning and libration damping. We prove that when the inclination is low and phase lags are linear in frequency, the Kaula series is equivalent to a corrected version of the MacDonald method. The correction to MacDonald's approach would be to set the phase lag of the integral bulge proportional to the instantaneous frequency. The equivalence of descriptions gets violated by a nonlinear frequency-dependence of the lag. We explain that both the MacDonald- and Darwin-torque-based derivations of the popular formula for the tidal despinning rate are limited to low inclinations and to the phase lags being linear in frequency. The Darwin-torque-based derivation, though, is general enough to accommodate both a finite inclination and the actual rheology. Although rheologies with Q scaling as the frequency to a positive power make the torque diverge at a zero frequency, this reveals not the impossible nature of the rheology, but a flaw in mathematics, i.e., a common misassumption that damping merely provides lags to the terms of the Fourier series for the tidal potential. A hydrodynamical treatment (Darwin 1879) had demonstrated that the magnitudes of the terms, too, get changed. Reinstating of this detail tames the infinities and rehabilitates the "impossible" scaling law (which happens to be the actual law the terrestrial planets obey at low frequencies).Comment: arXiv admin note: sections 4 and 9 of this paper contain substantial text overlap with arXiv:0712.105

An Investigation into the Radial Velocity Variations of CoRoT-7

CoRoT-7b, the first transiting ``superearth'' exoplanet, has a radius of 1.7 R_Earth and a mass of 4.8 M_Earth. Ground-based radial velocity measurements also detected an additional companion with a period of 3.7 days (CoRoT-7c) and a mass of 8.4 M_Earth. The mass of CoRoT-7b is a crucial parameter for planet structure models, but is difficult to determine because CoRoT-7 is a modestly active star and there is at least one additional companion. A Fourier analysis was performed on spectral data for CoRoT-7 taken with the HARPS spectrograph. These data include RV measurements, spectral line bisectors, the full width at half maximum of the cross-correlation function, and Ca II emission. The latter 3 quantities vary due to stellar activity and were used to assess the nature of the observed RV variations. An analysis of a sub-set of the RV measurements where multiple observations were made per night was also used to estimate the RV amplitude from CoRoT-7b that was less sensitive to activity variations. Our analysis indicates that the 0.85-d and 3.7-d RV signals of CoRoT-7b and CoRoT-7c are present in the spectral data with a high degree of statistical significance. We also find evidence for another significant RV signal at 9 days. An analysis of the activity indicator data reveals that this 9-d signal most likely does not arise from activity, but possibly from an additional companion. If due to a planetary companion the mass is m = 19.5 M_Earth, assuming co-planarity with CoRoT-7b. A dynamical study of the three planet system shows that it is stable over several hundred millions of years. Our analysis yields a RV amplitude of 5.04 +/- 1.09 m/s for CoRoT-7b which corresponds to a planet mass of m = 6.9 +/- 1.4 M_Earth. This increased mass would make the planet CoRoT-7b more Earth-like in its internal structure.Comment: 20 pages, 20 figure

Cosmic ray short burst observed with the Global Muon Detector Network (GMDN) on June 22, 2015

We analyze the short cosmic ray intensity increase ("cosmic ray burst": CRB) on June 22, 2015 utilizing a global network of muon detectors and derive the global anisotropy of cosmic ray intensity and the density (i.e. the omnidirectional intensity) with 10-minute time resolution. We find that the CRB was caused by a local density maximum and an enhanced anisotropy of cosmic rays both of which appeared in association with Earth's crossing of the heliospheric current sheet (HCS). This enhanced anisotropy was normal to the HCS and consistent with a diamagnetic drift arising from the spatial gradient of cosmic ray density, which indicates that cosmic rays were drifting along the HCS from the north of Earth. We also find a significant anisotropy along the HCS, lasting a few hours after the HCS crossing, indicating that cosmic rays penetrated into the inner heliosphere along the HCS. Based on the latest geomagnetic field model, we quantitatively evaluate the reduction of the geomagnetic cut-off rigidity and the variation of the asymptotic viewing direction of cosmic rays due to a major geomagnetic storm which occurred during the CRB and conclude that the CRB is not caused by the geomagnetic storm, but by a rapid change in the cosmic ray anisotropy and density outside the magnetosphere.Comment: accepted for the publication in the Astrophysical Journa

Eccentricities of Planets in Binary Systems

The most puzzling property of the extrasolar planets discovered by recent radial velocity surveys is their high orbital eccentricities, which are very difficult to explain within our current theoretical paradigm for planet formation. Current data reveal that at least 25% of these planets, including some with particularly high eccentricities, are orbiting a component of a binary star system. The presence of a distant companion can cause significant secular perturbations in the orbit of a planet. At high relative inclinations, large-amplitude, periodic eccentricity perturbations can occur. These are known as "Kozai cycles" and their amplitude is purely dependent on the relative orbital inclination. Assuming that every planet host star also has a (possibly unseen, e.g., substellar) distant companion, with reasonable distributions of orbital parameters and masses, we determine the resulting eccentricity distribution of planets and compare it to observations? We find that perturbations from a binary companion always appear to produce an excess of planets with both very high (e>0.6) and very low (e<0.1) eccentricities. The paucity of near-circular orbits in the observed sample implies that at least one additional mechanism must be increasing eccentricities. On the other hand, the overproduction of very high eccentricities observed in our models could be combined with plausible circularization mechanisms (e.g., friction from residual gas) to create more planets with intermediate eccentricities (e=0.1-0.6).Comment: 8 pages, to appear in "Close Binaries in the 21st Century: New Opportunities and Challenges", ed. A. Gimenez et al. (Springer

National Geodetic Satellite Program, Part II: Smithsonian Astrophysical Observatory

A sequence of advances in the determination of geodetic parameters presented by the Smithsonian Astrophysical Observatory are described. A Baker-Nunn photographic system was used in addition to a ruby-laser ranging system to obtain data for refinement of geodetic parameters. A summary of the data employed to: (1) derive coordinates for the locations of various tracking stations; and (2) determine the gravitational potential of the earth, is presented

Mapping the ν⊙\nu_\odot Secular Resonance for Retrograde Irregular Satellites

Constructing dynamical maps from the filtered output of numerical integrations, we analyze the structure of the $\nu_\odot$ secular resonance for fictitious irregular satellites in retrograde orbits. This commensurability is associated to the secular angle $\theta = \varpi - \varpi_\odot$, where $\varpi$ is the longitude of pericenter of the satellite and $\varpi_\odot$ corresponds to the (fixed) planetocentric orbit of the Sun. Our study is performed in the restricted three-body problem, where the satellites are considered as massless particles around a massive planet and perturbed by the Sun. Depending on the initial conditions, the resonance presents a diversity of possible resonant modes, including librations of $\theta$ around zero (as found for Sinope and Pasiphae) or 180 degrees, as well as asymmetric librations (e.g. Narvi). Symmetric modes are present in all giant planets, although each regime appears restricted to certain values of the satellite inclination. Asymmetric solutions, on the other hand, seem absent around Neptune due to its almost circular heliocentric orbit. Simulating the effects of a smooth orbital migration on the satellite, we find that the resonance lock is preserved as long as the induced change in semimajor axis is much slower compared to the period of the resonant angle (adiabatic limit). However, the librational mode may vary during the process, switching between symmetric and asymmetric oscillations. Finally, we present a simple scaling transformation that allows to estimate the resonant structure around any giant planet from the results calculated around a single primary mass.Comment: 11 pages, 13 figure