Parallel Monte Carlo sampling scheme for sphere and hemisphere

The sampling of certain solid angle is a fundamental operation in realistic image synthesis, where the rendering equation describing the light propagation in closed domains is solved. Monte Carlo methods for solving the rendering equation use sampling of the solid angle subtended by unit hemisphere or unit sphere in order to perform the numerical integration of the rendering equation. In this work we consider the problem for generation of uniformly distributed random samples over hemisphere and sphere. Our aim is to construct and study the parallel sampling scheme for hemisphere and sphere. First we apply the symmetry property for partitioning of hemisphere and sphere. The domain of solid angle subtended by a hemisphere is divided into a number of equal sub-domains. Each sub-domain represents solid angle subtended by orthogonal spherical triangle with fixed vertices and computable parameters. Then we introduce two new algorithms for sampling of orthogonal spherical triangles. Both algorithms are based on a transformation of the unit square. Similarly to the Arvo's algorithm for sampling of arbitrary spherical triangle the suggested algorithms accommodate the stratified sampling. We derive the necessary transformations for the algorithms. The first sampling algorithm generates a sample by mapping of the unit square onto orthogonal spherical triangle. The second algorithm directly compute the unit radius vector of a sampling point inside to the orthogonal spherical triangle. The sampling of total hemisphere and sphere is performed in parallel for all sub-domains simultaneously by using the symmetry property of partitioning. The applicability of the corresponding parallel sampling scheme for Monte Carlo and Quasi-D/lonte Carlo solving of rendering equation is discussed.

Insights into redox cycling of sulfur and iron in peatlands using high-resolution diffusive equilibrium thin film (DET) gel probe sampling

High spatial resolution vertical profiles of pore-water chemistry have been obtained for a peatland using diffusive equilibrium in thin films (DET) gel probes. Comparison of DET pore-water data with more traditional depth-specific sampling shows good agreement and the DET profiling method is less invasive and less likely to induce mixing of pore-waters. Chloride mass balances as water tables fell in the early summer indicate that evaporative concentration dominates and there is negligible lateral flow in the peat. Lack of lateral flow allows element budgets for the same site at different times to be compared. The high spatial resolution of sampling also enables gradients to be observed that permit calculations of vertical fluxes. Sulfate concentrations fall at two sites with net rates of 1.5 and 5.0nmol cm− 3 day− 1, likely due to a dominance of bacterial sulfate reduction, while a third site showed a net gain in sulfate due to oxidation of sulfur over the study period at an average rate of 3.4nmol cm− 3 day− 1. Behaviour of iron is closely coupled to that of sulfur; there is net removal of iron at the two sites where sulfate reduction dominates and addition of iron where oxidation dominates. The profiles demonstrate that, in addition to strong vertical redox related chemical changes, there is significant spatial heterogeneity. Whilst overall there is evidence for net reduction of sulfate within the peatland pore-waters, this can be reversed, at least temporarily, during periods of drought when sulfide oxidation with resulting acid production predominates.

Excited states of nitro-polypyridine metal complexes and their ultrafast decay. Time-resolved IR absorption, spectroelectrochemistry, and TD-DFT calculations of fac-[Re(Cl)(CO)3(5-Nitro-1,10-phenanthroline)]

The lowest absorption band of fac-[Re(Cl)(CO)(3)(5-NO2-phen)] encompasses two close-lying MLCT transitions. The lower one is directed to LUMO, which is heavily localized on the NO2 group. The UV-vis absorption spectrum is well accounted for by TD-DFT (G03/PBEPBE1/CPCM), provided that the solvent, MeCN, is included in the calculations. Near-UV excitation of fac-[Re(Cl)(CO)(3)(5-NO2-phen)] populates a triplet metal to ligand charge-transfer excited state, (MLCT)-M-3, that was characterized by picosecond time-resolved IR spectroscopy. Large positive shifts of the v(CO) bands upon excitation (+70 cm(-1) for the A'(1) band) signify a very large charge separation between the Re(Cl)(CO)3 unit and the 5-NO2-phen ligand. Details of the excited-state character are revealed by TD-DFT calculated changes of electron density distribution. Experimental excited-state v(CO) wavenumbers agree well with those calculated by DFT. The (MLCT)-M-3 state decays with a ca. 10 ps lifetime (in MeCN) into another transient species, that was identified by TRIR and TD-DFT calculations as an intraligand (3)n pi* excited state, whereby the electron density is excited from the NO2 oxygen lone pairs to the pi* system of 5-NO2-phen. This state is short-lived, decaying to the ground state with a similar to 30 ps lifetime. The presence of an n pi* state seems to be the main factor responsible for the lack of emission and the very short lifetimes of 3 MLCT states seen in all d(6)-metal complexes of nitro-polypyridyl ligands. Localization of the excited electron density in the lowest (MLCT)-M-3 states parallels localization of the extra electron in the reduced state that is characterized by a very small negative shift of the v(CO) IR bands (-6 cm(-1) for A'(1)) but a large downward shift of the v(s)(NO2) IR band. The Re-Cl bond is unusually stable toward reduction, whereas the Cl ligand is readily substituted upon oxidation.

Rhenium-to-benzoylpyridine and rhenium-to-bipyridine MLCT excited states offac-[Re(Cl)(4-benzoylpyridine)2(CO)3] andfac-[Re(4-benzoylpyridine)(CO)3(bpy)]+: a time-resolved spectroscopic and spectroelectrochemical study

The lowest allowed electronic transition of fac-[Re(Cl)(CO)(3)(bopy)(2)] (bopy = 4-benzoylpyridine) has a Re --> bopy MLCT character, as revealed by UV-vis and stationary resonance Raman spectroscopy. Accordingly, the lowest-lying, long-lived, excited state is Re --> bopy (MLCT)-M-3. Electronic depopulation of the Re(CO)(3) unit and population of a bopy pi* orbital upon excitation are evident by the upward shift of v(Cequivalent toO) vibrations and a downward shift of the ketone v(C=O) vibration, respectively, seen in picosecond time-resolved IR spectra. Moreover, reduction of a single bopy ligand in the (MLCT)-M-3 excited state is indicated by time-resolved visible and resonance Raman (TR3) spectra that show features typical of bopy(.-). In contrast, the lowest allowed electronic transition and lowest-lying excited state of a new complex fac-[Re(bopy)(CO)(3)(bpy)](+) (bpy = 2,2'-bipyridine) have been identified as Re --> bpy MLCT with no involvement of the bopy ligand, despite the fact that the first reduction of this complex is bopy-localized, as was proven spectroelectrochemically. This is a rare case in which the localizations of the lowest MLCT excitation and the first reduction are different. (MLCT)-M-3 excited states of both fac-[Re(Cl)(CO)(3)(bopy)(2)] and fac-[Re(bopy)(CO)(3)(bpy)](+) are initially formed vibrationally hot. Their relaxation is manifested by picosecond dynamic shifts of v(Cequivalent toO) IR bands. The X-ray structure of fac-[Re(bopy)(CO)(3)(bpy)](PF6CH3CN)-C-. has been determined.

Perturbation, extraction and refinement of invariant pairs for matrix polynomials

Generalizing the notion of an eigenvector, invariant subspaces are frequently used in the context of linear eigenvalue problems, leading to conceptually elegant and numerically stable formulations in applications that require the computation of several eigenvalues and/or eigenvectors. Similar benefits can be expected for polynomial eigenvalue problems, for which the concept of an invariant subspace needs to be replaced by the concept of an invariant pair. Little has been known so far about numerical aspects of such invariant pairs. The aim of this paper is to fill this gap. The behavior of invariant pairs under perturbations of the matrix polynomial is studied and a first-order perturbation expansion is given. From a computational point of view, we investigate how to best extract invariant pairs from a linearization of the matrix polynomial. Moreover, we describe efficient refinement procedures directly based on the polynomial formulation. Numerical experiments with matrix polynomials from a number of applications demonstrate the effectiveness of our extraction and refinement procedures.

De-agrarianization, re-agrarianization and local economic development: re-orientating livelihoods in African artisanal mining communities

This article contributes to the debate on livelihood diversification in rural sub-Saharan Africa, focusing specifically on the growing economic importance of artisanal and small-scale mining (ASM) in the region. The precipitous decline in the value of many export crops and the removal of subsidies on crucial inputs such as fertilizers have made smallholder production unviable, forcing many farmers to ‘branch out’ into non-farm activities to supplement their incomes. One of the more popular destinations for poor farmers is the low-tech ASM sector which, because of its low barriers to entry, has absorbed millions of rural Africans over the past two decades, the majority of whom are engaged in the extraction of near-surface mineral deposits located on concessions that have been demarcated to multinational corporations. The efforts made hitherto to control this illegal mining activity, both through force and regulation, however, have had little effect, forcing many of the region’s governments and private sector partners to ‘re-think’ their approaches. One strategy that has gained considerable attention throughout the region is intensified support for agrarian-orientated activities, many of which, despite the problems plaguing smallholder agricultural sector and challenges with making it more economically sustainable, are being lauded as appropriate ‘alternative’ sources of employment to artisanal mining. After examining where artisanal mining fits into the de-agrarianization ‘puzzle’ in sub-Saharan Africa, the article critiques the efficacy of ‘re-agrarianization’ as a strategy for addressing the region’s illegal mining problem. A case study of Ghana is used to shed further light on these issues.

Non-intrusive monitoring of atmospheric CO2 in analogue models of terrestrial carbon cycle

1. Closed Ecological Systems (CES) are small manmade ecosystems which do not have any material exchange with the surrounding environment. Recent ecological and technological advances enable successful establishment and maintenance of CES, making them a suitable tool for detecting and measuring subtle feedbacks and mechanisms. 2. As a part of an analogue (physical) C cycle modelling experiment, we developed a non-intrusive methodology to control the internal environment and to monitor atmospheric CO2 concentration inside 16 replicated CES. Whilst maintaining an air-tight seal of all CES, this approach allowed for access to the CO2 measuring equipment for periodic re-calibration and repairs. 3. To ensure reliable cross-comparison of CO2 observations between individual CES units and to minimise the cost of the system, only one CO2 sampling unit was used. An ADC BioScientific OP-2 (open-path) analyser mounted on a swinging arm was passing over a set of 16 measuring cells. Each cell was connected to an individual CES with air continuously circulating between them. 4. Using this setup, we were able to continuously measure several environmental variables and CO2 concentration within each closed system, allowing us to study minute effects of changing temperature on C fluxes within each CES. The CES and the measuring cells showed minimal air leakage during an experimental run lasting, on average, 3 months. The CO2 analyser assembly performed reliably for over 2 years, however an early iteration of the present design proved to be sensitive to positioning errors. 5. We indicate how the methodology can be further improved and suggest possible avenues where future CES based research could be applied.

Assessment of the RE(OH)3 Ising magnetic materials as possible candidates for the study of transverse-field-induced quantum phase transitions

The LiHoxY1−xF4 Ising magnetic material subject to a magnetic field perpendicular to the Ho3+ Ising direction has shown over the past 20 years to be a host of very interesting thermodynamic and magnetic phenomena. Unfortunately, the availability of other magnetic materials other than LiHoxY1−xF4 that may be described by a transverse-field Ising model remains very much limited. It is in this context that we use here a mean-field theory to investigate the suitability of the Ho(OH)3, Dy(OH)3, and Tb(OH)3 insulating hexagonal dipolar Ising-type ferromagnets for the study of the quantum phase transition induced by a magnetic field, Bx, applied perpendicular to the Ising spin direction. Experimentally, the zero-field critical (Curie) temperatures are known to be Tc≈2.54, 3.48, and 3.72 K, for Ho(OH)3, Dy(OH)3, and Tb(OH)3, respectively. From our calculations we estimate the critical transverse field, Bxc, to destroy ferromagnetic order at zero temperature to be Bxc=4.35, 5.03, and 54.81 T for Ho(OH)3, Dy(OH)3, and Tb(OH)3, respectively. We find that Ho(OH)3, similarly to LiHoF4, can be quantitatively described by an effective S=1/2 transverse-field Ising model. This is not the case for Dy(OH)3 due to the strong admixing between the ground doublet and first excited doublet induced by the dipolar interactions. Furthermore, we find that the paramagnetic (PM) to ferromagnetic (FM) transition in Dy(OH)3 becomes first order for strong Bx and low temperatures. Hence, the PM to FM zero-temperature transition in Dy(OH)3 may be first order and not quantum critical. We investigate the effect of competing antiferromagnetic nearest-neighbor exchange and applied magnetic field, Bz, along the Ising spin direction ẑ on the first-order transition in Dy(OH)3. We conclude from these preliminary calculations that Ho(OH)3 and Dy(OH)3 and their Y3+ diamagnetically diluted variants, HoxY1−x(OH)3 and DyxY1−x(OH)3, are potentially interesting systems to study transverse-field-induced quantum fluctuations effects in hard axis (Ising-type) magnetic materials.

Epitaxial growth of ultrathin palladium films on Re{0001}

Ultrathin bimetallic layers create unusual magnetic and surface chemical effects through the modification of electronic structure brought on by low dimensionality, polymorphism, reduced screening, and epitaxial strain. Previous studies have related valence and core-level shifts to surface reactivity through the d-band model of Hammer and Nørskov, and in heteroepitaxial films this band position is determined by competing effects of coordination, strain, and hybridization of substrate and overlayer states. In this study we employ the epitaxially matched Pd on Re{0001} system to grow films with no lateral strain. We use a recent advancement in low-energy electron diffraction to expand the data range sufficiently for a reliable determination of the growth sequence and out-of-plane surface relaxation as a function of film thickness. The results are supported by scanning tunneling microscopy and X-ray photoemission spectroscopy, which show that the growth is layer-by-layer with significant core-level shifts due to changes in film structure, morphology, and bonding.

Análise de dados coletados por respondent-driven sampling(RDS): um estudo da prevalência de HIV e fatores associados entre mulheres trabalhadoras do sexo em 10 cidades brasileiras

O grupo das mulheres trabalhadoras do sexo (MTS) é reconhecido como uma populaçãode maior risco à infecção pelo HIV, tanto pela prevalência elevada, como por suavulnerabilidade social como pelos fatores relacionados à própria atividade profissional. Porém, arealização de estudos nos subgrupos de maior risco ao HIV mediante estratégias convencionaisde amostragem é, em geral, problemática por essas populações possuírem pequena magnitudeem termos populacionais e por estarem vinculados a comportamentos estigmatizados ouatividades ilegais. Em 1997, foi proposto um método de amostragem probabilística parapopulações de difícil acesso denominado Respondent-Driven Sampling (RDS). O método éconsiderado como uma variante da amostragem em cadeia e possibilita a estimação estatísticados parâmetros de interesse. Na literatura internacional, para análise de dados coletados porRDS, muitos autores têm utilizado técnicas estatísticas multivariadas tradicionais, sem levar emconta a estrutura de dependência das observações, presente nos dados coletados por RDS.A presente tese tem por objetivo contribuir para suprir informações sobre as práticas derisco relacionadas ao HIV entre as mulheres trabalhadoras do sexo (MTS) com odesenvolvimento de método estatístico para análise de dados coletados com o método deamostragem RDS. Com tal finalidade, foram utilizadas as informações coletadas na PesquisaCorrente da Saúde realizada em dez cidades brasileiras, com 2.523 MTS recrutadas por RDS,entre os anos de 2008 e 2009. O questionário foi autopreenchido e incluiu módulos sobrecaracterísticas da atividade profissional, práticas sexuais, uso de drogas, testes periódicos deHIV, e acesso aos serviços de saúde.Primeiramente, foram descritos alguns pressupostos do RDS e todas as etapas deimplantação da pesquisa. Em seguida, foram propostos métodos de análise multivariada, considerando o RDS como um desenho complexo de amostragem.