Permanent electric dipole moment of molybdenum carbide

High resolution optical spectroscopy has been used to study a molecular beam of molybdenum monocarbide (MoC). The Stark effect of the R(e)(0) and Q(fe)(1) branch features of the [18.6] (3)Pi(1)-X (3)Sigma(-)(0,0) band system of (98)MoC were analyzed to determine the permanent electric dipole moments mu(e) of 2.68(2) and 6.07(18) D for the [18.6] (3)Pi(1)(nu=0) and X (3)Sigma(-)(nu=0) states, respectively. The dipole moments are compared with the experimental value for ruthenium monocarbide [T. C. Steimle et al., J. Chem. Phys. 118, 2620 (2003)] and with theoretical predictions. A molecular orbital correlation diagram is used to interpret the observed and predicted trends of ground state mu(e) values for the 4d-metal monocarbides series.     

Differential virial theorem in relation to a sum rule for the exchange-correlation force in density-functional theory

Holas and March [Phys. Rev. A. 51, 2040 (1995)] gave a formally exact theory for the exchange-correlation (xc) force F(xc)(r)= -inverted Deltaupsilon(xc)(r) associated with the xc potential upsilon(xc)(r) of the density-functional theory in terms of low-order density matrices. This is shown in the present study to lead, rather directly, to the determination of a sum rule nF(xc)=0 relating the xc force with the ground-state density n(r). Some connection is also made with an earlier result relating to the external potential by Levy and Perdew [Phys. Rev. A. 32, 2010 (1985)] and with the quite recent study of Joubert [J. Chem. Phys. 119, 1916 (2003)] relating to the separation of the exchange and correlation contributions.     

Ab initio potential energy surfaces, bound states, and electronic spectrum of the Ar-SH complex

New ab initio potential energy surfaces for the (2)Pi ground electronic state of the Ar-SH complex are presented, calculated at the RCCSD(T)/aug-cc-pV5Z level. Weakly bound rotation-vibration levels are calculated using coupled-channel methods that properly account for the coupling between the two electronic states. The resulting wave functions are analyzed and a new adiabatic approximation including spin-orbit coupling is proposed. The ground-state wave functions are combined with those obtained for the excited (2)Sigma(+) state [D. M. Hirst, R. J. Doyle, and S. R. Mackenzie, Phys. Chem. Chem. Phys. 6, 5463 (2004)] to produce transition dipole moments. Modeling the transition intensities as a combination of these dipole moments and calculated lifetime values [A. B. McCoy, J. Chem. Phys. 109, 170 (1998)] leads to a good representation of the experimental fluorescence excitation spectrum [M.-C. Yang, A. P. Salzberg, B.-C. Chang, C. C. Carter, and T. A. Miller, J. Chem. Phys. 98, 4301 (1993)].     

Soft sticky dipole-quadrupole-octupole potential energy function for liquid water: an approximate moment expansion

A new, efficient potential energy function for liquid water is presented here. The new model, which is referred here as the soft sticky dipole-quadrupole-octupole (SSDQO) model, describes a water molecule as a Lennard-Jones sphere with point dipole, quadrupole, and octupole moments. It is a single-point model and resembles the hard-sphere sticky dipole potential model for water by Bratko et al. [J. Chem. Phys. 83, 6367 (1985)] and the soft sticky dipole model by Ichiye and Liu [J. Phys. Chem. 100, 2723 (1996)] except now the sticky potential consists of an approximate moment expansion for the dimer interaction potential, which is much faster than the true moment expansion. The object here is to demonstrate that the SSDQO potential energy function can accurately mimic the potential energy function of a multipoint model using the moments of that model. First, the SSDQO potential energy function using the dipole, quadruple, and octupole moments from SPC/E, TIP3P, or TIP5P is shown to reproduce the dimer potential energy functions of the respective multipoint model. In addition, in Monte Carlo simulations of the pure liquid at room temperature, SSDQO reproduces radial distribution functions of the respective model. However, the Monte Carlo simulations using the SSDQO model are about three times faster than those using the three-point models and the long-range interactions decay faster for SSDQO (1/r(3) and faster) than for multipoint models (1/r). Moreover, the contribution of each moment to the energetics and other properties can be determined. Overall, the simplicity, efficiency, and accuracy of the SSDQO potential energy function make it potentially very useful for studies of aqueous solvation by computer simulations.     

Experimental measurement of the induced dipole moment of an isolated molecule in its ground and electronically excited states: indole and indole-H2O

Reported here are measurements of the magnitude and orientation of the induced dipole moment that is produced when an indole molecule in its ground S(0) and electronically excited S(1) states is polarized by the attachment of a hydrogen bonded water molecule in the gas phase complex indole-H(2)O. For the complex, we find the permanent dipole moment values mu(IW)(S(0)) = 4.4 D and mu(IW)(S(1)) = 4.0 D, values that are substantially different from calculated values based on vector sums of the dipole moments of the component parts. From this result, we derive the induced dipole moment values mu(I) (*)(S(0)) = 0.7 D and mu(I) (*)(S(1)) = 0.5 D. The orientation of the induced moment also is significantly different in the two electronic states. These results are quantitatively reproduced by a purely electrostatic calculation based on ab initio values of multipole moments.     

Nonlinear optical property calculations of polyynes with long-range corrected hybrid exchange-correlation functionals

Polarizabilities (alpha), second-hyperpolarizabilities (gamma), and the gamma scaling factors (c) of polyynes [H-(C[triple bond]C)(n)-H, n = 1-8] were evaluated using the long-range corrected (LC) density functional theory (DFT) and LC-DFT with a short-range Gaussian attenuation (LCgau), as well as high quality wavefunction methods. We show that the c values obtained from LC- and LCgau-DFT are consistent with those from CCSD(T) calculations. Furthermore, the polyyne c values we obtained are seen to be smaller than the c values derived from previously calculated polyene gamma values [Sekino et al., J. Chem. Phys. 126, 014107 (2007)] in all the methods we consider. We compare our results with those obtained experimentally [Shepkov et al., J. Chem. Phys. 120, 6807 (2004).] from end-capped polyynes [i-Pr(3)Si-(C[triple bond]C)(n)-Sii-Pr(3)], which show larger c values for polyynes than polyenes. Our alpha and gamma calculations with i-Pr(3)Si-(C[triple bond]C)(n)-Sii-Pr(3) (n = 4, 5, 6, and 8) show that i-Pr(3)Si- may participate in pi molecular orbital delocalization, which can unexpectedly affect the c value. We also confirm the importance of molecular geometry in these nonlinear optical calculations. We find that while LC- and LCgau-DFT excellently reproduce experimental geometries and bond length alternation (BLA), MP2 optimized geometries have a BLA that is too short to be used for accurate alpha and gamma calculations.     

Molecular beam optical Stark study of rhodium mononitride

The optical Stark effect in the Q(1) and R(0) lines of the [15.1]1-X (1)Sigma+ (1,0) band of rhodium mononitride (RhN) were recorded and analyzed to determine the permanent electric dipole moments mu for the X (1)Sigma+(upsilon=0) and [15.1]1(upsilon=1) states to be 2.43(5) and 1.75(1) D, respectively. The determined dipole moments are compared to predicted values obtained from density functional theory [Stevens et al., Chem. Phys. Lett. 421, 281 (2006)] and an all-electron ab initio calculation [Shim et al., J. Mol. Struct. THEOCHEM 393, 127 (1997)]. A simple single configuration molecular orbital correlation diagram is used to rationalize the relative values of mu for the 4d mononitrides and RhO. An electronic configuration for the [15.1]1 state is proposed based on the interpretation of the (103)Rh and (14)N magnetic hyperfine interactions.     

Atomic dipole moments calculated using analytical molecular second-moment gradients

We have implemented analytical second-moment gradients for Hartree-Fock and multiconfigurational self-consistent-field wave functions. The code is used to calculate atomic dipole moments based on the generalized atomic polar tensor (GAPT) formalism [Phys. Rev. Lett. 62, 1469 (1989)], and the proposal of Dinur and Hagler (DH) for the calculation of atomic multipoles [J. Chem. Phys. 91, 2949 (1989)]. Both approaches display smooth basis-set convergence toward a well-defined basis-set limit and give reasonable electron correlation effects on the calculated atomic properties. However, the atomic charges and atomic dipole moments obtained from the GAPT partitioning scheme are unable to provide even qualitatively meaningful molecular quadrupole moments for some molecules, and thus the atomic multipole moments calculated in this scheme cannot be considered well suited for analyzing the electron density in molecules and for calculating intermolecular interaction energies. In contrast, the DH approach gives atomic charges and dipole moments that by definition exactly reproduce the molecular quadrupole moments. The approach of DH is, however, restricted to planar molecules and thus suffers from not being applicable to molecules of arbitrary shape. Both the GAPT and DH approaches give rather poor results for octupole and hexadecapole moments, indicating that at least atomic quadrupole moments are required for an accurate representation of the molecular charge distribution in terms of atomic electric moments.     

Amine-hydrogen halide complexes: experimental electric dipole moments and a theoretical decomposition of dipole moments and binding energies

The Stark effect has been observed in the rotational spectra of several gas-phase amine-hydrogen halide complexes and the following electric dipole moments have been determined: H(3)(15)N-H(35)Cl (4.05865 +/- 0.00095 D), (CH(3))(3)(15)N-H(35)Cl (7.128 +/- 0.012 D), H(3)(15)N-H(79)Br (4.2577 +/- 0.0022 D), and (CH(3))(3)(15)N-H(79)Br (8.397 +/- 0.014 D). Calculations of the binding energies and electric dipole moments for the full set of complexes R(n)()(CH(3))(3)(-)(n)()N-HX (n = 0-3; X = F, Cl, Br) at the MP2/aug-cc-pVDZ level are also reported. The block localized wave function (BLW) energy decomposition method has been used to partition the binding energies into contributions from electrostatic, exchange, distortion, polarization, and charge-transfer terms. Similarly, the calculated dipole moments have been decomposed into distortion, polarization, and charge-transfer components. The complexes studied range from hydrogen-bonded systems to proton-transferred ion pairs, and the total interaction energies vary from 7 to 17 kcal/mol across the series. The individual energy components show a much wider variation than this, but cancellation of terms accounts for the relatively narrow range of net binding energies. For both the hydrogen-bonded complexes and the proton-transferred ion pairs, the electrostatic and exchange terms have magnitudes that increase with the degree of proton transfer but are of opposite sign, leaving most of the net stabilization to arise from polarization and charge transfer. In all of the systems studied, the polarization terms contribute the most to the induced dipole moment, followed by smaller but still significant contributions from charge transfer. A significant contribution to the induced moment of the ion pairs also arises from distortion of the HX monomer.     

Molecular nanopolaritonics: cross manipulation of near-field plasmons and molecules. I. Theory and application to junction control

Near-field interactions between plasmons and molecules are treated in a simple unified approach. The density matrix of a molecule is treated with linear-response random phase approximation and the plasmons are treated classically. The equations of motion for the combined system are linear, governed by a simple Liouvillian operator for the polariton (plasmon+molecule excitation) dynamics. The dynamics can be followed in time or directly in frequency space where a trace formula for the transmission is presented. A model system is studied, metal dots in a forklike arrangement, coupled to a two level system with a large transition-dipole moment. A Fano-type resonance [Phys. Rev. 103, 1202 (1956)] develops when the molecular response is narrower than the width of the absorption spectrum for the plasmons. We show that the direction of the dipole of the molecule determines the direction the polariton chooses. Further, the precise position of the molecule has a significant effect on the transfer.     

The performance of semilocal and hybrid density functionals in 3d transition-metal chemistry

We investigate the performance of contemporary semilocal and hybrid density functionals for bond energetics, structures, dipole moments, and harmonic frequencies of 3d transition-metal (TM) compounds by comparison with gas-phase experiments. Special attention is given to the nonempirical metageneralized gradient approximation (meta-GGA) of Tao, Perdew, Staroverov, and Scuseria (TPSS) [Phys. Rev. Lett. 91, 146401 (2003)], which has been implemented in TURBOMOLE for the present work. Trends and error patterns for classes of homologous compounds are analyzed, including dimers, monohydrides, mononitrides, monoxides, monofluorides, polyatomic oxides and halogenides, carbonyls, and complexes with organic pi ligands such as benzene and cyclopentadienyl. Weakly bound systems such as Ca(2), Mn(2), and Zn(2) are discussed. We propose a reference set of reaction energies for benchmark purposes. Our all-electron results with quadruple zeta valence basis sets validate semilocal density-functional theory as the workhorse of computational TM chemistry. Typical errors in bond energies are substantially larger than in (organic) main group chemistry, however. The Becke-Perdew'86 [Phys. Rev. A 38, 3098 (1988); Phys. Rev. B 33, 8822 (1986)] GGA and the TPSS meta-GGA have the best price/performance ratio, while the TPSS hybrid functional achieves a slightly lower mean absolute error in bond energies. The popular Becke three-parameter hybrid B3LYP underbinds significantly and tends to overestimate bond distances; we give a possible explanation for this. We further show that hybrid mixing does not reduce the width of the error distribution on our reference set. The error of a functional for the s-d transfer energy of a TM atom does not predict its error for TM bond energies and bond lengths. For semilocal functionals, self-interaction error in one- and three-electron bonds appears to be a major source of error in TM reaction energies. Nevertheless, TPSS predicts the correct ground-state symmetry in the vast majority of cases and rarely fails qualitatively. This further confirms TPSS as a general purpose functional that works throughout the periodic table. We also give workstation timing comparisons for the 645-atom protein crambin.     

Molecular properties from variational reduced-density-matrix theory with three-particle N-representability conditions

Molecular ground-state energies and two-electron reduced density matrices (2-RDMs) have recently been computed without the many-electron wave function by constraining the 2-RDM to satisfy a complete set of three-positivity conditions for N representability [D. A. Mazziotti, Phys. Rev. A 74, 032501 (2006)]. Energies at both equilibrium and nonequilibrium geometries are obtained within 0.3% of the correlation energy. In this paper the authors extend this work to examine the accuracy of molecular properties, including multipole moments and components of the ground-state energy, relative to full configuration interaction (FCI). Comparisons are also made with 2-RDM methods with two-positivity conditions and two-positivity plus the generalized T1T2 conditions as well as several approximate wave function methods. Using the 2-RDM method with three-positivity conditions, the authors obtain dipole, quadrupole, and octupole moments for BeH2, BH, H2O, CO, and NH3 at equilibrium geometries that are within 0.04% of their FCI values. In addition, for the potential energy surface of N2, the 2-RDM method with three-positivity yields not only accurate total ground-state energies but also accurate expectation values of the kinetic energy operator, the electron-nuclei potential, and electron-electron repulsion.     

Dipole active vibrations and dipole moments of N2 and O2 physisorbed on a metal surface

We have, in infrared reflection absorption measurements, observed narrow dipole active absorption lines associated with the fundamental internal vibrational transitions of N(2) and O(2) physisorbed at 30 K on the chemically inert Pt(111)(1 x 1)H surface. Such transitions are forbidden for free homonuclear molecules and become dipole active at a metal surface due to polarization induced surface dipole moments. The measurements show that the internal stretch vibration frequencies are lowered by 7-8 cm(-1) relative to the gas phase values. The measured static and dynamic dipole moments are in the ranges of 0.06-0.07 and 0.001-0.002 D, respectively. We find that good estimates of the induced dynamic as well as the static dipole moments can in general be obtained from a van der Waals model but that the ratios of the measured static and dynamic moments indicates a need for a refinement of the dipole moment function.     

A second-quantization framework for the unified treatment of relativistic and nonrelativistic molecular perturbations by response theory

A formalism is presented for the calculation of relativistic corrections to molecular electronic energies and properties. After a discussion of the Dirac and Breit equations and their first-order Foldy-Wouthuysen [Phys. Rev. 78, 29 (1950)] transformation, we construct a second-quantization electronic Hamiltonian, valid for all values of the fine-structure constant alpha. The resulting alpha-dependent Hamiltonian is then used to set up a perturbation theory in orders of alpha(2), using the general framework of time-independent response theory, in the same manner as for geometrical and magnetic perturbations. Explicit expressions are given to second order in alpha(2) for the Hartree-Fock model. However, since all relativistic considerations are contained in the alpha-dependent Hamiltonian operator rather than in the wave function, the same approach may be used for other wave-function models, following the general procedure of response theory. In particular, by constructing a variational Lagrangian using the alpha-dependent electronic Hamiltonian, relativistic corrections can be calculated for nonvariational methods as well.     

Correlation energy functionals dependent on an effective number of electrons: charged species and equilibrium geometries

Recently proposed spin-dependent and spin-independent correlation energy functionals [Perez-Jimenez et al., J. Chem. Phys. 116, 10571 (2002)] based on an effective number of electrons N are extended to deal with charged systems. By introducing the concept of an effective atomic number Z analogous to N, the spin-dependent functional in combination with Becke's exchange [Becke, Phys. Rev. A 38, 3098 (1988)] yields a mean absolute error (MAE) of 5.4 kcal/mol for the 88 ionization potentials and 58 electron affinities included in the extended G2 set, and a MAE of 4.1 kcal/mol for the 312 data comprising the above plus the 148 enthalpies of formation of the extended G2 set and the 18 total energies of the neutral atoms H through Ar. Geometry optimizations performed on the 53 molecules of the G2-1 test set with the above combination of exchange and correlation functionals yield MAEs of 0.017 A and 1.5 degrees for the 68 bond lengths and 29 angles analyzed as compared with the experimental estimates.     

An experimental and theoretical study of dipole moments of N-[4-(9-acridinylamino)-3-methoxyphenyl]methanesulfonamide

The excited state (S1) dipole moment of m-AMSA (1), an acridine derivative with antitumor activity, was determined from solvatochromic shifts of the lowest energy absorption band in several organic solvents. The effect of the solute shape and the values of polarizability on the determined change of dipole moment between ground and excited state was discussed. The dipole moments in S0 and S1 state were calculated in gas phase with semiempirical quantum-chemical and DFT and CIS methods and in solvents with SM5.4A solvation model and compared with values obtained experimentally. All the results show that the dipole moment of compound 1 in the excited state is higher than that in the ground state. These methods quite well predict the values of Deltamicro between two states of an investigated compound.     

Retarded boson-fermion interaction in atomic systems

The retarded interaction between an electron and a spin-0 nucleus is derived from electrodynamical perturbation theory. The contribution of retardation at order v(2)c(2) mimics the Breit interaction [Phys. Rev. 34, 553 (1929); 36, 388 (1930); 39, 616 (1932)] with the Dirac matrix alpha(2) being replaced by p(2)m(2)c where p(2) is the linear momentum operator for the nucleus. An effective one-electron retardation operator is obtained in relative coordinates, and this can be used through all orders in perturbation theory without any problem of infinite degeneracy. A few steps of unitary transformation lead to the nonrelativistic limit. The leading terms in retardation corrections to energy are of order (m(e)m(n))alpha(2)Z(4)(alpha(2)m(e)c(2)). The implications for atomic systems are discussed.     

Effect of nitro groups on the photo physical properties of benzimidazolone: a solvatochromic study

Electronic absorption and fluorescence spectra of mono, di, and tri-nitro benzimidazolones are measured at room temperature (298 K) in nine solvents with different polarities and the observed shifts are compared with benzimidazolone. Ground and excited state electric dipole moments are determined using the solvatochromic method based on the bulk solvent properties, F(1)(ε, n) and F(2)(ε, n). A reasonable agreement is observed between the experimental and ab initio dipole moments. Change in dipole moment is also determined using the solvatochromic method based on the microscopic solvent polarity parameter, (E(T)(N)), which considers the polarization changes due to hydrogen bonding in different solvents. It has been observed that the correlation of the solvatochromic Stokes shifts with the parameter (E(T)(N)), is superior to that derived using bulk solvent polarity functions for all the benzimidazolones reported in the present study. Calculated difference between excited state and ground state dipole moments seems to be a good measure of the effect of nitro group when correlated with (E(T)(N)).     

Dielectric Relaxation Around a Charged Colloidal Cylinder in an Electrolyte

The polarizability and corresponding dielectric relaxation of the Debye-Hückel (DH) atmosphere surrounding a charged rod-like polyelectrolyte immersed in an ionic solution of a symmetrical electrolyte is determined following the method developed by J. A. Fornés [Phys. Rev. E 57, 2110 (1998)]. Several formulas are given to estimate the DH atmosphere parameters, namely, the polarizability at zero frequency, alpha(0), the relaxation time, tau, the cloud capacitance, C, the average displacement of the ionic cloud, delta, the square root dipole moment quadratic fluctuation, (1/2), and the thermal fluctuating field, (1/2). The Poisson-Boltzmann equation is solved numerically to apply the theory to a highly charged polyelectrolyte such as DNA in solution, although formulas valid for the DH approximation are also given. A dispersion in the polarizability and correspondingly in the dielectric constant of these solutions in the microwave region is predicted. For instance, considering a DNA length of 1000 ?, with its reduced linear charge density xi(0)=4.25 and ionization factor gamma=0.5, immersed in a NaCl solution (40 mM), we predict a polarizability of the DH atmosphere at zero frequency alpha(0) of 1x10(-33) Fm(2) ( approximately 6.1x10(6)) times greater than the mean value of the polarizability of water) and the corresponding fluctuating dipole moment p of 2.1x10(-27) Cm ( approximately 600 times greater than the permanent dipole moment of water molecule). The relaxation time and the average displacement of the ionic cloud are tau=1.6 ns and delta=14. ?, respectively. This displacement is produced by the thermal fluctuating field, which, in this case, at room temperature is (1/2)=2 x10(6) V/m. Copyright 2000 Academic Press.