Ab initio benchmark study for the oxidative addition of CH4 to Pd: importance of basis-set flexibility and polarization

To obtain a state-of-the-art benchmark potential energy surface (PES) for the archetypal oxidative addition of the methane C-H bond to the palladium atom, we have explored this PES using a hierarchical series of ab initio methods (Hartree-Fock, second-order Møller-Plesset perturbation theory, fourth-order Møller-Plesset perturbation theory with single, double and quadruple excitations, coupled cluster theory with single and double excitations (CCSD), and with triple excitations treated perturbatively [CCSD(T)]) and hybrid density functional theory using the B3LYP functional, in combination with a hierarchical series of ten Gaussian-type basis sets, up to g polarization. Relativistic effects are taken into account either through a relativistic effective core potential for palladium or through a full four-component all-electron approach. Counterpoise corrected relative energies of stationary points are converged to within 0.1-0.2 kcal/mol as a function of the basis-set size. Our best estimate of kinetic and thermodynamic parameters is -8.1 (-8.3) kcal/mol for the formation of the reactant complex, 5.8 (3.1) kcal/mol for the activation energy relative to the separate reactants, and 0.8 (-1.2) kcal/mol for the reaction energy (zero-point vibrational energy-corrected values in parentheses). This agrees well with available experimental data. Our work highlights the importance of sufficient higher angular momentum polarization functions, f and g, for correctly describing metal-d-electron correlation and, thus, for obtaining reliable relative energies. We show that standard basis sets, such as LANL2DZ+ 1f for palladium, are not sufficiently polarized for this purpose and lead to erroneous CCSD(T) results. B3LYP is associated with smaller basis set superposition errors and shows faster convergence with basis-set size but yields relative energies (in particular, a reaction barrier) that are ca. 3.5 kcal/mol higher than the corresponding CCSD(T) values

Notations and conventions in molecular spectroscopy: part 1. General spectroscopic notation

The field of Molecular Spectroscopy was surveyed in order to determine a set of conventions and symbols which are in common use in the spectroscopic literature. This document, which is Part I in a series, establishes the notations and conventions used for general spectroscopic notations and deals with quantum mechanics, quantum numbers (vibrational states, angular momentum and energy levels), spectroscopic transitions, and miscellaneous notations (e.g. spectroscopic terms). Further parts will follow, dealing inter alia with symmetry notation, permutation and permutation-inversion symmetry notation, vibration-rotation spectroscopy and electronic spectroscopy.

Malacosporean-like spores in urine of rainbow trout react with antibody and DNA probes to Tetracapsuloides bryosalmonae

Tetracapsuloides bryosalmonae is the myxozoan parasite causing proliferative kidney disease (PKD) of salmonid fishes in Europe and North America. The complete life cycle of the parasite remains unknown despite recent discoveries that the stages infectious for fish develop in freshwater bryozoans. During the course of examinations of the urine of rainbow trout (Oncorhynchus mykiss) with or recovering from PKD we identified spores with features similar to those of T. bryosalmonae found in the bryozoan host. Spores found in the urine were subspherical, with a width of 16 mum and height of 14 mum, and possessed two soft valves surrounding two spherical polar capsules (2 mum in diameter) and a single sporoplasm. The absence of hardened valves is a distinguishing characteristic of the newly established class Malacosporea that includes T. bryosalmonae as found in the bryozoan host. The parasite in the urine of rainbow trout possessed only two polar capsules and two valve cells compared to the four polar capsules and four valves observed in the spherical spores of 19 mum in diameter from T. bryosalmonae from the bryozoan host. Despite morphological differences, a relationship between the spores in the urine of rainbow trout and T. bryosalmonae was demonstrated by binding of monoclonal and polyclonal antibodies and DNA probes specific to T. bryosalmonae.

Tests of MULTIMODE calculations of rovibrational energies of CH4

We report variational calculations of rovibrational energies of CH4 using the code MULTIMODE and an ab initio force field of Schwenke and Partridge. The systematic convergence of the energies with respect to the level of mode coupling is presented. Converged vibrational energies calculated using the five-mode representation of the potential for zero total angular momentum are compared with previous, benchmark calculations based on Radau coordinates using this force field for zero total angular momentum and for J = 1. Very good agreement with the previous benchmark calculations is found. (c) 2006 Elsevier B.V. All rights reserved.

Atmospheric conservation properties in ERA-Interim

We study the global atmospheric budgets of mass, moisture, energy and angular momentum in the latest reanalysis from the European Centre for Medium-Range Weather Forecasts (ECMWF), ERA-Interim, for the period 1989–2008 and compare with ERA-40. Most of the measures we use indicate that the ERA-Interim reanalysis is superior in quality to ERA-40. In ERA-Interim the standard deviation of the monthly mean global dry mass of 0.7 kg m−2 (0.007%) is slightly worse than in ERA-40, and long time-scale variations in dry mass originate predominately in the surface pressure field. The divergent winds are improved in ERA-Interim: the global standard deviation of the time-averaged dry mass budget residual is 10 kg m−2 day−1 and the quality of the cross-equatorial mass fluxes is improved. The temporal variations in the global evaporation minus precipitation (E − P) are too large but the global moisture budget residual is 0.003 kg m−2 day−1 with a spatial standard deviation of 0.3 kg m−2 day−1. Both the E − P over ocean and P − E over land are about 15% larger than the 1.1 Tg s−1 transport of water from ocean to land. The top of atmosphere (TOA) net energy losses are improved, with a value of 1 W m−2, but the meridional gradient of the TOA net energy flux is smaller than that from the Clouds and the Earth's Radiant Energy System (CERES) data. At the surface the global energy losses are worse, with a value of 7 W m−2. Over land however, the energy loss is only 0.5 W m−2. The downwelling thermal radiation at the surface in ERA-Interim of 341 W m−2 is towards the higher end of previous estimates. The global mass-adjusted energy budget residual is 8 W m−2 with a spatial standard deviation of 11 W m−2, and the mass-adjusted atmospheric energy transport from low to high latitudes (the sum for the two hemispheres) is 9.5 PW

Assessing wind profile effects on the global atmospheric torque

The impact of a new approach to the evaluation of surface gravity wave drag (GWD) is assessed. This approach uses linear theory, but incorporates the effects of wind profile shear and curvature, by means of a second-order WKB approximation. While the theory predicts the possibility of either drag enhancement or reduction, depending on the wind profile, results obtained with the ERA-40 reanalysis data clearly indicate the predominance of local drag enhancement. However, the global impact of shear on the atmospheric axial GWD torque comes mostly from regions with predominantly easterly flow, contributing to a slight reduction of the bias found in different studies of the global angular momentum budget. The relative correction due to shear on linear GWD is found not to depend too strongly on the levels chosen for the computation of the low-level wind derivatives.

Symmetric stability of compressible zonal flows on a generalized equatorial β plane

Sufficient conditions are derived for the linear stability with respect to zonally symmetric perturbations of a steady zonal solution to the nonhydrostatic compressible Euler equations on an equatorial � plane, including a leading order representation of the Coriolis force terms due to the poleward component of the planetary rotation vector. A version of the energy–Casimir method of stability proof is applied: an invariant functional of the Euler equations linearized about the equilibrium zonal flow is found, and positive definiteness of the functional is shown to imply linear stability of the equilibrium. It is shown that an equilibrium is stable if the potential vorticity has the same sign as latitude and the Rayleigh centrifugal stability condition that absolute angular momentum increase toward the equator on surfaces of constant pressure is satisfied. The result generalizes earlier results for hydrostatic and incompressible systems and for systems that do not account for the nontraditional Coriolis force terms. The stability of particular equilibrium zonal velocity, entropy, and density fields is assessed. A notable case in which the effect of the nontraditional Coriolis force is decisive is the instability of an angular momentum profile that decreases away from the equator but is flatter than quadratic in latitude, despite its satisfying both the centrifugal and convective stability conditions.

On the available energy of an axisymmetric vortex

A theory of available energy for axisymmetric circulations is presented. The theory is a generalization of the classical theory of available potential energy, in that it accounts for both thermal and angular momentum constraints on the circulation. The generalization relies on the Hamiltonian structure of the (conservative) dynamics, is exact at finite amplitude, and has a local form. Application of the theory is presented for the case of an axisymmetric vortex on an f -plane in the context of the Boussinesq equations.

Issues in stratosphere-troposphere coupling

An overview is given of current issues concerning the coupling between the stratosphere and troposphere. The tropopause region, more generally the upper troposphere/lower stratosphere, is the region of direct contact where exchange of material takes place. Dynamical coupling through angular momentum transfer by waves occurs nonlocally, and provides a generally negative torque on the stratosphere which drives an equator-to-pole circulation (i.e., towards the Earth’s axis of rotation). This wave-driven circulation is the principal mechanism for intraseasonal and interannual variability in the extratropical stratosphere. Although such variability is generally dynamical in origin, there are important chemical and radiative feedbacks. The location of the tropopause has implications for radiative forcing of climate, through its effect on the distribution of relatively short-lived greenhouse gases (ozone and water vapour). Some outstanding puzzles in our current understanding are identified. Attention is focused on possible climate sensitivities, and how these may be tested and constrained. Results from the Canadian Middle Atmosphere Model (CMAM), a fully interactive radiative-chemical-dynamical general circulation model, are used to illustrate some of the points.

Mechanisms for tropical upwelling in the stratosphere

The dynamics of the tropical upwelling branch of the stratospheric Brewer–Dobson circulation are examined, with a particular focus on the role of the middle-atmosphere Hadley circulation. Upwelling is examined in terms of both the diabatic circulation and Lagrangian trajectories using a zonally symmetric balance model. The behavior of the wave-driven circulation in the presence of angular momentum redistribution by the Hadley circulation is also considered. The results of the zonally symmetric model are compared with fields from a middle-atmosphere GCM. It is found that the Hadley circulation makes a significant contribution to annual mean tropical upwelling at the upwelling maximum in the vicinity of the stratopause, and can account for most of the annual mean upwelling seen in the GCM. In the mid- to lower stratosphere, the role of the Hadley circulation is much weaker and wave drag appears to be required to explain the observed upwelling, although the Hadley circulation makes a nonnegligible contribution to the annual cycle of the upwelling. Subtropical wave drag can produce annual mean upwelling through a nonlinear mechanism; viscosity is not required. However, the magnitude of the observed upwelling suggests that wave drag must penetrate quite close to the equator.

The middle-atmosphere Hadley circulation and equatorial inertial adjustment

In the tropical middle atmosphere the climatological radiative equilibrium temperature is inconsistent with gradient-wind balance and the available angular momentum, especially during solstice seasons. Adjustment toward a balanced state results in a type of Hadley circulation that lies outside the “downward control” view of zonally averaged dynamics. This middle-atmosphere Hadley circulation is reexamined here using a zonally symmetric balance model driven through an annual cycle. It is found that the inclusion of a realistic radiation scheme leads to a concentration of the circulation near the stratopause and to its closing off in the mesosphere, with no need for relaxational damping or a rigid lid. The evolving zonal flow is inertially unstable, leading to a rapid process of inertial adjustment, which becomes significant in the mesosphere. This short-circuits the slower process of angular momentum homogenization by the Hadley circulation itself, thereby weakening the latter. The effect of the meridional circulation associated with extratropical wave drag on the Hadley circulation is considered. It is shown that the two circulations are independent for linear (quasigeostrophic) zonal-mean dynamics, and interact primarily through the advection of temperature and angular momentum. There appears to be no significant coupling in the deep Tropics via temperature advection since the wave-driven circulation is unable to alter meridional temperature gradients in this region. However, the wave-driven circulation can affect the Hadley circulation by advecting angular momentum out of the Tropics. The validity of the zonally symmetric balance model with parameterized inertial adjustment is tested by comparison with a three-dimensional primitive equations model. Fields from a middle-atmosphere GCM are also examined for evidence of these processes. While many aspects of the GCM circulation are indicative of the middle-atmosphere Hadley circulation, particularly in the upper stratosphere, it appears that the circulation is obscured in the mesosphere and lower stratosphere by other processes.

The middle atmosphere

The last 50 years have seen enormous advances in our knowledge and understanding of the stratosphere and mesosphere, which together comprise the middle atmosphere. Beginning from a phase of basic discovery, we have now reached the stage where most observed phenomena can be modelled from first principles with a reasonable degree of fidelity, and where there is an overall theoretical framework which can be tested against measurements and models. This review surveys a number of major surprises in middle atmosphere science over the past 50 years. A phenomenological and historical approach is adopted in each case, leading up to the current literature. Along the way, a common thread emerges: the central role of waves, of various types, in redistributing angular momentum within the atmosphere, and the global nature of the atmospheric response to such redistribution

On the 'downward control' of extratropical diabatic circulations by eddy-induced mean zonal forces

The situation considered is that of a zonally symmetric model of the middle atmosphere subject to a given quasi-steady zonal force F̄, conceived to be the result of irreversible angular momentum transfer due to the upward propagation and breaking of Rossby and gravity waves together with any other dissipative eddy effects that may be relevant. The model's diabatic heating is assumed to have the qualitative character of a relaxation toward some radiatively determined temperature field. To the extent that the force F̄ may be regarded as given, and the extratropical angular momentum distribution is realistic, the extratropical diabatic mass flow across a given isentropic surface may be regarded as controlled exclusively by the F̄ distribution above that surface (implying control by the eddy dissipation above that surface and not, for instance, by the frequency of tropopause folding below). This “downward control” principle expresses a critical part of the dynamical chain of cause and effect governing the average rate at which photochemical products like ozone become available for folding into, or otherwise descending into, the extratropical troposphere. The dynamical facts expressed by the principle are also relevant, for instance, to understanding the seasonal-mean rate of upwelling of water vapor to the summer mesopause, and the interhemispheric differences in stratospheric tracer transport. The robustness of the principle is examined when F̄ is time-dependent. For a global-scale, zonally symmetric diabatic circulation with a Brewer-Dobson-like horizontal structure given by the second zonally symmetric Hough mode, with Rossby height HR = 13 km in an isothermal atmosphere with density scale height H = 7 km, the vertical partitioning of the unsteady part of the mass circulation caused by fluctuations in F̄ confined to a shallow layer LF̄ is always at least 84% downward. It is 90% downward when the force fluctuates sinusoidally on twice the radiative relaxation timescale and 95% if five times slower. The time-dependent adjustment when F̄ is changed suddenly is elucidated, extending the work of Dickinson (1968), when the atmosphere is unbounded above and below. Above the forcing, the adjustment is characterized by decay of the meridional mass circulation cell at a rate proportional to the radiative relaxation rate τr−1 divided by {1 + (4H2/HR2)}. This decay is related to the boundedness of the angular momentum that can be taken up by the finite mass of air above LF̄ without causing an ever-increasing departure from thermal wind balance. Below the forcing, the meridional mass circulation cell penetrates downward at a speed τr−1 HR2/H. For the second Hough mode, the time for downward penetration through one density scale height is about 6 days if the radiative relaxation time is 20 days, the latter being representative of the lower stratosphere. At any given altitude, a steady state is approached. The effect of a rigid lower boundary on the time-dependent adjustment is also considered. If a frictional planetary boundary layer is present then a steady state is ultimately approached everywhere, with the mass circulation extending downward from LF̄ and closing via the boundary layer. Satellite observations of temperature and ozone are used in conjunction with a radiative transfer scheme to estimate the altitudes from which the lower stratospheric diabatic vertical velocity is controlled by the effective F̄ in the real atmosphere. The data appear to indicate that about 80% of the effective control is usually exerted from below 40 km but with significant exceptions up to 70 km (in the high latitude southern hemispheric winter). The implications for numerical modelling of chemical transport are noted.

The effect of the magnetic field of the secondary star in dwarf novae

We investigate the effect of a secondary star magnetic field on the accretion disc dynamics of dwarf novae. Simulations have been carried out with a particle code and a dipolar magnetic field structure. The magnetic field acts to remove angular momentum from the disc material, increasing the inward mass flow. This makes the accretion disc more centrally condensed, causing a reduction in the recurrence time for dwarf nova outbursts. We have produced Doppler tomograms and light curves which may be compared with observations. These tomograms are significantly different from those produced in the absence of a magnetic field on the secondary. We derive an upper limit to the magnetic moment of the secondary star in UGem of mu_2<2x10^32 A m^2. The magnetic truncation of the accretion disc produces resonance phenomena similar to those seen in the superoutbursts of SUUMa systems. While these have not been observed for systems like UGem, observations of the SUUMa systems provide us with a useful diagnostic of the disc-field interaction. We are able to place an upper limit on the magnetic moment of the secondary in ZCha of mu_2<1x10^30 A m^2.

Sensitivity of resolved and parameterized surface drag to changes in resolution and parameterization

The relative contribution of resolved and parameterized surface drag towards balancing the atmospheric angular momentum flux convergence (AMFC), and their sensitivity to horizontal resolution and parameterization, are investigated in an atmospheric model. This sensitivity can be difficult to elucidate in free-running climate models, in which the AMFC varies with changing climatologies and, as a result, the relative contributions of surface terms balancing the AMFC also vary. While the sensitivity question has previously been addressed using short-range forecasts, we demonstrate that a nudging framework is an effective method for constraining the AMFC. The Met Office Unified Model is integrated at three horizontal resolutions ranging from 130 km (N96) to 25 km (N512) while relaxing the model’s wind and temperature fields towards the ERAinterim reanalysis within the altitude regions of maximum AMFC. This method is validated against short range forecasts and good agreement is found. These experiments are then used to assess the fidelity of the exchange between parameterized and resolved orographic torques with changes in horizontal resolution. Although the parameterized orographic torque reduces substantially with increasing horizontal resolution, there is little change in resolved orographic torque over 20N to 50N. The tendencies produced by the nudging routine indicate that the additional drag at lower horizontal resolution is excessive. When parameterized orographic blocking is removed at the coarsest of these resolutions, there is a lack of compensation, and even compensation of the opposite sense, by the boundary layer and resolved torques which is particularly pronounced over 20N to 50N. This study demonstrates that there is strong sensitivity in the behaviour of the resolved and parameterized surface drag over this region.