Automatic generation of active coordinates for quantum dynamics calculations: application to the dynamics of benzene photochemistry

A new practical method to generate a subspace of active coordinates for quantum dynamics calculations is presented. These reduced coordinates are obtained as the normal modes of an analytical quadratic representation of the energy difference between excited and ground states within the complete active space self-consistent field method. At the Franck-Condon point, the largest negative eigenvalues of this Hessian correspond to the photoactive modes: those that reduce the energy difference and lead to the conical intersection; eigenvalues close to 0 correspond to bath modes, while modes with large positive eigenvalues are photoinactive vibrations, which increase the energy difference. The efficacy of quantum dynamics run in the subspace of the photoactive modes is illustrated with the photochemistry of benzene, where theoretical simulations are designed to assist optimal control experiments

Nuclear relaxation contribution to static and dynamic (infinite frequency approximation) nonlinear optical properties by means of electrical property expansions: Application to HF, CH4, CF4, and SF6

Electrical property derivative expressions are presented for the nuclear relaxation contribution to static and dynamic (infinite frequency approximation) nonlinear optical properties. For CF4 and SF6, as opposed to HF and CH4, a term that is quadratic in the vibrational anharmonicity (and not previously evaluated for any molecule) makes an important contribution to the static second vibrational hyperpolarizability of CF4 and SF6. A comparison between calculated and experimental values for the difference between the (anisotropic) Kerr effect and electric field induced second-harmonic generation shows that, at the Hartree-Fock level, the nuclear relaxation/infinite frequency approximation gives the correct trend (in the series CH4, CF4, SF6) but is of the order of 50% too small

The Influence of Body Condition on the Stopover Ecology of Least Sandpipers in the Lower Mississippi Alluvial Valley during Fall Migration

Many shorebirds are long-distance migrants and depend on the energy gained at stopover sites to complete migration. Competing hypotheses have described strategies used by migrating birds; the energy-selection hypothesis predicts that shorebirds attempt to maximize energy gained at stopover sites, whereas the time-selection hypothesis predicts that shorebirds attempt to minimize time spent at stopover sites. The energy- and time-selection hypotheses both predict that birds in better condition will depart sites sooner. However, numerous studies of stopover duration have found little support for this prediction, leading to the suggestion that migrating birds operate under energy and time constraints for only a small portion of the migratory season. During fall migration 2002, we tested the prediction that birds in better condition depart stopover sites sooner by examining the relationship between stopover duration and body condition for migrating Least Sandpipers (Calidris minutilla) at three stopover sites in the Lower Mississippi Alluvial Valley. We also tested the assumption made by the Lower Mississippi Alluvial Valley Migratory Bird Science Team that shorebirds stay in the Mississippi Valley for 10 d. The assumption of 10 d was used to estimate the amount of habitat required by shorebirds in the Mississippi Valley during fall migration; a period longer than 10 d would increase the estimate of the amount habitat required. We used multiple-day constancy models of apparent survival and program MARK to estimate stopover duration for 293 individually color-marked and resighted Least Sandpipers. We found that a four-day constancy interval and a site x quadratic time trend interaction term best modeled apparent survival. We found only weak support for body condition as a factor explaining length of stopover duration, which is consistent with findings from similar work. Stopover duration estimates were 4.1 d (95% CI = 2.8–6.1) for adult Least Sandpipers at Bald Knob National Wildlife Refuge, Arkansas, 6.5 d (95% CI = 4.9–8.7) for adult and 6.1 d (95% CI =4.2–9.1) for juvenile Least Sandpipers at Yazoo National Wildlife Refuge, Mississippi, and 6.9 d (95% CI = 5.5–8.7) for juvenile Least Sandpipers at Morgan Brake National Wildlife Refuge, Mississippi. Based on our estimates of stopover duration and the assumption made by the Lower Mississippi Alluvial Valley Migratory Bird Science Team, there is sufficient habitat in the lower Mississippi Valley to support shorebirds during fall migration.

Stand and within-stand factors influencing Golden-winged Warbler use of regenerating stands in the central Appalachian Mountains

The Golden-winged Warbler (Vermivora chrysoptera) is currently being considered for protected status under the U.S. Endangered Species Act. The creation of breeding habitat in the Appalachian Mountains is considered a conservation priority for this songbird, which is dependent on extensively forested landscapes with adequate availability of young forest. We modeled abundance of Golden-winged Warbler males in regenerating harvested forest stands that were 0-17 years postharvest at both mid-Appalachian and northeast Pennsylvania regional scales using stand and within-stand characteristics of 222 regenerating stands, 2010-2011. Variables that were most influential at the mid-Appalachian scale were different than those in the northeast region. Across the mid-Appalachian ecoregion, the proportion of young forest cover, i.e., shrub/scrub cover, within 1 km of regenerating stands best explained abundance of Golden-winged Warblers. Golden-winged Warbler response was best explained by a concave quadratic relationship in which abundance was highest with 5-15% land in young forest cover. We also found evidence that the amount of herbaceous cover, i.e., the amount of grasses and forbs, within a regenerating stand positively influenced abundance of Golden-winged Warblers. In northeastern Pennsylvania, where young forest cover is found in high proportions, the distance to the nearest regenerating stand best explained variation in abundance of Golden-winged Warblers. Abundance of Golden-winged Warblers was <1 male per survey when another regenerating stand was >1500 m away. When modeling within-stand features in the northeast region, many of the models were closely ranked, indicating that multiple variables likely explained Golden-winged Warbler response to within-stand conditions. Based on our findings, we have proposed several management guidelines for land managers interested in creating breeding habitat for Golden-winged Warblers using commercial timber operations. For example, we recommend when managing for Golden-winged Warblers in the central Appalachian Mountains that managers should strive for 15% young forest in a heavily forested landscape (>70% forest cover) and cluster stands within 1-2 km of other young forest habitats.

A collocation method for high frequency scattering by convex polygons

We consider the problem of scattering of a time-harmonic acoustic incident plane wave by a sound soft convex polygon. For standard boundary or finite element methods, with a piecewise polynomial approximation space, the computational cost required to achieve a prescribed level of accuracy grows linearly with respect to the frequency of the incident wave. Recently Chandler–Wilde and Langdon proposed a novel Galerkin boundary element method for this problem for which, by incorporating the products of plane wave basis functions with piecewise polynomials supported on a graded mesh into the approximation space, they were able to demonstrate that the number of degrees of freedom required to achieve a prescribed level of accuracy grows only logarithmically with respect to the frequency. Here we propose a related collocation method, using the same approximation space, for which we demonstrate via numerical experiments a convergence rate identical to that achieved with the Galerkin scheme, but with a substantially reduced computational cost.

A Galerkin boundary element method for high frequency scattering by convex polygons

In this paper we consider the problem of time-harmonic acoustic scattering in two dimensions by convex polygons. Standard boundary or finite element methods for acoustic scattering problems have a computational cost that grows at least linearly as a function of the frequency of the incident wave. Here we present a novel Galerkin boundary element method, which uses an approximation space consisting of the products of plane waves with piecewise polynomials supported on a graded mesh, with smaller elements closer to the corners of the polygon. We prove that the best approximation from the approximation space requires a number of degrees of freedom to achieve a prescribed level of accuracy that grows only logarithmically as a function of the frequency. Numerical results demonstrate the same logarithmic dependence on the frequency for the Galerkin method solution. Our boundary element method is a discretization of a well-known second kind combined-layer-potential integral equation. We provide a proof that this equation and its adjoint are well-posed and equivalent to the boundary value problem in a Sobolev space setting for general Lipschitz domains.

QUAGMIRE v1.3: a quasi-geostrophic model for investigating rotating fluids experiments

QUAGMIRE is a quasi-geostrophic numerical model for performing fast, high-resolution simulations of multi-layer rotating annulus laboratory experiments on a desktop personal computer. The model uses a hybrid finite-difference/spectral approach to numerically integrate the coupled nonlinear partial differential equations of motion in cylindrical geometry in each layer. Version 1.3 implements the special case of two fluid layers of equal resting depths. The flow is forced either by a differentially rotating lid, or by relaxation to specified streamfunction or potential vorticity fields, or both. Dissipation is achieved through Ekman layer pumping and suction at the horizontal boundaries, including the internal interface. The effects of weak interfacial tension are included, as well as the linear topographic beta-effect and the quadratic centripetal beta-effect. Stochastic forcing may optionally be activated, to represent approximately the effects of random unresolved features. A leapfrog time stepping scheme is used, with a Robert filter. Flows simulated by the model agree well with those observed in the corresponding laboratory experiments.

Voronoi, Delaunay, and Block-Structured Mesh Refinement for Solution of the Shallow-Water Equations on the Sphere

Alternative meshes of the sphere and adaptive mesh refinement could be immensely beneficial for weather and climate forecasts, but it is not clear how mesh refinement should be achieved. A finite-volume model that solves the shallow-water equations on any mesh of the surface of the sphere is presented. The accuracy and cost effectiveness of four quasi-uniform meshes of the sphere are compared: a cubed sphere, reduced latitude–longitude, hexagonal–icosahedral, and triangular–icosahedral. On some standard shallow-water tests, the hexagonal–icosahedral mesh performs best and the reduced latitude–longitude mesh performs well only when the flow is aligned with the mesh. The inclusion of a refined mesh over a disc-shaped region is achieved using either gradual Delaunay, gradual Voronoi, or abrupt 2:1 block-structured refinement. These refined regions can actually degrade global accuracy, presumably because of changes in wave dispersion where the mesh is highly nonuniform. However, using gradual refinement to resolve a mountain in an otherwise coarse mesh can improve accuracy for the same cost. The model prognostic variables are height and momentum collocated at cell centers, and (to remove grid-scale oscillations of the A grid) the mass flux between cells is advanced from the old momentum using the momentum equation. Quadratic and upwind biased cubic differencing methods are used as explicit corrections to a fast implicit solution that uses linear differencing.

Mersenne numbers

These notes have been issued on a small scale in 1983 and 1987 and on request at other times. This issue follows two items of news. First, WaIter Colquitt and Luther Welsh found the 'missed' Mersenne prime M110503 and advanced the frontier of complete Mp-testing to 139,267. In so doing, they terminated Slowinski's significant string of four consecutive Mersenne primes. Secondly, a team of five established a non-Mersenne number as the largest known prime. This result terminated the 1952-89 reign of Mersenne primes. All the original Mersenne numbers with p < 258 were factorised some time ago. The Sandia Laboratories team of Davis, Holdridge & Simmons with some little assistance from a CRAY machine cracked M211 in 1983 and M251 in 1984. They contributed their results to the 'Cunningham Project', care of Sam Wagstaff. That project is now moving apace thanks to developments in technology, factorisation and primality testing. New levels of computer power and new computer architectures motivated by the open-ended promise of parallelism are now available. Once again, the suppliers may be offering free buildings with the computer. However, the Sandia '84 CRAY-l implementation of the quadratic-sieve method is now outpowered by the number-field sieve technique. This is deployed on either purpose-built hardware or large syndicates, even distributed world-wide, of collaborating standard processors. New factorisation techniques of both special and general applicability have been defined and deployed. The elliptic-curve method finds large factors with helpful properties while the number-field sieve approach is breaking down composites with over one hundred digits. The material is updated on an occasional basis to follow the latest developments in primality-testing large Mp and factorising smaller Mp; all dates derive from the published literature or referenced private communications. Minor corrections, additions and changes merely advance the issue number after the decimal point. The reader is invited to report any errors and omissions that have escaped the proof-reading, to answer the unresolved questions noted and to suggest additional material associated with this subject.

A wavenumber independent boundary element method for an acoustic scattering problem

In this paper we consider the impedance boundary value problem for the Helmholtz equation in a half-plane with piecewise constant boundary data, a problem which models, for example, outdoor sound propagation over inhomogeneous. at terrain. To achieve good approximation at high frequencies with a relatively low number of degrees of freedom, we propose a novel Galerkin boundary element method, using a graded mesh with smaller elements adjacent to discontinuities in impedance and a special set of basis functions so that, on each element, the approximation space contains polynomials ( of degree.) multiplied by traces of plane waves on the boundary. We prove stability and convergence and show that the error in computing the total acoustic field is O( N-(v+1) log(1/2) N), where the number of degrees of freedom is proportional to N logN. This error estimate is independent of the wavenumber, and thus the number of degrees of freedom required to achieve a prescribed level of accuracy does not increase as the wavenumber tends to infinity.

A high-wavenumber boundary-element method for an acoustic scattering problem

In this paper we show stability and convergence for a novel Galerkin boundary element method approach to the impedance boundary value problem for the Helmholtz equation in a half-plane with piecewise constant boundary data. This problem models, for example, outdoor sound propagation over inhomogeneous flat terrain. To achieve a good approximation with a relatively low number of degrees of freedom we employ a graded mesh with smaller elements adjacent to discontinuities in impedance, and a special set of basis functions for the Galerkin method so that, on each element, the approximation space consists of polynomials (of degree $\nu$) multiplied by traces of plane waves on the boundary. In the case where the impedance is constant outside an interval $[a,b]$, which only requires the discretization of $[a,b]$, we show theoretically and experimentally that the $L_2$ error in computing the acoustic field on $[a,b]$ is ${\cal O}(\log^{\nu+3/2}|k(b-a)| M^{-(\nu+1)})$, where $M$ is the number of degrees of freedom and $k$ is the wavenumber. This indicates that the proposed method is especially commendable for large intervals or a high wavenumber. In a final section we sketch how the same methodology extends to more general scattering problems.

Evaluating a new LISFLOOD-FP formulation with data from the summer 2007 floods in Tewkesbury, UK

Two-dimensional flood inundation modelling is a widely used tool to aid flood risk management. In urban areas, where asset value and population density are greatest, the model spatial resolution required to represent flows through a typical street network (i.e. < 10m) often results in impractical computational cost at the whole city scale. Explicit diffusive storage cell models become very inefficient at such high resolutions, relative to shallow water models, because the stable time step in such schemes scales as a quadratic of resolution. This paper presents the calibration and evaluation of a recently developed new formulation of the LISFLOOD-FP model, where stability is controlled by the Courant–Freidrichs–Levy condition for the shallow water equations, such that, the stable time step instead scales linearly with resolution. The case study used is based on observations during the summer 2007 floods in Tewkesbury, UK. Aerial photography is available for model evaluation on three separate days from the 24th to the 31st of July. The model covered a 3.6 km by 2 km domain and was calibrated using gauge data from high flows during the previous month. The new formulation was benchmarked against the original version of the model at 20 m and 40 m resolutions, demonstrating equally accurate performance given the available validation data but at 67x faster computation time. The July event was then simulated at the 2 m resolution of the available airborne LiDAR DEM. This resulted in a significantly more accurate simulation of the drying dynamics compared to that simulated by the coarse resolution models, although estimates of peak inundation depth were similar.

Anharmonic force field of acetylene

The harmonic and anharmonic force field of acetylene has been determined in a least-squares calculation from recently determined data on the spectroscopic constants of various isotopic species (including the vibrational l-doubling constant). A general quadratic and cubic force field was used, but a constrained quartic force field containing only 8 of the 23 possible quartic constants. The results are discussed and compared with earlier work.

The anharmonic force field of HCN and HCP

The quadratic, cubic, and quartic force field of HCN has been calculated by a least squares refinement to fit the most recent observed data on the vibration-rotation constants of HCN, DCN and H13CN. All of the observed parameters are fitted within their standard errors of observation. The corresponding parameters for other isotopic species are calculated. For HCP and DCP the more limited data available have been fitted to an anharmonic force field using constraints based on comparison with HCN. Using this force field the zero-point rotational constants B0 have been corrected to obtain the equilibrium constants Be, and hence the equilibrium structure has been determined to be re(CH) = 1•0692(7)A, and re(CP) = 1•5398(2)A.

A potential energy surface for the ground state of formaldehyde, H2CO(1A1)

A model potential energy function for the ground state of H2CO has been derived which covers the whole space of the six internal coordinates. This potential reproduces the experimental energy, geometry and quadratic force field of formaldehyde, and dissociates correctly to all possible atom, diatom and triatom fragments. Thus there are good reasons for believing it to be close to the true potential energy surface except in regions where both hydrogen atoms are close to the oxygen. It leads to the prediction that there should be a metastable singlet hydroxycarbene HCOH which has a planar trans structure and an energy of 2•31 eV above that of equilibrium formaldehyde. The reaction path for dissociation into H2 + CO is predicted to pass through a low symmetry transition state with an activation energy of 4•8 eV. Both of these predictions are in good agreement with recently published ab initio calculations.