A comparative analysis by means of quantum molecular similarity measures of density distributions derived from conventional ab initio and density functional methods

A procedure based on quantum molecular similarity measures (QMSM) has been used to compare electron densities obtained from conventional ab initio and density functional methodologies at their respective optimized geometries. This method has been applied to a series of small molecules which have experimentally known properties and molecular bonds of diverse degrees of ionicity and covalency. Results show that in most cases the electron densities obtained from density functional methodologies are of a similar quality than post-Hartree-Fock generalized densities. For molecules where Hartree-Fock methodology yields erroneous results, the density functional methodology is shown to yield usually more accurate densities than those provided by the second order Møller-Plesset perturbation theory

A new all-round density functional based on spin states and SN2 barriers

We report here a new empirical density functional that is constructed based on the performance of OPBE and PBE for spin states and SN 2 reaction barriers and how these are affected by different regions of the reduced gradient expansion. In a previous study [Swart, Sol̀, and Bickelhaupt, J. Comput. Methods Sci. Eng. 9, 69 (2009)] we already reported how, by switching between OPBE and PBE, one could obtain both the good performance of OPBE for spin states and reaction barriers and that of PBE for weak interactions within one and the same (SSB-sw) functional. Here we fine tuned this functional and include a portion of the KT functional and Grimme's dispersion correction to account for π- π stacking. Our new SSB-D functional is found to be a clear improvement and functions very well for biological applications (hydrogen bonding, π -π stacking, spin-state splittings, accuracy of geometries, reaction barriers)

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

Variational calculation of vibrational linear and nonlinear optical properties

A variational approach for reliably calculating vibrational linear and nonlinear optical properties of molecules with large electrical and/or mechanical anharmonicity is introduced. This approach utilizes a self-consistent solution of the vibrational Schrödinger equation for the complete field-dependent potential-energy surface and, then, adds higher-level vibrational correlation corrections as desired. An initial application is made to static properties for three molecules of widely varying anharmonicity using the lowest-level vibrational correlation treatment (i.e., vibrational Møller-Plesset perturbation theory). Our results indicate when the conventional Bishop-Kirtman perturbation method can be expected to break down and when high-level vibrational correlation methods are likely to be required. Future improvements and extensions are discussed

Basis set and electron correlation effects on initial convergence for vibrational nonlinear optical properties of conjugated organic molecules

The level of ab initio theory which is necessary to compute reliable values for the static and dynamic (hyper)polarizabilities of three medium size π-conjugated organic nonlinear optical (NLO) molecules is investigated. With the employment of field-induced coordinates in combination with a finite field procedure, the calculations were made possible. It is stated that to obtain reasonable values for the various individual contributions to the (hyper)polarizability, it is necessary to include electron correlation. Based on the results, the convergence of the usual perturbation treatment for vibrational anharmonicity was examined

Formulació matricial de la teoria de pertorbacions de Rayleigh-Schrödringer: monodimensional

Es repassa la formulació de la Teoria de Pertorbacions en notació matricial i s'exposa una aplicació senzilla com és la solució del problema de la partícula sotmesa a un potencial d'atracció dins la caixa quàntica monodimensional

Estructura i dinàmica del sistema bentònic en llacunes costaneres dels aiguamolls de l'Empordà

S'ha estudiat l'estructura i dinàmica del sistema bentònic de llacunes costaneres de la maresma dels aiguamolls de l'Empordà (NE de la península Ibèrica) i els factors que les determinen. Amb aquesta finalitat es van prendre mostres d'organismes del bentos (mensualment), nutrients i pigments de l'aigua (setmanalment) i el sediment (mensualment), així com paràmetres físics de l'aigua (setmanalment), durant dos cicles d'inundació (1997-1998 i 1998-1999). En aquesta maresma els factors determinants per als organismes bentònics i el contingut en nutrients del sediment són principalment físics: el grau de permanència de l'aigua, el confinament i la granulometria. De manera que els factors tròfics tenen menys pes que els físics tal com es posa de manifest amb les associacions de nematodes. El règim hídric de les llacunes integra des del grau de permanència de l'aigua fins a la salinitat i el confinament, per la qual cosa a partir de la caracterització dels diferents tipus de règim hídric es poden discriminar les tres tipologies d'ambients amb comunitats bentòniques diferenciades en l'estructura, la composició i la dinàmica temporal. Així, les tres comunitats diferenciades són: (1) la d'aigües permanents caracteritzada per una riquesa i diversitat elevades i constants en el temps, i per un nombre més elevat de tàxons sense adaptacions per passar situacions adverses; (2) la d'aigües semipermanents i temporànies d'inundació pulsativa, caracteritzada per una menor riquesa i diversitat, i per la dominància d'organismes amb estructures de resistència per passar la fase seca i aguantar condicions desfavorables, i (3) la d'aigües temporànies d'inundació continuada, caracteritzada per una major variabilitat de la riquesa i diversitat al llarg del cicle i per una singularitat taxonòmica més elevada, ja que s'observen cladòcers i una major representació taxonòmica dels insectes. Durant la realització d'aquest estudi els dos cicles d'inundació van suposar entrades d'aigua en el sistema diferents, per això es va diferenciar un cicle humit (1997-1998) i un de sec (1998-1999). En relació amb la major importància de les entrades d'aigua en el sistema, es va observar un contingut en nutrients de l'aigua significativament superior durant el cicle humit. Pel que fa al sediment, la seva major inèrcia explicaria la manca de diferències significatives en el contingut en nutrients del sediment del cicle humit i del sec. Ara bé, sí que s'observa diferència en la composició d'organismes. Així, la variabilitat interanual en l'estructura de la comunitat bentònica de la maresma es pot atribuir a una situació d'estrès provocada pel cicle sec. Aquesta situació va comportar un augment de la diversitat per disminució de la dominància a totes les llacunes estudiades, determinada, en gran manera, per la disminució de l'abundància original de poblacions dominants i característiques d'ambients permanents, Corophium orientale, o temporànies, Gammarus aequicauda. No existeix un patró estacional en les comunitats estudiades, sinó que els principals canvis en el plàncton i el bentos són deguts a pertorbacions hídriques. Aquestes pertorbacions provoquen que s'estableixi una situació en què dominen uns tàxons característics en el plàncton (Synchaeta i Eutintinnus) i, en canvi, que en el bentos es perdi l'estructura pròpia de les comunitats característiques dels diferents ambients (comunitats dominades per Corophium orinetale en aigües permanents i per Chironomus salinarius en temporànies). La severitat és diferent en funció del grau de permanència de l'aigua, i és més alta en ambients permanents i més baixa en els temporanis. La major severitat en ambients permanents ve determinada per la menor adaptació del bentos d'aigües permanents a les fluctuacions ambientals i a l'efecte de la predació sobre el plàncton. A més, tant en el bentos com en el plàncton es donen fenòmens d'histèresi, segons els quals l'efecte de la pertorbació no solament depèn de la causa pertorbadora sinó també de l'estat en què es troba la comunitat en el moment de la pertorbació.

Implementation and application of basis set superposition error-correction schemes to the theoretical modeling of weak intermolecular interactions

This thesis deals with the so-called Basis Set Superposition Error (BSSE) from both a methodological and a practical point of view. The purpose of the present thesis is twofold: (a) to contribute step ahead in the correct characterization of weakly bound complexes and, (b) to shed light the understanding of the actual implications of the basis set extension effects in the ab intio calculations and contribute to the BSSE debate. The existing BSSE-correction procedures are deeply analyzed, compared, validated and, if necessary, improved. A new interpretation of the counterpoise (CP) method is used in order to define counterpoise-corrected descriptions of the molecular complexes. This novel point of view allows for a study of the BSSE-effects not only in the interaction energy but also on the potential energy surface and, in general, in any property derived from the molecular energy and its derivatives A program has been developed for the calculation of CP-corrected geometry optimizations and vibrational frequencies, also using several counterpoise schemes for the case of molecular clusters. The method has also been implemented in Gaussian98 revA10 package. The Chemical Hamiltonian Approach (CHA) methodology has been also implemented at the RHF and UHF levels of theory for an arbitrary number interacting systems using an algorithm based on block-diagonal matrices. Along with the methodological development, the effects of the BSSE on the properties of molecular complexes have been discussed in detail. The CP and CHA methodologies are used for the determination of BSSE-corrected molecular complexes properties related to the Potential Energy Surfaces and molecular wavefunction, respectively. First, the behaviour of both BSSE-correction schemes are systematically compared at different levels of theory and basis sets for a number of hydrogen-bonded complexes. The Complete Basis Set (CBS) limit of both uncorrected and CP-corrected molecular properties like stabilization energies and intermolecular distances has also been determined, showing the capital importance of the BSSE correction. Several controversial topics of the BSSE correction are addressed as well. The application of the counterpoise method is applied to internal rotational barriers. The importance of the nuclear relaxation term is also pointed out. The viability of the CP method for dealing with charged complexes and the BSSE effects on the double-well PES blue-shifted hydrogen bonds is also studied in detail. In the case of the molecular clusters the effect of high-order BSSE effects introduced with the hierarchical counterpoise scheme is also determined. The effect of the BSSE on the electron density-related properties is also addressed. The first-order electron density obtained with the CHA/F and CHA/DFT methodologies was used to assess, both graphically and numerically, the redistribution of the charge density upon BSSE-correction. Several tools like the Atoms in Molecules topologycal analysis, density difference maps, Quantum Molecular Similarity, and Chemical Energy Component Analysis were used to deeply analyze, for the first time, the BSSE effects on the electron density of several hydrogen bonded complexes of increasing size. The indirect effect of the BSSE on intermolecular perturbation theory results is also pointed out It is shown that for a BSSE-free SAPT study of hydrogen fluoride clusters, the use of a counterpoise-corrected PES is essential in order to determine the proper molecular geometry to perform the SAPT analysis.

A educação inclusiva e a hiperactividade no jardim de infância

RESUMO: Propomos com este trabalho conhecer e aprofundar conhecimentos, através duma investigação-acção, relativamente à problemática em estudo - Perturbação de Hiperactividade com Défice de Atenção (PHDA). Também procurámos mostrar, através da intervenção, que o educador pode auxiliar no desenvolvimento e aprendizagem de uma criança com um diagnóstico de PHDA. Realizámos uma pesquisa bibliográfica sobre esta Perturbação, com uma breve perspectiva histórica, sobre o que é a PHDA, o diagnóstico, os subtipos e problemas associados e a intervenção em crianças com esta problemática. Pesquisámos, também, estratégias educativas de como actuar a nível pedagógico com esta criança e com o grupo onde se encontra incluída. Demos uma grande centralidade ao trabalho empírico, através da pesquisa documental, observações, entrevistas e sociometria. De seguida, elaborámos um projecto de investigação-acção, para o grupo onde a criança está incluída, utilizando estratégias de educação inclusiva, de modo a colmatar problemas comportamentais e de aprendizagem, derivados de défices de atenção e de concentração do grupo, através do grupo. Todo este processo deu resultados positivos, a nível dos comportamentos sociais e das aprendizagens académicas do grupo e do aluno desencadeador do Trabalho Projecto. Obrigou a uma acção, reflexão, acção que desencadeou uma procura e experimentação sistemática para dar resposta aos problemas surgidos, o que ajudou ao meu desenvolvimento pessoal e profissional e a uma vontade de investir mais e melhor no desencadeamento de práticas de educação cada vez mais inclusivas. ABSTRACT: With this paper, we want to learn, using the investigation-action, about the team we are working with, that is Perturbation of Hyperactivity and Deficit of Attention (PHDA). We also want to show that it is possible to help the children with the PHDA in their learning process using the right intervention. We searched for different opinions on this team, its history, how to diagnostic it and the different problems that this children can have. We also search for different strategies to teach the students with this problem as well as the others of the group. We gave special attention to the empirical work. For this reason we used the documental research, the observation, interviews and sociometry. After this we prepared the action research project for the group of children that work with the PHDA child. To do so, we used different inclusion strategies so that we could help him improving his behavior and learning process using the help of other students. On the end we had positive results in the behavior and also in the learning process with the group and the student who we based this work on. We used action/reflection/action that helped us to find out about the PHDA and the different strategies to try to solve this problem. All of this helped me as a teacher and as a person. It gave me the will to keep looking for better and different inclusive strategies.

Lesser Scaup Population Dynamics: What Can Be Learned from Available Data?

Populations of Lesser Scaup (Aythya affinis) have declined markedly in North America since the early 1980s. When considering alternatives for achieving population recovery, it would be useful to understand how the rate of population growth is functionally related to the underlying vital rates and which vital rates affect population growth rate the most if changed (which need not be those that influenced historical population declines). To establish a more quantitative basis for learning about life history and population dynamics of Lesser Scaup, we summarized published and unpublished estimates of vital rates recorded between 1934 and 2005, and developed matrix life-cycle models with these data for females breeding in the boreal forest, prairie-parklands, and both regions combined. We then used perturbation analysis to evaluate the effect of changes in a variety of vital-rate statistics on finite population growth rate and abundance. Similar to Greater Scaup (Aythya marila), our modeled population growth rate for Lesser Scaup was most sensitive to unit and proportional change in adult female survival during the breeding and non-breeding seasons, but much less so to changes in fecundity parameters. Interestingly, population growth rate was also highly sensitive to unit and proportional changes in the mean of nesting success, duckling survival, and juvenile survival. Given the small samples of data for key aspects of the Lesser Scaup life cycle, we recommend additional research on vital rates that demonstrate a strong effect on population growth and size (e.g., adult survival probabilities). Our life-cycle models should be tested and regularly updated in the future to simultaneously guide science and management of Lesser Scaup populations in an adaptive context.

Evaluating the Small Population Paradigm for Rare Large-Bodied Woodpeckers, with Implications for the Ivory-billed Woodpecker

Six large-bodied, ≥ 120 g, woodpecker species are listed as near-threatened to critically endangered by the International Union for Conservation of Nature (IUCN). The small population paradigm assumes that these populations are likely to become extinct without an increase in numbers, but the combined influences of initial population size and demographic rates, i.e., annual adult survival and fecundity, may drive population persistence for these species. We applied a stochastic, stage-based single-population model to available demographic rates for Dryocopus and Campephilus woodpeckers. In particular, we determined the change in predicted extinction rate, i.e., proportion of simulated populations that went extinct within 100 yr, to concomitant changes in six input parameters. To our knowledge, this is the first study to evaluate the combined importance of initial population size and demographic rates for the persistence of large-bodied woodpeckers. Under a worse-case scenario, the median time to extinction was 7 yr (range: 1–32). Across the combinations of other input values, increasing initial population size by one female induced, on average, 0.4%–3.2% (range: 0%–28%) reduction in extinction rate. Increasing initial population size from 5–30 resulted in extinction rates < 0.05 under limited conditions: (1) all input values were intermediate, or (2) Allee effect present and annual adult survival ≥ 0.8. Based on our model, these species can persist as rare, as few as five females, and thus difficult-to-detect, populations provided they maintain ≥ 1.1 recruited females annually per adult female and an annual adult survival rate ≥ 0.8. Athough a demographic-based population viability analysis (PVA) is useful to predict how extinction rate changes across scenarios for life-history attributes, the next step for modeling these populations should incorporate more easily acquired data on changes in patch occupancy to make predictions about patch colonization and extinction rates.

Solar influences on dynamical coupling between the stratosphere and troposphere

We use a simplified atmospheric general circulation model (AGCM) to investigate the response of the lower atmosphere to thermal perturbations in the lower stratosphere. The results show that generic heating of the lower stratosphere tends to weaken the sub-tropical jets and the tropospheric mean meridional circulations. The positions of the jets, and the extent of the Hadley cells, respond to the distribution of the stratospheric heating, with low latitude heating displacing them poleward, and uniform heating displacing them equatorward. The patterns of response to the low latitude heating are similar to those found to be associated with solar variability in previous observational data analysis, and to the effects of varying solar UV radiation in sophisticated AGCMs. In order to investigate the chain of causality involved in converting the stratospheric thermal forcing to a tropospheric climate signal we conduct an experiment which uses an ensemble of model spin-ups to analyse the time development of the response to an applied stratospheric perturbation. We find that the initial effect of the change in static stability at the tropopause is to reduce the eddy momentum flux convergence in this region. This is followed by a vertical transfer of the momentum forcing anomaly by an anomalous mean circulation to the surface, where it is partly balanced by surface stress anomalies. The unbalanced part drives the evolution of the vertically integrated zonal flow. We conclude that solar heating of the stratosphere may produce changes in the circulation of the troposphere even without any direct forcing below the tropopause. We suggest that the impact of the stratospheric changes on wave propagation is key to the mechanisms involved.

Chemical composition observed over the mid-Atlantic and the detection of pollution signatures far from source regions

The atmospheric composition of the central North Atlantic region has been sampled using the FAAM BAe146 instrumented aircraft during the Intercontinental Transport of Ozone and Precursors (ITOP) campaign, part of the wider International Consortium for Atmospheric Research on Transport and Transformation (ICARTT). This paper presents an overview of the ITOP campaign. Between late July and early August 2004, twelve flights comprising 72 hours of measurement were made in a region from approximately 20 to 40°W and 33 to 47°N centered on Faial Island, Azores, ranging in altitude from 50 to 9000 m. The vertical profiles of O3 and CO are consistent with previous observations made in this region during 1997 and our knowledge of the seasonal cycles within the region. A cluster analysis technique is used to partition the data set into air mass types with distinct chemical signatures. Six clusters provide a suitable balance between cluster generality and specificity. The clusters are labeled as biomass burning, low level outflow, upper level outflow, moist lower troposphere, marine and upper troposphere. During this summer, boreal forest fire emissions from Alaska and northern Canada were found to provide a major perturbation of tropospheric composition in CO, PAN, organic compounds and aerosol. Anthropogenic influenced air from the continental boundary layer of the USA was clearly observed running above the marine boundary layer right across the mid-Atlantic, retaining high pollution levels in VOCs and sulfate aerosol. Upper level outflow events were found to have far lower sulfate aerosol, resulting from washout on ascent, but much higher PAN associated with the colder temperatures. Lagrangian links with flights of other aircraft over the USA and Europe show that such signatures are maintained many days downwind of emission regions. Some other features of the data set are highlighted, including the strong perturbations to many VOCs and OVOCs in this remote region.

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.

How well do we understand and evaluate climate change feedback processes?

Processes in the climate system that can either amplify or dampen the climate response to an external perturbation are referred to as climate feedbacks. Climate sensitivity estimates depend critically on radiative feedbacks associated with water vapor, lapse rate, clouds, snow, and sea ice, and global estimates of these feedbacks differ among general circulation models. By reviewing recent observational, numerical, and theoretical studies, this paper shows that there has been progress since the Third Assessment Report of the Intergovernmental Panel on Climate Change in (i) the understanding of the physical mechanisms involved in these feedbacks, (ii) the interpretation of intermodel differences in global estimates of these feedbacks, and (iii) the development of methodologies of evaluation of these feedbacks (or of some components) using observations. This suggests that continuing developments in climate feedback research will progressively help make it possible to constrain the GCMs’ range of climate feedbacks and climate sensitivity through an ensemble of diagnostics based on physical understanding and observations.