Pairs of planets in a system may end up close to their host star on eccentric orbits as a consequence of planet planet scattering, Kozai or secular migration. In this scenario, general relativity and secular perturbations have comparable timescales and may interfere with each other with relevant effects on the eccentricity and pericenter evolution of the two planets. We explore, both analytically and via numerical integration, how the secular evolution is changed by general relativity for a wide range of different initial conditions. We find that when the faster secular frequency approaches the general relativity precession rate, which tipically occurs when the outer planet moves away from the inner one, it relaxes to it and a significant damping of the proper eccentricity of the inner planet occurs. The proper eccentricity of the outer planet is reduced as well due to the changes in the secular interaction of the bodies. The lowering of the peak eccentricities of the two planets during their secular evolution has important implications on their stability. A significant number of two planet systems, otherwise chaotic because of the mutual secular perturbations, are found stable when general relativity is included.

Planet--Planet scattering is an efficient and robust dynamical mechanism for producing eccentric exoplanets. Coupled to tidal interactions with the central star, it can also explain close--in giant planets on circularized and potentially misaligned orbits. We explore scattering events occurring close to the star and test if they can reproduce the main features of the observed orbital distribution of giant exoplanets on tight orbits.In our modeling we exploit a numerical integration code based on the Hermite algorithm and including the effects of general relativity, dynamical tides and two--body collisions.We find that P--P scattering events occurring in systems with three giant planets initially moving on circular orbits close to their star produce a population of planets similar to the presently observed one, including eccentric and misaligned close--in planets. The contribution of tides and general relativity is relevant in determining the final outcome of the chaotic phase. Even if two--body collisions dominate the chaotic evolution of three planets in crossing orbits close to their star, the final distribution shows a significant number of planets on eccentric orbits. The highly misaligned close--in giant planets are instead produced by systems where the initial semi--major axis of the inner planet was around 0.2 au or beyond.

We examined the excitations of planetesimals caused by the resonances of a giant planet in a protoplanetary gas disk. The highly excited planetesimals generate bow shocks, the mechanism of which results in chondrule formation, crystallization of silicate dust, and evaporation of icy planetesimals. The planetesimals beyond 2:1 resonance migrate owing to the gas drag and obtain the maximum eccentricity around 3:1 resonance, which is located at approximately half the planetary distance. The eccentricity depends on the parameters of the planetesimals and the Jovian planet, such as size and location, and gas density of the disk. The maximum relative velocity of a 100-km-sized planetesimal with respect to the gas disk reaches up to ~12 km/s in the case of Jupiter owing to secular resonance, which occurs because of the disk's gravity. We find that if a Jovian mass planet is located within 10 au, the planetesimals larger than 100 km gain sufficient velocity to cause the melting of chondrule precursors and crystallization of the silicate. The maximum velocity is higher for large planetesimals and eccentric planets. Planetesimals are trapped temporarily in the resonances and continue to have high speed over >1 Myr after the formation of a Jovian planet. This duration fits into the timescale of chondrule formation suggested by the isotopic data. The evaporation of icy planetesimals occurs when a Jovian planet is located within 15 au. This mechanism can be a new indicator of planet formation in exosystems if some molecules ejected from icy planetesimals are detected.

We report the detection of planets around two evolved giant stars from radial velocity measurements at Okayama Astrophysical observatory. 24 Boo (G3IV) has a mass of $0.99\,M_{\odot}$, a radius of $10.64\,R_{\odot}$, and a metallicity of ${\rm [Fe/H]}=-0.77$. The star hosts one planet with a minimum mass of $0.91\,M_{\rm Jup}$ and an orbital period of $30.35{\rm d}$. The planet has one of the shortest orbital periods among those ever found around evolved stars by radial-veloocity methods. The stellar radial velocities show additional periodicity with $150{\rm d}$, which are probably attributed to stellar activity. The star is one of the lowest-metallicity stars orbited by planets currently known. $\gamma$ Lib (K0III) is also a metal-poor giant with a mass of $1.47\,M_{\odot}$, a radius of $11.1\,R_{\odot}$, and ${\rm [Fe/H]}=-0.30$. The star hosts two planets with minimum masses of $1.02M_{\rm Jup}$ and $4.58\,M_{\rm Jup}$, and periods of $415{\rm d}$ and $964{\rm d}$, respectively. The star has the second lowest metallicity among the giant stars hosting more than two planets. Dynamical stability analysis for the $\gamma$ Lib system sets a minimum orbital inclination angle to be about $70^{\circ}$ and suggests that the planets are in 7:3 mean-motion resonance, though the current best-fitted orbits to the radial-velocity data are not totally regular.

We report the detection of a double planetary system around the evolved intermediate-mass star HD 47366 from precise radial-velocity measurements at Okayama Astrophysical Observatory, Xinglong Station, and Australian Astronomical Observatory. The star is a K1 giant with a mass of 1.81+-0.13M_sun, a radius of 7.30+-0.33R_sun, and solar metallicity. The planetary system is composed of two giant planets with minimum mass of 1.75^{+0.20}_{-0.17}Mjup and 1.86^{+0.16}_{-0.15}Mjup, orbital period of 363.3^{+2.5}_{-2.4} d and 684.7^{+5.0}_{-4.9} d, and eccentricity of 0.089^{+0.079}_{-0.060} and 0.278^{+0.067}_{-0.094}, respectively, which are derived by a double Keplerian orbital fit to the radial-velocity data. The system adds to the population of multi-giant-planet systems with relatively small orbital separations, which are preferentially found around evolved intermediate-mass stars. Dynamical stability analysis for the system revealed, however, that the best-fit orbits are unstable in the case of a prograde configuration. The system could be stable if the planets were in 2:1 mean-motion resonance, but this is less likely considering the observed period ratio and eccentricity. A present possible scenario for the system is that both of the planets have nearly circular orbits, namely the eccentricity of the outer planet is less than ~0.15, which is just within 1.4sigma of the best-fit value, or the planets are in a mutually retrograde configuration with a mutual orbital inclination larger than 160 degree.

Most of super-Earths detected by the radial velocity (RV) method have significantly smaller eccentricities than the eccentricities corresponding to velocity dispersion equal to their surface escape velocity ("escape eccentricities"). If orbital instability followed by giant impacts among protoplanets that have migrated from outer region is considered, it is usually considered that eccentricities of the merged bodies become comparable to those of orbital crossing bodies, which are excited up to their escape eccentricities by close scattering. However, the eccentricity evolution in the {\it in situ} accretion model has not been studied in detail. Here, we investigate the eccentricity evolution through {\it N}-body simulations. We have found that the merged planets tend to have much smaller eccentricities than the escape eccentricities due to very efficient collision damping. If the protoplanet orbits are initially well separated and their eccentricities are securely increased, an inner protoplanet collides at its apocenter with an outer protoplanet at its pericenter. The eccentricity of the merged body is the smallest for such configuration. Orbital inclinations are also damped by this mechanism and planets tend to share a same orbital plane, which is consistent with {\it Kepler} data. Such efficient collision damping is not found when we start calculations from densely packed orbits of the protoplanets. If the protoplanets are initially in the mean-motion resonances, which corresponds to well separated orbits, the {\it in situ} accretion model well reproduces the features of eccentricities and inclinations of multiple super-Earths/Earth systems discovered by RV and {\it Kepler} surveys.

We discuss the detectability of gravitationally bounded pairs of gas-giant planets (which we call "binary planets") in extrasolar planetary systems that are formed through orbital instability followed by planet-planet dynamical tides during their close encounters, based on the results of N-body simulations by Ochiai, Nagasawa and Ida (Paper I). Paper I showed that the formation probability of a binary is as much as $\sim 10\%$ for three giant planet systems that undergo orbital instability, and after post-capture long-term tidal evolution, the typical binary separation is 3--5 times the sum of physical radii of the planets. The binary planets are stable during main sequence lifetime of solar-type stars, if the stellarcentric semimajor axis of the binary is larger than 0.3 AU. We show that detecting modulations of transit light curves is the most promising observational method to detect binary planets. Since the likely binary separations are comparable to the stellar diameter, the shape of the transit light curve is different from transit to transit, depending on the phase of the binary's orbit. The transit durations and depth for binary planet transits are generally longer and deeper than those for the single planet case. We point out that binary planets could exist among the known inflated gas giant planets or objects classified as false positive detections at orbital radii > 0.3 AU, propose a binary planet explanation for the CoRoT candidate SRc01 E2 1066, and show that binary planets are likely to be present in, and could be detected using Kepler-quality data.

It is proposed that planetesimals perturbed by Jovian mean-motion resonances are the source of shock waves that form chondrules. It is considered that this shock-induced chondrule formation requires the velocity of the planetesimal relative to the gas disk to be on the order of > 7 km/s at 1 AU. In previous studies on planetesimal excitation, the effects of Jovian mean-motion resonance together with the gas drag were investigated, but the velocities obtained were at most 8 km/s in the asteroid belt, which is insufficient to account for the ubiquitous existence of chondrules. In this paper, we reexamine the effect of Jovian resonances and take into account the secular resonance in the asteroid belt caused by the gravity of the gas disk. We find that the velocities relative to the gas disk of planetesimals a few hundred kilometers in size exceed 12 km/s, and that this is achieved around the 3:1 mean-motion resonance. The heating region is restricted to a relatively narrow band between 1.5 AU and 3.5 AU. Our results suggest that chondrules were produced effectively in the asteroid region after Jovian formation. We also find that many planetesimals are scattered far beyond Neptune. Our findings can explain the presence of crystalline silicate in comets if the scattered planetesimals include silicate dust processed by shock heating.

We have investigated i) the formation of gravitationally bounded pairs of gas-giant planets (which we call "binary planets") from capturing each other through planet-planet dynamical tide during their close encounters and ii) the following long-term orbital evolution due to planet-planet and planet-star {\it quasi-static} tides. For the initial evolution in phase i), we carried out N-body simulations of the systems consisting of three jupiter-mass planets taking into account the dynamical tide. The formation rate of the binary planets is as much as 10% of the systems that undergo orbital crossing and this fraction is almost independent of the initial stellarcentric semi-major axes of the planets, while ejection and merging rates sensitively depend on the semi-major axes. As a result of circularization by the planet-planet dynamical tide, typical binary separations are a few times the sum of the physical radii of the planets. After the orbital circularization, the evolution of the binary system is governed by long-term quasi-static tide. We analytically calculated the quasi-static tidal evolution in later phase ii). The binary planets first enter the spin-orbit synchronous state by the planet-planet tide. The planet-star tide removes angular momentum of the binary motion, eventually resulting in a collision between the planets. However, we found that the binary planets survive the tidal decay for main-sequence life time of solar-type stars (~10Gyrs), if the binary planets are beyond ~0.3 AU from the central stars. These results suggest that the binary planets can be detected by transit observations at >0.3AU.

The KOI-94 system is a closely-packed, multi-transiting planetary system discovered by the Kepler space telescope. It is known as the first system that exhibited a rare event called a "planet-planet eclipse (PPE)," in which two planets partially overlap with each other in their double-transit phase. In this paper, we constrain the parameters of the KOI-94 system with an analysis of the transit timing variations (TTVs). Such constraints are independent of the radial velocity (RV) analysis recently performed by Weiss and coworkers, and valuable in examining the reliability of the parameter estimate using TTVs. We numerically fit the observed TTVs of KOI-94c, KOI-94d, and KOI-94e for their masses, eccentricities, and longitudes of periastrons, and obtain the best-fit parameters including $m_{\rm c} = 9.4_{-2.1}^{+2.4} M_{\oplus}$, $m_{\rm d} = 52.1_{-7.1}^{+6.9} M_{\oplus}$, $m_{\rm e} = 13.0_{-2.1}^{+2.5} M_{\oplus}$, and $e \lesssim 0.1$ for all the three planets. While these values are mostly in agreement with the RV result, the mass of KOI-94d estimated from the TTV is significantly smaller than the RV value $m_{\rm d} = 106 \pm 11 M_{\oplus}$. In addition, we find that the TTV of the outermost planet KOI-94e is not well reproduced in the current modeling. We also present analytic modeling of the PPE and derive a simple formula to reconstruct the mutual inclination of the two planets from the observed height, central time, and duration of the brightening caused by the PPE. Based on this model, the implication of the results of TTV analysis for the time of the next PPE is discussed.

The ubiquity of planets and diversity of planetary systems reveal planet formation encompass many complex and competing processes. In this series of papers, we develop and upgrade a population synthesis model as a tool to identify the dominant physical effects and to calibrate the range of physical conditions. Recent planet searches leads to the discovery of many multiple-planet systems. Any theoretical models of their origins must take into account dynamical interaction between emerging protoplanets. Here, we introduce a prescription to approximate the close encounters between multiple planets. We apply this method to simulate the growth, migration, and dynamical interaction of planetary systems. Our models show that in relatively massive disks, several gas giants and rocky/icy planets emerge, migrate, and undergo dynamical instability. Secular perturbation between planets leads to orbital crossings, eccentricity excitation, and planetary ejection. In disks with modest masses, two or less gas giants form with multiple super-Earths. Orbital stability in these systems is generally maintained and they retain the kinematic structure after gas in their natal disks is depleted. These results reproduce the observed planetary mass-eccentricity and semimajor axis-eccentricity correlations. They also suggest that emerging gas giants can scatter residual cores to the outer disk regions. Subsequent in situ gas accretion onto these cores can lead to the formation of distant (> 30AU) gas giants with nearly circular orbits.

The orbital distributions of currently observed extrasolar giant planets allow marginally stable orbits for hypothetical, terrestrial planets. In this paper, we propose that many of these systems may not have additional planets on these "stable" orbits, since past dynamical instability among giant planets could have removed them. We numerically investigate the effects of early evolution of multiple giant planets on the orbital stability of the inner, sub-Neptune-like planets which are modeled as test particles, and determine their dynamically unstable region. Previous studies have shown that the majority of such test particles are ejected out of the system as a result of close encounters with giant planets. Here, we show that secular perturbations from giant planets can remove test particles at least down to 10 times smaller than their minimum pericenter distance. Our results indicate that, unless the dynamical instability among giant planets is either absent or quiet like planet-planet collisions, most test particles down to ~0.1 AU within the orbits of giant planets at a few AU may be gone. In fact, out of ~30 % of survived test particles, about three quarters belong to the planet-planet collision cases. We find a good agreement between our numerical results and the secular theory, and present a semi-analytical formula which estimates the dynamically unstable region of the test particles just from the evolution of giant planets. Finally, our numerical results agree well with the observations, and also predict the existence of hot rocky planets in eccentric giant planet systems.

We report the detection of a double planetary system orbiting around the evolved intermediate-mass star HD 4732 from precise Doppler measurements at Okayama Astrophysical Observatory (OAO) and Anglo-Australian Observatory (AAO). The star is a K0 subgiant with a mass of 1.7 M_sun and solar metallicity. The planetary system is composed of two giant planets with minimum mass of msini=2.4 M_J, orbital period of 360.2 d and 2732 d, and eccentricity of 0.13 and 0.23, respectively. Based on dynamical stability analysis for the system, we set the upper limit on the mass of the planets to be about 28 M_J (i>5 deg) in the case of coplanar prograde configuration.

We report that HAT-P-7 has a common proper motion stellar companion. The companion is located at $\sim3.9$ arcsec to the east and estimated as an M5.5V dwarf based on its colors. We also confirm the presence of the third companion, which was first reported by Winn et al. (2009), based on long-term radial velocity measurements. We revisit the migration mechanism of HAT-P-7b given the presence of those companions, and propose sequential Kozai migration as a likely scenario in this system. This scenario may explain the reason for an outlier in the discussion of the spin-orbit alignment timescale for HAT-P-7b by Albrecht et al. (2012).

Many extrasolar planetary systems containing multiple super-Earths have been discovered. N-body simulations taking into account standard type-I planetary migration suggest that protoplanets are captured into mean-motion resonant orbits near the inner disk edge at which the migration is halted. Previous N-body simulations suggested that orbital stability of the resonant systems depends on number of the captured planets. In the unstable case, through close scattering and merging between planets, non-resonant multiple systems are finally formed. In this paper, we investigate the critical number of the resonantly trapped planets beyond which orbital instability occurs after disk gas depletion. We find that when the total number of planets ($N$) is larger than the critical number ($N_{\rm crit}$), crossing time that is a timescale of initiation of the orbital instability is similar to non-resonant cases, while the orbital instability never occurs within the orbital calculation time ($10^8$ Kepler time) for $N\leq N_{\rm crit}$. Thus, the transition of crossing time across the critical number is drastic. When all the planets are trapped in 7:6 resonance of adjacent pairs, $N_{\rm crit} = 4$. We examine the dependence of the critical number of 4:3, 6:5 and 8:7 resonance by changing the orbital separation in mutual Hill radii and planetary mass. The critical number increases with increasing the orbital separation in mutual Hill radii with fixed planetary mass and increases with increasing planetary mass with fixed the orbital separation in mutual Hill radii. We also calculate the case of a system which is not composed of the same resonance. The sharp transition of the stability can be responsible for the diversity of multiple super-Earths (non-resonant or resonant), that is being revealed by $Kepler$ mission.

We present results of direct imaging observations for HAT-P-7 taken with the Subaru HiCIAO and the Calar Alto AstraLux. Since the close-in transiting planet HAT-P-7b was reported to have a highly tilted orbit, massive bodies such as giant planets, brown dwarfs, or a binary star are expected to exist in the outer region of this system. We show that there are indeed two candidates for distant faint stellar companions around HAT-P-7. We discuss possible roles played by such companions on the orbital evolution of HAT-P-7b. We conclude that as there is a third body in the system as reported by Winn et al. (2009, ApJL, 763, L99), the Kozai migration is less likely while planet-planet scattering is possible.

The extrasolar planets (EPs) so far detected are very different to the planets in our own Solar System. Many of them have Jupiter-like masses and close-in orbits (the so-called hot planets, HPs), with orbital periods of only a few days. In this paper, we present a new statistical analysis of the observed EPs, focusing on the origin of the HPs. Among the several HP formation mechanisms proposed so far, the two main formation mechanisms are type II migration and scattering. In both cases, planets form beyond the so-called snow-line of the protoplanetary disk and then migrate inward due to angular momentum and energy exchange with either the protoplanetary disk or with companion planets. Although theoretical studies produce a range of observed features, no firm correspondence between the observed EPs and models has yet been established. In our analysis, by means of principal component analysis and hierarchical cluster analysis, we find convincing indications for the existence of two types of HPs, whose parameters reflect physical mechanisms of type II migration and scattering.

We revisit the "dynamical shakeup" model of Solar System terrestrial planet formation, wherein the whole process is driven by the sweeping of Jupiter's secular resonance as the gas disk is removed. Using a large number of 0.5 Gyr-long N-body simulations, we investigate the different outcomes produced by such a scenario. We confirm that in contrast to existing models, secular resonance sweeping combined with tidal damping by the disk gas can reproduce the low eccentricities and inclinations, and high radial mass concentration, of the Solar System terrestrial planets. At the same time, this also drives the final assemblage of the planets on a timescale of several tens of millions of years, an order of magnitude faster than inferred from previous numerical simulations which neglected these effects, but possibly in better agreement with timescales inferred from cosmochemical data. In addition, we find that significant delivery of water-rich material from the outer Asteroid Belt is a natural byproduct.

We have investigated the formation of close-in extrasolar giant planets through a coupling effect of mutual scattering, Kozai mechanism, and tidal circularization, by orbital integrations. We have carried out orbital integrations of three planets with Jupiter-mass, directly including the effect of tidal circularization. We have found that in about 30% runs close-in planets are formed, which is much higher than suggested by previous studies. We have found that Kozai mechanism by outer planets is responsible for the formation of close-in planets. During the three-planet orbital crossing, the Kozai excitation is repeated and the eccentricity is often increased secularly to values close enough to unity for tidal circularization to transform the inner planet to a close-in planet. Since a moderate eccentricity can remain for the close-in planet, this mechanism may account for the observed close-in planets with moderate eccentricities and without nearby secondary planets. Since these planets also remain a broad range of orbital inclinations (even retrograde ones), the contribution of this process would be clarified by more observations of Rossiter-McLaughlin effects for transiting planets.

The formation of gas giant planets is assumed to be preceded by the emergence of solid cores in the conventional sequential-accretion paradigm. This hypothesis implies that the presence of earth-like planets can be inferred from the detection of gas giants. A similar prediction cannot be made with the gravitational instability (hereafter GI) model which assumes that gas giants (hereafter giants) formed from the collapse of gas fragments analogous to their host stars. We propose an observational test for the determination of the dominant planet-formation channel. Based on the sequential-accretion (hereafter SA) model, we identify several potential avenues which may lead to the prolific formation of a population of close-in earth-mass ($M_\oplus$) planets (hereafter close-in earths) around stars with 1) short-period or 2) solitary eccentric giants and 3) systems which contain intermediate-period resonant giants. In contrast, these close-in earths are not expected to form in systems where giants originated rapidly through GI. As a specific example, we suggest that the SA processes led to the formation of the 7.5 $M_\oplus$ planet around GJ 876 and predict that it may have an atmosphere and envelope rich in O$_2$ and liquid water. Assessments of the ubiquity of these planets will lead to 1) the detection of the first habitable terrestrial planets, 2) the verification of the dominant mode of planet formation, 3) an estimate of the fraction of earth-harboring stars, and 4) modification of bio-marker signatures.

We study the possible signals of a pion string associated with the QCD chiral phase transition in LHC Pb - Pb collision at energy $ \sqrt{s}=5.5$ TeV. In terms of the Kibble-Zurek mechanism we discuss the production and evolution of the pion string. The pion string is not topologically stable, it decays into neutral pions and sigma mesons which in turn decay into pions. Our results show that all the neutral pions from the pion string are distributed at the low momentum and the ratio of neutral to charged pions from the pion string violates the isospin symmetry. For the momentum spectra of the total pions, the signal from the sigma particle decay which is from the pion string will be affected by the large decay width of the sigma significantly.

Most extrasolar planets are observed to have eccentricities much larger than those in the solar system. Some of these planets have sibling planets, with comparable masses, orbiting around the same host stars. In these multiple planetary systems, eccentricity is modulated by the planets' mutual secular interaction as a consequence of angular momentum exchange between them. For mature planets, the eigenfrequencies of this modulation are determined by their mass and semi-major axis ratios. But, prior to the disk depletion, self gravity of the planets' nascent disks dominates the precession eigenfrequencies. We examine here the initial evolution of young planets' eccentricity due to the apsidal libration or circulation induced by both the secular interaction between them and the self gravity of their nascent disks. We show that as the latter effect declines adiabatically with disk depletion, the modulation amplitude of the planets' relative phase of periapse is approximately invariant despite the time-asymmetrical exchange of angular momentum between planets. However, as the young planets' orbits pass through a state of secular resonance, their mean eccentricities undergo systematic quantitative changes. For applications, we analyze the eccentricity evolution of planets around Upsilon Andromedae and HD168443 during the epoch of protostellar disk depletion. We find that the disk depletion can change the planets' eccentricity ratio. However, the relatively large amplitude of the planets' eccentricity cannot be excited if all the planets had small initial eccentricities.

The formation probability of multiple winding topological defects is calculated by phase and flux distribution analysis based on the Kibble mechanism. The core size of defects is taken into account so that when it is much larger than the correlation length of the Higgs field, high winding configuration can be realized. For example, if the self coupling constant of the Higgs field is smaller than $4\times10^{-4}$, the topological inflation may occur at the grand unification energy scale by the multiple winding string produced during the GUT phase transition.

In models in which some of the scalar field are charged and some uncharged, interactions with a finite temperature plasma will lead to corrections to the effective potential of the charged fields which may stabilize embedded defects made up of uncharged fields. In these models embedded defects, solutions of the field equations which are unstable at zero temperature, may thus be stable in the early Universe, and may then play an important role in cosmology. Two prototypical examples are the pion string in the theory of strong interactions, and the electroweak Z-string in the standard electroweak theory.

The Jeans scale is estimated during the coexistence epoch of quark-gluon and hadron phases in the first-order QCD phase transition. It is shown that, contrary to the previous claims, reduction of the sound speed is so little that the phase transition does not affect evolution of cosmological density fluctuations appreciably.

We examine the collapse of an axion domain wall bounded by an axionic string. It is found that the collapse proceeds quickly and axion domain walls disappear. However axions are emitted in the collapse and its energy density increases during radiation dominated era and contributes significantly to the present mass density of the universe. In particular the axion emitted from the wall can account for the dark matter in the universe for $F_a\gsim 10^{10}\GeV$.

Time evolution of a circular cosmic string loop is investigated by numerically solving the field equations for the scalar and the gauge fields consisting of the vortex. It is shown that the result agrees with an analytic estimate based on the Nambu-Goto action, which supports its validity in analyzing nonstraight and rapidly moving strings.