Monte Carlo numerical treatment of large linear algebra problems

In this paper we deal with performance analysis of Monte Carlo algorithm for large linear algebra problems. We consider applicability and efficiency of the Markov chain Monte Carlo for large problems, i.e., problems involving matrices with a number of non-zero elements ranging between one million and one billion. We are concentrating on analysis of the almost Optimal Monte Carlo (MAO) algorithm for evaluating bilinear forms of matrix powers since they form the so-called Krylov subspaces. Results are presented comparing the performance of the Robust and Non-robust Monte Carlo algorithms. The algorithms are tested on large dense matrices as well as on large unstructured sparse matrices.

Modeling the structure of amorphous MoS3: a neutron diffraction and reverse Monte Carlo study

A model for the structure of amorphous molybdenum trisulfide, a-MoS3, has been created using reverse Monte Carlo methods. This model, which consists of chains Of MoS6 units sharing three sulfurs with each of its two neighbors and forming alternate long, nonbonded, and short, bonded, Mo-Mo separations, is a good fit to the neutron diffraction data and is chemically and physically realistic. The paper identifies the limitations of previous models based on Mo-3 triangular clusters in accounting for the available experimental data.

The influence of the carbon surface chemical composition on Dubinin-Astakhov equation parameters calculated from SF(6) adsorption data-grand canonical Monte Carlo simulation

Using grand canonical Monte Carlo simulation we show, for the first time, the influence of the carbon porosity and surface oxidation on the parameters of the Dubinin-Astakhov (DA) adsorption isotherm equation. We conclude that upon carbon surface oxidation, the adsorption decreases for all carbons studied. Moreover, the parameters of the DA model depend on the number of surface oxygen groups. That is why in the case of carbons containing surface polar groups, SF(6) adsorption isotherm data cannot be used for characterization of the porosity.

Population genetic diversity of the endemic Sardinian newt Euproctus platycephalus: implications for conservation

The Sardinian mountain newt Euproctus platycephalus, endemic to the island of Sardinia, (Italy), is considered a rare and threatened species and is classed as critically endangered by IUCN. It inhabits streams, small lakes and pools on the main mountain systems of the island. Threats from climatic and anthropogenic factors have raised concerns for the long-term survival of newt populations on the island. MtDNA sequencing was used to investigate the genetic population structure and phylogeography of this endemic species. Patterns of genetic variation were assessed by sequencing the complete Dloop region and part of the 12SrRNA, from 74 individuals representing four different populations. Analyses of molecular variance suggest that populations are significantly differentiated, and the distribution of haplotypes across the island shows strong geographical structuring. However, phylogenetic analyses also suggest that the Sardinian population consists of two distinct mtDNA groups, which may reflect ancient isolation and expansion events. Population structure, evolutionary history of the species and implications for the conservation of newt populations are discussed.

A novel parallel approach to the likelihood-based estimation of admixture in population genetics

Inferring population admixture from genetic data and quantifying it is a difficult but crucial task in evolutionary and conservation biology. Unfortunately state-of-the-art probabilistic approaches are computationally demanding. Effectively exploiting the computational power of modern multiprocessor systems can thus have a positive impact to Monte Carlo-based simulation of admixture modeling. A novel parallel approach is briefly described and promising results on its message passing interface (MPI)-based C implementation are reported.

Modelling heterotachy in phylogenetic inference by reversible-jump Markov chain Monte Carlo

The rate at which a given site in a gene sequence alignment evolves over time may vary. This phenomenon-known as heterotachy-can bias or distort phylogenetic trees inferred from models of sequence evolution that assume rates of evolution are constant. Here, we describe a phylogenetic mixture model designed to accommodate heterotachy. The method sums the likelihood of the data at each site over more than one set of branch lengths on the same tree topology. A branch-length set that is best for one site may differ from the branch-length set that is best for some other site, thereby allowing different sites to have different rates of change throughout the tree. Because rate variation may not be present in all branches, we use a reversible-jump Markov chain Monte Carlo algorithm to identify those branches in which reliable amounts of heterotachy occur. We implement the method in combination with our 'pattern-heterogeneity' mixture model, applying it to simulated data and five published datasets. We find that complex evolutionary signals of heterotachy are routinely present over and above variation in the rate or pattern of evolution across sites, that the reversible-jump method requires far fewer parameters than conventional mixture models to describe it, and serves to identify the regions of the tree in which heterotachy is most pronounced. The reversible-jump procedure also removes the need for a posteriori tests of 'significance' such as the Akaike or Bayesian information criterion tests, or Bayes factors. Heterotachy has important consequences for the correct reconstruction of phylogenies as well as for tests of hypotheses that rely on accurate branch-length information. These include molecular clocks, analyses of tempo and mode of evolution, comparative studies and ancestral state reconstruction. The model is available from the authors' website, and can be used for the analysis of both nucleotide and morphological data.

Habitat correlates of distribution and local population decline of the endemic Sardinian newt Euproctus platycephalus

Habitat-based statistical models relating patterns of presence and absence of species to habitat variables could be useful to resolve conservation-related problems and highlight the causes of population declines. In this paper, we apply such a modelling approach to an endemic amphibian, the Sardinian mountain newt Euproctus platycephalus, considered by IUCN a critically endangered species. Sardinian newts inhabit freshwater habitat in streams, small lakes and pools on the island of Sardinia (Italy). Reported declines of newt populations are not yet supported by quantitative data, however, they are perceived or suspected across the species' historical range. This study represents a first attempt trying to statistically relate habitat characteristics to Sardinian newt occurrence and persistence. Linear regression analysis revealed that newts are more likely to be found in sites with colder water temperature, less riparian vegetation and, marginally, absence of fish. The implications of the results for the conservation of the species are discussed, and suggestions for the short-term management of newt inhabited sites suggested. (C) 2003 Elsevier Ltd. All rights reserved.

Geographical distribution of the endemic Sardinian brook salamander, Euproctus platycephalus, and implications for its conservation

The Sardinian brook salamander, Euproctus platycephalus, is a cryptically coloured urodele found in streams, springs and pools in the main mountain systems of Sardinia, and is classified as critically endangered by IUCN. General reviews of the mountainous range where salamanders occur are numerous, but very few field-based distribution studies exist on this endemic species. Through a field and questionnaire survey, conducted between 1999 and 2001, we report a first attempt to increase data on the present distribution of E. platycephalus. A total of 14 localities where Sardinian salamanders are represented by apparently stable and in some cases abundant populations have been identified, as well as 30 sites where species presence has been recorded after 1991. Some 11 historical sites were identified which are no longer inhabited by the species. The implications of this distributional study for the conservation of the species and for the realization of an updated atlas are discussed.

Testing isotherm models and recovering empirical relationships for adsorption in microporous carbons using virtual carbon models and grand canonical Monte Carlo simulations

Using the plausible model of activated carbon proposed by Harris and co-workers and grand canonical Monte Carlo simulations, we study the applicability of standard methods for describing adsorption data on microporous carbons widely used in adsorption science. Two carbon structures are studied, one with a small distribution of micropores in the range up to 1 nm, and the other with micropores covering a wide range of porosity. For both structures, adsorption isotherms of noble gases (from Ne to Xe), carbon tetrachloride and benzene are simulated. The data obtained are considered in terms of Dubinin-Radushkevich plots. Moreover, for benzene and carbon tetrachloride the temperature invariance of the characteristic curve is also studied. We show that using simulated data some empirical relationships obtained from experiment can be successfully recovered. Next we test the applicability of Dubinin's related models including the Dubinin-Izotova, Dubinin-Radushkevich-Stoeckli, and Jaroniec-Choma equations. The results obtained demonstrate the limits and applications of the models studied in the field of carbon porosity characterization.

The Spanish Ambassador and Samuel Daniel's Vision of the Twelve Goddesses: A New Document [with text]

The performance of Samuel Daniel's masque The Vision of the Twelve Goddesses at court on January 8, 1604 took place in the midst of the preliminary negotiations that would lead to the signing of the Anglo-Spanish peace at Somerset House the following August. Philip III sent a special ambassador to England to congratulate James on his accession, and a series of tussles between Juan de Tassis and his French counterpart ensued. As a recently-discovered document in the Archivo General de Simancas reveals, Anna of Denmark intervened personally to insure that de Tassis, and not the Frenchman, attended the masque. This was a clear signal of James and Anna's peace aims, which de Tassis conveyed to the King of Spain; moreover, he enclosed in his dispatch a text of Daniel's masque which he clearly considered both political intelligence and of interest to the theater-loving Hapsburg monarch. The Simancas text of the Daniel masque is a new version, hitherto unknown, which adds to our knowledge of the circumstances in which the first Stuart masque was performed. Here we present a transcription and annotated translation of both de Tassis' letter and the text of the masque he had compiled for Philip III. (B. C.-E. and M. H.)

Euthydemoss of Bactria

Euthydemos I (ca. 260–200 bce) was king of Bactria from around 230. He founded a dynasty which, most notably under his son Demetrios I, extended the control of the Greco-Bactrian kings south of the Hindu Kush into Arachosia and India.

The credit relationship between Henry III and merchants of Douai and Ypres, 1247-70

This article looks at an important but neglected aspect of medieval sovereign debt, namely ‘accounts payable’ owed by the Crown to merchants and employees. It focuses on the unusually well-documented relationship between Henry III, King of England between 1216 and 1272, and Flemish merchants from the towns of Douai and Ypres, who provided cloth on credit to the royal wardrobe. From the surviving royal documents, we reconstruct the credit advanced to the royal wardrobe by the merchants of Ypres and Douai for each year between 1247 and 1270, together with the king's repayment history. The interactions between the king and the merchants are then analysed. The insights from this analysis are applied to the historical data to explain the trading decisions made by the merchants during this period, as well as why the strategies of the Yprois sometimes differed from those of the Douaissiens.

Quantifying uncertainty in critical loads: (B) acidity mass balance critical loads on a sensitive site

This paper reports an uncertainty analysis of critical loads for acid deposition for a site in southern England, using the Steady State Mass Balance Model. The uncertainty bounds, distribution type and correlation structure for each of the 18 input parameters was considered explicitly, and overall uncertainty estimated by Monte Carlo methods. Estimates of deposition uncertainty were made from measured data and an atmospheric dispersion model, and hence the uncertainty in exceedance could also be calculated. The uncertainties of the calculated critical loads were generally much lower than those of the input parameters due to a "compensation of errors" mechanism - coefficients of variation ranged from 13% for CLmaxN to 37% for CL(A). With 1990 deposition, the probability that the critical load was exceeded was > 0.99; to reduce this probability to 0.50, a 63% reduction in deposition is required; to 0.05, an 82% reduction. With 1997 deposition, which was lower than that in 1990, exceedance probabilities declined and uncertainties in exceedance narrowed as deposition uncertainty had less effect. The parameters contributing most to the uncertainty in critical loads were weathering rates, base cation uptake rates, and choice of critical chemical value, indicating possible research priorities. However, the different critical load parameters were to some extent sensitive to different input parameters. The application of such probabilistic results to environmental regulation is discussed.

Landscape typology in the Mediterranean context: a tool for habitat restoration

Despite the wide use of Landscape Character Assessment (LCA) as a tool for landscape planning in NW Europe, there are few examples of its application in the Mediterranean. This paper reports on the results from the development of a typology for LCA in a study area of northern Sardinia, Italy to provide a spatial framework for the analysis of current patterns of cork oak distribution and future restoration of this habitat. Landscape units were derived from a visual interpretation of map data stored within a GIS describing the physical and cultural characteristics of the study area. The units were subsequently grouped into Landscape Types according to the similarity of shared attributes using Two Way Indicator Species Analysis (TWINSPAN). The preliminary results showed that the methodology classified distinct Landscape Types but, based on field observations, there is a need for further refinement of the classification. The distribution and properties of two main cork oak habitats types was examined within the identified Landscape Types namely woodlands and wood pastures using Patch Analyst. The results show very clearly a correspondence between the distribution of cork oak pastures and cork oak woodland and landscape types. This forms the basis of the development of strategies for the maintenance, restoration and recreation of these habitat types within the study area, ultimately for the whole island of Sardinia. Future work is required to improve the landscape characterisation , particularly with respect to cultural factors, and to determine the validity of the landscape spatial framework for the analysis of cork oak distribution as part of a programme of habitat restoration and re-creation.