Äîêóìåíò âçÿò èç êýøà ïîèñêîâîé ìàøèíû. Àäðåñ îðèãèíàëüíîãî äîêóìåíòà : http://sn.sai.msu.ru/~sil/preprints/counter.pdf Äàòà èçìåíåíèÿ: Fri Jul 19 13:52:26 2002 Äàòà èíäåêñèðîâàíèÿ: Mon Oct 1 19:39:41 2012 Êîäèðîâêà: Ïîèñêîâûå ñëîâà: earth's atmosphere

The Astronomical Journal, 124: 706­721, 2002 August
# 2002. The American Astronomical Society. All rights reserved. Printed in U.S.A.

COMPLEX FORMATION HISTORY OF THE LENTICULAR GALAXIES WITH STELLAR COUNTERROTATION: NGC 4138 AND NGC 45501 V. L. Afanasiev
Special Astrophysical Observatory, Nizhnij Arkhyz, 369167 Russia; vafan@sao.ru

and O. K. Sil`chenko

2

Sternberg Astronomical Institute, Moscow, 119992 Russia; and Isaac Newton Institute of Chile, Moscow Branch; olga@sai.msu.su Received 2001 June 14; accepted 2002 March 24

ABSTRACT Two lenticular galaxies with counterrotating stellar components in their disks have been studied with the Multi-Pupil Fiber Spectrograph of the 6 m telescope (at the Special Astrophysical Observatory). In NGC 4138 the nucleus is certainly chemically distinct, with the central concentration of magnesium enhancement marginally exceeding that of iron. The ionized gas within R = 400 of the nucleus rotates circularly in an inclined ring, perhaps a kind of a circumnuclear polar ring surrounding a compact stellar minibar. NGC 4550 exhibits complex extended structures in its center that can be identified by enhanced levels of magnesium and iron indices; the stellar population of the unresolved nucleus is younger than that of the circumnuclear region. We conclude that although NGC 4550 does not possess a chemically distinct nucleus in the ordinary sense of the term, there are nevertheless clear signatures of a secondary star formation burst confined perhaps to the counterrotating stellar gaseous disk. We argue that the inclination of the counterrotating disk differs from that of the main stellar disk, implying that the two disks are not coplanar. Both galaxies seem to have possessed bars, either actually present or already dissolved. Key words: galaxies: evolution -- galaxies: individual (NGC 4550, NGC 4138) -- galaxies: nuclei -- galaxies: structure

1. INTRODUCTION

It was back in 1992 that we first reported the existence of chemically distinct (metal enriched) stellar nuclei in earlytype disk galaxies (Sil`chenko, Afanasiev, & Vlasyuk 1992). A preliminary statistical analysis implied that as much as half of all Sa or Sb galaxies and at least a quarter of all SOs may possess such nuclei (Sil`chenko 1994). We now know about two dozen well-studied cases. Recently a new FrenchDutch-British project of the two-dimensional spectrograph SAURON has started at the 4.2 m William Herschel Telescope (WHT), and one of the first attempts to map absorption-line indices in the centers of nearby early-type galaxies has led to a discovery of a chemically distinct nucleus in NGC 3384 (de Zeeuw et al. 2002). What may be an origin of the chemically distinct nuclei? We noted that the phenomenon of chemically distinct nuclei is often accompanied by that of multiple kinematical subsystems. For example, we found chemically distinct nuclei in NGC 7331 (Sil`chenko 1999b) and NGC 7217 (Sil`chenko & Afanasiev 2000) where counterrotating stellar subsystems had been reported earlier by Prada et al. (1996) and Merrifield & Kuijken (1994), respectively, as well as in NGC 4826 (Sil`chenko 1996) and NGC 7332 (Sil`chenko 1999a), which contain counterrotating gas (Rubin 1994; Rix et al. 1995; Fisher, Illingworth, & Franx 1994; Plana & Boulesteix 1996). As for a very regular Sb galaxy NGC 2841, we began by finding there a chemi1 Partly based on observations collected with the 6 m telescope (BTA) at the Special Astrophysical Observatory of the Russian Academy of Sciences. 2 Guest investigator of the UK Astronomy Data Centre.

cally distinct nucleus and a circumnuclear polar gaseous disk (Sil`chenko, Vlasyuk, & Burenkov 1997b), followed by a counterrotating stellar component, which we later discovered in its bulge (Afanasiev & Sil`chenko 1999). Circumnuclear polar gaseous disks in noninteracting spiral galaxies are quite a new phenomenon, which might also be related to chemically distinct nuclei; besides NGC 2841, we found them in NGC 6340 (Sil`chenko 2000) and NGC 7217 (Zasov & Sil`chenko 1997; Sil`chenko & Afanasiev 2000). Chemically distinct nuclei might be in some way genetically related to multiple kinematical subsystems. To further explore this idea, we undertook a special study of stellar population properties in the centers of two early-type disk galaxies with counterrotating global stellar disks: NGC 4138 and NGC 4550. Table 1 summarizes the global parameters of these galaxies. The most comprehensive analysis of the regular S0/a galaxy NGC 4138 has been performed by Jore, Broeils, & Haynes (1996), who investigated the H i distribution and kinematics, broadband (I filter) morphology, and differential photometry, and made long-slit kinematic observations of stars and ionized gas, which make up the most impressive part of their work. The above authors found that gas--both neutral and ionized--counterrotates with respect to the major stellar component; however, there is also a counterrotating stellar component, which exhibits disklike velocity dispersion and is concentrated toward the gaseous ring of radius 2200 (1.7 kpc). Also noteworthy are line-of-sight velocity variations along the minor axis shown by the ionized gas within ô500 of the center: even though the whole counterrotating matter is confined to the global galactic disk as Jore et al. (1996) argue, there is a violation of circular 706


NGC 4138 AND NGC 4550
TABLE 1 Global Parameters of the Galaxies Parameter Type (NED) ..................................... R25 (kpc) (LEDA) ............................ B0 (RC3) ......................................... T MB ................................................... ÏB þ V ÷0 (RC3) ............................... T ÏU þ B÷0 (RC3) ............................... T Vr(radio) (LEDA) (km sþ1) .............. Distance, Mpc (LEDA, H0 = 75 km sþ1 Mpcþ1) ............... Inclination (LEDA) (deg)................. P.A.phot (LEDA) (deg) ..................... NGC 4138 SA(r)0+ 6.4 12.08 þ18.94 0.83 0.31 907 14.6 61 150 NCG 4550 SB0 8.1 12.36 þ18.7 0.84 0.38 381 17 (Virgo) 90 178

707

dominated; if the latter is true, the question arises immediately: Why do Rix et al. (1992) see no kinematical signatures of the bulge along the major axis? The very center of the galaxy within R = 300 (Rix et al. 1992) also exhibits a rather unusual kinematical behavior: the stellar component with a zero rotation velocity and a velocity dispersion of only 80 km sþ1 combined with rather rapidly rotating ionized gas. Where is the kinetic energy of stars? Is the system not virialized?
2. OBSERVATIONS AND DATA REDUCTION

character of rotation or of a well-settled planar distribution in the very center. It also remains to carefully verify the results of the photometric analysis by Jore et al. (1996), who claim NGC 4138 (B/D = 0.13) to be a globally diskdominated galaxy despite its very early morphological type. The I-band brightness profile of NGC 4138 reported by Jore et al. (1996) fits well a single exponential law at R > 1000 , and the bulge of the galaxy is very compact and also exponential. At the same time Baggett, Baggett, & Anderson (1998) analyzed the V-band brightness profile of NGC 4138 and found a much more prominent de Vaucouleurs bulge with re = 2200 and an inner cutoff of the exponential disk at R = 1900 . The morphology of the galaxy remains quite unclear. We must add only that the galaxy possesses a ring of H ii regions at R % 2200 (Pogge & Eskridge 1987) while being morphologically unbarred. NGC 4550 is an S0 galaxy seen almost edge-on. Rubin, Graham, & Kenney (1992) have noted in it two counterrotating stellar disks of comparable strength. Rix et al. (1992) carefully analyzed a long-slit major-axis spectrum of this galaxy and showed its counterrotating stellar components to have identical exponential radial brightness distributions at R > 600 and uniformly low stellar velocity dispersion of less than 60 km sþ1. Though the disky nature of both counterrotating stellar components is beyond doubt, the morphology of this galaxy remains ambiguous, as is that of NGC 4138. Rix et al. (1992) and Fisher, Franx, & Illingworth (1996) noted a small de Vaucouleurs bulge seen only inside R % 500 , with B/T estimates ranging from 0.16 to 0.3; Gavazzi et al. (2000) confirmed the predominance of the global disk component (B/T = 0.21) but reported an exponential radial brightness distribution in the bulge in the H filter. However, the sophisticated two-dimensional decomposition of the V-image of NGC 4550 by Scorza et al. (1998) yielded an exponential disk fully embedded into a de Vaucouleurs bulge and everywhere fainter than the latter (the total B/D % 5); similar results were obtained by Baggett et al. (1998), who analyzed a onedimensional V-band brightness profile in terms of exponential and de Vaucouleurs law for the disk and bulge, respectively: the disk was found to be everywhere fainter than the bulge. Finally Simien & Michard (1990) and Caon, Capaccioli, & D'Onofrio (1993) concluded that the brightness profile of the dominant component in NGC 4550, be it bulge or disk, is neither exponential nor de Vaucouleurs; it can be best fitted by the Sersic law with n = 1.7. After all, we do not yet know whether this galaxy is disk- or bulge-

The spectral data that we analyze in this work were obtained with three different integral field spectrographs. Integral field spectroscopy is a rather new approach, which was first suggested by G. Courtes some 15 years ago (see Bacon et al. 1995 for a description of the underlying idea of the instrument). This technique allows one to obtain simultaneously a set of spectra covering a wide spectral interval from an extended sky area, for example, from a central part of a galaxy. A two-dimensional array of microlenses provides a set of micropupils, which are put onto to the entry of a spectrograph. A reduction of a full set of spectra corresponding to individual spatial elements yields a list of continuum and emission-line fluxes, line-of-sight velocities for stars and ionized gas, and equivalent widths of absorption lines, which are usually expressed as indices in the welldefined Lick system (Worthey et al. 1994). The resulting list can be transformed into two-dimensional maps of the parameters mentioned above for the central part of the galaxy under study. This approach offers, besides the wellknown benefits of field observations, a unique opportunity to superpose various two-dimensional distributions without bothering about positioning problems. Two spectral ranges ° are usually exposed, green (4800­5600 A) and red (6200­ ° ). The green spectra are used to calculate the Lick 6900 A indices H, Mg b, Fe 5270, and Fe 5335, which are suitable for determining metallicity, age, and the Mg/Fe ratio of old stellar populations (Worthey 1994), and to cross-correlate the galaxy spectra with the spectrum of a template star, usually of a K0 III­K3 III spectral type, to derive the line-ofsight velocity field for the stellar component and in some cases the map of stellar velocity dispersion. The red spectral range contains the strongest optical emission lines H and [N ii] 6583, which are used to derive the line-of-sight velocity field of ionized gas. The detailed log observations of NGC 4138 and NGC 4550 galaxies made with three different two-dimensional spectrographs is given in Table 2. The observations in 1998 January were made with the old variant of the Multi-Pupil Field Spectrograph (MPFS) of the 6 m telescope (Afanasiev et al. 1996). The instrument operated in two modes with 8 á 12 or 8 á 16 micropupil arrays, supported by 520 á 580 and 1040 á 1160 CCD detectors, respectively, with a spatial scale of 1 > 3 per micro° lens and a reciprocal dispersion of 1.6 A per pixel (the spec° tral resolution was 4­6 A and varied slightly over the frame). We performed wavelength calibration by using separate exposures for the spectrum of the hollow-cathode lamp filled with helium, neon, and argon; the internal accu° racy of linearization was typically 0.25 and 0.1 A in the green and red ranges, respectively. However, we checked the accuracy and absence of systematic shifts in the velocity scale by measuring strong night-sky emission lines [O i] 5577 and [O i] 6300. In this variant of two-dimensional


708

AFANASIEV & SIL`CHENKO
TABLE 2 Two-dimensional Spectroscopy of NGC 4138 and NGC 4550 Exposure (min) 60 40 60 30 45 45 120 Spectral Range ° (A)
00

Date 1998 1998 1999 1999 1999 1999 1999 Jan 21 ....... Jan 23 ....... Jun 6 ......... Jun 9 ......... Dec 12 ....... Dec 15 ....... Feb 21 .......

Galaxy NGC NGC NGC NGC NGC NGC NGC 4138 4550 4550 4550 4138 4138 4550

Configuration BTA/MPFS×CCD 520 á 580 BTA/MPFS×CCD 1040 á 1160 BTA/MPFS×CCD 1024 á 1024 BTA/MPFS×CCD 1024 á 1024 BTA/MPFS×CCD 1024 á 1024 BTA/MPFS×CCD 1024 á 1024 WHT/SAURON×CCD 2000 á 4000

Field 10 1100 1600 1600 1600 1600 3300
00

Seeing 2 2 1 2 3 2 3 > > > > > > > 4 2 7 5 0 1 0

á á á á á á á

16 2100 1500 1500 1500 1500 4100

4700­5400 4200­5600 4200­5600 5800­7200 4200­5600 5800­7200 4800­5400

spectroscopy, we also made a separate exposure of a blank sky area several arcminutes from the galaxy, because we needed sky subtraction to calculate properly the equivalent widths of the absorption lines. We checked the consistency of our index system with the Lick system by observing nine standard stars from the list of Worthey et al. (1994) and found our indices to be in excellent agreement with those tabulated by the above authors, with the mean deviations ° being less than 0.05 A. Later in 1998, a new variant of the MPFS became operational in the prime focus of the 6 m telescope.3 Compared with the old variant, the new MPFS has a larger field of view and wider total spectral range. This became possible because of the use of optical fibers transmitting to the spectrograph slit the light from 16 á 15 square elements of the galaxy image and another 16 fibers transmitting sky background light taken outside the galaxy, thus eliminating the need for separate sky exposures. The size of one spatial element was approximately 100 á 100 . The detector used was 1024 á 1024 TK CCD with a reciprocal dispersion of 1.35 ° ° A per pixel and a rather stable spectral resolution of 5 A. The wavelength calibration was performed in the same way as in the previous variants of the MPFS observations. To calibrate the new MPFS index system onto the standard Lick one, we observed 15 stars from the list of Worthey et al. (1994) during four observational runs. Figure 1a compares the indices calculated directly from the new MPFS spectra with those tabulated in Worthey et al. (1994). The data points can be seen to form tight linear relations, which, however, deviate slightly from bisectors. We have calculated the linear regression formulae to transform our index measurements into the Lick system; the rms scatter of data points ° in Figure 1a about these linear relations is $0.2 A for all four indices considered and is therefore within the observational errors quoted by Worthey et al. (1994). We corrected our index measurements for the broadening due to stellar velocity dispersion, which usually differs substantially from zero in the centers of early-type galaxies. We determined the correction values by artificial Gaussian broadening of the spectrum of the standard star HD 97907 by using various values. We then fitted fourth-order polynomials to the derived dependencies of index corrections on and applied the resulting corrections to the measured index values prior to reducing them into the Lick system. A typical correction ° for stellar velocity dispersion is about 0.4 A in the spectral index at = 200 km sþ1. The third two-dimensional spectrograph whose data we use in this work is the new SAURON instrument operated
3

See http://www.sao.ru/~gafan/devices/mpfs/.

at the 4.2 m William Herschel Telescope of La Palma Observatory (see Bacon et al. 2001 for a detailed description). We retrieved the data for NGC 4550 obtained with this instrument in 1999 February from the open Isaac Newton Group (ING) Archive of the UK Astronomy Data Centre. This instrument has a field of view of 4100 á 3300 and each spatial element has size 0 > 94 á 0 > 94. Sky background is exposed simultaneously with the target object 20 from the center of the galaxy. The data cover the ° spectral range 4800­5400 A with a reciprocal dispersion of ° 1.11­1.21 A, which varies from the left to the right edge of the frame. Neon spectrum is used for comparison and linea° rization is performed to an accuracy of 0.1 A by using a second-order polynomial. The consistency of the resulting index system with the Lick system is checked using seven stars from the list of Worthey et al. (1994). Figure 1b compares the indices measured directly from the SAURON spectra with the values tabulated by Worthey et al. (1994). The two index systems agree to within an rms scatter of 0.3 ° A, and therefore we apply no correction to the index system, based on the SAURON data. This example demonstrates that index systems usually depend on the shape of the continuum rather than on spectral resolution: the spectral reso° lution of SAURON, which is equal to 3.5 A, differs more ° from the standard Lick resolution of 8 A than does that of ° the new MPFS of 5 A, whereas the SAURON index system is closer to the standard one, perhaps because of the flat spectral response of the instrument. The exposure times for the galaxies observed with both variants of the MPFS were chosen to be long enough to achieve signal-to-noise ratios of no less than 60 (per angstrom unit) in the nuclei and 20 near the frame edges; according to Cardiel et al. (1998), this level of signal-tonoise ratio provides an accuracy of line indices ranging from ° ° 0.15 A in the center to 0.6 A for the individual spatial elements in the outer part of the galaxy. To keep the accuracy constant along the radius, we also co-added the element spectra obtained with every two-dimensional spectrograph within concentric circular rings centered on the nuclei and then studied the radial dependencies of the stellar population properties by comparing our azimuthally averaged absorption-line indices to evolutionary synthesis models of old stellar populations by Worthey (1994). We estimate the typical accuracy of our azimuthally averaged indices to be ° better than 0.1 A. For NGC 4138 we used in addition to the twodimensional spectral data some photometric data, namely, the 200 + 160 s exposures made with the Hubble Space Telescope Wide Field Planetary Camera 2 (HST WFPC2) on 1997 July 7, using the F547M, which we retrieved from


Fig. 1a

Fig. 1b Fig. 1.--Calibration of our measured indices into the Lick system by using observations of standard stars: (a) new MPFS (starting from 1999). Different symbols correspond to different observing seasons. (b) 1999 February SAURON data. Dashed lines show quadrant bisectors.


710

AFANASIEV & SIL`CHENKO

Vol. 124

the HST Archive, and broadband BRIK images of the galaxy obtained within the framework of the Ursa Major Cluster Survey of Tully, which we adopted from the CD-ROM collection of the Canadian Astronomy Data Center. All the data, spectroscopic and photometric, except those obtained with the new MPFS, were reduced with the software developed by V. V. Vlasyuk in the Special Astrophysical Observatory (Vlasyuk 1993). Primary reduction of the data obtained with the new MPFS was performed using IDL-based software developed by one of us (V. L. A.). The Lick indices were calculated with our own FORTRAN program and a FORTRAN program of A. Vazdekis.

3. NGC 4138

3.1. Nuclear and Circumnuclear Stellar Populations We measured radial variations of the Lick indices H, Mg b, Fe 5270, and Fe 5335 in the central part of NGC 4138 twice, in 1998 January with the old MPFS operating with a 8 á 12 microlens array and in 1999 December with the new MPFS and a 16 á 15 microlens array. Figure 2 compares two sets of azimuthally averaged measurements within R = 900 . No significant systematic discrepancies between index measurements are in evidence here despite quite different spectrographs and slightly different spectral resolutions. Both sets of measurements yield rather accurate indices,
Fig. 3.--Iron vs. magnesium diagram for azimuthally averaged index measurements in NGC 4138. The data points are taken in galactocentric radial steps of 1 > 3 and 100 for the 1998 January and 1999 December sets, respectively. The models of Worthey (1994) for old stellar populations with [Mg/Fe] = 0 are plotted as a reference; the metallicities of Worthey's models are (top to bottom) +0.50, +0.25, 0.00, þ0.22, þ0.50, þ1.00,þ1.50, and þ2.00.

Fig. 2.--Comparison of two MPFS sets of index measurements for NGC 4138. Two-dimensional measurements are azimuthally averaged over concentric rings centered onto the nucleus.

and one can even note the poorer spatial resolution in 1999 December to show up in the Mg b profile. Figure 3 shows a diagram of hFei versus Mg b [hFei (Fe 5270 + Fe 5335)/2], where we compare our measurements to simple models of old stellar populations by Worthey (1994). The indices are already corrected for the broadening due to stellar velocity dispersion as measured by Jore et al. (1996). One can see that the radial variations of metal indices in NGC 4138 are roughly parallel to model sequences, implying that the stellar populations in the center of NGC 4138 are mildly magnesium-overabundant ([Mg/Fe] % +0.2), and this magnesium-to-iron ratio remains almost constant along the radius. The behavior is similar to that of giant elliptical galaxies (Worthey, Faber, & Gonzalez 1992), although NGC 4138 is more often classified as an S0­Sa galaxy and its bulge is of a very moderate luminosity. Unfortunately, the strong contamination by emission prevents the use of the H index for determining the mean age of stellar population and disentangling age and metallicity effects in this particular case. Nonetheless, strong nucleus-to-bulge variations of mean metallicity show up in the hFei versus Mg b diagram: the nucleus and the bulge are closer to the models with above-solar and subsolar metallicities, respectively. If we assume that the ring at R = 500 ­700 , taken as a representative bulge area, lies beyond the region affected by the nucleus under our moderate seeing quality, the nucleus-to-bulge metal index differences are DMg b = 1.44 ô 0.20 and DhFei = 0.92 ô 0.13, corresponding to a total metallicity difference of 0.6 dex, or a factor of 4, if equal mean ages are assumed for both stellar populations. A relatively younger nucleus compared with the bulge, which is typical for chemically distinct nuclei (Sil`chenko 1999a, 1999b; Sil`chenko & Afanasiev 2000), should imply an even greater metallicity difference. Figure 4 demonstrates two-dimensional maps of the Lick indices H, Mg b, and hFei for the small central region of NGC 4138 according to our MPFS measurements of 1999 December smoothed by a Gaussian with FWHM = 2 > 3.


No. 2, 2002

NGC 4138 AND NGC 4550

711

° Fig. 4.--MPFS index maps for NGC 4138. Green (5000 A) continuum is superposed by isophotes; the second brightness peak northwest of the nucleus is due to a foreground star.

The map of H absorption-line index tells us little about the properties of stellar populations: throughout the entire central region of NGC 4138 it is contaminated by Balmer emission, and the strongest contamination is observed near the major axis of the continuum isophotes. But as far as metal line indices are concerned, both Mg b and hFei maps reveal the chemically distinct nature of the nucleus (the second continuum brightness peak on the major axis northwest of the nucleus and the corresponding minima in all indices are due to a bright blue foreground star). However, the Mg b and hFei indices can be seen to exhibit different behaviors: magnesium is centrally peaked, whereas the maximum of the iron index shows up as a resolved elongated feature. We have already seen similar `` iron-rich circumnuclear disks '' in galaxies with chemically distinct nuclei, e.g., in the S0 galaxy NGC 1023 (Sil`chenko 1999a) or in the Sb galaxy NGC 7331 (Sil`chenko 1999b), but it is for the first time that this

structure is elongated in a direction close to the minor axis of the continuum isophotes. 3.2. Central Kinematics of Stars and Ionized Gas Jore et al. (1996) performed a detailed analysis of the kinematics of stars and ionized (as well as neutral) gas in this galaxy, but this work is based only on observations made with a long slit oriented along the major and minor axes. We used a two-dimensional field spectrograph, and therefore our results shed additional light onto the peculiar kinematics of NGC 4138. Figure 5 shows two-dimensional velocity fields for stars and ionized gas (the latter is based on [N ii] 6583 emissionline measurements). Both velocity fields look quite regular: stars demonstrate a rigid rotation pattern, and the rotation of ionized gas shows marginal local maxima of relative

Fig. 5.--Line-of-sight velocities of the stellar and gaseous components in the center of NGC 4138 (isolines), superposed on the green and red continuum images.


712

AFANASIEV & SIL`CHENKO

Vol. 124

Fig. 6.--[N ii] 6583 emission-line surface brightness distribution (isolines) superposed on the red (6500) continuum image in the center of NGC 4138 (left) and for comparison an HST WFPC2 image at the same scale and with the same orientation (right).

velocity 500 ­600 from the center. In a rough sense, stars and ionized gas counterrotate, as reported by Jore et al. (1996): the receding side of the stellar rotation curve is the one northwest of the center, and that of the gas rotation curve, southeast of the center. But we now report a feature that is of greatest, we would even say crucial, importance: the kinematical major axis (the direction of maximum line-of-sight velocity gradient) of the stellar subsystem does not coincide with this of the ionized gas subsystem, implying that the two systems may rotate in different planes. Averaging over two measurement sets and three radial bins yields for stars a kinematical major-axis orientation of P.A.0,* = 327 ô 3 in the 2 > 5­500 radial interval. The fact that the line of nodes of the galactic global disk inferred from the orientation of the outermost isophotes does not deviate strongly from this direction (330 and 332 according to the Lyon-Meudon Extragalactic Database [LEDA] and to our own photometric measurements reported in x 3.3., respectively) leads us to conclude that stars within R = 300 ­500 rotate axisymmetrically in the global symmetry plane of NGC 4138. As for ionized gas, its kinematical major axis in the same radius range has a position angle of P.A.0,g = 116 ô 4 (according to the line-of-sight velocity fields based on [N ii] 6583 and [O iii] 5007 emission-line measurements) and thus deviates by 35 from the line of nodes. To what plane can we attribute this circumnuclear ionized gas? Or does it experience strong noncircular motions? Figure 6 (left) shows the surface brightness distribution in the [N ii] 6583 emission line. The major axis of the innermost emission-line isophotes has a position angle of P.A. % 113 , which is very close to that of the kinematical major axis [the emission isophote deviation from the line of nodes is even more pronounced in the more outer part at R > 700 , where P.A.(major axis) % 100 ]. The coincidence of the kinematical and photometric major axes is indicative of planar axisymmetric rotation. The central angular rotation velocity of stars, as estimated from the velocity field shown in Figure 5, is ! sin i = 11.3 ô 0.7 km sþ1 arcsecþ1; the maximum central line-of-sight velocity gradient for ionized gas is 12.3 ô 1.6 km sþ1 arcsecþ1. The comparable velocity dispersions of stars and ionized gas ($150 and $120 km sþ1,

respectively; Jore et al. 1996) in the center of NGC 4138 and the fact that the inclination of the gaseous disk implied by the ellipticity of the emission-line isophotes shown in Figure 6(left) does not differ strongly from that of the global stellar rotation plane lead us to interpret the coincidence of the maximum central line-of-sight velocity gradients for stellar and ionized gas subsystems as conclusive evidence of circumnuclear ionized gas being confined to and rotating circularly in a plane inclined with respect to the global symmetry plane of the galaxy. This point can be illustrated by an F547W-band image of the central part of NGC 4138 obtained with the HST WFPC2 (Fig. 6, right). One can see a straight dust lane crossing the center of the galaxy and extending to a radius of at least 300 ; the orientation of this dust lane is P.A. = 115 --the same as that of the kinematical and photometric major axes of the gaseous emission. The circumnuclear gaseous disk appears to contain a significant amount of dust. The overall appearance of Figure 6 (right) gives the impression that the dust lane crosses the central elongated structure at a right angle. If so, continuum isophotes should be strongly twisted in the center. Or perhaps this might be a false impression caused by the complex distribution of dust? To answer this question, in the next subsection we analyze the photometric structure of NGC 4138 3.3. Photometric Structure To investigate the photometric structure of NGC 4138, we used BRIK images of the galaxy taken from the survey of the Ursa Major Cluster of galaxies by Tully (CD-ROM collection at Harvard). The images were taken under good seeing conditions, FWHM* = 0 > 8/0 > 9, and their analysis can reveal a number of interesting details. Figure 7 shows the variations of isophotal characteristics along the radius. One can see that both the position angle of the major axis and the ellipticity remain more or less constant at R > 2500 , and we therefore may adopt these asymptotic values (the position angle of the line of nodes of P.A.0 % 150 and the inclination of about 50 , which corresponds to an ellipticity of 0.36) for the entire global stellar disk. Inside R = 2500 meas-


No. 2, 2002

NGC 4138 AND NGC 4550

713

Fig. 8.--Residual K-band image of NGC 4138 after subtraction of the global exponential disk. The image is deprojected to the face-on view with the line of nodes aligned horizontally.

Fig. 7.--Radial variations of isophote parameters in NGC 4138

urements made in different passbands are somewhat discrepant in the systematic sense; these data must be obviously affected by the large amount of dust, with the K-band measurements being the least sensitive. Strong variations of the isophote parameters are observed inside R = 300 . The blue-band isophote ellipticity increases toward the nucleus and we may interpret this rise as a dust effect. However, the major axis verges to the north-south direction consistently in all filters, implying that this turn cannot be interpreted in terms of clumpy absorption. The position angle of the circumnuclear continuum isophotes is close to P.A. % 13 , leading us to conclude that the galaxy should have a compact circumnuclear minibar or a minidisk around which a quasi-polar gaseous ring rotates (