Monte Carlo studies of phase transitions and adsorption : application to n-6 Lennard-Jones, C60 and zeolite systems

The solid-fluid transition properties of the n - 6 Lennard-Jones system are studied by means of extensive free energy calculations. Different values of the parameter n which regulates the steepness of the short-range repulsive interaction are investigated. Furthermore, the free energies of the n < 12 systems are calculated using the n = 12 system as a reference. The method relies on a generalization of the multiple histogram method that combines independent canonical ensemble simulations performed with different Hamiltonians and computes the free energy difference between them. The phase behavior of the fullerene C60 solid is studied by performing NPT simulations using atomistic models which treat each carbon in the molecule as a separate interaction site with additional bond charges. In particular, the transition from an orientationally frozen phase at low temperatures to one where the molecules are freely rotating at higher temperatures is studied as a function of applied pressure. The adsorption of molecular hydrogen in the zeolite NaA is investigated by means of grand-canonical Monte Carlo, in a wide range of temperatures and imposed gas pressures, and results are compared with available experimental data. A potential model is used that comprises three main interactions: van der Waals, Coulomb and induced polarization by the permanent electric field in the zeolite.

Phase transitions and nonlinear phenomena in neuronal network models

Communication and cooperation between billions of neurons underlie the power of the brain. How do complex functions of the brain arise from its cellular constituents? How do groups of neurons self-organize into patterns of activity? These are crucial questions in neuroscience. In order to answer them, it is necessary to have solid theoretical understanding of how single neurons communicate at the microscopic level, and how cooperative activity emerges. In this thesis we aim to understand how complex collective phenomena can arise in a simple model of neuronal networks. We use a model with balanced excitation and inhibition and complex network architecture, and we develop analytical and numerical methods for describing its neuronal dynamics. We study how interaction between neurons generates various collective phenomena, such as spontaneous appearance of network oscillations and seizures, and early warnings of these transitions in neuronal networks. Within our model, we show that phase transitions separate various dynamical regimes, and we investigate the corresponding bifurcations and critical phenomena. It permits us to suggest a qualitative explanation of the Berger effect, and to investigate phenomena such as avalanches, band-pass filter, and stochastic resonance. The role of modular structure in the detection of weak signals is also discussed. Moreover, we find nonlinear excitations that can describe paroxysmal spikes observed in electroencephalograms from epileptic brains. It allows us to propose a method to predict epileptic seizures. Memory and learning are key functions of the brain. There are evidences that these processes result from dynamical changes in the structure of the brain. At the microscopic level, synaptic connections are plastic and are modified according to the dynamics of neurons. Thus, we generalize our cortical model to take into account synaptic plasticity and we show that the repertoire of dynamical regimes becomes richer. In particular, we find mixed-mode oscillations and a chaotic regime in neuronal network dynamics.

Interacções intermoleculares fracas em materiais funcionais avançados: um estudo espectroscópico e computacional

A utilização combinada de espectroscopia vibracional e de cálculos envolvendo a teoria do funcional de densidade (DFT) possibilita o estudo de ligações de hidrogénio em fase condensada, assim como a análise da estrutura molecular dos sistemas em estudo. Por um lado, a espectroscopia vibracional permite a detecção de associações moleculares, enquanto os métodos computacionais auxiliam na obtenção de informação referente aos mecanismos de associação, nomeadamente no que diz respeito à possível estrutura de dímeros e compostos de inclusão em ciclodextrinas e às energias de interacção e de inclusão. O estudo que originou a presente dissertação pretende contribuir para o reforço da aplicação de estudos espectroscópicos e computacionais na elucidação de diversos fenómenos químicos, com especial destaque para o papel desempenhado por interacções intermoleculares fracas na estrutura e propriedades de materiais moleculares. No âmbito desta tese foram investigados os seguintes tópicos: polimorfismo e pseudopolimorfismo em sólidos farmacêuticos, transições de fase em misturas binárias de ácidos gordos, inclusão em ciclodextrinas, interacção de compostos farmacêuticos com superfícies metálicas e formação de agregados de água em materiais híbridos orgânicos-inorgânicos. Os sistemas foram analisados utilizando a espectroscopia vibracional – particularmente a espectroscopia de difusão de Raman – como técnica fundamental. Para uma melhor caracterização de processos envolvendo transições de fase, foram efectuados estudos com variação de temperatura, variação de humidade relativa e substituição isotópica. O estudo da interacção com superfícies metálicas foi realizado por espectroscopia de Raman intensificada à superfície. Dada a complexidade dos sistemas em estudo, a informação obtida por espectroscopia vibracional foi complementada por resultados de cálculos mecânico-quânticos. Em particular, os cálculos DFT foram utilizados para a optimização de geometrias e previsão de frequências vibracionais de moléculas e associações moleculares, permitindo assim a análise e interpretação de espectros vibracionais e a caracterização da estrutura de materiais.

Piezoresponse force microscopy of ferroelectric relaxors

Nesta tese, ferroeléctricos relaxor (I dont know uf the order is correct) de base Pb das familias (Pb,La)(Zr,Ti)O3 (PLZT), Pb(Mg1/3,Nb2/3)O3-PbTiO3 (PMN-PT), Pb(Zn1/3,Nb2/3)O3-PbTiO3 (PZN-PT) foram investigados e analisados. As propriedades ferroeléctricas e dieléctricas das amostras foram estudadas por métodos convencionais de macro e localmente por microscopia de força piezoeléctrica (PFM). Nos cerâmicos PLZT 9.75/65/35 o contraste da PFM à escala nanometrica _ foi investigado em função do tamanho e orientação dos grãos. Apurou-se que a intensidade do sinal piezoeléctrico das nanoestruturas diminui com o aumento da temperatura e desaparece a 490 K (La mol. 8%) e 420 K (9,5%). Os ciclos de histerese locais foram obtidos em função da temperatura. A evolução dos parâmetros macroscópicos e locais com a temperatura de superfície sugere um forte efeito de superfície nas transições de fase ferroeléctricas do material investigado. A rugosidade da parede de domínio é determinada por PFM para a estrutura de domínio natural existente neste ferroeléctrico policristalino. Além disso, os domínios ferroeléctricos artificiais foram criados pela aplicação de pulsos eléctricos à ponta do condutor PFM e o tamanho de domínio in-plane foi medido em função da duração do pulso. Todas estas experiências levaram à conclusão de que a parede de domínio em relaxors do tipo PZT é quase uma interface unidimensional. O mecanismo de contraste na superfície de relaxors do tipo PLZT é medido por PFMAs estruturas de domínio versus evolução da profundidade foram estudadas em cristais PZN-4,5%PT, com diferentes orientações através da PFM. Padrões de domínio irregulares com tamanhos típicos de 20-100 nm foram observados nas superfícies com orientação <001> das amostras unpoled?. Pelo contrário, os cortes de cristal <111> exibem domínios regulares de tamanho mícron normal, com os limites do domínio orientados ao longo dos planos cristalográficos permitidos. A existência de nanodomínios em cristais com orientação <001> está provisoriamente (wrong Word) atribuída à natureza relaxor de PZN-PT, onde pequenos grupos polares podem formar-se em coindições de zero-field-cooling (ZFC). Estes nanodomínios são considerados como os núcleos do estado de polarização oposta e podem ser responsáveis pelo menor campo coercitivo para este corte de cristal em particular. No entanto, a histerese local piezoelétrica realizada pelo PFM à escala nanométrica indica uma mudança de comportamento de PZN-PT semelhante para ambas as orientações cristalográficas investigadas. A evolução das estruturas de domínio com polimento abaixo da superfície do cristal foi investigada. O domínio de ramificações e os efeitos de polarização de triagem após o polimento e as medições de temperatura têm sido estudados pela PFM e pela análise SEM. Além disso, verificou-se que a intensidade do sinal piezoeléctrico a partir das estruturas de nanodomínio diminui com o aumento da temperatura, acabando por desaparecer aos 430 K (orientaçáo <111>) e 470 K (orientação <100>). Esta diferença de temperatura nas transições de fase local em cristais de diferentes orientações é explicada pelo forte efeito de superfície na transição da fase ferroelétrica em relaxors.A comutação da polarização em relaxor ergódico e nas fases ferroeléctricas do sistema PMN-PT foram realizadas pela combinação de três métodos, Microscopia de Força Piezoeléctrica, medição de um único ponto de relaxamento eletromecânico e por ultimo mapeamento de espectroscopia de tensão. A dependência do comportamento do relaxamento na amplitude e tempo da tensão de pulso foi encontrada para seguir um comportamento logarítmico universal com uma inclinação quase constante. Este comportamento é indicativo da progressiva população dos estados de relaxamento lento, ao contrário de uma relaxação linear na presença de uma ampla distribuição do tempo de relaxamento. O papel do comportamento de relaxamento, da não-linearidade ferroeléctrica e da heterogeneidade espacial do campo na ponta da sonda de AFM sobre o comportamento do ciclo de histerese é analisada em detalhe. Os ciclos de histerese para ergódica PMN- 10%PT são mostrados como cineticamente limitados, enquanto que no PMN, com maior teor de PT, são observados verdadeiros ciclos de histerese ferroeléctrica com viés de baixa nucleação.

Gas liquid equilibria at high pressures

Este trabalho tem como objectivo a medição da solubilidade de gases em líquidos a altas pressões. Para a realização experimental foi usada uma célula de altas pressões de volume variável e um método de observação directa para a detecção das transições de fase. As condições usadas neste trabalho experimental foram temperaturas entre 276.00 e 373.15 K e pressões ate 100 MPa. Numa primeira fase foram efectuadas medições de solubilidade de metano em anilina pura e para uma solução aquosa de composição mássica igual a 5% em anilina. A anilina é fundamentalmente usada em processos como a produção de borracha, de poliuretanos, pigmentos e tintas, fármacos, herbicidas e fungicidas. É usualmente produzida através de processos de redução do nitrobenzeno por reacção com o HCl. De facto, esta síntese é um processo químico complexo onde um grande número de processos compete entre si. Alterações nas condições do processo com a possível formação de produtos intermediários podem afectar a eficiência do processo. Para a sua melhoria foi sugerido que o hidrogénio usado fosse simultaneamente produzido e gasto no reactor principal. Neste caso e para a produção do mesmo, era necessário adicionar metano e água ao reactor. Tendo por base a ideia de que as reacções onde o hidrogénio é reagente e produto ocorrem em simultâneo, torna-se clara a importância do estudo da solubilidade do metano em anilina pura e em soluções aquosas desta. Numa segunda fase foi estudada a solubilidade do dióxido de carbono em soluções aquosas de tri-iso-butil(metil)fosfónio tosilato, com composições molares são de 4, 8 e 12% em líquido iónico. Este pertencente à família dos fosfónios. Possui uma viscosidade e densidade elevadas, é térmica e quimicamente estável e ainda possui uma elevada polaridade. Apresenta uma miscibilidade completa em água e nos solventes mais usuais, como o diclorometano e tolueno, não sendo no entanto míscivel em hexano. O tri-isobutil( metil)fosfónium tosilato é usado como solvente nos processos de hidroformilação de olefinas e ainda em processos de captura e conversão de dióxido de carbono. Neste trabalho experimental, a temperatura e a pressão foram inicialmente aumentadas até o sistema atingir o equilíbrio. A pressão é diminuída lentamente até se verificar o aparecimento/desaparecimento da última bolha de gás. A pressão à qual a última bolha de gás desaparece representa a pressão de equilíbrio para aquela temperatura. Este procedimento foi efectuado para vários sistemas e várias temperaturas.

Novel critical phenomena in optimization driven processes on networks

The work presented in this Ph.D thesis was developed in the context of complex network theory, from a statistical physics standpoint. We examine two distinct problems in this research field, taking a special interest in their respective critical properties. In both cases, the emergence of criticality is driven by a local optimization dynamics. Firstly, a recently introduced class of percolation problems that attracted a significant amount of attention from the scientific community, and was quickly followed up by an abundance of other works. Percolation transitions were believed to be continuous, until, recently, an 'explosive' percolation problem was reported to undergo a discontinuous transition, in [93]. The system's evolution is driven by a metropolis-like algorithm, apparently producing a discontinuous jump on the giant component's size at the percolation threshold. This finding was subsequently supported by number of other experimental studies [96, 97, 98, 99, 100, 101]. However, in [1] we have proved that the explosive percolation transition is actually continuous. The discontinuity which was observed in the evolution of the giant component's relative size is explained by the unusual smallness of the corresponding critical exponent, combined with the finiteness of the systems considered in experiments. Therefore, the size of the jump vanishes as the system's size goes to infinity. Additionally, we provide the complete theoretical description of the critical properties for a generalized version of the explosive percolation model [2], as well as a method [3] for a precise calculation of percolation's critical properties from numerical data (useful when exact results are not available). Secondly, we study a network flow optimization model, where the dynamics consists of consecutive mergings and splittings of currents flowing in the network. The current conservation constraint does not impose any particular criterion for the split of current among channels outgoing nodes, allowing us to introduce an asymmetrical rule, observed in several real systems. We solved analytically the dynamic equations describing this model in the high and low current regimes. The solutions found are compared with numerical results, for the two regimes, showing an excellent agreement. Surprisingly, in the low current regime, this model exhibits some features usually associated with continuous phase transitions.

Random path sampling applied to vector models of magnetism

The recently reported Monte Carlo Random Path Sampling method (RPS) is here improved and its application is expanded to the study of the 2D and 3D Ising and discrete Heisenberg models. The methodology was implemented to allow use in both CPU-based high-performance computing infrastructures (C/MPI) and GPU-based (CUDA) parallel computation, with significant computational performance gains. Convergence is discussed, both in terms of free energy and magnetization dependence on field/temperature. From the calculated magnetization-energy joint density of states, fast calculations of field and temperature dependent thermodynamic properties are performed, including the effects of anisotropy on coercivity, and the magnetocaloric effect. The emergence of first-order magneto-volume transitions in the compressible Ising model is interpreted using the Landau theory of phase transitions. Using metallic Gadolinium as a real-world example, the possibility of using RPS as a tool for computational magnetic materials design is discussed. Experimental magnetic and structural properties of a Gadolinium single crystal are compared to RPS-based calculations using microscopic parameters obtained from Density Functional Theory.