Higher-order equation-of-motion coupled-cluster methods for ionization processes

Compact algebraic equations defining the equation-of-motion coupled-cluster (EOM-CC) methods for ionization potentials (IP-EOM-CC) have been derived and computer implemented by virtue of a symbolic algebra system largely automating these processes. Models with connected cluster excitation operators truncated after double, triple, or quadruple level and with linear ionization operators truncated after two-hole-one-particle (2h1p), three-hole-two-particle (3h2p), or four-hole-three-particle (4h3p) level (abbreviated as IP-EOM-CCSD, CCSDT, and CCSDTQ, respectively) have been realized into parallel algorithms taking advantage of spin, spatial, and permutation symmetries with optimal size dependence of the computational costs. They are based on spin-orbital formalisms and can describe both alpha and beta ionizations from open-shell (doublet, triplet, etc.) reference states into ionized states with various spin magnetic quantum numbers. The application of these methods to Koopmans and satellite ionizations of N2 and CO (with the ambiguity due to finite basis sets eliminated by extrapolation) has shown that IP-EOM-CCSD frequently accounts for orbital relaxation inadequately and displays errors exceeding a couple of eV. However, these errors can be systematically reduced to tenths or even hundredths of an eV by IP-EOM-CCSDT or CCSDTQ. Comparison of spectroscopic parameters of the FH+ and NH+ radicals between IP-EOM-CC and experiments has also underscored the importance of higher-order IP-EOM-CC treatments. For instance, the harmonic frequencies of the A 2Sigma- state of NH+ are predicted to be 1285, 1723, and 1705 cm(-1) by IP-EOM-CCSD, CCSDT, and CCSDTQ, respectively, as compared to the observed value of 1707 cm(-1). The small adiabatic energy separation (observed 0.04 eV) between the X 2Pi and a 4Sigma- states of NH+ also requires IP-EOM-CCSDTQ for a quantitative prediction (0.06 eV) when the a 4Sigma- state has the low-spin magnetic quantum number (s(z) = 1/2). When the state with s(z) = 3/2 is sought, the energy separations converge much more rapidly with the IP-EOM-CCSD value (0.03 eV) already being close to the observed (0.04 eV).     

Metal-modified silica adsorbents for removal of humic substances in water

As novel adsorbents for humic substances, Fe-, Mg-, and Ca-modified silica gels SiO(2)Fe, SiO(2)Mg, and SiO(2)Ca were prepared, and their adsorbabilities to humic and fulvic acids were evaluated in water at 25 degrees C for 20 h. Among these adsorbents, SiO(2)Fe indicated the highest adsorbability, in which removing humic substances in water was accomplished to 80-97%. By calcination at 600 degrees C and washing with water, adsorbabilities of the silica adsorbents deteriorated except for the case of calcination of SiO(2)Ca. This is due to changing metal modification structures and losing chlorine, judging from elemental analysis, TG-DTA, and XPS. Especially, the modifier Fe was tightly fixed on the silica surface of SiO(2)Fe, since the metal content was almost constant even after the calcination and water washing unlike the other adsorbents. Therefore, we found that SiO(2)Fe is the most useful adsorbent among the silica adsorbents and that its modification structure is composed of SiOFe, FeCl, and FeOH. In addition, the adsorption mechanism is explained by an interaction between Fe and humic substance molecule having carboxylate and phenolate groups, accompanied with anion exchange of chloride ion Cl(-) for the humic substance molecule via the silica pores. SiO(2)Fe may be applied to an adsorbent alternative for charcoal in water purification plants, and the used SiO(2)Fe may be further reused as a fertilizer since humic substances have plant-glowing ability and silica strengthens plant parts.     

Regional self-interaction correction of density functional theory

We propose a new simple scheme for self-interaction correction (SIC) of exchange functionals in the density functional theory. In the new scheme, exchange energies are corrected by substituting exchange self-interactions for exchange functionals in regions of self-interaction. To classify the regions of self-interaction, we take advantage of the property of the total kinetic energy density approaching the Weizs?cker density in the case of electrons in isolated orbitals. The scheme differs from conventional SIC methods in that it produces optimized molecular structures. Applying the scheme to the calculation of reaction energy barriers showed that it provides a clear improvement in cases where the barriers are underestimated by conventional "pure" functionals. In particular, we found that this scheme even reproduces a transition state that is not given by pure functionals.     

Dynamics of coherent anharmonic phonons in bismuth using high density photoexcitation

We have investigated the dynamical properties of the coherent anharmonic phonons generated in Bi under high density excitation. The time-resolved reflectivity in the intensely photoexcited Bi film is modulated by the coherent A(1g) phonon oscillation with a time-dependent oscillation period. As the pump power density is increased, the line shape of the A(1g) mode in the Fourier transformed spectra becomes asymmetric, and the redshift of the phonon frequency is observed. Analysis of the transient redshift with a wavelet transform reveals that the frequency of the A(1g) mode depends on the squared amplitude of the oscillation, which is attributed to an anharmonicity of the lattice potential.