A Be star with a low nitrogen abundance in the SMC cluster NGC330

High- resolution UVES/ VLT spectra of B 12, an extreme pole- on Be star in the SMC cluster NGC 330, have been analysed using non-LTE model atmospheres to obtain its chemical composition relative to the SMC standard star AV304. We find a general underabundance of metals which can be understood in terms of an extra contribution to the stellar continuum due to emission from a disk which we estimate to be at the similar to 25% level. When this is corrected for, the nitrogen abundance for B12 shows no evidence of enhancement by rotational mixing as has been found in other non-Be B-type stars in NGC 330, and is inconsistent with evolutionary models which include the effects of rotational mixing. A second Be star, NGC330-B 17, is also shown to have no detectable nitrogen lines. Possible explanations for the lack of rotational mixing in these rapidly rotating stars are discussed, one promising solution being the possibility that magnetic fields might inhibit rotational mixing.

Photometry of cometary nuclei: Rotation rates, colours and a comparison with Kuiper Belt Objects

We present time-series data on Jupiter Family Comets (JFCs) 17P/Holmes, 47P/Ashbrook-Jackson and 137P/Shoemaker-Levy 2. In addition we also present results from `snap-shot' observations of comets 43P/Wolf-Harrington, 44P/Reinmuth 2, 103P/Hartley 2 and 104P/Kowal 2 taken during the same run. The comets were at heliocentric distances of between 3 and 7 au at this time. We present measurements of size and activity levels for the snap-shot targets. The time-series data allow us to constrain rotation periods and shapes, and thus bulk densities. We also measure colour indices (V - R) and (R - I) and reliable radii for these comets. We compare all of our findings to date with similar results for other comets and Kuiper Belt Objects (KBOs). We find that the rotational properties of nuclei and KBOs are very similar, that there is evidence for a cut-off in bulk densities at ~0.6 g cm-3 in both populations, and the colours of the two populations show similar correlations. For JFCs, there is no observational evidence for the optical colours being dependent on either position in the orbit or orbital parameters.

Early-type stars in the Galactic halo from the Palomar-Green survey - III. Completion of a magnitude range limited sample

High-resolution (R approximate to 40 000) echelle spectroscopic observations of 13 high-latitude early-type stars are presented. These stars comprise the final part of a complete magnitude range limited sample based on low-resolution spectroscopy of targets drawn from the Palomar-Green survey. The magnitude range under consideration is 13 less than or equal to B-PG less than or equal to 14.6, corresponding to an approximate distance limit for main-sequence B-type objects of 5 less than or equal to d less than or equal to 40 kpc. Three stars are found to be apparently normal, young stars, based on their positions on the (T-eff, log g) diagram, normal abundance patterns and relatively large projected rotational velocities. A further star, PG 1209+263, was found to belong to the chemically peculiar (CP) silicon star class of objects. The remainder are evolved subluminous stars lying on post- horizontal branch (post-HB) tracks, with the exception of PG 2120+062, which appears to be in a post-asymptotic giant branch evolutionary stage. For the young stars in the sample, we have derived distance and age estimates through comparison of the atmospheric parameters with recent theoretical evolutionary models. We discuss formation scenarios by comparing times-of- flight and evolutionary time-scales. It is found that all stars could have formed in the Galactic disc and been ejected from there soon after their birth, with the exception of PG 1209+263. The adopted proper motion is found to be a crucial factor in the kinematical analysis. We also present some number densities for young B-type halo stars, which indicate that they are extremely scarce objects.

Computational challenges in atomic, molecular and optical physics

Six challenges are discussed. These are the laser-driven helium atom; the laser-driven hydrogen molecule and hydrogen molecular ion: electron scattering (with ionization) from one-electron atoms; the vibrational and rotational structure of molecules such as H-3(+) and water at their dissociation limits; laser- heated clusters; and quantum degeneracy and Bose-Einstein condensation. The first four concern fundamental few-body systems where use of high-performance computing (HPC) is currently making possible accurate modelling from first principles. This leads to reliable predictions and support for laboratory experiment as well as true understanding of the dynamics. Important aspects of these challenges addressable only via a terascale facility are set out. Such a facility makes the last two challenges in the above list meaningfully accessible for the first time, and the scientific interest together with the prospective role for HPC in these is emphasized.

Escape times for rigid Brownian rotators in a bistable potential from the time evolution of the Green function and the characteristic time of the probability evolution

The greatest relaxation time for an assembly of three- dimensional rigid rotators in an axially symmetric bistable potential is obtained exactly in terms of continued fractions as a sum of the zero frequency decay functions (averages of the Legendre polynomials) of the system. This is accomplished by studying the entire time evolution of the Green function (transition probability) by expanding the time dependent distribution as a Fourier series and proceeding to the zero frequency limit of the Laplace transform of that distribution. The procedure is entirely analogous to the calculation of the characteristic time of the probability evolution (the integral of the configuration space probability density function with respect to the position co-ordinate) for a particle undergoing translational diffusion in a potential; a concept originally used by Malakhov and Pankratov (Physica A 229 (1996) 109). This procedure allowed them to obtain exact solutions of the Kramers one-dimensional translational escape rate problem for piecewise parabolic potentials. The solution was accomplished by posing the problem in terms of the appropriate Sturm-Liouville equation which could be solved in terms of the parabolic cylinder functions. The method (as applied to rotational problems and posed in terms of recurrence relations for the decay functions, i.e., the Brinkman approach c.f. Blomberg, Physica A 86 (1977) 49, as opposed to the Sturm-Liouville one) demonstrates clearly that the greatest relaxation time unlike the integral relaxation time which is governed by a single decay function (albeit coupled to all the others in non-linear fashion via the underlying recurrence relation) is governed by a sum of decay functions. The method is easily generalized to multidimensional state spaces by matrix continued fraction methods allowing one to treat non-axially symmetric potentials, where the distribution function is governed by two state variables. (C) 2001 Elsevier Science B.V. All rights reserved.

The role of Mittag-Leffler functions in anomalous relaxation

The inertia-corrected Debye model of rotational Brownian motion of polar molecules was generalized by Coffey et al. [Phys. Rev. E, 65, 32 102 (2002)] to describe fractional dynamics and anomalous rotational diffusion. The linear- response theory of the normalized complex susceptibility was given in terms of a Laplace transform and as a function of frequency. The angular-velocity correlation function was parametrized via fractal Mittag-Leffler functions. Here we apply the latter method and complex-contour integral- representation methods to determine the original time-dependent amplitude as an inverse Laplace transform using both analytical and numerical approaches, as appropriate. (C) 2004 Elsevier B.V. All rights reserved.

An exactly solvable four transition point model

Heavy particle collisions, in particular low-energy ion-atom collisions, are amenable to semiclassical JWKB phase integral analysis in the complex plane of the internuclear separation. Analytic continuation in this plane requires due attention to the Stokes phenomenon which parametrizes the physical mechanisms of curve crossing, non-crossing, the hybrid Nikitin model, rotational coupling and predissociation. Complex transition points represent adiabatic degeneracies. In the case of two or more such points, the Stokes constants may only be completely determined by resort to the so-called comparison- equation method involving, in particular, parabolic cylinder functions or Whittaker functions and their strong-coupling asymptotics. In particular, the Nikitin model is a two transition-point one-double-pole problem in each half-plane corresponding to either ingoing or outgoing waves. When the four transition points are closely clustered, new techniques are required to determine Stokes constants. However, such investigations remain incomplete, A model problem is therefore solved exactly for scattering along a one-dimensional z-axis. The energy eigenvalue is b(2)-a(2) and the potential comprises -z(2)/2 (parabolic) and -a(2) + b(2)/2z(2) (centrifugal/centripetal) components. The square of the wavenumber has in the complex z-plane, four zeros each a transition point at z = +/-a +/- ib and has a double pole at z = 0. In cases (a) and (b), a and b are real and unitarity obtains. In case (a) the reflection and transition coefficients are parametrized by exponentials when a(2) + b(2) > 1/2. In case (b) they are parametrized by trigonometrics when a(2) + b(2) <1/2 and total reflection is achievable. In case (c) a and b are complex and in general unitarity is not achieved due to loss of flux to a continuum (O'Rourke and Crothers, 1992 Proc. R. Sec. 438 1). Nevertheless, case (c) coefficients reduce to (a) or (b) under appropriate limiting conditions. Setting z = ht, with h a real constant, an attempt is made to model a two-state collision problem modelled by a pair of coupled first-order impact parameter equations and an appropriate (T) over tilde-tau relation, where (T) over tilde is the Stueckelberg variable and tau is the reduced or scaled time. The attempt fails because (T) over tilde is an odd function of tau, which is unphysical in a real collision problem. However, it is pointed out that by applying the Kummer exponential model to each half-plane (O'Rourke and Crothers 1994 J. Phys. B: At. Mel. Opt. Phys. 27 2497) the current model is in effect extended to a collision problem with four transition points and a double pole in each half-plane. Moreover, the attempt in itself is not a complete failure since it is shown that the result is a perfect diabatic inelastic collision for a traceless Hamiltonian matrix, or at least when both diagonal elements are odd and the off-diagonal elements equal and even.

Interpolation formula between very low and intermediate-to-high damping Kramers escape rates for single-domain ferromagnetic particles

It is shown that the Mel'nikov-Meshkov formalism for bridging the very low damping (VLD) and intermediate-to-high damping (IHD) Kramers escape rates as a function of the dissipation parameter for mechanical particles may be extended to the rotational Brownian motion of magnetic dipole moments of single-domain ferromagnetic particles in nonaxially symmetric potentials of the magnetocrystalline anisotropy so that both regimes of damping, occur. The procedure is illustrated by considering the particular nonaxially symmetric problem of superparamagnetic particles possessing uniaxial anisotropy subject to an external uniform field applied at an angle to the easy axis of magnetization. Here the Mel'nikov-Meshkov treatment is found to be in good agreement with an exact calculation of the smallest eigenvalue of Brown's Fokker-Planck equation, provided the external field is large enough to ensure significant departure from axial symmetry, so that the VLD and IHD formulas for escape rates of magnetic dipoles for nonaxially symmetric potentials are valid.

Microscopic models for dielectric relaxation in disordered systems

It is shown how the Debye rotational diffusion model of dielectric relaxation of polar molecules (which may be described in microscopic fashion as the diffusion limit of a discrete time random walk on the surface of the unit sphere) may be extended to yield the empirical Havriliak-Negami (HN) equation of anomalous dielectric relaxation from a microscopic model based on a kinetic equation just as in the Debye model. This kinetic equation is obtained by means of a generalization of the noninertial Fokker-Planck equation of conventional Brownian motion (generally known as the Smoluchowski equation) to fractional kinetics governed by the HN relaxation mechanism. For the simple case of noninteracting dipoles it may be solved by Fourier transform techniques to yield the Green function and the complex dielectric susceptibility corresponding to the HN anomalous relaxation mechanism.

More nitrogen rich B-type stars in the SMC cluster, NGC 330

High resolution spectra of seven early B-type giant/supergiant stars in the SMC cluster NGC330 are analysed to obtain their chemical compositions relative to SMC field and Galactic B-type stars. It is found that all seven stars are nitrogen rich with an abundance approximately 1.3 dex higher than an SMC main- sequence field B-type star, AV304. They also display evidence for deficiencies in carbon, but other metals have abundances typical of the SMC. Given the number of B-type stars with low projected rotational velocities in NGC330 (all our targets have v sin i <50 km s(-1)), we suggest that it is unlikely that the stars in our sample are seen almost pole-on, but rather that they are intrinsically slow rotators. Furthermore, none of our objects displays any evidence of significant Balmer emission excluding the possibility that these are Be stars observed pole-on. Comparing these results with the predictions of stellar evolution models including the effects of rotationally induced mixing, we conclude that while the abundance patterns may indeed be reproduced by these models, serious discrepancies exist. Most importantly, models including the effects of initially large rotational velocities do not reproduce the observed range of effective temperatures of our sample, nor the currently observed rotational velocities. Binary models may be able to produce stars in the observed temperature range but again may be incapable of producing suitable analogues with low rotational velocities. We also discuss the clear need for stellar evolution calculations employing the correct chemical mix of carbon, nitrogen and oxygen for the SMC.

A chemical analysis of five hot stars towards the Galactic centre

High resolution echelle spectroscopy is presented for thirteen stars lying in the direction of the Galactic centre which, on the basis of photographic photometry and low dispersion spectroscopy, have been classified as early-B-type. Eight of these stars have large rotational velocities which preclude a detailed analysis. The five stars with moderate to low projected rotational velocities have been analysed using model atmosphere techniques to determine atmospheric parameters and chemical compositions. Two of these stars appear to be evolved blue horizontal branch objects on the basis of their chemical compositions and small projected rotational velocity. The evolutionary status of a third is ambiguous but it is probably a post-asymptotic-giant branch star. The remaining two objects are probably young massive stars and show enhanced abundances of N, C, Mg and Si, consistent with their formation in the inner part of the Galactic disk. However their O abundances are normal, confirming results found previously for other early- type stars, which would imply a flat abundance gradient for this element in the inner region of our Galaxy.

Understanding B-type supergiants in the low metallicity environment of the SMC

Spectroscopic analyses of 7 SMC B-type supergiants and 1 giant have been undertaken using high resolution optical data obtained on the VLT with UVES. FASTWIND, a non-LTE, spherical, line-blanketed model atmosphere code was used to derive atmospheric and wind parameters of these stars as well as their absolute abundances. Mass-loss rates, derived from H-alpha profiles, are in poor agreement with metallicity dependent theoretical predictions. Indeed the wind-momenta of the SMC stars appear to be in good agreement with the wind-momentum luminosity relationship (WLR) of Galactic B-type stars, a puzzling result given that line-driven wind theory predicts a metallicity dependence. However the galactic stars were analysed using unblanketed model atmospheres which may mask any dependence on metallicity. A mean nitrogen enhancement of a factor of 14 is observed in the supergiants whilst only an enrichment of a factor of 4 is present in the giant, AV216. Similar excesses in nitrogen are observed in O-type dwarfs and supergiants in the same mass range, suggesting that the additional nitrogen is produced while the stars are still on the main-sequence. These nitrogen enrichments can be reproduced by current stellar evolution models, which include rotationally induced mixing, only if large initial rotational velocities of 300 kin s(-1) are invoked. Such large rotational velocities appear to be inconsistent with observed v sin i distributions for O-type stars and B-type supergiants. Hence it is suggested that the currently available stellar evolution models require more efficient mixing for lower rotational velocities.

The VLT-FLAMES survey of massive stars: mass loss and rotation of early-type stars in the SMC

We have studied the optical spectra of a sample of 31 O- and early B-type stars in the Small Magellanic Cloud, 21 of which are associated with the young massive cluster NGC 346. Stellar parameters are determined using an automated fitting method (Mokiem et al. 2005, A&A, 441, 711), which combines the stellar atmosphere code FASTWIND (Puls et al. 2005, A&A, 435, 669) with the genetic algorithm based optimisation routine PIKAIA (Charbonneau 1995, ApJS, 101, 309). Comparison with predictions of stellar evolution that account for stellar rotation does not result in a unique age, though most stars are best represented by an age of 1-3 Myr. The automated method allows for a detailed determination of the projected rotational velocities. The present day v(r) sin i distribution of the 21 dwarf stars in our sample is consistent with an underlying rotational velocity (v(r)) distribution that can be characterised by a mean velocity of about 160-190 km s(-1) and an effective half width of 100-150 km s(-1). The vr distribution must include a small percentage of slowly rotating stars. If predictions of the time evolution of the equatorial velocity for massive stars within the environment of the SMC are correct (Maeder & Meynet 2001, A&A, 373, 555), the young age of the cluster implies that this underlying distribution is representative for the initial rotational velocity distribution. The location in the Hertzsprung-Russell diagram of the stars showing helium enrichment is in qualitative agreement with evolutionary tracks accounting for rotation, but not for those ignoring vr. The mass loss rates of the SMC objects having luminosities of log L-star/L-circle dot greater than or similar to 5.4 are in excellent agreement with predictions by Vink et al. (2001, A&A, 369, 574). However, for lower luminosity stars the winds are too weak to determine. M accurately from the optical spectrum. Three targets were classified as Vz stars, two of which are located close to the theoretical zero-age main sequence. Three lower luminosity targets that were not classified as Vz stars are also found to lie near the ZAMS. We argue that this is related to a temperature effect inhibiting cooler from displaying the spectral features required for the Vz luminosity class.

Boundary conditions in the simplest model of linear and second harmonic magneto-optical effects

This paper is concerned with linear and nonlinear magneto- optical effects in multilayered magnetic systems when treated by the simplest phenomenological model that allows their response to be represented in terms of electric polarization, The problem is addressed by formulating a set of boundary conditions at infinitely thin interfaces, taking into account the existence of surface polarizations. Essential details are given that describe how the formalism of distributions (generalized functions) allows these conditions to be derived directly from the differential form of Maxwell's equations. Using the same formalism we show the origin of alternative boundary conditions that exist in the literature. The boundary value problem for the wave equation is formulated, with an emphasis on the analysis of second harmonic magneto-optical effects in ferromagnetically ordered multilayers. An associated problem of conventions in setting up relationships between the nonlinear surface polarization and the fundamental electric field at the interfaces separating anisotropic layers through surface susceptibility tensors is discussed. A problem of self- consistency of the model is highlighted, relating to the existence of resealing procedures connecting the different conventions. The linear approximation with respect to magnetization is pursued, allowing rotational anisotropy of magneto-optical effects to be easily analyzed owing to the invariance of the corresponding polar and axial tensors under ordinary point groups. Required representations of the tensors are given for the groups infinitym, 4mm, mm2, and 3m, With regard to centrosymmetric multilayers, nonlinear volume polarization is also considered. A concise expression is given for its magnetic part, governed by an axial fifth-rank susceptibility tensor being invariant under the Curie group infinityinfinitym.

Detection of C8H- and Comparison with C8H toward IRC+10216

We report the detection of new transitions of octatetraynyl (C8H) toward the circumstellar envelope IRC +10 216 using data taken with the 100 m Green Bank Telescope (GBT). In addition, we report five features from the Ku, K, and Q bands that have been identified as transitions of the octatetraynyl anion (C8H-). From a rotational temperature diagram and an assumed source size of 30", we find a total C8H column density of 8(3)×10^12 cm-2 and a rotational temperature of ~13 K. From the five detected transitions of C8H-, we find a total C8H- column density of ~2.1×10^12 cm-2 consistent with a rotational temperature of ~34 K for a total C8H/C8H- column density ratio of ~3.8. This observed C8H/C8H- column density ratio is similar to the theoretical prediction of 3.6, while the observed column densities were lower than that predicted by a factor of ~30. This prompted us to reinvestigate the initial conditions of the circumstellar envelope (CSE) model. The new model results are presented, and they more closely match the C8H and C8H- abundances observed with the GBT. Finally, we use the new CSE model results to predict the abundance of decapentaynyl (C10H), and we compare them with the measured upper limit found from the GBT observations.