Relations among several nuclear and electronic density functional reactivity indexes

A set of connections among several nuclear and electronic indexes of reactivity in the framework of the conceptual Density Functional Theory by using an expansion ofthe energy functional in terms of the total number of electrons and the normal coordinates within a canonical ensemble was derived. The relations obtained provided explicit links between important quantities related to the chemical reactivity of a system. This paper particularly demonstrates that the derivative of the electronic energy with respect to the external potential of a system in its equilibrium geometry was equal to the negative of the nuclear repulsion derivative with respect to the external potential

Two complementary molecular energy decomposition schemes: the Mayer and Ziegler-Rauk methods in comparison

In the present paper we discuss and compare two different energy decomposition schemes: Mayer's Hartree-Fock energy decomposition into diatomic and monoatomic contributions [Chem. Phys. Lett. 382, 265 (2003)], and the Ziegler-Rauk dissociation energy decomposition [Inorg. Chem. 18, 1558 (1979)]. The Ziegler-Rauk scheme is based on a separation of a molecule into fragments, while Mayer's scheme can be used in the cases where a fragmentation of the system in clearly separable parts is not possible. In the Mayer scheme, the density of a free atom is deformed to give the one-atom Mulliken density that subsequently interacts to give rise to the diatomic interaction energy. We give a detailed analysis of the diatomic energy contributions in the Mayer scheme and a close look onto the one-atom Mulliken densities. The Mulliken density ρA has a single large maximum around the nuclear position of the atom A, but exhibits slightly negative values in the vicinity of neighboring atoms. The main connecting point between both analysis schemes is the electrostatic energy. Both decomposition schemes utilize the same electrostatic energy expression, but differ in how fragment densities are defined. In the Mayer scheme, the electrostatic component originates from the interaction of the Mulliken densities, while in the Ziegler-Rauk scheme, the undisturbed fragment densities interact. The values of the electrostatic energy resulting from the two schemes differ significantly but typically have the same order of magnitude. Both methods are useful and complementary since Mayer's decomposition focuses on the energy of the finally formed molecule, whereas the Ziegler-Rauk scheme describes the bond formation starting from undeformed fragment densities

Basis set superposition error effects, excited-state potential energy surface and photodynamics of thymine

En aquesta tesi he estudiat l'efecte de l'error de superposició de base (BSSE) en la planaritat d'algunes molècules. He observat que l'ús d'alguns mètodes de càlcul amb determinades funcions de base descriuen mínims d'energia no planars per les bases nitrogenades de l'ADN. He demostrat que aquests problemes es poden arreglar utilitzant el mètode Counterpoise per corregir el BSSE en els càlculs. En aquesta tesi també he estudiat la fotofísica de la timina i els resultats mostren que existeixen dos camins de relaxació des de l'estat excitat que permeten la regeneració de l'estructura inicial de forma ultraràpida.

Slopes of Avian Species-Area Relationships, Human Population Density, and Environmental Factors

There is increasing interest in how humans influence spatial patterns in biodiversity. One of the most frequently noted and marked of these patterns is the increase in species richness with area, the species–area relationship (SAR). SARs are used for a number of conservation purposes, including predicting extinction rates, setting conservation targets, and identifying biodiversity hotspots. Such applications can be improved by a detailed understanding of the factors promoting spatial variation in the slope of SARs, which is currently the subject of a vigorous debate. Moreover, very few studies have considered the anthropogenic influences on the slopes of SARs; this is particularly surprising given that in much of the world areas with high human population density are typically those with a high number of species, which generates conservation conflicts. Here we determine correlates of spatial variation in the slopes of species–area relationships, using the British avifauna as a case study. Whilst we focus on human population density, a widely used index of human activities, we also take into account (1) the rate of increase in habitat heterogeneity with increasing area, which is frequently proposed to drive SARs, (2) environmental energy availability, which may influence SARs by affecting species occupancy patterns, and (3) species richness. We consider environmental variables measured at both local (10 km × 10 km) and regional (290 km × 290 km) spatial grains, but find that the former consistently provides a better fit to the data. In our case study, the effect of species richness on the slope SARs appears to be scale dependent, being negative at local scales but positive at regional scales. In univariate tests, the slope of the SAR correlates negatively with human population density and environmental energy availability, and positively with the rate of increase in habitat heterogeneity. We conducted two sets of multiple regression analyses, with and without species richness as a predictor. When species richness is included it exerts a dominant effect, but when it is excluded temperature has the dominant effect on the slope of the SAR, and the effects of other predictors are marginal.

An interpretation of baroclinic initial value problems: results for simple basic states with nonzero interior PV gradients

In the Eady model, where the meridional potential vorticity (PV) gradient is zero, perturbation energy growth can be partitioned cleanly into three mechanisms: (i) shear instability, (ii) resonance, and (iii) the Orr mechanism. Shear instability involves two-way interaction between Rossby edge waves on the ground and lid, resonance occurs as interior PV anomalies excite the edge waves, and the Orr mechanism involves only interior PV anomalies. These mechanisms have distinct implications for the structural and temporal linear evolution of perturbations. Here, a new framework is developed in which the same mechanisms can be distinguished for growth on basic states with nonzero interior PV gradients. It is further shown that the evolution from quite general initial conditions can be accurately described (peak error in perturbation total energy typically less than 10%) by a reduced system that involves only three Rossby wave components. Two of these are counterpropagating Rossby waves—that is, generalizations of the Rossby edge waves when the interior PV gradient is nonzero—whereas the other component depends on the structure of the initial condition and its PV is advected passively with the shear flow. In the cases considered, the three-component model outperforms approximate solutions based on truncating a modal or singular vector basis.

Aggregation of parasitism risk in an aphid-parasitoid system: Effects of plant patch size and aphid density

Few studies have linked density dependence of parasitism and the tritrophic environment within which a parasitoid forages. In the non-crop plant-aphid, Centaurea nigra-Uroleucon jaceae system, mixed patterns of density-dependent parasitism by the parasitoids Aphidius funebris and Trioxys centaureae were observed in a survey of a natural population. Breakdown of density-dependent parasitism revealed that density dependence was inverse in smaller colonies but direct in large colonies (>20 aphids), suggesting there is a threshold effect in parasitoid response to aphid density. The CV2 of searching parasitoids was estimated from parasitism data using a hierarchical generalized linear model, and CV2>1 for A. funebris between plant patches, while for T. centaureae CV2>1 within plant patches. In both cases, density independent heterogeneity was more important than density-dependent heterogeneity in parasitism. Parasitism by T. centaureae increased with increasing plant patch size. Manipulation of aphid colony size and plant patch size revealed that parasitism by A. funebris was directly density dependent at the range of colony sizes tested (50-200 initial aphids), and had a strong positive relationship with plant patch size. The effects of plant patch size detected for both species indicate that the tritrophic environment provides a source of host density independent heterogeneity in parasitism, and can modify density-dependent responses. (c) 2007 Gessellschaft fur Okologie. Published by Elsevier GmbH. All rights reserved.

The joint effects of larval density and C-14-cypermethrin on the life history and population growth rate of the midge Chironomus riparius

1. Chemical effects on organisms are typically assessed using individual-level endpoints or sometimes population growth rate (PGR), but such measurements are generally made at low population densities. In contrast most natural populations are subject to density dependence and fluctuate around the environmental carrying capacity as a result of individual competition for resources. As ecotoxicology aims to make reliable population projections of chemical impacts in the field, an understanding of how high-density or resource-limited populations respond to environmental chemicals is essential. 2. Our objective was to determine the joint effects of population density and chemical stress on the life history and PGR of an important ecotoxicological indicator species, Chironomus riparius, under controlled laboratory conditions. Populations were fed the same ration but initiated at different densities and exposed to a solvent control and three concentrations of C-14-cypermethrin in a sediment-water test system for 67 days at 20 +/- 1 degreesC. 3. Density had a negative effect on all the measured life-history traits, and PGR declined with increasing density in the controls. Exposure to C-14-cypermethrin had a direct negative effect on juvenile survival, presumably within the first 24 h because the chemical rapidly dissipated from the water column. Reductions in the initial larval densities resulted in an increase in the available resources for the survivors. Subsequently, exposed populations emerged sooner and started producing offspring earlier than the controls. C-14-cypermethrin had no effect on estimated fecundity and adult body weight but interacted with density to reduce the time to first emergence and first reproduction. As a result, PGR increased with cypermethrin concentration when populations were initiated at high densities. 4. Synthesis and applications. The results showed that the effects of C-14-cypermethrin were buffered at high density, so that the joint effects of density and chemical stress on PGR were less than additive. Low levels of chemical stressors may increase carrying capacity by reducing juvenile competition for resources. More and perhaps fitter adults may be produced, similar to the effects of predators and culling; however, toxicant exposure may result in survivors that are less tolerant to changing conditions. If less than additive effects are typical in the field, standard regulatory tests carried out at low density may overestimate the effects of environmental chemicals. Further studies over a wide range of chemical stressors and organisms with contrasting life histories are needed to make general recommendations.

The influence of larval density, food availability and habitat longevity on the life history and population growth rate of the midge Chironomus riparius

Despite long-standing interest in the forms and mechanisms of density dependence, these are still imperfectly understood. However, in a constant environment an increase in density must reduce per capita resource availability, which in turn leads to reduced survival, fecundity and somatic growth rate. Here we report two population experiments examining the density dependent responses under controlled conditions of an important indicator species, Chironomus riparius. The first experiment was run for 35 weeks and was started at low density with replicate populations being fed three different rations. Increased ration reduced generation time and increased population growth rate (pgr) but had no effect on survival, fecundity and female body weight in the first generation. In the second generation there was a six-fold increase in generation time, presumably due to the greatly reduced per capita resource availability as the estimated initial densities of the second generation were 300 times greater than the first. Juvenile survival to emergence, fecundity, adult body weight and pgr declined by 90%, 75%, 35% and 99%, respectively. These large between-generation effects may have obscured the effects of the threefold variation in ration, as only survival to emergence significantly increased with ration in the second generation. These results suggest that some chironomid larvae survive a reduction in resource availability by growing more slowly. In the ephemeral habitats sometimes occupied by C. riparius, the effects of population density may depend crucially on the longevity of the environment. A second experiment was therefore performed to measure pgr from six different starting densities over an eight-week period. The relationship between pgr and density was concave, viewed from above. At densities above 16 larvae per cm(2), less than 1% of the population emerged and no offspring were produced. Under the conditions of experiment 2 - an 8-week habitat lifespan carrying capacity was estimated as 8 larvae per cm(2).

A case of four-coordinate silver(I) adopting a high energy planar geometry. Experiment and theory

The ligands PhL and MeL are obtained by condensing 2-formylpyridine with benzil dihydrazone and diacetyl dihydrazone, respectively, in 2: 1 molar proportion. With silver( I), PhL yields a double-stranded dinuclear cationic helicate 1 in which the metal is tetrahedral but MeL gives a cationic one-dimensional polymeric complex 2 where silver( I) is distorted square planar and the ligand backbone is nearly planar. In both complexes, metal: ligand ratio is 1: 1. Ab initio calculations on the ligands at the HF/6-31+G* level reveal that while PhL strongly prefers a helical conformation, MeL has a natural inclination to remain in a planar conformation. Density functional theory calculations on model silver( I) complexes show that formation of the linear polymer in the case of MeL is also an important factor in imposing the planar geometry of Ag(I) in 2.

Solvent influences on metastable polymorph lifetimes: real-time interconversions using energy dispersive X-ray diffractometry

Solvent influences on the crystallization of polymorph and hydrate forms of the nootropic drug piracetam (2-oxo-pyrrolidineacetamide) were investigated from water, methanol, 2-propanol, isobutanol, and nitromethane. Crystal growth profiles of piracetam polymorphs were constructed using time-resolved diffraction snapshots collected for each solvent system. Measurements were performed by in situ energy dispersive X-ray diffraction recorded in Station 16.4 at the synchrotron radiation source (SRS) at Daresbury Laboratory, CCLRC UK. Crystallizations from methanol, 2-propanol, isobutanol, and nitromethane progressed in a similar fashion with the initial formation of form I which then converted relatively quickly to form II with form III being generated upon further cooling. However, considerable differences were observed for the polymorphs lifetime and both the rate and temperature of conversion using the different solvents. The thermodynamically unstable form I was kinetically favored in isobutanol and nitromethane where traces of this polymorph were observed below 10 degrees C. In contrast, the transformation of form II and subsequent growth of form III were inhibited in 2-propanol and nitromethane solutions. Aqueous solutions produced hydrate forms of piracetam which are different from the reported monohydrate; this crystallization evolved through successive generation of transient structures which transformed upon exchange of intramolecular water between the liquid and crystalline phases. (c) 2007 Wiley-Liss, Inc. and the American Pharmacists Association J Pharm Sci 96:1069-1078, 2007.

Regression based D-optimality experimental design for sparse kernel density estimation

This paper derives an efficient algorithm for constructing sparse kernel density (SKD) estimates. The algorithm first selects a very small subset of significant kernels using an orthogonal forward regression (OFR) procedure based on the D-optimality experimental design criterion. The weights of the resulting sparse kernel model are then calculated using a modified multiplicative nonnegative quadratic programming algorithm. Unlike most of the SKD estimators, the proposed D-optimality regression approach is an unsupervised construction algorithm and it does not require an empirical desired response for the kernel selection task. The strength of the D-optimality OFR is owing to the fact that the algorithm automatically selects a small subset of the most significant kernels related to the largest eigenvalues of the kernel design matrix, which counts for the most energy of the kernel training data, and this also guarantees the most accurate kernel weight estimate. The proposed method is also computationally attractive, in comparison with many existing SKD construction algorithms. Extensive numerical investigation demonstrates the ability of this regression-based approach to efficiently construct a very sparse kernel density estimate with excellent test accuracy, and our results show that the proposed method compares favourably with other existing sparse methods, in terms of test accuracy, model sparsity and complexity, for constructing kernel density estimates.

Low-lying excited states and primary photoproducts of [Os-3(CO)(10)-(s-cis-L)] (L = cyclohexa-1,3-diene, buta-1,3-diene)] clusters studied by picosecond time-resolved UV/Vis and IR spectroscopy and by density functional theory

Combined picosecond transient absorption and time-resolved infrared studies were performed, aimed at characterising low-lying excited states of the cluster [Os-3(CO)(10)(s-cis-L)] (L= cyclohexa-1,3-diene, 1) and monitoring the formation of its photoproducts. Theoretical (DFT and TD-DFT) calculations on the closely related cluster with L=buta-1,3-diene (2') have revealed that the low-lying electronic transitions of these [Os-3(CO)(10)(s-cis-1,3-diene)] clusters have a predominant sigma(core)pi*(CO) character. From the lowest sigmapi* excited state, cluster 1 undergoes fast Os-Os(1,3-diene) bond cleavage (tau=3.3 ps) resulting in the formation of a coordinatively unsaturated primary photoproduct (1a) with a single CO bridge. A new insight into the structure of the transient has been obtained by DFT calculations. The cleaved Os-Os(1,3-diene) bond is bridged by the donor 1,3-diene ligand, compensating for the electron deficiency at the neighbouring Os centre. Because of the unequal distribution of the electron density in transient la, a second CO bridge is formed in 20 ps in the photoproduct [Os-3(CO)(8)(mu-CO)(2)- (cyclohexa-1,3-diene)] (1b). The latter compound, absorbing strongly around 630 nm, mainly regenerates the parent cluster with a lifetime of about 100 ns in hexane. Its structure, as suggested by the DFT calculations, again contains the 1,3-diene ligand coordinated in a bridging fashion. Photoproduct 1b can therefore be assigned as a high-energy coordination isomer of the parent cluster with all Os-Os bonds bridged.

Kinetic energy analysis of the response of the Atlantic meridional overturning circulation to CO2-forced climate change

Atmosphere–ocean general circulation models (AOGCMs) predict a weakening of the Atlantic meridional overturning circulation (AMOC) in response to anthropogenic forcing of climate, but there is a large model uncertainty in the magnitude of the predicted change. The weakening of the AMOC is generally understood to be the result of increased buoyancy input to the north Atlantic in a warmer climate, leading to reduced convection and deep water formation. Consistent with this idea, model analyses have shown empirical relationships between the AMOC and the meridional density gradient, but this link is not direct because the large-scale ocean circulation is essentially geostrophic, making currents and pressure gradients orthogonal. Analysis of the budget of kinetic energy (KE) instead of momentum has the advantage of excluding the dominant geostrophic balance. Diagnosis of the KE balance of the HadCM3 AOGCM and its low-resolution version FAMOUS shows that KE is supplied to the ocean by the wind and dissipated by viscous forces in the global mean of the steady-state control climate, and the circulation does work against the pressure-gradient force, mainly in the Southern Ocean. In the Atlantic Ocean, however, the pressure-gradient force does work on the circulation, especially in the high-latitude regions of deep water formation. During CO2-forced climate change, we demonstrate a very good temporal correlation between the AMOC strength and the rate of KE generation by the pressure-gradient force in 50–70°N of the Atlantic Ocean in each of nine contemporary AOGCMs, supporting a buoyancy-driven interpretation of AMOC changes. To account for this, we describe a conceptual model, which offers an explanation of why AOGCMs with stronger overturning in the control climate tend to have a larger weakening under CO2 increase.

Kinetic energy analysis of the response of the Atlantic meridional overturning circulation to CO2-forced climate change

Atmosphere–ocean general circulation models (AOGCMs) predict a weakening of the Atlantic meridional overturning circulation (AMOC) in response to anthropogenic forcing of climate, but there is a large model uncertainty in the magnitude of the predicted change. The weakening of the AMOC is generally understood to be the result of increased buoyancy input to the north Atlantic in a warmer climate, leading to reduced convection and deep water formation. Consistent with this idea, model analyses have shown empirical relationships between the AMOC and the meridional density gradient, but this link is not direct because the large-scale ocean circulation is essentially geostrophic, making currents and pressure gradients orthogonal. Analysis of the budget of kinetic energy (KE) instead of momentum has the advantage of excluding the dominant geostrophic balance. Diagnosis of the KE balance of the HadCM3 AOGCM and its low-resolution version FAMOUS shows that KE is supplied to the ocean by the wind and dissipated by viscous forces in the global mean of the steady-state control climate, and the circulation does work against the pressure-gradient force, mainly in the Southern Ocean. In the Atlantic Ocean, however, the pressure-gradient force does work on the circulation, especially in the high-latitude regions of deep water formation. During CO2-forced climate change, we demonstrate a very good temporal correlation between the AMOC strength and the rate of KE generation by the pressure-gradient force in 50–70°N of the Atlantic Ocean in each of nine contemporary AOGCMs, supporting a buoyancy-driven interpretation of AMOC changes. To account for this, we describe a conceptual model, which offers an explanation of why AOGCMs with stronger overturning in the control climate tend to have a larger weakening under CO2 increase

Oxidation of low density lipoprotein by iron at lysosomal pH: implications for atherosclerosis

Low density lipoprotein (LDL) has recently been shown to be oxidised by iron within the lysosomes of macrophages and this is a novel potential mechanism for LDL oxidation in atherosclerosis. Our aim was to characterise the chemical and physical changes induced in LDL by iron at lysosomal pH and to investigate the effects of iron chelators and α-tocopherol on this process. LDL was oxidised by iron at pH 4.5 and 37°C and its oxidation monitored by spectrophotometry and HPLC. LDL was oxidised effectively by FeSO4 (5-50 µM) and became highly aggregated at pH 4.5, but not at pH 7.4. Cholesteryl esters decreased and after a pronounced lag 7-ketocholesterol increased greatly. Total hydroperoxides (measured by tri-iodide assay) increased up to 24 h and then decreased only slowly. The lipid composition after 12 h at pH 4.5 and 37°C was similar to that of LDL oxidised by copper at pH 7.4 and 4°C, i.e. rich in hydroperoxides but low in oxysterols. Previously oxidised LDL aggregated rapidly and spontaneously at pH 4.5, but not at pH 7.4. Ferrous was much more effective than ferric iron at oxidising LDL when added after the oxidation was already underway. The iron chelators diethylenetriaminepentaacetic acid and, to a lesser extent, desferrioxamine inhibited LDL oxidation when added during its initial stages, but were unable to prevent LDL aggregating after it had been partially oxidised. Surprisingly, desferrioxamine increased the rate of LDL modification when added late in the oxidation process. α-Tocopherol enrichment of LDL initially increased the oxidation of LDL, but inhibited it later. The presence of oxidised and highly aggregated lipid within lysosomes has the potential to perturb the function of these organelles and to promote atherosclerosis.