The Discovery of the Most Metal-Rich White Dwarf: Composition of a Tidally Disrupted Extrasolar Dwarf Planet

Abstract: Cool white dwarf stars are usually found to have an outer atmosphere that is practically pure in hydrogen or helium. However, a small fraction have traces of heavy elements that must originate from the accretion of extrinsic material, most probably circumstellar matter. Upon examining thousands of Sloan Digital Sky Survey spectra, we discovered that the helium-atmosphere white dwarf SDSS J073842.56+183509.6 shows the most severe metal pollution ever seen in the outermost layers of such stars. We present here a quantitative analysis of this exciting star by combining high S/N follow-up spectroscopic and photometric observations with model atmospheres and evolutionary models. We determine the global structural properties of our target star, as well as the abundances of the most significant pollutants in its atmosphere, i.e., H, O, Na, Mg, Si, Ca, and Fe. The relative abundances of these elements imply that the source of the accreted material has a composition similar to that of Bulk Earth. We also report the signature of a circumstellar disk revealed through a large infrared excess in JHK photometry. Combined with our inferred estimate of the mass of the accreted material, this strongly suggests that we are witnessing the remains of a tidally disrupted extrasolar body that was as large as Ceres.

\section{INTRODUCTION}

The vast majority of stars, some 97\% of them, start their lives with
masses less than about 8 \msun. After going through mass loss episodes
in the red giant phases, these stars ultimately run out of
thermonuclear fuel and end up as compact objects known as white
dwarfs. The white dwarf stars are characterized by masses around
0.6 \msun~ and dimensions comparable to that of the Earth, and they
slowly cool off over periods of billions of years by radiating away
the remaining thermal energy left in their core. Most of them are
composed mainly --- more than $99\%$ of their mass --- of carbon and
oxygen, the products of hydrogen and helium nuclear burning. However,
the processes of nuclear fusion and mass loss do not destroy or
evaporate 100$\%$ of the hydrogen and helium that were initially
present in the star at birth. Since the surface gravity of a white
dwarf is extremely high (log $g \sim$ 8), light elements such as
hydrogen and helium float to the surface in enough quantities to form
an optically thick photosphere, while heavier elements sink rapidly
out of sight. It is this efficient gravitational separation mechanism
that is responsible for the high purity of the hydrogen or helium
found in the atmospheres of most white dwarfs.

Traces of heavy elements are sometimes detected spectroscopically in
the atmospheres of cool white dwarfs where radiative levitation or
residual stellar winds no longer operate. Since their gravitational
settling timescales are much shorter than the evolutionary cooling
time, these elements cannot be primordial and must have accreted
relatively recently \citep[metals sink on timescales of at most a few
million years, while white dwarfs cool for several billion years,
][]{paquette}.

Passage through clouds in the interstellar medium (ISM) was, until
recently, the commonly accepted scenario to explain the atmospheric
pollution of these white dwarfs. However, two weaknesses in this
hypothesis have been 1) the amount of hydrogen accreted onto
helium-atmosphere stars is found to be several orders of magnitude
lower than that of metals \citep{wolff,dufour07} despite the fact that
the ISM is mostly constituted of hydrogen, and 2) investigations of
the galactic positions and kinematics of white dwarfs with traces of
metals all failed to show any evidence of an interaction with the ISM
\citep{aannestad,kilic07,farihi10}.

The last decade has seen the rise of a new paradigm as observations at
infrared wavelengths revealed the presence of dusty disks orbiting the
most metal-rich white dwarfs
\citep{becklin,kilic05,kilic06,vonhippel,jura07a,jura07b,farihi09}. It
now appears, in all likelihood, that the heavy elements in most, if
not all, polluted cool white dwarfs are accreted from this orbiting
reservoir whose origin is explained by the tidal disruption of one or
many orbitally perturbed asteroids
\citep{debes,jura03,jura06,jura08}. This alternative model also
naturally explains the lack of correlation with the ISM and the small
amount of hydrogen in the accreting material.

Traces of metals in the atmospheres of cool helium-rich white dwarfs
are usually revealed only by the presence of the broad Ca II H and K
absorption lines in the optical \citep[see][and references
therein]{dufour07}, and very few stars show other heavy elements
\citep [note that observations in the ultraviolet part of the
electromagnetic spectrum, unfortunately available only for a few
objects, see][reveal traces of a few more elements, with the strongest
transitions being due to C, Mg, Si and Fe]{wolff,friedrich,koester00}.

In a few rare cases, however, lines from a handful of elements can be
observed in the optical. GD 362 and GD 40 (see Figure 1) are two of
the most heavily polluted helium-rich white dwarfs currently
known. Both of these stars also possess infrared emitting debris disks
\citep{kilic05,jura07a,jura07b} which are believed to be the remains
of a single tidally disrupted minor body
\citep{jura06,jura08}. Pioneering analyses of these two stars, based
on Keck high resolution observations, have led to the determination of
the relative abundances of 16 and 9 elements respectively
\citep{zuckerman,klein}. The abundance patterns for these two objects
are remarkably consistent with the idea that the source of the
material is asteroids with a composition similar to that of Bulk Earth
\citep{zuckerman,klein}. Studies of heavily polluted white dwarfs can
thus offer a unique opportunity to learn about the chemistry of rocky
extrasolar planet, asteroid belts, and their subsequent dynamical
evolution following the red giant phases. As such, it is of utmost
importance to increase the sample of white dwarfs with a comprehensive
set of measured element abundances.

\begin{figure}[!ht]
\plotone{f1.ps}
\caption{Medium resolution spectra of
    J0738+1835, GD 40, and GD 362. The two spectra of J0738+1835 and
    that of GD 40 are from SDSS, while the spectrum of GD 362 is taken
    from Gianninas et al. (2004).}
\end{figure}

In what follows, we present a detailed analysis of SDSS
J073842.56+183509.6 (hereafter J0738+1835), the most metal-rich white
dwarf discovered so far. In \S~\ref{observation}, we describe the
observations. Our detailed analysis follows in \S~\ref{analysis}, and
the results are discussed and summarized in \S~\ref{conclusion}.

\section{TARGET SELECTION AND OBSERVATIONS}\label{observation}

Recently, thanks to the Sloan Digital Sky Survey (SDSS), the number of
spectroscopically identified white dwarfs has increased spectacularly
\citep[more than tenfold, and still counting, Kleinman et al. 2010, in
preparation, see also][]{Eisenstein06}, thus providing a large sample
of potential candidates amenable to detailed chemical analysis. In
order to find such candidates, we selected from the SDSS 7th data
release all spectroscopic objects with $g$ $<$ 19.5 mag that fall
within the color space for white dwarf stars. These spectra were then
visually inspected and classified into the various known white dwarf
spectral types. It is while completing this exercise that the
remarkable spectrum of SDSS J073842.56+183509.6 came to our attention
(see Figure 1). J0738+1835 is a DBZ white dwarf that shows not only
the usual Ca II H and K lines, but also several exceptionally strong
lines of O, Fe, Mg and Si. Remarkably, note that the Mg II
$\lambda$4481 line is deeper than He I $\lambda$4471, a unique feature
among known polluted white dwarfs. The severity of the metallic
pollution of J0738+1835 is particularly striking when its spectrum is
compared to those of GD 362 and GD 40, two objects already known to be
among the most polluted helium-rich white dwarfs currently identified
(see Figure 1). Hence, J0738+1835 is a white dwarf star that offers
great potential for a precise measurement of circumstellar material
composition, and it is only the third helium-rich white dwarf (after
GD 362 and GD 40) for which such a detailed analysis is possible from
optical spectroscopy.

The SDSS photometric magnitudes, in the $ugriz$ system, are $u$ =
17.527 $\pm$ 0.014, $g$ = 17.577 $\pm$ 0.009, $r$ = 17.822 $\pm$
0.009, $i$ = 18.077 $\pm$ 0.010 and $z$ = 18.321 $\pm$ 0.025,
suggesting an effective temperature in the 13,000-15,000 K range,
depending on the amount of reddenning that is assumed. Two spectra of
this object, covering the 3800-9200 \AA~ region at a resolution of
$\sim$ 3\AA~FWHM are available from the SDSS archive. The
SDSS-Modified Julian Date, plate and fiber ID number for these two
spectra are respectively 53431-2054-346 and 54495-2890-354. These two
original SDSS spectra of J0738+1835 presented in Figure 1 are too
noisy for a precise determination of the atmospheric parameters
(effective temperature, surface gravity, and abundances of the various
elements). We thus, as a first step, secured a new medium resolution
high S/N spectrum on UT 2009, November 19, using the 6.5 m MMT
telescope on Mount Hopkins, Arizona, equipped with the Blue Channel
Spectrograph. We used a 1$\arcsec$ slit and the 500 line mm$^{-1}$
grating in first order to obtain spectra with a wavelength coverage of
3500 $-$ 6630 \AA\ and a dispersion of 1.2 \AA~ per pixel. All spectra
were obtained at the parallactic angle. We used He-Ne-Ar comparison
lamp exposures and blue spectrophotometric standards \citep{massey88}
for wavelength and flux calibration, respectively.

In order to verify if an emitting debris disk is present (see below),
infrared observations are needed. We obtained $JHK$ photometry of
J0738+1835 using the Near Infra-Red Imager and Spectrometer (NIRI) on
Gemini-North. These observations were performed on 2010, January 03 as
part of the Director's discretionary time program GN-2009B-DD-8.  We
used a 13 position dither pattern with 21-30 s exposures. We do not
detect any sources above the background within 10" of the target. We
used the Gemini/NIRI package in IRAF to reduce the data and the UKIRT
faint standards \citep{leggett06} to calibrate the photometry. The
derived magnitudes, in the Mauna Kea photometric system, are $J$ =
17.965 $\pm$ 0.031, $H$ = 17.758 $\pm$ 0.033 and $K$ = 17.327 $\pm$
0.034 (or 0.1034 $\pm$ 0.0031 mJy, 0.0791 $\pm$ 0.0026 mJy and 0.0752
$\pm$ 0.0025 mJy respectively).

\section{DETAILED ANALYSIS}\label{analysis}

\subsection{Model Atmospheres and Fitting Technique}

The standard technique to determine the atmospheric parameters for
helium-rich white dwarfs is to perform a $\chi^2$ minimization of the
normalized He I line profiles using a grid of synthetic spectra.
However, in the case of J0738+1835, most of the He I lines are
contaminated by heavy element absorption and finding good
normalization points in the continuum is difficult due to the presence
of numerous metallic lines. Moreover, metallic absorptions have a
significant impact on the thermodynamic structure of the models,
making estimations of the effective temperature and gravity from the
He I lines alone unreliable. For example, in Figure 2, we show the
influence of metals on the temperature and pressure structure of a
\logg = 8, \Te = 14,000 K model atmosphere. This effect is explained
mostly by the flux redistribution due to line blanketing in the UV. We
find that the iron lines provide the dominant contribution in that
regard: for example, the thermodynamic structure of a model calculated
without the iron line opacity is much closer to that of a pure helium
model.

\begin{figure}[!ht]
\plotone{f2.ps}
\caption{Temperature and pressure as a function
    of the Rosseland optical depth for white dwarf model atmosphere
    with \Te = 14,000 K and \logg = 8.0. Thick line is for a pure
    helium atmosphere while the thin line is for a helium-rich
    atmosphere that includes traces of heavy elements (log (Ca/He) =
    -7.0 with the abundances of other elements relative to Ca set as
    described in the text).}
\end{figure}

In this connection, it is interesting to note that our own fit of an
optical SDSS spectrum of GD 40 \citep[another DBZ star analysed in
depth by][] {klein}, using only the He I lines as in the standard
approach and assuming \logg = 8.0 and a pure helium atmosphere, gives
\Te = 15,180 K, very similar to the 15,300 K of \citet{voss07}. In
contrast, a fit with a grid that includes heavy elements yields an
effective temperature up to 1000 K cooler, depending on the exact
amount of metals assumed. This should serve as a warning concerning
the use of pure helium atmospheres for the analysis of such stars.

Given those circumstances, we prefer to fit simultaneously the He I
lines and the numerous iron lines. This is performed by fitting the
whole spectrum with the solid angle, effective temperature, surface
gravity, and the abundance of iron as free parameters. In order to
take into account, to first order, the effect of other heavy elements,
we simply assume that the proportions of all these elements relative
to iron are the same as CI chondrites \citep{lodders}. This is a
reasonable assumption, as demonstrated by the study of GD 362 and GD
40 \citep{zuckerman,klein}, and has the advantage of minimizing the
consequences of unaccounted absorption from heavy elements other than
iron. In other words, it is better to include absorption from
additional metals rather than neglect their contribution completely.
This was chosen over a time-consuming strategy involving navigation in
a N(Z)+2 dimensinoal parameter space (N(Z) elements, effective
temperature, and surface gravity).

We use state-of-the-art model atmospheres based on a code similar to
that described in \citet{dufour07} and \citet{dufourphd}. We use
calcium as a reference element and set the abundances of all elements
with the proportions given by \citep{lodders}. This assumption needs
to be verified a posteriori (see below). Our model grid covers a range
from \Te = 13,000 to 16,000 K in steps of 500 K, \logg = 7.5 to 9.0 in
steps of 0.5 dex, and from log(Ca/He) = $-$5.5 to $-$7.0 in steps of
0.5 dex. Additional models with traces of hydrogen have also been
calculated. All models are calculated with the standard ML2
parametrization of the convective efficiency, but with a mixing length
$\alpha =1.25$ \citep{beauchamp}.

We first start our fitting procedure by assuming the canonical value
\logg = 8.0, but other values of the surface gravity ultimately had
to be explored (see below). We evaluated, as explained above, the
effective temperature and the iron abundance by fitting the optical
spectra. If the lines of a given heavy element are found to be in
disagreement with the observations, we readjust its abundance
accordingly and a new model is recalculated with the new abundances in
a self-consistent way. It is found that small variations of the
abundances of elements other than iron have only a marginal impact on
the thermodynamic structure of the star. Our inferred parameters are
thus not sensitive to our initial assumptions.

However, we noticed that it was not possible to get a satisfactory
simultaneous fit to the Mg I and Mg II lines under the assumption of
log $g$ = 8.0. For example, a good fit to the various Mg II lines
could be achieved only at the expense of a very bad fit to the Mg I
doublet ($\lambda\lambda$ 3832, 3838), the latter being not strong
enough for the assumed atmospheric parameters. This could have been an
indication that the effective temperature was overestimated, but the
much lower value needed to reconcile the Mg I and Mg II line strengths
was clearly incompatible with the He I and Fe I/Fe II line
profiles. Our first idea was that the thermodynamic structure of our
model was wrong due to one of the assumptions about the composition,
but models calculated by including only the observed elements were
essentially identical to the ones with all the elements included. We
also explored the parameter space to see if another combination of
abundances and effective temperature could solve that problem, but no
solutions were found.

Another way to favor the formation of Mg I relative to Mg II is to
increase the atmospheric pressure. This can be achieved if the surface
gravity is allowed to increase in our fitting procedure. In a manner
similar to that described above, we find that the Mg I and Mg II lines
can be fitted simultaneously (as well as the He I, Fe I, and Fe II
lines) if we increase the surface gravity to \logg = 8.5. Slightly
larger/lower values of \logg~ (compensated with an accordingly
larger/lower effective temperature) can also accommodate
simultaneously the Mg I and Mg II lines, but then the He I lines fit
is not as good. Our derived mass determination thus relies heavily on
the magnesium lines and the somewhat uncertain He I line profiles
(usually attributed to the treatment of van der Waals broadening).

Trends of higher mass at low effective temperature for DB stars has
raised concerns in the literature about the validity of the derived
masses from He I lines. However, as clearly shown in \citet[][see
their Figure 7]{limoges}, it seems that the spread in mass is probably
real after all. Since a higher than average surface gravity best fits
both the He I lines and the Mg lines, we have a high confidence that
J0738+1835 is indeed relatively massive. The debate will be settled
when a trigonometric parallax measurement and better broadening
theories become available, but until that happens, we consider our
solution as the best compromise.

There appears to be large discrepancies in the predicted line strength
of a few iron lines, the most noticeable being that of Fe II
$\lambda$4173.46 and near 5200 \AA (see Figure 3 and 4). These
discrepancies cannot be eliminated even by varying significantly from
our optimal solution the surface gravity and the effective temperature
of our model. It is most probable that this is a manifestation of the
uncertainties in the atomic data, $\log$ gf and Stark broadening
parameters in particular, used in our calculation (we used the well
known Kurucz linelists). We calculated a posteriori a synthetic
spectrum using the Vienna Atomic Line Database (VALD) lists. Using
these data helps to partly reduce the discrepancies for a few iron
lines (for example, the depth of the Fe II $\lambda$4173.46 is reduced
by half) but does not alter the abundance determinations for any
elements on average. Since our fit is based on numerous iron lines,
these few discrepant lines have a relatively weak weight on the iron
abundance determination and we present our final solution using the
Kurucz linelist for consistency with the rest of our analysis.

\begin{figure}[!ht]
\plotone{f3.ps}
\caption{MMT spectrum (in black) and our best
    fit model (red line). The strongest lines are indicated by tick
    marks. All the other numerous non-identified lines are from
    iron. They have not been marked for clarity. In order to take into
    account uncertainties on the slope of the continuum due to
    extinction and flux calibration, a small linear and quadratic term
    were also included in the fitting procedure.}
\end{figure}

\begin{figure}[!ht]
\plotone{f4.ps}
\caption{Same as Figure 3, but for the red
    part of the spectrum. Note that due to second order contamination
    longward of 5900 \AA, a second-order polynomial fit has been
    applied to the model spectrum continuum.}
\end{figure}

We note that the assumption on the abundance of heavy elements
described above also predicts weak lines of sodium (the Na D lines)
that are blended with the He I $\lambda$5876 line. There appear to be a
feature in the red wing of the He I line that is consistent with the
expected lines strength but given the signal-to-noise ratio of our
observations and the weakness of the presumable lines, we shy away
from claiming a detection. Two weak lines of Cr and Ti are also barely
predicted at a level comparable to the signal-to-noise ratio of our
observations. Future observations with higher signal-to-noise and
higher resolution (as well as in other part of the electromagnetic
spectrum) should eventually confirm their presence.

Our final solution is presented in in Table~1 while our best fit model
spectra are shown in Figures 3 and 4. Uncertainties obtained from the
covariance matrix are unrealistically small. These represent only
internal errors and thus underestimate the real uncertainties. We thus
evaluate the uncertainties on the model parameters by calculating a
series of models with various parameters near our optimal
solution. The constraint from fitting simultaneously the Mg I, Mg II
and He I lines dominates the errors on the effective temperature and
surface gravity, which we evaluate to be about 300 K and 0.2 dex
respectively. The uncertainties on the element abundances must also
take into account these uncertainties. In a similar manner to that
described in \citet{klein}, we varied the model \Te~ and \logg~ by
their uncertainties, one at a time, to evaluate the corresponding
errors on the various element abundances. The final uncertainties on
the abundances presented in Table~1 are the sum of the errors added in
quadrature. Those uncertainties are then properly propagated for the
determination of the errors on the cooling age, luminosity, mass,
radius and distance of the star (obtained by determining the solid
angle from the $ugriz$ photometry). The evolutionary models used are
similar to those described in \citep{FBB01} but with C/O cores, and
thickness of the helium and hydrogen layers of respectively q(He) =
10$^{-2}$ and q(H) = 10$^{-10}$, which are representative of
helium-rich atmosphere white dwarfs.

\input{tab1}

Finally, we calculated a posteriori a model with the final abundances
of only the observed elements. We find that the resulting synthetic
spectrum is practically identical to the one with all elements
included. This is not surprising given that Fe, Mg, Si, O, and Ca account
for more than $95\%$ of the contaminants by mass for Bulk Earth
\citep[a similar proportion was also found for GD 40,][]{klein}. A
similar proportion will probably be retrieved for J0738+1835 when high
resolution observations become available (ultraviolet observations
should also be helpful).

\subsection{Convection Zone Models}

In order to evaluate the total amount of heavy material mixed in the
outer layers of J0738+1835, we need to know the mass of the helium
convection zone into which the metals are diluted. Estimates of this
mass can readily be obtained using the Montr\'eal white dwarf building
codes that were developed over the years \citep[see,
e.g.,][]{BF94,FBB01}, and which have been maintained at the
state-of-the-art level. One possible option offered by these codes is
to compute complete, but static, stellar structures with a luminosity
profile closely following the mass profile, as appropriate for cool
white dwarfs shining through the loss of thermal energy. This is what
we adopted here.

To ease the comparison with the spectroscopic observations, we further
selected the option of fixing the surface gravity and the effective
temperature (as well as the compositional stratification and the core
composition) in these stellar model calculations. We also assumed
relatively thick helium envelopes in these models, log $(1 - M_{\rm
  He}/M_*) = -3.0$, much thicker than the superficial outer convection
zone due to helium partial ionization. Since the presence of heavy
elements does affect somewhat the extent of the convection zone,
models with an envelope composition consisting of helium and
homogeneous traces of metals were considered along with more standard
models having pure helium envelopes. We retained the results for a
metallicity of Z = 0.001 in what follows.

We investigated also how $M_{\rm He}/M_*$ depends on the assumed
convective efficiency. This is because, as was first shown by
\citet{ber95}, the calibration of the mixing-length derived in the
atmospheric layers (by comparing optical and UV observations) does not
apply at the base of the convection zone, where the flux of settling
heavy elements determines the observable abundances. In that case, the
calibration is best obtained by comparing the effective temperatures
at the blue edge of an instability strip as inferred from nonadiabatic
pulsation calculations with the values derived from spectroscopy. In
the case of pulsating H-rich (DA) white dwarfs, \citet{ber95} found
that the convective efficiency increases with depth, from the
atmospheric layers to the bottom of the convection zone, so that a
variable mixing length is required to account for both the
spectroscopic and pulsational measurements. Unfortunately, the
calibration based on pulsational properties is still uncertain,
particularly for the He-atmosphere white dwarfs \citep [see,
e.g.,][]{FB08}, so we decided to consider three different plausible
versions of the mixing-length theory for constraining the mass of the
outer convection zone in J0738+1835.

The results of our calculations are presented in Table~2 for values of
log $g$ and $T_{\rm eff}$ in the vicinity of those found
spectroscopically for J0738+1835. Our calculations indicate that the
mass of the helium convection zone in that regime is, fortunately,
only weakly dependent on the effective temperature and on the
convective theory flavor used in the stellar models. However, a strong
dependency on the surface gravity is found, $M_{\rm He}/M_*$
decreasing significantly as the surface gravity increases. Using our
atmospheric parameters, we estimate that the mass of the helium
convection zone lies somewhere between $10^{-5.70}M_*$ and
$10^{-6.75}M_*$. Correspondingly, the lower limits on the mass of the
polluting body (i.e. the sum of the masses of the detected elements in
Table~1) are between $2.7 \times 10^{23}$ g and $2.4 \times 10^{24}$
g, with a value of $ 4.3\times 10^{23}$ g corresponding to \logg =
8.5.

\input{tab2}

This is close to the mass of Ceres ($9.4 \times 10^{23}$ g) and is
about an order of magnitude more than the amount of metals found in
the convection zones of GD 40 \citep{klein} and GD 362
\citep{koester09} ($3.6 \times 10^{22}$ g and $1.8 \times 10^{22}$ g,
respectively), and comparable to the amount found in HS 2253+8023
\citep{jura06} which is, however, based solely on the derived Fe
abundance from an $IUE$ spectrum \citep{friedrich}.  Because an
unknown period of time has elapsed since the accretion process has
started, we do not know what fraction of the material has already sunk
out of sight below the base of the convection zone and, therefore,
this estimate for the total mass of heavy elements represents a lower
limit on the mass of the body responsible for this large
pollution. Given this, and the fact that another yet unknown quantity
of material is still orbiting the star in the form of a debris disk
(see below), it is safe to presume that the object that was tidally
destroyed by J0738+1835 was at least as large as the dwarf planet
Ceres.

This lower limit is sensitive to the determined value of the surface
gravity.  It will thus be of utmost importance to determine the
surface gravity with greater accuracy than is currently possible in
order to precisely constrain the mass of the asteroid/dwarf planet
that causes the pollution in J0738+1835. A good trigonometric parallax
measurement should help to reduce this uncertainty significantly.

\subsection{Infrared Photometry and Disk Model}

The next step was to determine if an accretion disk is still present
around J0738+1835, which would be revealed by an infrared excess. We
thus obtained infrared $JHK$ photometry with the 8 m Gemini North
telescope in Hawaii. Our observations, presented in Figure 5, clearly
reveal a near-IR excess.

\begin{figure}[!ht]
\plotone{f5.ps}
\caption{Photometric measurements in the GALEX,
    $ugriz$ and $JHK$ bands compared to models of the star and the
    debris disk.}
\end{figure}

Given the fact that most heavily polluted white dwarfs have debris
disks, our target is likely to have a disk as well. If the
near-infrared excess is due to a companion star, this star would have
$M_{\rm K}$ = 12.35 mag; an L6 dwarf \citep{leggett02}. Such a
companion would have $H-K$ = 0.6-0.9 mag. The near-infrared excess
around our target (after subtracting the contribution from the stellar
photosphere) has a $H-K$ color of 1.33 mag, incompatible with an L6
dwarf or any other M, L, and T dwarfs studied by
\citet{leggett02}. Therefore, the infrared excess around our target is
best explained by the presence of a dusty debris disk rather than a
brown dwarf companion.

We fit the infrared excess around J0738+1835 using the optically thick
flat-disk models of \citet{jura03}. Using the effective temperature,
radius, and distance estimates for the white dwarf (see Table~1), we
created a grid of disk models with inner temperatures 1000-1800 K,
outer temperatures 300-900 K (in steps of 100 K), and inclination
angles 15$^\circ$-75$^\circ$ (in steps of 15$^\circ$). The model that
best fits the infrared photometry has $T_{\rm in}=$ 1700 K, $T_{\rm
  out}=$ 700 K, and an inclination angle of 60$^\circ$. However, since
we only have near-IR data, the outer temperature is not well
constrained, and these parameters change depending on the assumed
inclination. For example, if the inclination angle is assumed to be
75$^\circ$, then the best fit model has $T_{\rm in}$=1600 K and
$T_{\rm out}$=300 K. Nevertheless, the near-IR data can be explained
with a reasonable set of parameters for a dusty disk, indicating that
the accretion process is most probably occurring at present. An inner
temperature of 1700 K is somewhat unusual for disks around white
dwarfs. However, the near-infrared flux excess around J0738+1835 is
similar to the excess seen around the white dwarf GD 56
\citep{kilic06}. A flat-disk model with an inner temperature of 1700 K
is required to explain the observed excess around GD 56
\citep{jura07a}. A better fit to the near-infrared data for
GD 56 can be obtained if the disk is substantially warped or puffed up
by the gravitational field of a planet \citep{jura09}. Similar
processes may be responsible for the observed near-infrared flux
distribution of the disk around J0738+1835. Mid-infrared observations
will be useful to check the flat-disk models for this target.

Unfortunately, the total mass of the disk cannot be determined for an
opaque ring. Masses for the dust disks around G29-38 and GD 362 have
been estimated from models of the mid-infrared emission
\citep{reach,becklin}. Based on the choice of optically thin or thick
models, these disks hold $\sim 10^{19}-10^{24}$ g of material. The
disk around J0738+1835 is likely to hold a similar amount of material
as well. Hence, the combined mass of the metals in the surface
convection zone and the dust disk is on the order of 4-14 $\times
10^{23}$ g, equivalent to the mass of a dwarf planet.

\section{DISCUSSION AND CONCLUSION}\label{conclusion}

According to the accretion/diffusion scenario
\citep{dupuisa,dupuisb,dupuisc}, a steady state abundance for a given
heavy element is rapidly reached after the accretion process
starts. Under those circumstances, the ratio of abundances measured
spectroscopically in the photosphere of J0738+1835 can be related,
unlike the case of HS 2253+802 which shows no disk \citep{farihi09},
to that of the incoming disk material.

It is interesting to note that the mass ratios of the most abundant
heavy elements found in J0738+1835, namely O, Mg, Si and Fe, are very
similar to that found in GD 40, which has an atmosphere polluted with
debris coming from a body with an Earth-like composition
\citep{klein}. However, calcium is depleted in J0738+1835 by a factor
of 4, perhaps indicating that mantle/crust differentiation has
occurred before the crust was lost. Moreover, if we assume that the
total amount of hydrogen present in the atmosphere of J0738+1835
originates from the last polluting event only, we can place an upper
limit of $\sim$1$\%$ by mass on the water or ice content of the
polluting body. This is a surprisingly low amount considering that
most asteroids in our solar system contain large quantities of ice on
their surface. This might indicate that the parent body was
sufficiently close to the progenitor red giant star for thermal
processes to have evaporated most of the ice or water initially
present \citep{jura10}. Thus it is possible that the surface of the
orbiting body was significantly altered in the late stages of stellar
evolution, but at this point it is premature to quantify this
effect. Given the uncertainties on the atmospheric parameters, we
shall refrain from speculating further on the true fate of the body
that polluted the atmosphere of J0738+1835. However, when parallax and
higher resolution spectroscopy observations become available,
J0738+1835 will become a perfect testbed for a precise study of not
only the abundances of the polluting extrasolar body, but also of the
effects of processes such as crust differentiation and thermal
heating.

It is worth noting that medium resolution observations of GD 362 and
GD 40 revealed only relatively weak absorptions of Ca, Mg, and Fe in
the optical (Figure 1). It was only when Keck high
resolution spectra became available that a myriad of elements were
uncovered, providing extremely valuable information about the accreted
material.  Given that J0738+1835's medium resolution spectrum is
showing even stronger metallic absorption lines than in the cases of
GD 362 and GD 40, further high resolution observations in the optical
and ultraviolet to study elemental abundances as well as at longer
wavelengths to characterize the debris disk are highly
desirable. These observations will undoubtedly offer an opportunity to
study the composition of an extrasolar dwarf planet with unprecedented
accuracy.

\acknowledgements This work was supported in part by NSERC Canada and
FQRNT Qu\'ebec. P.D is a CRAQ postdoctoral fellow, and P.B is a
Cottrell Scholar of Research for Science
Advancement. G.F. acknowledges the contribution of the Canada Research
Chair Program. MK acknowledges support from NASA through the Spitzer
Space Telescope Fellowship Program, under an award from Caltech. Based
on observations obtained at the MMT and Gemini Observatory. The MMT is
a joint facility of the Smithsonian Institution and the University of
Arizona. The Gemini observatory is operated by the Association of
Universities for Research in Astronomy, Inc., under a cooperative
agreement with the NSF on behalf of the Gemini partnership: the
National Science Foundation (United States), the Science and
Technology Facilities Council (United Kingdom), the National Research
Council (Canada), CONICYT (Chile), the Australian Research Council
(Australia), Ministerio da Ciencia, Tecnologia (Brazil) and Ministerio
de Ciencia, Tecnologia e Innovacion Productiva (Argentina).

Solar and Stellar Astrophysics (astro-ph.SR)


Polar-bulge galaxies

Abstract: Based on SDSS data, we have selected a sample of nine edge-on spiral galaxies with bulges whose major axes show a high inclination to the disk plane. Such objects are called polar-bulge galaxies. They are similar in their morphology to polar-ring galaxies, but the central objects in them have small size and low luminosity. We have performed a photometric analysis of the galaxies in the g and r bands and determined the main characteristics of their bulges and disks. We show that the disks of such galaxies are typical for the disks of spiral galaxies of late morphological types. The integrated characteristics of their bulges are similar to the parameters of normal bulges. The stellar disks of polar-bulge galaxies often show large-scale warps, which can be explained by their interaction with neighboring galaxies or external accretion from outside.

\section{Introduction}

The bulge is one of the main subsystems of spiral
galaxies. Bulges distinguish from the other galactic subsystems
(disks, bars, spiral arms, etc.) by their brightness
distribution, shape, stellar population,
the pattern of stellar motion, and other characteristics.
For a long time the bulges were considered
as a kind of small elliptical galaxies surrounded
by stellar disks. In the last couple of decades this
simplified picture has changed to a more complex one:
it turned out that the bulges are not a homogeneous
type of objects. They break up into at least two different
groups: the classical bulges that resemble elliptical
galaxies in many respects (see, e.g., Renzini, 1999)
and the pseudo-bulges that are closer to galactic disks
in a number of characteristics (Kormendy and Kennicutt
2004; Kormendy 2015). Athanassoula (2005)
introduced the third type,
boxy/peanut bulges, that are actually edge-on bars. The
situation is further complicated by the fact that different
types of bulges can coexist in the same galaxy (see,
e.g., Erwin et al. 2015).

In this paper we consider yet another type of
bulges, the bulges whose apparent major axis is highly
inclined to the galactic midplane. The
objects with such bulges are similar in their morphology
to the polar-ring galaxies whose polar component
dominates in luminosity. Corsini et al. (2012)
called such bulges {\it polar}. The polar bulges are extremely
rare. In a recent review devoted to the shape
of galactic bulges, M\'endez-Abreu (2015) mentions
only three such galaxies: NGC 4698 (Bertola et al.
1999; Corsini et al. 2012), NGC 4672 (Sarzi et al.
2000), and UGC 10043 (Matthews and de Grijs
2004). One of the reasons for the small number
of known polar bulges is that they can be identified
only in the case where the galactic disk is close to the
edge-on orientation to the line of sight and therefore,
the bulge inclination to the disk plane can be easily
detected. If this is not the case, then the bulge
looks like a bar and does not attract the attention of
researchers.

Given the small number of known objects,
the systematics of properties and the origin of the polar
bulges remain puzzles. In two of the three galaxies
(NGC 4698 and NGC 4672) the morphologically
decoupled bulges are known to be related to
the kinematically decoupled subsystems in their cores
that  rotate almost orthogonally to the galactic disks
(Bertola et al. 1999; Sarzi et al. 2000). No kinematically
decoupled core was detected in UGC 10043
(Matthews and de Grijs 2004), although as the authors
note, this could be hindered by the limited spatial
resolution of the spectra they used and by the
powerful absorption stripe that shields the core.

The existence of kinematically and morphologically
decoupled structures in galaxies is usually associated
with a ``secondary'' event in their history. To
explain the formation of polar bulges, Matthews and
de Grijs (2004) considered three possible ``secondary''
events: the capture of a disk by a pre-existing ``naked''
bulge (this scenario resembles the formation of polar-ring
galaxies with the external accretion; Reshetnikov
and Sotnikova 1997); the capture of a
small elliptical galaxy by a spiral galaxy (this mechanism
resembles the formation of collisional rings, but
the relative velocity of the galaxies must be small;
Appleton and Struck-Marcell 1996); and merging
of two disk galaxies (Bekki 1998). Qualitatively, all
these scenarios are capable of describing a number
of peculiarities of polar-bulge galaxies, and extensive
observational data are needed to choose between
them. First of all, it is necessary to
increase the sample of known objects of this type
and to perform their detailed observational study and
simulations.

In this paper, we present a small sample
of new candidates to the polar-bulge galaxies and some
results of their photometric analysis.
All numerical values in the paper
are given for the cosmological model with the Hubble
constant of 70 km s$^{-1}$ Mpc$^{-1}$ and
$\Omega_m=0.3$, $\Omega_{\Lambda}=0.7$.

\section{The sample of polar-bulge galaxies and their analysis}

\subsection{The sample of galaxies}

In order to search for galaxies with morphologically decoupled
bulges we examined images of galaxies
from the EGIS (Edge-on Galaxies In SDSS) catalog
(Bizyaev et al. 2014), which includes almost
six thousand edge-on galaxies. Our examination
revealed more than twenty objects whose bulges visually
appeared elongated along the galaxy minor
axis. As the next step, we performed a photometric decomposition
of the images of all selected galaxies and
estimated parameters of the
disks and bulges, including their apparent flattening
and position angles of the apparent major axes (see
the next section for the details of our analysis).
Our final sample of nine candidates to the polar-bulge galaxies
includes only those objects whose bulges in the $g$
and $r$ bands showed an apparent flattening $b/a < 0.9$,
while the difference between the position angles of the
major axes for the disks and bulges exceeded 30$^{\rm o}$.
The accuracy of the model parameters (see below) does
not allow us to reliably judge whether the bulges with
$b/a \geq 0.9$ can belong to the polar ones.

Table 1 is based of the SDSS\footnote{http://www.sdss.org}
and NED\footnote{http://ned.ipac.caltech.edu} data, and summarizes
the main characteristics of the final sample
of the polar-bulge galaxies. Their $r$-band contour maps are shown in
Fig. 1. As it can be seen from the table, two galaxies
(\#1 and \#5) were previously included in the catalog
of candidates to the polar-ring galaxies (SPRC; Moiseev
et al. 2011) due to their unusual morphology. The
galaxy UGC 10043 (no. 7 in Table 1) was studied by
Matthews and de Grijs (2004) (see the Introduction).

\begin{table*}
\caption{Candidates for the polar-bulge galaxies}
\begin{center}
\begin{tabular}{|c|c|c|c|c|c|}
\hline
         &           &              &      &            &     \\
\#        & SDSS name & Alternative name & Type & $r$ (mag)  &  $z$  \\
         &           &              &      &            &     \\
\hline
         &           &              &      &            &     \\
         1  &  SDSS J015858.39-002923.2     &  SPRC-77 &     & 17.53 & 0.08099 \\
         2  &  ~~~~~~~~~~J095346.69+351456.8&  &     & 16.62 & 0.03928 \\
         3  &  ~~~~~~~~~~J102343.67+421917.9&  &     & 15.59 & 0.04611 \\
         4  &  ~~~~~~~~~J124331.21-142015.8 &  FGC~1494, RFGC~2365 &  Sc & 16.72 & \\
         5  &  ~~~~~~~~~~J133904.58+020949.5&  UGC~8634, FGC~1649, &  Sb & 15.10 & 0.02337 \\
            &                               &  RFGC~2619, SPRC-42 & & & \\
         6  &  ~~~~~~~~~~J151223.37+013823.9&  &     & 16.58 & 0.02850 \\
         7  &  ~~~~~~~~~~J154841.10+215208.7&  UGC~10043, FGC~1953, & Sbc & 14.34 & 0.00721 \\
            &                               &  RFGC~3043 & & & \\
         8  &  ~~~~~~~~~~J160028.24+143201.2&  &     & 15.96 & 0.03413 \\
         9  &  ~~~~~~~~~~J222006.54+253621.0&  &     & 16.96 & 0.04164 \\
            &        &              &       &           &  \\
\hline
\end{tabular}
\end{center}
\end{table*}

\begin{figure}
\centering
\includegraphics[width=7.4cm, angle=0, clip=]{gal1.eps}\\
\includegraphics[width=7.4cm, angle=0, clip=]{gal2.eps}\\
\includegraphics[width=7.4cm, angle=0, clip=]{gal3.eps}\\
\includegraphics[width=7.4cm, angle=0, clip=]{gal4.eps}\\
\includegraphics[width=7.4cm, angle=0, clip=]{gal5.eps}\\
\includegraphics[width=7.4cm, angle=0, clip=]{gal6.eps}\\
\includegraphics[width=7.4cm, angle=0, clip=]{gal7.eps}\\
\includegraphics[width=7.4cm, angle=0, clip=]{gal8.eps}\\
\includegraphics[width=7.4cm, angle=0, clip=]{gal9.eps}\\
\caption{The $r$-band contour maps of the galaxies with morphologically
decoupled bulges. The contours were selected in a way to
highlight the global structure of each galaxy. The galaxy number from
Table 1 is indicated in the upper left corner of each panel.
The length of the bar in the lower left corner is 10 arcsec. }
\end{figure}

\subsection{Photometric analysis}

To determine the photometric and structural parameters
of the considered galaxies, we decomposed
their $g$- and $r$-band images within the framework
of the two-component ``bulge + disk'' model.
The galaxy images were taken from the SDSS archive (SDSS
DR12; see Alam et al. 2015).

To describe the surface brightness distribution in
the bulge, we used the Sersic (1963) law:

\begin{equation}
I(r) =
I_{0,\mathrm{b}} \exp
  \left[-\nu_n\left(\frac{r}{r_e}\right)^{\frac{1}{n}}\right],

\label{eqSersicMu}
\end{equation}
where $I_{0,\mathrm{b}}$, $r_\mathrm{e}$ and $n$ are the central
surface brightness, effective radius, and Sersic parameter, respectively;
$\nu_n$ is an $n$-dependent constant chosen in such
a way that a half of the total bulge luminosity was
encompassed within the $r_e$. The standard model of an edge-on
disk (van der Kruit and Searle 1981) was used as
the disk model:

\begin{equation}
I(r, z) =
I_{0,\mathrm{d}}
  \left(\frac{r}{h}\right) K_1
    \left(\frac{r}{h}\right) \sech^2\left(\frac{z}{z_0}\right),
\label{eqDiscEOn}
    \end{equation}
where $I_{0,\mathrm{d}}$, $h$ and $z_0$
are the central surface brightness,
radial scale length and vertical scale height of the disk, and
$K_1$ is a modified first-order Bessel function.

The decomposition, i.e., the search for optimal parameters
that ensured the minimal residual between
the model and observations, was performed according
to the following scheme. At the preliminary stage,
the preparation of images was made, including the
image rotation in such a way that the galaxy major
axis was oriented along the x-axis of the image (the
method from Martin-Navarro et al. (2012) was used
to find the galaxy position angle), the cropping of
the images, and the masking of background objects and
dust lanes in the galaxies (if they were visible). At
the next stage we obtained the first approximation
of the galaxy model using the Differential Evolution
algorithm from the "Imfit" package (Erwin 2015). The
advantage of this algorithm is that it does not require
specific initial conditions for the optimization process
and only the range of possible parameters is needed
for its operation. The parameters obtained in this
way were then improved using the gradient descent
algorithm implemented in the "galfit" package (Peng
et al. 2010). The output parameters from the "galfit" were
assumed as the final galaxy model parameters.

Since the galaxies in our sample have rather
small angular sizes, atmospheric seeing
strongly affects the parameters being determined
(e.g., Trujillo et al. 2001). To compensate for this
effect, we convolved the photometric model with the
point spread function (PSF) simulating atmospheric
seeing and compared the convolved model with the
observed image.
Here, we used the Moffat function
(Moffat 1969) as the PSF. We determined the parameters
of the Moffat function by fitting isolated field
stars. As a rule, there were 10 - 20 such PSF stars in
the field of each galaxy, which allowed us to reliably
determine the PSF parameters and to estimate the
dispersion of these parameters.

To estimate the errors of the decomposition results,
we performed Monte Carlo simulations by assuming
that the main sources of errors were the image
noise and the discrepancy between the model and
real PSFs. For each galaxy, we constructed a model
corresponding to the parameters determined during
the decomposition, convolved it with the observed
PSF, and added noise whose value corresponded to
the real parameters.
The PSF parameters differed from their
true values according to a normal law with a dispersion
equal to the dispersion of the parameters of individual
PSF stars obtained at the PSF construction
stage. For each galaxy, we repeated the decomposition procedure
250 times and obtained the mean values and the
uncertainties of the decomposition
parameters.

The mean difference between the total magnitudes
(defined as the sum of the model bulge and disk)
and the {\em modelMag} magnitudes from the SDSS is
$-0\fm37 \pm 0\fm29$  in the $g$ band and
$-0\fm23 \pm 0\fm25$ in
the $r$ band. If we take into account the fact that
the SDSS photometric data are marked as unreliable
for four galaxies, then the agreement becomes better
for the remaining five galaxies:
$-0\fm21 \pm 0\fm10$ and $-0\fm10 \pm 0\fm14$
in the same bands. In the subsequent
discussion, we use the results of our photometry.

The results of our photometric decomposition of the galaxies
from the sample in the $r$ band are shown in Table~2.
For each object, the first row of the table gives the
measured parameter and the second row gives the
measurement undertainty of this parameter estimated as
described above. The first column in Table 2 contains
the galaxy number according to Table 1. The
next four columns present the bulge parameters: the
effective surface brightness, effective radius, Sersic
index and apparent flattening. The next columns show the following characteristics of
the galactic disks: the observed central
surface brightness, exponential scale length, and vertical
scale height of the brightness distribution. The
observed difference between the position angles of the
major axes for the bulge and disk is given in the ninth
column of Table 2. (Note that $\Delta$P.A. in the table is
only the model-averaged difference of the position angles.
In practice, the position angles of the major axes
for the disks and bulges exhibit noticeable variations
for a number of galaxies.) The last three columns
give the galaxy absolute magnitude, its color index,
and the bulge-to-total
luminosity ratio. The surface brightnesses, absolute magnitudes,
and color indices in Table 2 were corrected
for the foreground extinction in the Milky Way according to Schlafly and
Finkbeiner (2011).

\begin{table*}
\caption{Photometric characteristics of the galaxies in the $r$ band}
\begin{center}
\begin{tabular}{|c|c|c|c|c|c|c|c|c|c|c|c|}
\hline
N & $\mu_e$ &  $r_e$ & $n$ & $q$  & $\mu_0$ & $h$   & $z_0$ & $\Delta$P.A. & $M_r$ & $g - r$ & $B/T$ \\
  &         & ($''$) &     &      &         & ($''$)& ($''$)& ($^{\rm o}$) &       &       \\
  \hline
  1 & 22.57   & 2.17    & 6.84& 0.85 &  22.16  & 4.25  & 0.88  &    73.0      & -20.51 & 0.75 &0.80  \\
    &  0.33   & 0.32    & 2.40& 0.04 &   0.10  & 0.17  & 0.06  &     3.6      &       &   &    \\
\hline
2 & 22.03   & 2.73    & 1.82& 0.87 &  20.93  & 5.09  & 1.21  &    74.7      & -19.85& 1.02 &0.36 \\
      &  0.44   & 0.49    & 2.23& 0.10 &   0.23  & 0.32  & 0.15  &    10.3      &       &  &    \\
  \hline
  3 & 21.46   & 2.17    & 5.44& 0.69 &  21.22  & 6.37  & 2.77  &    71.3      & -20.86& 0.56  &0.36 \\
    &  0.17   & 0.15    & 0.86& 0.03 &   0.02  & 0.07  & 0.03  &     1.8      &       &  &    \\
\hline
4 & 21.94   & 2.55    & 0.39& 0.47 &  21.31  & 6.48  & 1.60  &    89.5      &       & 0.58 &0.18 \\
  &  0.03   & 0.03    & 0.02& 0.01 &   0.01  & 0.06  & 0.02  &     0.6      &       &  &    \\
  \hline
  5 & 21.38   & 4.12    & 2.53& 0.42 &  21.39  & 12.00 & 3.10  &    72.2      & -20.07& 0.57 &0.36 \\
& 0.04    & 0.06    & 0.16& 0.02 &   0.01  &  0.09 & 0.02  &     0.4      &       &  &    \\
\hline
6 & 20.19   & 1.41    & 1.50& 0.66 &  22.47  & 8.22  & 2.43  &    49.8      & -19.09& 0.44 &0.56 \\
  & 0.13    & 0.07    & 0.08& 0.02 &   0.06  & 0.37  & 0.29  &     3.2      &       &  &    \\
  \hline
  7 & 20.89   & 5.95    & 0.86& 0.80 &  21.61  & 23.67 & 5.43  &    89.7      & -18.75& 0.64 &0.44 \\
    &  0.01   & 0.04    & 0.02& 0.01 &   0.01  &  0.08 & 0.02  &     0.1      &       &  &    \\
\hline
8 & 21.40   & 3.04    & 4.13& 0.42 &  21.88  & 8.28  & 2.97  &    39.1      & -20.09& 0.65 &0.47 \\
  &  0.06   & 0.08    & 0.48& 0.03 &   0.03  & 0.16  & 0.02  &     0.7      &       &  &    \\
  \hline
 9 & 22.07   & 2.69    & 0.72& 0.68 &  21.37  & 7.54  & 2.25  &    87.4      & -20.12& 0.57 &0.20 \\
&  0.36   & 0.23    & 0.08& 0.03 &   0.08  & 0.21  & 0.12  &     6.7      &       & &     \\
\hline
\end{tabular}
\end{center}
\end{table*}

\section{Results and discussion}

\subsection{General characteristics of the galaxies}

The integrated characteristics of the polar-bulge galaxies
in the $r$ band are as follows: the mean
absolute magnitude $< M > = -19.9 \pm 0.7$, the mean
color index $< g - r > = 0.64 \pm 0.16$, and the exponential
scale length of the disk $< h > = 5.1 \pm 1.0$ kpc. In
all these characteristics, our objects are typical bright
spiral galaxies of late morphological types (see, e.g.,
Bizyaev et al. 2014). The mean vertical scale height of
the stellar disk in the $r$ band, $< z_0 > = 1.50 \pm 0.54$ kpc,
is also near the peak of the distribution for the galaxies
from the EGIS catalog. If we exclude the two galaxies
(\#3 and \#6 in Table 1) whose inclination visually
deviates most strongly from 90$^\mathrm{o}$, then this value decreases
to $< z_0 > = 1.36 \pm 0.45$ kpc.

The mean relative disk thickness for the galaxies of
our sample is $< z_0 / h > = 0.285 \pm 0.07$; without galaxies
\#3 and \#6, it is $< z_0 / h > = 0.265 \pm 0.06$. These values
are close to the galactic disk thickness in well-defined
samples of edge-on spiral galaxies (Mosenkov et al. 2015).

The mean ratio of the disk exponential scale length in
the $g$ and $r$ bands for the considered sample is 1.09$\pm$0.10,
which implyes the existence
of a color gradient: the galactic disks become bluer
towards the periphery. This feature is typical for
most of spiral galaxies.

An interesting structural peculiarity of the galaxies
in our sample is that their disks exhibit prominent
integral-shaped warps (see Fig. 1) in approximately
half of the cases. Since the warp is easier
to detect in the case where its line of nodes is close
to the line of sight, it can be concluded that the
stellar disks in most of the polar-bulge galaxies are warped.
Strong optical warps are observed in galaxies relatively
rarely. As a rule, their existence is associated
with the gravitational interaction of galaxies and
external accretion (Reshetnikov and Combes 1998,
1999; Ann and Park 2006).

We did not specially study the spatial environment
of the galaxies from our sample, but a simple analysis
of the SDSS images (in particular, deeper images
from Stripe 82 are available for galaxy 1 from our
list; Abazajian et al. 2009) showed that relatively
close galaxies of comparable luminosity and a number
of fainter companions are observed near most of the
objects (six out of nine). Unfortunately, there are
no available redshifts for the neighboring galaxies,
which does not allow us to decide about their spatial
association with our objects.

More reliable data are available for previously
known galaxies with decoupled bulges. UGC 10043
(\#7 in Table 1) is a member of a group, and it
interacts with the galaxy MCG+04-37-035
(Aguirre et al. 2009). Two other galaxies
(NGC 4672 and NGC 4698; see Introduction)
are not isolated either. NGC 4672 is a member
of a group of galaxies (Garcia 1993); NGC 4698 is a
member of the Virgo cluster of galaxies.

\subsection{Photometric characteristics of bulges}

The distributions of the main characteristics of
the decoupled bulges in the $r$ band are shown in Fig. 2.
By comparing these distributions with the characteristics
of the central objects in polar-ring galaxies (Fig. 2
in Reshetnikov and Combes 2015) it can be noted
that, as expected, the bulges are more
compact and faint, on average. For example, the mean absolute
$r$-band magnitude for the bulges of our sample is
$< M > = -18.95 \pm 0.77$ (compare with -20.34 for the centers
of the polar-ring galaxies), while the mean effective radius
in the same band is $< r_e > = 1.86 \pm 0.74$ kpc (3.2 kpc
for polar-ring galaxies). The distribution of the Sersic
index for the bulges appears almost flat (Fig. 2b) and
does not show a distinct peak typical for the bulges
of spiral galaxies of early morphological types (Tasca
and White 2011). In general, the observed distribution
is close to that for the bulges of the late-type
spiral galaxies, with caution from the limited size of our sample.

\begin{figure}
\centering
\includegraphics[width=16cm, angle=-90, clip=]{fig2.ps}\\
\caption{Observed distributions of the polar-bulge characteristics in
the $r$ band: the absolute magnitudes (a), Sersic index (b), and
effective radius (c).}
\end{figure}

The mean color index of the decoupled bulges,
$< g - r > = 0.88 \pm 0.13$, is considerably redder than
that of the surrounding disks, for which $< g - r > =
0.49 \pm 0.22$. Given that the colors of the edge-on
disks were not corrected for the internal extinction, the
color difference between the bulges and disks must
be even greater. This large color difference between
the bulges and disks is typical for galaxies of late
morphological types (Mollenhoff 2004).

Fig. 3 shows the photometric scaling relations
for the bulges. The relation between the luminosity
of morphologically decoupled bulges and their Sersic
index is displayed in Fig. 3a. The dashed line
in Fig. 3 corresponds to mean relation
for the bulges of
almost 1000 galaxies constructed from the data of
Gadotti (2009) (see Eq. (25) and Fig. 17 in Mosenkov
et al. (2014)).
As it can be seen from Fig. 3a, the bulges
we consider roughly follow the trend typical for the regular
galaxies.

\begin{figure}
\centering
\includegraphics[width=16.5cm, angle=-90, clip=]{fig3.ps}\\
\caption{The scaling relations for the polar bulges in the $r$ band (circles).
(a) The luminosity vs Sersic index relation; the dashed
line indicates the mean relation for the bulges of regular galaxies (see text).
(b) The luminosity vs effective radius (in kpc) relation;
the mean relation for normal bulges is indicated by the dashed line.
(c) The Kormendy relation; the dashed line indicates the
relation for the E/S0 galaxies in the Coma cluster (Houghton et al. 2012);
the dots are the characteristics of the polar-ring central
galaxies (Reshetnikov and Combes 2015).}
\end{figure}

Fig. 3b illustrates the luminosity -- size relation
for the decoupled bulges in the comparison with the mean
relation for the bulges plotted from the data of
Gadotti (2009) (Eq. (31) in Mosenkov et al. (2014)).
The polar bulges follow the relation drawn for the regular bulges,
with a slight shift towards the larger sizes. It is difficult
to judge the significance of this shift because our
sample is small and in addition, the normal bulges
show a large scatter on the luminosity -- size diagram.

The position of the bulges in the Kormendy relation
are shown in Fig. 3c. The bulges are located
on this plane approximately along the canonical relation
for the E/S0 galaxies with a shift from it and from
the central polar-ring galaxies towards the smaller size
and down to the lower central brightness.
There are faint early-type
galaxies and bulges (Capaccioli et al. 1992) in
this region (approximately at log\,$r_\mathrm{e} \leq 0.5$). The polar
bulges are at the edge of the region where the polar-ring
galaxies are located and are comparable in their
characteristics to the faintest and smallest host
objects in the polar rings (Fig. 3c).

\section{Conclusions}

We presented a sample of nine edge-on spiral
galaxies whose bulges were decoupled from their
disks in the sense that the major axes of the bulges are highly inclined
to the disk midplane. Note that we selected the objects
based on a purely morphological analysis; therefore,
there can be chance projections of two galaxies
among our objects. The galaxies with such decoupled,
or polar, bulges are extremely rare. Based on
the number of galaxies we examined when searching
for the polar bulges, we can roughly estimate their
relative fraction among the spiral galaxies as $\sim 10^{-3}$.

Our analysis of the photometric characteristics for
the galactic disks in the polar-bulge galaxies showed them to be similar to the
disks of ordinary late-type (Sc--Sd) spiral galaxies.
If we disregard the elongation along
the minor axes of the galaxies, the bulges also
appear quite typical. The polar
bulges are similar to the faintest and smallest of
the central galaxies in the polar-ring galaxies.

The peculiarity of the polar-bulge galaxies is that
galaxies with strong stellar disk warps are encountered
relatively frequently among them. Since the warp
is often interpreted as a consequence of an external perturbation
and external accretion, these factors may
also be responsible for the formation of polar bulges.
As has been noted in Introduction, kinematically
decoupled subsystems rotating almost orthogonally
to the galactic disks are observed in two out of three
previously studied polar bulges. The existence of such
subsystems is usually associated with the external
accretion. Consequently, the detection of
such kinematically decoupled structures in the polar
bulges of other galaxies would give a good argument
for the hypothesis of the formation of decoupled
bulges under external accretion.

However, the question of what is the secondary
structure in the polar-bulge galaxies, their disks or bulges,
remains unclear. Their similarity to the polar-ring galaxies,
where the formation of polar disks/rings is well
reproduced in the models with the external accretion
from another galaxy or from a cloud of intergalactic gas
(Reshetnikov and Sotnikova 1997; Bournaud and
Combes 2003; Maccio et al. 2006), may suggest that
it is the disk in such objects that is secondary. From
this viewpoint, the polar-bulge galaxies are simply
the polar-ring galaxies with faint central galaxies. On
the other hand, as was discussed by Matthews and
de Grijs (2004), the opposite scenario, whereby a
pre-existing spiral galaxy captures a small elliptical
galaxy in the polar plane to become its bulge, is also
possible. Less exotic scenarios can also be
considered. In general, it is clear that new observing
data and realistic numerical simulations are necessary to
clarify the nature of these unique objects.

\bigskip
\section*{Acknowledgments}
This study was financially supported by the Russian
Foundation for Basic Research (project nos. 14-02-810
and 13-02-00416). DB is supported by grant
RSCF-14-22-00041.

\section*{REFERENCES}

\hspace{0.4cm}1.\,K.N.\,Abazajian, J.K.\,Adelman-McCarthy, M.A.\,Agueros, et al., Astrophys. J. Suppl. Ser.
182, 543 (2015).

2.\,P.\,Aguirre, J.M.\,Uson, and L.D.\,Matthews, Rev. Mex. Astron. Astrofis. 35, 201 (2009).

3.\,S.\,Alam, F.D.\,Albareti, C.A.\,Prieto, F.\,Anders,
F.S.\,Anderson, T.\,Anderton, B.H.\,Andrews, E.\,Armengaut,
et al., Astrophys. J. Suppl. Ser. 219, 12A (2015).

4.\,H.B.\,Ann and J.-C.\,Park, New Astron. 11, 293 (2006).

5.\,P.N.\,Appleton and C.\,Struck-Marcell, Fundam. Cosm. Phys. 16, 111 (1996).

6.\,E.\,Athanassoula, Mon. Not. R. Astron. Soc. 358, 1477 (2005).

7.\,K.\,Bekki, Astrophys. J. 499, 635 (1998).

8.\,F.\,Bertola, E.M.\,Corsini, J.C.\,Vega\,\,Beltran,
A.\,Pizzella, M.\,Sarzi, M.\,Cappellari, and S.J.\,Funes,
Astrophys. J. 519, L127 (1999).

9.\,D.V.\,Bizyaev, S.J.\,Kautsch, A.V.\,Mosenkov,
V.P.\,Reshetnikov, N.Ya.\,Sotnikova, N.V.\,Yablokova,
and R.W.\,Hillyer, Astrophys. J. 787, ID 24 (2014).

10.\,F.\,Bournaud and F.\,Combes, Astron. Astrophys. 401, 817 (2003).

11.\,M.\,Capaccioli, N.\,Caon, and M.\,D’Onofrio, Mon. Not. R. Astron. Soc. 259, 323 (1992).

12.\,E.M.\,Corsini, J.\,M\'endez-Abreu, N.\,Pastorello,
E.\,dalla Bonta, L.\,Morelli, A.\,Beifiori, A.\,Pizzella,
and F.\,Bertola, Mon. Not. R. Astron. Soc. 423, L79 (2012).

13.\,P.\,Erwin, Astrophys. J. 799, id. 226 (2015).

14.\,P.\,Erwin, R.P.\,Saglia, M.\,Fabricius, J.\,Thomas,
N.\,Nowak, S.\,Rusli, R.\,Bender, J.C.\,Vega Beltran,
and J.E.\,Beckman,Mon. Not. R. Astron. Soc. 446, 4039 (2015).

15.\,D.A.\,Gadotti, Mon. Not. R. Astron. Soc. 393, 1531 (2009).

16.\,A.M.\,Garcia, Astron. Astrophys. Suppl. Ser. 100, 47 (1993).

17.\,R.C.W.\,Houghton, R.L.\,Davies, E.\,dalla Bonita,
and R.\,Masters, Mon. Not. R. Astron. Soc. 423, 256 (2012).

18.\,J.\,Kormendy, Galactic Bulges, Ed. by E.\,Laurikainen,
R.F.\,Peletier, and D.A.\,Gadotti (Springer, New York, 2015).

19.\,J.\,Kormendy and R.C.\,Kennicutt, Ann. Rev. Astron. Astrophys. 42, 603 (2004).

20.\,P.C.\,van der Kruit and L.\,Searl, Astron. Astrophys. 95, 105 (1981).

21.\,A.V.\,Maccio, B.\,Moore, and J.\,Stadel, Astrophys. J. 636, L25 (2006).

22.\,I.\,Martin-Navarro, J.\,Bakos, I.\,Trujillo, J.H.\,Knapen,
E.\,Athanasoula, A.\,Bosma, S.\,Comeron, B.G.\,Elmegreen, et al., Mon. Not. R. Astron.
Soc. 427, 1102 (2012).

23.\,L.D.\,Matthews and R.\,de Grijs, Astron. J. 128, 137 (2004).

24.\,J.\,M\'endez-Abreu, Galactic Bulges, Ed. by E.\,Laurikainen,
R.F.\,Peletier, and D.A.\,Gadotti (Springer, New York, 2015).

25.\,A.F.J.\,Moffat, Astron. Astrophys. 3, 455 (1969).

26.\,A.V.\,Moiseev, K.I.\,Smirnova, A.A.\,Smirnova, and V.P.\,Reshetnikov,
Mon.Not. R. Astron. Soc. 418, 244 (2011).

27.\,C.\,Mollenhoff, Astron. Astrophys. 415, 63 (2004).

28.\,A.V.\,Mosenkov, N.Ya.\,Sotnikova, and V.P.\,Reshetnikov,
Mon. Not. R. Astron. Soc. 441, 1066 (2014).

29.\,A.V.\,Mosenkov, N.Ya.\,Sotnikova, V.P.\,Reshet- nikov,
D.V.\,Bizyaev, and S.J.\,Kautsch, Mon. Not. R. Astron. Soc. 451, 2376 (2015).

30.\,Ch.Y.\,Peng, L.C.\,Ho, Ch.D.\,Impey, and H.-W.\,Rix,
Astron. J. 139, 2097 (2010).

31.\,A.\,Renzini, The Formation of Galactic Bulges, Ed. by C.M.\,Carollo,
H.C.\,Ferguson, and R.F.G.\,Wyse (Cambridge Univ. Press, Cambridge, 1999).

32.\,V.\,Reshetnikov and F.\,Combes, Astron. Astrophys. 337, 9 (1998).

33.\,V.\,Reshetnikov and F.\,Combes, Astron. Astrophys. Suppl. Ser. 138, 101 (1999).

34.\,V.\,Reshetnikov and F.\,Combes, Mon. Not. R. Astron. Soc. 447, 2287 (2015).

35.\,V.\,Reshetnikov and N.\,Sotnikova, Astron. Astrophys. 325, 933 (1997).

36.\,M.\,Sarzi, E.M.\,Corsini, A.\,Pizzella, J.C.\,Vega Beltran, M.\,Cappellari,
J.G.\,Funes, and F.\,Bertola, Astron. Astrophys. 360, 439 (2000).

37.\,E.F.\,Schlafly and D.P.\,Finkbeiner, Astrophys. J. 737, ID 103 (2011).

38.\,J.L.\,Sersic, Bol. Asoc. Argentina Astron. 6, 41 (1963).

39.\,L.A.M.\,Tasca and S.D.M.\,White, Astron. Astrophys. 530, A106 (2011).

40.\,I.\,Trujillo, J.A.L.\,Aguerri, J.\,Cepa, and C.M.\,Gutierrez,
Mon. Not. R. Astron. Soc. 328, 977 (2001).

\end{document}

Astrophysics of Galaxies (astro-ph.GA); Cosmology and Nongalactic Astrophysics (astro-ph.CO)

Polar-bulge galaxies