Gamma–hadron separation methods for the VERITAS array of four imaging atmospheric Cherenkov telescopes

Ground-based arrays of imaging atmospheric Cherenkov telescopes have emerged as the most sensitive γ-ray detectors in the energy range of about 100 GeV and above. The strengths of these arrays are a very large effective collection area on the order of 10⁵ m², combined with excellent single photon angular and energy resolutions. The sensitivity of such detectors is limited by statistical fluctuations in the number of Cosmic-ray initiated air showers that resemble γ-ray air showers in many ways. In this paper, we study the performance of simple event reconstruction methods when applied to simulated data of the Very Energetic Radiation Imaging Telescope Array System (VERITAS) experiment. We review methods for reconstructing the arrival direction and the energy of the primary photons, and examine means to improve on their performance. For a software threshold energy of 300 GeV (100 GeV), the methods achieve point source angular and energy resolutions of σ_63% = 0.1° (0.2°) and σ_68% = 15% (22%), respectively. The main emphasis of the paper is the discussion of γ–hadron separation methods for the VERITAS experiment. We find that the information from several methods can be combined based on a likelihood ratio approach and the resulting algorithm achieves a γ–hadron suppression with a quality factor that is substantially higher than that achieved with the standard methods used so far.     

A photometric study of the detached binary WY Hydrae with twin components

We present a multicolor photometry for the eclipsing binary WY Hydrae, observed on four nights of 2008 December. From our new observations and Carr’s data, the photometric solutions were deduced by using the updated W–D program. The results show that WY Hya is a detached binary with a mass ratio of q = 0.970(±0.005).By analyzing the O–C curve, it is found that there exists either a continuous period increase or a cyclic variation. From Eq. (2), the orbital period of WY Hya secularly increases at a rate of dP/dt = +3.56(±0.37) × 10^?7 days/yr, which may be interpreted by some mass transfer for the near-contact configuration or tidal dissipation. From Eq. (3), the period and semi-amplitude of the periodic oscillation are P₃ = 95.4(±4.2) yr and A = 0^d.0087(±0^d.0003), respectively. This may be likely attributed by light-time effect via the presence of the assumed third body. Assumed in the coplanar orbit with the binary, the mass of the third body should be M₃ = 0.18 M_⊙. If the unseen additional companion exists, it will extract angular momentum from the binary system. Finally, WY Hya with high fill-out factors (i.e., f_1,2 > 80%), may evolve into a semi-detached configuration.     

Observation of DCl and upper limit to NH3 on Venus

To search for DCl in the Venus atmosphere, a spectrum near the D³⁵Cl (1–0) R4 line at 2141.54 cm^?1 was observed using the CSHELL spectrograph at NASA IRTF. Least square fitting to the spectrum by a synthetic spectrum results in a DCl mixing ratio of 17.8 ± 6.8 ppb. Comparing to the HCl abundance of 400 ± 30 ppb (Krasnopolsky [2010a] Icarus, 208, 314–322), the DCl/HCl ratio is equal to 280 ± 110 times the terrestrial D/H = 1.56 × 10^?4. This ratio is similar to that of HDO/H₂O = 240 ± 25 times the terrestrial HDO/H₂O from the VEX/SOIR occultations at 70–110 km. Photochemistry in the Venus mesosphere converts H from HCl to that in H₂O with a rate of 1.9 × 10⁹ cm^?2 s^?1 (Krasnopolsky [2012] Icarus, 218, 230–246). The conversion involves photolysis of HCl; therefore, the photochemistry tends to enrich D/H in HCl and deplete in H₂O. Formation of the sulfuric acid clouds may affect HDO/H₂O as well. The enriched HCl moves down by mixing to the lower atmosphere where thermodynamic equilibriums for H₂ and HCl near the surface correspond to D/H = 0.71 and 0.74 times that in H₂O, respectively. Time to establish these equilibriums is estimated at ～3 years and comparable to the mixing time in the lower atmosphere. Therefore, the enriched HCl from the mesosphere gives D back to H₂O near the surface. Comparison of chemical and mixing times favors a constant HDO/H₂O up to ～100 km and DCl/HCl equal to D/H in H₂O times 0.74.Ammonia is an abundant form of nitrogen in the reducing environments. Thermodynamic equilibriums with N₂ and NO near the surface of Venus give its mixing ratio of 10^?14 and 6 × 10^?7, respectively. A spectrum of Venus near the NH₃ line at 4481.11 cm^?1 was observed at NASA IRTF and resulted in a two-sigma upper limit of 6 ppb for NH₃ above the Venus clouds. This is an improvement of the previous upper limit by a factor of 5. If ammonia exists at the ppb level or less in the lower atmosphere, it quickly dissociates in the mesosphere and weakly affects its photochemistry.     

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.     

Outer Main Belt asteroids: Identification and distribution of four 3-μm spectral groups

This paper examines the distribution and the abundance of hydrated minerals (any mineral that contains H₂O or OH) on outer Main Belt asteroids spanning the 2.5 < a < 4.0 AU region. The hypothesis we are testing is whether planetesimals that accreted closer to the Sun experienced a higher degree of aqueous alteration. We would expect then to see a gradual decline of the abundance of hydrated minerals among the outer Main Belt asteroids with increasing heliocentric distance (2.5 < a < 4.0 AU). We measured spectra (0.8–2.5 μm and 1.9–4.1 μm) of 28 outer Main Belt asteroids using the SpeX spectrograph/imager at the NASA Infrared Telescope Facility (IRTF). We identified four groups on the basis of the shape and the band center of the 3-μm feature. The first group, which we call “sharp”, exhibits a sharp 3-μm feature, attributed to hydrated minerals (phyllosilicates). Most asteroids in this group are located in the 2.5 < a < 3.3 AU region. The second group, which we call “Ceres-like”, consists of 10 Hygiea and 324 Bamberga. Like Asteroid Ceres, these asteroids exhibit a 3-μm feature with a band center of 3.05 ± 0.01 μm that is superimposed on a broader absorption feature from ～2.8 to 3.7 μm. The third group, which we call “Europa-like”, includes 52 Europa, 31 Euphrosyne, and 451 Patientia. Objects in this group exhibit a 3-μm feature with a band center of 3.15 ± 0.01 μm. Both the Ceres-like and Europa-like groups are concentrated in the 2.5 < a < 3.3 AU region. The fourth group, which we call “rounded”, is concentrated in the 3.4 < a < 4.0 AU region. Asteroids in this group are characterized by a rounded 3-μm feature, attributed to H₂O ice. A similar rounded 3-μm feature was also identified in 24 Themis and 65 Cybele. Unlike the sharp group, the rounded group did not experience aqueous alteration. Of the asteroids observed in this study, 140 Siwa, a P-type, is the only one that does not exhibit a 3-μm feature. These results are important to constrain the nature and the degree of aqueous alteration in outer Main Belt asteroids.     

The 2009–2010 fade of Jupiter’s South Equatorial Belt: Vertical cloud structure models and zonal winds from visible imaging

The South Equatorial Belt (SEB) of Jupiter is known to alternate its appearance at visible wavelengths from a classical belt-like band most of the time to a short-lived zone-like aspect which is called a “fade” of the belt, hereafter SEBF. The albedo change of the SEB is due to a change in the structure and properties of the clouds and upper hazes. Recent works based on infrared observations of the last SEBF have shown that the aerosol density below 1 bar increased in parallel with the reflectivity change. However, the nature of the change in the upper clouds and hazes that produces the visible reflectivity change and whether or not this reflectivity change is accompanied by a change in the winds at the upper cloud level remained unknown. In this paper we focus in the near ultraviolet to near infrared reflected sunlight (255–953 nm) to address these two issues. We characterize the vertical cloud structure above the ammonia condensation level from Hubble Space Telescope images, and the zonal wind velocities from long-term high-quality images coming from the International Outer Planet Watch database, both during the SEB and SEBF phases. We show that reflectivity changes do not happen simultaneously in this wavelength range, but they start earlier in the most deep-sensing filters and end in 2010 with just minor changes in those sensing the highest particle layers. Our models require a substantial increase of the optical thickness of the cloud deck at 1.0 ± 0.4 bar from τ_cloud = 6 ± 2 in July 2009 (SEB phase) to semiinfinite at visual wavelengths in 2010 (SEBF). Upper tropospheric particles (～240–1400 mbar) are also required to become substantially reflectant and their single scattering albedo in the blue changes from ?₀ = 0.905 ± 0.005 in November 2009 to ?₀ = 0.95 ± 0.01 in June 2010. No significant changes were found in the cloud top heights or in the particle density of the tropospheric haze. The disturbance travels from the levels below ～3 bar to a level about 400 ± 100 mbar. We derive an upward velocity of 0.15 ± 0.05 cm/s, in agreement with a diffusive process in Jupiter’s upper troposphere requiring a mean eddy coefficient K ～ 8 × 10⁵ cm² s^?1. On the other hand, cloud tracking on the IOPW imaging showed no significant changes in the zonal wind profile between the SEB and SEBF stages. As in other visually huge changes in Jupiter’s cloud morphology and structure, the wind profile remains robust, possibly indicating a deeply rooted dynamical regime.     

CCD photometric analysis of the W UMa-type binary V376 Andromeda

This study presents the absolute parameters of the contact binary system V376 And. CCD photometric observations were made at the Çanakkale Onsekiz Mart University Observatory in 2004. The instrumental magnitudes of all observed stars were converted into standard magnitudes. New BV light curves of the system were analysed using the Wilson–Devinney method supplemented with a Monte Carlo type algorithm. Since there are large asymmetries between maxima (i.e., O’Connell effect) in these light curves, two different models (one with a cool spot and one with a hot spot) were applied to the photometric data. The best fit, which was obtained with a large hot spot on the secondary component, gives V376 And as an A sub-type contact binary in poor thermal contact and a small value of the filling factor (f ≈ 0.07). Combining the solutions of our light curves and Rucinski et al. (2001)’s radial velocity curves, the following absolute parameters of the components were determined: M₁ = 2.44 ± 0.04 M_⊙, M₂ = 0.74 ± 0.03 M_⊙, R₁ = 2.60 ± 0.03 R_⊙, R₂ = 1.51 ± 0.02 R_⊙, L₁ = 40 ± 4 L_⊙ and L₂ = 5 ± 1 L_⊙. We also discuss the evolution of the system, which appears to have an age of 1.6 Gyr. The distance to V376 And was calculated as 230 ± 20 pc from this analysis, taking into account interstellar extinction.     

Optical constants from 370 nm to 900 nm of Titan tholins produced in a low pressure RF plasma discharge

Determining the optical constants of Titan aerosol analogues, or tholins, has been a major concern for the last three decades because they are essential to constrain the numerical models used to analyze Titan’s observational data (albedo, radiative transfer, haze vertical profile, surface contribution, etc.). Here we present the optical constant characterization of tholins produced with an RF plasma discharge in a (95%N₂–5%CH₄) gas mixture simulating Titan’s main atmospheric composition, and deposited as a thin film on an Al–SiO₂ substrate. The real and imaginary parts, n and k, of the tholin complex refractive index have been determined from 370 nm to 900 nm wavelength using spectroscopic ellipsometry. The values of n decrease from n = 1.64 (at 370 nm) to n = 1.57 (at 900 nm) as well as the values of k which feature two behaviors: an exponential decay from 370 nm to 500 nm, with k = 12.4 × e^?0.018λ (where λ is expressed in nm), followed by a plateau, with k = (1.8 ± 0.2) × 10^?3. The trends observed for the PAMPRE tholins optical constants are compared to those determined for other Titan tholins, as well as to the optical constants of Titan’s aerosols retrieved from observational data.     

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.     

Chandra X-ray observations of two unusual BAL quasars

We report sensitive Chandra X-ray non-detections of two unusual, luminous Iron Low-Ionization Broad Absorption Line Quasars (FeLoBALs). The observations do detect a non-BAL, wide-binary companion quasar to one of the FeLoBAL quasars. We combine X-ray-derived column density lower limits (assuming solar metallicity) with column densities measured from ultraviolet spectra and CLOUDY photoionization simulations to explore whether constant-density slabs at broad-line region densities can match the physical parameters of these two BAL outflows, and find that they cannot. In the “overlapping-trough” object SDSS J0300+0048, we measure the column density of the X-ray absorbing gas to be N_H ? 1.8 × 10²⁴ cm^?2. From the presence of Fe ii UV78 absorption but lack of Fe ii UV195/UV196 absorption, we infer the density in that part of the absorbing region to be n_e ? 10⁶ cm^?3. We do find that a slab of gas at that density might be able to explain this object’s absorption. In the Fe iii-dominant object SDSS J2215–0045, the X-ray absorbing column density of N_H ? 3.4 × 10²⁴ cm^?2 is consistent with the Fe iii-derived N_H ? 2 × 10²² cm^?2 provided the ionization parameter is log U > 1.0 for both the n_e = 10¹¹ cm^?3 and n_e = 10¹² cm^?3 scenarios considered (such densities are required to produce Fe iii absorption without Fe iiabsorption). However, the velocity width of the absorption rules out its being concentrated in a single slab at these densities. Instead, this object’s spectrum can be explained by a low density, high ionization and high temperature disk wind that encounters and ablates higher density, lower ionization Fe iii-emitting clumps.     

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.     

Ensemble properties of comets in the Sloan Digital Sky Survey

We present the ensemble properties of 31 comets (27 resolved and 4 unresolved) observed by the Sloan Digital Sky Survey (SDSS). This sample of comets represents about 1 comet per 10 million SDSS photometric objects. Five-band (u, g, r, i, z) photometry is used to determine the comets’ colors, sizes, surface brightness profiles, and rates of dust production in terms of the Afρ formalism. We find that the cumulative luminosity function for the Jupiter Family Comets in our sample is well fit by a power law of the form N(<H) ∝ 10^{(0.49±0.05)H} for H < 18, with evidence of a much shallower fit N(<H) ∝ 10^{(0.19±0.03)H} for the faint (14.5 < H < 18) comets. The resolved comets show an extremely narrow distribution of colors (0.57 ± 0.05 in g ? r for example), which are statistically indistinguishable from that of the Jupiter Trojans. Further, there is no evidence of correlation between color and physical, dynamical, or observational parameters for the observed comets.     

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.     

A photometric study of the W UMa-type binary DF Canum Venaticorum

We present new photometric observations for the eclipsing binary DF CVn, and determined five light minimum times. By using the Wilson–Devinney code, two sets of photometric solutions were deduced from our observations in 2009. The asymmetric light curves obtained on 2009 March 5 were modeled by a dark spot on the more massive component. The results indicate that DF CVn is a W-type weak-contact binary, with a mass ratio of q ～ 0.28 and an overcontact degree of f ～ 20%. From the O ? C curve of minimum times, it is found that there exists a cyclic variation, whose period and semi-amplitude are P₃ = 17.2(±0.9) year and A = 0.^d0070(±0.^d0008), respectively. This kind of cyclic oscillation may possibly result from the light-time effect due to the presence of an unseen third body. This kind of additional body may extract angular momentum from the central system. The low-amplitude changes of the light curves on a short-time scale (e.g., half a month) may be attributed to the dark spot activity, which may result in angular momentum loss via magnetic breaking. With angular momentum loss, the weak-contact binary DF CVn will evolve into a deep-contact configuration.     

BV photometry of quadrupole system ET Bootis

In this study, we present the first Johnson BV photometry of the eclipsing binary star ET Bootis, which is member of a physically connected visual pair. Analysis of times of light minima enables us to calculate accurate ephemeris of the system via O–C analysis and observed an increase in period which we believe is a result of the light-time effect in the outer visual orbit. Secondly, we determined the total brightness and color of the system in light maxima and minima. Photometric solution of the system indicates that the contribution of the visual pair to the total light is about 40% in Johnson V band. Furthermore, photometric analysis shows that the primary star in the eclipsing binary has F8 spectral type while it confirms the G5 spectral type for the visual pair. Masses of the components in eclipsing binary are M₁ = 1.109 ± 0.014 M_⊙ and M₂ = 1.153 ± 0.011 M_⊙. Absolute radii of the components are R₁ = 1.444 ± 0.007 R_⊙ and R₂ = 1.153 ± 0.007 R_⊙. Physical properties of the components leads 176 ± 7 pc distance for the system and suggests an age of 6.5 billion years.     

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.     

A new absolute magnitude calibration with 2MASS for cataclysmic variables

Using reliable trigonometric measurements, we find that the absolute magnitude of cataclysmic variables depends on the orbital period and de-reddened (J − H)₀ and (H − K _s)₀ colours of 2MASS (Two Micron All Sky Survey) photometric system. The calibration equation covers the ranges 0.032^d < P_orb ⩽ 0.454^d, −0.08 < (J − H)₀ ⩽ 1.54, −0.03 < (H − K_s)₀ ⩽ 0.56 and 2.0 < M_J < 11.7; It is based on trigonometric parallaxes with relative errors of (σ_π/π) ⩽ 0.4. By using the period-luminosity-colours (PLCs) relation, we estimated the distances of cataclysmic variables with orbital periods and 2MASS observations and compared them with distances found from other methods. We suggest that the PLCs relation can be a useful statistical tool to estimate the distances of cataclysmic variables.     

Revealing the parameters of the open cluster Ruprecht 58

We present results of a study that combines UBVI photometry, MK spectral classification and proper motions in the area of the, up to now unknown, open cluster Ruprecht 58 at the Puppis region. Star counts from the 2MASS data catalog together with the analysis of CCD UBVI photometry demonstrate that it is a real open cluster with 9′ size approximately. The cluster is placed at a distance of 3.9 kpc and is about 250 Myr old with mean reddening E_(B−V) = 0.33 mag. Proper motions confirm Ruprecht 58 is a real cluster with mean absolute proper motions μ_αcosδ = −2.77 ± 0.45 mas/yr and μ_δ = 4.54 ± 0.45 mas/yr in the magnitude range 13.5 < V < 14.5 and μ_αcosδ = −2.70 ± 0.32 mas/yr and μ_δ = 3.19 ± 0.32 mas/yr in the range 14.5 < V < 16.0. The computation of the cluster mass spectrum slope yielded x = 1.8 in the mass range from ≈1.4 to ≈4m_⊙.     

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.