The expulsion and retention of dust grains by galactic discs

We have constructed a numerical model of a galaxy that consists of a stellar, gas and dust disc imbedded within a dark halo. We have used this model to assess the radiation, gravitational and viscous forces on dust grains and to trace their motion through the interstellar medium over a period of 10⁹ yr. We conclude that the disc opacity is a crucial factor in understanding the motion of the grains. Large grains (≈0.1 μm) with low disc opacity will lead to dust expulsion from the stellar disc, while high opacity leads to dust retention. Reasonable disc opacities lead to the recycling of the larger grains from the outer to the inner regions of the galaxy. The larger grains travel at higher velocities than small grains (0.01−0.001 μm), and so the smaller grains remain relatively close to their formation sites. Dust can 'leak' out over the entire surface of the disc because of the imbalance of radiation and gravitational forces. The dust is dynamically coupled to the gas and so although the gas lags behind the dust it is carried along with it. This explains the close correlation between the far-infrared emission from dust and the gas column density. We use a simple analytical model to show how the dust mass of a galaxy may evolve with time and how a significant fraction (90 per cent) of the total dust mass produced may have been expelled into the intergalactic medium.     

An excursion set model for the distribution of dark matter and dark matter haloes

A model of the gravitationally evolved dark matter distribution, in the Eulerian space, is developed. It is a simple extension of the excursion set model that is commonly used to estimate the mass function of collapsed dark matter haloes. In addition to describing the evolution of the Eulerian space distribution of the haloes, the model allows one to describe the evolution of the dark matter itself. It can also be used to describe density profiles, on scales larger than the virial radius of these haloes, and to quantify the way in which matter flows in and out of Eulerian cells. When the initial Lagrangian space distribution is white noise Gaussian, the model suggests that the Inverse Gaussian distribution should provide a reasonably good approximation to the evolved Eulerian density field, in agreement with numerical simulations. Application of this model to clustering from more general Gaussian initial conditions is discussed at the end.     

Numerical simulations of protostellar encounters — II. Coplanar disc–disc encounters

It is expected that an average protostar will undergo at least one impulsive interaction with a neighbouring protostar whilst a large fraction of its mass is still in a massive, extended disc. Such interactions must have a significant impact upon the evolution of the protostars and their discs.   We have carried out a series of simulations of coplanar encounters between two stars, each possessing a massive circumstellar disc, using an SPH code that models gravitational, hydrodynamic and viscous forces. We find that during a coplanar encounter, disc material is swept up into a shock layer between the two interacting stars, and the layer then fragments to produce new protostellar condensations. The truncated remains of the discs may subsequently fragment; and the outer regions of the discs may be thrown off to form circumbinary disc-like structures around the stars. Thus coplanar disc–disc encounters lead efficiently to the formation of multiple star systems and small- N clusters, including substellar objects.     

Non-stationary pair plasma in a pulsar magnetosphere and the two-stream Instability

We test a new emission mechanism in pulsar magnetospheres, eventually responsible in part for the high level of observed radio radiation. This is carried out by comparing the efficiency of the two-stream instability of Langmuir waves in a pulsar emission region, where the stationary and non-stationary characters of pair plasma outflows produced in the gap region are characterized by two different time-scales. On the shorter time-scale, the Ruderman &38; Sutherland 'sparking' phenomenon leads to the creation of pair plasma clouds, in motion along magnetic field lines, that contain particles with a large spectrum of momenta. The overlapping of particles with different energies produced in successive clouds results in an efficient 'two stream'-like instability. This effect is a consequence of the non-stationary character of the pair plasma produced in the gap region, just above the magnetic poles of the neutron star. On a long time-scale, resulting pair plasma outflows in pulsar magnetospheres can be treated as stationary. In this case, the instability which results from interaction between existing primary beam particles and the pair plasma is negligible, whereas the instability owing to interaction between electrons and positrons of the pair plasma itself, and more precisely to their relative drift motion along curved magnetic field lines, is effective. We derive characteristic features of the triggered instability, using specific distribution functions to describe either particles in the assembly of clouds or relative drifting of electrons and positrons in these same plasma clouds. Although linear and local, our treatment suggests that non-stationary effects may compete with, or even dominate over, drifting effects in parts of pulsar emission regions.     

Double detonations at the core–envelope boundary in Type Ia supernovae

An off-centre detonation propagating near the interface between a C–O core and a He envelope in a Type Ia supernova explosion is modelled as a steady two-dimensional similarity solution at a plane interface. We assume that in both regions the energy release occurs in an infinitely thin detonation, which produces material in nuclear statistical equilibrium (NSE) in He and in nuclear statistical quasi-equilibrium in C–O. An α-network is then used to determine the effect of the associated rarefaction wave in the C–O on the final abundance of intermediate elements. We find that, although there is a significant effect, the rarefaction is not strong enough to quench the reactions and prevent the C–O from burning to NSE.     

A timing formula for main-sequence star binary pulsars

In binary radio pulsars with a main-sequence star companion, the spin-induced quadrupole moment of the companion gives rise to a precession of the binary orbit. As a first approximation one can model the secular evolution caused by this classical spin-orbit coupling by linear-in-time changes of the longitude of periastron and the projected semi-major axis of the pulsar orbit. This simple representation of the precession of the orbit neglects two important aspects of the orbital dynamics of a binary pulsar with an oblate companion. First, the quasiperiodic effects along the orbit, owing to the anisotropic 1/ r ³ nature of the quadrupole potential. Secondly, the long-term secular evolution of the binary orbit, which leads to an evolution of the longitude of periastron and the projected semi-major axis, which is non-linear in time.   In this paper a simple timing formula for binary radio pulsars with a main-sequence star companion is presented which models the short-term secular and most of the short-term periodic effects caused by the classical spin-orbit coupling. I also give extensions of the timing formula that account for long-term secular changes in the binary pulsar motion. It is shown that the short-term periodic effects are important for the timing observations of the binary pulsar PSR B1259–63. The long-term secular effects are likely to become important in the next few years of timing observations of the binary pulsar PSR J0045–7319. They could help to restrict or even determine the moments of inertia of the companion star and thus probe its internal structure.   Finally, I reinvestigate the spin-orbit precession of the binary pulsar PSR J0045–7319 since the analysis given in the literature is based on an incorrect expression for the precession of the longitude of periastron. A lower limit of 20° for the inclination of the B star with respect to the orbital plane is derived.     

武汉台重力长期变化分析

通过对武汉台超导重力仪９年多连续观测的重力长期变化资料分析,得到长周期重力潮汐参数、该频段气压回归因子Ｃ、极移重力效应因子Ｐ和长周期气压回归因子Ｃ１分别为：Ｃ：－０．３８３（±０．０１４）×１０－８ｍ·ｓ－２／１００Ｐａ;Ｍｍ：１．１１７３（±０．０７２０）,－１°．１７８７（±３°．６９４３）;Ｍｆ：１．１４３２（±０．０４８４）,５°．２２３５（±２°．４２５３）;Ｍｔｍ：１．２２７６（±０．２２４５）,－１７°．５６４８（±１０°．４７５８）;Ｃ１：－０．３１２（±０．０１５）×１０－８ｍ·ｓ－２／１００Ｐａ;Ｐ：１．９０４７（±０．０６９５）．重力长期变化在Ｃｈａｎｄｌｅｒ周期附近的振幅约为３．２６×１０－８ｍ·ｓ－２。     

On the solar cycle dependence of winds and planetary waves as seen from mid-latitude D1 LF mesopause region wind measurements

At the Collm Observatory of the University of Leipzig LF D1 low-frequency total reflection nighttime wind measurements have been carried out continuously for more than two decades. Using a multiple regression analysis to derive prevailing winds, tides and the quasi-2-day wave from the half-hourly mean values of the horizontal wind components, monthly mean values of mesopause wind parameters are obtained that can be analysed with respect to long-term trends and influences of solar variability. The response of the prevailing wind to the 11-year solar cycle differs throughout the year. While in winter no significant correlation between the zonal prevailing wind and solar activity is found, in spring and summer a negative correlation between the TWC can be seen from the measurements. This is connected with stronger vertical gradients of the zonal prevailing wind during solar maximum than during solar minimum. Since the amplitude of the quasi-2-day wave is dependent on the zonal mean wind vertical gradient, this is connected with a positive correlation between solar activity and quasi-two-day wave activity.Paper Presented at the Second IAGA/ICMA (IAMAS) Workshop on Solar Activity Forcing of the Middle Atmosphere, Prague, August 1997     

Parametrization of the precipitation in the Northern Hemisphere and its verification in Mexico

To improve results in monthly rainfall prediction, a parametrization of precipitation has been developed. The thermodynamic energy equation used in the Adem thermodynamic model (ATM) and the Clausius and Clapeyron equation, were used to obtain a linear parametrization of the precipitation anomalies as a function of the surface temperature and the 700 mb temperature anomalies. The observed rainfall in Mexico over 36 months, from January 1981 to December 1983, was compared with the results obtained of the heat released by condensation, which is proportional to precipitation, using our theoretical formula, and those obtained using a statistical formula, which was derived for the ATM using 12 years of hemispheric real data. The verification using our formula in Mexico, showed better results than the one using the statistical formula.