Desenvolupament computacional de la semblança molecular quàntica

La present Tesi Doctoral, titulada desenvolupament computacional de la semblança molecular quàntica, tracta, fonamentalment, els aspectes de càlcul de mesures de semblança basades en la comparació de funcions de densitat electrònica.El primer capítol, Semblança quàntica, és introductori. S'hi descriuen les funcions de densitat de probabilitat electrònica i llur significança en el marc de la mecànica quàntica. Se n'expliciten els aspectes essencials i les condicions matemàtiques a satisfer, cara a una millor comprensió dels models de densitat electrònica que es proposen. Hom presenta les densitats electròniques, mencionant els teoremes de Hohenberg i Kohn i esquematitzant la teoria de Bader, com magnituds fonamentals en la descripció de les molècules i en la comprensió de llurs propietats.En el capítol Models de densitats electròniques moleculars es presenten procediments computacionals originals per l'ajust de funcions densitat a models expandits en termes de gaussianes 1s centrades en els nuclis. Les restriccions físico-matemàtiques associades a les distribucions de probabilitat s'introdueixen de manera rigorosa, en el procediment anomenat Atomic Shell Approximation (ASA). Aquest procediment, implementat en el programa ASAC, parteix d'un espai funcional quasi complert, d'on se seleccionen variacionalment les funcions o capes de l'expansió, d'acord als requisits de no negativitat. La qualitat d'aquestes densitats i de les mesures de semblança derivades es verifica abastament. Aquest model ASA s'estén a representacions dinàmiques, físicament més acurades, en quant que afectades per les vibracions nuclears, cara a una exploració de l'efecte de l'esmorteïment dels pics nuclears en les mesures de semblança molecular. La comparació de les densitats dinàmiques respecte les estàtiques evidencia un reordenament en les densitats dinàmiques, d'acord al que constituiria una manifestació del Principi quàntic de Le Chatelier. El procediment ASA, explícitament consistent amb les condicions de N-representabilitat, s'aplica també a la determinació directe de densitats electròniques hidrogenoides, en un context de teoria del funcional de la densitat.El capítol Maximització global de la funció de semblança presenta algorismes originals per la determinació de la màxima sobreposició de les densitats electròniques moleculars. Les mesures de semblança molecular quàntica s'identifiquen amb el màxim solapament, de manera es mesuri la distància entre les molècules, independentment dels sistemes de referència on es defineixen les densitats electròniques. Partint de la solució global en el límit de densitats infinitament compactades en els nuclis, es proposen tres nivells de aproximació per l'exploració sistemàtica, no estocàstica, de la funció de semblança, possibilitant la identificació eficient del màxim global, així com també dels diferents màxims locals. Es proposa també una parametrització original de les integrals de recobriment a través d'ajustos a funcions lorentzianes, en quant que tècnica d'acceleració computacional. En la pràctica de les relacions estructura-activitat, aquests avenços possibiliten la implementació eficient de mesures de semblança quantitatives, i, paral·lelament, proporcionen una metodologia totalment automàtica d'alineació molecular. El capítol Semblances d'àtoms en molècules descriu un algorisme de comparació dels àtoms de Bader, o regions tridimensionals delimitades per superfícies de flux zero de la funció de densitat electrònica. El caràcter quantitatiu d'aquestes semblances possibilita la mesura rigorosa de la noció química de transferibilitat d'àtoms i grups funcionals. Les superfícies de flux zero i els algorismes d'integració usats han estat publicats recentment i constitueixen l'aproximació més acurada pel càlcul de les propietats atòmiques. Finalment, en el capítol Semblances en estructures cristal·lines hom proposa una definició original de semblança, específica per la comparació dels conceptes de suavitat o softness en la distribució de fonons associats a l'estructura cristal·lina. Aquests conceptes apareixen en estudis de superconductivitat a causa de la influència de les interaccions electró-fonó en les temperatures de transició a l'estat superconductor. En aplicar-se aquesta metodologia a l'anàlisi de sals de BEDT-TTF, s'evidencien correlacions estructurals entre sals superconductores i no superconductores, en consonància amb les hipòtesis apuntades a la literatura sobre la rellevància de determinades interaccions.Conclouen aquesta tesi un apèndix que conté el programa ASAC, implementació de l'algorisme ASA, i un capítol final amb referències bibliogràfiques.

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.

Alfredo Pareja: ensayista y crítico (Del Historiador y el Ensayista)

La autora revisa la ensayística de Alfredo Pareja sobre temas de arte colonial, cultura y literatura. Da cuenta de la evolución de su noción de lo artístico desde 1933, cuando Pareja marca distancia con la generación anterior y con el arielismo vigente, al hacer suya la defensa del realismo social, más adelante, sin embargo, matizaría estos juicios en otros ensayos. Su postura crítica, en términos generales, no defiende la originalidad ni la forma del texto por sí mismas, sino el valor del contenido, el humanismo empieza a ser el hilo conductor de su narrativa e iluminador de su estética. Por otro lado, siente la necesidad de resignificar lo local en su contacto con lo universal, reconociendo el afán de contemporaneidad de su generación, su defensa de lo mestizo está presente en sus reflexiones sobre el barroquismo latinoamericano. Finalmente, resulta interesante que Pareja aproxime dos grandes momentos del arte ecuatoriano: la pintura y escultura colonial de la llamada escuela de Quito, y la propuesta estética de su generación.

California's Endangered Belding’s Savannah Sparrow (Passerculus sandwichensis beldingi): Tolerance of Pedestrian Disturbance

Salt marshes constitute habitat islands for many endemic animal species, particularly along the California coast, where urban sprawl has fragmented this habitat. Recreational activities in salt marshes have increased recently, posing an interesting problem: how do endemic species lacking alternative habitat modify their tolerance to humans? We assessed seasonal and site variations in three tolerance parameters (distances at which animals became alert, fled, and moved after fleeing) of California's endangered Belding’s Savannah Sparrow ((Passerculus sandwichensis beldingi). We approached individuals on trails in three salt marshes with different levels of vehicle and pedestrian traffic. Belding’s Savannah Sparrows became aware and fled at shorter distances in the salt marsh coincident with greater levels of recreational activity as a result of habituation or visual obstruction effects. Seasonal effects in tolerance varied between sites. Alert and flight initiation distances were higher in the pre-nesting than in the non-breeding season in the site with the highest levels of recreational use likely due to greater exposure of breeding individuals; however, the opposite seasonal trend was found in each of the two sites with relatively lower human use, probably because individuals were less spatially attached in the non-breeding season when they foraged in aggregations. Distance fled was greater in the non-breeding than in the breeding season. Our findings call for dynamic management of recreational activities in different salt marshes depending on the degree of exposure to humans and seasonal variations in tolerance. We recommend a minimum approaching distance of 63 m and buffer areas of 1.3 ha around Belding's Savannah Sparrows.

Declines of Aerial Insectivores in North America Follow a Geographic Gradient

North American birds that feed on aerial insects are experiencing widespread population declines. An analysis of the North American Breeding Bird Survey trend estimates for 1966 to 2006 suggests that declines in this guild are significantly stronger than in passerines in general. The pattern of decline also shows a striking geographical gradient, with aerial insectivore declines becoming more prevalent towards the northeast of North America. Declines are also more acute in species that migrate long distances compared to those that migrate short distances. The declines become manifest, almost without exception, in the mid 1980s. The taxonomic breadth of these downward trends suggests that declines in aerial insectivore populations are linked to changes in populations of flying insects, and these changes might be indicative of underlying ecosystem changes.

Do the Golden-winged Warbler and Blue-winged Warbler Exhibit Species-specific Differences in their Breeding Habitat Use?

We compared habitat features of Golden-winged Warbler (Vermivora chrysoptera) territories in the presence and absence of the Blue-winged Warbler (V. cyanoptera) on reclaimed coal mines in southeastern Kentucky, USA. Our objective was to determine whether there are species specific differences in habitat that can be manipulated to encourage population persistence of the Golden-winged Warbler. When compared with Blue-winged Warblers, Golden-winged Warblers established territories at higher elevations and with greater percentages of grass and canopy cover. Mean territory size (minimum convex polygon) was 1.3 ha (se = 0.1) for Golden-winged Warbler in absence of Blue-winged Warbler, 1.7 ha (se = 0.3) for Golden-winged Warbler coexisting with Blue-winged Warbler, and 2.1 ha (se = 0.3) for Blue-winged Warbler. Territory overlap occurred within and between species (18 of n = 73 territories, 24.7%). All Golden-winged and Blue-winged Warblers established territories that included an edge between reclaimed mine land and mature forest, as opposed to establishing territories in open grassland/shrubland habitat. The mean distance territories extended from a forest edge was 28.0 m (se = 3.8) for Golden-winged Warbler in absence of Blue-winged Warbler, 44.7 m (se = 5.7) for Golden-winged Warbler coexisting with Blue-winged Warbler, and 33.1 m (se = 6.1) for Blue-winged Warbler. Neither territory size nor distances to forest edges differed significantly between Golden-winged Warbler in presence or absence of Blue-winged Warbler. According to Monte Carlo analyses, orchardgrass (Dactylis glomerata), green ash (Fraxinus pennsylvanica) seedlings and saplings, and black locust (Robinia pseudoacacia) saplings were indicative of sites with only Golden-winged Warblers. Sericea lespedeza, goldenrod (Solidago spp.), clematis vine (Clematis spp.), and blackberry (Rubus spp.) were indicative of sites where both species occurred. Our findings complement recent genetic studies and add another factor for examining Golden-winged Warbler population decline. Further, information from our study will aid land managers in manipulating habitat for the Golden-winged Warbler.

Bachman’s Sparrow (Peucaea aestivalis) response to variation in the extent of burns conducted during the nesting season

Bachman’s Sparrow (Peucaea aestivalis), an endemic North American passerine, requires frequent (≤ 3 yr) prescribed fires to maintain preferred habitat conditions. Prescribed fires that coincide with the sparrow’s nesting season are increasingly used to manage sparrow habitat, but concerns exist regarding the effects that nesting-season fires may pose to this understory-dwelling species. Previous studies suggested that threats posed by fires might be lessened by reducing the extent of prescribed fires, thereby providing unburned areas close to the areas where fires eliminate ground-cover vegetation. To assess this hypothesis, we monitored color-marked male Bachman’s Sparrows on 2 sites where the extent of nesting-season fires differed 5-fold (> 70 ha vs. < 15 ha). Monthly survival for males did not differ between the large- and small-extent treatments, and survival rates exceeded 90% for all months except one during the second year of our study when fires were applied later in the season. Male densities also did not differ between treatments, but treatment-by-year interactions pointed to effects relating to the specific time that fires were applied. The distances separating observations of marked males before and after burns were smaller on small-extent treatments in the first year of study but larger on the small-extent treatments in the second year of study. Burn extents also had no consistent effect on postburn reproductive status. The largest extent we examined could have been too small to affect sparrow populations, but responses may also reflect sustainable metapopulation dynamics in a setting where a large sparrow population is maintained at a regional scale (> 100,000 ha) using frequent prescribed fire (≤ 2-yr return intervals). Additional research is needed regarding the effects that nesting-season fires may have on small, isolated populations as well as sites where much larger burn extents (> 100 ha) or longer burn intervals (> 2 yr) are used.

Biased representation of disturbance rates in the roadside sampling frame in boreal forests: implications for monitoring design

The North American Breeding Bird Survey (BBS) is the principal source of data to inform researchers about the status of and trend for boreal forest birds. Unfortunately, little BBS coverage is available in the boreal forest, where increasing concern over the status of species breeding there has increased interest in northward expansion of the BBS. However, high disturbance rates in the boreal forest may complicate roadside monitoring. If the roadside sampling frame does not capture variation in disturbance rates because of either road placement or the use of roads for resource extraction, biased trend estimates might result. In this study, we examined roadside bias in the proportional representation of habitat disturbance via spatial data on forest “loss,” forest fires, and anthropogenic disturbance. In each of 455 BBS routes, the area disturbed within multiple buffers away from the road was calculated and compared against the area disturbed in degree blocks and BBS strata. We found a nonlinear relationship between bias and distance from the road, suggesting forest loss and forest fires were underrepresented below 75 and 100 m, respectively. In contrast, anthropogenic disturbance was overrepresented at distances below 500 m and underrepresented thereafter. After accounting for distance from road, BBS routes were reasonably representative of the degree blocks they were within, with only a few strata showing biased representation. In general, anthropogenic disturbance is overrepresented in southern strata, and forest fires are underrepresented in almost all strata. Similar biases exist when comparing the entire road network and the subset sampled by BBS routes against the amount of disturbance within BBS strata; however, the magnitude of biases differed. Based on our results, we recommend that spatial stratification and rotating panel designs be used to spread limited BBS and off-road sampling effort in an unbiased fashion and that new BBS routes be established where sufficient road coverage exists.

A Newton method for pose refinement of 3D models

Different optimization methods can be employed to optimize a numerical estimate for the match between an instantiated object model and an image. In order to take advantage of gradient-based optimization methods, perspective inversion must be used in this context. We show that convergence can be very fast by extrapolating to maximum goodness-of-fit with Newton's method. This approach is related to methods which either maximize a similar goodness-of-fit measure without use of gradient information, or else minimize distances between projected model lines and image features. Newton's method combines the accuracy of the former approach with the speed of convergence of the latter.

Electrospinning atactic polystyrene: A neutron scattering study

Electrospinning is a method used to produce nanoscale to microscale sized polymer fibres. In this study we electrospin 1:1 blends of deuterated and hydrogenated atactic-Polystyrene from N,N-Dimethylformamide for small angle neutron scattering experiments in order to analyse the chain conformation in the electrospun fibres. Small angle neutron scattering was carried out on randomly orientated fibre mats obtained using applied voltages of 10kV-15kV and needle tip to collector distances of 20cm and 30cm. Fibre diameters varied from 3mm - 20mm. Neutron scattering data from fibre samples were compared with bulk samples of the same polymer blend. The scattering data indicates that there are pores and nanovoiding present in the fibres; this was confirmed by scanning electron microscopy. A model that combines the scattering from the pores and the labelled polymer chains was used to extract values for the radius of gyration. The radius of gyration in the fibres is found to vary little with the applied voltage, but varies with the initial solution concentration and fibre diameter. The values for the radius of gyration in the fibres are broadly equivalent to that of the bulk state.

Microscopic Mechanisms of Strain Hardening in Glassy Polymers

The mechanisms underlying the increase in stress for large mechanical strains of a polymer glass, quantified by the strain-hardening modulus, are still poorly understood. In the present paper we aim to elucidate this matter and present new mechanisms. Molecular-dynamics simulations of two polymers with very different strain-hardening moduli (polycarbonate and polystyrene) have been carried out. Nonaffine displacements occur because of steric hindrances and connectivity constraints. We argue that it is not necessary to introduce the concept of entanglements to understand strain hardening, but that hardening is rather coupled with the increase in the rate of nonaffine particle displacements. This rate increases faster for polycarbonate, which has the higher strain-hardening modulus. Also more nonaffine chain stretching is present for polycarbonate. It is shown that the inner distances of such a nonaffinely deformed chain can be well described by the inner distances of the worm-like chain, but with an effective stiffness length (equal to the Kuhn length for an infinite worm-like chain) that increases during deformation. It originates from the finite extensibility of the chain. In this way the increase in nonaffine particle displacement can be understood as resulting from an increase in the effective stiffness length of the perturbed chain during deformation, so that at larger strains a higher rate of plastic events in terms of nonaffine displacement is necessary, causing in turn the observed strain hardening in polymer glasses.

A field study of factors influencing the concentrations of a traffic-related pollutant in the vicinity of a complex urban junction

The paper describes a field study focused on the dispersion of a traffic-related pollutant within an area close to a busy intersection between two street canyons in Central London. Simultaneous measurements of airflow, traffic flow and carbon monoxide concentrations ([CO]) are used to explore the causes of spatial variability in [CO] over a full range of background wind directions. Depending on the roof-top wind direction, evidence of both flow channelling and recirculation regimes were identified from data collected within the main canyon and the intersection. However, at the intersection, the merging of channelled flows from the canyons increased the flow complexity and turbulence intensity. These features, coupled with the close proximity of nearby queuing traffic in several directions, led to the highest overall time-average measured [CO] occurring at the intersection. Within the main street canyon, the data supported the presence of a helical flow regime for oblique roof-top flows, leading to increased [CO] on the canyon leeward side. Predominant wind directions led to some locations having significantly higher diurnal average [CO] due to being mostly on the canyon leeward side during the study period. For all locations, small changes in the background wind direction could cause large changes in the in-street mean wind angle and local turbulence intensity, implying that dispersion mechanisms would be highly sensitive to small changes in above roof flows. During peak traffic flow periods, concentrations within parallel side streets were approximately four times lower than within the main canyon and intersection which has implications for controlling personal exposure. Overall, the results illustrate that pollutant concentrations can be highly spatially variable over even short distances within complex urban geometries, and that synoptic wind patterns, traffic queue location and building topologies all play a role in determining where pollutant hot spots occur.

Imagining migration: Placing children's understanding of 'moving house' in Malawi and Lesotho

Research pertaining to children's geographies has mainly focused on children's physical experiences of space, with their 'imagined geographies' receiving far less attention. The few studies of children's imagined geographies that exist tend to focus on children's national identities and their understanding of distant places. However, children's lives are not necessarily static and they often move between places. Research has not so far considered children's images of these transitional spaces or how such images are constructed. Through an examination of over 800 thematic drawings and stories, regarding 'moving house, produced by children aged 10-17 years in urban and rural communities of Lesotho and Malawi, this paper explores southern African children's representations of migration. The research considers how ideas of migration are culturally-constructed based on notions of family, home and kinship, particularly in relation to the fluid family structure characteristic of most southern African societies. The results suggest that most children imagine migration as a household rather than an individual process.. rarely including micro -migrations between extended family households in their drawings. Further, children's images of migration are place-rooted in everyday life experiences. Their representations concentrate on the reasons for migration, both negative and positive, which are specifically related to their local social and environmental situations and whether house moves take place locally or over longer distances. The paper concludes by exploring the implications of these conceptualisations of moving house for children's contemporary migration experiences, particularly in light of changing family structures due to the effects of the HIV/AIDS pandernic. (c) 2005 Elsevier Ltd. All rights reserved

Estimating the spatial scales of regionalized variables by nested sampling, hierarchical analysis of variance and residual maximum likelihood

The variogram is essential for local estimation and mapping of any variable by kriging. The variogram itself must usually be estimated from sample data. The sampling density is a compromise between precision and cost, but it must be sufficiently dense to encompass the principal spatial sources of variance. A nested, multi-stage, sampling with separating distances increasing in geometric progression from stage to stage will do that. The data may then be analyzed by a hierarchical analysis of variance to estimate the components of variance for every stage, and hence lag. By accumulating the components starting from the shortest lag one obtains a rough variogram for modest effort. For balanced designs the analysis of variance is optimal; for unbalanced ones, however, these estimators are not necessarily the best, and the analysis by residual maximum likelihood (REML) will usually be preferable. The paper summarizes the underlying theory and illustrates its application with data from three surveys, one in which the design had four stages and was balanced and two implemented with unbalanced designs to economize when there were more stages. A Fortran program is available for the analysis of variance, and code for the REML analysis is listed in the paper. (c) 2005 Elsevier Ltd. All rights reserved.

Distribution and mobility of trace elements in soils and vegetation around the mining and smelting areas of Tharsis, Riotinto and Huelva, Iberian Pyrite Belt, SW Spain

Trace elements may present an environmental hazard in the vicinity of mining and smelting activities. However, the factors controlling their distribution and transfer within the soil and vegetation systems are not always well defined. Total concentrations of up to 15,195 mg center dot kg (-1) As, 6,690 mg center dot kg(-1) Cu, 24,820 mg center dot kg(-1) Pb and 9,810 mg center dot kg(-1) Zn in soils, and 62 mg center dot kg(-1) As, 1,765 mg center dot kg(-1) Cu, 280 mg center dot kg(-1) Pb and 3,460 mg center dot kg (-1) Zn in vegetation were measured. However, unusually for smelters and mines of a similar size, the elevated trace element concentrations in soils were found to be restricted to the immediate vicinity of the mines and smelters (maximum 2-3 km). Parent material, prevailing wind direction, and soil physical and chemical characteristics were found to correlate poorly with the restricted trace element distributions in soils. Hypotheses are given for this unusual distribution: (1) the contaminated soils were removed by erosion or (2) mines and smelters released large heavy particles that could not have been transported long distances. Analyses of the accumulation of trace elements in vegetation (median ratios: As 0.06, Cu 0.19, Pb 0.54 and Zn 1.07) and the percentage of total trace elements being DTPA extractable in soils (median percentages: As 0.06%, Cu 15%, Pb 7% and Zn 4%) indicated higher relative trace element mobility in soils with low total concentrations than in soils with elevated concentrations.