Search

In this paper, we present a new approach to treat the electronic structure of a molecule in solution. Unlike the hybrid-type method, such as the reference interaction site model self-consistent-field theory, the new approach describes not only the electronic structure of solute but also solute–solvent interactions in terms of the quantum chemistry based on the Hartree–Fock frozen density formulation. In the treatment, the quantum effect due to solvent, including exchange repulsion, is projected on to the solute Hamiltonian using the spectral representation method. The solvent distribution around the solute is handled by the integral equation theory of liquids. As illustrative applications of the approach, the electronic and solvation structure of noble atoms, neon and argon, in liquid neon are studied. We also investigate the electronic structure of an excess electron in liquid helium. The preliminary results demonstrate that the quantum-mechanical effect on the electronic and solvation structure of the solute due to solvent molecules is successfully represented by the new method.In this paper, we present a new approach to treat the electronic structure of a molecule in solution. Unlike the hybrid-type method, such as the reference interaction site model self-consistent-field theory, the new approach describes not only the electronic structure of solute but also solute–solvent interactions in terms of the quantum chemistry based on the Hartree–Fock frozen density formulation. In the treatment, the quantum effect due to solvent, including exchange repulsion, is projected on to the solute Hamiltonian using the spectral representation method. The solvent distribution around the solute is handled by the integral equation theory of liquids. As illustrative applications of the approach, the electronic and solvation structure of noble atoms, neon and argon, in liquid neon are studied. We also investigate the electronic structure of an excess electron in liquid helium. The preliminary results demonstrate that the quantum-mechanical effect on the electronic and solvation structure of the s...

In a previous paper of this series [Paper III: Nakatsuji, J. Chem. Phys. 105, 2465 (2001)], the author showed a high potentiality of the extended coupled cluster (ECC) method to calculate the exact wave function of the ground state. In this paper, we propose ECC-configuration interaction (CI) method, which is an accurate useful method to calculate the excited states from the ECC wave function of the ground state. In contrast to the ECC method, the standard ECC-CI method is approximate, but we can make it exact by generalizing its excitation operator (ECC-CI general). The ECC-CI method is applicable not only to the excited states having the same spin-space symmetry as the ground state, but also to those having different spin-space symmetries and to the ionized and electron-attached states. The theoretical framework of the ECC-CI method is similar to that of the symmetry-adapted-cluster (SAC)-CI method proposed in 1978 by the present author. Next in this paper, we examine the performance of the methods proposed in this series of papers for a simple one-dimensional harmonic oscillator. The iterative configuration interaction (ICI) and ECC methods are examined for the ground state and the ICI-CI and ECC-CI methods for the excited states. The ICI method converges well to the exact ground state and the excited states are calculated nicely by the ICI-CI method in both the standard and general active spaces. In contrast to the simplest (S)ECC examined in Paper III, the ECC2 method shows quite a rapid convergence to the exact ground state, which enables us to calculate the true exact wave function in the ECC form. The ECC-CI methods in both the standard and general active spaces also work well to calculate the excited states. Thus, we conclude that the ICI and ECC approaches have a potentiality to provide useful method to calculate accurate wave functions of the ground and excited states. A merit of ECC is that it provides the exact wave function in a simple explicit form.

We continue to study exponential ansatz as a candidate of the structure of the exact wave function. We divide the Hamiltonian into ND (number of divisions) parts and extend the concept of the coupled cluster (CC) theory such that the cluster operator is made of the divided Hamiltonian. This is called extended coupled cluster (ECC) including ND variables (ECCND). It is shown that the S(simplest)ECC, including only one variable (ND=1), is exact in the sense that it gives an explicit solution of the Schrodinger equation when its single variable is optimized by the variational or H-nijou method. This fact further implies that the ECCND wave function with ND⩾2 should also have a freedom of the exact wave function. Therefore, by applying either the variational equation or the H-nijou equation, ECCND would give the exact wave function. Though these two methods give different expressions, the difference between them should vanish for the exact wave function. This fact solves the noncommuting problem raised in Paper I [H. Nakatsuji, J. Chem. Phys. 113, 2949 (2000)]. Further, ECCND may give more rapidly converging solution than SECC because of its non-linear character, ECCND may give the exact wave function at the sets of variables different from SECC. Thus, ECCND is exact not only for ND=1, but also for ND⩾2. The operator of the ECC, exp(S), is an explicit expression of the wave operator that transforms a reference function into the exact wave function. The coupled cluster including general singles and doubles (CCGSD) proposed in Paper I is an important special case of the ECCND. We have summarized the method of solution for the SECC and ECCND truncated at order n. The performance of SECC and ECC2 is examined for a simple example of harmonic oscillator and the convergence to the exact wave function is confirmed for both cases. Quite a rapid convergence of ECC2 encourages an application of the ECCND to more general realistic cases.

We report results from molecular dynamics simulations of the freezing transition of TIP5P water molecules confined between two parallel plates under the influence of a homogeneous external electric field, with magnitude of 5 V/nm, along the lateral direction. For water confined to a thickness of a trilayer we find two different phases of ice at a temperature of T=280 K. The transformation between the two, proton-ordered, ice phases is found to be a strong first-order transition. The low-density ice phase is built from hexagonal rings parallel to the confining walls and corresponds to the structure of cubic ice. The high-density ice phase has an in-plane rhombic symmetry of the oxygen atoms and larger distortion of hydrogen bond angles. The short-range order of the two ice phases is the same as the local structure of the two bilayer phases of liquid water found recently in the absence of an electric field [J. Chem. Phys. 119, 1694 (2003)]. These high- and low-density phases of water differ in local ordering at the level of the second shell of nearest neighbors. The results reported in this paper, show a close similarity between the local structure of the liquid phase and the short-range order of the corresponding solid phase. This similarity might be enhanced in water due to the deep attractive well characterizing hydrogen bond interactions. We also investigate the low-density ice phase confined to a thickness of 4, 5, and 8 molecular layers under the influence of an electric field at T=300 K. In general, we find that the degree of ordering decreases as the distance between the two confining walls increases. (C) 2004 American Institute of Physics.

In this paper we investigate the influence of self-affine roughness on the friction coefficient mu(f) of a rubber body under incomplete contact onto a solid surface. The roughness is characterized by the rms amplitude w, the correlation length xi, and the roughness exponent H. It is shown that with increasing surface roughening at short and/or long length scales (decreasing H and/or increasing ratio w/xi, respectively), the maximum of the friction coefficient mu(f) shifts to lower sliding velocities. The latter occurs only for conditions of incomplete contact for small contact length scales lambda (xi). In all cases, the friction coefficient mu(f) increases monotonically with decreasing roughness exponent H and/or increasing roughness ratio w/xi and attains its maximum value for sufficiently large contact length scales (xi).

Quasirelativistic theory of magnetic shielding constants based on the Douglas–Kroll–Hess transformation of the magnetic potential presented in a previous paper is extended to molecular systems that contain heavy elements. The gauge-including atomic orbital method is adapted to the quasirelativistic Hamiltonian to allow origin-independent calculations. The present theory is applied to the proton and halogen magnetic shielding constants of hydrogen halides and the 199Hg magnetic shielding constants and chemical shifts of mercury dihalides and methyl mercury halides. While the relativistic correction to the magnetic interaction term has little effect on the proton magnetic shielding constants, this correction is a dominant origin of the heavy atom shifts of the magnetic shielding constants of heavy halogens and mercury. The basis set-dependence of mercury shielding constants is quite large in the relativistic calculation; it is important to use the basis functions that are optimized by the relativistic metho...

In this paper a combined experimental and quantum chemical study of the geometry and opto-electronic properties of unsubstituted and dialkoxy-sustituted phenylene–vinylene oligomers (PV’s) is presented. The optical absorption spectra for PV cations with different chain lengths and substitution patterns were measured using pulse radiolysis with time-resolved spectrophotometric detection from 1380 to 500 nm (0.9 to 2.5 eV). The geometries of the PV’s studied were optimized using density functional theory (DFT) for both the neutral and singly charged molecule. The spectra for the PV radical cations were then calculated using singly excited configuration interaction with an intermediate neglect of differential overlap reference wave function method together with the DFT geometry. The agreement between experimental and theoretical absorption energies is excellent; most of the calculated radical cation absorption energies are within 0.15 eV of the experimental values. The pattern of dialkoxy-substitution is found to have a large effect on the optical absorption spectrum of the cation. Using the calculated charge distribution it is shown that the degree of delocalization of the charge correlates with the energy of the lowest absorption band. If alkoxy side chains are present on some of the rings the positive charge tends to localize at those sites.

In molecular J-aggregates one often observes an increase of the fluorescence decay time when increasing the temperature from 0 K. This phenomenon is usually attributed to the thermal population of the dark Frenkel exciton states that lie above the superradiant bottom state of the exciton band. In this paper, we study this effect for a homogeneous one-dimensional aggregate in a host medium and we model the scattering between different exciton states as arising from their coupling to the host vibrations. A Pauli master equation is used to describe the redistribution of excitons over the band. The rates entering this equation are calculated within the framework of first-order perturbation theory, assuming a linear on-site interaction between excitons and acoustic phonons. Solving the master equation numerically for aggregates of up to 100 molecules, we calculate the temperature dependence of the fluorescence kinetics in general and the decay time scale in particular. The proper definition of the fluorescence decay time is discussed in detail. We demonstrate that, even at a quantum yield of unity, the possibility to directly interpret fluorescence experiments in terms of a simple radiative time scale depends crucially on the initial excitation conditions in combination with the competition between spontaneous emission and intraband phonon-assisted relaxation.

In this paper we present the implementation of the two-component scaled zeroth-order regular approximation (ZORA) method in the molecular electronic structure package GAMESS-UK. It is the first application of this method, which was earlier investigated in the context of density functional theory, in molecular ab initio basis set calculations. The performance of the method is tested in atomic calculations, which we can compare with numerical results, on xenon and radon and in molecular calculations on the molecules AgH, HI, I2, AuH, TlH, and Bi2. In calculations on the I2 molecule we investigated the effect of the different approaches regarding the internal Coulomb matrix used in the ZORA method. For the remaining molecules we computed harmonic frequencies and bond lengths. It is shown that the scaled ZORA approach is a cost-effective alternative to the Dirac–Fock method.

In this paper a detailed theoretical analysis is presented of the temperature dependent radiative decay in aggregates of pseudoisocyanine (PIC). Our approach extends the original linear exciton-phonon coupling model including static disorder and second order exciton-phonon interactions. It is shown that for a one-dimensional exciton model neither of these additional effects alone or in combination with linear electron–phonon coupling can explain the steep rise in radiative lifetime at 40 K observed in the J-aggregate of PIC. However, when the aggregate assembles into a two-dimensional bricklike structure its radiative dynamics can be simulated, with linear exciton-optical phonon coupling as the only source for exciton scattering. Exciton-phonon scattering transfers oscillator strength from the k=0 state to other band states and also generates a nonequilibrium population among the exciton states, which persists during the superradiant decay. These effects together explain the marked temperature dependence of the radiative lifetime of the PIC J aggregate. When disorder limits the coherence length at low temperatures to a few molecules, as seems the case in several light harvesting complexes, the exciton population can equilibrate on the time scale of the superradiance. This situation pertains to the strong collision limit of the master equation, where the radiative decay is insensitive to details of the electron–phonon coupling, but only senses change in the thermal population among the exciton states.

The visible absorption spectrum of NO2 is very dense and irregular, and shows signs of a chaotic frequency and intensity distribution in the higher energy region. The complexity of the spectrum is related to a conical intersection of the potential energy surfaces of the two lowest electronic states. Above the conical intersection strong vibronic interactions lead to hybrid eigenstates, which can be viewed as mixtures of low vibrational levels of the electronically excited state and high vibrational levels of the electronic ground state. As a contribution to the elucidation of the nature of the vibronic bands of NO2 we have measured high-resolution spectra of a number of vibronic bands in the region between 10 000 and 14 000 cm−1 by exciting a supersonically cooled beam of NO2 molecules with a narrow-band Ti:Sapphire ring laser. The energy absorbed by the molecules was detected by a bolometer, and in some cases, laser-induced fluorescence was detected. The hyperfine structure is dominated by the Fermi-contact interaction and the magnitude of this interaction is a direct measure of the (electronic) composition of the hybrid eigenstates. In the present paper we have restricted our analysis to transitions of K⎯=0 stacks. The fine- and hyperfine structure of each rotational transition can be analyzed by using an effective Hamiltonian approach. The very good agreement that is found between the calculated transition strengths and the measured line intensities is evidence that in the spectral region studied, rovibronic interactions play a minor role. The composition of the hybrid eigenstates is compared with ab initio calculations reported in the literature, leading to the conclusion that measurements of the hyperfine structure are a helpful tool in characterizing vibronic bands.

In the present paper we consider nonlinear elasticity and friction of grafted persistent chains, which are highly stretched in the normal to the surface direction due to orientational interactions. We examine the normal and the lateral forces both in equilibrium and under shear sliding when the monolayer is confined by a bare surface. We show that in the confined monolayer in equilibrium the tilted orientation of the director becomes stable. In the sliding regime the friction force passes through a maximum value. The additional normal force in the sliding regime, when the distance between the surfaces is fixed, is also considered. We show that this force is attractive for small velocities and becomes repulsive for high velocities after the friction force passes through the maximum value.

The flux-flux correlation function formalism is a standard and widely used approach for the computation of reaction rates. In this paper we introduce a method to compute the classical and quantum flux-flux correlation functions for anharmonic barriers essentially analytically through the use of the classical and quantum normal forms. In the quantum case we show that the quantum normal form reduces the computation of the flux-flux correlation function to that of an effective one dimensional anharmonic barrier. The example of the computation of the quantum flux-flux correlation function for a fourth order anharmonic barrier is worked out in detail, and we present an analytical expression for the quantum mechanical microcanonical flux-flux correlation function. We then give a discussion of the short-time and harmonic limits., 9 pages, 2 figures

Low-frequency intermolecular dynamics in liquids is studied by ultrafast four- and six-wave mixing. The theory of these nonlinear optical processes is given for electronically nonresonant optical interactions up to fifth order in the electric field. The Born-Oppenheimer approximation is used to separate the motional part of the response functions from coordinate independent electronic hyperpolarizabilities. A large variety of experiments, involving far-infrared absorption, ordinary Rayleigh-Raman or hyper Rayleigh-Raman scattering is covered by this theory. The response in nonresonant six-wave mixing comprises four dynamically different processes. It is shown that one of the terms contains information on the time scale(s) of intermolecular dynamics, that is not available from lower-order nonresonant experiments. For instance, homogenous and inhomogeneous contributions to line broadening can be distinguished. The optical response of harmonic nuclear motion is calculated for nonlinear coordinate dependence of the polarizabilities. Results for level-dependent and level-independent damping of the motion are compared. It is shown that level-dependent damping destroys the interference between different quantum mechanical pathways, yielding an extra contribution to the fifth-order response that has not been discussed before. When two or more nuclear modes determine the optical response, their relative contributions to the four- and six-wave mixing signals are in general different. These contributions are determined by the coordinate dependence of the electronic polarizability, which is usually not fully known, Model calculations are presented for the dynamic parameters of liquid CS2. The theory of this paper will be employed in Part II, to analyze experimental results on femtosecond four- and six-wave mixing.

In quantum-mechanical/molecular-mechanical (QM/MM) treatment of chemical reactions in condensed phases, one solves the electronic Schrodinger equation for the solute (or an active site) under the electrostatic field from the environment. This Schrodinger equation depends parametrically on the solute nuclear coordinates R and the external electrostatic potential V. This fact suggests that one may use R and V as natural collective coordinates for describing the entire system, where V plays the role of collective solvent variables. In this paper such an (R,V) representation of the QM/MM canonical ensemble is described, with particular focus on how to treat charge transfer processes in this representation. As an example, the above method is applied to the proton-coupled electron transfer of a ubiquinol analog with phenoxyl radical in acetonitrile solvent. Ab initio free-energy surfaces are calculated as functions of R and V using the reference interaction site model self-consistent field method, the equilibrium points and the minimum free-energy crossing point are located in the (R,V) space, and then the kinetic isotope effects (KIEs) are evaluated approximately. The results suggest that a stiffer proton potential at the transition state may be responsible for unusual KIEs observed experimentally for related systems.

In this paper we present a comprehensive lattice Monte Carlo study of long chain conformations at solid-polymer melt interfaces. Segmental scale interfacial features, like the bond orientational distribution were found to be independent of surface-segment energetics, and statistically identical with Helfand’s predictions for the full-occupancy, infinite chain length limit. Conformational statistics of chains longer than 5-6 statistical segments followed the predictions of the Scheutjens-Fleer theory and the same power laws as a single ideal chain at the critical value of the surface-segment adsorption free-energy. Our simulations tested the predictions of random walk next to a “reflective” surface statistics for the spatial variations of chain dimensions and chain center of mass density. It was found that these statistics furnish the correct long chain limit, independently of surface-segment energetics. The random walk next to a “reflective” boundary predictions for the “adsorbed” amount and the distributions of tail, loop, and train number, as well as tail and loop size were in quantitative agreement with the simulation data. The correspondence between random walks and real (or simulated) chains required the knowledge of a single, microscopic parameter, the number of chemical segments per statistical segment, a∞. This quantity was very close to the average length of adsorbed sequences (trains), in the long chain limit. Our simulations tested thoroughly and established firmly the validity of “reflective” boundary statistics in the melt. The inevitability of these statistics has broad implications on “desorption” kinetics, chain mobility, and chain relaxation, which are currently under study.

The fifth-order two-dimensional (2D) Raman signals have been calculated from the equilibrium and nonequilibrium (finite field) molecular dynamics simulations. The equilibrium method evaluates response functions with equilibrium trajectories, while the nonequilibrium method calculates a molecular polarizability from nonequilibrium trajectories for different pulse configurations and sequences. In this paper, we introduce an efficient algorithm which hybridizes the existing two methods to avoid the time-consuming calculations of the stability matrices which are inherent in the equilibrium method. Using nonequilibrium trajectories for a single laser excitation, we are able to dramatically simplify the sampling process. With this approach, the 2D Raman signals for liquid xenon, carbon disulfide, water, acetonitrile, and formamide are calculated and discussed. Intensities of 2D Raman signals are also estimated and the peak strength of formamide is found to be only five times smaller than that of carbon disulfide.

The ground state structure for the CCCCCl radical was computed by using symmetry-adapted cluster configuration-interaction (SAC-CI) theory along with density functional theory to overcome the differences raised in the recently published paper [Y. Sumiyoshi et al., Chem. Phys. Lett. 414, 82 (2005)] between the theory and the experiment. SAC-CI results clearly support the earlier experimental conclusion that the radical has the bent ground state structure corresponding to 2Pi symmetry. Contrarily, probably due to spin contamination, mixing of a bent doublet ground state with the quartet components of a linear structure, coupled-cluster singles and doubles (CCSD) calculations were unable to provide reliable results. Results obtained using density functional theory also show that the radical has a bent structure. Some low-lying doublet excited states were also studied using the SAC-CI theory. The energy difference between the ground Pi state and the nearby Sigma state is around 0.2 eV. The excitation energy for the transition with the largest oscillator strength agrees with the strongest absorption peak.

In this paper we examine the effect of radiative and radiationless processes on optical free induction decay (OFID) in singlet–singlet transitions of molecules in the condensed phase. The main conclusion is that the applicability of OFID to molecular singlet–singlet transitions is severely limited by the presence of an intermediate triplet state, which quite often forms a bottleneck in the optical pumping cycle. Using the theoretically derived formulas for the OFID signal, we examine the feasibility of OFID experiments on the mixed crystals of porphin in n‐octane and naphthalene in durene and conclude that such experiments, if at all fruitful, would be difficult to execute. We re‐evaluate the reported OFID results of pentacene in p‐terphenyl and obtain for the lowest 1B2u←1A1g transition a dipole moment of 0.9±0.3 D.

In this paper the recently reported fluorescence transient, appearing during the onset of cw dye laser excitation of the energetically lower sites of pentacene in p-terphenyl at 2 K, is analyzed quantitatively. The role of the lowest triplet state in causing the long (compared to the fluorescence lifetime) transient, is firmly established. For the lower sites a kinetic analysis yields a T1←S1 intersystem crossing (ISC) yield of about 0.4% and an S0 ←T1 ISC rate of (2.2 ± 0.1) × 10^4 s–1. Using these data in calculating the S1 ←S0 transition dipole moment from our earlier OFID experiments yields for the O1 site (proto) μ=0.71 ±0.24 D. The S0 ←S1 internal conversion yield is estimated to be 22%. The fluorescence transients of the energetically higher sites of pentacene in p-terphenyl exhibit a much faster decay, which indicates a drastic increase in ISC yield from S1. Optical free induction decay (OFID) experiments performed on these sites confirm this and imply that the T1 ←S1 ISC yield for these sites is ~60%. We further discuss results of an OFID experiment on pentecene in naphthalene and derive in this case at T1←S1 ISC yield of ≤2.8% and an S0←S1 internal conversion yield of 31.5 ± 1.5%. The off-resonance pumping effect in a 3-level system was also studied, using the parameter values for pentacene in p-terphenyl. For the lower sites, with the excitation intensities used (~1 MW/m2) the part of the inhomogeneous absorption line that is pumped has a width of about 2.5 GHz (effective excitation bandwidth

Pyrazine is shown by its Molecular Eigenstate Spectra and by its quantum beats to be a small molecule. The effect of Coriolis coupling and non-resonant light scattering are discussed. This paper further treats the so-called "fast component". Its existence with a decay time of about 100 psec has by now been firmly established. It finds its cause in the "bunching" of the molecular eigenstates in frequency space. For broad (i.e. psec) lasers it is therefore still useful to speak of pyrazine as an ILS-molecule. The characterization depends on the laser one uses in the experiment.