Preconditioning and iterative solution of symmetric indefinite linear systems arising from interior point methods for linear programming

We study the preconditioning of symmetric indefinite linear systems of equations that arise in interior point solution of linear optimization problems. The preconditioning method that we study exploits the block structure of the augmented matrix to design a similar block structure preconditioner to improve the spectral properties of the resulting preconditioned matrix so as to improve the convergence rate of the iterative solution of the system. We also propose a two-phase algorithm that takes advantage of the spectral properties of the transformed matrix to solve for the Newton directions in the interior-point method. Numerical experiments have been performed on some LP test problems in the NETLIB suite to demonstrate the potential of the preconditioning method discussed.

C-H···O H-bonded complexes: how does basis set superposition error change their potential-energy surfaces?

Geometries, vibrational frequencies, and interaction energies of the CNH⋯O3 and HCCH⋯O3 complexes are calculated in a counterpoise-corrected (CP-corrected) potential-energy surface (PES) that corrects for the basis set superposition error (BSSE). Ab initio calculations are performed at the Hartree-Fock (HF) and second-order Møller-Plesset (MP2) levels, using the 6-31G(d,p) and D95++(d,p) basis sets. Interaction energies are presented including corrections for zero-point vibrational energy (ZPVE) and thermal correction to enthalpy at 298 K. The CP-corrected and conventional PES are compared; the unconnected PES obtained using the larger basis set including diffuse functions exhibits a double well shape, whereas use of the 6-31G(d,p) basis set leads to a flat single-well profile. The CP-corrected PES has always a multiple-well shape. In particular, it is shown that the CP-corrected PES using the smaller basis set is qualitatively analogous to that obtained with the larger basis sets, so the CP method becomes useful to correctly describe large systems, where the use of small basis sets may be necessary

How does basis set superposition error change the potential surfaces for hydrogen-bonded dimers?

We describe a simple method to automate the geometric optimization of molecular orbital calculations of supermolecules on potential surfaces that are corrected for basis set superposition error using the counterpoise (CP) method. This method is applied to the H-bonding complexes HF/HCN, HF/H2O, and HCCH/H2O using the 6-31G(d,p) and D95 + + (d,p) basis sets at both the Hartree-Fock and second-order Møller-Plesset levels. We report the interaction energies, geometries, and vibrational frequencies of these complexes on the CP-optimized surfaces; and compare them with similar values calculated using traditional methods, including the (more traditional) single point CP correction. Upon optimization on the CP-corrected surface, the interaction energies become more negative (before vibrational corrections) and the H-bonding stretching vibrations decrease in all cases. The extent of the effects vary from extremely small to quite large depending on the complex and the calculational method. The relative magnitudes of the vibrational corrections cannot be predicted from the H-bond stretching frequencies alone

Linear amplification of marginally neutral baroclinic waves

Baroclinic wave development is investigated for unstable parallel shear flows in the limit of vanishing normal-mode growth rate. This development is described in terms of the propagation and interaction mechanisms of two coherent structures, called counter-propagating Rossby waves (CRWs). It is shown that, in this limit of vanishing normal-mode growth rate, arbitrary initial conditions produce sustained linear amplification of the marginally neutral normal mode (mNM). This linear excitation of the mNM is subsequently interpreted in terms of a resonance phenomenon. Moreover, while the mathematical character of the normal-mode problem changes abruptly as the bifurcation point in the dispersion diagram is encountered and crossed, it is shown that from an initial-value viewpoint, this transition is smooth. Consequently, the resonance interpretation remains relevant (albeit for a finite time) for wavenumbers slightly different from the ones defining cut-off points. The results are further applied to a two-layer version of the classic Eady model in which the upper rigid lid has been replaced by a simple stratosphere.

Comparison of EDTA and EDDS as potential soil amendments for enhanced phytoextraction of heavy metals

Phytoextraction has been proposed as an alternative remediation technology for soils polluted with heavy metals or radionuclides, but is generally conceived as too slow working. Enhancing the accumulation of trace pollutants in harvestable plant tissues is a prerequisite for the technology to be practically applicable. The chelating aminopolycarboxylic acid, ethylene diamine tetraacetate (EDTA), has been found to enhance shoot accumulation of heavy metals. However, the use of EDTA in phytoextraction may not be suitable due to its high environmental persistence, which may lead to groundwater contamination. This paper aims to assess whether ethylene diamine disuccinate (EDDS), a biodegradable chelator, can be used for enhanced phytoextraction purposes. A laboratory experiment was conducted to examine mobilisation of Cd, Cu, Cr, Ni, Pb and Zn into the soil solution upon application of EDTA or EDDS. The longevity of the induced mobilisation was monitored for a period of 40 days after application. Estimated effect half lives ranged between 3.8 and 7.5 days for EDDS, depending on the applied dose. The minimum observed effect half life of EDTA was 36 days, while for the highest applied dose no decrease was observed throughout the 40 day period of the mobilisation experiment. Performance of EDTA and EDDS for phytoextraction was evaluated by application to Helianthus annuus. Two other potential chelators, known for their biodegradability in comparison to EDTA, were tested in the plant experiment: nitrilo acetic acid (NTA) and citric acid. Uptake of heavy metals was higher in EDDS-treated pots than in EDTA-treated pots. The effects were still considered insufficiently high to consider efficient remediation. This may be partly due to the choice of timing for application of the soil amendment. Fixing the time of application at an earlier point before harvest may yield better results. NTA and citric acid induced no significant effects on heavy metal uptake. (C) 2004 Elsevier Ltd. All rights reserved.

Wetland restoration within agricultural watersheds: Balancing water quality protection with habitat conservation

Peat wetlands that have been restored from agricultural Land have the potential to act as Long term sources of phosphorus (P) and, therefore have to potenital to accelerate freshwater eutrophication. During a two-year study the water table in a eutrophic fen peat that was managed by pump drainage fluctuated annually between +20 cm and -60 cm relative to ground Level. This precise management was facilitated by the high hydraulic conductivity (K) of the humified peat (1.1 x 10(-5) m s(-1)) below around 60 cm depth. However, during one week of intermittent pumping, as much as 50 g ha(-1) dissolved P entered the pumped ditch. Summer. rainfall events and autumn reflooding also triggered P losses. The P Losses were attributed to the low P sorption capacity (217 mg kg(-1)) of the saturated peat below 60 cm, combined with its high K and the reductive dissolution of Fe bound P.

A potential energy surface for the ground state of acetylene, C2H2

A potential energy function has been derived for the ground state surface of C2H2 as a many-body expansion. The 2- and 3-body terms have been obtained by preliminary investigation of the ground state surfaces of CH2( 3B1) and C2H( 2Σ+). A 4-body term has been derived which reproduces the energy, geometry and harmonic force field of C2H2. The potential has a secondary minimum corresponding to the vinylidene structure and the geometry and energy of this are in close agreement with predictions from ab initio calculations. The saddle point for the HCCH-H2CC rearrangement is predicted to lie 2•530 eV above the acetylene minimum.

In situ degradability of dry matter and neutral-detergent fibre of thorn scrubland forage consumed by goats in the semi-arid region of north Mexico

Three goats provided with oesophageal and ruminal cannulae were used to determine variations in dry matter (DM) and neutral-detergent fibre (NDF) degradability of the forage consumed when grazing thorn scrubland in the semi-arid region of north Mexico, during two consecutive dry and wet periods. Ingesta samples were incubated intraruminally, the data were fitted to the exponential equation P = a + b (1-e(-ct)) and statistically analysed using a randomized-block design. Organic matter and crude protein (CP) contents were higher (P < 0.05) in the wet seasons. Values of NDF were similar in dry and wet season of both years whereas higher numerical values of acid-detergent fibre (ADF), lignin and cellulose were registered in the dry seasons. DM and NDF degradabilities after 24 and 48 h of ruminal incubation were higher (P < 0.05) in the wet seasons. Higher values (P < 0.05) in DM and NDF bag losses at zero time (A fraction) were registered in the two wet seasons. The insoluble but fermentable DM and NDF (B fractions) were higher (P < 0.05) in the 1999 wet season and variable in the rest of the studied period. Numerically higher values of DM and NDF c fraction were found in wet periods, whereas DM and NDF potential degradabilities were higher (P < 0.05) in the wet season in 1999 and similar across seasons in 2000. Lowest (P < 0.05) contents of CP in grazed forage, DM and NDF degradabilities after 48 h of ruminal incubation, and A, and B, and c fractions were observed in the dry seasons. Thus, these results may be related to both the lower feeding value of forage consumed by the animals and lower performance of livestock during this period. Then, the DM and NDF degradability after 48 h, together with the insoluble but fermentable matter and the c fraction permit the nutritive value of the forage consumed by grazing goats to be accurately described.

Full-dimensional vibrational calculations for H5O2+ using an ab initio potential energy surface

We report quantum diffusion Monte Carlo (DMC) and variational calculations in full dimensionality for selected vibrational states of H5O2+ using a new ab initio potential energy surface [X. Huang, B. Braams, and J. M. Bowman, J. Chem. Phys. 122, 044308 (2005)]. The energy and properties of the zero-point state are focused on in the rigorous DMC calculations. OH-stretch fundamentals are also calculated using "fixed-node" DMC calculations and variationally using two versions of the code MULTIMODE. These results are compared with infrared multiphoton dissociation measurements of Yeh [L. I. Yeh, M. Okumura, J. D. Myers, J. M. Price, and Y. T. Lee, J. Chem. Phys. 91, 7319 (1989)]. Some preliminary results for the energies of several modes of the shared hydrogen are also reported.

Full-dimensional quantum calculations of ground-state tunneling splitting of malonaldehyde using an accurate ab initio potential energy surface

Quantum calculations of the ground vibrational state tunneling splitting of H-atom and D-atom transfer in malonaldehyde are performed on a full-dimensional ab initio potential energy surface (PES). The PES is a fit to 11 147 near basis-set-limit frozen-core CCSD(T) electronic energies. This surface properly describes the invariance of the potential with respect to all permutations of identical atoms. The saddle-point barrier for the H-atom transfer on the PES is 4.1 kcal/mol, in excellent agreement with the reported ab initio value. Model one-dimensional and "exact" full-dimensional calculations of the splitting for H- and D-atom transfer are done using this PES. The tunneling splittings in full dimensionality are calculated using the unbiased "fixed-node" diffusion Monte Carlo (DMC) method in Cartesian and saddle-point normal coordinates. The ground-state tunneling splitting is found to be 21.6 cm(-1) in Cartesian coordinates and 22.6 cm(-1) in normal coordinates, with an uncertainty of 2-3 cm(-1). This splitting is also calculated based on a model which makes use of the exact single-well zero-point energy (ZPE) obtained with the MULTIMODE code and DMC ZPE and this calculation gives a tunneling splitting of 21-22 cm(-1). The corresponding computed splittings for the D-atom transfer are 3.0, 3.1, and 2-3 cm(-1). These calculated tunneling splittings agree with each other to within less than the standard uncertainties obtained with the DMC method used, which are between 2 and 3 cm(-1), and agree well with the experimental values of 21.6 and 2.9 cm(-1) for the H and D transfer, respectively. (C) 2008 American Institute of Physics.

MULTIMODE quantum calculations of intramolecular vibrational energies of the water dimer and trimer using ab initio-based potential energy surfaces

We report vibrational configuration interaction calculations of the monomer fundamentals of (H2O)(2), (D2O)(2), (H2O)(3), and (D2O)(3) using the code MULTIMODE and full dimensional ab initio-based global potential energies surfaces (PESs). For the dimer the HBB PES [Huang , J. Chem. Phys 128, 034312 (2008)] is used and for the trimer a new PES, reported here, is used. The salient properties of the new trimer PES are presented and compared to previous single-point calculations and the vibrational energies are compared with experiments. (C) 2008 American Institute of Physics.

Case study: object and ground point separation on a river bank using airborne based LIDAR data

In the past decade, airborne based LIght Detection And Ranging (LIDAR) has been recognised by both the commercial and public sectors as a reliable and accurate source for land surveying in environmental, engineering and civil applications. Commonly, the first task to investigate LIDAR point clouds is to separate ground and object points. Skewness Balancing has been proven to be an efficient non-parametric unsupervised classification algorithm to address this challenge. Initially developed for moderate terrain, this algorithm needs to be adapted to handle sloped terrain. This paper addresses the difficulty of object and ground point separation in LIDAR data in hilly terrain. A case study on a diverse LIDAR data set in terms of data provider, resolution and LIDAR echo has been carried out. Several sites in urban and rural areas with man-made structure and vegetation in moderate and hilly terrain have been investigated and three categories have been identified. A deeper investigation on an urban scene with a river bank has been selected to extend the existing algorithm. The results show that an iterative use of Skewness Balancing is suitable for sloped terrain.

Bio-energy retains its mitigation potential under elevated CO2

Background If biofuels are to be a viable substitute for fossil fuels, it is essential that they retain their potential to mitigate climate change under future atmospheric conditions. Elevated atmospheric CO2 concentration [CO2] stimulates plant biomass production; however, the beneficial effects of increased production may be offset by higher energy costs in crop management. Methodology/Main findings We maintained full size poplar short rotation coppice (SRC) systems under both current ambient and future elevated [CO2] (550 ppm) and estimated their net energy and greenhouse gas balance. We show that a poplar SRC system is energy efficient and produces more energy than required for coppice management. Even more, elevated [CO2] will increase the net energy production and greenhouse gas balance of a SRC system with 18%. Managing the trees in shorter rotation cycles (i.e. 2 year cycles instead of 3 year cycles) will further enhance the benefits from elevated [CO2] on both the net energy and greenhouse gas balance. Conclusions/significance Adapting coppice management to the future atmospheric [CO2] is necessary to fully benefit from the climate mitigation potential of bio-energy systems. Further, a future increase in potential biomass production due to elevated [CO2] outweighs the increased production costs resulting in a northward extension of the area where SRC is greenhouse gas neutral. Currently, the main part of the European terrestrial carbon sink is found in forest biomass and attributed to harvesting less than the annual growth in wood. Because SRC is intensively managed, with a higher turnover in wood production than conventional forest, northward expansion of SRC is likely to erode the European terrestrial carbon sink.

Threshold-free object and ground point separation in LIDAR data

Light Detection And Ranging (LIDAR) is an important modality in terrain and land surveying for many environmental, engineering and civil applications. This paper presents the framework for a recently developed unsupervised classification algorithm called Skewness Balancing for object and ground point separation in airborne LIDAR data. The main advantages of the algorithm are threshold-freedom and independence from LIDAR data format and resolution, while preserving object and terrain details. The framework for Skewness Balancing has been built in this contribution with a prediction model in which unknown LIDAR tiles can be categorised as “hilly” or “moderate” terrains. Accuracy assessment of the model is carried out using cross-validation with an overall accuracy of 95%. An extension to the algorithm is developed to address the overclassification issue for hilly terrain. For moderate terrain, the results show that from the classified tiles detached objects (buildings and vegetation) and attached objects (bridges and motorway junctions) are separated from bare earth (ground, roads and yards) which makes Skewness Balancing ideal to be integrated into geographic information system (GIS) software packages.

Potential energy surface and moleculardynamics of MbNO: existence of an unsuspected FeON minimum

Ligands such as CO, O2, or NO are involved in the biological function of myoglobin. Here we investigate the energetics and dynamics of NO interacting with the Fe(II) heme group in native myoglobin using ab initio and molecular dynamics simulations. At the global minimum of the ab initio potential energy surface (PES), the binding energy of 23.4 kcal/mol and the Fe-NO structure compare well with the experimental results. Interestingly, the PES is found to exhibit two minima: There exists a metastable, linear Fe-O-N minimum in addition to the known, bent Fe-N-O global minimum conformation. Moreover, the T-shaped configuration is found to be a saddle point, in contrast to the corresponding minimum for NO interacting with Fe(III). To use the ab initio results for finite temperature molecular dynamics simulations, an analytical function was fitted to represent the Fe-NO interaction. The simulations show that the secondary minimum is dynamically stable up to 250 K and has a lifetime of several hundred picoseconds at 300 K. The difference in the topology of the heme-NO PES from that assumed previously (one deep, single Fe-NO minimum) suggests that it is important to use the full PES for a quantitative understanding of this system. Why the metastable state has not been observed in the many spectroscopic studies of myoglobin interacting with NO is discussed, and possible approaches to finding it are outlined.