The orbital distribution of Near-Earth Objects inside Earth’s orbit

Canada’s Near-Earth Object Surveillance Satellite (NEOSSat), set to launch in early 2012, will search for and track Near-Earth Objects (NEOs), tuning its search to best detect objects with a < 1.0 AU. In order to construct an optimal pointing strategy for NEOSSat, we needed more detailed information in the a < 1.0 AU region than the best current model (Bottke, W.F., Morbidelli, A., Jedicke, R., Petit, J.M., Levison, H.F., Michel, P., Metcalfe, T.S. [2002]. Icarus 156, 399–433) provides. We present here the NEOSSat-1.0 NEO orbital distribution model with larger statistics that permit finer resolution and less uncertainty, especially in the a < 1.0 AU region. We find that Amors = 30.1 ± 0.8%, Apollos = 63.3 ± 0.4%, Atens = 5.0 ± 0.3%, Atiras (0.718 < Q < 0.983 AU) = 1.38 ± 0.04%, and Vatiras (0.307 < Q < 0.718 AU) = 0.22 ± 0.03% of the steady-state NEO population. Vatiras are a previously undiscussed NEO population clearly defined in our integrations, whose orbits lie completely interior to that of Venus. Our integrations also uncovered the unexpected production of retrograde orbits from main-belt asteroid sources; this retrograde NEA population makes up ?0.1% of the steady-state NEO population. The relative NEO impact rate onto Mercury, Venus, and Earth, as well as the normalized distribution of impact speeds, was calculated from the NEOSSat-1.0 orbital model under the assumption of a steady-state. The new model predicts a slightly higher Mercury impact flux.     

Radiative heating of carbonaceous near-Earth objects as a cause of thermal metamorphism for CK chondrites

Metamorphic CK carbonaceous chondrites display matrix textures that are best explained by a transient thermal event with temperatures in the 550–950 K range and durations in the order of days to years, longer than what is commonly admitted for shock events but shorter than what is required for nuclide decay. We propose that radiative heating of small carbonaceous meteoroids with perihelia close to the Sun could account for the petrological features observed in CK chondrites. Numerical thermal modeling, using favorable known NEOs orbital parameters (perihelion distances between 0.07 and 0.15 AU) and physical properties of CV and CK chondrites (albedo in the range 0.01–0.1, 25% porosity, thermal diffusivity of 0.5–1.5 W m^?1 K^?1), shows that radiative heating can heat carbonaceous meteoroids in the meter size range to core temperatures up to 1050 K, consistent with the metamorphic temperatures estimated for CK chondrites. Sizes of known CV and CK chondrites indicate that all these objects were small meteoroids (radii from a few cm to 2.5 m) prior to their atmospheric entry. Simulations of dynamic orbits for NEO objects suggest that there are numerous such bodies with suitable orbits and properties, even if they are only a small percentage of all NEOs. Radiative heating would be a secondary process (superimposed on parent-body processes) affecting meteoroids formed by the disruption of an initially homogeneous CV3-type parent body. Different petrologic types can be accounted for depending on the sizes and heliocentric distances of the objects in such a swarm.     

Searching for superhump period of cataclysmic variable AY Psc

The results obtained from unfiltered photometric CCD observations of AY Psc made during 17, 20 and 12 nights, respectively, in 2003, 2004 and 2005 are presented. A period of 0.21732 ± 0.00001 d was detected in the data. This period is consistent with the previously proposed orbital period of P_orb = 0.2173209 d ((Diaz and Steiner, 1990)). Since this period was present in the light curves taken in all three years, with no apparent change in its value or amplitude, it is interpreted as the orbital period of this binary system. In addition, quasi-periodicities of 0.2057 ± 0.0001 d, 0.2063 ± 0.0001 d, 0.2072 ± 0.0001 d for the years 2003, 2004 and 2005, respectively, were also discovered. These periods were interpreted as negative superhump periods and it was seen that they changes from year to year. Therefore AY Psc is then classified as a negative superhump system.     

Cyclic period changes and the light-time effect in eclipsing binaries: A low-mass companion around the system VV Ursae Majoris

In this article, a period analysis of the late-type eclipsing binary VV UMa is presented. This work is based on the periodic variation of eclipse timings of the VV UMa binary. We determined the orbital properties and mass of a third orbiting body in the system by analyzing the light-travel time effect. The O−C diagram constructed for all available minima times of VV UMa exhibits a cyclic character superimposed on a linear variation. This variation includes three maxima and two minima within approximately 28,240 orbital periods of the system, which can be explained as the light-travel time effect (LITE) because of an unseen third body in a triple system that causes variations of the eclipse arrival times. New parameter values of the light-time travel effect because of the third body were computed with a period of 23.22 ± 0.17 years in the system. The cyclic-variation analysis produces a value of 0.0139 day as the semi-amplitude of the light-travel time effect and 0.35 as the orbital eccentricity of the third body. The mass of the third body that orbits the eclipsing binary stars is 0.787 ± 0.02 M_⊙, and the semi-major axis of its orbit is 10.75 AU.     

CCD Photometry of the W UMa type star QX Andromeda

We present new B- and V-band photometry of the W UMa-type binary system QX And, which is a member of the open cluster NGC 752. Revised orbital period and new ephemerides were given for the binary system based on the data of times of light minima. The result of a period analysis reveals that the system is undergoing a continuous orbital period increase during the past decades. The rate of period increasing turns out to be about 2.7 × 10^?7 d yr^?1. With the Wilson–Devinney code, a photometric solution is computed. It yields a contact configuration for the system with a filling factor of 0.361. Combining the results from the photometric solution along with that from the radial-velocity observations, we have determined the absolute parameters for the two components of the system. The masses, radii and luminosity of the primary and secondary stars are calculated as 1.43 ± 0.04 M_⊙, 1.45 ± 0.09 R_⊙, 2.87 ± 0.40 L_⊙ and 0.44 ± 0.02 M_⊙, 0.87 ± 0.05 R_⊙, 0.99 ± 0.13 L_⊙, respectively. The evolutionary status and physical nature of the contact binary system were discussed compared with the theoretical models.     

The orbital and superhump periods of the dwarf nova HS 0417+7445 in Camelopardalis

We present the 2005–2010 outburst history of the SU UMa-type dwarf HS 0417+7445, along with a detailed analysis of extensive time-series photometry obtained in March 2008 during the second recorded superoutburst of the system. The mean outburst interval is 197 ± 59 d, with a median of 193 d. The March 2008 superoutburst was preceded by a precursor outburst, had an amplitude of 4.2 magnitudes, and the whole event lasted about 16 days. No superhumps were detected during the decline from the precursor outburst, and our data suggests instead that orbital humps were present during that phase. Early superhumps detected during the rise to the superoutburst maximum exhibited an unusually large fractional period excess of ? = 0.137 (P_sh = 0.0856(88) d). Following the maximum, a linear decline in brightness followed, lasting at least 6 days. During this decline, a stable superhump period of P_sh = 0.07824(2) d was measured. Superimposed on the superhumps were orbital humps, which allowed us to accurately measure the orbital period of HS 0417+7445, P_orb = 0.07531(8) d, which was previously only poorly estimated. The fractional superhump period excess during the main phase of the outburst was ? = 0.037, which is typical for SU UMa dwarf novae with similar orbital period. Our observations are consistent with the predictions of the thermal-tidal instability model for the onset of superoutbursts, but a larger number of superoutbursts with extensive time-series photometry during the early phases of the outburst would be needed to reach a definite conclusion on the cause of superoutbursts.     

The near-contact binary star RZ Dra revisited

This paper presents the absolute parameters of RZ Dra. New CCD observations were made at the Mt. Suhora Observatory in 2007. Two photometric data sets (1990 BV and 2007 BVRI) were analysed using modern light-curve synthesis methods. Large asymmetries in the light curves may be explained in terms of a dark starspot on the primary component, an A6 type star. Due to this magnetic activity, the primary component would appear to belong to the class of Ap-stars and would show small amplitude with δ Scuti-type pulsations. With this in mind, a time-series analysis of the residual light curves was made. However, we found no evidence of pulsation behaviour in RZ Dra. Combining the solutions of our light curves and Rucinski et al. (2000)’s radial velocity curves, the following absolute parameters of the components were determined: M₁ = 1.63 ± 0.03 M_⊙, M₂ = 0.70 ± 0.02 M_⊙, R₁ = 1.65 ± 0.02R_⊙, R₂ = 1.15 ± 0.02 R_⊙, L₁ = 9.72 ± 0.30 L_⊙ and L₂ = 0.74 ± 0.10 L_⊙. The distance to RZ Dra was calculated as 400 ± 25 pc, taking into account interstellar extinction. The orbital period of the system was studied using updated O–C information. It was found that the orbital period varied in its long-period sinusoidal form, superimposed on a downward parabola. The parabolic term shows a secular period decrease at a slow rate of 0.06 ± 0.02 s per century and is explained by the mass loss via magnetized wind of the Ap-star primary. The tilted sinusoidal form of the period variation may be considered as an apparent change and may be interpreted in terms of the light-time effect due to the presence of a third body.     

Main Belt Comet P/2008 R1 Garradd: Duration of activity

In this work we investigated changes of the water emission from a model comet of the size and orbital elements of Comet P/2008 R1 (Garradd). We performed simulations for model cometary nuclei of different compositions and two different orientations in space. Our simulations indicate, that the emission of water decreases from one orbital period to another one, but in some cases slowly. When the rotation axis of the nucleus lies in the orbital plane the seasonal maximum of water production during the first two orbital periods can be as high as about 10²⁶ mol s^?1, but decreases by two orders of magnitude during only 50 orbital periods. The highest rate of water production after many orbital periods is expected when the rotation axis is perpendicular to the orbital plane – the seasonal maximum of water production can be about 5 × 10²⁵ mol s^?1 during the first two orbital periods after activation of the comet and no more than 0.8 × 10²⁵ mol s^?1 500 orbital periods later. The upper estimate for the production of water derived from observations of P/2008 R1 (Garradd) by Jewitt et al. (Jewitt, D., Yang, B., Haghighipour, N. [2009]. Astron. J. 137, 4313–4321) is 5 × 10²⁵.     

Absolute properties of the overcontact binary HH Boo

We obtained multi-colour light curves of the overcontact binary system HH Boo and analysed the orbital period variation of the system. Our analysis tentatively indicates either mass transfer from the secondary to the primary or mass loss from the system at a rate of -5.04 × 10⁻⁷ M_⊙ per year. Through a combined analysis of the published radial velocity curve and light curves, we determined an inclination (i) of 69°.71 ± 0°.16 and a semi-major axis (a) of 2.246 ± 0.064 R_⊙ for HH Boo. The masses of the primary and secondary components were found to be 0.92 ± 0.08 M_⊙ and 0.58 ± 0.06 M_⊙, respectively. The radius determined for the primary was 0.98 ± 0.03 R_⊙, while that determined for the secondary was 0.80 ± 0.02 R_⊙. We demonstrated that HH Boo is most likely a member of the A-type subclass of W UMa binaries.     

The mass,orbit, and tidal evolution of the Quaoar–Weywot system

Here we present new adaptive optics observations of the Quaoar–Weywot system. With these new observations we determine an improved system orbit. Due to a 0.39 day alias that exists in available observations, four possible orbital solutions are available with periods of ～11.6, ～12.0, ～12.4, and ～12.8 days. From the possible orbital solutions, system masses of 1.3–1.5 ± 0.1 × 10²¹ kg are found. These observations provide an updated density for Quaoar of 2.7–5.0 g cm^?3. In all cases, Weywot’s orbit is eccentric, with possible values ～0.13–0.16. We present a reanalysis of the tidal orbital evolution of the Quaoar–Weywot system. We have found that Weywot has probably evolved to a state of synchronous rotation, and has likely preserved its initial inclination over the age of the Solar System. We find that for plausible values of the effective tidal dissipation factor tides produce a very slow evolution of Weywot’s eccentricity and semi-major axis. Accordingly, it appears that Weywot’s eccentricity likely did not tidally evolve to its current value from an initially circular orbit. Rather, it seems that some other mechanism has raised its eccentricity post-formation, or Weywot formed with a non-negligible eccentricity.     

The outcome of the protoplanetary disk of very massive stars

We suggest that planets, brown dwarfs, and even low mass stars can be formed by fragmentation of protoplanetary disks around very massive stars (M ? 100 M_⊙). We discuss how fragmentation conditions make the formation of very massive planetary systems around very massive stars favorable. Such planetary systems are likely to be composed of brown dwarfs and low mass stars of ～0.1–0.3 M_⊙, at orbital separations of ～ few × 100–10⁴ AU. In particular, scaling from solar-like stars suggests that hundreds of Mercury-like planets might orbit very massive stars at ～10³ AU where conditions might favor liquid water. Such fragmentation objects can be excellent targets for the James Webb Space Telescope and other large telescopes working in the IR bands. We predict that deep observations of very massive stars would reveal these fragmentation objects, orbiting in the same orbital plane in cases where there are more than one object.     

Stardust-NExT,Deep Impact,and the accelerating spin of 9P/Tempel 1

The evolution of the spin rate of Comet 9P/Tempel 1 through two perihelion passages (in 2000 and 2005) is determined from 1922 Earth-based observations taken over a period of 13 year as part of a World-Wide observing campaign and from 2888 observations taken over a period of 50 days from the Deep Impact spacecraft. We determine the following sidereal spin rates (periods): 209.023 ± 0.025°/dy (41.335 ± 0.005 h) prior to the 2000 perihelion passage, 210.448 ± 0.016°/dy (41.055 ± 0.003 h) for the interval between the 2000 and 2005 perihelion passages, 211.856 ± 0.030°/dy (40.783 ± 0.006 h) from Deep Impact photometry just prior to the 2005 perihelion passage, and 211.625 ± 0.012°/dy (40.827 ± 0.002 h) in the interval 2006–2010 following the 2005 perihelion passage. The period decreased by 16.8 ± 0.3 min during the 2000 passage and by 13.7 ± 0.2 min during the 2005 passage suggesting a secular decrease in the net torque. The change in spin rate is asymmetric with respect to perihelion with the maximum net torque being applied on approach to perihelion. The Deep Impact data alone show that the spin rate was increasing at a rate of 0.024 ± 0.003°/dy/dy at JD2453530.60510 (i.e., 25.134 dy before impact), which provides independent confirmation of the change seen in the Earth-based observations.The rotational phase of the nucleus at times before and after each perihelion and at the Deep Impact encounter is estimated based on the Thomas et al. (Thomas et al. [2007]. Icarus 187, 4–15) pole and longitude system. The possibility of a 180° error in the rotational phase is assessed and found to be significant. Analytical and physical modeling of the behavior of the spin rate through of each perihelion is presented and used as a basis to predict the rotational state of the nucleus at the time of the nominal (i.e., prior to February 2010) Stardust-NExT encounter on 2011 February 14 at 20:42.We find that a net torque in the range of 0.3–2.5 × 10⁷ kg m² s^?2 acts on the nucleus during perihelion passage. The spin rate initially slows down on approach to perihelion and then passes through a minimum. It then accelerates rapidly as it passes through perihelion eventually reaching a maximum post-perihelion. It then decreases to a stable value as the nucleus moves away from the Sun. We find that the pole direction is unlikely to precess by more than ～1° per perihelion passage. The trend of the period with time and the fact that the modeled peak torque occurs before perihelion are in agreement with published accounts of trends in water production rate and suggests that widespread H₂O out-gassing from the surface is largely responsible for the observed spin-up.     

Chronology and sources of lunar impact bombardment

The Moon has suffered intense impact bombardment ending at 3.9 Gyr ago, and this bombardment probably affected all of the inner Solar System. Basin magnetization signatures and lunar crater size-distributions indicate that the last episode of bombardment at about 3.85 Gyr ago was less extensive than previously thought. We explore the contribution of the primordial Mars-crosser population to early lunar bombardment. We find that Mars-crosser population initially decays with a 80-Myr half-life, with the long tail of survivors clustering on temporarily non-Mars-crossing orbits between 1.8 and 2 AU. These survivors decay with half-life of about 600 Myr and are progenitors of the extant Hungaria asteroid group in the same region. We estimate the primordial Mars-crosser population contained about 0.01–0.02 Earth masses. Such initial population is consistent with no lunar basins forming after 3.8 Gya and the amount of mass in the Hungaria group. As they survive longer and in greater numbers than other primordial populations, Mars-crossers are the best candidate for forming the majority of lunar craters and basins, including most of the Nectarian system. However, this remnant population cannot produce Imbrium and Orientale basins, which formed too late and are too large to be part of a smooth bombardment. We propose that the Imbrian basins and craters formed in a discrete event, consistent with the basin magnetization signatures and crater size-distributions. This late “impactor shower” would be triggered by a collisional disruption of a Vesta-sized body from this primordial Mars-crossing population (Wetherill, G.W. [1975]. Proc. Lunar Sci. Conf. 6, 1539–1561) that was still comparable to the present-day asteroid belt a 3.9 Gya. This tidal disruption lead to a short-lived spike in bombardment by non-chondritic impactors with a non-asteroidal size–frequency distribution, in agreement with available evidence. This body (“Wetherill’s object”) also uniquely matches the constraints for the parent body of mesosiderite meteorites. We propose that the present-day sources of mesosiderites are multi-km-sized asteroids residing in the Hungaria group, that have been implanted there soon after the original disruption of their parent 3.9 Gyr ago.     

Periodicity of flare index revisited using the Hilbert–Huang transform method

Using the Hilbert–Huang Transform (HHT) method, we investigate the periodicity in the monthly mean flare indices from 1966 to 2007 (corresponding to almost four complete solar cycles), calculated by T. Atac and A. Ozguc. The results show as following. (1) The periods of 9.37 ± 2.53, 11.8 ± 0.172 and 23.6 ± 0.316 years are found to be statistically significant in the flare index. The most eminent period is 9.37 ± 2.53 years. (2) Other periods of 0.237 ± 0.196 years (86.6 ± 71.6 days), 0.525 ± 0.0508 years (191 ± 18.5 days), 1.05 ± 0.478 years (383 ± 174 days) and 2.37 ± 0.395 years are below the 99% confidence level line, suggesting they are due to stochastic random noise.     

A CCD photometric study of the late type contact binary EK Comae Berenices

We present CCD photometric observations of the W UMa type contact binary EK Comae Berenices using the 2 m telescope of IUCAA Girawali Observatory, India. The star was classified as a W UMa type binary of subtype-W by Samec et al. (1996). The new V band photometric observations of the star reveal that shape of the light curve has changed significantly from the one observed by Samec et al. (1996). A detailed analysis of the light curve obtained from the high-precision CCD photometric observations of the star indicates that EK Comae Berenices is not a W-type but an A-type totally eclipsing W UMa contact binary. The photometric mass ratio is determined to be 0.349 ± 0.005. A temperature difference of ΔT = 141 ± 10 K between the components and an orbital inclination of i[°] = 89.800 ± 0.075 were obtained for the binary system. Absolute values of masses, radii and luminosities are estimated by means of the standard mass-luminosity relation for zero age main-sequence stars. The star shows O’Connell effect, asymmetries in the light curve shape around the primary and secondary maximum. The observed O’Connell effect is explained by the presence of a hot spot on the primary component.     

Efficiency of a wide-area survey in achieving short- and long-term warning for small impactors

We consider a network of telescopes capable of scanning all the observable sky each night and targeting Near-Earth objects (NEOs) in the size range of the Tunguska-like asteroids, from 160 m down to 10 m. We measure the performance of this telescope network in terms of the time needed to discover at least 50% of the impactors in the considered population with a warning time large enough to undertake proper mitigation actions. The warning times are described by a trimodal distribution and the telescope network has a 50% probability of discovering an impactor of the Tunguska class with at least one week of advance already in the first 10 yr of operations of the survey. These results suggest that the studied survey would be a significant addition to the current NEO discovery efforts.     

The detached eclipsing binary TX Her revisited

This paper presents new CCD Bessell BVRI light curves and photometric analysis of the Algol-type binary star TX Her. The CCD observations were carried out at Çanakkale Onsekiz Mart University Observatory in 2010. New BVRI light curves from this study and radial velocity curves from Popper (1970) were solved simultaneously using modern light and radial velocity curves synthesis methods. The general results show that TX Her is a well-detached eclipsing binary, however, both component stars fill at least half of their Roche lobes. A significant third light contribution to the total light of the system could not be determined. Using O–C residuals formed by the updated minima times, an orbital period study of the system was performed. It was confirmed that the tilted sinusoidal O–C variation corresponds to an apparent period variation caused by the light travel time effect due to an unseen third body. The following absolute parameters of the components were derived: M₁ = 1.62 ± 0.04 M_⊙, M₂ = 1.45 ± 0.03 M_⊙, R₁ = 1.69 ± 0.03 R_⊙, R₂ = 1.43 ± 0.03 R_⊙, L₁ = 8.21 ± 0.90 L_⊙ and L₂ = 3.64 ± 0.60 L_⊙. The distance to TX Her was calculated as 155 ± 10 pc, taking into account interstellar extinction. The position of the components of TX Her in the HR diagram are also discussed. The components are young stars with an age of ～500 Myr.     

A local interstellar spectrum for galactic electrons

A heliopause spectrum at 122 AU from the Sun is presented for galactic electrons over an energy range from 1 MeV to 50 GeV that can be considered the lowest possible local interstellar spectrum (LIS). The focus of this work is on the spectral shape of the LIS below ∼1.0 GeV. The study is done by using a comprehensive numerical model for solar modulation in comparison with Voyager 1 observations at ∼112 AU from the Sun and PAMELA data at Earth. Below ∼1.0 GeV, this LIS exhibits a power law with E^{−(1.55 ± 0.05)}, where E is the kinetic energy of these electrons. However, reproducing the PAMELA electron spectrum averaged for 2009, requires a LIS with a different power law of the form E^{−(3.15 ± 0.05)} above ∼5 GeV. Combining the two power laws with a smooth transition from low to high energies yields a LIS over the full energy range that is relevant and applicable to the modulation of cosmic ray electrons in the heliosphere. The break occurs between ∼800 MeV and ∼2 GeV as a characteristic feature of this LIS. The power-law form below ∼1 GeV produces a challenge to the origin of these low energy galactic electrons. On the other hand, the results of this study can be used as a gauge for astrophysical modeling of the local interstellar spectrum for electrons.     

Rotation period of Venus estimated from Venus Express VIRTIS images and Magellan altimetry

The 1.02 μm wavelength thermal emission of the nightside of Venus is strongly anti-correlated to the elevation of the surface. The VIRTIS instrument on Venus Express has mapped this emission and therefore gives evidence for the orientation of Venus between 2006 and 2008. The Magellan mission provided a global altimetry data set recorded between 1990 and 1992. Comparison of these two data sets reveals a deviation in longitude indicating that the rotation of the planet is not fully described by the orientation model recommended by the IAU. This deviation is sufficiently large to affect estimates of surface emissivity from infrared imaging. A revised period of rotation of Venus of 243.023 ± 0.002 d aligns the two data sets. This period of rotation agrees with pre-Magellan estimates but is significantly different from the commonly accepted value of 243.0185 ± 0.0001 d estimated from Magellan radar images. It is possible that this discrepancy stems from a length of day variation with the value of 243.023 ± 0.002 d representing the average of the rotation period over 16 years.     

A photometric study of a weak-contact binary: V873 Per

We present a photometric study of a weak-contact binary V873 Per. New observations in BVR filter bands showed asymmetric light curves to be a negative type of the O’Connell effect, which can be described by magnetic activity of a cool spot on the more massive component. Our photometric solutions showed that V873 Per is a W-type with a mass ratio of q = 2.504(±0.0029), confirming the results of Samec et al. (2009). The derived contact degree was found to be f = 18.10%(±1.36%). Moreover, our analysis found the cyclic variation with the period of about 4 yr that could be due to existence of the third companion in the system or the mechanism of magnetic activity cycle in the binary. While available data indicated that the long-term orbital period tends to be stable rather than decreasing.