Monte Carlo radiation hydrodynamics:Methods, tests and application to Type Ia supernova ejecta

In astrophysical systems, radiation-matter interactions are important in transferring energy and momentum between the radiation field and the surrounding material. This coupling often makes it necessary to consider the role of radiation when modelling the dynamics of astrophysical fluids. During the last few years, there have been rapid developments in the use of Monte Carlo methods for numerical radiative transfer simulations. Here, we present an approach to radiation hydrodynamics that is based on coupling Monte Carlo radiative transfer techniques with finite-volume hydrodynamical methods in an operator-split manner. In particular, we adopt an indivisible packet formalism to discretize the radiation field into an ensemble of Monte Carlo packets and employ volume-based estimators to reconstruct the radiation field characteristics. In this paper the numerical tools of this method are presented and their accuracy is verified in a series of test calculations. Finally, as a practical example, we use our approach to study the influence of the radiation-matter coupling on the homologous expansion phase and the bolometric light curve of Type Ia supernova explosions. © 2012 The Authors Monthly Notices of the Royal Astronomical Society © 2012 RAS.

Empirical performance of RSSI-based Monte Carlo localisation for active RFID patient tracking systems

The range of potential applications for indoor and campus based personnel localisation has led researchers to create a wide spectrum of different algorithmic approaches and systems. However, the majority of the proposed systems overlook the unique radio environment presented by the human body leading to systematic errors and inaccuracies when deployed in this context. In this paper RSSI-based Monte Carlo Localisation was implemented using commercial 868 MHz off the shelf hardware and empirical data was gathered across a relatively large number of scenarios within a single indoor office environment. This data showed that the body shadowing effect caused by the human body introduced path skew into location estimates. It was also shown that, by using two body-worn nodes in concert, the effect of body shadowing can be mitigated by averaging the estimated position of the two nodes worn on either side of the body. © Springer Science+Business Media, LLC 2012.

Monte Carlo simulation of complex cohesive fracture in random heterogeneous quasi-brittle materials

A numerical method is developed to simulate complex two-dimensional crack propagation in quasi-brittle materials considering random heterogeneous fracture properties. Potential cracks are represented by pre-inserted cohesive elements with tension and shear softening constitutive laws modelled by spatially varying Weibull random fields. Monte Carlo simulations of a concrete specimen under uni-axial tension were carried out with extensive investigation of the effects of important numerical algorithms and material properties on numerical efficiency and stability, crack propagation processes and load-carrying capacities. It was found that the homogeneous model led to incorrect crack patterns and load–displacement curves with strong mesh-dependence, whereas the heterogeneous model predicted realistic, complicated fracture processes and load-carrying capacity of little mesh-dependence. Increasing the variance of the tensile strength random fields with increased heterogeneity led to reduction in the mean peak load and increase in the standard deviation. The developed method provides a simple but effective tool for assessment of structural reliability and calculation of characteristic material strength for structural design.

Strong electronic correlation in the hydrogen chain: A variational Monte Carlo study

In this paper, we report a fully ab initio variational Monte Carlo study of the linear and periodic chain of hydrogen atoms, a prototype system providing the simplest example of strong electronic correlation in low dimensions. In particular, we prove that numerical accuracy comparable to that of benchmark density-matrix renormalization-group calculations can be achieved by using a highly correlated Jastrow-antisymmetrized geminal power variational wave function. Furthermore, by using the so-called "modern theory of polarization" and by studying the spin-spin and dimer-dimer correlations functions, we have characterized in detail the crossover between the weakly and strongly correlated regimes of this atomic chain. Our results show that variational Monte Carlo provides an accurate and flexible alternative to highly correlated methods of quantum chemistry which, at variance with these methods, can be also applied to a strongly correlated solid in low dimensions close to a crossover or a phase transition.

Body shadowing mitigation using differentiated LOS/NLOS channel models for RSSI-based Monte Carlo personnel localization

Research into localization has produced a wealth of algorithms and techniques to estimate the location of wireless network nodes, however the majority of these schemes do not explicitly account for non-line of sight conditions. Disregarding this common situation reduces their accuracy and their potential for exploitation in real world applications. This is a particular problem for personnel tracking where the user's body itself will inherently cause time-varying blocking according to their movements. Using empirical data, this paper demonstrates that, by accounting for non-line of sight conditions and using received signal strength based Monte Carlo localization, meter scale accuracy can be achieved for a wrist-worn personnel tracking tag in a 120 m indoor office environment. © 2012 IEEE.

Quantum annealing of an ising spin-glass by Green's function Monte Carlo

We present an implementation of quantum annealing (QA) via lattice Green's function Monte Carlo (GFMC), focusing on its application to the Ising spin glass in transverse field. In particular, we study whether or not such a method is more effective than the path-integral Monte Carlo- (PIMC) based QA, as well as classical simulated annealing (CA), previously tested on the same optimization problem. We identify the issue of importance sampling, i.e., the necessity of possessing reasonably good (variational) trial wave functions, as the key point of the algorithm. We performed GFMC-QA runs using such a Boltzmann-type trial wave function, finding results for the residual energies that are qualitatively similar to those of CA (but at a much larger computational cost), and definitely worse than PIMC-QA. We conclude that, at present, without a serious effort in constructing reliable importance sampling variational wave functions for a quantum glass, GFMC-QA is not a true competitor of PIMC-QA.

Monte Carlo studies of quantum and classical annealing on a double well

We present results for a variety of Monte Carlo annealing approaches, both classical and quantum, benchmarked against one another for the textbook optimization exercise of a simple one-dimensional double well. In classical (thermal) annealing, the dependence upon the move chosen in a Metropolis scheme is studied and correlated with the spectrum of the associated Markov transition matrix. In quantum annealing, the path integral Monte Carlo approach is found to yield nontrivial sampling difficulties associated with the tunneling between the two wells. The choice of fictitious quantum kinetic energy is also addressed. We find that a "relativistic" kinetic energy form, leading to a higher probability of long real-space jumps, can be considerably more effective than the standard nonrelativistic one.

Monte Carlo calculations of quantum yield in inhomogeneous PtSi/p-Si Schottky barriers

Monte Carlo calculations of quantum yield in PtSi/p-Si infrared detectors are carried out taking into account the presence of a spatially distributed barrier potential. In the 1-4 mu m wavelength range it is found that the spatial inhomogeneity of the barrier has no significant effect on the overall device photoresponse. However, above lambda = 4.0 mu m and particularly as the cut-off wavelength (lambda approximate to 5.5 mu m) is approached, these calculations reveal a difference between the homogeneous and inhomogeneous barrier photoresponse which becomes increasingly significant and exceeds 50% at lambda = 5.3 mu m. It is, in fact, the inhomogeneous barrier which displays an increased photoyield, a feature that is confirmed by approximate analytical calculations assuming a symmetric Gaussian spatial distribution of the barrier. Furthermore, the importance of the silicide layer thickness in optimizing device efficiency is underlined as a trade-off between maximizing light absorption in the silicide layer and optimizing the internal yield. The results presented here address important features which determine the photoyield of PtSi/Si Schottky diodes at energies below the Si absorption edge and just above the Schottky barrier height in particular.

Root Cause Analysis by a Combined Sparse Classification and Monte Carlo Approach

Classification methods with embedded feature selection capability are very appealing for the analysis of complex processes since they allow the analysis of root causes even when the number of input variables is high. In this work, we investigate the performance of three techniques for classification within a Monte Carlo strategy with the aim of root cause analysis. We consider the naive bayes classifier and the logistic regression model with two different implementations for controlling model complexity, namely, a LASSO-like implementation with a L1 norm regularization and a fully Bayesian implementation of the logistic model, the so called relevance vector machine. Several challenges can arise when estimating such models mainly linked to the characteristics of the data: a large number of input variables, high correlation among subsets of variables, the situation where the number of variables is higher than the number of available data points and the case of unbalanced datasets. Using an ecological and a semiconductor manufacturing dataset, we show advantages and drawbacks of each method, highlighting the superior performance in term of classification accuracy for the relevance vector machine with respect to the other classifiers. Moreover, we show how the combination of the proposed techniques and the Monte Carlo approach can be used to get more robust insights into the problem under analysis when faced with challenging modelling conditions.

Self-consistent modelling of line-driven hot-star winds with Monte Carlo radiation hydrodynamics

Radiative pressure exerted by line interactions is a prominent driver of outflows in astrophysical systems, being at work in the outflows emerging from hot stars or from the accretion discs of cataclysmic variables, massive young stars and active galactic nuclei. In this work, a new radiation hydrodynamical approach to model line-driven hot-star winds is presented. By coupling a Monte Carlo radiative transfer scheme with a finite volume fluid dynamical method, line-driven mass outflows may be modelled self-consistently, benefiting from the advantages of Monte Carlo techniques in treating multiline effects, such as multiple scatterings, and in dealing with arbitrary multidimensional configurations. In this work, we introduce our approach in detail by highlighting the key numerical techniques and verifying their operation in a number of simplified applications, specifically in a series of self-consistent, one-dimensional, Sobolev-type, hot-star wind calculations. The utility and accuracy of our approach are demonstrated by comparing the obtained results with the predictions of various formulations of the so-called CAK theory and by confronting the calculations with modern sophisticated techniques of predicting the wind structure. Using these calculations, we also point out some useful diagnostic capabilities our approach provides. Finally, we discuss some of the current limitations of our method, some possible extensions and potential future applications.

Es teatro vuestra oración: Santa Teresa’s performed confession in Juan Mayorga’s La lengua en pedazos

Juan Mayorga’s La Lengua en Pedazos (2010) strikes at the heart of the compositional circumstances of St Teresa's Libro de la Vida– staging, and arguably heightening the origins of her rhetorical strategies, the sense of awareness of readership and potential censure we encounter within the Libro de la Vida. His inquisitor refuses to be complicit in the tacit agreement that the word spoken in the theatrical space can conjure new realities –insistent on underscoring the textual origin of the visions painfully and partially offered up for his and our scrutiny. I will suggest that the persistent undertow towards a meta-commentary on the unmaking and remaking of the autobiographical text creates an unresolved tension between Teresa’s eloquent ability to take the spectator to a place beyond language, and our awareness that we are in the presence of a consummate performer, the textual source for the script itself produced with a supreme awareness of audience scrutiny. The play reflects ongoing lines of inquiry in our evolving understanding of the cultural production of Teresa and other holy women of the Early Modern period.

Técnica de canto lírico e de teatro musical: práticas de crossover

O presente trabalho tem como objectivo estudar as práticas de crossover entre as técnicas vocais que os cantores líricos e de teatro musical utilizam, assim como as técnicas vocais subjacentes às referidas práticas. Descrevem-se as técnicas vocais utilizadas no canto lírico e no teatro musical por pedagogos que fundamentam o seu trabalho com as descobertas da investigação na área da voz, e comparam-se as referidas técnicas para entender quais os pontos comuns e quais os pontos divergentes. Devido à elevada percentagem de pontos comuns às duas técnicas concluiu-se que são muito próximas entre si, o que faz sentido por serem executadas pela mesma estrutura fisiológica. Apresentam-se as entrevistas efectuadas a cantores profissionais de canto lírico e de teatro musical sobre os aspectos fundamentais das técnicas vocais que utilizam e como fazem o crossover entre as mesmas. Dos resultados dos inquéritos concluiu-se que a maioria dos cantores utiliza habitualmente práticas de crossover na sua performance. A segunda conclusão retirada dos resultados do inquérito foi que a execução das referidas práticas é intuitiva na maioria dos casos, e não conscientemente efectuada. Apresenta-se um caso de aplicação em contexto performativo das práticas de crossover: o papel do soprano na cantata cénica "Moby Dick - Aos Peixes". A utilização tecnicamente consciente das práticas de crossover permitiu estabilizar a execução vocal desde o início dos ensaios e obter posteriormente uma performance consistente mas versátil, sem fixar a execução vocal ao longo da carreira do espectáculo. Os apêndices incluem informação anatomofisiológica útil para este estudo, um resumo dos métodos de estudo científico da voz, o questionário utilizado no inquérito, as tabelas dos dados obtidos, a partitura anotada de "Aos Peixes" e o DVD do espectáculo realizado no Centro Cultural de Belém, em Lisboa. ABSTRACT: The present work aims to understand the crossover mecanisms that the classical singers and the musical theatre singers use, and the vocal techniques underlying those practices. This work describes the vocal techniques taught in the lyrical singing and in the musical theatre singing by teachers who base their pedagogy on the findings of scientific investigation of voice and singing. The techniques used in both fields are compared to understand their similarities and diferences. This process led to the conclusion that both techniques are very close, due to the high percentage of common items found, and this makes sense since both techniques are produced by the same physiologic struture: the vocal system. Professional singers from the lyrical and the musical theatre scene were interviewed to explain the basic foundations of their vocal technique and how do they do the crossover between those styles, from a technical point of view. The results of these interviews led to the conclusion that the majority of the singers performs crossover actions in their singing. The second conclusion is that for the majority of singers these crossover actions are intuitive, inspired by the music, the text or the dramatic context, and not informed by technically conscious actions. It is presented a case study of how these crossover methods were used in a staged cantata: the soprano role in "Moby Dick - Aos Peixes". The use of conscious technically informed crossover practice allowed to stabilize the vocal execution right from the beginning of rehearsals and obtain afterwards a performance which was both consistent and versatile, not fixed, during the running of the show. The appendixes include useful anatomical and physiological information, a summary of the methods of the scientific study of voice, the formulary used in the enquiry, the data tables of the field work, the annotated score of "Moby Dick - Aos Peixes" and the DVD from the play filmed at Centro Cultural de Belém, in Lisbon.

Para estudar o teatro do Renascimento italiano. A questão da tragédia (explorações)

Em Para estudar o teatro do Renascimento italiano. A questão da tragédia (explorações) estabelecem-se as pistas para o reconhecimento do fracasso italiano na elaboração quinhentista de textos trágicos, equacionando as causas desse fracasso e assinalando alguns dos textos e autores cuja produção marcou esse trajecto.

Estudo de Monte Carlo da electroluminescência em detectores gasosos

Uma descrição detalhada do processo de electroluminescência é um prérequisito na optimização de detectores gasosos para sistemas de imagiologia, astrofísica, física de altas energias e experiências de eventos raros. Neste trabalho, é apresentada e caracterizada uma nova e versátil plataforma de simulação da emissão de luz durante a deriva de electrões em gases nobres, desenvolvida usando os programas Magboltz e Garfield. Propriedades intrínsecas da electroluminescência em gases nobres são calculadas e apresentadas em função do campo eléctrico aplicado, nomeadamente eficiências, rendimento e flutuações estatísticas associadas. São obtidos resultados em grande concordância com dados experimentais e simulações Monte Carlo anteriores. A plataforma é usada para determinar as condições óptimas de funcionamento de detectores como o NEXT (Neutrino Experiment with a Xenon TPC) e outros baseados nas micro-estruturas GEM (Gas Electron Multiplier) e MHSP (Micro- Hole & Strip Plate).