The mechanisms of anisotropic near-IR tunnel ionization and high-order harmonic generation (HHG) in a CO molecule are theoretically investigated by using the multiconfiguration timedependent Hartree-Fock (MCTDHF) method developed for the simulation of multielectron dynamics of molecules. The multielectron dynamics obtained by numerically solving the equations of motion (EOMs) in the MCTDHF method is converted to a single orbital picture in the natural orbital representation where the first-order reduced density matrix is diagonalized. The ionization through each natural orbital is examined and the process of HHG is classified into different optical paths designated by a combinations of initial, intermediate and final natural orbitals. The EOMs for natural spin-orbitals are also derived within the framework of the MCTDHF, which maintains the first-order reduced density matrix to be a diagonal one throughout the time propagation of a many-electron wave function. The orbital dependent, timedependent effective potentials that govern the dynamics of respective time-dependent natural orbitals are deduced from the derived EOMs, of which the temporal variation can be used to interpret the motion of the electron density associated with each natural spin-orbital. The roles of the orbital shape, multiorbital ionization, linear Stark effect and multielectron interaction in the ionization and HHG of a CO molecule are revealed by the effective potentials obtained. When the laser electric field points to the nucleus O from C, tunnel ionization from the C atom side is enhanced; a hump structure originating from multielectron interaction is then formed on the top of the field-induced distorted barrier of the HOMO effective potential. This hump formation, responsible for the directional anisotropy of tunnel ionization, restrains the influence of the linear Stark effect on the energy shifts of bound states.

To theoretically demonstrate the binding of a positron to a nonpolar or small dipole molecule, we have calculated the vibrational averaged positron affinity (PA) values along the harmonic asymmetric stretching vibrational coordinate with the configuration interaction level of multi-component molecular orbital method for CXY (X, Y = O, S, and Se) molecules. For CSe2 and CSSe molecules, a positron can even be attached at the equilibrium structures, due to the effect of the induced dipole moment with large polarizability values. For a CS2 molecule, the positive PA value is obtained at the lowest vibrational excited state in our scheme. Although there is no direct experimental evidence for the positron-binding to CO2, COS, and COSe molecules, we have predicted positron-binding for these molecules at higher vibrational excited states.

A systematic quantum Monte Carlo study of 2p atoms (C, N, O) and 3p atoms (Si, P, S) is performed to investigate the influence of correlation on the interpretation of Hund's multiplicity rule, which is an extension of our previous study of the carbon atom [J. Chem. Phys. 121, 7144 (2004)] to heavier atoms. The accuracy in the present study is significantly improved as compared with the previous study. A detailed analysis of the correlation contribution to individual energy components of the total energy is given beyond the self-consistent Hartree-Fock calculation. The stability of the highest spin-multiplicity state of all the atoms is ascribed to the greater electron-nucleus attraction energy that is gained at the cost of increasing the electron-electron repulsion energy as well as the kinetic energy. The present study demonstrates that correlation does not change the above conclusion due to the Hartree-Fock theory to support Boyd's less screening mechanism.

Using the framework of multiconfiguration theory, where the wavefunction Φ(t) of a many-electron system at time t is expanded as Φ(t)=Σ(I)C(I)(t)Φ(I)(t) in terms of electron configurations {Φ(I)(t)}, we divided the total electronic energy E(t) as E(t)=Σ(I)|C(I)(t)|(2)E(I)(t) . Here E(I)(t) is the instantaneous phase changes of C(I)(t) regarded as a configurational energy associated with Φ(I)(t). We then newly defined two types of time-dependent states: (i) a state at which the rates of population transfer among configurations are all zero; (ii) a state at which {E(I)(t)} associated with the quantum phases of C(I)(t) are all the same. We call the former time-dependent state a classical stationary state by analogy with the stationary (steady) states of classical reaction rate equations and the latter one a quantum stationary state. The conditions (i) and (ii) are satisfied simultaneously for the conventional stationary state in quantum mechanics. We numerically found for a LiH molecule interacting with a near-infrared (IR) field ε(t) that the condition (i) is satisfied whenever the average velocity of electrons is zero and the condition (ii) is satisfied whenever the average acceleration is zero. We also derived the chemical potentials μ(j)(t) for time-dependent natural orbitals ϕ(j)(t) of a many-electron system. The analysis of the electron dynamics of LiH indicated that the temporal change in Δμ(j)(t) ≡ μ(j)(t) + ε(t) · d(j)(t) - μ(j)(0) correlates with the motion of the dipole moment of ϕ(j)(t), d(j)(t). The values Δμ(j)(t) are much larger than the energy ζ(j)(t) directly supplied to ϕ(j)(t) by the field, suggesting that valence electrons exchange energy with inner shell electrons. For H2 in an intense near-IR field, the ionization efficiency of ϕ(j)(t) is correlated with Δμ(j)(t). Comparing Δμ(j)(t) to ζ(j)(t), we found that energy accepting orbitals of Δμ(j)(t) > ζ(j)(t) indicate high ionization efficiency. The difference between Δμ(j)(t) and ζ(j)(t) is significantly affected by electron-electron interactions in real time.

A unified interpretation of Hund's first and second rules for 2p (C, N, O) and 3p (Si, P, S) atoms is given by Hartree-Fock (HF) and multiconfiguration Hartree-Fock (MCHF) methods. Both methods exactly satisfy the virial theorem, in principle, which enables one to analyze individual components of the total energy E(=T+V(en)+V(ee)), where T, V(en), and V(ee) are the kinetic, the electron-nucleus attraction, and the electron-electron repulsion energies, respectively. The correct interpretation for each of the two rules can only be achieved under the condition of the virial theorem 2T+V=0 by investigating how V(en) and V(ee) interplay to attain the lower total potential energy V(=V(en)+V(ee)). The stabilization of the more stable states for all the 2p and 3p atoms is ascribed to a greater V(en) that is caused by contraction of the valence orbitals accompanied with slight expansion of the core orbitals. The contraction of the valence orbitals for the two rules is a consequence of reducing the Hartree screening of the nucleus at short interelectronic distances. The reduced screening in the first rule is due to a greater amount of Fermi hole contributions in the state with the highest total spin-angular momentum S. The reduced screening in the second rule is due to the fact that two valence electrons are more likely to be on opposite sides of the nucleus in the state with the highest total orbital-angular momentum L. For each of the two rules, the inclusion of correlation does not qualitatively change the HF interpretation, but HF overestimates the energy difference ∣ΔE∣ between two levels being compared. The magnitude of the correlation energy is significantly larger for the lower L states than for the higher L states since two valence electrons in the lower L states are less likely to be on opposite sides of the nucleus. The MCHF evaluation of ∣ΔE∣ is in excellent agreement with experiment. The present HF and MCHF calculations demonstrate the above statements that were originally given by Katriel [Theor. Chem. Acta 23, 309 (1972); 26, 163 (1972)]. We have, for the first time, analyzed the correlation-induced changes in the radial density distribution for the excited LS terms of the 2p and 3p atoms as well as for the ground LS term.