Simulações computacionais são ferramentas essenciais para a compreensão num nível microscópico de diversos fenômenos que ocorrem na natureza. Em particular, simulações ab initio, isto é, de primeiros princípios, podem predizer novas propriedades e auxiliar na interpretação de resultados experimentais, sem a necessidade de potenciais empíricos os quais são ajustados para uma específica configuração do sistema. No entanto, mesmo para simulações ab initio são necessárias aproximações tanto para o cálculo de estrutura eletrônica quanto para a descrição dos movimentos nucleares. Nesta tese, novas aproximações para o funcional de energia de troca e correlação da Teoria do Funcional da Densidade (DFT) são propostas e testadas no cálculo de estrutura eletrônica de átomos, moléculas e sólidos. Para a descrição dos movimentos nucleares a técnica de Monte Carlo é utilizada, porém as energias totais são obtidas com a DFT. Propomos também uma nova metodologia que permite descrever movimentos intramoleculares de sistemas cujas frequências de vibração não podem ser tratadas classicamente. Como aplicações estudamos aglomerados de água e água líquida, mostrando a relevância dessa nova metodologia na descrição de propriedades estruturais, vibracionais e de momento de dipolo desses sistemas. / Computer simulations are essential tools for a microscopic understanding of many processes that occur in nature. In particular, ab initio simulations, i.e., first principles simulations, can predict new properties and support experimental results, without the need to use empirical potentials which are fitted for a specific configuration of the system. However, even in ab initio simulations it is necessary to do approximations for the electronic structure calculations and for the description of the nuclear movements. In this thesis, new approximations for the exchange-correlation energy functional of Density Functional Theory (DFT) are proposed and tested for atoms, molecules and solids. The description of the nuclear movements was done using the Monte Carlo technique, however the total energy calculations were obtained by DFT. We also propose a new methodology which allows the description of intramolecular movements whose vibrational frequencies can not be treated classically. As an illustration, we have studied water clusters and liquid water, showing the relevance of this new methodology on the description of structural, vibrational and of the dipole moment of these systems.

Simulações computacionais são ferramentas essenciais para a compreensão num nível microscópico de diversos fenômenos que ocorrem na natureza. Em particular, simulações ab initio, isto é, de primeiros princípios, podem predizer novas propriedades e auxiliar na interpretação de resultados experimentais, sem a necessidade de potenciais empíricos os quais são ajustados para uma específica configuração do sistema. No entanto, mesmo para simulações ab initio são necessárias aproximações tanto para o cálculo de estrutura eletrônica quanto para a descrição dos movimentos nucleares. Nesta tese, novas aproximações para o funcional de energia de troca e correlação da Teoria do Funcional da Densidade (DFT) são propostas e testadas no cálculo de estrutura eletrônica de átomos, moléculas e sólidos. Para a descrição dos movimentos nucleares a técnica de Monte Carlo é utilizada, porém as energias totais são obtidas com a DFT. Propomos também uma nova metodologia que permite descrever movimentos intramoleculares de sistemas cujas frequências de vibração não podem ser tratadas classicamente. Como aplicações estudamos aglomerados de água e água líquida, mostrando a relevância dessa nova metodologia na descrição de propriedades estruturais, vibracionais e de momento de dipolo desses sistemas. / Computer simulations are essential tools for a microscopic understanding of many processes that occur in nature. In particular, ab initio simulations, i.e., first principles simulations, can predict new properties and support experimental results, without the need to use empirical potentials which are fitted for a specific configuration of the system. However, even in ab initio simulations it is necessary to do approximations for the electronic structure calculations and for the description of the nuclear movements. In this thesis, new approximations for the exchange-correlation energy functional of Density Functional Theory (DFT) are proposed and tested for atoms, molecules and solids. The description of the nuclear movements was done using the Monte Carlo technique, however the total energy calculations were obtained by DFT. We also propose a new methodology which allows the description of intramolecular movements whose vibrational frequencies can not be treated classically. As an illustration, we have studied water clusters and liquid water, showing the relevance of this new methodology on the description of structural, vibrational and of the dipole moment of these systems.

Neste trabalho, nós investigamos o comportamento crítico dinâmico de três modelos estatísticos utilizando simulações Monte Carlo em tempos curtos. Inicialmente, estudamos os modelos tridimensionais de dupla-troca e de Heisenberg. O expoente dinâmico de persistência global, bem como o expoente z são estimados através de duas técnicas. Para obter o expoente de persistência global, aplicamos diretamente a lei de potência obtida para a probabilidade de persistência global e em seguida fizemos o colapso de uma função universal para duas redes de tamanhos diferentes. Para estimar o valor de z, nós usamos uma função mista que combina resultados de simulações realizadas com diferentes condições iniciais e o cumulante de Binder de quarta ordem dependente do tempo. O expoente dinâmico que governa o comportamento tipo lei de potência da magnetização inicial, é estimado através da correlação temporal da magnetização (modelos de dupla-troca e Heisenberg) e da aplicação direta de uma lei de potência (modelo de Heisenberg). Os expoentes estáticos da magnetização e comprimento de correlação são estimados seguindo o comportamento de escala do parâmetro de ordem e sua derivada, respectivamente. Os resultados confirmam que esses dois modelos pertencem à mesma classe de universalidade. Em seguida, alguns expoentes críticos dinâmicos e estáticos são estimados no ponto de bifurcação do modelo de spin com simetria Z(5) bidimensional. Neste ponto, o modelo apresenta dois parâmetros de ordem diferentes, cada um possuindo um conjunto diferente de índices críticos. Os valores dos expoentes críticos estáticos estão em boa concordância com os resultados exatos. Até onde sabemos, está é a primeira tentativa de se obter os expoentes críticos dinâmicos para os modelos de dupla troca, Heisenberg e para o modelo Z(5). / In this work, we investigate the dynamic critical behavior of three statistical models by using short-time Monte Carlo simulations. At first, we study the three-dimensional double-exchange and Heisenberg models. The global persistence exponent, as well as the exponent z are estimated through two techniques. The dynamical exponent of global persistence is obtained by using the straight application of the power law obtained for the global persistence probability and by following the scaling collapse of a universal function for two diferent lattice sizes. To estimate the value of z, we use a mixed function which combines results obtained from samples submitted to diferent initial configurations and the time dependent fourth-order Binder cumulant. The dynamical exponent which governs the power law behavior of the initial magnetization, is estimated through the time correlation of the magnetization (double-exchange and Heisenberg models) and through the straight application of a power law(Heisenberg model). The statical exponents of the magnetization and correlation length are estimated through the scaling behavior of the order parameter and its derivative, respectively. The results confirm which those models belong to the same universality class. Following, the dynamical exponents and the statical exponents are estimated at the bifurcation point of the two-dimensional Z(5)-symmetric spin model. In this point, the model presents two diferent order parameters, each one possessing a diferent set of critical indices. The values of the static critical exponents are in good agreement with the exact results. Our study is, to the best of our knowledge, the first attempt to obtain the dynamic critical exponents of the double-exchange, Heisenberg, and Z(5) models.

Conselho Nacional de Desenvolvimento Científico e Tecnológico / Nesta dissertação são apresentados resultados de simulações Monte Carlo de fluorescência de raios X (XRF), utilizando o programa GEANT4, para medidas de espessura de revestimento metálico (Ni e Zn) em base metálica (Fe). As simulações foram feitas para dois tamanhos de espessura para cada metal de revestimento, (5μm e 10μm), com passos de 0,1 μm e 0,001 μm e com 106 histórias. No cálculo da espessura do revestimento foram feitas as aproximações de feixe de raios X monoenegético, com a análise da transmissão apenas da energia do K-alfa e para uma geometria compatível com um sistema real de medição (ARTAX-200). Os resultados mostraram a eficiência da metodologia de simulação e do cálculo da espessura do revestimento, o que permitirá futuros cálculos, inclusive para multirevestimentos metálicos em base metálica. / This dissertation is presented results of x-ray fluorescence (XRF) Monte Carlo simulations) using GEANT4 for metallic coating (Ni and Zn) thickness determination on metallic base (Fe). The simulations were made for two values of coating thickness for each covering (5 μm and 10 μm), with steps of 0.1 μm and 0.001 μm, and with 106 stories. The monoenergetic x-ray approach was used assuming the transmission of only K-alpha line. The experiment and geometry was similar to the real system of measurement (ARTAX-200 ). The results showed the efficiency of the simulation methodology, as well as the ability to measure the thickness of coating by XRF. Thus, a continuation of this study for multilayer metal coating could be promising.

Conselho Nacional de Desenvolvimento Científico e Tecnológico / Nesta dissertação são apresentados resultados de simulações Monte Carlo de fluorescência de raios X (XRF), utilizando o programa GEANT4, para medidas de espessura de revestimento metálico (Ni e Zn) em base metálica (Fe). As simulações foram feitas para dois tamanhos de espessura para cada metal de revestimento, (5μm e 10μm), com passos de 0,1 μm e 0,001 μm e com 106 histórias. No cálculo da espessura do revestimento foram feitas as aproximações de feixe de raios X monoenegético, com a análise da transmissão apenas da energia do K-alfa e para uma geometria compatível com um sistema real de medição (ARTAX-200). Os resultados mostraram a eficiência da metodologia de simulação e do cálculo da espessura do revestimento, o que permitirá futuros cálculos, inclusive para multirevestimentos metálicos em base metálica. / This dissertation is presented results of x-ray fluorescence (XRF) Monte Carlo simulations) using GEANT4 for metallic coating (Ni and Zn) thickness determination on metallic base (Fe). The simulations were made for two values of coating thickness for each covering (5 μm and 10 μm), with steps of 0.1 μm and 0.001 μm, and with 106 stories. The monoenergetic x-ray approach was used assuming the transmission of only K-alpha line. The experiment and geometry was similar to the real system of measurement (ARTAX-200 ). The results showed the efficiency of the simulation methodology, as well as the ability to measure the thickness of coating by XRF. Thus, a continuation of this study for multilayer metal coating could be promising.

Neste trabalho, nós investigamos o comportamento crítico dinâmico de três modelos estatísticos utilizando simulações Monte Carlo em tempos curtos. Inicialmente, estudamos os modelos tridimensionais de dupla-troca e de Heisenberg. O expoente dinâmico de persistência global, bem como o expoente z são estimados através de duas técnicas. Para obter o expoente de persistência global, aplicamos diretamente a lei de potência obtida para a probabilidade de persistência global e em seguida fizemos o colapso de uma função universal para duas redes de tamanhos diferentes. Para estimar o valor de z, nós usamos uma função mista que combina resultados de simulações realizadas com diferentes condições iniciais e o cumulante de Binder de quarta ordem dependente do tempo. O expoente dinâmico que governa o comportamento tipo lei de potência da magnetização inicial, é estimado através da correlação temporal da magnetização (modelos de dupla-troca e Heisenberg) e da aplicação direta de uma lei de potência (modelo de Heisenberg). Os expoentes estáticos da magnetização e comprimento de correlação são estimados seguindo o comportamento de escala do parâmetro de ordem e sua derivada, respectivamente. Os resultados confirmam que esses dois modelos pertencem à mesma classe de universalidade. Em seguida, alguns expoentes críticos dinâmicos e estáticos são estimados no ponto de bifurcação do modelo de spin com simetria Z(5) bidimensional. Neste ponto, o modelo apresenta dois parâmetros de ordem diferentes, cada um possuindo um conjunto diferente de índices críticos. Os valores dos expoentes críticos estáticos estão em boa concordância com os resultados exatos. Até onde sabemos, está é a primeira tentativa de se obter os expoentes críticos dinâmicos para os modelos de dupla troca, Heisenberg e para o modelo Z(5). / In this work, we investigate the dynamic critical behavior of three statistical models by using short-time Monte Carlo simulations. At first, we study the three-dimensional double-exchange and Heisenberg models. The global persistence exponent, as well as the exponent z are estimated through two techniques. The dynamical exponent of global persistence is obtained by using the straight application of the power law obtained for the global persistence probability and by following the scaling collapse of a universal function for two diferent lattice sizes. To estimate the value of z, we use a mixed function which combines results obtained from samples submitted to diferent initial configurations and the time dependent fourth-order Binder cumulant. The dynamical exponent which governs the power law behavior of the initial magnetization, is estimated through the time correlation of the magnetization (double-exchange and Heisenberg models) and through the straight application of a power law(Heisenberg model). The statical exponents of the magnetization and correlation length are estimated through the scaling behavior of the order parameter and its derivative, respectively. The results confirm which those models belong to the same universality class. Following, the dynamical exponents and the statical exponents are estimated at the bifurcation point of the two-dimensional Z(5)-symmetric spin model. In this point, the model presents two diferent order parameters, each one possessing a diferent set of critical indices. The values of the static critical exponents are in good agreement with the exact results. Our study is, to the best of our knowledge, the first attempt to obtain the dynamic critical exponents of the double-exchange, Heisenberg, and Z(5) models.

Estudamos o diagrama de fases e o ponto tricrítico da transição de colapso em um modelo de animais na rede quadrada, a partir da expansão em série da compressibilidade isotérmica KT do sistema. Como função das variáveis x (fugacidade) e y = e1/T (T é a temperatura reduzida), a serie KT é analisada utilizando-se a técnica dos aproximantes diferenciais parciais. Determinamos o padrão de fluxo das trajetórias características de um típico aproximante diferencial parcial com ponto fixo estável. Obtemos estimativas satisfatórias para a fugacidade tricrítica Xt = 0.024 &#177 0.005 e a temperatura tricritica Tt = 0.54 &#177 0.04. Considerando somente campos de escala lineares, obtemos também o expoente de escala &#947 = 1.4 &#177 0.2 e o expoente \"crossover\" &#934 = 0.66 &#177 0.08. Nossos resultados estão em boa concordância com estimativas prévias obtidas por outros métodos. Também estudamos um processo de polimerização ramifIcada através de simulações computacionais na rede quadrada baseadas em um modelo de crescimento cinético generalizado para se incorporar ramifIcações e impurezas. A configuração do polímero e identificada com uma árvore-ligação (\"bond tree\") a fim de se examinar os aspectos topológicos. As dimensões fractais dos aglomerados (\"clusters\") são obtidas na criticalidade. As simulações também permitem o estudo da evolução temporal dos aglomerados bem como a determinação das auto-correlações temporais e expoentes críticos dinâmicos. Com relação aos efeitos de tamanho finito, uma técnica de cumulantes de quarta ordem e empregada para se estimar a probabilidade de ramificação critica bc e os expoentes críticos v e &#946. Na ausência de impurezas, a rugosidade da superfície e descrita em termos dos expoentes de Hurst. Finalmente, simulamos este modelo de crescimento cinético na rede quadrada utilizando um método de Monte Carlo para estudar a polimerização ramificada com interações atrativas de curto alcance entre os monômeros. O diagrama de fases que separa os regimes de crescimento finito e infinito e obtido no plano (T,b) (T é a temperatura reduzida e b é a probabilidade de ramificação). No limite termodinâmico, extrapolamos a temperatura T&#8727 = 0.102 &#177 0.005 abaixo da qual a fase e sempre infinita. Observamos também a ocorrência de uma transição de rugosidade na superfície do polímero. / The phase diagram and the tricritical point of a collapsing lattice animal are studied through an extended series expansion of the isothermal compressibility KT on a square lattice. As a function of the variables x (fugacity) and y = e1/T (T is the reduced temperature), this series KT is investigated using the partial differential approximants technique. The characteristic flow pattern of partial differential approximant trajectories is determined for a typical stable fixed point. We obtain satisfactory estimates for the tricritical fugacity Xt = 0.024 &#177 0.005and temperature Tt = 0.54 &#177 0.04.Taking into account only linear scaling fields we are also able to get the scaling exponent &#947 = 1.4 &#177 0.2 and the crossover exponent &#934 = 0.66 &#177 0.08. Our results are in good agreement with previous estimates from other methods. We also study ramified polymerization through computational simulations on the square lattice of a kinetic growth model generalized to incorporate branching and impurities. The polymer configuration is identified with a bond tree in order to examine its topology. The fractal dimensions of clusters are obtained at criticality. Simulations also allow the study of time evolution of clusters as well as the determination of time autocorrelations and dynamical critical exponents. In regard to finite size effects, a fourth-order cumulant technique is employed to estimate the critical branching probability be and the critical exponents v and &#946. In the absence of impurities, the surface roughness is described in terms of the Hurst exponents. Finally we simulate this kinetic growth model on the square lattice using a Monte Carlo approach in order to study ramified polymerization with short distance attractive interactions between monomers. The phase boundary separating finite from infinite growth regimes is obtained in the (T,b) space (T is the reduced temperature and b is the branching probability). In the thermodynamic limit, we extrapolate the temperature T = 0.102 &#177 0.005 below which the phase is found to be always infinite. We also observe the occurrence of a roughening transition at the polymer surface.

Fundação Carlos Chagas Filho de Amparo a Pesquisa do Estado do Rio de Janeiro / Nessa tese foi feito um estudo preliminar, destinado à elaboração do programa experimental inicial para a primeira instalação da tomografia com prótons (pCT) brasileira por meio de modelagem computacional. A terapia com feixe de prótons é uma forma bastante precisa de tratamento de câncer. Atualmente, o planejamento de tratamento é baseado na tomografia computadorizada com raios X, alternativamente, a tomografia com prótons pode ser usada. Algumas questões importantes, como efeito de escala e a Curva de Calibração (fonte de dados iniciais para planejamento de terapia com prótons), foram estudados neste trabalho. A passagem de prótons com energias iniciais de 19,68MeV; 23MeV; 25MeV; 49,10MeV e 230MeV pelas camadas de materiais variados (água, alumínio, polietileno, ouro) foi simulada usando códigos Monte Carlo populares como SRIM e GEANT4. Os resultados das simulações foram comparados com a previsão teórica (baseada na solução aproximada da equação de transporte de Boltzmann) e com resultados das simulações feitas com outro popular código Monte Carlo MCNPX. Análise comparativa dos resultados das simulações com dados experimentais publicados na literatura científica para alvos grossos e na faixa de energias de prótons usada em medidas em pCT foi feita. Foi observado que apesar de que todos os códigos mostram os resultados parecidos alguns deslocamentos não sistemáticos podem ser observados. Foram feitas observações importantes sobre a precisão dos códigos e uma necessidade em medidas sistemáticas de frenagem de prótons em alvos grossos foi declarada. / In the present work a preliminary research via computer simulations was made in order to elaborate a prior program for the first experimental pCT setup in Brazil. Proton therapy is a high precise form of a cancer treatment. Treatment planning nowadays is performed basing on X ray Computer Tomography data (CT), alternatively the same procedure could be performed using proton Computer Tomography (pCT). Some important questions, as a scale effect and so called Calibration Curve (as a source of primary data for pCT treatment planning) were studied in this work. The 19.68MeV; 23MeV; 25MeV; 49.10MeV e 230MeV protons passage through varied absorbers (water, aluminum, polyethylene, gold) were simulated by such popular Monte Carlo packages as SRIM and GEANT4. The simulation results were compared with a theoretic prevision based on approximate solution of the Boltzmann transport equation and with simulation results of the other popular Monte Carlo code MCNPX. The comparative analysis of the simulations results with the experimental data published in scientific literature for thick absorbers and within the energy range used in the pCT measurements was made. It was noted in spite of the fact that all codes showed similar results some nonsystematic displacements can be observed. Some important observations about the codes precision were made and a necessity of the systematic measurements of the proton stopping power in thick absorbers was declared.

Fundação Carlos Chagas Filho de Amparo a Pesquisa do Estado do Rio de Janeiro / Nessa tese foi feito um estudo preliminar, destinado à elaboração do programa experimental inicial para a primeira instalação da tomografia com prótons (pCT) brasileira por meio de modelagem computacional. A terapia com feixe de prótons é uma forma bastante precisa de tratamento de câncer. Atualmente, o planejamento de tratamento é baseado na tomografia computadorizada com raios X, alternativamente, a tomografia com prótons pode ser usada. Algumas questões importantes, como efeito de escala e a Curva de Calibração (fonte de dados iniciais para planejamento de terapia com prótons), foram estudados neste trabalho. A passagem de prótons com energias iniciais de 19,68MeV; 23MeV; 25MeV; 49,10MeV e 230MeV pelas camadas de materiais variados (água, alumínio, polietileno, ouro) foi simulada usando códigos Monte Carlo populares como SRIM e GEANT4. Os resultados das simulações foram comparados com a previsão teórica (baseada na solução aproximada da equação de transporte de Boltzmann) e com resultados das simulações feitas com outro popular código Monte Carlo MCNPX. Análise comparativa dos resultados das simulações com dados experimentais publicados na literatura científica para alvos grossos e na faixa de energias de prótons usada em medidas em pCT foi feita. Foi observado que apesar de que todos os códigos mostram os resultados parecidos alguns deslocamentos não sistemáticos podem ser observados. Foram feitas observações importantes sobre a precisão dos códigos e uma necessidade em medidas sistemáticas de frenagem de prótons em alvos grossos foi declarada. / In the present work a preliminary research via computer simulations was made in order to elaborate a prior program for the first experimental pCT setup in Brazil. Proton therapy is a high precise form of a cancer treatment. Treatment planning nowadays is performed basing on X ray Computer Tomography data (CT), alternatively the same procedure could be performed using proton Computer Tomography (pCT). Some important questions, as a scale effect and so called Calibration Curve (as a source of primary data for pCT treatment planning) were studied in this work. The 19.68MeV; 23MeV; 25MeV; 49.10MeV e 230MeV protons passage through varied absorbers (water, aluminum, polyethylene, gold) were simulated by such popular Monte Carlo packages as SRIM and GEANT4. The simulation results were compared with a theoretic prevision based on approximate solution of the Boltzmann transport equation and with simulation results of the other popular Monte Carlo code MCNPX. The comparative analysis of the simulations results with the experimental data published in scientific literature for thick absorbers and within the energy range used in the pCT measurements was made. It was noted in spite of the fact that all codes showed similar results some nonsystematic displacements can be observed. Some important observations about the codes precision were made and a necessity of the systematic measurements of the proton stopping power in thick absorbers was declared.

Estudamos o diagrama de fases e o ponto tricrítico da transição de colapso em um modelo de animais na rede quadrada, a partir da expansão em série da compressibilidade isotérmica KT do sistema. Como função das variáveis x (fugacidade) e y = e1/T (T é a temperatura reduzida), a serie KT é analisada utilizando-se a técnica dos aproximantes diferenciais parciais. Determinamos o padrão de fluxo das trajetórias características de um típico aproximante diferencial parcial com ponto fixo estável. Obtemos estimativas satisfatórias para a fugacidade tricrítica Xt = 0.024 &#177 0.005 e a temperatura tricritica Tt = 0.54 &#177 0.04. Considerando somente campos de escala lineares, obtemos também o expoente de escala &#947 = 1.4 &#177 0.2 e o expoente \"crossover\" &#934 = 0.66 &#177 0.08. Nossos resultados estão em boa concordância com estimativas prévias obtidas por outros métodos. Também estudamos um processo de polimerização ramifIcada através de simulações computacionais na rede quadrada baseadas em um modelo de crescimento cinético generalizado para se incorporar ramifIcações e impurezas. A configuração do polímero e identificada com uma árvore-ligação (\"bond tree\") a fim de se examinar os aspectos topológicos. As dimensões fractais dos aglomerados (\"clusters\") são obtidas na criticalidade. As simulações também permitem o estudo da evolução temporal dos aglomerados bem como a determinação das auto-correlações temporais e expoentes críticos dinâmicos. Com relação aos efeitos de tamanho finito, uma técnica de cumulantes de quarta ordem e empregada para se estimar a probabilidade de ramificação critica bc e os expoentes críticos v e &#946. Na ausência de impurezas, a rugosidade da superfície e descrita em termos dos expoentes de Hurst. Finalmente, simulamos este modelo de crescimento cinético na rede quadrada utilizando um método de Monte Carlo para estudar a polimerização ramificada com interações atrativas de curto alcance entre os monômeros. O diagrama de fases que separa os regimes de crescimento finito e infinito e obtido no plano (T,b) (T é a temperatura reduzida e b é a probabilidade de ramificação). No limite termodinâmico, extrapolamos a temperatura T&#8727 = 0.102 &#177 0.005 abaixo da qual a fase e sempre infinita. Observamos também a ocorrência de uma transição de rugosidade na superfície do polímero. / The phase diagram and the tricritical point of a collapsing lattice animal are studied through an extended series expansion of the isothermal compressibility KT on a square lattice. As a function of the variables x (fugacity) and y = e1/T (T is the reduced temperature), this series KT is investigated using the partial differential approximants technique. The characteristic flow pattern of partial differential approximant trajectories is determined for a typical stable fixed point. We obtain satisfactory estimates for the tricritical fugacity Xt = 0.024 &#177 0.005and temperature Tt = 0.54 &#177 0.04.Taking into account only linear scaling fields we are also able to get the scaling exponent &#947 = 1.4 &#177 0.2 and the crossover exponent &#934 = 0.66 &#177 0.08. Our results are in good agreement with previous estimates from other methods. We also study ramified polymerization through computational simulations on the square lattice of a kinetic growth model generalized to incorporate branching and impurities. The polymer configuration is identified with a bond tree in order to examine its topology. The fractal dimensions of clusters are obtained at criticality. Simulations also allow the study of time evolution of clusters as well as the determination of time autocorrelations and dynamical critical exponents. In regard to finite size effects, a fourth-order cumulant technique is employed to estimate the critical branching probability be and the critical exponents v and &#946. In the absence of impurities, the surface roughness is described in terms of the Hurst exponents. Finally we simulate this kinetic growth model on the square lattice using a Monte Carlo approach in order to study ramified polymerization with short distance attractive interactions between monomers. The phase boundary separating finite from infinite growth regimes is obtained in the (T,b) space (T is the reduced temperature and b is the branching probability). In the thermodynamic limit, we extrapolate the temperature T = 0.102 &#177 0.005 below which the phase is found to be always infinite. We also observe the occurrence of a roughening transition at the polymer surface.