CDM, Feedback and the Hubble Sequence

We have performed TreeSPH simulations of galaxy formation in a standard ΛCDM cosmology, including effects of star formation, energetic stellar feedback processes and a meta-galactic UV field, and obtain a mix of disk, lenticular and elliptical galaxies. The disk galaxies are deficient in angular momentum by only about a factor of two compared to observed disk galaxies. The stellar disks have approximately exponential surface density profiles, and those of the bulges range from exponential to r ^1/4, as observed. The bulge-to-disk ratios of the disk galaxies are consistent with observations and likewise are their integrated B-V colours, which have been calculated using stellar population synthesis techniques. Furthermore, we can match the observed I-band Tully-Fisher (TF) relation, provided that the mass-to-light ratio of disk galaxies is (M/L _I) ≃ 0.6–0.7. The ellipticals and lenticulars have approximately r ^1/4 stellar surface density profiles, are dominated by non-disklike kinematics and flattened due to non-isotropic stellar velocity distributions, again consistent with observations. This revised version was published online in September 2006 with corrections to the Cover Date.     

Evolutionary processes in protoplanetary accretion disks: the propagation of axisymmetric shock waves

In this paper we investigate both the global and the local hydrodynamics of axisymmetric accretion disks around young stellar objects under the simultaneous action of viscosity, self-gravity and pressure forces. For simplicity, we take for the global model a polytropic equation of state, make the infinitely thin disk approximation and characterize the surface density and temperature profiles in the disk as power laws in the radial distance r from the protostar. We solve the problem of the general density profile of a Keplerian disk showing that self-gravity could not be an important factor for the fast formation of the rocky cores of giant gaseous planets in our solar system. Under the hypothesis that the unperturbed rotation curve of the disk is nearly Keplerian throughout the radial extent, we can estimate with our polytropic model a lower limit for the resulting masses M_d(r) of stable disks up to 100 AU. These masses are in the range of the so-called minimum mass solar nebular (_d/M_s ≈ 0.01–0.02).By adopting a simplified viscosity model, where the height-integrated turbulent dynamical viscosity ν is a function of the surface density σ like η ∝ σ^Γ, we derive in the local shearing sheet model linearized evolution equations for small density perturbations describing both a diffusion process and the propagation of acoustic density waves. We solve a special initial value problem and calculate the appropriate Green's function. The analytical solutions so obtained describe in the case Γ < 0 the successive formation of quasi-stationary ring-shaped density structures in a disk with a definite mode of maximum instability, whereas in the case Γ > Γ_c the density wave equation describes the propagation of an “overstable” ring-shaped acoustic density wavelet to the outer ranges of the accretion disk. Whereas the group velocity of the wave packet is subsonic, the phase velocities of individual wave crests in the wave packet are supersonic. The mode of maximum instability, the growth rate and the number of growing waves in the wavelet are controlled by Γ and α. Our present knowledge concerning turbulent viscosity in protoplanetary disks is not sufficient to decide whether or not the case Γ > Γ_c is realized.The suggested structuring processes in the linear theory should initiate in the non-linear regime the formation of narrow ring-shaped density shock waves moving through the protoplanetary disk. These non-linear waves could produce extremely spatially and temporally heterogeneous temperature regions in the disk. We speculate that ring-shaped density waves, excited by inner boundary conditions and which have dominated the disk's evolution at early times, are responsible both for the fast growth of dust to planetesimals and at least for the rapid accretion of the rocky cores of giant gaseous planets in the protoplanetary accretion disk (shock wave trigger hypothesis). We derive provisional scaling rules for planetary systems regarding the spacing of orbits as a function of the mass ratio of the protoplanetary disk to the protostar. However, further analytical work and linear as well as nonlinear numerical simulations of density waves excited by inner boundary conditions are needed to consolidate the results and speculations of our linear wave mechanics in the future.     

Spectroscopic study of the peculiar galaxy UGC 5119

We present the results of our study of the stellar kinematics in the elliptical galaxy UGC 5119, which has previously been suspected to be a polar-ring galaxy. We have detected a rapidly rotating disk in the central region (r ≤ 3.2 kpc) of the galaxy’s main body and found a radial velocity gradient along its minor axis (in the putative ring). We conclude that UGC 5119 is a medium-luminosity elliptical galaxy with a rapidly rotating disk component and a stellar (probably polar) ring. We have calculated the Lick indices of the Hβ, Mggb, Fe 5270, and Fe 5335 absorption lines and compared them with evolutionary synthesis models. Differences in the [Mg/Fe] ratios, metallicities, and ages of the stars have been found: the young stellar population with a solar [Mg/Fe] ratio and a high metallicity dominates in the circumnuclear region (r ≤ 1 kpc), while the old one with a low metal abundance dominates in the ring.     

Properties of spherical galaxies and clusters with an NFW density profile

Using the standard dynamical theory of spherical systems, we calculate the properties of spherical galaxies and clusters whose density profiles obey the universal form first obtained in high-resolution cosmological N -body simulations by Navarro, Frenk & White (NFW). We adopt three models for the internal kinematics: isotropic velocities, constant anisotropy and increasingly radial OsipkovMerritt anisotropy. Analytical solutions are found for the radial dependence of the mass, gravitational potential, velocity dispersion, energy and virial ratio and we test their variability with the concentration parameter describing the density profile and amount of velocity anisotropy. We also compute structural parameters, such as half-mass radius, effective radius and various measures of concentration. Finally, we derive projected quantities, the surface mass density and line-of-sight as well as aperture-velocity dispersion, all of which can be directly applied in observational tests of current scenarios of structure formation. On the mass scales of galaxies, if constant mass-to-light is assumed, the NFW surface density profile is found to fit HubbleReynolds laws well. It is also well fitted by Sérsic R ^{1/ m} laws, for     but in a much narrower range of m and with much larger effective radii than are observed. Assuming in turn reasonable values of the effective radius, the mass density profiles imply a mass-to-light ratio that increases outwards at all radii.     

The red giant branch in the Tycho-2 catalogue

Based on multicolor photometry from the 2MASS and Tycho-2 catalogues, we have produced a sample of 38 368 branch red giants that has less than 1% of admixtures and is complete within 500 pc of the Sun. The sample includes 30 671 K giants, 7544Mgiants, 49 C giants, and 104 suspected supergiants or S stars. The photometric distances have been calculated for K, M, and C stars with an accuracy of 40%. Tycho-2 proper motions and PCRV radial velocities are used to analyze the stellar kinematics. The decrease in the stellar distribution density with distance from the Galactic equator approximated by the barometric law, contrary to the Besanconmodel of the Galaxy, and the kinematic parameters calculated using the Ogorodnikov-Milne model characterize the overwhelming majority of the selected K and M giants as disk stars with ages of more than 3 Gyr. A small number of K and M giants are extremely young or, conversely, thick-disk ones. The latter show a nonuniform distribution in the phase space of coordinates and velocities, arguing against isothermality and full relaxation of the disk and for the theory of dynamical streams or superclusters. The spatial distribution and kinematics of the selected C stars force us to consider them as asymptotic branch giants with masses of more than 2M _⊙ and ages of less than 2 Gyr probably associated with the Gould Belt. The offset of the Sun above the Galactic equator has been found from the distribution of stars to be 13 ± 2 pc, which coincides with the previously obtained value for the clump red giants.     

Gravitational stability and dynamical overheating of galactic stellar disks

Stellar velocity dispersion data at galactocentric distance of two disk radial scale lengths (R = 2h), available in the literature allowed us to determine the upper limits of disk local surface densities at a given R and (by extrapolation) total masses of disks proceeding from the marginal gravitational stability condition. A comparison of the obtained disk masses with the photometric estimates based on the stellar population models indicates the absence of strong dynamical overheating inmost spiral galaxies and hence the absence of significant major merging events, which were able to heat dynamically the inner parts of disks. The same conclusion is valid for some of S0 galaxies. However, a significant part of the latter possesses stellar velocity dispersion, which exceeds the threshold value needed for gravitational stability. Dynamically overheated disks occur both among paired and isolated galaxies. Disk to total mass ratios within R = 4h found for marginally stable disks in most cases lie in the range 0.5–0.8 with the absence of the clearly defined correlation of this ratio with color index or morphological type.     

Structure and stability of rotating fluid disks around massive objects. I. Newtonian formulation

In this paper we have presented a very general class of solutions for rotating fluid disks around massive objects (neglecting the self gravitation of the disk) with density as a function of the radial coordinate only and pressure being nonzero. Having considered a number of cases with different density and velocity distributions, we have analysed the stability of such disks under both radial and axisymmetric perturbations. For a perfect gas disk with γ= 5/3 the disk is stable with frequency (MG/r³)^1/2 for purely radial pulsation with expanding and contracting boundary. In the case of axisymmetric perturbation the critical γ_c for neutral stability is found to be much less than 4/3 indicating that such disks are mostly stable under such perturbations. On leave of absence from Government College, Jagdalpur 494005.     

旋涡星系颜色梯度的研究进展

旋涡星系的颜色梯度反映了其星族构成沿径向的分布,包含了星系恒星形成历史的信息.因此,对旋涡星系颜色梯度的研究有助于理解星系的形成和演化过程.大部分旋涡星系存在负的颜色梯度,其主要原因是旋涡星系存在星族梯度.颜色梯度与星系的面亮度之间存在内禀的相关,表明质量面密度在星系的形成和演化过程中具有重要作用.     

Coronal densities and magnetic fields from K-coronameter and type IV radio burst data

From K-coronameter data we have obtained an electron density profile above the active region responsible for the Type IV burst observed on 14 September 1966. If the observed frequency cutoff in the burst's spectrum is caused by the Razin effect, then the coronal electron density may be derived from the intensity variation in the burst as it propagates outwards from the Sun. We show that the electron density profiles obtained from K-coronameter data (appropriate to 1.125 <r/R < 2.0) and from the radio data (2.2< r/R < 2.5) form a continuous distribution. We conclude that the cutoff is due to the Razin effect, and that radiation in the burst is due to relativistic electrons having a steep inverse power-law energy distribution. From the electron density profile derived from the radio data, we find that the coronal magnetic field was 0.26 G at r/R = 2.2.     

The extended rotation curve and the dark matter halo of M33

We present the 21-cm rotation curve of the nearby galaxy M33 out to a galactocentric distance of 16 kpc (13 disc scalelengths). The rotation curve keeps rising out to the last measured point and implies a dark halo mass ≳5×10¹⁰ M_⊙. The stellar and gaseous discs provide virtually equal contributions to the galaxy gravitational potential at large galactocentric radii, but no obvious correlation is found between the radial distribution of dark matter and the distribution of stars or gas.
Results of the best fit to the mass distribution in M33 picture a dark halo which controls the gravitational potential from 3 kpc outward, with a matter density which decreases radially as R ^−1.3. The density profile is consistent with the theoretical predictions for structure formation in hierarchical clustering cold dark matter (CDM) models, and favours lower mass concentrations than those expected in the standard cosmogony.     

A calculation of structural parameters of spherical stellar systems by means of characteristic functions method. I. Theory

To calculate structural parameters of stellar systems such as an effective radius and central space (or surface) density, the method of characteristic functions is suggested. The characteristic function of the system is a Fourier image of their normalized space density profile f₃(r). In the case of spherical symmetry the probability distribution of r (Q₃(r) = (3/a³)r²f₃(r)) and its orthogonal projections have the same characteristic functions. This fact is used to calculate the effective radii of a few star cluster models (King law, Plummer model and Gausian profile). It is shown, that the characteristic function for King law clusters tends to a finite generalised function if the concentration parameter c is large. The expression for the effective radius (at c ≫ 1) is given. The formula of the effective radius in the Plummer model as well as the relation between the one-dimensional central velocity dispersion and the root mean square velocity are obtained. It is shown, that in the Gaussian model and for King law clusters the effective radius (half-mass visual radius) can differ from the effective (harmonic) radius a few times. This fact should be taken into account in estimating the mass-to-light ratio from the virial mass of such systems using the King radius.     

Report on the physical characteristics of Vaidya-Tikekar's exact relativistic model for a superdense star

We report here the results of our examination of the physical properties of Vaidya-Tikekar's model for a relativistic star. Full details will be published elsewhere. The analysis yields a strong indication that the model is stable with respect to infinitesimal radial oscillations. We find that the adiabatic speed of sound is smaller than the speed of light everywhere inside the fluid sphere if the radius of the sphere is larger than 1.46 times its Schwarzschild radius. We also find that the fluid must necessarily be supraluminal somewhere if the speed of sound is decreasing outwards close to the center. We further find that the strong energy condition is fulfilled everywhere if it is fulfilled at the center. Since the ratio of the pressure p and the density ⋅ is decreasing outwards, this indicates that the temperature gradient is negative. We also find that the relativistic adiabatic index is larger than two. Demanding the fluid to be causal, and taking the pressure and the density to be somewhere given by 7.4 ⋅ 10³³ dynes/cm³ and 5.1 ⋅ 10¹⁴ g/cm³, we calculate the maximum mass of the fluid sphere to be 3 solar masses.     

Estimating the dark halo mass from the relative thickness of stellar disks

We analyze the relationship between the mass of a spherical component and the minimum possible thickness of stable stellar disks. This relationship for real galaxies allows the lower limit on the dark halo mass to be estimated (the thinner the stable stellar disk is, the more massive the dark halo must be). In our analysis, we use both theoretical relations and numerical N-body simulations of the dynamical evolution of thin disks in the presence of spherical components with different density profiles and different masses. We conclude that the theoretical relationship between the thickness of disk galaxies and the mass of their spherical components is a lower envelope for the model data points. We recommend using this theoretical relationship to estimate the lower limit for the dark halo mass in galaxies. The estimate obtained turns out to be weak. Even for the thinnest galaxies, the dark halo mass within four exponential disk scale lengths must be more than one stellar disk mass.     

Do disk galaxies with abnormally low mass-to-light ratios exist?

We have performed photometric B, V, and R observations of nine disk galaxies that presumably have abnormally low total mass-to-light (M/L) ratios for given color indices. Our data on surface photometry are used to analyze the possible causes of anomalous M/L estimates. In many cases, these can be the result of errors in photometry or rotational velocity determination but can also reflect the real peculiarities of the stellar composition of galaxies. Comparison of the photometric and dynamical disk mass estimates obtained by analyzing the rotational velocities shows that low M/L values for a given color index are probably real for some of the galaxies. This is primarily true of NGC 4826 (Sab), NGC 5347 (Sab), and NGC 6814 (Sb). The small number of such galaxies suggests that the stellar initial mass function is universal. However, a small fraction of galaxies probably may have a non-typical mass function “depleted” in low-mass stars. Such galaxies require a more careful study.     

The pressure of an equilibrium interstellar medium in galactic disks

Based on an axisymmetric galactic disk model, we estimate the equilibrium gas pressure P/k in the disk plane as a function of the galactocentric distance R for several galaxies (MW, M33, M51, M81, M100, M101, M106, and the SMC). For this purpose, we solve a self-consistent system of equations by taking into account the gas self-gravity and the presence of a dark pseudo-isothermal halo. We assume that the turbulent velocity dispersions of the atomic and molecular gases are fixed and that the velocity dispersion of the old stellar disk corresponds to its marginal stability (except for the Galaxy and the SMC). We also consider a model with a constant disk thickness. Of the listed galaxies, the SMC and M51 have the highest pressure at a given relative radius R/R ₂₅, while M81 and the Galaxy has the lowest pressure. The pressure dependence of the relative molecular gas fraction confirms the existence of a positive correlation between these quantities, but it is not so distinct as that obtained previously when the pressure was estimated very roughly. This dependence breaks down for the inner regions of M81 and M106, probably because the gas pressure has been underestimated in the bulge region. We discuss the possible effects of factors other than the pressure affecting the relative content of the molecular gas in the galaxies under consideration.     

Evolution of the Radial Abundance Gradient and Cold Gas of the Galactic Disk

In order to understand the forming mechanism of the radial abun- dance gradient of the Galactic disk and the evolution of cold gas, we have con- structed a chemical evolution model of the Galactic disk, in which the star for- mation law concerned with molecular hydrogens is adopted, and the evolution of mass surface density is calculated for the molecular and atomic hydrogens separately, then the model predictions and the observed radial distributions of some physical quantities are compared. The result indicates that the model prediction is sensitive to the adopted infall timescale, the model which adopts the star formation law concerned with the molecular hydrogens can agree well with the major observed properties of the Galactic disk, especially can obtain naturally the radial oxygen abundance gradient of the Galactic disk, and the radial surface density profile of cold gas. The assumption of instantaneous or non-instantaneous recycling approximation has a small effect on the evolution of cold gas, especially in the case of rather low gas density.     

Mechanisms of the vertical secular heating of a stellar disk

We investigate the nonlinear growth stages of the bending instability in stellar disks with exponential radial density profiles. We found that the unstable modes are global (the wavelengths are larger than the disk scale lengths) and that the instability saturation level is much higher than that following from a linear criterion. The instability saturation time scales are of the order of one billion years or more. For this reason, the bending instability can play an important role in the secular heating of a stellar disk in the z direction. In an extensive series of numerical N-body simulations with a high spatial resolution, we were able to scan in detail the space of key parameters (the initial disk thickness z₀, the Toomre parameter Q, and the ratio of dark halo mass to disk mass M_h/M_d). We revealed three distinct mechanisms of disk heating in the z direction: bending instability of the entire disk, bending instability of the bar, and heating on vertical inhomogeneities in the distribution of stellar matter.     

Magnetic star-disk interaction in classical T Tauri systems

We simulate the impact of a dipolar stellar magnetic field rooted in a classical T Tauri star on the accretion disk and the halo above using a 2.5D finite difference code. The gas is assumed resistive, and inside the disk accretion is driven by a Shakura-Sunyaev-type eddy viscosity. The rotational shear between the star and the Keplerian disk causes the magnetic field to be wound up and stretched outwards, away from the star. Part of the field lines open and an outflow is launched. Direct disk disruption by the Lorentz force only occurs for sufficient field strength. For our model system with a solar-mass central star, an accretion rate of 10^-7M⊙/a, and a viscosity parameter α_SS=0.01, a field strength of 1 kG, measured at the poles on the surface of the star, was found insufficient for disk disruption.     

Minor Mergers of Galaxies: Theory and Observations

We compare results from numerical simulations with observations of edge-on galaxies interacting/merging with a small companion (Schwarzkopf and Dettmar,2000), hereafter S&D00). Observations show a clear influence of the merging and interacting process on disk scale parameters h (radial scalelength), z ₀ (vertical scalelength) and their ratio (h/z ₀), leading to a heating and thickening of the stellar disk. Our numerical simulations show the same behaviour but differ significantly in the magnitude of the change of the disk scale parameters. This revised version was published online in September 2006 with corrections to the Cover Date.