We present a realistic study for electronic and magnetic properties in dilute magnetic semiconductor (Ga,Mn)As. A multi-orbital Haldane-Anderson model parameterized by density-functional calculations is presented and solved with the Hirsch-Fye quantum Monte Carlo algorithm. Results well reproduce experimental results in the dilute limit. When the chemical potential is located between the top of the valence band and an impurity bound state, a long-range ferromagnetic correlations between the impurities, mediated by antiferromagnetic impurity-host couplings, are drastically developed. We observe an anisotropic character in local density of states at the impurity-bound-state energy, which is consistent with the STM measurements. The presented combined approach thus offers a firm starting point for realistic calculations of the various family of dilute magnetic semiconductors. The discovery of ferromagnetism in dilute magnetic semiconductor (DMS) materials, represented by a (Ga,Mn)-As system, has created much activity in the field of spintronics. [1][2][3][4] In the low-doping regime (( 1%), (Ga,Mn)-As is insulating with a clear experimental evidence for the presence of an Mn-induced impurity band located at 110 meV above the valence band.5) The position of the impurity level and the impurity-induced carriers are key quantities in generating the ferromagnetism. Various model studies [6][7][8][9][10] have clarified that this ferromagnetic correlations between the impurities are based on the polarization mechanism of the impurity-induced carriers.However, it is still an open question whether this mechanism really works in real materials. In order to clarify this problem and step forward to the challenge for predicting new functional properties of the DMS materials and their realistic designs, many tasks such as nonempirical electronic-band structure of the host compound, reliable estimations for impurity-host hybridizations, and rigorous treatments for local electron correlations at the impurity site are required. A combined approach of density functional theory (DFT) 11) and quantum Monte Carlo (QMC) technique 12) is clearly useful for this purpose. The hybridization or mixing parameters between the host and impurity orbitals is a crucial parameter characterizing material differences and therefore this parameter should be estimated on the basis of ab initio density-functional calculations. On the other hand, density-functional calculations often fail to describe local electron correlations, so this problem must be treated with more accurate solvers such as the QMC technique. Critical comparisons with the experimental data are needed for a check whether this approach is really reliable in describing the details of the ''low-energy'' electronic and magnetic properties and/or whether it has an extended applicability for many DMS materials.In this letter, we present a comprehensive study for the generation mechanism of the ferromagnetism of a dilute (Ga,Mn)As system. A multi-orbital Haldane-Anderson model 13) is employed...