During planet formation gravitational interaction between a planetary embryo and the protoplanetary gas disc causes orbital migration of the planetary embryo, which plays an important role in shaping the final planetary system. While migration sometimes occurs in the supersonic regime, wherein the relative velocity between the planetary embryo and the gas is higher than the sound speed, migration prescriptions proposed thus far describing the planet-disc interaction force and the timescales of orbital change in the supersonic regime are inconsistent with one another. Here we discuss the details of existing prescriptions in the literature and derive a new simple and intuitive formulation for planet-disc interactions based on dynamical friction that can be applied in both supersonic and subsonic cases. While the existing prescriptions assume particular disc models, ours include the explicit dependence on the disc parameters; hence it can be applied to discs with any radial surface density and temperature dependence (except for the local variations with radial scales less than the disc scale height). Our prescription will reduce the uncertainty originating from different literature formulations of planet migration and will be an important tool to study planet accretion processes, especially when studying the formation of close-in low-mass planets that are commonly found in exoplanetary systems.

The formation of multiple close-in low-mass exoplanets is still a mystery. The challenge is to build a system wherein the outermost planet is beyond 0.2 AU from the star. Here we investigate how the prescription for type I planet migration affects the ability to trap multiple planets in a resonant chain near the inner edge of the protostellar disc. A sharp edge modelled as a hyperbolic tangent function coupled with supersonic corrections to the classical type I migration torques results in the innermost planets being pushed inside the cavity through resonant interaction with farther planets because migration is starward at slightly supersonic eccentricities. Planets below a few Earth masses are generally trapped in a resonant chain with the outermost planet near the disc edge, but long-term stability is not guaranteed. For more massive planets the migration is so fast that the eccentricity of the innermost resonant pair is excited to highly supersonic levels due to decreased damping on the innermost planet as it is pushed inside the cavity; collisions frequently occur and the system consists one or two intermediate-mass planets residing closer to the star than the disc's inner edge. We found a neat pileup of resonant planets outside the disc edge only if the corotation torque does not rapidly diminish at high eccentricity. We call for detailed studies on planet migration near the disc's inner edge, which is still uncertain, and for an improved understanding of eccentricity damping and disc torques in the supersonic regime.

Dynamical models of planet formation coupled with cosmochemical data from martian meteorites show that Mars' isotopic composition is distinct from that of Earth. Reconciliation of formation models with meteorite data require that Mars grew further from the Sun than its present position. Here, we evaluate this compositional difference in more detail by comparing output from two $N$-body planet formation models. The first of these planet formation models simulates what is termed the "Classical" case wherein Jupiter and Saturn are kept in their current orbits. We compare these results with another model based on the "Grand Tack", in which Jupiter and Saturn migrate through the primordial asteroid belt. Our estimate of the average fraction of chondrite assembled into Earth and Mars assumes that the initial solid disk consists of only sources of enstatite chondrite composition in the inner region, and ordinary chondrite in the outer region. Results of these analyses show that both models tend to yield Earth and Mars analogues whose accretion zones overlap. The Classical case fares better in forming Mars with its documented composition (29% to 68% enstatite chondrite plus 32% to 67% ordinary chondrite) though the Mars analogues are generally too massive. However, if we include the restriction of mass on the Mars analogues, the Classical model does not work better. We also further calculate the isotopic composition of $^{17} \rm O$, $^{50} \rm Ti$, $^{54} \rm Cr$, $^{142} \rm Nd$, $^{64} \rm Ni$, and $^{92} \rm Mo$ in the martian mantle from the Grand Tack simulations. We find that it is possible to match the calculated isotopic composition of all the above elements in Mars' mantle with their measured values, but the resulting uncertainties are too large to place good restriction on the early dynamical evolution and birth place of Mars.

Highly volcanic exoplanets, which can be variously characterized as 'lava worlds', 'magma ocean worlds', or 'super-Ios' are high priority targets for investigation. The term 'lava world' may refer to any planet with extensive surface lava lakes, while the term 'magma ocean world' refers to planets with global or hemispherical magma oceans at their surface. 'Highly volcanic planets', including super-Ios, may simply have large, or large numbers of, active explosive or extrusive volcanoes of any form. They are plausibly highly diverse, with magmatic processes across a wide range of compositions, temperatures, activity rates, volcanic eruption styles, and background gravitational force magnitudes. Worlds in all these classes are likely to be the most characterizable rocky exoplanets in the near future due to observational advantages that stem from their preferential occurrence in short orbital periods and their bright day-side flux in the infrared. Transit techniques should enable a level of characterization of these worlds analogous to hot Jupiters. Understanding processes on highly volcanic worlds is critical to interpret imminent observations. The physical states of these worlds are likely to inform not just geodynamic processes, but also planet formation, and phenomena crucial to habitability. Volcanic and magmatic activity uniquely allows chemical investigation of otherwise spectroscopically inaccessible interior compositions. These worlds will be vital to assess the degree to which planetary interior element abundances compare to their stellar hosts, and may also offer pathways to study both the very young Earth, and the very early form of many silicate planets where magma oceans and surface lava lakes are expected to be more prevalent. We suggest that highly volcanic worlds may become second only to habitable worlds in terms of both scientific and public long-term interest.

The geophysics of extrasolar planets is a scientific topic often regarded as standing largely beyond the reach of near-term observations. This reality in no way diminishes the central role of geophysical phenomena in shaping planetary outcomes, from formation, to thermal and chemical evolution, to numerous issues of surface and near-surface habitability. We emphasize that for a balanced understanding of extrasolar planets, it is important to look beyond the natural biases of current observing tools, and actively seek unique pathways to understand exoplanet interiors as best as possible during the long interim prior to a time when internal components are more directly accessible. Such pathways include but are not limited to: (a) enhanced theoretical and numerical modeling, (b) laboratory research on critical material properties, (c) measurement of geophysical properties by indirect inference from imprints left on atmospheric and orbital properties, and (d) the purpose-driven use of Solar System object exploration expressly for its value in comparative planetology toward exoplanet-analogs. Breaking down barriers that envision local Solar System exploration, including the study of Earth's own deep interior, as separate from and in financial competition with extrasolar planet research, may greatly improve the rate of needed scientific progress for exoplanet geophysics. As the number of known rocky and icy exoplanets grows in the years ahead, we expect demand for expertise in 'exogeoscience' will expand at a commensurately intense pace. We highlight key topics, including: how water oceans below ice shells may dominate the total habitability of our galaxy by volume, how free-floating nomad planets may often attain habitable subsurface oceans supported by radionuclide decay, and how deep interiors may critically interact with atmospheric mass loss via dynamo-driven magnetic fields.

Context. Planet formation with pebbles has been proposed to solve a couple of long-standing issues in the classical formation model. Some sophisticated simulations have been done to confirm the efficiency of pebble accretion. However, there has not been any global N-body simulations that compare the outcomes of planet formation via pebble accretion with observed extrasolar planetary systems. Aims. In this paper, we study the effects of a range of initial parameters of planet formation via pebble accretion, and present the first results of our simulations. Methods. We incorporate the pebble accretion model by Ida et al. (2016) in the N-body code SyMBA (Duncan et al. 1998), along with the effects of gas accretion, eccentricity and inclination damping and planet migration in the disc. Results. We confirm that pebble accretion leads to a variety of planetary systems, but have difficulty in reproducing observed properties of exoplanetary systems, such as planetary mass, semimajor axis, and eccentricity distributions. The main reason behind this is a too-efficient type I migration, which sensitively depends on the disc model. However, our simulations also lead to a few interesting predictions. First, we find that formation efficiencies of planets depend on the stellar metallicities, not only for giant planets, but also for Earths (Es) and Super-Earths (SEs). The dependency for Es/SEs is subtle. Although higher metallicity environments lead to faster formation of a larger number of Es/SEs, they also tend to be lost later via dynamical instability. Second, our results indicate that a wide range of bulk densities observed for Es and SEs is a natural consequence of dynamical evolution of planetary systems. Third, the ejection trend of our simulations suggest that one free-floating E/SE may be expected for two smaller-mass planets.

Overcoming type I migration and preventing low-mass planets from spiralling into the central star is a long-studied topic. It is well known that outward migration is possible in viscous-heated discs relatively close to the central star because the entropy gradient can be sufficiently steep that the positive corotation torque overcomes the negative Lindblad torque. Yet efficiently trapping planets in this region remains elusive. Here we study disc conditions that yield outward migration for low-mass planets under specific planet migration prescriptions. Using a steady-state disc model with a constant alpha-viscosity, outward migration is only possible when the negative temperature gradient exceeds ~0.87. We derive an implicit relation for the maximum mass at which outward migration is possible as a function of viscosity and disc scale height. We apply these criteria, using a simple power-law disc model, to planets that have reached their pebble isolation mass after an episode of rapid accretion. It is possible to trap planets with the pebble isolation mass farther than the inner edge of the disc provided that alpha_crit >~ 0.004 for discs older than 1 Myr. In very young discs the high temperature causes the planets to grow to masses exceeding the maximum for outward migration. As the disc evolves, these more massive planets often reach the central star, generally only towards the end of the disc's lifetime. Saving super-Earths is therefore a delicate interplay between disc viscosity, the opacity profile and the temperature gradient in the viscously heated inner disc.

With approximately one ninth of Earth's mass, Mars is widely considered to be a stranded planetary embryo that never became a fully-grown planet. A currently popular planet formation theory predicts that Mars formed near Earth and Venus and was subsequently scattered outwards to its present location. In such a scenario, the compositions of the three planets are expected to be similar to each other. However, bulk elemental and isotopic data for martian meteorites demonstrate that key aspects of Mars' composition are markedly different from that of Earth. This suggests that Mars formed outside of the terrestrial feeding zone during primary accretion. It is therefore probable that Mars always remained significantly farther from the Sun than Earth; its growth was stunted early and its mass remained relatively low. Here we identify a potential dynamical pathway that forms Mars in the asteroid belt and keeps it outside of Earth's accretion zone while at the same time accounting for strict age and compositional constraints, as well as mass differences. Our uncommon pathway (approximately 2% probability) is based on the Grand Tack scenario of terrestrial planet formation, in which the radial migration by Jupiter gravitationally sculpts the planetesimal disc at Mars' current location. We conclude that Mars' formation requires a specific dynamical pathway, while this is less valid for Earth and Venus. We further predict that} Mars' volatile budget is most likely different from Earth's and that Venus formed close enough to our planet that it is expected to have a nearly identical composition from common building blocks.

It is generally accepted that silicate-metal (`rocky') planet formation relies on coagulation from a mixture of sub-Mars sized planetary embryos and (smaller) planetesimals that dynamically emerge from the evolving circum-solar disc in the first few million years of our Solar System. Once the planets have, for the most part, assembled after a giant impact phase, they continue to be bombarded by a multitude of planetesimals left over from accretion. Here we place limits on the mass and evolution of these planetesimals based on constraints from the highly siderophile element (HSE) budget of the Moon. Outcomes from a combination of N-body and Monte Carlo simulations of planet formation lead us to four key conclusions about the nature of this early epoch. First, matching the terrestrial to lunar HSE ratio requires either that the late veneer on Earth consisted of a single lunar-size impactor striking the Earth before 4.45 Ga, or that it originated from the impact that created the Moon. An added complication is that analysis of lunar samples indicates the Moon does not preserve convincing evidence for a late veneer like Earth. Second, the expected chondritic veneer component on Mars is 0.06 weight percent. Third, the flux of terrestrial impactors must have been low ( <=10^(-6) M_earth/Myr) to avoid wholesale melting of Earth's crust after 4.4~Ga, and to simultaneously match the number of observed lunar basins. This conclusion leads to an Hadean eon which is more clement than assumed previously. Last, after the terrestrial planets had fully formed, the mass in remnant planetesimals was ~10^(-3) M_earth, lower by at least an order of magnitude than most previous models suggest. Our dynamically and geochemically self-consistent scenario requires that future N-body simulations of rocky planet formation either directly incorporate collisional grinding or rely on pebble accretion.

The Grand Tack model of terrestrial planet formation has emerged in recent years as the premier scenario used to account for several observed features of the inner solar system. It relies on early migration of the giant planets to gravitationally sculpt and mix the planetesimal disc down to ~1 AU, after which the terrestrial planets accrete from material left in a narrow circum-solar annulus. Here we have investigated how the model fares under a range of initial conditions and migration course-change (`tack') locations. We have run a large number of N-body simulations with a tack location of 1.5 AU and 2 AU and tested initial conditions using equal mass planetary embryos and a semi-analytical approach to oligarchic growth. We make use of a recent model of the protosolar disc that takes account of viscous heating, include the full effect of type 1 migration, and employ a realistic mass-radius relation for the growing terrestrial planets. Results show that the canonical tack location of Jupiter at 1.5 AU is inconsistent with the most massive planet residing at 1 AU at greater than 95% confidence. This favours a tack farther out at 2 AU for the disc model and parameters employed. Of the different initial conditions, we find that the oligarchic case is capable of statistically reproducing the orbital architecture and mass distribution of the terrestrial planets, while the equal mass embryo case is not.

Recent observations started revealing the compositions of protostellar discs and planets beyond the Solar System. In this paper, we explore how the compositions of terrestrial planets are affected by dynamical evolution of giant planets. We estimate the initial compositions of building blocks of these rocky planets by using a simple condensation model, and numerically study the compositions of planets formed in a few different formation models of the Solar System. We find that the abundances of refractory and moderately volatile elements are nearly independent of formation models, and that all the models could reproduce the abundances of these elements of the Earth. The abundances of atmophile elements, on the other hand, depend on the scattering rate of icy planetesimals into the inner disc as well as the mixing rate of the inner planetesimal disc. For the classical formation model, neither of these mechanisms are efficient and the accretion of atmophile elements during the final assembly of terrestrial planets appears to be difficult. For the Grand Tack model, both of these mechanisms are efficient, which leads to a relatively uniform accretion of atmophile elements in the inner disc. It is also possible to have a "hybrid" scenario where the mixing is not very efficient but the scattering is efficient. The abundances of atmophile elements in this case increases with orbital radii. Such a scenario may occur in some of the extrasolar planetary systems which are not accompanied by giant planets or those without strong perturbations from giants. We also confirm that the Grand Tack scenario leads to the distribution of asteroid analogues where rocky planetesimals tend to exist interior to icy ones, and show that their overall compositions are consistent with S-type and C-type chondrites, respectively.

We present new $L^\prime$ (3.8 $\mu m$) and Br-$\alpha$ (4.05 $\mu m$) data and reprocessed archival $L^\prime$ data for the young, planet-hosting star HR 8799 obtained with Keck/NIRC2, VLT/NaCo and Subaru/IRCS. We detect all four HR 8799 planets in each dataset at a moderate to high signal-to-noise (SNR $\gtrsim$ 6-15). We fail to identify a fifth planet, "HR 8799 f", at $r$ $<$ 15 $AU$ at a 5-$\sigma$ confidence level: one suggestive, marginally significant residual at 0.2" is most likely a PSF artifact. Assuming companion ages of 30 $Myr$ and the Baraffe (Spiegel \& Burrows) planet cooling models, we rule out an HR 8799 f with mass of 5 $M_{J}$ (7 $M_{J}$), 7 $M_{J}$ (10 $M_{J}$), and 12 $M_{J}$ (13 $M_{J}$) at $r_{proj}$ $\sim$ 12 $AU$, 9 $AU$, and 5 $AU$, respectively. All four HR 8799 planets have red early T dwarf-like $L^\prime$ - [4.05] colors, suggesting that their SEDs peak in between the $L^\prime$ and $M^\prime$ broadband filters. We find no statistically significant difference in HR 8799 cde's colors. Atmosphere models assuming thick, patchy clouds appear to better match HR 8799 bcde's photometry than models assuming a uniform cloud layer. While non-equilibrium carbon chemistry is required to explain HR 8799 bc's photometry/spectra, evidence for it from HR 8799 de's photometry is weaker. Future, deep IR spectroscopy/spectrophotometry with the Gemini Planet Imager, SCExAO/CHARIS, and other facilities may clarify whether the planets are chemically similar or heterogeneous.

We have performed unbiased spectral line surveys at 3 mm band toward the Galactic circumnuclear disk (CND) and Sgr A* using the Nobeyama Radio Observatory (NRO) 45 m radio telescope. The target positions are two tangential points of the CND and the direction of Sgr A*. We have obtained three wide-band spectra which cover the frequency range from 81.3 GHz to 115.8 GHz, detecting 46 molecular lines from 30 species including 10 rare isotopomers and four hydrogen recombination lines. Each line profile consists of multiple velocity components which arise from the CND, +50 km/s and +20 km/s clouds (GMCs), and the foreground spiral arms. We define the specific velocity ranges which represent the CND and the GMCs toward each direction, and classify the detected lines into three categories: the CND-/GMC-/HBD-types, based on the line intensities integrated over the defined velocity ranges. The CND- and GMC-types are the lines which mainly trace the CND and the GMCs, respectively. The HBD-type possesses the both characteristics of the CND-/GMC-types. We also present the lists of line intensities and other parameters, as well as intensity ratios, which must be useful to investigate the difference between nuclear environments of our Galaxy and of others.

We report the results of millimeter-wave molecular line observations of the Tornado Nebula (G357.7--0.1), which is a bright radio source behind the Galactic Center region. A 15'x15' area was mapped in the J=1--0 lines of CO, 13CO, and HCO+ with the Nobeyama Radio Observatory 45-m telescope. The VLA archival data of OH at 1720 MHz were also reanalyzed. We found two molecular clouds with separate velocities, V_LSR=-14 km/s and +5 km/s. These clouds show rough spatial anti-correlation. Both clouds are associated with OH 1720 MHz emissions in the area overlapping with the Tornado Nebula. The spatial and velocity coincidence indicates violent interaction between the clouds and the Tornado nebula. Modestly excited gas prefers the position of the Tornado "head" in the -14 km/s cloud, also suggesting the interaction. Virial analysis shows that the +5 km/s cloud is more tightly bound by self-gravity than the -14 km/s cloud. We propose a formation scenario for the Tornado Nebula; the +5 km/s cloud collided into the -14 km/s cloud, generating a high-density layer behind the shock front, which activates a putative compact object by Bondi-Hoyle-Lyttleton accretion to eject a pair of bipolar jets.

Out of the handful of asteroids that have been imaged, some have distributions of blocks that are not easily explained. In this paper, we investigate the possibility that seismic shaking leads to the size sorting of particles in asteroids. In particular, we focus on the so-called Brazil Nut Effect (BNE) that separates large particles from small ones under vibrations. We study the BNE over a wide range of parameters by using the N-body code PKDGRAV, and find that the effect is largely insensitive to the coefficients of restitution, but sensitive to friction constants and oscillation speeds. Agreeing with the previous results, we find that convection drives the BNE, where the intruder rises to the top of the particle bed. For the wide-cylinder case, we also observe a "whale" effect, where the intruder follows the convective current and does not stay at the surface. We show that the non-dimensional critical conditions for the BNE agree well with previous studies. We also show that the BNE is scalable for low-gravity environments and that the rise speed of an intruder is proportional to the square root of the gravitational acceleration. Finally, we apply the critical conditions to observed asteroids, and find that the critical oscillation speeds are comparable to the seismic oscillation speeds that are expected from non-destructive impacts.

We report the detection of extremely broad emission toward two molecular clumps in the Galactic central molecular zone. We have mapped the Sagittarius C complex ($-0^\circ.61 < l < -0^\circ.27$, $-0^\circ.29 < b < 0^\circ.04$) in the HCN $J$ = 4--3, $\mathrm{^{13}CO}$ $J$ = 3--2, and $\mathrm{H^{13}CN}$ $J$ = 1--0 lines with the ASTE 10 m and NRO 45 m telescopes, detecting bright emission with $80\mbox{--}120$\ $\mathrm{km\,s^{-1}}$ velocity width (in full-width at zero intensity) toward CO$-0.30$$-0.07$ and CO$-0.40$$-0.22$, which are high velocity compact clouds (HVCCs) identified with our previous CO $J$ = 3--2 survey. Our data reveal an interesting internal structure of CO$-0.30$$-0.07$ comprising a pair of high velocity lobes. The spatial-velocity structure of CO$-0.40$$-0.22$ can be also understood as multiple velocity component, or a velocity gradient across the cloud. They are both located on the rims of two molecular shells of about 10 pc in radius. Kinetic energies of CO$-0.30$$-0.07$ and CO$-0.40$$-0.22$ are $\left(0.8\mbox{--}2\right)\times10^{49}$ erg and $\left(1\mbox{--}4\right)\times10^{49}$ erg, respectively. We propose several interpretations of their broad emission: collision between clouds associated with the shells, bipolar outflow, expansion driven by supernovae (SNe), and rotation around a dark massive object. There scenarios cannot be discriminated because of the insufficient angular resolution of our data, though the absence of visible energy sources associated with the HVCCs seems to favor the cloud--cloud collision scenario. Kinetic energies of the two molecular shells are $1\times10^{51}$ erg and $0.7\times10^{51}$ erg, which can be furnished by multiple SN or hypernova explosions in $2\times10^5$ yr. These shells are candidates of molecular superbubbles created after past active star formation.

Recent observations have shown that at least some close-in exoplanets maintain eccentric orbits despite tidal circularization timescales that are typically shorter than stellar ages. We explore gravitational interactions with a distant planetary companion as a possible cause of these non-zero eccentricities. For simplicity, we focus on the evolution of a planar two-planet system subject to slow eccentricity damping and provide an intuitive interpretation of the resulting long-term orbital evolution. We show that dissipation shifts the two normal eigenmode frequencies and eccentricity ratios of the standard secular theory slightly, and that each mode decays at its own rate. Tidal damping of the eccentricities drives orbits to transition between periods of pericenter circulation and libration, and the planetary system settles into a locked state where the pericenters are nearly aligned or anti-aligned. Once in the locked state, the eccentricities of the two orbits decrease very slowly due to tides rather than at the much more rapid single-planet rate, and thus eccentric orbits, even for close-in planets, can often survive much longer than the age of the system. Assuming that an observed close-in planet on an elliptical orbit is apsidally-locked to a more distant, and perhaps unseen companion, we provide a constraint on the mass, semi-major axis, and eccentricity of the companion. We find the observed two-planet system HAT-P-13 might be in just such an apsidally-locked state, with parameters that obey our constraint well. We also survey close-in single planets, and found that none provide compelling evidence for unseen companions. Instead, we suspect that (1) orbits are circular, (2) tidal damping rates are slower than our assumption, or (3) a recent event has excited these eccentricities. Our method should prove useful for interpreting the results of current and future planet searches.

We analyze new/archival VLT/NaCo and Gemini/NICI high-contrast imaging of the young, self-luminous planet $\beta$ Pictoris b in seven near-to-mid IR photometric filters, using advanced image processing methods to achieve high signal-to-noise, high precision measurements. While $\beta$ Pic b's near-IR colors mimick that of a standard, cloudy early-to-mid L dwarf, it is overluminous in the mid-infrared compared to the field L/T dwarf sequence. Few substellar/planet-mass objects -- i.e. $\kappa$ And b and 1RXJ 1609B -- match $\beta$ Pic b's $JHK_{s}L^\prime$ photometry, and its 3.1 $\mu m$ and 5 $\mu m$ photometry are particularly difficult to reproduce. Atmosphere models adopting cloud prescriptions and large ($\sim$ 60 $\mu m$) dust grains fail to reproduce the $\beta$ Pic b spectrum. However, models incorporating thick clouds similar to those found for HR 8799 bcde but also with small (a few microns) modal particle sizes yield fits consistent with the data within uncertainties. Assuming solar abundance models, thick clouds, and small dust particles ($<a>$ = 4 $\mu m$) we derive atmosphere parameters of log(g) = 3.8 $\pm$ 0.2 and $T_{eff}$ = 1575--1650 $K$, an inferred mass of 7$^{+4}_{-3}$ $M_{J}$, and a luminosity of log(L/L$_{\odot}$) $\sim$ -3.80 $\pm$ 0.02. The best-estimated planet radius, $\approx$ 1.65 $\pm$ 0.06 $R_{J}$, is near the upper end of allowable planet radii for hot-start models given the host star's age and likely reflects challenges with constructing accurate atmospheric models. Alternatively, these radii are comfortably consistent with hot-start model predictions if $\beta$ Pic b is younger than $\approx$ 7 Myr, consistent with a late formation, well after its host star's birth $\sim$ 12$^{+8}_{-4}$ Myr ago.

We mapped molecular gas toward the supernova remnant W44 in the HCO+ J=1-0 line with the Nobeyama Radio Observatory 45 m telescope and in the CO J=3-2 line with the Atacama Submillimeter Telescope Experiment 10 m telescope. High-velocity emission wings were detected in both lines over the area where the radio shell of W44 overlaps the molecular cloud in the plane of the sky. We found that the average velocity distributions of the wing emission can be fitted by a uniform expansion model. The best-fit expansion velocities are 12.2+-0.3 km/s and 13.2+-0.2 km/s in HCO+ and CO, respectively. The non-wing CO J=3-2 component is also fitted by the same model with an expansion velocity of 4.7+-0.1 km/s . This component might be dominated by a post shock higher-density region where the shock velocity had slowed down. The kinetic energy of shocked molecular gas is estimated to be (3.5+-1.3)x10^{49} erg. Adding this and the energy of the previously identified HI shell, we concluded that (1.2+-0.2)x10^{50} erg has been converted into gas kinetic energy from the initial baryonic energy of the W44 supernova. We also found ultra-high-velocity CO J=3-2 wing emission with a velocity width of ~100 km/s at (l, b)=(+34.73d, -0.47d). The origin of this extremely high-velocity wing is a mystery.

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 present Subaru/IRCS J band data for Fomalhaut and a (re)reduction of archival 2004--2006 HST/ACS data first presented by Kalas et al. (2008). We confirm the existence of a candidate exoplanet, Fomalhaut b, in both the 2004 and 2006 F606W data sets at a high signal-to-noise. Additionally, we confirm the detection at F814W and present a new detection in F435W. Fomalhaut b's space motion may be consistent with it being in an apsidally-aligned, non debris ring-crossing orbit, although new astrometry is required for firmer conclusions. We cannot confirm that Fomalhaut b exhibits 0.7-0.8 mag variability cited as evidence for planet accretion or a semi-transient dust cloud. The new, combined optical SED and IR upper limits confirm that emission identifying Fomalhaut b originates from starlight scattered by small dust, but this dust is most likely associated with a massive body. The Subaru and IRAC/4.5 micron upper limits imply M < 2 Mj, still consistent with the range of Fomalhaut b masses needed to sculpt the disk. Fomalhaut b is very plausibly "a planet identified from direct imaging" even if current images of it do not, strictly speaking, show thermal emission from a directly imaged planet.

This paper reports the discovery of a helical molecular cloud in the central molecular zone (CMZ) of our Galaxy. This "pigtail" molecular cloud appears at (l, b, V_LSR) ~ (-0.7deg, +0.0deg, -70 to -30 km/s), with a spatial size of ~ (20 pc)^2 and a mass of (2-6) 10^5 solar masses. This is the third helical gaseous nebula found in the Galactic center region to date. Line intensity ratios indicate that the pigtail molecular cloud has slightly higher temperature and/or density than the other normal clouds in the CMZ. We also found a high-velocity wing emission near the footpoint of this cloud. We propose a formation model of the pigtail molecular cloud. It might be associated with a magnetic tube that is twisted and coiled because of the interaction between clouds in the innermost x_1 orbit and ones in the outermost x_2 orbit.

We present previously unpublished July 2005 $H$-band coronagraphic data of the young, planet-hosting star HR 8799 from the newly-released Keck/NIRC2 archive. Despite poor observing conditions, we detect three of the planetary companions (HR 8799 bcd), two of them (HR 8799 bc) without advanced image processing. Comparing these data with previously published 1998-2011 astrometry and that from re-reduced October 2010 Keck data constrains the orbits of the planets. Analyzing the planets' astrometry separately, HR 8799 d's orbit is likely inclined at least 25$^\circ$ from face-on and the others may be on in inclined orbits. For semimajor axis ratios consistent with a 4:2:1 mean-motion resonance, our analysis yields precise values for HR 8799 bcd's orbital parameters and strictly constrains the planets' eccentricities to be less than 0.18--0.3. However, we find no acceptable orbital solutions with this resonance that place the planets in face-on orbits; HR 8799 d shows the largest deviation from such orbits. Moreover, few orbits make HR 8799 d coplanar with b and c, whereas dynamical stability analyses used to constrain the planets' masses typically assume coplanar and/or face-on orbits. This paper illustrates the significant science gain enabled with the release of the NIRC2 archive.

This letter presents a Nyquist-sampled, high-resolution [CI] 3P1-3P0 map of the -0.2 deg < l < 1.2 deg x -0.1 deg < b < 0 deg region in the Central Molecular Zone (CMZ) taken with the Atacama Submillimeter Telescope Experiment (ASTE) 10 m telescope. We have found that molecular clouds in the CMZ can be classified into two groups according to their [CI]/13CO intensity ratios: a bulk component consisting with clouds with a low, uniform [CI]/13CO ratio (0.45) and another component consisting of clouds with high [CI]/13CO ratios (> 0.8). The [CI]-enhanced regions appear in M-0.02-0.07, the circumnuclear disk, the 180-pc ring and the high velocity compact cloud CO+0.02-0.02. We have carried out a large velocity gradient (LVG) analysis and have derived the C^0/CO column density ratio for M-0.02-0.07 as 0.47, which is approximately twice that of the bulk component of the CMZ (0.26). We propose several hypotheses on the origin of high C^0 abundance in M-0.02-0.07, including cosmic-ray/X-ray dissociation and mechanical dissociation of CO in the pre-existing molecular clouds. We also suggest the possibility that M-0.02-0.07 is a cloud at an early stage of chemical evolution from diffuse gas, which was possibly formed by the bar-induced mass inflow in the Galactic Center region.

We present and analyze a new M' detection of the young exoplanet beta Pictoris b from 2008 VLT/NaCo data at a separation of ~ 4 AU and a high signal-to-noise rereduction of L' data taken in December 2009. Based on our orbital analysis, the planet's orbit is viewed almost perfectly edge-on (i ~ 89 degrees) and has a Saturn-like semimajor axis of 9.50 (+3.93, -1.7) AU. Intriguingly, the planet's orbit is aligned with the major axis of the outer disk (Omega ~ 31 degrees) but probably misaligned with the warp/inclined disk at 80 AU often cited as a signpost for the planet's existence. Our results motivate new studies to clarify how $\beta$ Pic b sculpts debris disk structures and whether a second planet is required to explain the warp/inclined disk.

We present new 1--1.25 micron (z and J band) Subaru/IRCS and 2 micron (K band) VLT/NaCo data for HR 8799 and a rereduction of the 3--5 micron MMT/Clio data first presented by Hinz et al. (2010). Our VLT/NaCo data yields a detection of a fourth planet at a projected separation of ~ 15 AU -- "HR 8799e". We also report new, albeit weak detections of HR 8799b at 1.03 microns and 3.3 microns. Empirical comparisons to field brown dwarfs show that at least HR 8799b and HR8799c, and possibly HR 8799d, have near-to-mid IR colors/magnitudes significantly discrepant from the L/T dwarf sequence. Standard cloud deck atmosphere models appropriate for brown dwarfs provide only (marginally) statistically meaningful fits to HR 8799b and c for unphysically small radii. Models with thicker cloud layers not present in brown dwarfs reproduce the planets' SEDs far more accurately and without the need for rescaling the planets' radii. Our preliminary modeling suggests that HR 8799b has log(g) = 4--4.5, Teff = 900K, while HR 8799c, d, and (by inference) e have log(g) = 4--4.5, Teff = 1000--1200K. Combining results from planet evolution models and new dynamical stability limits implies that the masses of HR 8799b, c, d, and e are 6--7 Mj, 7--10 Mj, 7--10 Mj and 7--10 Mj. 'Patchy" cloud prescriptions may provide even better fits to the data and may lower the estimated surface gravities and masses. Finally, contrary to some recent claims, forming the HR 8799 planets by core accretion is still plausible, although such systems are likely rare.

Recent discoveries of several transiting planets with clearly non-zero eccentricities and some large inclinations started changing the simple picture of close-in planets having circular and well-aligned orbits. Two major scenarios to form such planets are planet migration in a disk, and planet--planet interactions combined with tidal dissipation. The former scenario can naturally produce a circular and low-obliquity orbit, while the latter implicitly assumes an initially highly eccentric and possibly high-obliquity orbit, which are then circularized and aligned via tidal dissipation. We investigate the tidal evolution of transiting planets on eccentric orbits. We show that the current and future orbital evolution of these systems is likely dominated by tidal dissipation, and not by a more distant companion. Although most of these close-in planets experience orbital decay all the way to the Roche limit, there are two characteristic evolution paths for them, depending on the relative efficiency of tidal dissipation inside the star and the planet. We point out that the current observations may be consistent with one of them. Our results suggest that at least some of the close-in planets with non-zero orbital eccentricity may have been formed by tidally circularizing an initially eccentric orbit. We also find that even when the stellar spin-orbit misalignment is observed to be small at present, some systems could have had a highly misaligned orbit in the past. Finally, we also re-examine the recent claim by Levrard et. al., who found that all orbital and spin parameters evolve on a similar timescale to orbital decay.

The discovery of over 400 extrasolar planets allows us to statistically test our understanding of formation and dynamics of planetary systems via numerical simulations. Traditional N-body simulations of multiple-planet systems without gas disks have successfully reproduced the eccentricity (e) distribution of the observed systems, by assuming that the planetary systems are relatively closely packed when the gas disk dissipates, so that they become dynamically unstable within the stellar lifetime. However, such studies cannot explain the small semi-major axes (a) of extrasolar planetary systems, if planets are formed, as the standard planet formation theory suggests, beyond the ice line. In this paper, we numerically study the evolution of three-planet systems in dissipating gas disks, and constrain the initial conditions that reproduce the observed semi-major axis and eccentricity distributions simultaneously. We adopt the initial conditions that are motivated by the standard planet formation theory, and self-consistently simulate the disk evolution, and planet migration by using a hybrid N-body and 1D gas disk code. We also take account of eccentricity damping, and investigate the effect of saturation of corotation resonances on the evolution of planetary systems. We find that the semi-major axis distribution is largely determined in a gas disk, while the eccentricity distribution is determined after the disk dissipation. We also find that there may be an optimum disk mass which leads to the observed a-e distribution. Our simulations generate a larger fraction of planetary systems trapped in mean-motion resonances (MMRs) than the observations, indicating that the disk's perturbation to the planetary orbits may be important to explain the observed rate of MMRs. We also find much lower occurrence of planets on retrograde orbits than the current observations of close-in planets suggest.

Large-scale CO surveys have revealed that the central molecular zone (CMZ) of our Galaxy is characterized by a number of expanding shells/arcs, and by a peculiar population of compact clouds with large velocity widths -- high-velocity compact clouds (HVCCs). We have performed a large-scale CO J=3-2 survey of the CMZ from 2005 to 2008 with the Atacama Submillimeter-wave Telescope Experiment (ASTE). The data cover almost the full extent of the CMZ and the Clump 2 with a 34 arcsec grid spacing. The CO J=3-2 data were compared with the CO J=1--0 data taken with the Nobeyama Radio Observatory 45 m telescope. Molecular gas in the CMZ exhibits higher J=3-2/J=1-0 intensity ratios [R(3-2/1-0)=0.7] than that in the Galactic disk [R(3-2/1-0)=0.4]. We extracted highly excited, optically thin gas by the criterion, R(3-2/1-0)>=1.5. Clumps of high R(3-2/1-0) gas were found in the Sgr A and Sgr C regions, near SNR G 0.9+0.1, and three regions with energetic HVCCs; CO 1.27+0.01, CO -0.41-0.23, and CO -1.21-0.12. We also found a number of small spots of high R(3-2/1-0) gas over the CMZ. Many of these high R(3-2/1-0) clumps and spots have large velocity widths, and some apparently coincides with HVCCs, suggesting that they are spots of shocked molecular gas. Their origin should be local explosive events, possibly supernova explosions. These suggest that the energetic HVCCs are associated with massive compact clusters, which have been formed by microbursts of star formation. Rough estimates of the energy flow from large to small scale suggest that the supernova shocks can make a significant contribution to turbulence activation and gas heating in the CMZ.

The growth of Jovian mass planets during migration in their protoplanetary disks is one of the most important problems that needs to be solved in light of observations of the exosolar planets. Studies of the migration of planets in standard gas disk models routinely show that migration is too fast to form Jovian planets, and that such migrating planetary cores generally plunge into the central stars in less than a Myr. In previous work, we have shown that a poorly ionized, less viscous region in a protoplanetary disk called a dead zone slows down the migration of fixed-mass planets. In this paper, we extend our numerical calculations to include dead zone evolution along with the disk, as well as planet formation via accretion of rocky and gaseous materials. Using our symplectic-integrator-gas dynamics code, we find that dead zones, even in evolving disks wherein migrating planets grow by accretion, still play a fundamental role in saving planetary systems. We demonstrate that Jovian planets form within 2.5 Myr for disks that are ten times more massive than a minimum mass solar nebula (MMSN) with an opacity reduction and without slowing down migration artificially. Our simulations indicate that protoplanetary disks with an initial mass comparable to the MMSN only produce Neptunian mass planets. We also find that planet migration does not help core accretion as much in the oligarchic planetesimal accretion scenario as it was expected in the runaway accretion scenario. Therefore we expect that an opacity reduction (or some other mechanisms) is needed to solve the formation timescale problem even for migrating protoplanets, as long as we consider the oligarchic growth. We also point out a possible role of a dead zone in explaining long-lived, strongly accreting gas disks.

Strong tidal interaction with the central star can circularize the orbits of close-in planets. With the standard tidal quality factor Q of our solar system, estimated circularization times for close-in extrasolar planets are typically shorter than the ages of the host stars. While most extrasolar planets with orbital radii a < 0.1 AU indeed have circular orbits, some close-in planets with substantial orbital eccentricities have recently been discovered. This new class of eccentric close-in planets implies that either their tidal Q factor is considerably higher, or circularization is prevented by an external perturbation. Here we constrain the tidal Q factor for transiting extrasolar planets by comparing their circularization times with accurately determined stellar ages. Using estimated secular perturbation timescales, we also provide constraints on the properties of hypothetical second planets exterior to the known ones.

The ensemble of now more than 250 discovered planetary systems displays a wide range of masses, orbits and, in multiple systems, dynamical interactions. These represent the end point of a complex sequence of events, wherein an entire protostellar disk converts itself into a small number of planetary bodies. Here, we present self-consistent numerical simulations of this process, which produce results in agreement with some of the key trends observed in the properties of the exoplanets. Analogs to our own solar system do not appear to be common, originating from disks near the boundary between barren and (giant) planet-forming.

Observations in the past decade have revealed extrasolar planets with a wide range of orbital semimajor axes and eccentricities. Based on the present understanding of planet formation via core accretion and oligarchic growth, we expect that giant planets often form in closely packed configurations. While the protoplanets are embedded in a protoplanetary gas disk, dissipation can prevent eccentricity growth and suppress instabilities from becoming manifest. However, once the disk dissipates, eccentricities can grow rapidly, leading to close encounters between planets. Strong planet--planet gravitational scattering could produce both high eccentricities and, after tidal circularization, very short-period planets, as observed in the exoplanet population. We present new results for this scenario based on extensive dynamical integrations of systems containing three giant planets, both with and without residual gas disks. We assign the initial planetary masses and orbits in a realistic manner following the core accretion model of planet formation. We show that, with realistic initial conditions, planet--planet scattering can reproduce quite well the observed eccentricity distribution. Our results also make testable predictions for the orbital inclinations of short-period giant planets formed via strong planet scattering followed by tidal circularization.

The tidal interaction between a disk and a planet leads to the planet's migration. A long-standing question regarding this mechanism is how to stop the migration before planets plunge into their central stars. In this paper, we propose a new, simple mechanism to significantly slow down planet migration, and test the possibility by using a hybrid numerical integrator to simulate the disk-planet interaction. The key component of the scenario is the role of low viscosity regions in protostellar disks known as dead zones, which affect planetary migration in two ways. First of all, it allows a smaller-mass planet to open a gap, and hence switch the faster type I migration to the slower type II migration. Secondly, a low viscosity slows down type II migration itself, because type II migration is directly proportional to the viscosity. We present numerical simulations of planetary migration by using a hybrid symplectic integrator-gas dynamics code. Assuming that the disk viscosity parameter inside the dead zone is (alpha=1e-4-1e-5), we find that, when a low-mass planet (e.g. 1-10 Earth masses) migrates from outside the dead zone, its migration is stopped due to the mass accumulation inside the dead zone. When a low-mass planet migrates from inside the dead zone, it opens a gap and slows down its migration. A massive planet like Jupiter, on the other hand, opens a gap and slows down inside the dead zone, independent of its initial orbital radius. The final orbital radius of a Jupiter mass planet depends on the dead zone's viscosity. For the range of alpha's noted above, this can vary anywhere from 7 AU, to an orbital radius of 0.1 AU that is characteristic of the hot Jupiters.

The in-orbit performance and calibration of the Hard X-ray Detector (HXD) on board the X-ray astronomy satellite Suzaku are described. Its basic performances, including a wide energy bandpass of 10-600 keV, energy resolutions of ~4 keV (FWHM) at 40 keV and ~11% at 511 keV, and a high background rejection efficiency, have been confirmed by extensive in-orbit calibrations. The long-term gains of PIN-Si diodes have been stable within 1% for half a year, and those of scintillators have decreased by 5-20%. The residual non-X-ray background of the HXD is the lowest among past non-imaging hard X-ray instruments in energy ranges of 15-70 and 150-500 keV. We provide accurate calibrations of energy responses, angular responses, timing accuracy of the HXD, and relative normalizations to the X-ray CCD cameras using multiple observations of the Crab Nebula.

Most low-mass protostellar disks evolve in clustered environments where they are affected by external radiation fields, while others evolve in more isolated star-forming regions. Assuming that the magneto-rotational instability (MRI) is the main source of viscosity, we calculate the size of a poorly ionized, MRI inactive, and hence low viscosity region - the "dead zone" - in these protostellar disks. We include disk ionization by X-rays, cosmic rays, radioactive elements and thermal collisions, recombination by molecules, metals, and grains, as well as the effect of turbulence stimulation in the dead zone by the active layers lying above it. We also calculate the gap-opening masses of planets, which are determined by a disk's viscosity and a disk aspect ratio, for disks in these environments and compare them with each other. We find that the dead zone is a robust feature of the protostellar disks that is largely independent of their environment, typically stretching out to ~ 15 AU. We analyze the possible effects of dead zones on planet formation, migration, and eccentricity evolution. We show that the gap-opening mass inside the dead zone is expected to be of the order of terrestrial and ice giant mass planets while that outside the dead zone is Jovian or super-Jovian mass planets, largely independent of the star-forming environment. We show that dead zones can significantly slow down both type I and type II planetary migration due to their lower viscosity. We also find that the growth of eccentricity of massive extrasolar planets is particularly favorable through the planet-disk interaction inside the dead zones due to the large gaps expected to be opened by planets.

Protoplanets accrete material from their natal protostellar disks until they are sufficiently massive to open a gap in the face of the disk's viscosity that arises from the magneto-rotational instability (MRI). By computing the ionization structure within observationally well-constrained disk models, we demonstrate that poorly ionized, low viscosity "dead zones" stretch out to 12 AU within typical disks. We find that planets of terrestrial mass robustly form within the dead zones while massive Jovian planets form beyond. Dead zones will also halt the rapid migration of planets into their central stars. Finally, we argue that the gravitational scattering of low mass planets formed in the dead zone, to larger radii by a rapidly accreting Jupiter beyond, can explain the distribution of planetary masses in our solar system.

The final masses of Jovian planets are attained when the tidal torques that they exert on their surrounding protostellar disks are sufficient to open gaps in the face of disk viscosity, thereby shutting off any further accretion. In sufficiently well-ionized disks, the predominant form of disk viscosity originates from the Magneto-Rotational Instability (MRI) that drives hydromagnetic disk turbulence. In the region of sufficiently low ionization rate -- the so-called dead zone -- turbulence is damped and we show that lower mass planets will be formed. We considered three ionization sources (X-rays, cosmic rays, and radioactive elements) and determined the size of a dead zone for the total ionization rate by using a radiative, hydrostatic equilibrium disk model developed by Chiang et al. (2001). We studied a range of surface mass density (Sigma_{0}=10^3 - 10^5 g cm^{-2}) and X-ray energy (kT_{x}=1 - 10 keV). We also compared the ionization rate of such a disk by X-rays with cosmic rays and find that the latter dominate X-rays in ionizing protostellar disks unless the X-ray energy is very high (5 - 10 keV). Among our major conclusions are that for typical conditions, dead zones encompass a region extending out to several AU -- the region in which terrestrial planets are found in our solar system. Our results suggest that the division between low and high mass planets in exosolar planetary systems is a consequence of the presence of a dead zone in their natal protoplanetary disks. We also find that the extent of a dead zone is mainly dependent on the disk's surface mass density. Our results provide further support for the idea that Jovian planets in exosolar systems must have migrated substantially inwards from their points of origin.