Dans cette thèse on étudie des schémas numériques pour des processus<br />/X/ à coefficients discontinus. Un premier schéma pour le cas<br />unidimensionnel utilise les Équations Différentielles Stochastiques<br />avec Temps Local. En effet en dimension un les processus /X/ sont<br />solutions de telles équations. On construit une grille sur la droite<br />réelle, qu'une bijection adéquate transforme en une grille uniforme<br />de pas /h/. Cette bijection permet de transformer /X/ en /Y/ qui se<br />comporte localement comme un Skew Brownian Motion, pour lequel on<br />connaît les probabilités de transition sur une grille uniforme, et le<br />temps moyen passé sur chaque cellule de cette grille. Une marche<br />aléatoire peut alors être construite, qui converge vers /X/ en racine<br />de /h/. Toujours dans le cas unidimensionnel on propose un deuxième<br />schéma plus général. On se donne une grille non uniforme sur la<br />droite réelle, dont les cellules ont une taille proportionnelle à<br />/h/. On montre qu'on peut relier les probabilités de transition de<br />/X/ sur cette grille, ainsi que le temps moyen passé par /X/ sur<br />chacune de ses cellules, à des solutions de problèmes d'EDP<br />elliptiques ad hoc. Une marche aléatoire en temps et en espace est<br />ainsi construite, qui permet d'approcher /X/ à nouveau en racine de<br />/h/. Ensuite on présente des pistes pour adapter cette dernière<br />approche au cas bidimensionnel et les problèmes que cela soulève.<br />Enfin on illustre par des exemples numériques les schémas étudiés.

In this report we will calculate the amplitude for a scalar-to-scalars (φ3 <img src="http://www.diva-portal.org/cgi-bin/mimetex.cgi?%5Crightarrow" /> φ2φ2) decay which involves a triangle loop. We compute the real and imaginary part of the amplitude separately and will argue that this is much more straightforward and practical in this case rather than having to deal with or worry about branch cuts of logarithms. We will derive simple cutting rules closely related to the imaginary part of the amplitude. In doing this, we derive a formula that deals with expressions of the form δ[f(x,y)]δ[g(x,y)], containing two Dirac delta functions.

We give a rather self-contained introduction to the rare pion to dielectron decay which in nontrivial leading order is given by a QED triangle loop. We work within the dispersive framework where the imaginary part of the amplitude is obtained via the Cutkosky rules. We derive these rules in detail. Using the twofold Mellin-Barnes representation for the pion transition form factor, we derive a simple expression for the branching ratio B(π0 <img src="http://www.diva-portal.org/cgi-bin/mimetex.cgi?%5Crightarrow" /> e+e-) which we then test for various models. In particular a more recent form factor derived from a Lagrangian for light pseudoscalars and vector mesons inspired by effective field theories. Comparison with the KTeV experiment at Fermilab is made and we find that we are more than 3σ below the KTeV experiment for some of the form factors. This is in agreement with other theoretical models, such as the Vector Meson Dominance model and the quark-loop model within the constituent-quark framework. But we also find that we can be in agreement with KTeV if we explore some freedom of the form factor not fixed by the low-energy Lagrangian.

>Magister Scientiae - MSc / We present various methods of pricing Asian options. The methods include Monte Carlo simulations designed using control and antithetic variates, numerical solution of partial differential equation and using lower bounds.The price of the Asian option is known to be a certain risk-neutral expectation. Using the Feynman-Kac theorem, we deduce that the problem of determining the expectation implies solving a linear parabolic partial differential equation. This partial differential equation does not admit explicit solutions due to the fact that the distribution of a sum of lognormal variables is not explicit. We then solve the partial differential equation numerically using finite difference and Monte Carlo methods.Our Monte Carlo approach is based on the pseudo random numbers and not deterministic sequence of numbers on which Quasi-Monte Carlo methods are designed. To make the Monte Carlo method more effective, two variance reduction techniques are discussed.Under the finite difference method, we consider explicit and the Crank-Nicholson’s schemes. We demonstrate that the explicit method gives rise to extraneous solutions because the stability conditions are difficult to satisfy. On the other hand, the Crank-Nicholson method is unconditionally stable and provides correct solutions. Finally, we apply the pricing methods to a similar problem of determining the price of a European-style arithmetic basket option under the Black-Scholes framework. We find the optimal lower bound, calculate it numerically and compare this with those obtained by the Monte Carlo and Moment Matching methods.Our presentation here includes some of the most recent advances on Asian options, and we contribute in particular by adding detail to the proofs and explanations. We also contribute some novel numerical methods. Most significantly, we include an original contribution on the use of very sharp lower bounds towards pricing European basket options.

Neste trabalho estudamos modelos teóricos que descrevem sistemas eletrônicos fortemente correlacionados, em especial o modelo t-J, e suas aplicações a compostos de óxidos de cobre, notadamente os compostos que apresentam supercondutividade de alta temperatura crítica e o composto Sr2CuO2Cl2. No primeiro capítulo do trabalho, fazemos uma exposição de três modelos que envolvem o tratamento das interações elétron-elétron, que são os modelos de Hubbard de uma banda, o modelo de Heisenberg e o modelo t-J. Na dedução deste último fazemos uma expansão canônica do hamiltoniano de Hubbard, no limite de acoplamento forte, levando-nos a obter um novo hamiltoniano que pode ser utilizado para descrever um sistema antiferromagnético bidimensional na presen- ça de lacunas, que é exatamente o que caracteriza os compostos supercondutores de alta temperatura crítica na sua fase de baixa dopagem.Após termos obtido o hamiltoniano que descreve o modelo t-J, aplicamos à este uma descrição de polarons de spin, numa representação de holons, que são férmions sem spin, e spinons, que são bósons que carregam somente os graus de liberdade de spin. Utilizando uma função de Green para descrever a propagação do polaron pela rede, obtemos uma equação para a sua autoenergia somando uma série de diagramas de Feynman, sendo que para este cálculo utilizamos a aproxima ção de Born autoconsistente[1]. Do ponto de vista numérico demonstramos que a equação integral de Dyson resultante do tratamento anterior não requer um procedimento iterativo para sua solução, e com isto conseguimos trabalhar com sistemas com grande número de partículas. Os resultados mostram, como um aspecto novo, que o tempo de vida média do holon tem um valor bastante grande no ponto (π,0 ) da rede recíproca, perto da singularidade de Van Hove mencionada na literatura[2]. Este aspecto, e suas implicações, é amplamente discutido neste capítulo. No capítulo 3 estudamos o modelo estendido t-t'-J, com tunelamento à segundos vizinhos e a incorporação dos termos de três sítios[3]. Fazemos a mesma formulação do capítulo anterior, e discutimos as aplicações dos nossos resultados ao óxido mencionado anteriormente. Finalmente, no último capítulo apresentamos uma aplicação original do modelo t-J à uma rede retangular, levemente distorcida, e demonstramos que os resultados do capítulo 3 são reproduzidos sem necessidade de introduzir termos de tunelamento adicionais no hamiltoniano. Esta aplicação pode se tornar relevante para o estudo das fases de tiras encontradas recentemente nesses materiais de óxidos de cobre.

Neste trabalho estudamos modelos teóricos que descrevem sistemas eletrônicos fortemente correlacionados, em especial o modelo t-J, e suas aplicações a compostos de óxidos de cobre, notadamente os compostos que apresentam supercondutividade de alta temperatura crítica e o composto Sr2CuO2Cl2. No primeiro capítulo do trabalho, fazemos uma exposição de três modelos que envolvem o tratamento das interações elétron-elétron, que são os modelos de Hubbard de uma banda, o modelo de Heisenberg e o modelo t-J. Na dedução deste último fazemos uma expansão canônica do hamiltoniano de Hubbard, no limite de acoplamento forte, levando-nos a obter um novo hamiltoniano que pode ser utilizado para descrever um sistema antiferromagnético bidimensional na presen- ça de lacunas, que é exatamente o que caracteriza os compostos supercondutores de alta temperatura crítica na sua fase de baixa dopagem.Após termos obtido o hamiltoniano que descreve o modelo t-J, aplicamos à este uma descrição de polarons de spin, numa representação de holons, que são férmions sem spin, e spinons, que são bósons que carregam somente os graus de liberdade de spin. Utilizando uma função de Green para descrever a propagação do polaron pela rede, obtemos uma equação para a sua autoenergia somando uma série de diagramas de Feynman, sendo que para este cálculo utilizamos a aproxima ção de Born autoconsistente[1]. Do ponto de vista numérico demonstramos que a equação integral de Dyson resultante do tratamento anterior não requer um procedimento iterativo para sua solução, e com isto conseguimos trabalhar com sistemas com grande número de partículas. Os resultados mostram, como um aspecto novo, que o tempo de vida média do holon tem um valor bastante grande no ponto (π,0 ) da rede recíproca, perto da singularidade de Van Hove mencionada na literatura[2]. Este aspecto, e suas implicações, é amplamente discutido neste capítulo. No capítulo 3 estudamos o modelo estendido t-t'-J, com tunelamento à segundos vizinhos e a incorporação dos termos de três sítios[3]. Fazemos a mesma formulação do capítulo anterior, e discutimos as aplicações dos nossos resultados ao óxido mencionado anteriormente. Finalmente, no último capítulo apresentamos uma aplicação original do modelo t-J à uma rede retangular, levemente distorcida, e demonstramos que os resultados do capítulo 3 são reproduzidos sem necessidade de introduzir termos de tunelamento adicionais no hamiltoniano. Esta aplicação pode se tornar relevante para o estudo das fases de tiras encontradas recentemente nesses materiais de óxidos de cobre.

Made available in DSpace on 2014-10-09T12:53:07Z (GMT). No. of bitstreams: 0 / Made available in DSpace on 2014-10-09T13:59:04Z (GMT). No. of bitstreams: 0 / Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP) / Tese (Doutoramento) / IPEN/T / Instituto de Pesquisas Energeticas e Nucleares - IPEN/CNEN-SP / FAPESP:03/01261-0

Neste trabalho estudamos modelos teóricos que descrevem sistemas eletrônicos fortemente correlacionados, em especial o modelo t-J, e suas aplicações a compostos de óxidos de cobre, notadamente os compostos que apresentam supercondutividade de alta temperatura crítica e o composto Sr2CuO2Cl2. No primeiro capítulo do trabalho, fazemos uma exposição de três modelos que envolvem o tratamento das interações elétron-elétron, que são os modelos de Hubbard de uma banda, o modelo de Heisenberg e o modelo t-J. Na dedução deste último fazemos uma expansão canônica do hamiltoniano de Hubbard, no limite de acoplamento forte, levando-nos a obter um novo hamiltoniano que pode ser utilizado para descrever um sistema antiferromagnético bidimensional na presen- ça de lacunas, que é exatamente o que caracteriza os compostos supercondutores de alta temperatura crítica na sua fase de baixa dopagem.Após termos obtido o hamiltoniano que descreve o modelo t-J, aplicamos à este uma descrição de polarons de spin, numa representação de holons, que são férmions sem spin, e spinons, que são bósons que carregam somente os graus de liberdade de spin. Utilizando uma função de Green para descrever a propagação do polaron pela rede, obtemos uma equação para a sua autoenergia somando uma série de diagramas de Feynman, sendo que para este cálculo utilizamos a aproxima ção de Born autoconsistente[1]. Do ponto de vista numérico demonstramos que a equação integral de Dyson resultante do tratamento anterior não requer um procedimento iterativo para sua solução, e com isto conseguimos trabalhar com sistemas com grande número de partículas. Os resultados mostram, como um aspecto novo, que o tempo de vida média do holon tem um valor bastante grande no ponto (π,0 ) da rede recíproca, perto da singularidade de Van Hove mencionada na literatura[2]. Este aspecto, e suas implicações, é amplamente discutido neste capítulo. No capítulo 3 estudamos o modelo estendido t-t'-J, com tunelamento à segundos vizinhos e a incorporação dos termos de três sítios[3]. Fazemos a mesma formulação do capítulo anterior, e discutimos as aplicações dos nossos resultados ao óxido mencionado anteriormente. Finalmente, no último capítulo apresentamos uma aplicação original do modelo t-J à uma rede retangular, levemente distorcida, e demonstramos que os resultados do capítulo 3 são reproduzidos sem necessidade de introduzir termos de tunelamento adicionais no hamiltoniano. Esta aplicação pode se tornar relevante para o estudo das fases de tiras encontradas recentemente nesses materiais de óxidos de cobre.

Made available in DSpace on 2014-10-09T12:53:07Z (GMT). No. of bitstreams: 0 / Made available in DSpace on 2014-10-09T13:59:04Z (GMT). No. of bitstreams: 0 / Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP) / Uma nova metodologia para a medida absoluta da fração efetiva de nêutrons atrasados &beta;eff, baseada em técnicas de análise de ruído microscópico e no modelo de Duas- Regiões, foi desenvolvida no reator IPEN/MB-01. Diferentemente das demais técnicas, tais como o Método de Bennet Modificado, o Método do Número de Nelson e o Método da fonte de 252Cf, a principal vantagem da metodologia proposta é a obtenção de &beta;eff de um modo puramente experimental, sem a necessidade de quaisquer outros parâmetros, sejam estes calculados ou provenientes de outros experimentos. Com a finalidade de validar este novo método, uma série de experimentos Rossi-&alpha; e Feynman-&alpha; foram realizados no reator IPEN/MB-01. De acordo com a metodologia proposta, &beta;eff foi estimado com uma incerteza de 0.67%, a qual atende aos requisitos de precisão almejados. Além disso, o tempo de geração de nêutrons prontos , dentre outros parâmetros, também foi obtido experimentalmente via esta metodologia. Em geral, os parâmetros medidos estão em acordo com resultados provenientes de experimentos de análise de ruído macroscópico. Nas comparações teoria-experimento, os valores de &beta;eff medidos neste trabalho mostram que a biblioteca JENDL3.3 apresenta a melhor performance (dentro de 1%). Esta concordância justifica a redução no yield de fissão do 235U proposta por Sakurai e Okajima. / Tese (Doutoramento) / IPEN/T / Instituto de Pesquisas Energeticas e Nucleares - IPEN/CNEN-SP / FAPESP:03/01261-0

Submitted by Pedro Guimarães (pedroengel@hotmail.com) on 2012-10-03T00:30:05Z No. of bitstreams: 1 Dissertação.pdf: 709819 bytes, checksum: 11a410ca7be737e1b89ed81a75cd3a0d (MD5) / Approved for entry into archive by Marcia Bacha (marcia.bacha@fgv.br) on 2012-12-20T16:42:02Z (GMT) No. of bitstreams: 1 Dissertação.pdf: 709819 bytes, checksum: 11a410ca7be737e1b89ed81a75cd3a0d (MD5) / Made available in DSpace on 2012-12-20T16:42:10Z (GMT). No. of bitstreams: 1 Dissertação.pdf: 709819 bytes, checksum: 11a410ca7be737e1b89ed81a75cd3a0d (MD5) Previous issue date: 2012-06-29 / I present here an approach that unify a variety of derivative pricing models that consists of attaining a Partial Differential Equation(PDE) by intuitive arguments and give its solution by Feynman-Kac method as a conditinal expectation of a markovian process. The expectation is taken in a risk neutral world(or risk neutral measure) where all the assets grow at the risk free rate. I also present how to make this specific change of measure, connecting the real world to the risk neutral world, and show that the relevant element for the measure change is the market price of factor risk. When the market is complete the market price of risk is unique and when the market is incomplete there is a variety of possible prices to the market price of factor risks that satisfy no arbitrage arguments. In the latter case the parameters are usually chosen to calibrate the model to market data. / Apresento aqui uma abordagem que unifica a literatura sobre os vários modelos de apreçamento de derivativos que consiste em obter por argumentos intuitivos de não arbitragem uma Equação Diferencial Parcial(EDP) e através do método de Feynman-Kac uma solução que é representada por uma esperança condicional de um processo markoviano do preço do derivativo descontado pela taxa livre de risco. Por este resultado, temos que a esperança deve ser tomada com relação a processos que crescem à taxa livre de risco e por este motivo dizemos que a esperança é tomada em um mundo neutro ao risco(ou medida neutra ao risco). Apresento ainda como realizar uma mudança de medida pertinente que conecta o mundo real ao mundo neutro ao risco e que o elemento chave para essa mudança de medida é o preço de mercado dos fatores de risco. No caso de mercado completo o preço de mercado do fator de risco é único e no caso de mercados incompletos existe uma variedade de preços aceitáveis para os fatores de risco pelo argumento de não arbitragem. Neste último caso, os preços de mercado são geralmente escolhidos de forma a calibrar o modelo com os dados de mercado.