Disentangling multi-level systems: averaging, correlations and memory

We consider two weakly coupled systems and adopt a perturbative approach based on the Ruelle response theory to study their interaction. We propose a systematic way of parameterizing the effect of the coupling as a function of only the variables of a system of interest. Our focus is on describing the impacts of the coupling on the long term statistics rather than on the finite-time behavior. By direct calculation, we find that, at first order, the coupling can be surrogated by adding a deterministic perturbation to the autonomous dynamics of the system of interest. At second order, there are additionally two separate and very different contributions. One is a term taking into account the second-order contributions of the fluctuations in the coupling, which can be parameterized as a stochastic forcing with given spectral properties. The other one is a memory term, coupling the system of interest to its previous history, through the correlations of the second system. If these correlations are known, this effect can be implemented as a perturbation with memory on the single system. In order to treat this case, we present an extension to Ruelle's response theory able to deal with integral operators. We discuss our results in the context of other methods previously proposed for disentangling the dynamics of two coupled systems. We emphasize that our results do not rely on assuming a time scale separation, and, if such a separation exists, can be used equally well to study the statistics of the slow variables and that of the fast variables. By recursively applying the technique proposed here, we can treat the general case of multi-level systems.

Approximate Gauss-Newton methods for nonlinear least squares problems

The Gauss–Newton algorithm is an iterative method regularly used for solving nonlinear least squares problems. It is particularly well suited to the treatment of very large scale variational data assimilation problems that arise in atmosphere and ocean forecasting. The procedure consists of a sequence of linear least squares approximations to the nonlinear problem, each of which is solved by an “inner” direct or iterative process. In comparison with Newton’s method and its variants, the algorithm is attractive because it does not require the evaluation of second-order derivatives in the Hessian of the objective function. In practice the exact Gauss–Newton method is too expensive to apply operationally in meteorological forecasting, and various approximations are made in order to reduce computational costs and to solve the problems in real time. Here we investigate the effects on the convergence of the Gauss–Newton method of two types of approximation used commonly in data assimilation. First, we examine “truncated” Gauss–Newton methods where the inner linear least squares problem is not solved exactly, and second, we examine “perturbed” Gauss–Newton methods where the true linearized inner problem is approximated by a simplified, or perturbed, linear least squares problem. We give conditions ensuring that the truncated and perturbed Gauss–Newton methods converge and also derive rates of convergence for the iterations. The results are illustrated by a simple numerical example. A practical application to the problem of data assimilation in a typical meteorological system is presented.

Relationships between lightning and properties of convective cloud clusters

Satellite observations of convective system properties and lightning flash rate are used to investigate the ability of potential lightning parameterizations to capture both the dominant land-ocean contrast in lightning occurrence and regional differences between Africa, the Amazon and the islands of the maritime continent. As found in previous studies, the radar storm height is tightly correlated with the lightning flash rate. A roughly second order power-law fit to the mean radar echo top height above the 0C isotherm is shown to capture both regional and land-ocean contrasts in lightning occurrence and flash rate using a single set of parameters. Recent developments should soon make it possible to implement a parameterization of this kind in global models. Parameterizations based on cloud top height, convective rain rate and convective rain fraction all require the use of separate fits over land and ocean and fail to capture observed differences between continental regions.

Application of the Lighthill-Ford theory of spontaneous imbalance to clear-air turbulence forecasting

A new method of clear-air turbulence (CAT) forecasting based on the Lighthill–Ford theory of spontaneous imbalance and emission of inertia–gravity waves has been derived and applied on episodic and seasonal time scales. A scale analysis of this shallow-water theory for midlatitude synoptic-scale flows identifies advection of relative vorticity as the leading-order source term. Examination of leading- and second-order terms elucidates previous, more empirically inspired CAT forecast diagnostics. Application of the Lighthill–Ford theory to the Upper Mississippi and Ohio Valleys CAT outbreak of 9 March 2006 results in good agreement with pilot reports of turbulence. Application of Lighthill–Ford theory to CAT forecasting for the 3 November 2005–26 March 2006 period using 1-h forecasts of the Rapid Update Cycle (RUC) 2 1500 UTC model run leads to superior forecasts compared to the current operational version of the Graphical Turbulence Guidance (GTG1) algorithm, the most skillful operational CAT forecasting method in existence. The results suggest that major improvements in CAT forecasting could result if the methods presented herein become operational.

The acoustic signature of fluid flow in complex porous media

Effective medium approximations for the frequency-dependent and complex-valued effective stiffness tensors of cracked/ porous rocks with multiple solid constituents are developed on the basis of the T-matrix approach (based on integral equation methods for quasi-static composites), the elastic - viscoelastic correspondence principle, and a unified treatment of the local and global flow mechanisms, which is consistent with the principle of fluid mass conservation. The main advantage of using the T-matrix approach, rather than the first-order approach of Eshelby or the second-order approach of Hudson, is that it produces physically plausible results even when the volume concentrations of inclusions or cavities are no longer small. The new formulae, which operates with an arbitrary homogeneous (anisotropic) reference medium and contains terms of all order in the volume concentrations of solid particles and communicating cavities, take explicitly account of inclusion shape and spatial distribution independently. We show analytically that an expansion of the T-matrix formulae to first order in the volume concentration of cavities (in agreement with the dilute estimate of Eshelby) has the correct dependence on the properties of the saturating fluid, in the sense that it is consistent with the Brown-Korringa relation, when the frequency is sufficiently low. We present numerical results for the (anisotropic) effective viscoelastic properties of a cracked permeable medium with finite storage porosity, indicating that the complete T-matrix formulae (including the higher-order terms) are generally consistent with the Brown-Korringa relation, at least if we assume the spatial distribution of cavities to be the same for all cavity pairs. We have found an efficient way to treat statistical correlations in the shapes and orientations of the communicating cavities, and also obtained a reasonable match between theoretical predictions (based on a dual porosity model for quartz-clay mixtures, involving relatively flat clay-related pores and more rounded quartz-related pores) and laboratory results for the ultrasonic velocity and attenuation spectra of a suite of typical reservoir rocks. (C) 2003 Elsevier B.V. All rights reserved.

Approximations to the scattering of water waves by steep topography

A new method is developed for approximating the scattering of linear surface gravity waves on water of varying quiescent depth in two dimensions. A conformal mapping of the fluid domain onto a uniform rectangular strip transforms steep and discontinuous bed profiles into relatively slowly varying, smooth functions in the transformed free-surface condition. By analogy with the mild-slope approach used extensively in unmapped domains, an approximate solution of the transformed problem is sought in the form of a modulated propagating wave which is determined by solving a second-order ordinary differential equation. This can be achieved numerically, but an analytic solution in the form of a rapidly convergent infinite series is also derived and provides simple explicit formulae for the scattered wave amplitudes. Small-amplitude and slow variations in the bedform that are excluded from the mapping procedure are incorporated in the approximation by a straightforward extension of the theory. The error incurred in using the method is established by means of a rigorous numerical investigation and it is found that remarkably accurate estimates of the scattered wave amplitudes are given for a wide range of bedforms and frequencies.

Vibrational intensities in methane

The absorption intensities of the two infra-red active vibrations in methane have been obtained from a perturbation calculation on the equilibrium wave functions derived in the preceding paper. The perturbation field is the change in the potential field due to the nuclei which results from moving the nuclei in the vibrational coordinate concerned, and a simplified form of second order perturbation theory, developed by Pople and Schofield, is used for the calculation. The main approximation involved is the neglect of f and higher harmonics in the spherical harmonic expansion of the nuclear field. The resulting dipole moment derivatives are approximately three times larger than the experimental values, but they show qualitative features and sign relationships which are significant.

Quartic anharmonic resonances in acetylene molecules

Formulas are derived for the quartic anharmonic resonance coefficients observed to be important between C–H stretching and the combination of one quantum of C≡C stretching and two quanta of H–C≡C bending in a number of acetylene molecules. Examples of this resonance are ν3 with ν2+ν4+ν5 in 12C2H2, ν1 with ν2+2ν5 in 13C2H2, and ν1 with ν2+2ν4 in monofluoroacetylene and monochloroacetylene. The coefficients characterizing the resonances in these examples, which we denote K3,245, K1,255, and K1,244, arise from cubic and quartic terms in the anharmonic force field, in the normal coordinate representation, through second order and first order perturbation treatments respectively, where the second order resonances are calculated by a Van Vleck resonance formalism. The experimentally determined values of these coefficients are compared with values calculated from model anharmonic force fields.

The v2, v5 Raman band system of 13CH3F

The anisotropic and isotropic components of the ν2, ν5 rotation-vibrational Raman bands of 13CH3F were obtained separately. The two upper states are coupled by a strong second-order Coriolis resonance. The anisotropic spectrum was analyzed by means of a program system due to R. Escribano. A contour simulation and a least-squares fit of 233 assigned transitions yielded values for ν5, ΔA5, ΔA2, and Aζ5a, 5b(z). The 13C shifts of ν2 and ν5 were obtained from the isotropic spectrum.

Simultaneous analysis of v1, v4, 2v2, v2+v5, and 2v5 bands of 12CH3F

The Fourier-transform spectrum of CH3F from 2800 to 3100 cm−1, obtained by Guelachvili in Orsay at a resolution of about 0.003 cm−1, was analyzed. The effective Hamiltonian used contained all symmetry allowed interactions up to second order in the Amat-Nielsen classification, together with selected third-order terms, amongst the set of nine vibrational basis functions represented by the states ν1(A1), ν4(E), 2ν2(A1), ν2 + ν5(E), 2ν50(A1), and 2ν5±2(E). A number of strong Fermi and Coriolis resonances are involved. The vibrational Hamiltonian matrix was not factorized beyond the requirements of symmetry. A total of 59 molecular parameters were refined in a simultaneous least-squares analysis to over 1500 upper-state energy levels for J ≤ 20 with a standard deviation of 0.013 cm−1. Although the standard deviation remains an order of magnitude greater than the precision of the measurements, this work breaks new ground in the simultaneous analysis of interacting symmetric top vibrational levels, in terms of the number of interacting vibrational states and the number of parameters in the Hamiltonian.

Detection of frailty in Weibull lifetime data using Outlier tests

Heterogeneity in lifetime data may be modelled by multiplying an individual's hazard by an unobserved frailty. We test for the presence of frailty of this kind in univariate and bivariate data with Weibull distributed lifetimes, using statistics based on the ordered Cox-Snell residuals from the null model of no frailty. The form of the statistics is suggested by outlier testing in the gamma distribution. We find through simulation that the sum of the k largest or k smallest order statistics, for suitably chosen k , provides a powerful test when the frailty distribution is assumed to be gamma or positive stable, respectively. We provide recommended values of k for sample sizes up to 100 and simple formulae for estimated critical values for tests at the 5% level.

Molecular structures of trimethylchlorogermane, (CH3)(3)GeCl, and trimethylbromogermane, (CH3)(3)GeBr, obtained by gas-phase electron diffraction and theoretical calculations

The structures of trimethylchlorogermane ((CH3)(3)GeCl) and trimethylbromogermane ((CH3)(3)GeBr) have been determined by gas-phase electron diffraction (GED), augmented by the results from ab initio calculations employing second-order Moller-Plesset (MP2) level of theory and the 6-311+G(d) basis set. All the electrons were included in the correlation calculation. The results from the ab initio calculations indicated that these molecules have C-3v symmetry, and models with this symmetry were used in the electron diffraction analysis. The results for the principal distances (r(g)) and angles (angle(alpha)) from the combined GED/ab initio study of trimethylchlorogermane (with estimated 2sigma uncertainties) are: r(Ge-C) = 1.950(4) Angstrom, r(Ge-Cl) = 2.173(4) Angstrom, r(C-H) = 1.090(9) Angstrom, angleCGeC = 112.7(7)degrees, angleCGeCl = 106.0(8)degrees, angleGeCH = 107.8(12)degrees. The results for the principal distances (r(g)) and angles (angle(alpha)) from the combined GED/ab initio study of trimethylbromogermane (with estimated 2sigma uncertainties) are: r(Ge-C) = 1.952(7) Angstrom, r(Ge-Br) = 2.325(4) Angstrom, r(C-H) = 1. 140(28) Angstrom, angleCGeC = 114.2(11)degrees, angleCGeBr = 104.2(13)degrees, angleGeCH 106.9(43)degrees. Local C-3v symmetry and staggered conformation were assumed for the methyl groups.

Time-resolved gas-phase kinetic and quantum chemical studies of the reaction of silylene with oxygen

Time-resolved kinetic studies of the reaction of silylene, SiH2, generated by laser flash photolysis of phenylsilane, have been carried out to obtain rate constants for its bimolecular reaction with O-2. The reaction was studied in the gas phase over the pressure range 1-100 Torr in SF6 bath gas, at five temperatures in the range 297-600 K. The second order rate constants at 10 Torr were fitted to the Arrhenius equation: log(k/cm(3) molecule(-1) s(-1)) = (-11.08 +/- 0.04) + (1.57 +/- 0.32 kJ mol(-1))/RT ln10 The decrease in rate constant values with increasing temperature, although systematic is very small. The rate constants showed slight increases in value with pressure at each temperature, but this was scarcely beyond experimental uncertainty. From estimates of Lennard-Jones collision rates, this reaction is occurring at ca. 1 in 20 collisions, almost independent of pressure and temperature. Ab initio calculations at the G3 level backed further by multi-configurational (MC) SCF calculations, augmented by second order perturbation theory (MRMP2), support a mechanism in which the initial adduct, H2SiOO, formed in the triplet state (T), undergoes intersystem crossing to the more stable singlet state (S) prior to further low energy isomerisation processes leading, via a sequence of steps, ultimately to dissociation products of which the lowest energy pair are H2O + SiO. The decomposition of the intermediate cyclo-siladioxirane, via O-O bond fission, plays an important role in the overall process. The bottleneck for the overall process appears to be the T -> S process in H2SiOO. This process has a small spin orbit coupling matrix element, consistent with an estimate of its rate constant of 1 x 10(9) s(-1) obtained with the aid of RRKM theory. This interpretation preserves the idea that, as in its reactions in general, SiH2 initially reacts at the encounter rate with O-2. The low values for the secondary reaction barriers on the potential energy surface account for the lack of an observed pressure dependence. Some comparisons are drawn with the reactions of CH2 + O-2 and SiCl2 + O-2.

Time-resolved gas-phase kinetic and quantum chemical studies of the reaction of silylene with nitric oxide

Time-resolved kinetic studies of the reaction of silylene, SiH2, generated by laser flash photolysis of phenylsilane, have been carried out to obtain rate constants for its bimolecular reaction with NO. The reaction was studied in the gas phase over the pressure range 1-100 Torr in SF6 bath gas at five temperatures in the range 299-592 K. The second-order rate constants at 10 Torr fitted the Arrhenius equation log(k/cm(3) molecule(-1) s(-1)) = (- 11.66 +/- 0.01) + (6.20 +/- 0.10 kJ mol(-1))IRT In 10 The rate constants showed a variation with pressure of a factor of ca. 2 over the available range, almost independent of temperature. The data could not be fitted by RRKM calculations to a simple third body assisted association reaction alone. However, a mechanistic model with an additional (pressure independent) side channel gave a reasonable fit to the data. Ab initio calculations at the G3 level supported a mechanism in which the initial adduct, bent H2SiNO, can ring close to form cyclo-H2SiNO, which is partially collisionally stabilized. In addition, bent H2SiNO can undergo a low barrier isomerization reaction leading, via a sequence of steps, ultimately to dissociation products of which the lowest energy pair are NH2 + SiO. The rate controlling barrier for this latter pathway is only 16 kJ mol(-1) below the energy of SiH2 + NO. This is consistent with the kinetic findings. A particular outcome of this work is that, despite the pressure dependence and the effects of the secondary barrier (in the side reaction), the initial encounter of SiH2 with NO occurs at the collision rate. Thus, silylene can be as reactive with odd electron molecules as with many even electron species. Some comparisons are drawn with the reactions of CH2 + NO and SiCl2 + NO.

The addition reaction between silylene and ethyne: further isotope studies, pressure dependence studies, and quantum chemical calculations

Time-resolved kinetic studies of the reaction of dideutero-silylene, SiD2, generated by laser flash photolysis of phenylsilane-d(3), have been carried out to obtain rate constants for its bimolecular reaction with C2H2. The reaction was studied in the gas phase over the pressure range 1-100 Torr in SF6 bath gas, at five temperatures in the range 297-600 K. The second-order rate constants obtained by extrapolation to the high-pressure limits at each temperature fitted the Arrhenius equation log(k(infinity)/cm(3) molecule(-1) s(-1)) = (-10.05 +/- 0.05) + (3.43 +/- 0.36 kJ mol(-1))/RT ln 10. The rate constants were used to obtain a comprehensive set of isotope effects by comparison with earlier obtained rate constants for the reactions of SiH2 with C2H2 and C2D2. Additionally, pressure-dependent rate constants for the reaction of SiH2 with C2H2 in the presence of He (1-100 Tort) were obtained at 300, 399, and 613 K. Quantum chemical (ab initio) calculations of the SiC2H4 reaction system at the G3 level support the initial formation of silirene, which rapidly isomerizes to ethynylsilane as the major pathway. Reversible formation of vinylsilylene is also an important process. The calculations also indicate the involvement of several other intermediates, not previously suggested in the mechanism. RRKM calculations are in semiquantitative agreement with the pressure dependences and isotope effects suggested by the ab initio calculations, but residual discrepancies suggest the possible involvement of the minor reaction channel, SiH2 + C2H2 - SWPO + C2H4. The results are compared and contrasted with previous studies of this reaction system.