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...

The complexity of the absorption spectrum of NO2NO2 can be attributed to a conical intersection of the potential energy surfaces of the two lowest electronic states, the electronic ground state of 2A12A1 symmetry and the first electronically excited state of 2B22B2 symmetry. In a previous paper we reported on the feasibility of using the hyperfine splittings, specifically the Fermi-contact interaction, to determine the electronic ground state character of the excited vibronic states in the region just above the conical intersection; 10 000 to 14 000 cm−114 000 cm−1 above the electronic ground state. High-resolution spectra of a number of vibronic bands in this region were measured by exciting a supersonically cooled beam of NO2NO2 molecules with a narrow-band Ti:Sapphire ring laser. The energy absorbed by the molecules was detected by the use of a bolometer. In the region of interest, rovibronic interactions play no significant role, with the possible exception of the vibronic band at 12 658 cm−1,12 658 cm−1, so that the fine- and hyperfine structure of each rotational transition could be analyzed by using an effective Hamiltonian. In the previous paper we restricted ourselves to an analysis of transitions of the K⎯=0K−=0 stack. In the present paper we extend the analysis to transitions of the K⎯=1K−=1 stack, from which, in addition to hyperfine coupling constants, values of the AA rotational constants of the excited NO2NO2 molecules can be determined. Those rotational constants also contain information about the electronic composition of the vibronic states, and, moreover, about the geometry of the NO2NO2 molecule in the excited state of interest. The results of our analyses are compared with those obtained by other authors. The conclusion arrived at in our previous paper that determining Fermi-constants is useful to help characterize the vibronic bands, is corroborated. In addition, the AA rotational constants correspond to geometries that are consistent with the electronic composition of the relevant excited states as expected from the Fermi-constants., The complexity of the absorption spectrum of NO2 can be attributed to a conical intersection of the potential energy surfaces of the two lowest electronic states, the electronic ground state of 2A1 symmetry and the first electronically excited state of 2B2 symmetry. In a previous paper we reported on the feasibility of using the hyperfine splittings, specifically the Fermi-contact interaction, to determine the electronic ground state character of the excited vibronic states in the region just above the conical intersection; 10 000 to 14 000 cm−1 above the electronic ground state. High-resolution spectra of a number of vibronic bands in this region were measured 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 the use of a bolometer. In the region of interest, rovibronic interactions play no significant role, with the possible exception of the vibronic band at 12 658 cm−1, so that the fine- and hyperfine structure of each rotational transition could be analyzed by using an effective Hamiltonian. In the previous paper we restricted ourselves to an analysis of transitions of the K⎯=0 stack. In the present paper we extend the analysis to transitions of the K⎯=1 stack, from which, in addition to hyperfine coupling constants, values of the A rotational constants of the excited NO2 molecules can be determined. Those rotational constants also contain information about the electronic composition of the vibronic states, and, moreover, about the geometry of the NO2 molecule in the excited state of interest. The results of our analyses are compared with those obtained by other authors. The conclusion arrived at in our previous paper that determining Fermi-constants is useful to help characterize the vibronic bands, is corroborated. In addition, the A rotational constants correspond to geometries that are consistent with the electronic composition of the relevant excited states as expected from the Fermi-constants.

Nucleation rate data, obtained from expansion wave tube experiments, are reported for several vapor–gas mixtures at high pressure. Results are given for water–vapor in the presence of helium and nitrogen gas, and for n-nonane in helium and methane. For all these mixtures, carrier gas pressures of 10, 25, and 40 bar have been applied, with temperatures ranging from 230 to 250 K. An extended form of the nucleation theorem (in terms of the derivative of the nucleation rate with respect to carrier gas pressure) is derived, which appears to be very helpful in the interpretation of high pressure data. It can be used to obtain the carrier gas content of the critical nucleus directly from the pressure dependence of experimental nucleation rates. Combining this method with the theoretical considerations of part I of this paper [J. Chem. Phys. 111, 8524 (1999), preceding paper]: the nucleation behavior of water at high pressures of both helium and nitrogen can quantitatively be understood. For n-nonane in helium our "pressure perturbation approach" is also valid. For n-nonane in methane, however, this approach fails because of the high methane solubility in the liquid phase.

The speed and angular distribution of O atoms arising from the photofragmentation of C(5)H(8)-O(2), the isoprene-oxygen van der Waals complex, in the wavelength region of 213-277 nm has been studied with the use of a two-color dissociation-probe method and the velocity map imaging technique. Dramatic enhancement in the O atoms photo-generation cross section in comparison with the photodissociation of individual O(2) molecules has been observed. Velocity map images of these "enhanced" O atoms consisted of five channels, different in their kinetic energy, angular distribution, and wavelength dependence. Three channels are deduced to be due to the one-quantum excitation of the C(5)H(8)-O(2) complex into the perturbed Herzberg III state ((3)Δ(u)) of O(2). This excitation results in the prompt dissociation of the complex giving rise to products C(5)H(8)+O+O when the energy of exciting quantum is higher than the complex photodissociation threshold, which is found to be 41740 ± 200 cm(-1) (239.6±1.2 nm). This last threshold corresponds to the photodissociation giving rise to an unexcited isoprene molecule. The second channel, with threshold shifted to the blue by 1480 ± 280 cm(-1), corresponds to dissociation with formation of rovibrationally excited isoprene. A third channel was observed at wavelengths up to 243 nm with excitation below the upper photodissociation threshold. This channel is attributed to dissociation with the formation of a bound O atom C(5)H(8)-O(2) + hv → C(5)H(8)-O(2)((3)Δ(u)) → C(5)H(8)O + O and/or to dissociation of O(2) with borrowing of the lacking energy from incompletely cooled complex internal degrees of freedom C(5)H(8)*-O(2) + hv → C(5)H(8)*-O(2)((3)Δ(u)) → C(5)H(8) + O + O. The kinetic energy of the O atoms arising in two other observed channels corresponds to O atoms produced by photodissociation of molecular oxygen in the excited a (1)Δ(g) and b (1)Σ(g)(+) singlet states as the precursors. This indicates the formation of singlet oxygen O(2)(a (1)Δ(g)) and O(2)(b (1)Σ(g)(+)) after excitation of the C(5)H(8)-O(2) complex. Cooperative excitation of the complex with a simultaneous change of the spin of both partners (1)X-(3)O(2) + hν → (3)X-(1)O(2) → (3)X + (1)O(2) is suggested as a source of singlet oxygen O(2)(a (1)Δ(g)) and O(2)(b (1)Σ(g)(+)). This cooperative excitation is in agreement with little or no vibrational excitation of O(2)(a (1)Δ(g)), produced from the C(5)H(8)-O(2) complex as studied in the current paper as well as from the C(3)H(6)-O(2) and CH(3)I-O(2) complexes reported in our previous paper [Baklanov et al., J. Chem. Phys. 126, 124316 (2007)]. The formation of O(2)(a (1)Δ(g)) from C(5)H(8)-O(2) was observed at λ(pump) = 213-277 nm with the yield going down towards the long wavelength edge of this interval. This spectral profile is interpreted as the red-side wing of the band of a cooperative transition (1)X-(3)O(2) + hν → (3)X(T(2))-(1)O(2)(a (1)Δ(g)) in the C(5)H(8)-O(2) complex.

We investigate the vibrational shift of beryllium oxide (BeO) in Xe matrix as well as in Ar matrix environments by mixed quantum-classical simulation and examine the origin of spectral shift in details. BeO is known to form strong chemical complex with single rare gas atom, and it is predicted from the gas phase calculations that vibrational frequencies are blueshifted by 78 cm^[-1] and 80 cm^[-1] upon formation of XeBeO and ArBeO, respectively. When the effects of other surrounding rare gas atoms are included by Monte Carlo simulations, it is found that the vibrational frequencies are redshifted by 21 cm^[-1] and 8 cm^[-1] from the isolated XeBeO and ArBeO complexes, respectively. The vibrational shift of XeBeO in Ar matrix is also calculated and compared with experimental data. In all simulations examined in this paper, the calculated vibrational frequency shifts from the isolated BeO molecule are in reasonable agreement with experimental values. The spectral shift due to the rare-gas-complex formation of RgBeO (Rg = Xe or Ar) is not negligible as seen in the previous studies, but it is shown in this paper that the effects of other surrounding rare gas atoms should be carefully taken into account for quantitative description of the spectral shifts and that these two effects are competing in vibrational spectroscopy of BeO in matrix environments.

Current efforts to characterize and study interstellar polycyclic aromatic hydrocarbons(PAHs) rely heavily on theoretically predicted infrared (IR) spectra. Generally, such studies use the scaled harmonic frequencies for band positions and double harmonic approximation for intensities of species, and then compare these calculated spectra with experimental spectra obtained under matrix isolation conditions. High-resolution gas-phase experimental spectroscopic studies have recently revealed that the double harmonic approximation is not sufficient for reliable spectra prediction. In this paper, we present the anharmonic theoretical spectra of three PAHs: naphthalene, anthracene, and tetracene, computed with a locally modified version of the SPECTRO program using Cartesian derivatives transformed from Gaussian 09 normal coordinate force constants. Proper treatments of Fermi resonances lead to an impressive improvement on the agreement between the observed and theoretical spectra, especially in the C–H stretching region. All major IR absorption features in the full-scale matrix-isolated spectra, the high-temperature gas-phase spectra, and the most recent high-resolution gas-phase spectra obtained under supersonically cooled molecular beam conditions in the CH-stretching region are assigned.

This paper presents a detailed study of a water adlayer adsorbed on Pt(111) and Rh(111) surfaces using periodic density functional theory methods. The interaction between the metal surface and the water molecules is assessed from molecular dynamics simulation data and single point electronic structure calculations of selected configurations. It is argued that the electron bands around the Fermi level of the metal substrate extend over the water adlayer. As a consequence in the presence of the water layer the surface as a whole still maintains its metallic conductivity-a result of a crucial importance for understanding the process of electron transfer through the water/metal interface and electrochemical reactions in particular. Our results also indicate that there exists a weak bond between the hydrogen of the water and the Rh metal atoms as opposed to the widespread (classical) models based on purely repulsive interaction. This suggests that the commonly used classical interactions potentials adopted for large scale molecular dynamics simulations of water/metal interfaces may need revision. Two adsorption models of water on transition metals with the OH bonds pointing towards or away of the surface are also examined. It is shown that due to the very close values of their adsorption energies one should consider the real structure of water on the surface as a mixture of these simple "up" and "down" models. A model for the structure of the adsorbed water layer on Rh(111) is proposed in terms of statistical averages from molecular dynamics simulations.

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.

This paper discusses a time-dependent density functional theory study of the effect of molecular structure on the excited state polarizability of conjugated molecules. A short phenylenevinylene oligomer containing three phenyl rings (PV2, distyryl benzene) is taken as a model system. Introduction of methyl substituents is shown to have only a small influence on the increase in polarizability upon excitation (the excess polarizability, Delta(alpha) over bar). Methoxy groups have a much larger effect but in this case Delta(alpha) over bar depends strongly on the dihedral angle between the side chain and the backbone of the molecule. If the central phenyl ring of PV2 has a meta-configuration rather than para, both the optical absorption spectrum and the excess polarizability change considerably. (C) 2003 American Institute of Physics.

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 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.

The key characteristics of electronic excitations of many-electron systems, the excitation energies ωα and the oscillator strengths fα, can be obtained from linear response theory. In one-electron models and within the adiabatic approximation, the zeros of the inverse response matrix, which occur at the excitation energies, can be obtained from a simple diagonalization. Particular cases are the eigenvalue equations of time-dependent density functional theory (TDDFT), time-dependent density matrix functional theory, and the recently developed phase-including natural orbital (PINO) functional theory. In this paper, an expression for the oscillator strengths fα of the electronic excitations is derived within adiabatic response PINO theory. The fα are expressed through the eigenvectors of the PINO inverse response matrix and the dipole integrals. They are calculated with the phaseincluding natural orbital functional for two-electron systems adapted from the work of Lowdin ¨ and Shull on two-electron systems (the phase-including Löwdin-Shull functional). The PINO calculations reproduce the reference fα values for all considered excitations and bond distances R of the prototype molecules H2 and HeH+ very well (perfectly, if the correct choice of the phases in the functional is made). Remarkably, the quality is still very good when the response matrices are severely restricted to almost TDDFT size, i.e., involving in addition to the occupied-virtual orbital pairs just (HOMO+1)-virtual pairs (R1) and possibly (HOMO+2)-virtual pairs (R2). The shape of the curves fα(R) is rationalized with a decomposition analysis of the transition dipole moments. peerReviewed

The discrete Kohn variational reactive scattering method presented in the preceding paper is applied to the reaction H+O2→ OH+O. The essential features of the method are the use of an interaction region grid, fully coupled open and closed distorted waves (DWs) and the use of a singular value decomposition technique to construct the optimal boundary conditions for the open DWs. A convergence test is presented at a total energy of 0.817 eV above the bottom of the H+O2 well. It is shown that very well converged results may be obtained in a calculation with a relatively small interaction region grid, when at least a few asymptotically closed DWs are included in the trial wave function. Furthermore, the number of open distorted waves (DWs) may be considerably smaller than the number of open channels. Six additional points are computed in an energy range of 1.2 meV, scanning through a narrow resonance around 0.817 eV. The results are in very good agreement with the hyperspherical coordinate propagation calculat...

In this paper we present time-dependent density functional calculations on frequency-dependent first (β) and second (γ) hyperpolarizabilities for the set of small molecules, N2, CO2, CS2, C2H4, NH3, CO, HF, H2O, and CH4, and compare them to Hartree–Fock and correlated ab initio calculations, as well as to experimental results. Both the static hyperpolarizabilities and the frequency dispersion are studied. Three approximations to the exchange-correlation (xc) potential are used: the widely used Local Density Approximation (LDA), the Becke–Lee–Yang–Parr (BLYP) Generalized Gradient Approximation (GGA), as well as the asymptotically correct Van Leeuwen–Baerends (LB94) potential. For the functional derivatives of the xc potential the Adiabatic Local Density Approximation (ALDA) is used. We have attempted to estimate the intrinsic quality of these methods by using large basis sets, augmented with several diffuse functions, yielding good agreement with recent numerical static LDA results. Contrary to claims which have appeared in the literature on the basis of smaller studies involving basis sets of lesser quality, we find that the static LDA results for β and γ are severely overestimated, and do not improve upon the (underestimated) Hartree–Fock results. No improvement is provided by the BLYP potential which suffers from the same incorrect asymptotic behavior as the LDA potential. The results are however clearly improved upon by the LB94 potential, which leads to underestimated results, slightly improving the Hartree–Fock results. The LDA and BLYP potentials overestimate the frequency dependence as well, which is once again improved by the LB94 potential. Future improvements are expected to come from improved models for asymptotically correct exchange-correlation potentials. Apart from the LB94 potential used in this work, several other asymptotically correct potentials have recently been suggested in the literature and can also be expected to improve considerably upon the relatively poor LDA and GGA results, for both the static properties and their frequency dependence.

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.

In this paper, we study the evolution of phase-separating binary mixtures which are subjected to alternate cooling and heating cycles. An initially homogeneous mixture is rapidly quenched to a temperature T(1)T(c), where T(c) is the critical temperature. The mixture undergoes phase separation for a while and is then suddenly heated to a temperature T(2)T(c). These cycles are repeated to create a domain morphology with multiple length scales, i.e., the structure factor is characterized by multiple peaks. For phase separation in d = 2 systems, we present numerical and analytical results for the emergence and growth of this multiple-scale morphology.

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.

The exchange-correlation potentials v xc which are currently fashionable in density functional theory ~DFT!, such as those obtained from the local density approximation ~LDA! or generalized gradient approximations ~GGAs!, all suffer from incorrect asymptotic behavior. In atomic calculations, this leads to substantial overestimations of both the static polarizability and the frequency dependence of this property. In the present paper, it is shown that the errors in atomic static dipole and quadrupole polarizabilities are reduced by almost an order of magnitude, if a recently proposed model potential with correct Coulombic long-range behavior is used. The frequency dependence is improved similarly. The model potential also removes the overestimation in molecular polarizabilities, leading to slight improvements for average molecular polarizabilities and their frequency dependence. For the polarizability anisotropy we find that the model potential results do not improve over the LDA and GGA results. Our method for calculating frequency-dependent molecular response properties within time-dependent DFT, which we described in more detail elsewhere, is summarized. © 1996 American Institute of Physics.@S0021-9606~96!02031-4#

In this paper, an extensive analysis of free volume related properties of a lipid membrane is given. Using percolation theory, and comparing the free volume properties to those of a soft polymer, additional insights are obtained. The analyses are discussed within the framework of the four region model. It is concluded that the four regions have very different free volume properties. The region containing the dense part of the lipid tails resembles a soft polymer membrane to a large extent. The middle part of the membrane is more similar to a low density alkane. The consequences of the computed free volume properties on the permeation process of small penetrants are discussed. (C) 1996 American Institute of physics.

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 describe a numerical method for the calculation of collective diffusion relaxation mechanisms in quenched block copolymer melts. The method entails the repeated calculation of two opposing fields-an external potential field U, conjugate to the density field rho, and an energetic interaction field E. The external field is calculated by numerical inversion of the density functionals and the energetic interaction field is calculated directly by integration over the density field. When the two fields are balanced U = E, we recover the self-consistent field solutions; when the two fields are off balance, the spatial gradient of E-U is the thermodynamic force which drives the collective diffusion. We introduce a simple local coupling approximation for the Onsager kinetic coefficients of short freely jointed chains in weakly ordered systems. Fluctuations are added by incorporation of a random Langevin force in the diffusion equation. Numerical results of decomposition in symmetric and asymmetric diblock copolymer melts indicate that the method is capable of describing extremely slow defect annihilation relaxation modes. We find that in the nonlinear regime, the density patterns evolve to metastable states, in which isolated defects separate relatively well-ordered crystalline microdomains. These final states are typical for many industrial applications of incompletely relaxed copolymer melts.

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.

In this paper, we report the first rotationally resolved one- and two-color resonance enhanced multiphoton ionization photoelectron spectroscopy (REMPI-PES) study of the HCl molecule. The agreement between our experimental branching ratios and theoretical investigations is excellent. We also report the first zero kinetic energy pulsed field ionization (ZEKE-PFI) experiments carried out in a "magnetic bottle'' electron spectrometer. A direct comparison is made between ZEKE-PFI and REMPI-PES spectra for ionization via several rotational levels of the F 12(v'=0) and f 32(v'=0) Rydberg states of HCl. Large differences in both the spin-orbit and rotational branching ratios are found between the ZEKE-PFI and REMPI-PES spectra. These differences can be understood qualitatively on the basis of rotational and spin-orbit autoionization mechanisms. The Journal of Chemical Physics is copyrighted by The American Institute of Physics.

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, a full account of an earlier report [Takahashi, Okuyama, and Kimura, J. Mol. Struct. 249, 47 (1991)] on the cation vibrational spectra of trans and gauche n‐propylbenzene, which were obtained by means of a two‐color (1+1’) multiphoton ionization threshold photoelectron technique is presented. The trans and gauche cations were separately produced by (1+1’) multiphoton ionization resonant through the S1 vibronic levels of n‐propylbenzene in supersonic free jets in which both the trans and gauche isomers exist. From the observed threshold photoelectron spectra, adiabatic ionization energies are determined as Ia(trans)=70 278 cm−1 (8.7134 eV) and Ia(gauche)=70 420 cm−1 (8.7311 eV) with an accuracy of ±8 cm−1 (±1 meV). Furthermore, four benzene ring modes (6b+, 1+, 12+, and 18a+) as well as several low‐frequency torsional and bending modes have been identified which are sensitive to the relative conformations of the n‐propyl group with respect to the benzene ring [(trans)+: 82, 212, and 300 cm−1...

This paper discusses the conformational phase change of an isolated polymer chain which is capable of forming physical bonds with solvent molecules. On the basis of the mean‐field theory of Flory type, we derive a formula for the temperature dependence of the expansion factor of the chain. Our theory takes account of the extra mixing entropy change induced by the solvation. We predict that the physical‐bond formation between polymer and solvent molecules causes a reentrant conformational change between coiled and globular state when temperature is varied. We also show that the polymer chain exhibits a sharp collapse near the Θ temperature as the temperature is increased to the lower critical solution temperature of the polymer solution. This behavior can be interpreted in terms of the breakup of the polymer–solvents complex. The result is compared with the observed coil–globule transition on poly (N‐isopropylacrylalmide) in water.

An examination of the finite system size effects on the calculated interfacial properties of the molten salt KI is described, as an extension of the work described in the preceding paper [J. Chem. Phys. 115, 8603 (2001)]. It is shown that, for system sizes beyond ∼1000 ions, the finite size effects are accounted for by the predictions of capillary wave (CW) theory, as in van der Waals liquids. Hence, there do not appear to be additional long-range effects arising from the Coulombic interactions. Techniques used to avoid truncation effects in the Coulombic and dispersion interactions appear to have eliminated system size dependent artefacts from these sources. The scaling behavior predicted from CW theory is used to predict the bulk surface tension for the chosen interaction potential, which is then compared with experimental values, and very good agreement found. Examination of the interfacial layers show that the local structure remains strongly charge ordered, but the surface is highly roughened by clus...

The hydrogen storage system LiH + NH(3) ↔ LiNH(2) + H(2) is one of the most promising hydrogen storage systems, where the reaction yield can be increased by replacing Li in LiH with other alkali metals (Na or K) in order of LiNaK. In this paper, we have studied the alkali metal M (M = Li, Na, K) dependence of the reactivity of MH with NH(3) by calculating the potential barrier of the H(2) desorption process from the reaction of an M(2)H(2) cluster with an NH(3) molecule based on the ab initio structure optimization method. We have shown that the height of the potential barrier becomes lower in order of Li, Na, and K, where the difference of the potential barrier in Li and Na is relatively smaller than that in Na and K, and this tendency is consistent with the recent experimental results. We have also shown that the H-H distance of the H(2) dimer at the transition state takes larger distance and the change of the potential energy around the transition state becomes softer in order of Li, Na, and K. There are almost no M dependence in the charge of the H atom in NH(3) before the reaction, while that of the H atom in M(2)H(2) takes larger negative value in order of Li, Na, and K. We have also performed molecular dynamics simulations on the M(2)H(2)-NH(3) system and succeeded to reproduce the H(2) desorption from the reaction of Na(2)H(2) with NH(3).

A new effective Hamiltonian is presented for the analysis of the high-resolution spectra of open-shell van der Waals complexes containing the O2 molecule. The effects of electron spin are included but the complications of nuclear spin and resultant nuclear spin splitting are neglected. The Hamiltonian is composed of the rotational, centrifugal distortion, and spin-spin interaction terms. The resulting energy levels are divided into two well-separated groups and the pattern is a complicated function of θ (the angle that the O2 molecule makes with the principal a axis of the complex) and φ (the azimuthal angle of the O2 out of the plane defined by the a and b axes of the complex). This model has been successfully applied to analyze the high-resolution spectrum of O2-N2O in the region of the N2O monomer ν3 vibrational band, which will be presented in a separate paper. © 1997 American Institute of Physics.

Molecular dynamics simulation has been performed on a liquid film that is sheared in between solid surfaces. As a shear is given to the liquid film, a Couette-like flow is generated in the liquid and energy conversion occurs from the macroscopic flow to the thermal energy, which is discharged back to the solid walls. In such a way, momentum and thermal energy fluxes are present simultaneously. And all these thermal and fluid phenomena take place in highly nonequilibrium state where thermal energy is not distributed equally to each degree of freedom of molecular motion in the vicinities of the solid-liquid interface. In the present paper, platinum and water are employed as solid and liquid, respectively. First, the structure and orientation of water molecules in the vicinities of the solid surfaces are analyzed and how these structure and orientation are influenced by the shear is considered. Based on this result, momentum and thermal energy transfer in the vicinities of and at the solid-liquid interfaces are investigated in detail. Results are compared with those of our previous study, in which monatomic and diatomic molecules are employed as liquid.

The reactivity and spectra of formaldehyde isomers and dimeric complexes between them are studied with ab initio methods. A large number of complexes between H2CO, trans-HCOH, cis-HCOH is calculated. Infrared and Raman spectra of (H2CO)2 are calculated with relatively simple methods using spectroscopic masses and scaled force constants. In this way, the structure of dimers in matrices can be deduced. Hydroxycarbene (HCOH) plays a key role in a model that explains a large number of experimental facts of the nonzero pressure photochemistry. Hydroxycarbene forms complexes with H2CO; the stabilization is due to classical hydrogen bonds. HCOH is a new example of an ambiphilic carbene. Addition products are formed from HCOH···H2CO complexes. The calculations show that, in agreement with matrix experiments, glycoaldehyde and methanol are easily formed. The formation of trans-HCOH occurs through a dimeric interaction with the shifting hydrogen. This bimolecular process is 9.6 kcal/mol (6–31G*) in favor of the unimolecular conversion. cis-HCOH might be formed via a nonplanar transition state, where also stabilization at the carbenic center is possible. When higher concentrations of HCOH are available, a hydrogen exchange mechanism easily transfers hydroxycarbene back to H2CO. Several experiments are suggested in this paper; notably, molecular beam and isotopic-mixture experiments will give useful information. The involvement of HCOH in the light-induced formose reaction is suggested.

This paper reports an ab initio CI calculation of the radiationless decay of the formaldehyde 1A2 state. First we derive quantitative conditions, which a basis set must satisfy if it is to be used for describing radiationless decay. We checked these conditions for the formaldehyde molecule and found them to be satisfied for the adiabatic Born–Oppenheimer set used in the calculation. We then derive a general equation for the coupling elements resulting from this basis set. With the method used, rotational coupling could be treated completely equivalent with vibrational coupling; this rotational coupling turned out to be not important in formaldehyde however. The coupling elements for D2CO are a factor 10 smaller than the corresponding ones in H2CO. The results of the calculation show, that the internal conversion in formaldehyde is an example of the so-called resonance case. Therefore the decay cannot be described by the model proposed by Yeung and Moore, where S1¿S0 internal conversion is the rate determining step in the photodissociation. Finally, we discuss the applicability of the ’’golden rule’’ in describing radiationless decay.

This paper reports an ab initio CI calculation of the radiative 1A1¿1A2 transition of H2CO and D2CO. Throughout the calculation the electronic wavefunctions and transition moments are explicitly calculated as functions of the nuclear geometry, contrary to the conventional Herzberg–Teller approach. The evaluation of the vibrational wavefunctions and integrals was made numerically. The results show that the excited state frequency for mode 3 has to be reassigned and that the calculated vibrational structure agrees well with the experimental intensities.

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.

This paper presents the exact solution of a lattice statistical problem posed by Makita and Miki in their theory of the anomalous specific heat of the ferroelectric sodium trihydrogen selenite. We conclude with a comparison of the experimental and calculated entropies.