Accurate measurements of stellar metallicity gradients in the radial and vertical directions of the disk and their temporal variations provide important constraints on the formation and evolution of the Milky Way disk...Accurate measurements of stellar metallicity gradients in the radial and vertical directions of the disk and their temporal variations provide important constraints on the formation and evolution of the Milky Way disk. We use 297 042 main sequence turn-off stars selected from the LAMOST Spectroscopic Survey of the Galactic Anticenter（LSS-GAC） to determine the radial and vertical gradients of stellar metallicity,△[Fe/H]/△R and △[Fe/H]/△｜Z ｜ of the Milky Way disk in the direction of the anticenter. We determine ages of those turn-off stars by isochrone fitting and measure the temporal variations of metallicity gradients. We have carried out a detailed analysis of the selection effects resulting from the selection, observation and data reduction of LSS-GAC targets and the potential biases of a magnitude limited sample on the determinations of metallicity gradients. Our results show that the gradients, both in the radial and vertical directions, exhibit significant spatial and temporal variations. The radial gradients yielded by stars with the oldest ages（ 11 Gyr） are essentially zero at all heights from the disk midplane, while those given by younger stars are always negative. The vertical gradients deduced from stars with the oldest ages（ 11 Gyr）are negative and only show very weak variations with Galactocentric distance in the disk plane, R, while those yielded by younger stars show strong variations with R.After being essentially flat at the earliest epochs of disk formation, the radial gradients steepen as age decreases, reaching a maximum（steepest） at age 7–8 Gyr, and then they flatten again. Similar temporal trends are also found for the vertical gradients. We infer that the assembly of the Milky Way disk may have experienced at least two distinct phases. The earlier phase is probably related to a slow, pressure-supported collapse of gas, when the gas settles down to the disk mainly in the vertical direction. In the later phase, there are significant radial flows of gas in the disk, and the rate of gas inflow near the solar neighborhood reaches a maximum around a lookback time of 7–8 Gyr.The transition between the two phases occurs around a lookback time between 8 and11 Gyr. The two phases may be responsible for the formation of the Milky Way＇s thick and thin disks, respectively. Also, as a consequence, we recommend that stellar age is a natural, physical criterion to distinguish stars from the thin and thick disks. From an epoch earlier than 11 Gyr to one between 8 and 11 Gyr, there is an abrupt, significant change in magnitude of both the radial and vertical metallicity gradients, suggesting that stellar radial migration is unlikely to play an important role in the formation of the thick disk.展开更多
In the tandem planet formation regime,planets form at two distinct sites where solid particles are densely accumulated due to the on/off state of the magnetorotational instability（MRI）.We found that tandem planet fo...In the tandem planet formation regime,planets form at two distinct sites where solid particles are densely accumulated due to the on/off state of the magnetorotational instability（MRI）.We found that tandem planet formation can reproduce the solid component distribution of the Solar System and tends to produce a smaller number of large planets through continuous pebble flow into the planet formation sites.In the present paper,we investigate the dependence of tandem planet formation on the vertical magnetic field of the protoplanetary disk.We calculated two cases of BZ 3.4 × 10^-3 G and BZ = 3.4 × 10^-5 G at 100 AU as well as the canonical case of BZ = 3.4 × 10^-4 G.We found that tandem planet formation holds up well in the case of the strong magnetic field（BZ 3.4 × 10^-3 G）.On the other hand,in the case of a weak magnetic field（BZ= 3.4 × 10^-5 G） at 100 AU,a new regime of planetary growth is realized：the planets grow independently at different places in the dispersed area of the MRl-suppressed region of r-8-30 AU at a lower accretion rate of M 〈 10^-7.4M⊙yr^-1.We call this the ＂dispersed planet formation＂ regime.This may lead to a system with a larger number of smaller planets that gain high eccentricity through mutual collisions.展开更多
We consider the existence of a neutron star magnetic field by the detected cyclotron lines. We collected data on nine sources of high-mass X-ray binaries with supergiant companions as a test case for our model, to dem...We consider the existence of a neutron star magnetic field by the detected cyclotron lines. We collected data on nine sources of high-mass X-ray binaries with supergiant companions as a test case for our model, to demonstrate their distribution and evolution. The wind velocity, spin period and magnetic field strength are studied under different mass loss rates. In our model, correlations between mass-loss rate and wind velocity are found and can be tested in further observations. We examine the parameter space where wind accretion is allowed, avoiding the barrier of rotating magnetic fields, with robust data on the magnetic field of neutron stars. Our model shows that most sources(six of nine systems) can be fed by the wind with relatively slow velocity, and this result is consistent with previous predictions. In a few sources,our model cannot fit the standard wind accretion scenario. In these peculiar cases, other scenarios(disk formation, partial Roche lobe overflow) should be considered. This would provide information about the evolutionary tracks of various types of binaries, and thus exhibit a clear dichotomy behavior in wind-fed X-ray binary systems.展开更多
We present the results of our recent study on the interactions between a giant planet and a self-gravitating gas disk. We investigate how the disk＇s self-gravity affects the gap formation process and the migration of...We present the results of our recent study on the interactions between a giant planet and a self-gravitating gas disk. We investigate how the disk＇s self-gravity affects the gap formation process and the migration of the giant planet. Two series of 1-D and 2-D hydrodynamic simulations are performed. We select several surface densities and focus on the gravitationally stable region. To obtain more reliable gravity torques exerted on the planet, a refined treatment of the disk＇s gravity is adopted in the vicinity of the planet. Our results indicate that the net effect of the disk＇s self- gravity on the gap formation process depends on the surface density of the disk. We notice that there are two critical values, ∑I and ∑n. When the surface density of the disk is lower than the first one,∑0 〈 ∑I, the effect of self-gravity suppresses the formation of a gap. When ∑0 〉 ∑I, the self-gravity of the gas tends to benefit the gap formation process and enlarges the width/depth of the gap. According to our 1-D and 2-D simulations, we estimate the first critical surface density to be ∑I ≈ 0.8 MMSN. This effect increases until the surface density reaches the second critical value ∑n- When ∑0 〉 ∑n, the gravitational turbulence in the disk becomes dominant and the gap formation process is suppressed again. Our 2-D simulations show that this critical surface density is around 3.5 MMSN. We also study the associated orbital evolution of a giant planet. Under the effect of the disk＇s self-gravity, the migration rate of the giant planet increases when the disk is dominated by gravitational turbulence. We show that the migration timescale correlates with the effective viscosity and can be up to 104 yr.展开更多
We investigate the behavior of the snowline in a protoplanetary disk and the relationship between the radius of the snowline and properties of molecular cloud cores.In our disk model,we consider mass influx from the g...We investigate the behavior of the snowline in a protoplanetary disk and the relationship between the radius of the snowline and properties of molecular cloud cores.In our disk model,we consider mass influx from the gravitational collapse of a molecular cloud core,irradiation from the central star,and thermal radiation from the ambient molecular cloud gas.As the protoplanetary disk evolves,the radius of the snowline increases first to a maximum value Rmax,and then decreases in the late stage of evolution of the protoplanetary disk.The value of Rmaxis dependent on the properties of molecular cloud cores(mass M;,angular velocity ω and temperature T;).Many previous works found that solid material tends to accumulate at the location of the snowline,which suggests that the snowline is the preferred location for giant planet formation.With these conclusions,we compare the values of R;with semimajor axes of giant planets in extrasolar systems,and find that Rmaxmay provide an upper limit for the locations of the formation of giant planets which are formed by the core accretion model.展开更多
We present a new united theory of planet formation,which includes magneto-rotational instability（MRl） and porous aggregation of solid particles in a consistent way.We show that the ＂tandem planet formation＂ regime...We present a new united theory of planet formation,which includes magneto-rotational instability（MRl） and porous aggregation of solid particles in a consistent way.We show that the ＂tandem planet formation＂ regime is likely to result in solar system-like planetary systems.In the tandem planet formation regime,planetesimals form at two distinct sites：the outer and inner edges of the MRl suppressed region.The former is likely to be the source of the outer gas giants,and the latter is the source for the inner volatile-free rocky planets.Our study spans disks with a various range of accretion rates,and we find that tandem planet formation can occur for M = 10^7.3- 10^-6.9Myr^-1.The rocky planets form between 0.4-2 AU,while the icy planets form between 6-30 All;no planets form in 2-6 AU region for any accretion rate.This is consistent with the gap in the solid component distribution in the solar system,which has only a relatively small Mars and a very small amount of material in the main asteroid belt from 2-6 AU.The tandem regime is consistent with the idea that the Earth was initially formed as a completely volatile-free planet.Water and other volatile elements came later through the accretion of icy material by occasional inward scattering from the outer regions.Reactions between reductive minerals,such as schreibersite（Fe-jP）,and water are essential to supply energy and nutrients for primitive life on Earth.展开更多
We studied the particle growth in a protoplanetary disk in a high-ionization environment and found that icy planet formation is inactive for a disk with an ionization rate 100 times higher than that of the present Sol...We studied the particle growth in a protoplanetary disk in a high-ionization environment and found that icy planet formation is inactive for a disk with an ionization rate 100 times higher than that of the present Solar System. In particular, in the case of M 〈 10^（-7.4）M_☉yr^（-1）, only rocky planet formation occurs. In such a case, all the solid materials in the disk drift inward, eventually reach the inner MRI front,and accumulate there. They form a dense, thin sub-disk of solid particles, which undergoes gravitational instability to form rocky planetesimals. The planetesimals rapidly grow into a planet through pebble accretion. Consequently, rocky planets tend to be much larger than planets formed through other regimes（tandem planet formation regime and dispersed planet formation regime）, in which icy planet formation actively takes place. These rocky planets may evolve into hot Jupiters if they grow fast enough to the critical core mass of the runaway gas accretion before the dispersal of the disk gas, or they may evolve into super-Earths if the gas dispersed sufficiently early.展开更多
Using a sample of over 70 000 red clump（RC） stars with 5%–10% distance accuracy selected from the LAMOST Spectroscopic Survey of the Galactic Anti-center（LSS-GAC）, we study the radial and vertical gradients of th...Using a sample of over 70 000 red clump（RC） stars with 5%–10% distance accuracy selected from the LAMOST Spectroscopic Survey of the Galactic Anti-center（LSS-GAC）, we study the radial and vertical gradients of the Galactic disk（s） mainly in the anti-center direction, covering a significant volume of the disk in the range of projected Galactocentric radius 7 ≤ RGC ≤ 14 kpc and height from the Galactic midplane 0 ≤ ｜Z ｜ ≤ 3 kpc. Our analysis shows that both the radial and vertical metallicity gradients are negative across much of the volume of the disk that is probed, and they exhibit significant spatial variations. Near the solar circle（7 ≤ RGC ≤ 11.5 kpc）, the radial gradient has a moderately steep, negative slope of-0.08 dex kpc-1near the midplane（｜Z ｜ 〈 0.1 kpc）, and the slope flattens with increasing ｜Z ｜. In the outer disk（11.5 〈 RGC ≤ 14 kpc）, the radial gradients have an essentially constant, much less steep slope of-0.01 dex kpc-1at all heights above the plane, suggesting that the outer disk may have experienced an evolutionary path different from that of the inner disk. The vertical gradients are found to flatten largely with increasing RGC. However, the vertical gradient of the lower disk（0 ≤ ｜Z ｜ ≤ 1 kpc）is found to flatten with RGC quicker than that of the upper disk（1 〈 ｜Z ｜ ≤ 3 kpc）.Our results should provide strong constraints on the theory of disk formation and evolution, as well as the underlying physical processes that shape the disk（e.g. gas flows,radial migration, and internal and external perturbations）.展开更多
More than a decade of dedicated experimental work on the collisional physics of protoplanetary dust has brought us to a point at which the growth of dust aggregates can - for the first time - be self-consistently and ...More than a decade of dedicated experimental work on the collisional physics of protoplanetary dust has brought us to a point at which the growth of dust aggregates can - for the first time - be self-consistently and reliably modeled. In this article, the emergent collision model for protoplanetery dust aggregates, as well as the numerical model for the evolution of dust aggregates in protoplanetary disks, is reviewed. It turns out that, after a brief period of rapid collisional growth of fluffy dust aggregates to sizes of a few centimeters, the protoplanetary dust particles are subject to bouncing collisions, in which their porosity is considerably decreased. The model results also show that low-velocity fragmentation can reduce the final mass of the dust aggregates but that it does not trigger a new growth mode as discussed previously. According to the current stage of our model, the direct formation of kilometer-sized planetesimals by collisional sticking seems unlikely, implying that collective effects, such as the streaming instability and the gravitational instability in dust-enhanced regions of the protoplanetary disk, are the best candidates for the processes leading to planetesimals.展开更多
The purpose of this paper is to explore the influences of cooling timescale on fragmentation of self-gravitating protoplanetary disks. We assume the cooling timescale, expressed in terms of the dynamical timescale Ω ...The purpose of this paper is to explore the influences of cooling timescale on fragmentation of self-gravitating protoplanetary disks. We assume the cooling timescale, expressed in terms of the dynamical timescale Ω tcool, has a power-law dependence on temperature and density, Ω toool ∝∑-aT-b, where a and b are con- stants. We use this cooling timescale in a simple prescription for the cooling rate, du/dt = -u/tcoll, where u is the internal energy. We perform our simulations using the smoothed particle hydrodynamics method. The simulations demonstrate that the disk is very sensitive to the cooling timescale, which depends on density and tem- perature. Under such a cooling timescale, the disk becomes gravitationally unstable and clumps form in the disk. This property even occurs for cooling timescales which are much longer than the critical cooling timescale, Ω toool≥ 7. We show that by adding the dependence of a cooling timescale on temperature and density, the number of clumps increases and the clumps can also form at smaller radii. The simulations im- ply that the sensitivity of a cooling timescale to density is more than to temperature, because even for a small dependence of the cooling timescale on density, clumps can still form in the disk. However, when the cooling timescale has a large dependence on temperature, clumps form in the disk. We also consider the effects of artificial viscos- ity parameters on fragmentation conditions. This consideration is performed in two cases, where Ω tcool is a constant and Ω tcool is a function of density and temperature. The simulations consider both cases, and results show the artificial viscosity param- eters have rather similar effects. For example, using too small of values for linear and quadratic terms in artificial viscosity can suppress the gravitational instability and consequently the efficiency of the clump formation process decreases. This property is consistent with recent simulations of self-gravitating disks. We perform simulations with and without the Balsara form of artificial viscosity. We find that in the cooling and self-gravitating disks without the Balsara switch, the clumps can form more easily than those with the Balsara switch. Moreover, in both cases where the Balsara switch is present or absent, the simulations show that the cooling timescale strongly depends on density and temperature.展开更多
The effects of viscosity on the circumplanetary disks residing in the vicinity of protoplanets are investigated through two-dimensional hydrodynamical simulations with the shearing sheet model. We find that viscosity ...The effects of viscosity on the circumplanetary disks residing in the vicinity of protoplanets are investigated through two-dimensional hydrodynamical simulations with the shearing sheet model. We find that viscosity can considerably affect properties of the circumplanetary disk when the mass of the protoplanet Mp ~ 33 Me, where Me is the Earth＇s mass. However, effects of viscosity on the circumplanetary disk are negligibly small when the mass of the protoplanet Mp 〉 33 Me. We find that when Mp ~ 33 Me, viscosity can markedly disrupt the spiral structure of the gas around the planet and smoothly distribute the gas, which weakens the torques exerted on the protoplanet. Thus, viscosity can slow the migration speed of a protoplanet. After including viscosity, the size of the circumplanetary disk can be decreased by a factor of 〉~ 20%. Viscosity helps to transport gas into the circumplanetary disk from the differentially rotating circumstellar disk. The mass of the circumplanetary disk can be increased by a factor of 50% after viscosity is taken into account when Mp ~ 33 Me. Effects of viscosity on the formation of planets and satellites are briefly discussed.展开更多
We investigate the effects of the cooling function in the formation of clumps of protoplanetary disks using two-dimensional smoothed particle hydrody- namic simulations. We use a simple prescription for the cooling ra...We investigate the effects of the cooling function in the formation of clumps of protoplanetary disks using two-dimensional smoothed particle hydrody- namic simulations. We use a simple prescription for the cooling rate of the flow, du/dt = -u/τcool, where u and %ool are the internal energy and cooling timeseale, respectively. We assume the ratio of local＇cooling to dynamical timescale, Ωτcool =β, to be a constant and also a function of the local temperature. We found that for the constantβ and γ = 5/3, fragmentation occurs only forβ ≤ 7. However, in the case ofβ having temperature dependence and γ = 5/3, fragmentation can also occur for larger values ofβ. By increasing the temperature dependence of the cooling timescale, the mass accretion rate decreases, the population of clumps/fragments increases, and the clumps/fragments can also form in the smaller radii. Moreover, we found that the clumps can form even in a low mass accretion rate, ≤10-7M⊙yr-1, in the case of temperature-dependentβ. However, clumps form with a larger mass accretion rate, 〉 10-7M⊙ yr-1, in the case of constantβ.展开更多
The effect of self-gravity on protoplanetary disks is investigated.The mechanisms of angular momentum transport and energy dissipation are assumed to be the viscosity due to turbulence in the accretion disk.The energy...The effect of self-gravity on protoplanetary disks is investigated.The mechanisms of angular momentum transport and energy dissipation are assumed to be the viscosity due to turbulence in the accretion disk.The energy equation is considered in a situation where the released energy by viscosity dissipation is balanced with cooling processes.The viscosity is obtained by equality of dissipation and cooling functions,and is used to derive the angular momentum equation.The cooling rate of the flow is calculated by a prescription,du/dt = u/τ cool,where u and τ cool are the internal energy and cooling timescale,respectively.The ratio of local cooling to dynamical timescales Ωτ cool is assumed to be a constant and also a function of the local temperature.The solutions for protoplanetary disks show that in the case of Ωτ cool = constant,the disk does not exhibit any gravitational instability over small radii for a typical mass accretion rate,˙ M = 10 6 M yr 1,but when choosing Ωτ cool to be a function of temperature,gravitational instability can occur for this value of mass accretion rate or even less in small radii.Also,by studying the viscosity parameter α,we find that the strength of turbulence in the inner part of self-gravitating protoplanetary disks is very low.These results are qualitatively consistent with direct numerical simulations of protoplanetary disks.Also,in the case of cooling with temperature dependence,the effect of physical parameters on the structure of the disk is investigated.These solutions demonstrate that disk thickness and the Toomre parameter decrease by adding the ratio of disk mass to central object mass.However,the disk thickness and the Toomre parameter increase by adding mass accretion rate.Furthermore,for typical input parameters such as mass accretion rate 10 6 M yr 1,the ratio of the specific heat γ = 5/3 and the ratio of disk mass to central object mass q = 0.1,gravitational instability can occur over the whole radius of the disk excluding the region very near the central object.展开更多
We analyze the radial and vertical metallicity and [α/Fe] gradients of the disk stars of a disk galaxy simulated in a fully cosmological setting with the chemo- dynamical galaxy evolution code GCD＋. We study how the...We analyze the radial and vertical metallicity and [α/Fe] gradients of the disk stars of a disk galaxy simulated in a fully cosmological setting with the chemo- dynamical galaxy evolution code GCD＋. We study how the radial abundance gradients vary as a function of height above the plane and find that the metallicity （[α/Fe]） gra- dient becomes more positive （negative） with increasing height, changing sign around 1.5 kpc above the plane. At the largest vertical height （2 〈 丨z丨 〈 3 kpc）, our simulated galaxy shows a positive radial metallicity gradient. We find that the positive metallicity gradient is caused by the age-metallicity and age-velocity dispersion relation, where younger stars have higher metallicity and lower velocity dispersion. Due to the age- velocity dispersion relation, a greater fraction of younger stars reaches 丨z丨 〉 2 kpc at the outer region, because of the lower gravitational restoring force of the disk, i.e. flaring. As a result, the fraction of younger stars with higher metallicity due to the age-metallicity relation becomes higher at the outer radii, which makes the median metallicity higher at the outer radii. Combining this result with the recently observed age-metallicity and age-velocity dispersion relation for the Milky Way thick disk stars suggested by Haywood et al., we argue that the observed （small） positive radial metal- licity gradient at large heights of the Milky Way disk stars can be explained by flaring of the younger thick and/or thin disk stars.展开更多
基金supported by the National Key Basic Research Program of China (2014CB845700)supported by the National Natural Science Foundation of China (Grant No.11473001)B.Q.C acknowledges partial funding from the China Postdoctoral Science Foundation (2014M560843)
文摘Accurate measurements of stellar metallicity gradients in the radial and vertical directions of the disk and their temporal variations provide important constraints on the formation and evolution of the Milky Way disk. We use 297 042 main sequence turn-off stars selected from the LAMOST Spectroscopic Survey of the Galactic Anticenter（LSS-GAC） to determine the radial and vertical gradients of stellar metallicity,△[Fe/H]/△R and △[Fe/H]/△｜Z ｜ of the Milky Way disk in the direction of the anticenter. We determine ages of those turn-off stars by isochrone fitting and measure the temporal variations of metallicity gradients. We have carried out a detailed analysis of the selection effects resulting from the selection, observation and data reduction of LSS-GAC targets and the potential biases of a magnitude limited sample on the determinations of metallicity gradients. Our results show that the gradients, both in the radial and vertical directions, exhibit significant spatial and temporal variations. The radial gradients yielded by stars with the oldest ages（ 11 Gyr） are essentially zero at all heights from the disk midplane, while those given by younger stars are always negative. The vertical gradients deduced from stars with the oldest ages（ 11 Gyr）are negative and only show very weak variations with Galactocentric distance in the disk plane, R, while those yielded by younger stars show strong variations with R.After being essentially flat at the earliest epochs of disk formation, the radial gradients steepen as age decreases, reaching a maximum（steepest） at age 7–8 Gyr, and then they flatten again. Similar temporal trends are also found for the vertical gradients. We infer that the assembly of the Milky Way disk may have experienced at least two distinct phases. The earlier phase is probably related to a slow, pressure-supported collapse of gas, when the gas settles down to the disk mainly in the vertical direction. In the later phase, there are significant radial flows of gas in the disk, and the rate of gas inflow near the solar neighborhood reaches a maximum around a lookback time of 7–8 Gyr.The transition between the two phases occurs around a lookback time between 8 and11 Gyr. The two phases may be responsible for the formation of the Milky Way＇s thick and thin disks, respectively. Also, as a consequence, we recommend that stellar age is a natural, physical criterion to distinguish stars from the thin and thick disks. From an epoch earlier than 11 Gyr to one between 8 and 11 Gyr, there is an abrupt, significant change in magnitude of both the radial and vertical metallicity gradients, suggesting that stellar radial migration is unlikely to play an important role in the formation of the thick disk.
基金supported by a Grant-in-Aid for Scientific Research on Innovative Areas(Grant No.26106006)
文摘In the tandem planet formation regime,planets form at two distinct sites where solid particles are densely accumulated due to the on/off state of the magnetorotational instability（MRI）.We found that tandem planet formation can reproduce the solid component distribution of the Solar System and tends to produce a smaller number of large planets through continuous pebble flow into the planet formation sites.In the present paper,we investigate the dependence of tandem planet formation on the vertical magnetic field of the protoplanetary disk.We calculated two cases of BZ 3.4 × 10^-3 G and BZ = 3.4 × 10^-5 G at 100 AU as well as the canonical case of BZ = 3.4 × 10^-4 G.We found that tandem planet formation holds up well in the case of the strong magnetic field（BZ 3.4 × 10^-3 G）.On the other hand,in the case of a weak magnetic field（BZ= 3.4 × 10^-5 G） at 100 AU,a new regime of planetary growth is realized：the planets grow independently at different places in the dispersed area of the MRl-suppressed region of r-8-30 AU at a lower accretion rate of M 〈 10^-7.4M⊙yr^-1.We call this the ＂dispersed planet formation＂ regime.This may lead to a system with a larger number of smaller planets that gain high eccentricity through mutual collisions.
基金the Abdul Hamed Shoman Foundation (Grant No. 6/2017) for supporting this projectsupported by the JSPS KAKENHI (Grant No. 18K03706)+1 种基金supported by the National Key R&D Program of China (2016YFA0400801)the National Natural Science Foundation of China (Grant No. U1838201)
文摘We consider the existence of a neutron star magnetic field by the detected cyclotron lines. We collected data on nine sources of high-mass X-ray binaries with supergiant companions as a test case for our model, to demonstrate their distribution and evolution. The wind velocity, spin period and magnetic field strength are studied under different mass loss rates. In our model, correlations between mass-loss rate and wind velocity are found and can be tested in further observations. We examine the parameter space where wind accretion is allowed, avoiding the barrier of rotating magnetic fields, with robust data on the magnetic field of neutron stars. Our model shows that most sources(six of nine systems) can be fed by the wind with relatively slow velocity, and this result is consistent with previous predictions. In a few sources,our model cannot fit the standard wind accretion scenario. In these peculiar cases, other scenarios(disk formation, partial Roche lobe overflow) should be considered. This would provide information about the evolutionary tracks of various types of binaries, and thus exhibit a clear dichotomy behavior in wind-fed X-ray binary systems.
基金Supported by the National Natural Science Foundation of China
文摘We present the results of our recent study on the interactions between a giant planet and a self-gravitating gas disk. We investigate how the disk＇s self-gravity affects the gap formation process and the migration of the giant planet. Two series of 1-D and 2-D hydrodynamic simulations are performed. We select several surface densities and focus on the gravitationally stable region. To obtain more reliable gravity torques exerted on the planet, a refined treatment of the disk＇s gravity is adopted in the vicinity of the planet. Our results indicate that the net effect of the disk＇s self- gravity on the gap formation process depends on the surface density of the disk. We notice that there are two critical values, ∑I and ∑n. When the surface density of the disk is lower than the first one,∑0 〈 ∑I, the effect of self-gravity suppresses the formation of a gap. When ∑0 〉 ∑I, the self-gravity of the gas tends to benefit the gap formation process and enlarges the width/depth of the gap. According to our 1-D and 2-D simulations, we estimate the first critical surface density to be ∑I ≈ 0.8 MMSN. This effect increases until the surface density reaches the second critical value ∑n- When ∑0 〉 ∑n, the gravitational turbulence in the disk becomes dominant and the gap formation process is suppressed again. Our 2-D simulations show that this critical surface density is around 3.5 MMSN. We also study the associated orbital evolution of a giant planet. Under the effect of the disk＇s self-gravity, the migration rate of the giant planet increases when the disk is dominated by gravitational turbulence. We show that the migration timescale correlates with the effective viscosity and can be up to 104 yr.
基金supported by the National Natural Science Foundation of China (Grant No. 11504150)
文摘We investigate the behavior of the snowline in a protoplanetary disk and the relationship between the radius of the snowline and properties of molecular cloud cores.In our disk model,we consider mass influx from the gravitational collapse of a molecular cloud core,irradiation from the central star,and thermal radiation from the ambient molecular cloud gas.As the protoplanetary disk evolves,the radius of the snowline increases first to a maximum value Rmax,and then decreases in the late stage of evolution of the protoplanetary disk.The value of Rmaxis dependent on the properties of molecular cloud cores(mass M;,angular velocity ω and temperature T;).Many previous works found that solid material tends to accumulate at the location of the snowline,which suggests that the snowline is the preferred location for giant planet formation.With these conclusions,we compare the values of R;with semimajor axes of giant planets in extrasolar systems,and find that Rmaxmay provide an upper limit for the locations of the formation of giant planets which are formed by the core accretion model.
基金supported by Grant-in-Aid for Scientific Research on Innovative Areas Number 26106006
文摘We present a new united theory of planet formation,which includes magneto-rotational instability（MRl） and porous aggregation of solid particles in a consistent way.We show that the ＂tandem planet formation＂ regime is likely to result in solar system-like planetary systems.In the tandem planet formation regime,planetesimals form at two distinct sites：the outer and inner edges of the MRl suppressed region.The former is likely to be the source of the outer gas giants,and the latter is the source for the inner volatile-free rocky planets.Our study spans disks with a various range of accretion rates,and we find that tandem planet formation can occur for M = 10^7.3- 10^-6.9Myr^-1.The rocky planets form between 0.4-2 AU,while the icy planets form between 6-30 All;no planets form in 2-6 AU region for any accretion rate.This is consistent with the gap in the solid component distribution in the solar system,which has only a relatively small Mars and a very small amount of material in the main asteroid belt from 2-6 AU.The tandem regime is consistent with the idea that the Earth was initially formed as a completely volatile-free planet.Water and other volatile elements came later through the accretion of icy material by occasional inward scattering from the outer regions.Reactions between reductive minerals,such as schreibersite（Fe-jP）,and water are essential to supply energy and nutrients for primitive life on Earth.
基金partly supported by Grant-in-Aid for Scientific Research on Innovative Areas Number 26106006
文摘We studied the particle growth in a protoplanetary disk in a high-ionization environment and found that icy planet formation is inactive for a disk with an ionization rate 100 times higher than that of the present Solar System. In particular, in the case of M 〈 10^（-7.4）M_☉yr^（-1）, only rocky planet formation occurs. In such a case, all the solid materials in the disk drift inward, eventually reach the inner MRI front,and accumulate there. They form a dense, thin sub-disk of solid particles, which undergoes gravitational instability to form rocky planetesimals. The planetesimals rapidly grow into a planet through pebble accretion. Consequently, rocky planets tend to be much larger than planets formed through other regimes（tandem planet formation regime and dispersed planet formation regime）, in which icy planet formation actively takes place. These rocky planets may evolve into hot Jupiters if they grow fast enough to the critical core mass of the runaway gas accretion before the dispersal of the disk gas, or they may evolve into super-Earths if the gas dispersed sufficiently early.
基金supported by the National Key Basic Research Program of China (2014CB845700)the National Natural Science Foundation of China (Grant No. 11473001)
文摘Using a sample of over 70 000 red clump（RC） stars with 5%–10% distance accuracy selected from the LAMOST Spectroscopic Survey of the Galactic Anti-center（LSS-GAC）, we study the radial and vertical gradients of the Galactic disk（s） mainly in the anti-center direction, covering a significant volume of the disk in the range of projected Galactocentric radius 7 ≤ RGC ≤ 14 kpc and height from the Galactic midplane 0 ≤ ｜Z ｜ ≤ 3 kpc. Our analysis shows that both the radial and vertical metallicity gradients are negative across much of the volume of the disk that is probed, and they exhibit significant spatial variations. Near the solar circle（7 ≤ RGC ≤ 11.5 kpc）, the radial gradient has a moderately steep, negative slope of-0.08 dex kpc-1near the midplane（｜Z ｜ 〈 0.1 kpc）, and the slope flattens with increasing ｜Z ｜. In the outer disk（11.5 〈 RGC ≤ 14 kpc）, the radial gradients have an essentially constant, much less steep slope of-0.01 dex kpc-1at all heights above the plane, suggesting that the outer disk may have experienced an evolutionary path different from that of the inner disk. The vertical gradients are found to flatten largely with increasing RGC. However, the vertical gradient of the lower disk（0 ≤ ｜Z ｜ ≤ 1 kpc）is found to flatten with RGC quicker than that of the upper disk（1 〈 ｜Z ｜ ≤ 3 kpc）.Our results should provide strong constraints on the theory of disk formation and evolution, as well as the underlying physical processes that shape the disk（e.g. gas flows,radial migration, and internal and external perturbations）.
基金funded by the German Space Agency (DLR) under grant Nos. 50WM0336, 50WM0636 and 50WM0936the Deutsche Forschungsgemeinschaft (DFG) under grant No. Bl298/7-1
文摘More than a decade of dedicated experimental work on the collisional physics of protoplanetary dust has brought us to a point at which the growth of dust aggregates can - for the first time - be self-consistently and reliably modeled. In this article, the emergent collision model for protoplanetery dust aggregates, as well as the numerical model for the evolution of dust aggregates in protoplanetary disks, is reviewed. It turns out that, after a brief period of rapid collisional growth of fluffy dust aggregates to sizes of a few centimeters, the protoplanetary dust particles are subject to bouncing collisions, in which their porosity is considerably decreased. The model results also show that low-velocity fragmentation can reduce the final mass of the dust aggregates but that it does not trigger a new growth mode as discussed previously. According to the current stage of our model, the direct formation of kilometer-sized planetesimals by collisional sticking seems unlikely, implying that collective effects, such as the streaming instability and the gravitational instability in dust-enhanced regions of the protoplanetary disk, are the best candidates for the processes leading to planetesimals.
基金Financial support from the research council of Damghan University with grant number 91/phys/108/204
文摘The purpose of this paper is to explore the influences of cooling timescale on fragmentation of self-gravitating protoplanetary disks. We assume the cooling timescale, expressed in terms of the dynamical timescale Ω tcool, has a power-law dependence on temperature and density, Ω toool ∝∑-aT-b, where a and b are con- stants. We use this cooling timescale in a simple prescription for the cooling rate, du/dt = -u/tcoll, where u is the internal energy. We perform our simulations using the smoothed particle hydrodynamics method. The simulations demonstrate that the disk is very sensitive to the cooling timescale, which depends on density and tem- perature. Under such a cooling timescale, the disk becomes gravitationally unstable and clumps form in the disk. This property even occurs for cooling timescales which are much longer than the critical cooling timescale, Ω toool≥ 7. We show that by adding the dependence of a cooling timescale on temperature and density, the number of clumps increases and the clumps can also form at smaller radii. The simulations im- ply that the sensitivity of a cooling timescale to density is more than to temperature, because even for a small dependence of the cooling timescale on density, clumps can still form in the disk. However, when the cooling timescale has a large dependence on temperature, clumps form in the disk. We also consider the effects of artificial viscos- ity parameters on fragmentation conditions. This consideration is performed in two cases, where Ω tcool is a constant and Ω tcool is a function of density and temperature. The simulations consider both cases, and results show the artificial viscosity param- eters have rather similar effects. For example, using too small of values for linear and quadratic terms in artificial viscosity can suppress the gravitational instability and consequently the efficiency of the clump formation process decreases. This property is consistent with recent simulations of self-gravitating disks. We perform simulations with and without the Balsara form of artificial viscosity. We find that in the cooling and self-gravitating disks without the Balsara switch, the clumps can form more easily than those with the Balsara switch. Moreover, in both cases where the Balsara switch is present or absent, the simulations show that the cooling timescale strongly depends on density and temperature.
基金Supported by the National Natural Science Foundation of Chinasupported in part by the Natural Science Foundation of China(Grant Nos.10833002,10825314,11103059,11121062 and 11133005)+1 种基金the National Basic Research Program of China(973 Program,2009CB824800)the CAS/SAFEA International Partnership Program for Creative Research Teams
文摘The effects of viscosity on the circumplanetary disks residing in the vicinity of protoplanets are investigated through two-dimensional hydrodynamical simulations with the shearing sheet model. We find that viscosity can considerably affect properties of the circumplanetary disk when the mass of the protoplanet Mp ~ 33 Me, where Me is the Earth＇s mass. However, effects of viscosity on the circumplanetary disk are negligibly small when the mass of the protoplanet Mp 〉 33 Me. We find that when Mp ~ 33 Me, viscosity can markedly disrupt the spiral structure of the gas around the planet and smoothly distribute the gas, which weakens the torques exerted on the protoplanet. Thus, viscosity can slow the migration speed of a protoplanet. After including viscosity, the size of the circumplanetary disk can be decreased by a factor of 〉~ 20%. Viscosity helps to transport gas into the circumplanetary disk from the differentially rotating circumstellar disk. The mass of the circumplanetary disk can be increased by a factor of 50% after viscosity is taken into account when Mp ~ 33 Me. Effects of viscosity on the formation of planets and satellites are briefly discussed.
文摘We investigate the effects of the cooling function in the formation of clumps of protoplanetary disks using two-dimensional smoothed particle hydrody- namic simulations. We use a simple prescription for the cooling rate of the flow, du/dt = -u/τcool, where u and %ool are the internal energy and cooling timeseale, respectively. We assume the ratio of local＇cooling to dynamical timescale, Ωτcool =β, to be a constant and also a function of the local temperature. We found that for the constantβ and γ = 5/3, fragmentation occurs only forβ ≤ 7. However, in the case ofβ having temperature dependence and γ = 5/3, fragmentation can also occur for larger values ofβ. By increasing the temperature dependence of the cooling timescale, the mass accretion rate decreases, the population of clumps/fragments increases, and the clumps/fragments can also form in the smaller radii. Moreover, we found that the clumps can form even in a low mass accretion rate, ≤10-7M⊙yr-1, in the case of temperature-dependentβ. However, clumps form with a larger mass accretion rate, 〉 10-7M⊙ yr-1, in the case of constantβ.
文摘The effect of self-gravity on protoplanetary disks is investigated.The mechanisms of angular momentum transport and energy dissipation are assumed to be the viscosity due to turbulence in the accretion disk.The energy equation is considered in a situation where the released energy by viscosity dissipation is balanced with cooling processes.The viscosity is obtained by equality of dissipation and cooling functions,and is used to derive the angular momentum equation.The cooling rate of the flow is calculated by a prescription,du/dt = u/τ cool,where u and τ cool are the internal energy and cooling timescale,respectively.The ratio of local cooling to dynamical timescales Ωτ cool is assumed to be a constant and also a function of the local temperature.The solutions for protoplanetary disks show that in the case of Ωτ cool = constant,the disk does not exhibit any gravitational instability over small radii for a typical mass accretion rate,˙ M = 10 6 M yr 1,but when choosing Ωτ cool to be a function of temperature,gravitational instability can occur for this value of mass accretion rate or even less in small radii.Also,by studying the viscosity parameter α,we find that the strength of turbulence in the inner part of self-gravitating protoplanetary disks is very low.These results are qualitatively consistent with direct numerical simulations of protoplanetary disks.Also,in the case of cooling with temperature dependence,the effect of physical parameters on the structure of the disk is investigated.These solutions demonstrate that disk thickness and the Toomre parameter decrease by adding the ratio of disk mass to central object mass.However,the disk thickness and the Toomre parameter increase by adding mass accretion rate.Furthermore,for typical input parameters such as mass accretion rate 10 6 M yr 1,the ratio of the specific heat γ = 5/3 and the ratio of disk mass to central object mass q = 0.1,gravitational instability can occur over the whole radius of the disk excluding the region very near the central object.
基金Supported by the National Natural Science Foundation of China
文摘We analyze the radial and vertical metallicity and [α/Fe] gradients of the disk stars of a disk galaxy simulated in a fully cosmological setting with the chemo- dynamical galaxy evolution code GCD＋. We study how the radial abundance gradients vary as a function of height above the plane and find that the metallicity （[α/Fe]） gra- dient becomes more positive （negative） with increasing height, changing sign around 1.5 kpc above the plane. At the largest vertical height （2 〈 丨z丨 〈 3 kpc）, our simulated galaxy shows a positive radial metallicity gradient. We find that the positive metallicity gradient is caused by the age-metallicity and age-velocity dispersion relation, where younger stars have higher metallicity and lower velocity dispersion. Due to the age- velocity dispersion relation, a greater fraction of younger stars reaches 丨z丨 〉 2 kpc at the outer region, because of the lower gravitational restoring force of the disk, i.e. flaring. As a result, the fraction of younger stars with higher metallicity due to the age-metallicity relation becomes higher at the outer radii, which makes the median metallicity higher at the outer radii. Combining this result with the recently observed age-metallicity and age-velocity dispersion relation for the Milky Way thick disk stars suggested by Haywood et al., we argue that the observed （small） positive radial metal- licity gradient at large heights of the Milky Way disk stars can be explained by flaring of the younger thick and/or thin disk stars.
