A treatment of atomic mean square displacement in higher order perturbation theory

A general derivation of the anharmonic coefficients for a periodic lattice invoking the special case of the central force interaction is presented. All of the contributions to mean square displacement (MSD) to order 14 perturbation theory are enumerated. A direct correspondance is found between the high temperature limit MSD and high temperature limit free energy contributions up to and including 0(14). This correspondance follows from the detailed derivation of some of the contributions to MSD. Numerical results are obtained for all the MSD contributions to 0(14) using the Lennard-Jones potential for the lattice constants and temperatures for which the Monte Carlo results were calculated by Heiser, Shukla and Cowley. The Peierls approximation is also employed in order to simplify the numerical evaluation of the MSD contributions. The numerical results indicate the convergence of the perturbation expansion up to 75% of the melting temperature of the solid (TM) for the exact calculation; however, a better agreement with the Monte Carlo results is not obtained when the total of all 14 contributions is added to the 12 perturbation theory results. Using Peierls approximation the expansion converges up to 45% of TM• The MSD contributions arising in the Green's function method of Shukla and Hubschle are derived and enumerated up to and including 0(18). The total MSD from these selected contributions is in excellent agreement with their results at all temperatures. Theoretical values of the recoilless fraction for krypton are calculated from the MSD contributions for both the Lennard-Jones and Aziz potentials. The agreement with experimental values is quite good.

The domain derivative in rough-surface scattering and rigorous estimates for first-order perturbation theory

We investigate Fréchet differentiability of the scattered field with respect to variation in the boundary in the case of time–harmonic acoustic scattering by an unbounded, sound–soft, one–dimensional rough surface. We rigorously prove the differentiability of the scattered field and derive a characterization of the Fréchet derivative as the solution to a Dirichlet boundary value problem. As an application of these results we give rigorous error estimates for first–order perturbation theory, justifying small perturbation methods that have a long history in the engineering literature. As an application of our rigorous estimates we show that a plane acoustic wave incident on a sound–soft rough surface can produce an unbounded scattered field.

Geocenter coordinates estimated from GNSS data as viewed by perturbation theory

Time series of geocenter coordinates were determined with data of two global navigation satellite systems (GNSSs), namely the U.S. GPS (Global Positioning System) and the Russian GLONASS (Global’naya Nawigatsionnaya Sputnikowaya Sistema). The data was recorded in the years 2008–2011 by a global network of 92 permanently observing GPS/GLONASS receivers. Two types of daily solutions were generated independently for each GNSS, one including the estimation of geocenter coordinates and one without these parameters. A fair agreement for GPS and GLONASS was found in the geocenter x- and y-coordinate series. Our tests, however, clearly reveal artifacts in the z-component determined with the GLONASS data. Large periodic excursions in the GLONASS geocenter z-coordinates of about 40 cm peak-to-peak are related to the maximum elevation angles of the Sun above/below the orbital planes of the satellite system and thus have a period of about 4 months (third of a year). A detailed analysis revealed that the artifacts are almost uniquely governed by the differences of the estimates of direct solar radiation pressure (SRP) in the two solution series (with and without geocenter estimation). A simple formula is derived, describing the relation between the geocenter z-coordinate and the corresponding parameter of the SRP. The effect can be explained by first-order perturbation theory of celestial mechanics. The theory also predicts a heavy impact on the GNSS-derived geocenter if once-per-revolution SRP parameters are estimated in the direction of the satellite’s solar panel axis. Specific experiments using GPS observations revealed that this is indeed the case. Although the main focus of this article is on GNSS, the theory developed is applicable to all satellite observing techniques. We applied the theory to satellite laser ranging (SLR) solutions using LAGEOS. It turns out that the correlation between geocenter and SRP parameters is not a critical issue for the SLR solutions. The reasons are threefold: The direct SRP is about a factor of 30–40 smaller for typical geodetic SLR satellites than for GNSS satellites, allowing it in most cases to not solve for SRP parameters (ruling out the correlation between these parameters and the geocenter coordinates); the orbital arc length of 7 days (which is typically used in SLR analysis) contains more than 50 revolutions of the LAGEOS satellites as compared to about two revolutions of GNSS satellites for the daily arcs used in GNSS analysis; the orbit geometry is not as critical for LAGEOS as for GNSS satellites, because the elevation angle of the Sun w.r.t. the orbital plane is usually significantly changing over 7 days.

Full-dimensional analytical ab initio potential energy surface of the ground state of HOI

Extensive ab initio calculations using a complete active space second-order perturbation theory wavefunction, including scalar and spin-orbit relativistic effects with a quadruple-zeta quality basis set were used to construct an analytical potential energy surface (PES) of the ground state of the [H, O, I] system. A total of 5344 points were fit to a three-dimensional function of the internuclear distances, with a global root-mean-square error of 1.26 kcal mol(-1). The resulting PES describes accurately the main features of this system: the HOI and HIO isomers, the transition state between them, and all dissociation asymptotes. After a small adjustment, using a scaling factor on the internal coordinates of HOI, the frequencies calculated in this work agree with the experimental data available within 10 cm(-1). (C) 2011 American Institute of Physics. [doi: 10.1063/1.3615545]

Electronic Structure of the Ground and Low-Lying Electronic States of MoB and MoB(+)

Multiconfigurational second-order perturbation theory (CASSCF//CASPT2) and quadruple-zeta ANO-RCC basis sets were employed to investigate the ground and low-lying electronic states of MoB and MoB(+). Spectroscopic constants, potential energy curves, wavefunctions, Mulliken population analyses, and ionization energies are given. The ground state of MoB is of X(6)Pi symmetry (R(e) = 1.968 angstrom, omega(e) = 664 cm(-1), and mu = 2.7 D), giving rise to a Omega = 7/2 ground state after including spin-orbit coupling. For MoB(+), the ground state is computed to be of X(7)Sigma(+) symmetry (R(e) = 2.224 angstrom, omega(e) = 141 cm(-1), and mu = 1.2 D), with an adiabatic ionization energy of 7.19 eV and a vertical one of 7.53 eV. (C) 2011 Wiley Periodicals, Inc. Int J Quantum Chem 111: 3362-3370, 2011

Solvatochromic Shift of the pi -> pi* Transition in All-Trans, Cis-13, Cis-11, Cis-9, and Cis-7 Retinal Isomers Induced by Water and Methanol

The solvatochromic shift of the lowest singlet it pi -> pi* electronic transition in the all-trans, cis-13, cis-11, cis-9, and cis-7 retinal isomers were computed under the influence of water, methanol, and benzene solvents. Excitation energies were calculated in gas phase and in solution. The calculations in solution were performed considering the sequential Monte Carlo (MC) /Quantum Mechanical approach. The MC simulations were performed considering the full retinal isomer molecules and 900 water molecules, 900 methanol, or 400 benzene ones. The OPLS/AA parametrization was chosen for retinal, methanol, and benzene molecules and the SPC model was used for water one. From the MC calculations 100 independent configurations were selected, with 100 solvent molecules in thermodynamical equilibrium at T = 298.15 K. Average point-charges were obtained from those independent configurations for water, methanol, and benzene solvent. TDDFT and CASSCF//CASPT2 methodologies were used to compute the vertical excitation energy of the retinal isomers in different environment. (C) 2010 Wiley Periodicals, Inc. Int J Quantum Chem 110: 2076-2087, 2010

Electronic structure and chemical bonding in W(2) molecule

The electronic structure of the lowest-lying electronic states of W(2) were investigated at the CASPT2 level. The ground state is a X(1)Sigma(+)(g) state, followed by the a(3)Delta(u), b(3)Sigma(+)(u) and A(1)Delta(u) electronic states. Seven low-lying Omega-states were computed: (1)0(g)(+), (2)3(u), (3)2(u), (4)1(u), (5)0(u)(-), (6)1(u), and (7)2(u), with the ground state corresponding to the (1)0(g)(+)(X(1)Sigma(+)(g)) state. Comparison with the other VIB transition metal group dimers indicates a common pattern of electronic structure and spectroscopic properties. (C) 2010 Elsevier B.V. All rights reserved.

Electronic Structure and Chemical Bonding in the Lowest Electronic States of TcN

Multiconfiguration second-order perturbation theory, with the inclusion of relativistic effects and spin-orbit Coupling, was employed to investigate the nature of the ground and low-lying Lambda-S and Omega states of the TcN molecule. Spectroscopic constants, effective bond order, and potential energy curves for 13 low-lying Lambda-S states and 5 Omega states are given, The computed ground state of TcN is of Omega = 3 symmetry (R(e) = 1.605 angstrom and omega(e) = 1085 cm(-1)), originating mainly from the (3)Delta Lambda-S ground state. This result is contrasted with the nature of the ground state for other VIIB transtion-metal mononitrides, including X(3)Sigma(-) symmetry for MnN and Omega = 0(+) symmetry for ReN, derived also from a X(3)Sigma(-) state.

The intramolecular charge transfer in a donor-pi-acceptor dianion probed by resonance Raman spectroscopy and quantum chemical calculations

The dideprotonation of 4-(4-nitrophenylazo)resorcinol generates an anionic species with substantial electronic pi delocalization. As compared to the parent neutral species, the anionic first excited electronic transition, characterized as an intramolecular charge transfer (ICT) from the CO(-) groups to the NO(2) moiety, shows a drastic red shift of ca. 200 nm in the lambda(max) in the UV-vis spectrum, leading to one of the lowest ICT energies observed (lambda(max) = 630 nm in dimethyl sulfoxide (DMSO)) in this class of push-pull molecular systems. Concomitantly, a threefold increase in the molar absorptivity (epsilon(max)) in comparison to the neutral species is observed. The resonance Raman enhancement profiles reveal that in the neutral species the chromophore involves several modes, as nu(C-N), nu(N=N), nu(C=C) and nu(s)(NO(2)), whereas in the dianion, there is a selective enhancement of the NO(2) vibrational modes. The quantum chemical calculations of the electronic transitions and vibrational wavenumbers led to a consistent analysis of the enhancement patterns observed in the resonance Raman spectra. Copyright (C) 2009 John Wiley & Sons, Ltd.

Electronic structure and chemical bonding in the ground states of Tc(2) and Re(2)

Multiconfiguration second-order perturbation theory, including relativistic effects and spin-orbit coupling, has been employed to investigate the nature of the chemical bonding in the ground state of Tc(2) and Re(2). The Tc(2) ground state is found to be a 0(g)(+) state, with an effective bond order (EBO) of 4.4, and a dissociation energy of 3.25 eV. The Re(2) ground state is a 1(g) state, with EBO = 4.3. Almost degenerate to it, is a 0(g)(+) state (T(e) = 77 cm(-1)), with EBO = 4.1. Experimental evidence also indicates that the ground state is of 1(g) nature. The dissociation energy is computed to be 5.0 eV in agreement with an experimental estimate of 4 +/- 1 eV.

New Relativistic Atomic Natural Orbital Basis Sets for Lanthanide Atoms with Applications to the Ce Diatom and LuF3

New basis sets of the atomic natural orbital (ANO) type have been developed for the lanthanide atoms La-Lu. The ANOs have been obtained from the average density matrix of the ground and lowest excited states of the atom, the positive ions, and the atom in an electric field. Scalar relativistic effects are included through the use of a Douglas-Kroll-Hess Hamiltonian. Multiconfigurational wave functions have been used with dynamic correlation included using second-order perturbation theory (CASSCF/CASPT2). The basis sets are applied in calculations of ionization energies and some excitation energies. Computed ionization energies have an accuracy better than 0.1 eV in most cases. Two molecular applications are inluded as illustration: the cerium diatom and the LuF3 molecule. In both cases it is shown that 4f orbitals are not involved in the chemical bond in contrast to an earlier claim for the latter molecule.

Electronic Structure and Chemical Bonding in the Ground and Low-Lying Electronic States of Ta(2)

The electronic structure and chemical bonding of the ground and low-lying Lambda - S and Omega states of Ta(2) were investigated at the multiconfiguration second-order perturbation theory (CASSCF//CASPT2) level. The ground state of Ta(2) is computed to be a X(3)Sigma(-)(g) state (R(e) = 2.120 angstrom, omega(e) = 323 cm(-1), and D(e) = 4.65 eV), with two low-lying singlet states close to it (a(1) Sigma(+)(g) : T(e) = 409 cm(-1), R(e) = 2.131 angstrom, and omega(e) = 313 cm(-1); b(1) Gamma(g): T(e) = 1, 038 cm(-1), R(e) = 2.127 angstrom, and omega(e) = 316 cm(-1)). These electronic states are derived from the same electronic configuration: vertical bar 13 sigma(2)(g)14 sigma(2)(g)7 delta(2)(g)13 pi(4)(u)>. The effective bond order of the X(3) Sigma(-)(g) state is 4.52, which indicates that the Ta atoms are bound by a quintuple chemical bond. The a(1) Sigma(+)(g) state interacts strongly with the X(3)Sigma(-)(g) g ground state by a second-order spin-orbit interaction, giving rise to the (1)0(g)(+) (ground state) (dominated by the X(3)Sigma(-)(g) Lambda - S ground state) and (9)0(g)(+) (dominated by the a(1) Sigma(+)(g) Lambda - S state) Omega states. These results are in line with those reported for the group 5B homonuclear transition metal diatomics. (C) 2010 Wiley Periodicals, Inc. Int J Quantum Chem 111: 1306-1315, 2011

Azo-hydrazone tautomerism in protonated aminoazobenzenes: Resonance Raman spectroscopy and quantum-chemical calculations

The protonation effect on the vibrational and electronic spectra of 4-aminoazobenzene and 4-(dimethylamino)azobenzene was investigated by resonance Raman spectroscopy, and the results were discussed on the basis of quantum-chemical calculations. Although this class of molecular systems has been investigated in the past concerning the azo-hydrazone tautomerism, the present work is the first to use CASSCF/CASPT2 calculations to unveil the structure of both tautomers as well the nature of the molecular orbitals involved in chromophoric moieties responsible for the resonance Raman enhancement patterns. More specifically both the resonance Raman and theoretical results show clearly that in the neutral species, the charge transfer transition involves mainly the azo moiety, whereas in the protonated forms there is a great difference, depending on the tautomer. In fact, for the azo tautomer the transition is similar to that observed in the corresponding neutral species, whereas in the hydrazone tautomer such a transition is much more delocalized due to the contribution of the quinoid structure. The characterization of protonated species and the understanding of the tautomerization mechanism are crucial for controlling molecular properties depending on the polarity and pH of the medium.

Low-lying singlet and triplet electronic states of RhB

The low-lying X-1 Sigma(+), a(3)Delta, A(1)Delta, b(3)Sigma(+), B-1 Pi, c(3)Pi, C-1 Phi, D-1 Sigma(+), E-1 Pi, d(1)Phi, and e(3)Pi electronic states of RhB have been investigated at the ab initio level, using the multistate multiconfigurational second-order perturbation (MS-CASPT2) theory, with extended atomic basis sets and inclusion of scalar relativistic effects. Among the eleven electronic states included in this work, only three (the X-1 Sigma(+), D-1 Sigma(+), and E-1 Pi states) have been investigated experimentally. Potential energy curves, spectroscopic constants, dipole moments, binding energies, and chemical bonding aspects are presented for all electronic states.

A theoretical study of the binding and electronic spectrum of the Mo-2 molecule

Multiconfigurational SCF and second-order perturbation theory have been employed to study seven low-lying singlet and triplet electronic states of the Mo-2 molecule. The bond order of the ground state has been analyzed based on the effective bond order (EBO), indicating that a fully developed sextuple bond is formed between the two Mo atoms. The experimentally observed excited states a(3)Sigma(+)(u) and A(1)Sigma(+)(u) have been determined and the so-called (3)Lambda excited state identified as the b(3)Sigma(+)(u) state, in agreement with experimental expectations. (C) 2007 Elsevier B.V. All rights reserved.