The ionization of two-active-electron systems by intense laser fields is investigated theoretically. In comparison with time-dependent Hartree-Fock and exact two electron simulation, we show that the ionization rate is overestimated in SAE approximation. A modified single-active-electron model is formulated by taking into account of the dynamical core polarization. Applying the new approach to Ca atoms, it is found that the polarization of the core can be considered instantaneous and the large polarizability of the cation suppresses the ionization by 50% while the photoelectron cut-off energy increases slightly. The existed tunneling ionization formulation can be corrected analytically by considering core polarization.PACS numbers: 33.80. Rv, 42.50.Hz, 42.65.Re Various of non-perturbative phenomena occurring during atom-laser interactions are started with single ionization, e.g., above threshold ionization (ATI) and high harmonic generation (HHG). Although they have been successfully interpreted by the rescattering model based on single active electron (SAE) approximation (see Reviews e.g., [1,2]), detailed examination showed that multielectron effects are embedded in the photon and electron spectra [3][4][5][6][7][8][9]. It is found that high-order harmonic generation (HHG) from molecules records interference of different channels suggesting more than one molecular orbitals are involved [3] and electron rearrangement is occuring [4], which is certainly beyond the scope of the SAE theory. On the other hand, two-electron events such as non-sequential double ionization can not be explained either without considering the electron-electron interaction [10]. It is thus desirable to examine in details the multielectron effects occurring in the ionization of atomic systems beyond SAE.The single ionization of atoms in strong laser fields can be pictured as tunneling of one electron through the barrier formed by the atomic potential and the laseratom dipole interaction. The Keldysh parameter measures the ratio of tunneling time to the optical period, γ = I p /2U p , where I p = κ 2 /2 is the ionization potential and U p = E 2 /4ω 2 is the ponderomotive energy of a free electron in a laser field of strength E and frequency of ω. When γ < 1, tunnel ionization occurs so rapid that the electric field can be considered as a static field at each instant. The so-called adiabatic approximation is the root of Ammosov-Delone-Krainov (ADK) -like theories [11] for obtaining ionization rates. Based on this picture, the rate is mainly determined by the unitless quantity κ 3 /E with κ 3 representing the atomic field strength at the classical radius of the electron motion.It is obvious that the adiabatic approximation will break down if the atomic potential acting on the tunneling electron is varing sooner than the tunneling time. For more than one electron systems, the core can be polarized by the laser fields, hence the atomic potential is time-varying. In the case of absence of resonant excitation, the polarization is instantaneou...