abstract: Non-orthogonal multiple access (NOMA) has attracted much recent attention owing to its capability for improving the system spectral efficiency in wireless communications. Deploying NOMA in heterogeneous network can satisfy users' explosive data traffic requirements, and NOMA will likely play an important role in the fifth-generation (5G) mobile communication networks. However, NOMA brings new technical challenges on resource allocation due to the mutual cross-tier interference in heterogeneous networks. In this article, to study the tradeoff between data rate performance and energy consumption in NOMA, we examine the problem of energy-efficient user scheduling and power optimization in 5G NOMA heterogeneous networks. The energy-efficient user scheduling and power allocation schemes are introduced for the downlink 5G NOMA heterogeneous network for perfect and imperfect channel state information (CSI) respectively. Simulation results show that the resource allocation schemes can significantly increase the energy efficiency of 5G NOMA heterogeneous network for both cases of perfect CSI and imperfect CSI.
author: Haijun Zhang, , Fang Fang, , Julian Cheng, , Keping Long, , Wei Wang, , and Victor C.M. Leung,  [^1]
title: Energy-Efficient Resource Allocation in NOMA Heterogeneous Networks

# Introduction

In order to meet the explosive growth of wireless traffic data requirement and overcome the shortage of the frequency resource, non-orthogonal multiple access (NOMA) has been proposed as a promising technique in the fifth generation (5G) of wireless communication systems . Different from orthogonal frequency division multiple access (OFDMA) technique using orthogonal subcarriers, NOMA enables multiple users to be multiplexed on the same frequency band by applying the successive interference cancelation (SIC) technique at the receivers . According to the NOMA protocol, the user who has worse channel state information (CSI) is assigned more power than the user who has better CSI . This protocol guarantees the weak users' access by successively decoding the signals from the received signals in the decreasing order of power. It was shown that the NOMA scheme can achieve higher spectral efficiency and higher data rate than the traditional OFDMA scheme .

Driven by the rapid increase of wireless terminal equipments and wide usage of mobile Internet, heterogeneous networks have emerged as one of the most promising network infrastructures to provide high system throughput and large coverage of indoor and cell edge scenarios in 5G wireless communication systems . In such an architecture, a macrocell is overlaid by several small cells, e.g., microcell, picocell and femtocell, to significantly improve the system throughput and the spectral efficiency. For heterogeneous networks, frequency band sharing between macrocell and small cells is viable, and it is also more efficient to reuse the frequency bands within a macrocell . However, the cross-tier interference can severely degrade the quality of the wireless transmission. The advantage of the heterogeneous network also comes with some fundamental challenges on cross-tier interference mitigation, which have been studied in previous research by some techniques, e.g. precoding technique and resource management . In this article, we apply NOMA with SIC technology to heterogeneous networks in order to alleviate the cross-tier interference and improve the system energy efficiency via resource optimization. In the traditional OFDMA heterogeneous networks, the frequency band can be divided into several sub frequency bands, and users in the macro cell and small cells are assigned to different sub frequency bands in order to avoid the cross-tier interference . However, in NOMA heterogeneous networks, SIC is applied at the receivers to allow multiple users to be multiplexed on the same sub frequency band and the strong users with higher CSI can cancel the interference from the weak user with lower CSI multiplexed on the same sub frequency band.

With the overwhelming increase of the traffic data and mobile devices, the energy cost has rapidly increased and become an important issue in the green cellular network because of the increasing amount of $\rm CO_2$ emission levels caused by energy consumption . Thus, fast growing energy consumption and limited global energy resources are the important motivation for the research on energy efficient communications in NOMA heterogeneous networks.

Different from the existing works in NOMA networks, in this article, we focus on the energy efficient wireless communications for the downlink 5G NOMA heterogeneous networks via resource optimization. More specifically, we study user scheduling and power allocation by considering both perfect CSI and imperfect CSI. We first present network architecture of a 5G NOMA heterogeneous network. Then we describe an effective resource allocation scheme in 5G NOMA heterogeneous networks to improve the system energy efficiency with perfect CSI. A user scheduling and power optimization scheme with imperfect CSI is also introduced. Finally, we conclude this article.

# Resource Optimization for NOMA Heterogeneous Networks

In NOMA heterogeneous networks, a macrocell with one macro base station (BS) is overlaid by several small cells where each small cell implements NOMA technique. In this system, macrocell users (MUEs) and small cell users (SUEs) can be served on the same sub frequency band at the same time to achieve high spectral efficiency. In other words, the overlay frequency band sharing model implies that macrocell users can access the sub frequency band freely and efficiently under the minimum data rate requirement determined by the quality of service (QoS). High data rate can be achieved by allocating more power to the UEs; however, this will incur high power consumption. Hence, finding the tradeoff between the data rate improvement and power consumption has become an important and inevitable trend in the next generation of wireless communication systems. For energy efficient wireless communications, we aim to design a resource scheduling scheme to balance the system sum rate improvement and the power consumption.

## NOMA Overview

We now present an overview of the resource allocation for NOMA systems. Different from the conventional OFDMA technique, multiple users in NOMA system can share the same frequency band instead of using orthogonal spectrum. By applying SIC technique at the receivers, a user in NOMA system can successively extract the intended signal from the received aggregated signals. In , SIC was exploited to decode the signals from the received signals in the decreasing order of power. Therefore, the user who has higher channel gain can cancel the interference signal from the user who has lower channel gain. NOMA can be combined with MIMO to further improve capacity . Particularly, the pioneering work of beamspace MIMO-NOMA was firstly proposed to integrate NOMA for millimeter wave massive MIMO systems in , where a dynamic power allocation scheme with low complexity was proposed by solving the challenging joint power optimization problem. The statistical channel information was applied at the transmitter to achieve the maximum of system energy efficiency. The energy-efficient resource allocation for a downlink NOMA single cell network has been investigated in . The authors in applied matching theory to assign the users to the different subchannels, and utilized difference of convex functions (DC) programming to allocate different levels of power to the multiplexed users and the subchannels. To reduce the complexity of the decoding at receivers, the number of users allocated on the same sub-band was limited to two. The NOMA system equipped with the proposed resource allocation scheme can achieve higher sum rate and higher energy efficiency than the traditional OFDMA system. In OFDMA heterogeneous networks, users in the same small cell are not allowed to share the same frequency band in order to avoid the inter-cell interference .

In the NOMA heterogeneous networks, NOMA can be implemented in each small cell, which means that the users in the same small cell are allowed to share the same sub frequency band. Different from NOMA single cell system, users in the same small cell can share the same sub-band with the users in macrocell. The interference signal from the users in the small cell can be detected with the superposition coding technique and removed by the strong users in the macrocell via SIC at the receivers.

## NOMA Heterogeneous Network Model

In this system, we assume each small cell occupies one sub frequency band. It means that all the users in the same small cell can be multiplexed on the same sub frequency band. The MUEs can be co-multiplexed on the same sub frequency band with SUEs. Each BS is responsible for frequency bands assignment (user scheduling), access procedures and power allocation.

In the following, we use sub-band to replace sub frequency band for the sake of simplicity. In NOMA heterogeneous networks, each UE can receive different signals from the macro BS and a small BS simultaneously. Furthermore, each BS transmits signals to UEs by using superposition coding. As shown in Fig. 1, we consider a downlink NOMA heterogeneous network where a primary macrocell hosts a macro base station (MBS) with $K$ overlaid small cells. Let $M$ denote the total number of active MUEs in the macrocell. Each small cell hosts one small BS (SBS) and $F$ SUEs that are randomly distributed within the small cell. The MBS transmits signals to $M$ MUEs through these $K$ sub-bands. The total bandwidth of the system is denoted as $B$, which is divided into $K$ sub-bands and each sub-band occupies bandwidth of $B/K$. The transmit power for the entire system is $P_s$. Because the transmit power for MBS should be higher than SBS, we assume that both the MBS and all the SBSs have transmit power $P_s/2$ and each SBS has transmit power $P_s/(2K)$. A block fading channel is considered in the system model, where the channel fading of each subcarrier is assumed to be the same within a sub-band, but it varies independently across different sub-bands. In this system, the MUEs can share the same frequency band with the small cells. The cross-tier interference caused by the primary macrocell will be effectively mitigated by applying the NOMA technique. Note that interference between small cells is neglected because different small cells will occupy different sub-bands and the sub-bands are independent with each other. The small cell is allocated lower power than the macrocell and the signal within the small cell can suffer severe wall penetration loss . Therefore, the co-tier inter-small cell interference is negligible when compared with cross-tier interference . We assume that MBS and SBS are connected by the wired links, thus the backhaul will not be considered in this article.

The resource optimization problem is formulated as an energy efficiency maximization problem for the downlink NOMA heterogeneous networks. Our goal is to maximize the total system sum rate with a unit power through resource scheduling. In other words, in order to guarantee amount of system data rate, we aim to reduce the power consumption as much as possible. In our study, energy efficiency is defined as the ratio of the system sum rate and total power consumption. The total power consumption includes both circuitry-consumed power of the mobile devices and transmission power. Therefore, we use bits per Joule to measure the system energy efficiency performance. The energy-efficient resource optimization problem can be formulated as the maximization of the energy efficiency. In this problem, user scheduling and power allocation are coupled each other. In order to solve this problem efficiently, the resource allocation for energy efficiency maximization problem can be divided into user scheduling and power allocation sub-problems. In the following, we will discuss the resource optimization for NOMA heterogeneous networks based on two cases: perfect CSI and imperfect CSI.

# Energy Efficient Resource Optimization for NOMA Heterogeneous Networks with Perfect CSI

We first consider the BSs have knowledge of channel state information. For user scheduling, matching theory can be utilized to match different UEs with different sub-bands. In the proposed energy efficient subchannel assignment framework for a downlink NOMA single cell network, the utility function of each user is defined as the summation of each subchannel energy efficiency which can be described as the ratio of the subchannel data rate and power consumption. The power consumption for each subchannel includes circuit power and transmission power on this subchannel. In the NOMA heterogeneous networks, we define our user scheduling optimization sub-problem as the summation of small cells energy efficiency. For user scheduling scheme design, we first assume equal power allocation on each sub-band. This user scheduling optimization is a non-convex mixed integer programming problem. The global optimal user scheduling scheme can be obtained by an exhaustive search. However, the complexity of the exhaustive search is extremely high and this approach is infeasible in practice. Assume we have $M$ MUEs and $N$ small cells. The scheduler needs to search $\frac{{(M)!}}{{{2^N}}}$ combinations. Therefore, a suboptimal user scheduling scheme can be designed for this system by using a matching theory. In the sub-optimal scheme, the complexity of the worst case is $O({N^2M})$. Taking natural logarithm of the complexity, the logarithm complexity is $O(\ln M)$. Since $O(\ln M)<O(M\ln M)$. The complexity of the sub-optimal scheme is much less than the exhaustive searching. In the proposed user scheduling scheme, we consider all the UEs as a user set, and consider $K$ sub-bands as a sub-band set. The user scheduling process can be described as the two-sided dynamic matching between the user set and the sub-band set. For each sub-band, we set a maximum number of users that can be multiplexed on one sub-band in order to reduce the decoding complexity at the receivers. Therefore, the energy efficiency of the overall system can be formulated as the maximization of energy efficiency under the constraints of maximal transmission power, the minimum user data rate as well as the maximum number of users can be allocated on one sub-band. Fig. 2 describes the resource optimization procedure including user scheduling and power allocation.

## User Scheduling

In user scheduling, we consider the SUEs scheduling and the MUEs scheduling separately. Since NOMA is implemented in each small cell, all the SUEs in one small cell can only multiplex on one sub-band. Then SUEs scheduling can also be called small cells scheduling. In the small cells scheduling, since the scheduler has the perfect CSI, each sub-band has preferred sequence of small cells to assign. This preferred sequence is generated according to the channel gains between SBS and SUEs within each small cell. Assume, for the time being, the equal power allocation among SUEs in the small cell, the energy efficiency of the small cell on each sub-band can be calculated by the scheduler. The scheduler chooses the sub-band with the highest channel gain and schedule the corresponding small cell on this sub-band. After all the small cells have occupied sub-bands, the scheduler schedules the MUEs on these occupied sub-bands, and the MUEs scheduling depends on the energy efficiency that the MUE can provide for the sub-bands. Therefore, the SUEs and MUEs co-multiplex on these sub-bands. In the MUEs scheduling scheme, each sub-band has a preferred allocation list of MUEs in a decreasing order of the channel gains between MBS and MUEs. For each sub-band, the scheduler finds the highest channel gain and schedule its corresponding MUE onto this sub-band. According to the condition of the maximum number of sub-band users, the MUE allocates on the sub-band directly when the number of the multiplexed users on this sub-band is less than the maximum number. However, when the number of the multiplexed users on this sub-band equals the maximum, the scheduler compares the energy efficiency of users and assign the sub-band to the users who can provide the highest energy efficiency. In the energy efficiency calculation, we assume the power of MUEs equals the power of small cell on this sub-band and each SUE has the equal power. The MUEs matching process terminates when all the UEs have been multiplexed on the sub-bands. The optimal sub-band matching can only be obtained by searching all matching combinations of users and sub-bands. This procedure comes with extremely high complexity and is not practical. The above proposed sub-optimal user scheduling has lower complexity than the optimal solution.

## Power Allocation

In user scheduling, we have assumed the equal power allocation across sub-bands and equal power for small cells and MUEs on the same sub-band. In order to further improve the system energy efficiency, sophisticated power allocation schemes can be proposed to replace the equal power allocation scheme. In the NOMA heterogeneous networks, the energy efficiency optimization problem is non-convex. To solve this problem, DC programming approach and time-sharing method can be used to transform the non-convex combinatorial optimization problem into a convex optimization problem. In this article, we exploit fractional transmit power allocation (FTPA) to further allocate powers among users and across sub-bands . Due to its low computational complexity, FTPA has been widely adopted in the OFDMA and NOMA systems.

Figure 3 shows the performance of the system energy efficiency versus the number of the small cells. The bandwidth is limited to 10 MHz, and we set the peak power of the entire system to be 20 W and circuit power to be 0.1 W on each sub-band. Figure 3 indicates that the total energy efficiency increases when the number of the small cells increases. From Fig. 3, we can observe that NOMA heterogeneous network (NOMA HetNet) outperforms the OFDM heterogeneous network (OFDM HetNet) because OFDM HetNet cannot fully utilize the spectrum resources.

Fig. 4 shows the cumulative distribution function (CDF) of user average data rate comparison between NOMA HetNet and OFDM HetNet with equal power allocation scheme. In Fig. 4, the CDF of user average data rate in NOMA HetNet is always below that in OFDM HetNet. This means that, given a target data rate, the number of users whose average data rate is larger than this target data rate in NOMA HetNet is less than that in NOMA HetNet system. This is because NOMA technique has higher spectrum efficiency than the OFDMA scheme.

# Energy Efficient Resource Optimization for NOMA Heterogeneous Network with Imperfect CSI

In this subsection, we study the energy efficient resource allocation for the downlink NOMA heterogeneous network with imperfect CSI. In practice, the perfect CSI at the transmitter is difficult to achieve due to channel estimation errors, feedback and quantization errors. To maximize the energy efficiency, we can formulate the energy efficient resource allocation as a probabilistic mixed non-convex optimization problem, under the constraints of outage probability limitation, the maximum transmitted power and the number of multiplexed users on one sub-band.

The resource allocation scheme is optimized based on the imperfect CSI. To solve this problem, we decouple the problem into user scheduling and power allocation sub-problems separately to maximize the system energy efficiency. We assume that the BSs have the estimated value of the CSI. In this situation, the user data rate may not meet the minimum data requirement determined by QoS. Therefore, an outage probability requirement is considered for the resource scheduling to maximize the system energy efficiency. The energy efficient optimization problem is formulated as a probabilistic mixed non-convex optimization problem. In order to efficiently solve the optimization problem, we can first transform the probabilistic mixed non-convex optimization problem to a non-probabilistic optimization problem. The outage probability constraint can be transformed to minimum power constraints for each UEs sharing the same sub-band. In this transformed problem, we can treat user scheduling and power allocation separately. The user scheduling starts with assigning equal power allocation on sub-bands. The optimal solution to the user scheduling sub-problem is challenging to obtain in practice as it requires to search all the possible combinations. A sub-optimal user scheduling algorithm can be obtained by using the estimated CSI. In order to cancel the interference from the other small cells, we assign each small cell a sub-band. MUEs multiplex on these sub-bands, and hence the small cell users not interfere with each other. The MUE with high channel gain is chosen to multiplex on its corresponding sub-band. For the each sub-band, we set the maximum number of multiplexed users to reduce the complexity of decoding at the receivers. The user scheduling terminates if there is no MUE left to be allocated. During the user scheduling, the power allocation for the multiplexed users on one sub-band can be determined by FTPA. The complexity of the proposed sub-optimal algorithm is less than the optimal solution that can only be obtained by the exhaustive search. Given the user scheduling scheme, the underlying optimization problem can be shown to be convex function respect to the power variable under the constraint of the maximum power of the system. Therefore, a unique optimal solution can be found by a gradient assisted binary search algorithm. Therefore, the new power allocation scheme can replace the equal power allocation to further improve the system energy efficiency.

Figure 5 shows the energy efficiency versus the number of the small cells with the channel gain estimation error variance 0.05. As shown, the energy efficiency increases when the number of the small cells grows. From Fig. 5, the NOMA HetNet outperforms OFDM HetNet in terms of energy efficiency, and the energy efficiency of NOMA HetNet with FTPA is higher than that of scheme with equal power allocation.

Figure 6 shows the energy efficiency of the NOMA HetNet performance with different estimation error variances. In both OFDM HetNet and NOMA HetNet schemes, the system energy efficiency deteriorates when the error variances increases. The NOMA HetNet always outperforms OFDM HetNet. For example, when the number of small cells is 20, the energy efficiency of NOMA HetNet with FTPA power allocation scheme is 27% more than that of the OFDM HetNet with the equal power (EQ) allocation scheme under the imperfect channel CSI.

# Open Issues, and Challenges

Resource allocation is one of the methods to improve the energy efficiency of 5G NOMA networks. Small cell in HetNet can also be energy efficient because low power is typically required due to the short distance between transmitter and receiver. However, there are still many challenges and open issues in energy-efficient 5G NOMA networks. 1) NOMA with MIMO: Not only user scheduling and power allocation, but also beamforming can be considered in the energy efficient resource allocation. Furthermore, different configuration of antennas for macrocell and small cell in HetNet makes the problem more complex. Finding effective ways to combine NOMA and massive MIMO is a challenging but important topic in future research works. 2) NOMA with energy harvesting: To save more energy, 5G NOMA HetNet is expected to harvest energy from solar, wind, thermoelectric, and so on. One open problem is how to manage those renewable energy source because they are intermittent over time and space. 3) NOMA with game theory: User pairing is key research direction in energy-efficient 5G NOMA networks. Matching game theory can be an effective tools to distributed model and optimize the user pairing of energy-efficient 5G NOMA networks. 4) NOMA with open/closed access small cells: There are three access modes in small cell networks, open access, closed access, and hybrid access. Different access mode of NOMA small cells requires different energy saving methods.

# Conclusion

We introduced the concept of NOMA heterogeneous network and examined the problem of energy efficient user scheduling and power allocation in 5G NOMA heterogeneous network by considering both perfect CSI and imperfect CSI. We solved the resource optimization problem with the help of convex optimization. Based on the obtained results, we proposed a two-step user scheduling and power optimization scheme. The effectiveness of the proposed schemes was compared to the existing scheme and verified by simulations in terms of energy efficiency. It is envisioned that a joint user scheduling and power allocation approach can further enhance the energy efficiency of the overall system.

[^1]: Haijun Zhang and Keping Long are with the Engineering and Technology Research Center for Convergence Networks, University of Science and Technology Beijing, Beijing, 100083, China (e-mail: firstname.lastname@example.com, email@example.com).

    Fang Fang and Julian Cheng are with School of Engineering, The University of British Columbia, Kelowna, BC, V1X 1V7, Canada (e-mail: email@example.com, email@example.com).

    Wei Wang is with the College of Information Science and Electronic Engineering, Zhejiang University, Hangzhou 310027, China, and also with the Zhejiang Provincial Key Laboratory of Information Processing, Communication and Networking, Hangzhou 310027, China (e-mail: firstname.lastname@example.com).

    Victor C.M. Leung is with the Department of Electrical and Computer Engineering, The University of British Columbia, Vancouver, BC V6T 1Z4 Canada (e-mail: firstname.lastname@example.com).


abstract: We derive the dependence of the fraction of passive central galaxies on the mass of their supermassive black holes for a sample of over 400,000 SDSS galaxies at z $<$ 0.2. Our large sample contains galaxies in a wide range of environments, with stellar masses $8 < {\rm log}(M_{*}/M_{\odot}) < 12$, spanning the entire morphological spectrum from pure disks to spheroids. We derive estimates for the black hole masses from measured central velocity dispersions and bulge masses, using a variety of published empirical relationships. We find a very strong dependence of the passive fraction on black hole mass, which is largely unaffected by the details of the black hole mass estimate. Moreover, the passive fraction relationship with black hole mass remains strong and tight even at fixed values of galaxy stellar mass ($M_{*}$), dark matter halo mass ($M_{\rm halo}$), and bulge-to-total stellar mass ratio ($B/T$). Whereas, the passive fraction dependence on $M_{*}$, $M_{\rm halo}$ and $B/T$ is weak at fixed $M_{BH}$. These observations show that, for central galaxies, $M_{BH}$ is the strongest correlator with the passive fraction, consistent with quenching from AGN feedback.
author: Asa F. L. Bluck$^{1,*}$, Sara L. Ellison$^{1}$, David R. Patton$^{2}$, Luc Simard$^{3}$, J. Trevor Mendel$^{4}$, Hossein Teimoorinia$^{1}$, Jorge Moreno$^{5,6}$, and Else Starkenburg$^{7,8}$
$^1$ Department of Physics and Astronomy, University of Victoria, 3800 Finnerty Road, Victoria, British Columbia, V8P 1A1, Canada
$^2$ Department of Physics and Astronomy, Trent University, 1600 West Bank Drive, Peterborough, Ontario, K9J 7B8, Canada
$^3$ National Research Council of Canada, Herzberg Institute of Astrophysics, 5071 West Saanich Road, Victoria, British Columbia, V9E 2E7, Canada
$^4$ Max-Planck-Institut für Extraterrestrische Physik (MPE), Giessenbachstrasse, D-85748 Garching, Germany
$^5$ Department of Physics and Astronomy, California State Polytechnic University Pomona, Pomona, CA 91768, USA
$^6$ TAPIR, Mailcode 350-17, California Institute of Technology, Pasadena, CA 91125, USA
$^7$ Schwarzschild Fellow, Leibniz-Institut für Astrophysik Potsdam (AIP), An der Sternwarte 16, 14482 Potsdam, Germany
$^8$ CIfAR Global Scholar
$*$ Email: email@example.com
title: Why do galaxies stop forming stars? I. The passive fraction - black hole mass relation for central galaxies

# Introduction

Understanding why galaxies stop forming stars is an important unresolved question in galaxy evolution. Only $\sim$`<!-- -->`{=html}10% of baryons reside within galaxies (e.g. Fukugita & Peebles 2004, Shull et al. 2012), yet since galaxies lie at nodes in the cosmic web corresponding to local minima in the gravitational potential well, naively one would expect far more baryons to collate in galaxies, ultimately forming more stars. Theoretical models offer a wide range of solutions to this problem, relying on the physics of gas, stars, and black hole accretion disks as so called 'baryonic feedback'. However, observational studies are required to test these models and provide evidence for their range and applicability.

Since the discovery that most galaxies contain a supermassive black hole (e.g. Kormendy & Richstone 1995), the energy released from forming these objects has become a popular mechanism for regulating gas flows and star formation in simulations, particularly for more massive galaxies (e.g. Croton et al. 2006, Bower et al. 2006, 2008, Somerville et al. 2008, Guo et al. 2011, Vogelsberger et al. 2014a,b, Schaye et al. 2015). In fact, substantial feedback from accretion around supermassive black holes is required in cosmological semi-analytic models, semi-empirical models, and hydrodynamical simulations to achieve the steep slope of the high-mass-end of the galaxy stellar mass function (e.g. Henriques et al. 2014, Schaye et al. 2015, Vogelsberger et al. 2014a,b).

Observationally, predominantly only very indirect means for linking the star formation in large populations of galaxies to their AGN have been found, e.g. the strong dependence of the passive fraction on galaxy stellar mass (e.g. Baldry et al. 2006, Peng et al. 2010, 2012) and galaxy structure (e.g. Driver et al. 2006, Cameron et al. 2009, Cameron & Driver 2009, Bluck et al. 2014). These galaxy properties both correlate (weakly) with black hole mass and, hence, the total energy released in forming the black hole (e.g. Silk & Rees 1998, Fabian 1999, Bluck et al. 2011, 2014). More direct measurements of AGN driven winds in galaxies and radio bubbles in galaxy haloes have provided evidence for the mechanisms of AGN feedback, but typically only for a very small number of galaxies (e.g. McNamara et al. 2000, Nulsen et al. 2005, McNamara et al. 2007, Dunn et al. 2010, Fabian 2012, Cicone et al. 2013).

More recent work has linked the passive fraction of large populations of galaxies to the central density within 1 kpc (Cheung et al. 2012, Fang et al. 2013, Woo et al. 2014) and to the mass of the galactic bulge (Bluck et al. 2014, Lang et al. 2014). Both of these galaxy properties are expected to correlate strongly with the mass of the central black hole (e.g. Haring & Rix 2004, McConnell & Ma 2013) and hence may provide qualitative support for AGN feedback driven quenching.

We find in Bluck et al. (2014) that 'bulge mass is king' in the sense that bulge mass is a tighter and steeper correlator to the passive fraction for centrals than any other variable considered (including stellar mass, halo mass, disk mass, local density, and galaxy structure). We argue for an explanation of this observation through AGN feedback, relying on the tight relationships between the bulge and central black hole (e.g. Haring & Rix 2004). The primary motivation for this Letter is to expand on that work by testing the role of AGN feedback in quenching central galaxies more directly.

In this Letter we estimate black hole masses from their bulge masses and central velocity dispersions, using a wide range of published empirical relationships with dynamical measurements of black hole masses (e.g. McConnell & Ma 2013). From this, we derive the relationship between passive fraction and estimated black hole mass. Our goal is to test whether AGN feedback is a viable route for quenching galaxies, and test whether or not it is dominant for a large sample of central galaxies. Throughout this Letter we assume a $\Lambda$CDM cosmology with H$_{0}$ = 70 km s$^{-1}$, $\Omega_{m}$ = 0.3, $\Omega_{\Lambda}$ = 0.7, and adopt AB magnitude units.

# Data Overview & Passive Fraction

We use the Sloan Digital Sky Survey Data Release 7 (SDSS DR7, Abazajian et al. 2009) spectroscopic sample as our data source. From this we collate a sample of 538046 galaxies (423480 centrals and 114566 satellites) with $8 < {\rm log}(M_{*}/M_{\odot}) < 12$ at z $<$ 0.2. In this Letter (the first paper in the series) we concentrate on the central galaxies. The star formation rates (SFR) for these galaxies are calculated from emission lines for non-AGN star forming galaxies and from the strength of the 4000 Å break for non-emission line galaxies and AGN (Brinchmann et al. 2004). A fibre correction is applied based on galaxy colour and magnitude outside the aperture. The stellar masses for the galaxies, and their component disks and spheroids, are derived in Mendel et al. (2014), based on SED fitting to a dual Sérsic fit of the $ugriz$ wavebands (Simard et al. 2011). An n = 4 bulge and n = 1 disk model is used, and we test the reliability of this approach in Bluck et al. (2014) via model data. We define the galaxy structure (or morphology) to be $B/T = M_{\rm bulge} / M_{*}$, where $M_{*}$ is the total stellar mass of the galaxy.

The velocity dispersions are derived from the widths of absorption lines taken from Bernardi et al. (2003) with an updated method implemented as in Bernardi et al. (2007) to the later data releases. Velocity dispersions from absorption lines with a S/N $<$ 3.5 are discarded from our sample, and those with $\sigma$ $<$ 70 km/s are deselected for some analyses, due to the instrumental resolution of the SDSS spectra. We also restrict our final sample to galaxies with an error on the velocity dispersion of $\sigma_{\rm err}$ $<$ 50 km/s. We recover $\sim$ 80 % of our parent sample which pass these data quality cuts.

The halo masses used in this Letter are derived from an abundance matching technique applied to the total stellar mass of the group or cluster (from Yang et al. 2007, 2008, 2009). Testing of the group finding algorithm on model galaxies from the Millennium Simulation (Springel et al. 2005) showed that over 90 % of galaxies are correctly assigned to groups at $M_{\rm halo} > 10^{12} M_{\odot}$. Using these group catalogues, centrals are defined as the most massive galaxy in the group, with satellites being any other galaxy within the group.

We follow the prescription for defining the passive fraction in Bluck et al. (2014). We define a galaxy to be passive if it is forming stars at a rate a factor of ten times lower than (emission line, non-AGN) star forming galaxies matched at the same stellar mass and redshift. This cleanly divides the two peaks of the bimodal distribution in SFR at all mass ranges. The passive fraction is the ratio of passive-to-total galaxies in each binning. We correct for the flux limit of the SDSS by weighting each galaxy in the passive fraction by the inverse of the volume over which its $ugriz$ magnitudes would pass all of the selection criteria. The errors on the passive fraction are computed in this work via the jack-knife technique. See Bluck et al. (2014) §2 & 3 for full details on these data and techniques.

# Black Hole Masses

Today there are less than 100 reliable dynamical measurements of central black hole masses in existence (see McConnell & Ma 2013 for a compilation). This is far too few to investigate the dependence of the passive fraction on black hole mass directly, due in part to the statistical nature of the passive fraction and the need for a wide range in galaxy properties. Fortunately, there exist tight relationships between dynamical measurements of central black hole masses and the properties of their host galaxies, particularly the bulge or spheroid in which the black hole resides (e.g. Magorrian et al. 1998; Ferrarese & Merritt 2000; Gebhardt et al. 2000; Haring & Rix 2004). On a case by case basis, using these relationships to estimate a given galaxy's black hole mass would lead to a high level of uncertainty, but for very large samples the average values should yield reliable measurements, in the regime where the relationships are tested.

In this work, we use a number of the galaxy - black hole relationships from the literature to estimate black hole masses. Specifically, we use central velocity dispersion ($\sigma_{c}$), stellar mass of the bulge or spheroid ($M_{\rm bulge}$), and combinations of these two. We also test separating the sample into early- and late-type systems and using separate scaling laws for each (as suggested in McConnell et al. 2011). The formulae used to estimate the black hole masses, and the references for these, are all presented in Table 1.

For the velocity dispersions, we first make an aperture correction, so that all measurements are made at the same effective aperture. We use the formula in Jorgensen et al. (1995), specifically calculating:

$$\sigma_{c} = \big(\frac{R_{c}}{R_{\rm app}}\big)^{-0.04} \sigma_{\rm app}$$

where the central radius, $R_{c}$, is chosen to be in line with the measurements made in the literature, with typical values of $R_{e}$ or $R_{e}/8$ used. $R_{e}$ is the bulge (or spheroid) effective radius, calculated in our sample through Sérsic index fitting (see Simard et al. 2011 for full details). We note that the aperture correction only affects the final mass estimate by typically $<$`<!-- -->`{=html}10%. For the bulge masses no such correction is necessary, and thus, this provides a natural test to the reliability of the method.

Since our goal is to present the passive fraction dependence on black hole mass for all galaxies, we do not restrict the sample to pure spheroids. Thus, there is a substantial population of disks present in our sample. This poses no significant difficulty for measuring the black hole masses via bulge mass (but see Bluck et al. 2014 appendices for a discussion on their reliability). However, for the velocity dispersions, contamination in their measurement from rotating disks could be an issue. We test this in the next section by looking at early- and late-type galaxies separately, and comparing estimates of $M_{BH}$ from velocity dispersion methods to bulge mass methods. We additionally test the impact of the choice of black hole mass calibration by comparing a number of different techniques, noting a high level of consistency between them (see next section).

```latex
\begin{table*}
\caption{Black Hole Mass Estimates}
 \label{tab1}
 \begin{tabular}{@{}cccccc}
  \hline
\hline
Source & log($M_{BH}/M_{\odot}$) =  & Morphology & Dispersion (dex) & $N_{\rm Gal}$ & Fig. 1 Colour \\
\hline
Hopkins et al. (2007) & 0.54 $\times$ log($M_{\rm bulge} / 10^{11} M_{\odot}$) + 2.18 $\times$ log($\sigma_{c}$ / 200 km/s) + 8.24 & No Cut & 0.22 & 38 & Blue \\
%\hline
Ferrarese \& Merritt (2000) & 4.80 $\times$ log($\sigma_{c}$) - 2.90 & No Cut & 0.31 & 12 & Green\\
%\hline
Haring \& Rix (2004) & 1.12 $\times$ log($M_{\rm bulge} / 10^{11} M_{\odot}$) + 8.2 & No Cut & 0.33 & 30 & Cyan\\
%\hline
McConnell \& Ma (2013) & 5.64 $\times$ log($\sigma_{c}$ / 200 km/s) + 8.32 & No Cut & 0.38 & 72 & Orange\\
%\hline
McConnell \& Ma (2013) & 5.20 $\times$ log($\sigma_{c}$ / 200 km/s) + 8.39 & ETGs & 0.34 & 53 & Red$^{*}$\\
%\hline
McConnell \& Ma (2013) & 5.06 $\times$ log($\sigma_{c}$ / 200 km/s) + 8.07 & LTGs & 0.46 & 19 & Red$^{*}$\\
\hline
\end{tabular}
$^{*}$ for the ETG/ LTG samples, galaxies defined as early types ($B/T >$ 0.5) are fit with the ETG formula and those defined as late types ($B/T <$ 0.5) are fit with the LTG formula, they are then recombined into the same sample, and shown in red in Fig. 1. $N_{\rm Gal}$ is the number of galaxies used in the dynamical fit.
\end{table*}
```

# Results

## The $f_{\rm passive} - M_{BH}$ Relationship for Central Galaxies

In this section we compute, for the first time, the relationship between the fraction of passive galaxies and the mass of their central supermassive black holes, for a population of over 400,000 central galaxies (see Fig. 1). The left panel shows the $f_{\rm passive} - M_{BH}$ relationship for all central galaxies, which are defined as the most massive group members in the catalogues of Yang et al. (2009). We show this relationship via a number of different black hole mass calibrations, summarized in Table 1. These include relationships from various authors based on central velocity dispersion, bulge mass, and a combination of the two. For the full morphological sample, there is remarkably good agreement between all of the methods used for estimating black hole masses.

The $f_{\rm passive} - M_{BH}$ relationship is very steep for central galaxies, with galaxies with low black hole masses being predominantly star forming and galaxies with high black hole masses being predominantly passive. The cross-over mass, where 50% of galaxies are passive occurs at $M_{BH} \sim 10^{7.5} M_{\odot}$. Values of black hole mass $< 10^{6} M_{\odot}$ are highly uncertain due to the $M_{BH} - \sigma$ , $M_{\rm bulge}$ relations not extending into this domain. Additionally, velocity dispersions $<$ 70 km/s are less reliable than for the rest of the sample (e.g. Bernardi et al. 2007), and this corresponds to approximately $M_{BH}$ $< 10^{6} M_{\odot}$ as well. For these reasons, we advise readers to treat the low $M_{BH}$ values more as upper limits (this threshold is indicated with a dashed vertical line in all relevant plots).

In the middle and right panels of Fig. 1, we show the passive fraction - black hole mass relationship for centrals split into early-types ($B/T >$ 0.5) and late-types ($B/T <$ 0.5) respectively (computed in Mendel et al. 2014). For the early-type sample, all of the black hole mass estimates lead to very similar passive fraction relationships, as with the full sample. However, for the late-type galaxies, there is some disagreement between methods which use velocity dispersion alone and those which incorporate bulge mass. Specifically, the passive fraction is lower for velocity dispersion estimates of black hole mass than with bulge mass estimates. This is most probably explained by the velocity dispersions being contaminated by disk rotation into the plane of the sky, and by measurements of the absorption line widths being contaminated by multiple emission lines obscuring the true continua. Both of these effects can lead to an over-estimate of the black hole mass, and hence lead to lower passive fractions for the measured black hole masses from velocity dispersion methods compared to bulge mass methods.

For bulge mass only methods contamination of the kinematic measurement in disks is clearly not an issue, but interestingly, in the case where we use $\sigma_{c}$ and $M_{\rm bulge}$ together (Hopkins et al. 2007), the contamination is also negligible (compare blue and cyan lines in Fig. 1, right panel). Thus, we posit that to estimate black hole masses for the full morphological sample (i.e. including disk galaxies) a measure of the mass of the bulge is *essential* to include. For the remainder of the analyses in this Letter we focus on the black hole mass estimates from bulge mass and velocity dispersion together (from Hopkins et al. 2007, see top row of Table 1), where kinematic contamination is not a significant issue. This technique draws on multiple complementary measurements of the central regions of galaxies, from spectroscopic kinematics and photometric profile and SED fitting. We note that all of our conclusions remain unchanged if we restrict to early-type galaxies only and use velocity dispersion methods to compute our black hole masses. The advantage of incorporating the bulge mass is to leave the sample complete.

In an upcoming paper in the series (Paper II, Bluck et al. in prep.) we show the passive fraction relationship with black hole mass for satellites. We find that, whereas centrals with low black hole masses are essentially all star forming, satellites with low black hole masses are more frequently passive. This hints at the need for additional quenching mechanisms to AGN feedback for satellites, but not necessarily for centrals.

## Is the $f_{\rm passive} - M_{BH}$ Relationship Fundamental for Central Galaxies?

The steep dependence of the passive fraction on supermassive black hole mass for centrals (Fig. 1, left panel) is consistent with the quenching of central galaxies being driven by AGN feedback. However, given that there are observationally well known, and theoretically expected, correlations between the black hole mass (derived from $M_{\rm bulge}$ and/or $\sigma_{c}$) and other galaxy properties (e.g. $M_{\rm halo}$, $M_{*}$, $B/T$), it is not immediately clear whether the $f_{\rm passive} - M_{BH}$ relationship is fundamental or not (i.e. whether it is derivative of some other, more fundamental, parameter).

To test the possibility that another galaxy property is more fundamental to passivity than black hole mass, we look at the $f_{\rm passive} - M_{BH}$ relationship at a fixed values of the following galaxies properties: $M_{\rm halo}$, $M_{*}$, and $B/T$ (presented in Fig. 2, top row). We find that the $f_{\rm passive} - M_{BH}$ relationship remains tight and steep even at a fixed halo mass, stellar mass and galaxy structure ($B/T$). Out of these three variables, we find the greatest variation at a fixed $M_{BH}$ with $B/T$, suggesting that the structure of a galaxy is an important secondary consideration for whether a galaxy will be forming stars or not. In Fig. 2 (bottom row) we present the $f_{\rm passive} - M_{\rm halo}, M_{*}$, and $B/T$ relationships, each at fixed $M_{BH}$. These are all remarkably flat compared to the $f_{\rm passive} - M_{BH}$ relationship, suggesting that it must be more fundamental than these alternatives.

Particularly, the $f_{\rm passive} - M_{\rm halo}$ relationship, at a fixed $M_{BH}$, is almost entirely flat, indicating essentially no dependence of galaxy quenching on the mass of the halo, once the correlation with black hole mass is accounted for. This has a profound implication for models of galaxy quenching. Methods that derive the energy to do work on the gas in a galaxy from the dark matter halo cannot be driving the quenching of centrals (e.g. halo mass quenching, Dekel & Birnboim 2006, Dekel et al. 2009, Woo et al. 2013). However, at high $M_{BH}$ values, there is some hint of a subtle positive correlation of the passive fraction with $M_{\rm halo}$, which may be explained if AGN feedback is more efficient in higher mass haloes (e.g. Dekel et al. 2014).

The $f_{\rm passive} - M_{*}$ relationship at fixed $M_{BH}$ is either flat, or in certain $M_{BH}$ ranges, *negative*. This striking result clearly indicates that stellar mass is not the primary driver to passivity in central galaxies (as advocated in Baldry et al. 2006, Peng et al. 2010, 2012). This result further suggests that the quenching of central galaxies cannot be driven by stellar or supernova feedback, which would both correlate with $M_{*}$ (as integrated star formation rate) more strongly than with $M_{BH}$. The negative trend is most probably explained by the fact that increasing $M_{*}$ at fixed $M_{BH}$ results in an increase in $M_{\rm disk}$, and hence a decrease in $B/T$. This follows directly from the $M_{BH} - M_{\rm bulge}$ relationship, and we have seen already that $B/T$ is an important *secondary* correlator to the passive fraction.

The $f_{\rm passive} - M_{BH}$ relationship is also steeper and tighter than the $f_{\rm passive} - B/T$ relationship (which was suggested to be primary in Driver et al. 2006, Cameron et al. 2009). This implies that major galaxy mergers (which can destroy in situ disks, e.g., Cole et al. 2000) are not necessary for quenching central galaxies, but substantial central black holes are. However, see Hopkins et al. (2013) for evidence that disks can regrow in the merger process complicating this issue. Additionally, morphological quenching (e.g. Martig et al. 2009) cannot be the primary driver to the quenching of star formation in centrals, since black hole mass correlates much stronger with the passive fraction than galaxy structure. See Bluck et al. (2014) §5 for a more thorough discussion on the role of various processes in quenching central galaxies.

We consider the relative role of the supermassive black hole in quenching central and satellites galaxies compared to many other galaxy properties (including bulge mass, disk mass, and local density) in an upcoming work in this series (Paper III, Teimoorinia et al., in prep.). Black hole mass remains the most important single correlator to the passive fraction in this larger analysis as well.

# Conclusions

In this Letter, we present for the first time the relationship between the fraction of passive galaxies and the mass of their supermassive black holes, for a sample of over 400,000 central galaxies at z $<$ 0.2 (Fig. 1). We find that the $f_{\rm passive} - M_{BH}$ relationship is very steep, with the vast majority of galaxies with low black hole masses being star forming and the vast majority of galaxies with high black hole masses being passive. We find a cross-over mass (where 50% of galaxies are passive) at $M_{BH} \sim 10^{7.5} M_{\odot}$. The steep dependence of the passive fraction on black hole mass is consistent with AGN feedback being the dominant mechanism for quenching central galaxies. We test splitting our sample into early- and late-types, finding generally very good accord between all of the black hole mass estimates, except for late-type galaxies. Here we find that incorporating bulge mass into the black hole mass estimate is essential because velocity dispersion alone methods are vulnerable to contamination from disk rotation.

Finally, we find that black hole mass is a more fundamental correlator to the passive fraction than total stellar mass, group halo mass, or galaxy structure (Fig. 2). This is highly suggestive of the quenching of central galaxies being driven by AGN feedback. However, it is not inconceivable that the mechanism which grows the central region of galaxies (and presumably their supermassive black holes as well) is directly responsible for the quenching of central galaxies. Yet this would have to occur in such a way as to leave the central density as the primary correlator to the passive fraction, with $M_{\rm halo}$, $M_{*}$, and $B/T$ being substantially less important. We conclude that this is less likely than AGN feedback being responsible.


abstract: Large amount of observational spectroscopic data are recently available from different observational projects, like Sloan Digital Sky Survey. It's become more urgent to identify white dwarfs stars based on data itself i.e. without modelling white dwarf atmospheres. In particular, existing methods of white dwarfs identification presented in Kleinman et al. (2004) and in Eisenstein et al. (2006) did not allow to find all the white dwarfs in examined data. We intend to test various criteria of searching for white dwarf candidates, based on photometric and spectral features.
address: $^1$University of Warsaw Astronomical Observatory, Al. Ujazdowskie 4, 00-478 Warsaw, Poland; $^2$The Andrzej Sołtan Institute for Nuclear Studies, Hoża 69, 00-681 Warsaw, Poland
author: Mirosław Należyty$^{1}$, Agnieszka Majczyna$^{1,2}$, Anna Ciechanowska$^{1}$ and Jerzy Madej$^{1}$
title: Searching for white dwarfs candidates in Sloan Digital Sky Survey Data

# Introduction

Because of their physical properties, white dwarfs stars are difficult to observe and identification. The first version of the McCook & Sion catalogue (1987) contains 1279 white dwarfs only. Recently, observational techniques was considerable improve, so the amount of photometric and spectral data rapidly grows. For this reason also the number of identified white dwarfs grows, but it is very needed to elaborate efficient identification and classification methods.

Until now, many of the white dwarf samples were selected, based on different selection criteria. For example Fleming et al. (1996) published a white dwarfs catalogue based on ROSAT data (Trumper 1983, 1992). A classification was made by searching for position coincidence between ROSAT sources and hot white dwarfs from McCook & Sion catalogue. Sources with no optical counterparts were marked as white dwarf if was visible only in the softest X-ray energy range. White dwarfs in EUV, selected also from the ROSAT all-sky survey by Marsh et al. (1997a,b), 129 white dwarf stars from Bergeron, Saffer & Liebert (1992), 200 white dwarfs from the Palomar Green survey (Green, Schmidt & Liebert (1986) analyzed by Liebert & Bergeron (1995) - it is only a few samples of these objects.

The Sloan Digital Sky Survey project (York et al. 2000) provide us large amount of very interesting data. Based on Data Release 1 and 4 (DR1 and DR4, respectively) there were made two catalogues by Kleinman et al. (2004) and Eisenstein et al. (2006), which contain 2551 and 9316 objects, respectively. Used identification methods did not allow to find all white dwarfs with whole range of effective temperatures and types in data sample. Kleinman et al. (2004) chose white dwarf candidates based on dereddened colour indexes $u-g$ and $g-r$ and a magnitude in $u$ filter. Selected white dwarf candidates was next verified by visual inspection and prescribed its type - DA, DB, DZ etc. After this procedure, spectra of selected objects were fitted by theoretical DA and DB spectra.

Eisenstein et al. (2006) disposed larger amount of data so they used more automatic selection procedure. They select stars fulfilled applicable criteria (see details in Eisenstein et al. 2006). In this manner they obtain sample of blue stars, which observational spectra were fitted by theoretical spectra. Stars with parameters characteristic for white dwarfs were next verified and classified during different tests including visual inspection. This method has some limitations, for example does not work for stars cooler than 8000K.

We intend to find a method of the white dwarf candidates selection - or, if it is possible, white dwarfs itself - based only on photometric and spectral features, chosen to obtain the best distinguishing between white dwarf stars and other type objects. Before we start fitting theoretical model atmospheres to observing spectra.

Not for the first time we have used data from the Sloan Digital Sky Survey (York et al. 2000). We chose Data Release 5 (DR5) and we selected spectral data for 10% randomly choosing objects. In this manner we obtained very imposing amount of nearly 16000 different type objects, including galaxies, quasars, stars, and of course white dwarfs. Additionally we also used white dwarfs identified in SDSS DR4 by Eisenstein et al. (2006), treating these stars on the one hand as a potential help to improve our white dwarfs selection method, on the other hand as a control data.

# Balmer lines and their parameters

The most characteristic feature of many DA white dwarf spectra is a presence of wide, hydrogen absorption Balmer lines (except for the hottest objects). It would appear that it should be a good selection criteria of degenerated stars with high surface gravity, but not every white dwarf have wide Balmer lines. Nevertheless we decided to check this

possibility. Because of large amount (15998) of objects in our sample we had to create proper software, which automatically finds hydrogen lines from Balmer series, and then calculates its various parameters like line widths at given depth, fluxes in lines etc. Fig. 1 shows one of the obtained results, a dependency of FWHM of Balmer $H_{\alpha}$ line on colour index $u-g$ for objects with $H_{\alpha}$, $H_{\beta}$, and $H_{\gamma}$ lines present in their spectra. Clearly seeing bimodal structure is not a by-product of our Balmer lines searching, but it is a result of a bimodal distribution of the colour index $u-g$. A location of the white dwarfs identified in SDSS DR4 (Kleinman et al. 2004), and belonged to our sample with detected Balmer lines, denoted by black dots, suggests, that all of the white dwarfs lies in the left part of this diagram, with colour indexes $u-g < 0.7$. In general, this is not a truth. Visual inspection shows, that some amount of white dwarf stars also lies in the right part of the diagram, and have larger values of colour indexes $u-g$. However, our requirement for the presence of hydrogen Balmer lines in spectra reduces our sample to practically star-like type objects only. In fact, it is decreases a number of white dwarf candidates from 15998 to 4613. Unfortunately, reduced in this way sample does not contain the hottest white dwarf stars, because there is no Balmer lines visible in their spectra.

# Fluxes in given ranges

Colour indexes based on $ugriz$ photometry are not very useful for us, generally because the filters are too wide. So we decide to define own artificial filters, i.e. own wavelength ranges, in which we calculated fluxes. Choosing our filters we intended to not contain any Balmer

lines. In Fig. 2 we present a dependency of two flux logarithms, calculated in two arbitrarily chosen wavelength ranges: $4550-4700$Å and $8000-8500$Å for all of 15998 objects in our sample. Objects with detected $H_{\alpha}$, $H_{\beta}$, and $H_{\gamma}$ lines (dark gray dots) occupy well defined area, so we decided to watch them closer. Fig. 3 shows logarithms of fluxes in the same wavelength ranges as in Fig. 2, but for 4613 objects with Balmer lines only. It is not very difficult to see some structure in this diagram. At least part of the objects divides to, not so bad, separate regions. This impression is additionally magnified by black dots showing positions of the white dwarfs from DR4, which occupy one of the mentioned above regions only. We decided to approximate a boundary between these two areas by the linear function $y = 0.94118x + 0.33235$, with empirically chosen constants. Preliminary results of the visual inspection shows, that spectra of almost all of the 2202 objects located in the lower-right region looks like white dwarf spectra. Because of relatively poor separation of these two areas, between 2411 objects belonging to the upper-left region we still could find relatively small amount of white dwarf stars, located close to the borderline.

# Conclusions

Our method of selecting white dwarf candidates is based on searching for objects with hydrogen Balmer lines visible in their spectra, and on flux calculations in well selected wavelength ranges. Preliminary results show that our method allows to select quite complete sample of white dwarf stars, under the above assumptions. Of course this method should be tested in detail, and it needs some improvements (for example increasing separation between two areas in Fig. 3, probably by choosing better wavelength ranges). It is quite possible, that using this method we shall be able to find also the hottest white dwarfs, which are too hot to show Balmer lines in their spectra, although we will need the other criteria to select star-like type objects. We plan to do this in future.

## Acknowledgments

This work has been supported by the Polish Ministry of Science and Higher Education grant No. N N203 4061 33. We also thank Institut d'Estudis Espacials de Catalunya (IEEC) for financial support.

Funding for creation and distribution of the SDSS Archive has been provided by the Alfred P. Sloan Foundation, the Participating Institutions, the National Aeronautics and Space Administration, the National Science Foundation, the US Department of Energy, the Japanese Monbukagakusho, and Max Planck Society. The SDSS Web site is http://www.sdss.org/ .

The SDSS is managed by the Astrophysical Research Consortium (ARC) for the Participanting Institutions. The Participanting Institutions are The University of Chicago, Femilab, the Institute for Advanced Study, the Japan Participation Group, the Johns Hopkins University, Los Alamos National Laboratory, the Max-Planck-Institute for Astronomy (MPIA), the Max-Planck-Institute for Astrophysics (MPA), New Mexico State University of Pittsburgh, Princeton University, the United States Naval Observatory, and the University of Washington.

# References