Applying many-body expansion aiming to make ab in initio molecular dynamics more efficient

In this present work we present a methodology that aims to apply the many-body expansion to decrease the computational cost of ab initio molecular dynamics, keeping acceptable accuracy on the results. We implemented this methodology in a program which we called ManBo. In the many-body expansion approach, we partitioned the total energy E of the system in contributions of one body, two bodies, three bodies, etc., until the contribution of the Nth body [1-3]: E = E1 + E2 + E3 + …EN. The E1 term is the sum of the internal energy of the molecules; the term E2 is the energy due to interaction between all pairs of molecules; E3 is the energy due to interaction between all trios of molecules; and so on. In Manbo we chose to truncate the expansion in the contribution of two or three bodies, both for the calculation of the energy and for the calculation of the atomic forces. In order to partially include the many-body interactions neglected when we truncate the expansion, we can include an electrostatic embedding in the electronic structure calculations, instead of considering the monomers, pairs and trios as isolated molecules in space. In simulations we made we chose to simulate water molecules, and use the Gaussian 09 as external program to calculate the atomic forces and energy of the system, as well as reference program for analyzing the accuracy of the results obtained with the ManBo. The results show that the use of the many-body expansion seems to be an interesting approach for reducing the still prohibitive computational cost of ab initio molecular dynamics. The errors introduced on atomic forces in applying such methodology are very small. The inclusion of an embedding electrostatic seems to be a good solution for improving the results with only a small increase in simulation time. As we increase the level of calculation, the simulation time of ManBo tends to largely decrease in relation to a conventional BOMD simulation of Gaussian, due to better scalability of the methodology presented. References [1] E. E. Dahlke and D. G. Truhlar; J. Chem. Theory Comput., 3, 46 (2007). [2] E. E. Dahlke and D. G. Truhlar; J. Chem. Theory Comput., 4, 1 (2008). [3] R. Rivelino, P. Chaudhuri and S. Canuto; J. Chem. Phys., 118, 10593 (2003).

Study of surface free energy on a lattice-gas model applied to Liquid bridge

The pulmonary crackling and the formation of liquid bridges are problems that for centuries have been attracting the attention of scientists. In order to study these phenomena, it was developed a canonical cubic lattice-gas­ like model to explain the rupture of liquid bridges in lung airways [A. Alencar et al., 2006, PRE]. Here, we further develop this model and add entropy analysis to study thermodynamic properties, such as free energy and force. The simulations were performed using the Monte Carlo method with Metropolis algorithm. The exchange between gas and liquid particles were performed randomly according to the Kawasaki dynamics and weighted by the Boltzmann factor. Each particle, which can be solid (s), liquid (l) or gas (g), has 26 neighbors: 6 + 12 + 8, with distances 1, √2 and √3, respectively. The energy of a lattice's site m is calculated by the following expression: Em = ∑k=126 Ji(m)j(k) in witch (i, j) = g, l or s. Specifically, it was studied the surface free energy of the liquid bridge, trapped between two planes, when its height is changed. For that, was considered two methods. First, just the internal energy was calculated. Then was considered the entropy. It was fond no difference in the surface free energy between this two methods. We calculate the liquid bridge force between the two planes using the numerical surface free energy. This force is strong for small height, and decreases as the distance between the two planes, height, is increased. The liquid-gas system was also characterized studying the variation of internal energy and heat capacity with the temperature. For that, was performed simulation with the same proportion of liquid and gas particle, but different lattice size. The scale of the liquid-gas system was also studied, for low temperature, using different values to the interaction Jij.

Dynamics of photoexcited carriers in the presence of disorder in radial heterostructured GaAs/AlGaAs/GaAs nanowires

Time-resolved photoluminescence was employed to study electron-hole dynamics in radial heterostructured GaAs/AlGaAs/GaAs core/inner shell/outer shell nanowires. It was found that impurity random potential results in a red shift of the recombination time maximum with respect to the photoluminescence peak energy.

Towards the construction of a validated numerical system to study the mesoscale dynamics of the North East Atlantic (2003-2006)

Máster en Oceanografía

Numerical simulations of magma chamber dynamics at Campi Flegrei, and associated seismicity, deformation and gravity changes

Understanding the complex relationships between quantities measured by volcanic monitoring network and shallow magma processes is a crucial headway for the comprehension of volcanic processes and a more realistic evaluation of the associated hazard. This question is very relevant at Campi Flegrei, a volcanic quiescent caldera immediately north-west of Napoli (Italy). The system activity shows a high fumarole release and periodic ground slow movement (bradyseism) with high seismicity. This activity, with the high people density and the presence of military and industrial buildings, makes Campi Flegrei one of the areas with higher volcanic hazard in the world. In such a context my thesis has been focused on magma dynamics due to the refilling of shallow magma chambers, and on the geophysical signals detectable by seismic, deformative and gravimetric monitoring networks that are associated with this phenomenologies. Indeed, the refilling of magma chambers is a process frequently occurring just before a volcanic eruption; therefore, the faculty of identifying this dynamics by means of recorded signal analysis is important to evaluate the short term volcanic hazard. The space-time evolution of dynamics due to injection of new magma in the magma chamber has been studied performing numerical simulations with, and implementing additional features in, the code GALES (Longo et al., 2006), recently developed and still on the upgrade at the Istituto Nazionale di Geofisica e Vulcanologia in Pisa (Italy). GALES is a finite element code based on a physico-mathematical two dimensional, transient model able to treat fluids as multiphase homogeneous mixtures, compressible to incompressible. The fundamental equations of mass, momentum and energy balance are discretised both in time and space using the Galerkin Least-Squares and discontinuity-capturing stabilisation technique. The physical properties of the mixture are computed as a function of local conditions of magma composition, pressure and temperature.The model features enable to study a broad range of phenomenologies characterizing pre and sin-eruptive magma dynamics in a wide domain from the volcanic crater to deep magma feeding zones. The study of displacement field associated with the simulated fluid dynamics has been carried out with a numerical code developed by the Geophysical group at the University College Dublin (O’Brien and Bean, 2004b), with whom we started a very profitable collaboration. In this code, the seismic wave propagation in heterogeneous media with free surface (e.g. the Earth’s surface) is simulated using a discrete elastic lattice where particle interactions are controlled by the Hooke’s law. This method allows to consider medium heterogeneities and complex topography. The initial and boundary conditions for the simulations have been defined within a coordinate project (INGV-DPC 2004-06 V3_2 “Research on active volcanoes, precursors, scenarios, hazard and risk - Campi Flegrei”), to which this thesis contributes, and many researchers experienced on Campi Flegrei in volcanological, seismic, petrological, geochemical fields, etc. collaborate. Numerical simulations of magma and rock dynamis have been coupled as described in the thesis. The first part of the thesis consists of a parametric study aimed at understanding the eect of the presence in magma of carbon dioxide in magma in the convection dynamics. Indeed, the presence of this volatile was relevant in many Campi Flegrei eruptions, including some eruptions commonly considered as reference for a future activity of this volcano. A set of simulations considering an elliptical magma chamber, compositionally uniform, refilled from below by a magma with volatile content equal or dierent from that of the resident magma has been performed. To do this, a multicomponent non-ideal magma saturation model (Papale et al., 2006) that considers the simultaneous presence of CO2 and H2O, has been implemented in GALES. Results show that the presence of CO2 in the incoming magma increases its buoyancy force promoting convection ad mixing. The simulated dynamics produce pressure transients with frequency and amplitude in the sensitivity range of modern geophysical monitoring networks such as the one installed at Campi Flegrei . In the second part, simulations more related with the Campi Flegrei volcanic system have been performed. The simulated system has been defined on the basis of conditions consistent with the bulk of knowledge of Campi Flegrei and in particular of the Agnano-Monte Spina eruption (4100 B.P.), commonly considered as reference for a future high intensity eruption in this area. The magmatic system has been modelled as a long dyke refilling a small shallow magma chamber; magmas with trachytic and phonolitic composition and variable volatile content of H2O and CO2 have been considered. The simulations have been carried out changing the condition of magma injection, the system configuration (magma chamber geometry, dyke size) and the resident and refilling magma composition and volatile content, in order to study the influence of these factors on the simulated dynamics. Simulation results allow to follow each step of the gas-rich magma ascent in the denser magma, highlighting the details of magma convection and mixing. In particular, the presence of more CO2 in the deep magma results in more ecient and faster dynamics. Through this simulations the variation of the gravimetric field has been determined. Afterward, the space-time distribution of stress resulting from numerical simulations have been used as boundary conditions for the simulations of the displacement field imposed by the magmatic dynamics on rocks. The properties of the simulated domain (rock density, P and S wave velocities) have been based on data from literature on active and passive tomographic experiments, obtained through a collaboration with A. Zollo at the Dept. of Physics of the Federici II Univeristy in Napoli. The elasto-dynamics simulations allow to determine the variations of the space-time distribution of deformation and the seismic signal associated with the studied magmatic dynamics. In particular, results show that these dynamics induce deformations similar to those measured at Campi Flegrei and seismic signals with energies concentrated on the typical frequency bands observed in volcanic areas. The present work shows that an approach based on the solution of equations describing the physics of processes within a magmatic fluid and the surrounding rock system is able to recognise and describe the relationships between geophysical signals detectable on the surface and deep magma dynamics. Therefore, the results suggest that the combined study of geophysical data and informations from numerical simulations can allow in a near future a more ecient evaluation of the short term volcanic hazard.

Comparison of model potentials for molecular dynamics simulation of crystalline silica

Simulationen von SiO2 mit dem von van Beest, Kramer und vanSanten (BKS) entwickelten Paarpotenzial erzeugen vielezufriedenstellende Ergebnisse, aber auch charakteristischeSchwachstellen. In dieser Arbeit wird das BKS-Potenzial mitzwei kürzlich vorgeschlagenen Potenzialen verglichen, dieeffektiv Mehrteilchen-Wechselwirkungen beinhalten. Der ersteAnsatz erlaubt dazu fluktuierende Ladungen, der zweiteinduzierbare Polarisierungen auf den Sauerstoffatomen. Die untersuchten Schwachstellen des BKS Potenzialsbeinhalten das Verhältnis der zwei Gitterkonstanten a und cim Quarzübergang, das von BKS falsch beschrieben wird.Cristobalit und Tridymit erscheinen instabil mit BKS.Weiterhin zeigt die BKS-Zustandsdichte charakteristischeAbweichungen von der wahren Zustandsdichte. DerÜbergangsdruck für den Stishovit I-II Übergang wird deutlichüberschätzt. Das Fluktuierende-Ladungs-Modell verbesserteinige der genannten Punkte, reproduziert aber viele andereEigenschaften schlechter als BKS. DasFluktierende-Dipol-Modell dagegen behebt alle genanntenArtefakte. Zusätzlich wird der druckinduzierte Phasenübergang imalpha-Quarz untersucht. Alle Potentiale finden die selbeStruktur für Quarz II. Bei anschliessender Dekompressionerzeugt BKS eine weitere Phase, während die beiden anderenPotentiale wieder zum alpha-Quarz zurückkehren. Weiterhinwerden zwei Methoden entwickelt, um die piezoelektrischenKonstanten bei konstantem Druck zu bestimmen. Die Ergebnissegeben Hinweise auf eine möglicherweisenicht-elektrostatische Natur der Polarisierungen imFluktuierende-Dipole-Modell. Mit dieser Interpretation scheint das Fluktuierende-DipolPotential alle verfügbaren experimentellen Daten am bestenvon allen drei untersuchten Ansätzen zu reproduzieren.

Electronic energy transfer processes in pi-conjugated polymers

Currently pi-conjugated polymers are considered as technologically interesting materials to be used as functional building elements for the development of the new generation of optoelectronic devices. More specifically during the last few years, poly-p-phenylene materials have attracted considerable attention for their blue photoluminescence properties. This Thesis deals with the optical properties of the most representative blue light poly-p-phenylene emitters such as poly(fluorene), oligo(fluorene), poly(indenofluorene) and ladder-type penta(phenylene) derivatives. In the present work, laser induced photoluminescence spectroscopy is used as a major tool for the study of the interdependence between the dynamics of the probed photoluminescence, the molecular structures of the prepared polymeric films and the presence of chemical defects. Complementary results obtained by two-dimensional wide-angle X-ray diffraction are reported. These findings show that the different optical properties observed are influenced by the intermolecular solid-state interactions that in turn are controlled by the pendant groups of the polymer backbone. A significant feedback is delivered regarding the positive impact of a new synthetic route for the preparation of a poly(indenofluorene) derivative on the spectral purity of the compound. The energy transfer mechanisms that operate in the studied systems are addressed by doping experiments. After the evaluation of the structure/property interdependence, a new optical excitation pathway is presented. An efficient photon low-energy up-conversion that sensitises the blue emission of poly(fluorene) is demonstrated. The observed phenomenon takes place in poly(fluorene) derivatives hosts doped with metallated octaethyl porphyrins, after quasi-CW photoexcitation of intensities in the order of kW/cm2. The up-conversion process is parameterised in terms of temperature, wavelength excitation and central metal cation in the porphyrin ring. Additionally the observation of the up-conversion is extended in a broad range of poly-p-phenylene blue light emitting hosts. The dependence of the detected up-conversion intensity on the excitation intensity and doping concentration is reported. Furthermore the dynamics of the up-conversion intensity are monitored as a function of the doping concentration. These experimental results strongly suggest the existence of triplet-triplet annihilation events into the porphyrin molecules that are subsequently followed by energy transfer to the host. After confirming the occurrence of the up-conversion in solutions, cyclic voltammetry is used in order to show that the up-conversion efficiency is partially determined from the energetic alignment between the HOMO levels of the host and the dopant.

N-Heterocyclic carbene complexes of silver, rhodium and iron: structures, dynamics and catalysis

The research performed in the framework of this Master Thesis has been directly inspired by the recent work of an organometallic research group led by Professor Maria Cristina Cassani on a topic related to the structures, dynamics and catalytic activity of N-heterocyclic carbene-amide rhodium(I) complexes1. A series of [BocNHCH2CH2ImR]X (R = Me, X = I, 1a’; R = Bz, X = Br, 1b’; R = trityl, X = Cl, 1c’) amide-functionalized imidazolium salts bearing increasingly bulky N-alkyl substituents were synthetized and characterized. Subsequently, these organic precursors were employed in the synthesis of silver(I) complexes as intermediate compounds on a way to rhodium(I) complexes [Rh(NBD)X(NHC)] (NHC = 1-(2-NHBoc-ethyl)-3-R-imidazolin-2-ylidene; X = Cl, R = Me (3a’), R = Bz (3b’), R = trityl (3c’); X = I, R = Me (4a’)). VT NMR studies of these complexes revealed a restricted rotation barriers about the metal-carbene bond. However, while the rotation barriers calculated for the complexes in which R = Me, Bz (3a’,b’ and 4a) matched the experimental values, this was not true in the trityl case 3c’, where the experimental value was very similar to that obtained for compound 3b’ and much smaller with respect to the calculated one. In addition, the energy barrier derived for 3c’ from line shape simulation showed a strong dependence on the temperature, while the barriers measured for 3a’,b’ did not show this effect. In view of these results and in order to establish the reasons for the previously found inconsistency between calculated and experimental thermodynamic data, the first objective of this master thesis was the preparation of a series of rhodium(I) complexes [Rh(NBD)X(NHC)] (NHC = 1-benzyl-3-R-imidazolin-2-ylidene; X = Cl, R = Me, Bz, trityl, tBu), containing the benzyl substituent as a chiral probe, followed by full characterization. The second objective of this work was to investigate the catalytic activity of the new rhodium compounds in the hydrosilylation of terminal alkynes for comparison purposes with the reported complexes. Another purpose of this work was to employ the prepared N-heterocyclic ligands in the synthesis of iron(II)-NHC complexes.

The coupling between electron transfer and Protein/Solvent Dynamics in Photosynthetic Reaction Centers: Spectroscopic Studies in Amorphous Matrices

We investigated at the molecular level protein/solvent interactions and their relevance in protein function through the use of amorphous matrices at room temperature. As a model protein, we used the bacterial photosynthetic reaction center (RC) of Rhodobacter sphaeroides, a pigment protein complex which catalyzes the light-induced charge separation initiating the conversion of solar into chemical energy. The thermal fluctuations of the RC and its dielectric conformational relaxation following photoexcitation have been probed by analyzing the recombination kinetics of the primary charge-separated (P+QA-) state, using time resolved optical and EPR spectroscopies. We have shown that the RC dynamics coupled to this electron transfer process can be progressively inhibited at room temperature by decreasing the water content of RC films or of RC-trehalose glassy matrices. Extensive dehydration of the amorphous matrices inhibits RC relaxation and interconversion among conformational substates to an extent comparable to that attained at cryogenic temperatures in water-glycerol samples. An isopiestic method has been developed to finely tune the hydration level of the system. We have combined FTIR spectral analysis of the combination and association bands of residual water with differential light-minus-dark FTIR and high-field EPR spectroscopy to gain information on thermodynamics of water sorption, and on structure/dynamics of the residual water molecules, of protein residues and of RC cofactors. The following main conclusions were reached: (i) the RC dynamics is slaved to that of the hydration shell; (ii) in dehydrated trehalose glasses inhibition of protein dynamics is most likely mediated by residual water molecules simultaneously bound to protein residues and sugar molecules at the protein-matrix interface; (iii) the local environment of cofactors is not involved in the conformational dynamics which stabilizes the P+QA-; (iv) this conformational relaxation appears to be rather delocalized over several aminoacidic residues as well as water molecules weakly hydrogen-bonded to the RC.

Influence of spatial constraints in non-crystalline regions on the polymer dynamics in semi-crystalline polyethylene: a solid state NMR study

A broad variety of solid state NMR techniques were used to investigate the chain dynamics in several polyethylene (PE) samples, including ultrahigh molecular weight PEs (UHMW-PEs) and low molecular weight PEs (LMW-PEs). Via changing the processing history, i.e. melt/solution crystallization and drawing processes, these samples gain different morphologies, leading to different molecular dynamics. Due to the long chain nature, the molecular dynamics of polyethylene can be distinguished in local fluctuation and long range motion. With the help of NMR these different kinds of molecular dynamics can be monitored separately. In this work the local chain dynamics in non-crystalline regions of polyethylene samples was investigated via measuring 1H-13C heteronuclear dipolar coupling and 13C chemical shift anisotropy (CSA). By analyzing the motionally averaged 1H-13C heteronuclear dipolar coupling and 13C CSA, the information about the local anisotropy and geometry of motion was obtained. Taking advantage of the big difference of the 13C T1 relaxation time in crystalline and non-crystalline regions of PEs, the 1D 13C MAS exchange experiment was used to investigate the cooperative chain motion between these regions. The different chain organizations in non-crystalline regions were used to explain the relationship between the local fluctuation and the long range motion of the samples. In a simple manner the cooperative chain motion between crystalline and non-crystalline regions of PE results in the experimentally observed diffusive behavior of PE chain. The morphological influences on the diffusion motion have been discussed. The morphological factors include lamellar thickness, chain organization in non-crystalline regions and chain entanglements. Thermodynamics of the diffusion motion in melt and solution crystallized UHMW-PEs is discussed, revealing entropy-controlled features of the chain diffusion in PE. This thermodynamic consideration explains the counterintuitive relationship between the local fluctuation and the long range motion of the samples. Using the chain diffusion coefficient, the rates of jump motion in crystals of the melt crystallized PE have been calculated. A concept of "effective" jump motion has been proposed to explain the difference between the values derived from the chain diffusion coefficients and those in literatures. The observations of this thesis give a clear demonstration of the strong relationship between the sample morphology and chain dynamics. The sample morphologies governed by the processing history lead to different spatial constraints for the molecular chains, leading to different features of the local and long range chain dynamics. The knowledge of the morphological influence on the microscopic chain motion has many implications in our understanding of the alpha-relaxation process in PE and the related phenomena such as crystal thickening, drawability of PE, the easy creep of PE fiber, etc.

N-Heterocyclic carbene complexes of rhodium: structures, dynamics and catalysis

A series of imidazolium salts of the type [BocNHCH2CH2ImR]X (Boc = t-Bu carbamates; Im = imidazole) (R = Me, X = I, 1a; R = Bn, X = Br, 1b; R = Trityl, X = Cl, 1c) and [BnImR’]X (R’ = Me, X = Br, 1d; R’ = Bn, X = Br, 1e; R’ = Trityl, X = Cl, 1g; R’ = tBu, X = Br, 1h) bearing increasingly bulky substituents were synthetized and characterized. Subsequently, these precursors were employed in the synthesis of silver(I)-N-heterocyclic (NHC) complexes as transmetallating reagents for the preparation of rhodium(I) complexes [RhX(NBD)(NHC)] (NHC = 1-(2-NHBoc-ethyl)-3-R-imidazolin-2-ylidene; X = Cl; R = Me, 4a; R = Bn, 4b; R = Trityl, 4c; X = I, R = Me, 5a; NHC = 1-Bn-3-R’-imidazolin-2-ylidene; X = Cl; R’ = Me, 4d, R’ = Bn, 4e, R’ = Trityl, 4g; R’ = tBu, 4h). VT NMR studies of these complexes revealed a restricted rotation barriers about the metal-carbene bond. While the rotation barriers calculated for the complexes in which R = Me, Bn (4a,b,d,e and 5a) matched the experimental values, this was not true for the complexes 4c,g, bearing a trityl group for which the values are much smaller than the calculated ones. Energy barriers for 4c,g, derived from a line shape simulation, showed a strong dependence on the temperature while for 4h the rotational energy barrier is stopped at room temperature. The catalytic activity of the new rhodium compounds was investigated in the hydrosilylation of terminal alkynes and in the addition of phenylboronic acid to benzaldehyde. The imidazolium salts 1d,e were also employed in the synthesis of new iron(II)-NHC complexes. Finally, during a six-months stay at the University of York a new ligand derived from Norharman was prepared and employed in palladium-mediated cross-coupling.

Magnetization dynamics in spin valves

Key technology applications like magnetoresistive sensors or the Magnetic Random Access Memory (MRAM) require reproducible magnetic switching mechanisms. i.e. predefined remanent states. At the same time advanced magnetic recording schemes push the magnetic switching time into the gyromagnetic regime. According to the Landau-Lifschitz-Gilbert formalism, relevant questions herein are associated with magnetic excitations (eigenmodes) and damping processes in confined magnetic thin film structures.rnObjects of study in this thesis are antiparallel pinned synthetic spin valves as they are extensively used as read heads in today’s magnetic storage devices. In such devices a ferromagnetic layer of high coercivity is stabilized via an exchange bias field by an antiferromagnet. A second hard magnetic layer, separated by a non-magnetic spacer of defined thickness, aligns antiparallel to the first. The orientation of the magnetization vector in the third ferromagnetic NiFe layer of low coercivity - the freelayer - is then sensed by the Giant MagnetoResistance (GMR) effect. This thesis reports results of element specific Time Resolved Photo-Emission Electron Microscopy (TR-PEEM) to image the magnetization dynamics of the free layer alone via X-ray Circular Dichroism (XMCD) at the Ni-L3 X-ray absorption edge.rnThe ferromagnetic systems, i.e. micron-sized spin valve stacks of typically deltaR/R = 15% and Permalloy single layers, were deposited onto the pulse leading centre stripe of coplanar wave guides, built in thin film wafer technology. The ferromagnetic platelets have been applied with varying geometry (rectangles, ellipses and squares), lateral dimension (in the range of several micrometers) and orientation to the magnetic field pulse to study the magnetization behaviour in dependence of these magnitudes. The observation of magnetic switching processes in the gigahertz range became only possible due to the joined effort of producing ultra-short X-ray pulses at the synchrotron source BESSY II (operated in the so-called low-alpha mode) and optimizing the wave guide design of the samples for high frequency electromagnetic excitation (FWHM typically several 100 ps). Space and time resolution of the experiment could be reduced to d = 100 nm and deltat = 15 ps, respectively.rnIn conclusion, it could be shown that the magnetization dynamics of the free layer of a synthetic GMR spin valve stack deviates significantly from a simple phase coherent rotation. In fact, the dynamic response of the free layer is a superposition of an averaged critically damped precessional motion and localized higher order spin wave modes. In a square platelet a standing spin wave with a period of 600 ps (1.7 GHz) was observed. At a first glance, the damping coefficient was found to be independent of the shape of the spin-valve element, thus favouring the model of homogeneous rotation and damping. Only by building the difference in the magnetic rotation between the central region and the outer rim of the platelet, the spin wave becomes visible. As they provide an additional efficient channel for energy dissipation, spin waves contribute to a higher effective damping coefficient (alpha = 0.01). Damping and magnetic switching behaviour in spin valves thus depend on the geometry of the element. Micromagnetic simulations reproduce the observed higher-order spin wave mode.rnBesides the short-run behaviour of the magnetization of spin valves Permalloy single layers with thicknesses ranging from 3 to 40 nm have been studied. The phase velocity of a spin wave in a 3 nm thick ellipse could be determined to 8.100 m/s. In a rectangular structure exhibiting a Landau-Lifschitz like domain pattern, the speed of the field pulse induced displacement of a 90°-Néel wall has been determined to 15.000 m/s.rn

Properties of cosmic voids in dark energy simulations

The last decade has witnessed the establishment of a Standard Cosmological Model, which is based on two fundamental assumptions: the first one is the existence of a new non relativistic kind of particles, i. e. the Dark Matter (DM) that provides the potential wells in which structures create, while the second one is presence of the Dark Energy (DE), the simplest form of which is represented by the Cosmological Constant Λ, that sources the acceleration in the expansion of our Universe. These two features are summarized by the acronym ΛCDM, which is an abbreviation used to refer to the present Standard Cosmological Model. Although the Standard Cosmological Model shows a remarkably successful agreement with most of the available observations, it presents some longstanding unsolved problems. A possible way to solve these problems is represented by the introduction of a dynamical Dark Energy, in the form of the scalar field ϕ. In the coupled DE models, the scalar field ϕ features a direct interaction with matter in different regimes. Cosmic voids are large under-dense regions in the Universe devoided of matter. Being nearby empty of matter their dynamics is supposed to be dominated by DE, to the nature of which the properties of cosmic voids should be very sensitive. This thesis work is devoted to the statistical and geometrical analysis of cosmic voids in large N-body simulations of structure formation in the context of alternative competing cosmological models. In particular we used the ZOBOV code (see ref. Neyrinck 2008), a publicly available void finder algorithm, to identify voids in the Halos catalogues extraxted from CoDECS simulations (see ref. Baldi 2012 ). The CoDECS are the largest N-body simulations to date of interacting Dark Energy (DE) models. We identify suitable criteria to produce voids catalogues with the aim of comparing the properties of these objects in interacting DE scenarios to the standard ΛCDM model, at different redshifts. This thesis work is organized as follows: in chapter 1, the Standard Cosmological Model as well as the main properties of cosmic voids are intro- duced. In chapter 2, we will present the scalar field scenario. In chapter 3 the tools, the methods and the criteria by which a voids catalogue is created are described while in chapter 4 we discuss the statistical properties of cosmic voids included in our catalogues. In chapter 5 the geometrical properties of the catalogued cosmic voids are presented by means of their stacked profiles. In chapter 6 we summarized our results and we propose further developments of this work.

Evaluation of energy level to be absorbed by tractor ROPS: Actual Tests, Simulation and Computation

A highly dangerous situations for tractor driver is the lateral rollover in operating conditions. Several accidents, involving tractor rollover, have indeed been encountered, requiring the design of a robust Roll-Over Protective Structure (ROPS). The aim of the thesis was to evaluate tractor behaviour in the rollover phase so as to calculate the energy absorbed by the ROPS to ensure driver safety. A Mathematical Model representing the behaviour of a generic tractor during a lateral rollover, with the possibility of modifying the geometry, the inertia of the tractor and the environmental boundary conditions, is proposed. The purpose is to define a method allowing the prediction of the elasto-plastic behaviour of the subsequent impacts occurring in the rollover phase. A tyre impact model capable of analysing the influence of the wheels on the energy to be absorbed by the ROPS has been also developed. Different tractor design parameters affecting the rollover behaviour, such as mass and dimensions, have been considered. This permitted the evaluation of their influence on the amount of energy to be absorbed by the ROPS. The mathematical model was designed and calibrated with respect to the results of actual lateral upset tests carried out on a narrow-track tractor. The dynamic behaviour of the tractor and the energy absorbed by the ROPS, obtained from the actual tests, showed to match the results of the model developed. The proposed approach represents a valuable tool in understanding the dynamics (kinetic energy) and kinematics (position, velocity, angular velocity, etc.) of the tractor in the phases of lateral rollover and the factors mainly affecting the event. The prediction of the amount of energy to be absorbed in some cases of accident is possible with good accuracy. It can then help in designing protective structures or active security devices.

Structure and dynamics of ionic liquid,electrode interfaces

In the early 20th century, Gouy, Chapman, and Stern developed a theory to describe the capacitance and the spatial ion distribution of diluted electrolytes near an electrode. After a century of research, considerable progress has been made in the understanding of the electrolyte/electrode interface. However, its molecular-scale structure and its variation with an applied potential is still under debate. In particular for room-temperature ionic liquids, a new class of solventless electrolytes, the classical theories for the electrical double layer are not applicable. Recently, molecular dynamics simulations and phenomenological theories have attempted to explain the capacitance of the ionic liquid/electrode interface with the molecular-scale structure and dynamics of the ionic liquid near the electrode. rnHowever, experimental evidence is very limited. rnrnIn the presented study, the ion distribution of an ionic liquid near an electrode and its response to applied potentials was examined with sub-molecular resolution. For this purpose, a new sample chamber was constructed, allowing in situ high energy X-ray reflectivity experiments under potential control, as well as impedance spectroscopy measurements. The combination of structural information and electrochmical data provided a comprehensive picture of the electric double layer in ionic liquids. Oscillatory charge density profiles were found, consisting of alternating anion- and cation-enriched layers at both, cathodic and anodic, potentials. This structure was shown to arise from the same ion-ion correlations dominating the liquid bulk structure that were observed as a distinct X-ray diffraction peak. Therefore, existing physically motivated models were refined and verified by comparison with independent measurements. rnrnThe relaxation dynamics of the interfacial structure upon potential variation were studied by time resolved X-ray reflectivity experiments with sub-millisecond resolution. The observed relaxation times during charging/discharging are consistent with the impedance spectroscopy data revealing three processes of vastly different characteristic time-scales. Initially, the ion transport normal to the interface happens on a millisecond-scale. Another 100-millisecond-scale process is associated with molecular reorientation of electrode-adsorbed cations. Further, a minute-scale relaxation was observed, which is tentatively assigned to lateral ordering within the first layer.