Search

We present new stellar evolutionary sequences of very metal-rich stars evolved with the Monash Stellar Structure code and with MESA. The Monash models include masses of $1-8M_{\odot}$ with metallicities $Z=0.04$ to $Z=0.1$ and are evolved from the main sequence to the thermally-pulsing asymptotic giant branch (AGB). These are the first $Z=0.1$ AGB models in the literature. The MESA models include intermediate-mass models with $Z=0.06$ to $Z=0.09$ evolved to the onset of the thermally-pulsing phase. Third dredge-up only occurs in intermediate-mass models $Z \le 0.08$. Hot bottom burning (HBB) shows a weaker dependence on metallicity, with the minimum mass increasing from 4.5$M_{\odot}$ for $Z=0.014$ to $\approx 5.5 M_{\odot}$ for Z = 0.04, $6M_{\odot}$ for $ 0.05 \le Z \le 0.07$ and above 6.5$ M_{\odot}$ for $Z\ge 0.08$. The behaviour of the $Z=0.1$ models is unusual; most do not experience He-shell instabilities owing to rapid mass-loss on the early part of the AGB. Turning off mass-loss produces He-shell instabilities, however thermal pulses are weak and result in no third dredge-up. The minimum mass for carbon ignition is reduced from 8$M_{\odot}$ for $Z=0.04$ to 7$M_{\odot}$ for $Z=0.1$, which implies a reduction in the minimum mass for core-collapse supernovae. MESA models of similarly high metallicity ($Z=0.06 - 0.09$) show the same lowering of the minimum mass for carbon ignition: carbon burning is detected in a $6 M_{\odot}$ model at the highest metallicity ($Z=0.09$) and in all $7 M_{\odot}$ models with $Z \ge 0.06$. This demonstrates robustness of the lowered carbon burning threshold across codes., 18 pages, 10 figures, accepted for publication in MNRAS

Ellipsoidal variables present light-curve modulations caused by stellar distortion, induced by tidal interaction with their companions. An analytical approximated model of the ellipsoidal modulation is given as a discrete Fourier series by Morris and Naftilan 1993 (MN93). Based on numerical simulations using the PHOEBE code we present here updated amplitudes of the first three harmonics of the model. The expected amplitudes are given as a function of the mass ratio and inclination of the binary system and the fillout factor of the primary---the ratio between the stellar radius and that of its Roche lobe. The corrections can get up to 30% relative to the MN93 model for fillout factors close to unity. The updated model can be instrumental in searching for short-period binaries with compact-object secondaries in large data sets of photometric light curves. As shown in one OGLE light-curve example, the minimum mass ratio can be obtained by using only the amplitudes of the three harmonics and an estimation of the stellar temperature. High enough amplitudes can help to identify binaries with mass ratios larger than unity, some of which might have compact companions., Comment: 13 pages, 6 figures

In protoplanetary discs (PPDs) consisting of gas and dust particles, fluid instabilities induced by the drag force, including secular gravitational instability (SGI) can facilitate planet formation. Although SGI subject to the axisymmetric perturbations was originally studied in the absence of gas feedback and it then generalized using a two-fluid approach, the fate of the nonaxisymmetric SGI, in either case, is an unexplored problem. We present a linear perturbation analysis of the nonaxisymmetric SGI in a PPD by implementing a two-fluid model. We explore the growth of the local, nonaxisymmetric perturbations using a set of linearized perturbation equations in a sheared frame. The nonaxisymmetric perturbations display significant growth during a finite time interval even when the system is stable against the axisymmetric perturbations. Furthermore, the surface density perturbations do not show continuous growth but are temporally amplified. We also study cases where the dust component undergoes amplification whereas the gas component remains stable. The amplitude amplification, however, strongly depends on the model parameters. In the minimum mass solar nebula (MMSN), for instance, the dust fluid amplification at the radial distance 100 au occurs when the Stokes number is about unity. But the amplification factor reduces as the dust and gas coupling becomes weaker. Furthermore, perturbations with a larger azimuthal wavelength exhibit a larger amplification factor., accepted for publication in MNRAS

We study the properties of distinct dark matter halos (i.e., those that are not subhalos) that have a final virial mass $M_{\mathrm{vir}}$ at $z = 0$ less than their peak mass ($M_{\mathrm{peak}}$) in the Bolshoi-Planck cosmological simulation. We identify two primary causes of halo mass loss: relaxation after a major merger and tidal stripping by a massive neighbouring halo. Major mergers initially boost $M_{\mathrm{vir}}$ and typically cause the final halo to become more prolate and less relaxed and to have higher spin and lower NFW concentration. As the halo relaxes, high energy material from the recent merger gradually escapes beyond the virial radius, temporarily resulting in a net negative accretion rate that reduces the halo mass by $5-15\%$ on average. Halos that experience a major merger around $z = 0.4$ typically reach a minimum mass near $z = 0$. Tidal stripping mainly occurs in dense regions, and it causes halos to become less prolate and have lower spins and higher NFW concentrations. Tidally stripped halos often lose a large fraction of their peak mass ($> 20\%$) and most never recover (or even reattain a positive accretion rate). Low mass halos can be strongly affected by both post-merger mass loss and tidal stripping, while high mass halos are predominantly influenced by post-merger mass loss and show few signs of significant tidal stripping.

We present new results from the Parallaxes of Southern Extremely Cool dwarfs program to measure parallaxes, proper motions and multi-epoch photometry of L and early T dwarfs. The observations were made on 108 nights over the course of 8 years using the Wide Field Imager on the ESO 2.2m telescope. We present 118 new parallaxes of L \& T dwarfs of which 52 have no published values and 24 of the 66 published values are preliminary estimates from this program. The parallax precision varies from 1.0 to 15.5 mas with a median of 3.8 mas. We find evidence for 2 objects with long term photometric variation and 24 new moving group candidates. We cross-match our sample to published photometric catalogues and find standard magnitudes in up to 16 pass-bands from which we build spectral energy distributions and H-R diagrams. This allows us to confirm the theoretically anticipated minimum in radius between stars and brown dwarfs across the hydrogen burning minimum mass. We find the minimum occurs between L2 and L6 and verify the predicted steep dependence of radius in the hydrogen burning regime and the gentle rise into the degenerate brown dwarf regime. We find a relatively young age of 2 Gyr from the kinematics of our sample., Comment: In press MNRAS

The formation of satellite systems in circum-planetary discs is considered to be similar to the formation of rocky planets in a proto-planetary disc, especially Super-Earths. Thus, it is possible to use systems with large satellites to test formation theories that are also applicable to extrasolar planets. Furthermore, a better understanding of the origin of satellites might yield important information about the environment near the growing planet during the last stages of planet formation. In this work we investigate the formation and migration of the Jovian satellites through N-body simulations. We simulated a massive, static, low viscosity, circum-planetary disc in agreement with the minimum mass sub-nebula model prescriptions for its total mass. In hydrodynamic simulations we found no signs of gaps, therefore type II migration is not expected. Hence, we used analytic prescriptions for type I migration, eccentricity and inclination damping, and performed N-body simulations with damping forces added. Detailed parameter studies showed that the number of final satellites is strong influenced by the initial distribution of embryos, the disc temperature, and the initial gas density profile. For steeper initial density profiles it is possible to form systems with multiple satellites in resonance while a flatter profile favours the formation of satellites close to the region of the Galilean satellites. We show that the formation of massive satellites such as Ganymede and Callisto can be achieved for hotter discs with an aspect ratio of H/r ~ 0.15 for which the ice line was located around 30 R_J., 17 pages, 21 figures

Recently revealed differences in planets around M dwarf vs. solar-type stars could arise from differences in their primordial disks, and surveys of T Tauri stars find a correlation between stellar mass and disk mass. "Minimum" disks have been reconstructed for the Solar System and solar-type stars and here this exercise is performed for M dwarfs using Kepler-detected planets. Distribution of planet mass between current orbits produces a disk with total mass of ~0.009Msun and a power-law profile with index 2.2. Disk reconstruction from the output of a forward model of planet formation indicates that the effect of detection bias on disk profile is slight and that the observed scatter in planet masses and semi-major axes is consistent with a universal disk profile. This nominal M dwarf disk is more centrally concentrated than those inferred around the solar-type stars observed by Kepler, and the mass surface density beyond 0.02 AU is sufficient for in situ accretion of planets as single embryos. The mass of refractory solids within 0.5 AU is 5.6Mearth compared to 4Mearth for solar-type stars, in contrast with the trend with total disk mass. The total solids beyond 0.5 AU is sufficient for the core of at least one giant planet., MNRAS, in press

We aim to understand under which conditions a low mass planet can open a gap in viscous dusty protoplanetary discs. For this purpose, we extend the theory of dust radial drift to include the contribution from the tides of an embedded planet and from the gas viscous forces. From this formalism, we derive i) a grain size-dependent criterion for dust gap opening in discs, ii) an estimate of the location of the outer edge of the dust gap and iii) an estimate of the minimum Stokes number above which low-mass planets are able to carve gaps which appear only in the dust disc. These analytical estimates are particularly helpful to appraise the minimum mass of an hypothetical planet carving gaps in discs observed at long wavelengths and high resolution. We validate the theory against 3D SPH simulations of planet-disc interaction in a broad range of dusty protoplanetary discs. We find a remarkable agreement between the theoretical model and the numerical experiments., Comment: 17 pages, 13 figures. Accepted for publication in MNRAS

Fractal dimensions can be used to characterize the clustering and lacunarities in density distributions. We use generalized fractal dimensions to study the neutral hydrogen distribution (HI) during the epoch of reionization. Using a semi-numeric model of ionized bubbles to generate the HI field, we calculate the fractal dimensions for length scales $\sim 10 h^{-1}$ cMpc. We find that the HI field displays significant multifractal behaviour and is not consistent with homogeneity at these scales when the mass averaged neutral fraction $\bar{x}_{\rm HI}^M \gtrsim 0.5$. This multifractal nature is driven entirely by the shapes and distribution of the ionized regions. The sensitivity of the fractal dimension to the neutral fraction implies that it can be used for constraining reionization history. We find that the fractal dimension is relatively less sensitive to the value of the minimum mass of ionizing haloes when it is in the range $\sim 10^9 - 10^{10} h^{-1} M_{\odot}$. Interestingly, the fractal dimension is very different when the reionization proceeds inside-out compared to when it is outside-in. Thus the multifractal nature of HI density field at high redshifts can be used to study the nature of reionization., 11 pages, 10 figures, accepted for publication in MNRAS

Recent timing analysis reveals that the orbital period of the first discovered accreting millisecond pulsar SAX J1808.4-3658 is increasing at a rate $\dot{P}_{\rm orb}=(3.89\pm0.15)\times 10^{-12}~\rm s\,s^{-1}$, which is at least one order of magnitude higher than the value arising from the conservative mass transfer. An ejection of mass loss rate of $10^{-9}~\rm M_{\odot}{\rm yr}^{-1}$ from the donor star at the inner Lagrangian point during the quiescence state could interpret the observed orbital period derivative. However, it is unknown whether this source can offer such a high mass loss rate. In this work, we attempt to investigate an evolutionary channel towards SAX J1808.4-3658. Once the accretion disk becomes thermally and viscously unstable, the spin-down luminosity of the millisecond pulsar and the X-ray luminosity during outbursts are assumed to evaporate the donor star, and the resulting winds carry away the specific orbital angular momentum at the inner Lagrangian point. Our scenario could yield the observed orbital period, the orbital period derivative, the peak X-ray luminosity during outbursts. Low-mass X-ray binaries with a $1.0~\rm M_{\odot}$ donor star, and an orbital period in the range of 0.8 - 1.5 d, may be the progenitor of SAX J1808.4-3658. Our numerical calculations propose that the current donor star mass is $0.044 ~\rm M_{\odot}$, which is in approximately agreement with the minimum mass of the donor star. In addition, our scenario can also account for the formation of black widows or the diamond planet like PSR J1719-1438., 7 pages, 5 figures, MNRAS in press

We present a new technique that allows us to compute ensemble statistics on a local basis, directly relating halo properties to their local environment. This is achieved by the use of a correlated ensemble in which the Large Scale Structure is common to all realizations while having each an independent halo population. The correlated ensemble can be stacked, effectively increasing the halo number density by an arbitrary factor, thus breaking the fundamental limit in the halo number density given by the haloe mass function. This technique allows us to compute \textit{local ensemble statistics} of the matter/haloe distribution at any position in the simulation box, while removing the intrinsic stochasticity in the halo formation process and directly relating halo properties to their environment. We introduce the \textit{Multum In Parvo} (MIP) correlated ensemble simulation consisting of 220 realizations on a 32 h$^{-1}$ Mpc box with $256^3$ particles each. This is equivalent in terms of effective volume and number of particles to a box of $\sim 193$ h$^{-1}$ Mpc of side with $\sim 1540^3$ particles containing $\sim 5\times 10^6$ haloes with a minimum mass of $3.25 \times 10^9$ h$^{-1}$ M$_{\odot}$. The potential of the technique presented here is illustrated by computing the local ensemble statistics of the halo ellipticity and halo shape-LSS alignment. We show that, while there are general trends in the ellipticity and alignment of haloes with their LSS, there are also significant spatial variations which has important implications for observational studies of galaxy shape and alignment., Published in MNRAS

The majority of discovered exoplanetary systems harbour a new class of planets, bodies typically several times more massive than Earth but orbiting their host stars well inside the orbit of Mercury. The origin of these close-in super-Earths and mini-Neptunes is a major unanswered question in planet formation. Unlike Earth, whose atmosphere contains $, 10 pages, 10 figures. Abstract abridged for submission. Fixed typographical error in Eq. 6 exponent

A comprehensive photometric study and an investigation of the orbital period variation of V53 in the globular cluster M4 are presented. The photometric study reveals that the mass ratio and the contact degree of V53 are q similar to 0.078 and f similar to 69%, respectively. The observed variation in the light curve can be explained by adjusting the spot parameters. V53 belongs to extreme mass ratio (q = 50%) binaries, and its mass ratio is close to the minimum mass ratio predicted by theoretical studies, making it a potential object for studying the evolution of binaries and the formation of blue stragglers and FK Com-type stars. The orbital period of V53 shows a long-term decrease at a rate of dp/dt = 5.89(+/- 0.02) x 10(-8) dyr(-1). This secular period decrease may be caused by the combination of mass transfer from the more massive component to the less massive component and an angular momentum loss via magnetic braking. As this mass transfer and angular momentum loss continues, V53 will ultimately evolve into a single fast-rotation star. By studying the statistics of all the contact binaries in globular clusters that have been analyzed, we found a possible correlation between the contact degree and whether or not a contact binary is a blue straggler. A contact binary is likely to become a blue straggler when its fill-out factor is more than 46.25(+/- 2.05)%. More samples should be introduced to confirm this preliminary result in the future.

We study torques on migrating low-mass planets in locally isothermal discs. Previous work on low-mass planets generally kept the planet on a fixed orbit, after which the torque on the planet was measured. In addition to these static torques, when the planet is allowed to migrate it experiences dynamical torques, which are proportional to the migration rate and whose sign depends on the background vortensity gradient. We show that in discs a few times more massive than the Minimum Mass Solar Nebula, these dynamical torques can have a profound impact on planet migration. Inward migration can be slowed down significantly, and if static torques lead to outward migration, dynamical torques can take over, taking the planet beyond zero-torque lines set by saturation of the corotation torque in a runaway fashion. This means the region in non-isothermal discs where outward migration is possible can be larger than what would be concluded from static torques alone., 14 pages, 13 figures, accepted for publication in MNRAS

We present the results of a multi-wavelength study of circumstellar discs around 44 young stellar objects in the 3 Myr old nearby Chamaeleon I star-forming region. In particular, we explore the far-infrared/submm regime using Herschel fluxes. We show that Herschel fluxes at 160-500$\,\mu$m can be used to derive robust estimates of the disc mass. The median disc mass is 0.005$M_{\odot}$ for a sample of 28 Class IIs and 0.006$M_{\odot}$ for 6 transition disks (TDs). The fraction of objects in Chamaeleon-I with at least the `minimum mass solar nebula' is 2-7%. This is consistent with previously published results for Taurus, IC348, $\rho$ Oph. Diagrams of spectral slopes show the effect of specific evolutionary processes in circumstellar discs. Class II objects show a wide scatter that can be explained by dust settling. We identify a continuous trend from Class II to TDs. Including Herschel fluxes in this type of analysis highlights the diversity of TDs. We find that TDs are not significantly different to Class II discs in terms of far-infrared luminosity, disc mass or degree of dust settling. This indicates that inner dust clearing occurs independently from other evolutionary processes in the discs., Comment: 16 pages, 13 figures. Accepted for publication in MNRAS

We carry out a new analysis of the published radial velocity data for the planet-hosting star HD82943. We include the recent Keck/HIRES measurements as well as the aged but much more numerous CORALIE data. We find that the CORALIE radial velocity measurements are polluted by a systematic annual variation which affected the robustness of many previous results. We show that after purging this variation, the residuals still contain a clear signature of an additional $\sim 1100$ days periodicity. The latter variation leaves significant hints in all three independent radial velocity subsets that we analysed: the CORALIE data, the Keck data acquired prior to a hardware upgrade, and the Keck data taken after the upgrade. We mainly treat this variation as a signature of a third planet in the system, although we cannot rule out other interpretations, such as long-term stellar activity. We find it easy to naturally obtain a stable three-planet radial-velocity fit close to the three-planet mean-motion resonance 1:2:5, with the two main planets (those in the 1:2 resonance) in an aligned apsidal corotation. The dynamical status of the third planet is still uncertain: it may reside in as well as slightly out of the 5:2 resonance. We obtain the value of $\sim 1075$ days for its orbital period and of $\sim 0.3 M_{\rm Jup}$ for its minimum mass, while the eccentric parameters are uncertain., 18 pages, 5 tables, 18 figures; accepted for publication in MNRAS

The Kepler-36 system consists of two planets that are spaced unusually close together, near the 7:6 mean motion resonance. While it is known that mean motion resonances can easily form by convergent migration, Kepler-36 is an extreme case due to the close spacing and the relatively high planet masses of 4 and 8 times that of the Earth. In this paper, we investigate whether such a system can be obtained by interactions with the protoplanetary disc. These discs are thought to be turbulent and exhibit density fluctuations which might originate from the magneto-rotational instability. We adopt a realistic description for stochastic forces due to these density fluctuations and perform both long term hydrodynamical and N-body simulations. Our results show that planets in the Kepler-36 mass range can be naturally assembled into a closely spaced planetary system for a wide range of migration parameters in a turbulent disc similar to the minimum mass solar nebula. The final orbits of our formation scenarios tend to be Lagrange stable, even though large parts of the parameter space are chaotic and unstable., Comment: 13 pages, 8 figures, accepted for publication in MNRAS

Two multifrequency campaigns were carried out on OJ287 in 2005: in April when it was in its pre-outburst state, and in November, during the main 12 yr cycle outburst. The wavelength coverage was from radio to X-rays. In the optical-to-UV range the differential spectrum between the observations has a bremsstrahlung spectral shape, consistent with gas at $3 \times 10^{5}K$ temperature. Our result supports the hydrogen column density of the OJ287 host galaxy of $\sim9.3\times 10^{20} cm^{-2}$, the average value found by Gosh & Soundararajaperumal. The $3 \times 10^{5}K$ bremsstrahlung radiation was predicted in the binary black hole model of OJ287, and it arises from a hot bubble of gas which is torn off the accretion disc by the impact of the secondary. As this radiation is not Doppler boosted, the brightness of the outburst provides an estimate for the mass of the secondary black hole, $\sim1.4\times10^{8}$ solar mass. In order to estimate the mass of the primary black hole, we ask what is the minimum mass ratio in a binary system which allows the stability of the accretion disc. By using particle simulations, we find that the ratio is $\sim1.3\times10^{2}$. This makes the minimum mass of the primary $\sim1.8\times10^{10}$ solar mass, in agreement with the mass determined from the orbit solution, $1.84 \times 10^{10}$ solar mass. With this mass value and the measured K-magnitude of the bulge of the host galaxy of OJ287, the system lies almost exactly on the previously established correlation in the black hole mass vs. K-magnitude diagramme. It supports the extension of this correlation to brighter magnitudes and to more massive black holes than has been done previously.

Using the model for (bottom-up) hierarchical halo growth recently developed by Salvador-Sol\'e et al. (2012), we derive the typical spherically averaged density profile for haloes with several relevant masses in the concordant warm dark matter (\Lambda WDM) cosmology with non-thermal sterile neutrinos of two different masses. The predicted density profiles become flat at small radii, as expected from the effects of the spectrum cutoff. The core cannot be resolved, however, because the non-null particle velocity yields the fragmentation of minimum mass protohaloes in small nodes, which invalidates the model at the corresponding radii.

Wide binaries are tracers of the gravity field of the Galaxy, but their study requires some caution. A large list of common proper motion stars selected from the third Astronomischen Gesellschaft Katalog (AGK3) was monitored with the CORAVEL (for COrrelation RAdial VELocities) spectrovelocimeter, in order to prepare a sample of physical binaries with very wide separations. 66 stars received special attention, since their radial velocities (RV) seemed to be variable. These stars were monitored over several years in order to derive the elements of their spectroscopic orbits. In addition, 10 of them received accurate RV measurements from the SOPHIE spectrograph of the T193 telescope at the Observatory of Haute-Provence. For deriving the orbital elements of double-lined spectroscopic binaries (SB2), a new method was applied, which assumed that the RV of blended measurements are linear combinations of the RV of the components. 13 SB2 orbits were thus calculated. The orbital elements were eventually obtained for 52 spectroscopic binaries (SB), two of them making a triple system. 40 SB received their first orbit and the orbital elements were improved for 10 others. In addition, 11 SB were discovered with very long periods for which the orbital parameters were not found. The median period of the 40 first orbits is 1 yr, and several SB should be resolved or should receive an astrometric orbit in future, providing the masses of the components. In addition, it appeared that HD 153252 has a close companion, which is a candidate brown dwarf with a minimum mass of 50 Jupiter masses. The final selection of wide binaries and the derivation of their statistical properties will be presented in a second paper.

The aim of the present paper is to investigate the spatial structure of a protoplanetary disc whose dynamics is governed by magnetorotational turbulence. We perform a series of local 3D chemo-radiative MHD simulations located at different radii of a disc which is twice as massive as the standard minimum mass solar nebula of Hayashi (1981). The ionisation state of the disc is calculated by including collisional ionisation, stellar X-rays, cosmic rays and the decay of radionuclides as ionisation sources, and by solving a simplified chemical network which includes the effect of the absorption of free charges by {\mu}m-sized dust grains. In the region where the ionisation is too low to assure good coupling between matter and magnetic fields, a non-turbulent central "dead zone" forms, which ranges approximately from a distance of 2 AU to 4 AU from the central star. The approach taken in the present work allows for the first time to derive the global spatial structure of a protoplanetary disc from a set of physically realistic numerical simulations.

We calculate radial migration rates of protoplanets in laminar minimum mass solar nebula discs using three-dimensional self-gravitating radiation hydrodynamical (RHD) models. The protoplanets are free to migrate, whereupon their migration rates are measured. For low mass protoplanets (10-50 M_\oplus) we find increases in the migration timescales of up to an order of magnitude between locally-isothermal and RHD models. In the high-mass regime the migration rates are changed very little. These results are arrived at by calculating migration rates in locally-isothermal models, before sequentially introducing self-gravity, and radiative transfer, allowing us to isolate the effects of the additional physics. We find that using a locally-isothermal equation of state, without self-gravity, we reproduce the migration rates obtained by previous analytic and numerical models. We explore the impact of different protoplanet models, and changes to their assumed radii, upon migration. The introduction of self-gravity gives a slight reduction of the migration rates, whilst the inertial mass problem, which has been proposed for high mass protoplanets with circumplanetary discs, is reproduced. Upon introducing radiative transfer to models of low mass protoplanets (\approx 10 M_\oplus), modelled as small radius accreting point masses, we find outward migration with a rate of approximately twice the analytic inward rate. However, when modelling such a protoplanet in a more realistic manner, with a surface which enables the formation of a deep envelope, this outward migration is not seen.

We consider the formation and migration of protoplanetary embryos in disks around the stars in tight binary systems (separations ~ 20 AU. In such systems, the initial stages of runaway embryo formation are expected to only take place within some critical disk radius a_{crit}, due to the perturbing effect of the binary companions (Thebault et al. 2009). We perform n-body simulations of the evolution of such a population of inner-disk embryos surrounded by an outer-disk of smaller planetesimals. Taking Alpha Centauri-B as our fiducial reference example in which a_{crit} ~ 0.7 AU, and using a Minimum Mass Nebular Model with $\Sigma \propto a^{-3/2}$, we find that within 10^6 yrs (10^7 yrs), systems will on average contain embryos which have migrated out to 0.9 AU (1.2 AU), with the average outer-most body having a mass of 0.2 M_{earth} 0.4 M_{earth}. Changes to increase the surface density of solids or to use a flatter profile both produce increased embryo migration and growth. At a given time, the relative change in semi-major axis of the outer-most embryo in these simulations is found to be essentially independent of a_{crit}, and we note that little further embryo migration takes place beyond 10^7 years. We conclude that the suppression of runaway growth outside a_{crit} does not mean that the habitable zones in such tight binary systems will be devoid of detectable, terrestrial mass planets, even if a_{crit} lies significantly interior to the inner edge of the habitable zone.

When the total angular momentum of a binary system is at a critical (minimum) value, a tidal instability occurs (Darwin’s instability), eventually forcing the stars to merge into a single, rapidly rotating object. The instability sets in at some critical separation which in the case of contact binaries corresponds to a minimum mass ratio depending on dimensionless gyration radius k1. If one considers n = 3 polytrope (fully radiative primary with � 1 = 4/3), k 2 = 0.075 and qmin ≈ 0.085–0.095. There appears to be, however, some W UMa-type binaries with q values very close, if not below these theoretical limits, implying that primary in these systems is probably more centrally condensed. We try to solve the discrepancy between theory and observations by considering rotating polytropes. We show by deriving and solving a modified Lane–Emden equation for n = 3 polytrope that including the effects of rotation does increase the central concentration and could reduce qmin to as low as 0.070–0.074, more consistent with the observed population.

We report on the detection of four extrasolar planets orbiting evolved intermediate-mass stars from a precise Doppler survey of G and K giants at Okayama Astrophysical Observatory. All of the host stars are considered to be formerly early F-type or A-type dwarfs when they were on the main sequence. 14 And (K0 III) is a clump giant with a mass of 2.2 M_solar and has a planet of minimum mass m_2sin i=4.8 M_Jup in a nearly circular orbit with a 186 day period. This is one of the innermost planets around evolved intermediate-mass stars and such planets have only been discovered in clump giants. 81 Cet (G5 III) is a clump giant with 2.4 M_solar hosting a planet of m_2sin i=5.3 M_Jup in a 953 day orbit with an eccentricity of e=0.21. 6 Lyn (K0 IV) is a less evolved subgiant with 1.7 M_solar and has a planet of m_2sin i=2.4 M_Jup in a 899 day orbit with e=0.13. HD 167042 (K1 IV) is also a less evolved star with 1.5 M_solar hosting a planet of m_2sin i=1.6 M_Jup in a 418 day orbit with e=0.10. This planet was independently announced by Johnson et al. (2008, ApJ, 675, 784). All of the host stars have solar or sub-solar metallicity, which supports the lack of metal-rich tendency in planet-harboring giants in contrast to the case of dwarfs., Comment: 23 pages, 6 figures, accepted for publication in PASJ

In this Letter, the radial velocity variability of the chemically peculiar star Epsilon Ursae Majoris ($\epsilon$ UMa) from the sharp cores of the hydrogen lines is investigated. This study is based on the ELODIE archival data obtained at different phases of the rotational cycle. The star exhibits low-amplitude radial velocity variations with a period of P=5.0887 d. The best Keplerian solution yields an eccentricity e=0.503 and a minimum mass ~14.7$M_{\rm Jup}$ on the hypothesis that the rotational axis of $\epsilon$ UMa is perpendicular to the orbital plane. This result indicate that the companion is the brown-dwarf with the projected semi-amplitude variation of the radial velocity $K_{\rm 2}$=135.9 km/sec and the sine of inclination times semi-major axis $a_{2}$sin(i)=0.055 au.

Disks in the 6 Myr old cluster eta Chamaeleontis were searched for emission from hot H2. Around the M3 star ECHAJ0843.3-7905 we detect circumstellar gas orbiting at ~2 AU. If the gas is UV-excited, the ro-vibrational line traces a hot gas layer supported by a disk of mass ~0.03Msolar, similar to the minimum mass solar nebula. Such a gas reservoir at 6 Myr would promote the formation and inwards migration of gas giant planets., Accepted for publication in MNRAS. 9 pages

We investigate the dependence of stellar properties on the opacity limit for fragmentation that is set by the metallicity of a molecular cloud. We compare the results from two large-scale hydrodynamical simulations of star cluster formation, which resolve the fragmentation process down to the opacity limit, the first of which was reported by Bate, Bonnell & Bromm. The initial conditions of the two calculations are identical, but in the new simulation the onset of the opacity limit occurs at a lower gas density, and this is expected to increase the minimum mass of a brown dwarf by a factor of 3 [to ≈ 9 Jupiter masses (M J)]. We find that the lowest mass object is a factor of 3 higher in the low-metallicity calculation, as expected. However, apart from this shift of the low-mass cut-off, the initial mass functions (IMFs) produced by the two calculations are indistinguishable. In particular, the median (characteristic) mass is unchanged. These results add support to the accretion‐ejection model proposed by Bate & Bonnell for the origin of the initial mass function (IMF), which predicts that the characteristic mass should vary in proportion to the mean thermal Jeans mass in the cloud. They also indicate that the form of the IMF above the low-mass cut-off should not display a strong metallicity dependence, assuming that the cooling is dominated by dust and that the overall mean thermal Jeans mass of a molecular cloud does not depend on its metallicity. However, if the mean thermal Jeans mass of a molecular cloud is set by the thermal behaviour of gas during the formation of the cloud, this should lead to an indirect dependence of the characteristic mass of the IMF on metallicity because of the link between the characteristic mass and the mean thermal Jeans mass of the cloud. Ke yw ords: accretion, accretion discs ‐ hydrodynamics ‐ binaries: general ‐ stars: formation ‐ stars: low-mass, brown dwarfs ‐ stars: luminosity function, mass function.

By comparing the PSCz galaxy power spectrum with the results of nested pure dark matter N-body simulations, we try to understand how infrared-selected galaxies populate dark-matter haloes, paying special attention to the method of halo identification in the simulations. We thus test the hypothesis that baryonic physics negligibly affects the distribution of galaxies down to the smallest scales yet observed. We are successful in reproducing the PSCz power spectrum on scales < ~40 h/Mpc, near our resolution limit, by imposing a central density cut-off on simulated haloes, which gives a rough minimum mass and circular velocity of haloes in which PSCz galaxies formed., 12 pages, 16 figures (one added), conforms to version in MNRAS

We report precise Doppler measurements of the stars HD 216437, HD 196050 and HD 160691 obtained with the Anglo-Australian Telescope using the UCLES spectrometer together with an iodine cell as part of the Anglo-Australian Planet Search. Our measurements reveal periodic Keplerian velocity variations that we interpret as evidence for planets in orbit around these solar type stars. HD 216437 has a period of 1294$\pm$250 d, a semi-amplitude of 38$\pm$4 m s$^{-1}$ and of an eccentricity of 0.33$\pm$0.09. The minimum (M sin $i$) mass of the companion is 2.1$\pm$0.3 M$_{\rm JUP}$ and the semi-major axis is 2.4$\pm$0.5 au. HD 196050 has a period of 1288$\pm$230 d, a semi-amplitude of 54$\pm$8 m s$^{-1}$ and an eccentricity of 0.28$\pm$0.15. The minimum mass of the companion is 3.0$\pm$0.5 M$_{\rm JUP}$ and the semi-major axis is 2.3$\pm$0.5 au. We also report further observations of the metal rich planet bearing star HD 160691. Our new solution confirms the previously reported planet and shows a trend indicating a second, longer-period companion. These discoveries add to the growing numbers of midly-eccentric, long-period extra-solar planets around Sun-like stars. As seems to be typical of stars with planets, both stars are metal-rich., Comment: 23 pages incl. 6 tables and 3 figures

We describe a method for predicting detection limits of minority elements in electron energy loss spectroscopy (EELS), and its implementation as a software package that gives quantitative predictions for user-specified materials and experimental conditions. The method is based on modeling entire energy loss spectra, including shot noise as well as instrumental noise, and taking into account all the relevant experimental parameters. We describe the steps involved in modeling the entire spectrum, from the zero loss up to inner shell edges, and pay particular attention to the contributions to the pre-edge background. The predicted spectra are used to evaluate the signal-to-noise ratios (SNRs) for inner shell edges from user-specified minority elements. The software also predicts the minimum detectable mass (MDM) and minimum mass fraction (MMF). It can be used to ascertain whether an element present at a particular concentration should be detectable for given experimental conditions, and also to quickly and quantitatively explore ways of optimizing the experimental conditions for a particular EELS analytical task. We demonstrate the usefulness of the software by confirming the recent empirical observation of single atom detection using EELS of phosphorus in thin carbon films, and show the effect on the SNR of varying the acquisition parameters. The case of delta-doped semiconductors is also considered as an important example from materials science where low detection limits and high spatial resolution are essential, and the feasibility of such characterization using EELS is assessed.

We discuss the evolution and the fragmentation of cylindrical gas clouds with primordial composition, using the virial equation. Comparing the timescale of collapse with that of fragmentation, we obtain the fragment mass scale for the various initial density. We estimate the minimum mass of a fragment that is formed from the primordial gas cloud. The mass scale is essentially determined by the Chandrasekhar mass. The effect of dark matter is also investigated. If the filamentary cloud consists not only of ordinary gas but also some dark matter, the fragmentation of the cylindrical cloud occurs at the lower density than the pure gas case. As a result, the fragment mass scale hardly becomes Chandrasekhar mass scale. In this case, massive star formation is strongly expected.

We have investigated numerically the orbital instability of a protoplanet system while taking account of the gasdrag force due to the solar nebula. In the present work, we considered an equally spaced five-protoplanet (with the same mass of 1 ×10 −7 M� ) system, in which their initial orbits are coplanar and circular, and assumed that the gas-drag force is proportional to the square of the relative velocity between the gas and a protoplanet. We first reexamined and confirmed that, under a gas-free condition, log 10 Tinst can be approximately written as a linear function of the initial orbital separation distance, ∆˜ a0, where Tinst is the time of the orbital instability (i.e., the time of the first orbital crossing between any two protoplanets). Next, we investigated the instability time under the gas-drag effect, T gas inst , and found that T gas inst suddenly becomes large compared with Tinst, when ∆˜ a0 is larger than a certain critical separation distance, (∆˜ a0)crit. Furthermore, we showed that (∆˜ a0)crit can be described semi-analytically as a function of the gaseous density. From a function extrapolated with a density in the minimum mass nebula model, we estimated (∆˜ a0)crit in the nebula as being about 10 Hill radius at 1 AU.

When the total angular momentum of a binary system J tot = J orb + J spin is at a certain critical (minimum) value, a tidal instability occurs which eventually forces the stars to merge into a single, rapidly rotating object. The instability occurs when J orb = 3J spin , which in the case of contact binaries corresponds to a minimum mass ratio q mm ≈ 0.071-0.078. The minimum mass ratio is obtained under the assumption that stellar radii are fixed and independent. This is not the case with contact binaries where, according to the Roche model, we have R 2 = R 2 (R 1, a, q). By finding a new criterion for contact binaries, which arises from dJ tot = 0, and assuming k 2 1 ≠ k 2 2 for the component's dimensionless gyration radii, a theoretical lower limit q min = 0.094-0.109 for overcontact degree f = 0-1 is obtained.

We run very large cosmological N-body hydrodynamical simulations in order to study statistically the baryon fractions in early dark matter haloes. We critically examine how differences in the initial conditions affect the gas fraction in the redshift range z= 11–21. We test three different linear power spectra for the initial conditions. (1) A complete heating model, which is our fiducial model; this model follows the evolution of overdensities correctly, according to Naoz & Barkana (2005), in particular including the spatial variation of the speed of sound of the gas due to Compton heating from the CMB. (2) An equal-δ model, which assumes that the initial baryon fluctuations are equal to those of the dark matter, while conserving σ8 of the total matter. (3) A mean cs model, which assumes a uniform speed of sound of the gas. The latter two models are often used in the literature. We calculate the baryon fractions for a large sample of haloes in our simulations. Our fiducial model implies that before reionization and significant stellar heating took place, the minimum mass needed for a minihalo to keep most of its baryons throughout its formation was ∼3 × 104 M⊙. However, the alternative models yield a wrong (higher by about 50 per cent) minimum mass, since the system retains a memory of the initial conditions. We also demonstrate this using the ‘filtering mass’ from linear theory, which accurately describes the evolution of the baryon fraction throughout the simulated redshift range.

We aim at understanding the massive star formation (MSF) limit m(r) = 870M_☉(r/pc)^(1.33) in the mass–size space of molecular structures recently proposed by Kauffmann & Pillai. As a first step, we build on the property that power-law density profiles for molecular clumps combined with a volume density threshold for the overall star formation naturally leads to mass–radius relations for molecular clumps containing given masses of star-forming gas. Specifically, we show that the mass m_(clump) and radius r_(clump) of molecular clumps whose density profile slope is −p and which contain a mass mth of gas denser than a density threshold ρ_(th) obeys the following: m_(clump) = m^(p/3)_(th) (^(4πpth) _(3-p))^((3-p)/3) r^(3-p)_(clump). In a second step, we use the relation between the mass of embedded clusters and the mass of their most massive star to estimate the minimum mass of the star-forming gas needed to form a 10-M_⊙ star. Assuming a star formation efficiency (SFE) of SFE ≃ 0.30, this gives m_(th,crit)≃150 M_⊙. In a third step, we demonstrate that, for sensible choices of the clump density index (p ≃ 1.7) and of the cluster formation density threshold (n_(th)≃ 10^4 cm^(−3)), the line of constant m_(th,crit)≃150 M_⊙ in the mass–radius space of molecular structures equates to the MSF limit for spatial scales larger than 0.3 pc. Hence, the observationally inferred MSF limit of Kauffmann & Pillai is consistent with a threshold in star-forming gas mass beyond which the star-forming gas reservoir is large enough to allow the formation of massive stars. For radii smaller than 0.3 pc, the MSF limit is shown to be consistent with the formation of a 10-M_⊙ star (m_(th,crit) ≃ 30 M_⊙ with SFE ≃ 0.3) out of its individual pre-stellar core of density threshold n_(th) ≃ 10^5 cm^(-3). The inferred density thresholds for the formation of star clusters and individual stars within star clusters match those previously suggested in the literature.

Using Eggleton's stellar evolution code, we study the minimum mass ratio ($q_{\rm min}$) of W Ursae Majoris (W UMa) binaries that have different primary masses. It is found that the minimum mass ratio of W UMa binaries decreases with increasing mass of the primary if the primary's mass is less than about 1.3$M_{\rm \odot}$, and above this mass the ratio is roughly constant. By comparing the theoretical minimum mass ratio with the observational data, it is found that the existence of low-$q$ systems can be explained by the different structure of the primaries with different masses. This suggests that the dimensionless gyration radius ($k_1^2$) and thus the structure of the primary is very important in determining the minimum mass ratio. In addition, we investigate the mass loss during the merging process of W UMa systems and calculate the rotation velocities of the single stars formed by the merger of W UMa binaries due to tidal instability. It is found that in the case of the conservation of mass and angular momentum, the merged single stars rotate with a equatorial velocity of about $588\sim819$ km s$^{-1}$, which is much larger than their break-up velocities ($v_{\rm b}$). This suggests that the merged stars should extend to a very large radius (3.7$\sim$5.3 times the radii of the primaries) or W UMa systems would lose a large amount of mass (21$\sim$33 per cent of the total mass) during the merging process. If the effect of magnetic braking is considered, the mass loss decreases to be 12$\sim$18 per cent of their total masses. This implies that the significant angular momentum and mass might be lost from W UMa systems in the course of the merging process, and this kind of mass and angular momentum loss might be driven by the release of orbital energy of the secondaries, which is similar to common-envelope evolution.

The lower limit of the mass of the first stars suggested recently may imply the formation of massive stars of mass greater than 8 solar mass irrespective of the details of the initial mass function. The production of heavy metals from the first stars will ensure a requisite for the existence of life without the anthropic principle., Comment: 3 Pages, to be published in Progress of Theoretical Physics 1997 January issue

We report the detections of substellar companions orbiting around seven evolved intermediate-mass stars from precise Doppler measurements at Okayama Astrophysical Observatory. o UMa (G4 II-III) is a giant with a mass of 3.1 M_sun and hosts a planet with minimum mass of m_2sini=4.1 M_J in an orbit with a period P=1630 d and an eccentricity e=0.13. This is the first planet candidate (< 13 M_J) ever discovered around stars more massive than 3 M_sun. o CrB (K0 III) is a 2.1 M_sun giant and has a planet of m_2sini=1.5 M_J in a 187.8 d orbit with e=0.19. This is one of the least massive planets ever discovered around ~2 M_sun stars. HD 5608 (K0 IV) is an 1.6 M_sun subgiant hosting a planet of m_2sini=1.4 M_J in a 793 d orbit with e=0.19. The star also exhibits a linear velocity trend suggesting the existence of an outer, more massive companion. 75 Cet (G3 III:) is a 2.5 M_sun giant hosting a planet of m_2sini=3.0 M_J in a 692 d orbit with e=0.12. The star also shows possible additional periodicity of about 200 d and 1880 d with velocity amplitude of ~7--10 m/s, although these are not significant at this stage. nu Oph (K0 III) is a 3.0 M_sun giant and has two brown-dwarf companions of m_2sini= 24 M_J and 27 M_J, in orbits with P=530.3 d and 3190 d, and e=0.126 and 0.17, respectively, which were independently announced by Quirrenbach et al. (2011). The ratio of the periods is close to 1:6, suggesting that the companions are in mean motion resonance. We also independently confirmed planets around k CrB (K0 III-IV) and HD 210702 (K1 IV), which had been announced by Johnson et al. (2008) and Johnson et al. (2007a), respectively. All of the orbital parameters we obtained are consistent with the previous results., 21 pages, 14 figures, accepted for publication in PASJ

The near contact binary BF Virginis was monitored by us from 2007 to 2011 with the 60-cm telescope of Yunnan Observatory. Eight new observed times of light minimum were derived. All of the available eclipse times, spanning more than 100 yr, were analyzed. It is found that there exists a cyclic variation with a semi-amplitude of 0.034 d and a period of 63.3 yr, superposed by an upward parabolic variation. The upward parabolic change can be explained as a continuous period increase due to mass transfer from the secondary to the primary component of BF Vir, which supports the secondary filling near contact configuration. The period oscillation can be explained either by the light-time effect due to the presence of a third body or by magnetic activity. If it is caused by the light-time effect, the minimum mass and the corresponding orbital radius of the additional body should be m(3) = 0.93 M-circle dot and a(3) = 19.4 AU, respectively. On the other hand, the period oscillation may be explained as a consequence of magnetic cycles in its secondary component. In this case, the required variation of the quadrupole moment is calculated to be Delta Q = 2.1 x 10(50) g cm(2).