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.

Study of multiferroic materials

The present PhD work aims the research and development of materials that exhibit multiferroic properties, in particular having a significant interaction between ferromagnetism and ferroelectricity; either directly within an intrinsic single phase or by combining extrinsic materials, achieving the coupling of properties through mechanic phenomena of the respective magnetostriction and piezoelectricity. These hybrid properties will allow the cross modification of magnetic and electric polarization states by the application of cross external magnetic and/or electric fields, giving way to a vast area for scientific investigation and potential technological applications in a new generation of electronic devices, such as computer memories, signal processing, transducers, sensors, etc. Initial experimental work consisted in chemical synthesis of nano powders oxides by urea pyrolysis method: A series of ceramic bulk composites with potential multiferroic properties comprised: of LuMnO3 with La0.7Sr0.3MnO3 and BaTiO3 with La0.7Ba0.3MnO3; and a series based on the intrinsic multiferroic LuMn1-zO3 phase modified with of Manganese vacancies. The acquisition of a new magnetron RF sputtering deposition system, in the Physics Department of Aveiro University, contributed to the proposal of an analogous experimental study in multiferroic thin films and multilayer samples. Besides the operational debut of this equipment several technical upgrades were completed like: the design and construction of the heater electrical contacts; specific shutters and supports for the magnetrons and for the substrate holder and; the addition of mass flow controllers, which allowed the introduction of N2 or O2 active atmosphere in the chamber; and the addition of a second RF generator, enabling co-deposition of different targets. Base study of the deposition conditions and resulting thin films characteristics in different substrates was made from an extensive list of targets. Particular attention was given to thin film deposition of magnetic phases La1-xSrxMnO3, La1-xBaxMnO3 and Ni2+x-yMn1-xGa1+y alloy, from the respective targets: La0.7Sr0.3MnO3, La0.7Ba0.3MnO3; and NiGa with NiMn. Main structural characterization of samples was performed by conventional and high resolution X-Ray Diffraction (XRD); chemical composition was determined by Electron Dispersion Spectroscopy (EDS); magnetization measurements recur to a Vibrating Sample Magnetometer (VSM) prototype; and surface probing (SPM) using Magnetic-Force (MFM) and Piezo-Response (PFM) Microscopy. Results clearly show that the composite bulk samples (LuM+LSM and BTO+LBM) feat the intended quality objectives in terms of phase composition and purity, having spurious contents below 0.5 %. SEM images confirm compact grain packaging and size distribution around the 50 nm scale. Electric conductivity, magnetization intensity and magneto impedance spreading response are coherent with the relative amount of magnetic phase in the sample. The existence of coupling between the functional phases is confirmed by the Magnetoelectric effect measurements of the sample “78%LuM+22%LSM” reaching 300% of electric response for 1 T at 100 kHz; while in the “78%BTO+22%LBM” sample the structural transitions of the magnetic phase at ~350 K result in a inversion of ME coefficient the behavior. A functional Magneto-Resistance measurement system was assembled from the concept stage until the, development and operational status; it enabled to test samples from 77 to 350 K, under an applied magnetic field up to 1 Tesla with 360º horizontal rotation; this system was also designed to measure Hall effect and has the potential to be further upgraded. Under collaboration protocols established with national and international institutions, complementary courses and sample characterization studies were performed using Magneto-Resistance (MR), Magneto-Impedance (MZ) and Magneto-Electric (ME) measurements; Raman and X-ray Photoelectron Spectroscopy (XPS); SQUID and VSM magnetization; Scanning Electron Microscopy (SEM) and Rutherford Back Scattering (RBS); Scan Probe Microscopy (SPM) with Band Excitation Probe Spectroscopy (BEPS); Neutron Powder Diffraction (NPD) and Perturbed Angular Correlations (PAC). Additional collaboration in research projects outside the scope of multiferroic materials provided further experience in sample preparation and characterization techniques, namely VSM and XPS measurements were performed in cubane molecular complex compounds and enable to identify the oxidation state of the integrating cluster of Ru ions; also, XRD and EDS/SEM analysis of the acquired targets and substrates implied the devolution of some items not in conformity with the specifications. Direct cooperation with parallel research projects regarding multiferroic materials, enable the assess to supplementary samples, namely a preliminary series of nanopowder Y1-x-yCaxØyMn1O3 and of Eu0.8Y0.2MnO3, a series of micropowder composites of LuMnO3 with La0.625Sr0.375MnO3 and of BaTiO3 with hexagonal ferrites; mono and polycrystalline samples of Pr1-xCaxMnO3, La1-xSrxMnO3 and La1-xCaxMnO3.