Molecular high-order harmonic generation (MHOHG) and molecular orbital ionization rates are calculated for the nonsymmetric OCS and symmetric CS 2 molecules using numerical solutions of Kohn-Sham (KS) equations of time-dependent density functional theory in the nonlinear nonperturbative regime of laser-molecule interactions for different laser-molecule orientations and intensities. It is found that the ionization of inner-shell KS molecular orbitals contributes significantly to the ionization and MHOHG processes for intensities I 3.5 × 10 14 W/cm 2 . Ionization rate maxima correspond to the alignment of maximum KS orbital densities with the laser pulse polarization instead of orbital ionization potentials. Furthermore, degeneracies of orbitals are removed as a function of laser-molecule angle, thus affecting ionization rates, the MHOHG spectra, and their polarizations, the latter allowing for identifying inner-orbital ionization.We report in this paper the nonlinear nonperturbative response to short (few cycles) intense laser pulses of the nonsymmetric OCS and symmetric CS 2 molecules. The timedependent density functional theory (TDDFT) Kohn-Sham (KS) equation for each occupied molecular orbital (MO) in these molecules was discretized in space using finitedifference (FD) grid techniques as reported earlier for CO 2 and other molecules [1,2]. We have used the Van Leeuwen and Baerends [3] potential (LB94), which introduces a gradient correction to the local density approximation exchange correlation so as to reproduce correctly the Coulomb asymptotic behavior of the potential. It reproduces KS MO negative energies [1,4,5] nearly equal to the ionization potentials (IP) obtained from photoelectron spectra, which are only accurately given normally by many-body Dyson orbitals [6]. The bond length has been set to its experimental [7-10] value, i.e., OCS bonds length are CS = 0.156 nm, CO = 0.1157 nm, and CS = 0.155 nm for CS 2, and these molecules are oriented along the z axis. The time-dependent ionization probability P i,σ (t) of an individual spin-orbital is calculated [2] aswhere N i,σ (t) = ψ i,σ (r,t)|ψ i,σ (r,t) is the time-dependent population (survival probability) of the i,σ th spin population. The number of electrons inside the grid of volume V and remaining population arewhere n(r,0) = 2 N σ i |ψ iσ (r,0)| 2 is the total number of electrons in the initial state (without external field). The total ionization yield probability γ (t f ), computed at the final time t f of the propagation, is then obtained from the difference between initial and final probability β(t f ) in Eq. (2) and * andre.bandrauk@usherbrooke.ca tabulated in Table I for different laser-molecule angles θ and intensities.The power spectrum of the dipole accelerationd z (t) in a given direction yields the predicted MHOHG spectra [11,12] Recent work has shown that the most accurate method for generating the MHOHG spectra when dealing with molecular systems [5,11,12] is to calculate the dipole acceleration d z (t) from the exact time-dependent fu...