Vibrational dynamics and structure of natural feldspars

The present study gives a contribution to the knowledge on the Na-feldspar and plagioclases, extending the database of the Raman spectra of plagioclases with different chemical compositions and structural orders. This information may be used for the future planetary explorations by “rovers”, for the investigation of ceramics nanocrystal materials and for the mineralogical phase identification in sediments. Na-feldspar and plagioclase solid solution have been investigated by Raman spectroscopy in order to determine the relationships between the vibrational changes and the plagioclase crystal chemistry and structure. We focused on the Raman micro-spectroscopy technique, being a non-destructive method, suited for contactless analysis with high spatial resolution. Chemical and structural analyses have been performed on natural samples to test the usefulness of Raman spectroscopy as a tool in the study of the pressure-induced structural deformations, the disordering processes due to change in the Al-Si distribution in the tetrahedral sites and, finally, in the determination of the anorthitic content (Anx) in plagioclase minerals. All the predicted 39 Ag Raman active modes have been identified and assigned to specific patterns of atomic vibrational motion. A detailed comparison between experimental and computed Raman spectra has been performed and previous assignments have been revised, solving some discrepancies reported in recent literature. The ab initio calculation at the hybrid HF/DFT level with the WC1LYP Hamiltonian has proven to give excellent agreement between calculated and experimentally measured Raman wavenumbers and intensities in triclinic minerals. A short digression on the 36 infrared active modes of Na-feldspar has been done too. The identification of all 39 computed Raman modes in the experimentally measured spectra of the fully ordered Na-feldspar, known as low albite, along with the detailed description of each vibrational mode, has been essential to extend the comparative analysis to the high pressure and high temperature structural forms of albite, which reflect the physical–chemical conditions of the hosting rocks. The understanding of feldspar structure response to pressure and temperature is crucial in order to constrain crustal behaviour. The compressional behaviour of the Na-feldspar has been investigated for the first time by Raman spectroscopy. The absence of phase transitions and the occurrence of two secondary compression mechanisms acting at different pressures have been confirmed. Moreover, Raman data suggest that the internal structural changes are confined to a small pressure interval, localized around 6 GPa, not spread out from 4 to 8 GPa as suggested by previous X-rays studies on elasticity. The dominant compression mechanisms act via tetrahedral tilting, while the T-O bond lengths remain nearly constant at moderate compressional regimes. At the spectroscopic level, this leads to the strong pressure dependencies of T-O-T bending modes, as found for the four modes at 478, 508, 578 and 815 cm-1. The Al-Si distribution in the tetrahedral sites affects also the Raman spectrum of Na-feldspar. In particular, peak broadening is more sensitive than peak position to changes in the degree of order. Raman spectroscopy is found to be a good probe for local ordering, in particular being sensitive to the first annealing steps, when the macroscopic order parameter is still high. Even though Raman data are scattered and there are outliers in the estimated values of the degree of order, the average peak linewidths of the Na-feldspar characteristic doublet band, labelled here as υa and υb, as a function of the order parameter Qod show interesting trends: both peak linewidths linearly increase until saturation. From Qod values lower than 0.6, peak broadening is no more affected by the Al-Si distribution. Moreover, the disordering process is found to be heterogeneous. SC-XRD and Raman data have suggested an inter-crystalline inhomogeneity of the samples, i.e., the presence of regions with different defect density on the micrometric scale. Finally, the influence of Ca-Na substitution in the plagioclase Raman spectra has been investigated. Raman spectra have been collected on a series of well characterized natural, low structural plagioclases. The variations of the Raman modes as a function of the chemical composition and the structural order have been determined. The number of the observed Raman bands at each composition gives information about the unit-cell symmetry: moving away from the C1 structures, the number of the Raman bands enhances, as the number of formula units in the unit cell increases. The modification from an “albite-like” Raman spectrum to a more “anorthite-like” spectrum occurs from sample An78 onwards, which coincides with the appearance of c reflections in the diffraction patterns of the samples. The evolution of the Raman bands υa and υb displays two changes in slope at ~An45 and ~An75: the first one occurs between e2 and e1 plagioclases, the latter separates e1 and I1 plagioclases with only b reflections in their diffraction patterns from I1 and P1 samples having b and c reflections too. The first variation represents exactly the e2→e1 phase transitions, whereas the second one corresponds in good approximation to the C1→I1 transition, which has been determined at ~An70 by previous works. The I1→P1 phase transition in the anorthite-rich side of the solid solution is not highlighted in the collected Raman spectra. Variations in peak broadening provide insights into the behaviour of the order parameter on a local scale, suggesting an increase in the structural disorder within the solid solution, as the structures have to incorporate more Al atoms to balance the change from monovalent to divalent cations. All the information acquired on these natural plagioclases has been used to produce a protocol able to give a preliminary estimation of the chemical composition of an unknown plagioclase from its Raman spectrum. Two calibration curves, one for albite-rich plagioclases and the other one for the anorthite-rich plagioclases, have been proposed by relating the peak linewidth of the most intense Raman band υa and the An content. It has been pointed out that the dependence of the composition from the linewidth can be obtained only for low structural plagioclases with a degree of order not far away from the references. The proposed tool has been tested on three mineralogical samples, two of meteoric origin and one of volcanic origin. Chemical compositions by Raman spectroscopy compare well, within an error of about 10%, with those obtained by elemental techniques. Further analyses on plagioclases with unknown composition will be necessary to validate the suggested method and introduce it as routine tool for the determination of the chemical composition from Raman data in planetary missions.

Noisy random Boolean formulae:a statistical physics perspective

Properties of computing Boolean circuits composed of noisy logical gates are studied using the statistical physics methodology. A formula-growth model that gives rise to random Boolean functions is mapped onto a spin system, which facilitates the study of their typical behavior in the presence of noise. Bounds on their performance, derived in the information theory literature for specific gates, are straightforwardly retrieved, generalized and identified as the corresponding macroscopic phase transitions. The framework is employed for deriving results on error-rates at various function-depths and function sensitivity, and their dependence on the gate-type and noise model used. These are difficult to obtain via the traditional methods used in this field.

On the robustness of random Boolean formulae

Random Boolean formulae, generated by a growth process of noisy logical gates are analyzed using the generating functional methodology of statistical physics. We study the type of functions generated for different input distributions, their robustness for a given level of gate error and its dependence on the formulae depth and complexity and the gates used. Bounds on their performance, derived in the information theory literature for specific gates, are straightforwardly retrieved, generalized and identified as the corresponding typical-case phase transitions. Results for error-rates, function-depth and sensitivity of the generated functions are obtained for various gate-type and noise models. © 2010 IOP Publishing Ltd.

Neural population coding is optimized by discrete tuning curves

The sigmoidal tuning curve that maximizes the mutual information for a Poisson neuron, or population of Poisson neurons, is obtained. The optimal tuning curve is found to have a discrete structure that results in a quantization of the input signal. The number of quantization levels undergoes a hierarchy of phase transitions as the length of the coding window is varied. We postulate, using the mammalian auditory system as an example, that the presence of a subpopulation structure within a neural population is consistent with an optimal neural code.

Dielectric and electro-optical properties of some cyanobiphenyl liquid-crystals

A variety of methods have been reviewed for obtaining parallel or perpendicular alignment in liquid-crystal cells. Some of these methods have been selected and developed and were used in polarised spectroscopy, dielectric and electro-optic studies. Also, novel dielectric and electro-optic cells were constructed for use over a range of temperature. Dielectric response of thin layers of E7 and E8 (eutectic mixture liquid-crystals) have been measured in the frequency range (12 Hz-100 kHz) and over a range of temperature (183-337K). Dielectric spectra were also obtained for supercooled E7 and E8 in the Hz and kHz range. When the measuring electric field was parallel to the nematic director, one loss peak (low-frequency relaxation process) was observed for E7 and for E8, that exhibits a Debye-type behaviour in the supercooled systems. When the measuring electric field was perpendicular to the nematic director, two resolved dielectric processes have been observed. The phase transitions, effective molecular polarisabilities, anisotropy of polarisabilities and order parameters of three liquid crystal homologs (5CB, 6CB, and 7CB), 60CB and three eutectic nematic mixtures E7, E8, and E607 were calculated using optical and density data measured at several temperatures. The order parameters calculated using the different methods of Vuks, Neugebauer, Saupe-Maier, and Palffy-Muhoray are nearly the same for the liquid crystals considered in the present study. Also, the interrelationship between density and refractive index and the molecular structure of these liquid crystals were established. Accurate dielectric and dipole results of a range of liquid-crystal forming molecules at several temperatures have reported. The role of the cyano-end group, biphenyl core, and flexible tail in molecular association, were investigated using the dielectric method for some molecules which have a structural relationship to the nematogens. Analysis of the dielectric data for solution of the liquid-crystals indicated a high molecular association, comparable to that observed in the nematic or isotropic phases. Electro-optic Kerr effect were investigated for some alkyl cyanobiphenyls, their nematic mixtures and the eutectic mixture liquid-crystals E7 and E8 in the isotropic phase and solution. The Kerr constant of these liquid crystals found to be very high at the nematic-isotropic transition temperatures as the molecules are expected to be highly ordered close to phase transition temperatures. Dynamic Kerr effect behaviour and transient molecular reorientation were also observed in thin layers of some alkyl cyanobiphenyls. Dichroic ratio R and order parameters of solutions containing some azo and anthraquinone dyes in the nematic solvent (E7 and E8), were investigated by the measurement of the intensity of the absorption bands in the visible region of parallel aligned samples. The effective factors on the dichroic ratio of the dyes dissolved in the nematic solvents were determined and discussed.

Statistical physics of error-correcting codes

In this thesis we use statistical physics techniques to study the typical performance of four families of error-correcting codes based on very sparse linear transformations: Sourlas codes, Gallager codes, MacKay-Neal codes and Kanter-Saad codes. We map the decoding problem onto an Ising spin system with many-spins interactions. We then employ the replica method to calculate averages over the quenched disorder represented by the code constructions, the arbitrary messages and the random noise vectors. We find, as the noise level increases, a phase transition between successful decoding and failure phases. This phase transition coincides with upper bounds derived in the information theory literature in most of the cases. We connect the practical decoding algorithm known as probability propagation with the task of finding local minima of the related Bethe free-energy. We show that the practical decoding thresholds correspond to noise levels where suboptimal minima of the free-energy emerge. Simulations of practical decoding scenarios using probability propagation agree with theoretical predictions of the replica symmetric theory. The typical performance predicted by the thermodynamic phase transitions is shown to be attainable in computation times that grow exponentially with the system size. We use the insights obtained to design a method to calculate the performance and optimise parameters of the high performance codes proposed by Kanter and Saad.

From the physics of interacting polymers to optimizing routes on the London Underground

Optimizing paths on networks is crucial for many applications, ranging from subway traffic to Internet communication. Because global path optimization that takes account of all path choices simultaneously is computationally hard, most existing routing algorithms optimize paths individually, thus providing suboptimal solutions. We use the physics of interacting polymers and disordered systems to analyze macroscopic properties of generic path optimization problems and derive a simple, principled, generic, and distributed routing algorithm capable of considering all individual path choices simultaneously. We demonstrate the efficacy of the algorithm by applying it to: (i) random graphs resembling Internet overlay networks, (ii) travel on the London Underground network based on Oyster card data, and (iii ) the global airport network. Analytically derived macroscopic properties give rise to insightful new routing phenomena, including phase transitions and scaling laws, that facilitate better understanding of the appropriate operational regimes and their limitations, which are difficult to obtain otherwise.

Structural Characterization of Metal Hydrides for Energy Applications

Hydrogen can be an unlimited source of clean energy for future because of its very high energy density compared to the conventional fuels like gasoline. An efficient and safer way of storing hydrogen is in metals and alloys as hydrides. Light metal hydrides, alanates and borohydrides have very good hydrogen storage capacity, but high operation temperatures hinder their application. Improvement of thermodynamic properties of these hydrides is important for their commercial use as a source of energy. Application of pressure on materials can have influence on their properties favoring hydrogen storage. Hydrogen desorption in many complex hydrides occurs above the transition temperature. Therefore, it is important to study the physical properties of the hydride compounds at ambient and high pressure and/or high temperature conditions, which can assist in the design of suitable storage materials with desired thermodynamic properties. ^ The high pressure-temperature phase diagram, thermal expansion and compressibility have only been evaluated for a limited number of hydrides so far. This situation serves as a main motivation for studying such properties of a number of technologically important hydrides. Focus of this dissertation was on X-ray diffraction and Raman spectroscopy studies of Mg2FeH6, Ca(BH4) 2, Mg(BH4)2, NaBH4, NaAlH4, LiAlH4, LiNH2BH3 and mixture of MgH 2 with AlH3 or Si, at different conditions of pressure and temperature, to obtain their bulk modulus and thermal expansion coefficient. These data are potential source of information regarding inter-atomic forces and also serve as a basis for developing theoretical models. Some high pressure phases were identified for the complex hydrides in this study which may have better hydrogen storage properties than the ambient phase. The results showed that the highly compressible B-H or Al-H bonds and the associated bond disordering under pressure is responsible for phase transitions observed in brorohydrides or alanates. Complex hydrides exhibited very high compressibility suggesting possibility to destabilize them with pressure. With high capacity and favorable thermodynamics, complex hydrides are suitable for reversible storage. Further studies are required to overcome the kinetic barriers in complex hydrides by catalytic addition. A comparative study of the hydride properties with that of the constituting metal, and their inter relationships were carried out with many interesting features.^

Exploring thermal and mechanical properties of selected transition elements under extreme conditions: Experiments at high pressures and high temperatures

Transition metals (Ti, Zr, Hf, Mo, W, V, Nb, Ta, Pd, Pt, Cu, Ag, and Au) are essential building units of many materials and have important industrial applications. Therefore, it is important to understand their thermal and physical behavior when they are subjected to extreme conditions of pressure and temperature. This dissertation presents: • An improved experimental technique to use lasers for the measurement of thermal conductivity of materials under conditions of very high pressure (P, up to 50 GPa) and temperature (T up to 2500 K). • An experimental study of the phase relationship and physical properties of selected transition metals, which revealed new and unexpected physical effects of thermal conductivity in Zr, and Hf under high P-T. • New phase diagrams created for Hf, Ti and Zr from experimental data. • P-T dependence of the lattice parameters in α-hafnium. Contrary to prior reports, the α-ω phase transition in hafnium has a negative dT/dP slope. • New data on thermodynamic and physical properties of several transition metals and their respective high P-T phase diagrams. • First complete thermodynamic database for solid phases of 13 common transition metals was created. This database has: All the thermochemical data on these elements in their standard state (mostly available and compiled); All the equations of state (EoS) formulated from pressure-volume-temperature data (measured as a part of this study and from literature); Complete thermodynamic data for selected elements from standard to extreme conditions. The thermodynamic database provided by this study can be used with available thermodynamic software to calculate all thermophysical properties and phase diagrams at high P-T conditions. For readers who do not have access to this software, tabulated values of all thermodynamic and volume data for the 13 metals at high P-T are included in the APPENDIX. In the APPENDIX, a description of several other high-pressure studies of selected oxide systems is also included. Thermophysical properties (Cp, H, S, G) of the high P-T ω-phase of Ti, Zr and Hf were determined during the optimization of the EoS parameters and are presented in this study for the first time. These results should have important implications in understanding hexagonal-close-packed to simple-hexagonal phase transitions in transition metals and other materials.

Numerical study of the enlarged O(5) symmetry of the 3D antiferromagnetic RP^(2) spin model

We investigate by means of Monte Carlo simulation and finite-size scaling analysis the critical properties of the three dimensional O (5) non-linear σ model and of the antiferromagnetic RP^(2) model, both of them regularized on a lattice. High accuracy estimates are obtained for the critical exponents, universal dimensionless quantities and critical couplings. It is concluded that both models belong to the same universality class, provided that rather non-standard identifications are made for the momentum-space propagator of the RP^(2) model. We have also investigated the phase diagram of the RP^(2) model extended by a second-neighbor interaction. A rich phase diagram is found, where most of the phase transitions are of the first order.

Physics of Hexagonal Limit-Periodic Phases: Thermodynamics, Formation and Vibrational Modes

Limit-periodic (LP) structures exhibit a type of nonperiodic order yet to be found in a natural material. A recent result in tiling theory, however, has shown that LP order can spontaneously emerge in a two-dimensional (2D) lattice model with nearest-and next-nearest-neighbor interactions. In this dissertation, we explore the question of what types of interactions can lead to a LP state and address the issue of whether the formation of a LP structure in experiments is possible. We study emergence of LP order in three-dimensional (3D) tiling models and bring the subject into the physical realm by investigating systems with realistic Hamiltonians and low energy LP states. Finally, we present studies of the vibrational modes of a simple LP ball and spring model whose results indicate that LP materials would exhibit novel physical properties.

A 2D lattice model defined on a triangular lattice with nearest- and next-nearest-neighbor interactions based on the Taylor-Socolar (TS) monotile is known to have a LP ground state. The system reaches that state during a slow quench through an infinite sequence of phase transitions. Surprisingly, even when the strength of the next-nearest-neighbor interactions is zero, in which case there is a large degenerate class of both crystalline and LP ground states, a slow quench yields the LP state. The first study in this dissertation introduces 3D models closely related to the 2D models that exhibit LP phases. The particular 3D models were designed such that next-nearest-neighbor interactions of the TS type are implemented using only nearest-neighbor interactions. For one of the 3D models, we show that the phase transitions are first order, with equilibrium structures that can be more complex than in the 2D case.

In the second study, we investigate systems with physical Hamiltonians based on one of the 2D tiling models with the goal of stimulating attempts to create a LP structure in experiments. We explore physically realizable particle designs while being mindful of particular features that may make the assembly of a LP structure in an experimental system difficult. Through Monte Carlo (MC) simulations, we have found that one particle design in particular is a promising template for a physical particle; a 2D system of identical disks with embedded dipoles is observed to undergo the series of phase transitions which leads to the LP state.

LP structures are well ordered but nonperiodic, and hence have nontrivial vibrational modes. In the third section of this dissertation, we study a ball and spring model with a LP pattern of spring stiffnesses and identify a set of extended modes with arbitrarily low participation ratios, a situation that appears to be unique to LP systems. The balls that oscillate with large amplitude in these modes live on periodic nets with arbitrarily large lattice constants. By studying periodic approximants to the LP structure, we present numerical evidence for the existence of such modes, and we give a heuristic explanation of their structure.

The origin of incipient ferroelectricity in lead telluride.

The interactions between electrons and lattice vibrations are fundamental to materials behaviour. In the case of group IV-VI, V and related materials, these interactions are strong, and the materials exist near electronic and structural phase transitions. The prototypical example is PbTe whose incipient ferroelectric behaviour has been recently associated with large phonon anharmonicity and thermoelectricity. Here we show that it is primarily electron-phonon coupling involving electron states near the band edges that leads to the ferroelectric instability in PbTe. Using a combination of nonequilibrium lattice dynamics measurements and first principles calculations, we find that photoexcitation reduces the Peierls-like electronic instability and reinforces the paraelectric state. This weakens the long-range forces along the cubic direction tied to resonant bonding and low lattice thermal conductivity. Our results demonstrate how free-electron-laser-based ultrafast X-ray scattering can be utilized to shed light on the microscopic mechanisms that determine materials properties.

Nonequilibrium critical scaling from quantum thermodynamics

The emerging field of quantum thermodynamics is contributing important results and insights into archetypal many-body problems, including quantum phase transitions. Still, the question whether out-of-equilibrium quantities, such as fluctuations of work, exhibit critical scaling after a sudden quench in a closed system has remained elusive. Here, we take a novel approach to the problem by studying a quench across an impurity quantum critical point. By performing density matrix renormalization group computations on the two-impurity Kondo model, we are able to establish that the irreversible work produced in a quench exhibits finite-size scaling at quantum criticality. This scaling faithfully predicts the equilibrium critical exponents for the crossover length and the order parameter of the model, and, moreover, implies a new exponent for the rescaled irreversible work. By connecting the irreversible work to the two-impurity spin correlation function, our findings can be tested experimentally.

Neutron Scattering Studies on Complex Magnetic Structures in Magnetoelectric Materials

In this thesis, the magnetic properties of four transition-metal oxides are presented. Their multiferroic and magnetoelectric phases have been investigated by means of different neutron scattering techniques. The materials TbMnO3 and MnWO4 belong to the group of spin-induced multiferroics. Their ferroelectric polarization can be explained by the inverse DzyaloshinskiiMoriya interaction. Another common feature of both materials is the presence of subsequent magnetic transitions from a spin-density wave to a spin spiral. The features of the phase transitions have been studied in both materials and it could be shown that diffuse magnetic scattering from the spin spiral is present even in the ordered spin-density wave phase. The excitation spectrum in the multiferroic phase of TbMnO3 was investigated in detail and a comprehensive dataset was obtained using time-of-flight spectroscopy. A spin-wave model could be obtained which can quantitatively describe the full dispersion. Furthermore, the polarization of the zone-center excitations could be derived which fit well to data from inelastic neutron spectroscopy and infrared spectroscopy. With the combination of spherical neutron polarimetry and a poling of the sample by an electric field, it was possible to observe the chiral magnetic component of the magnetic excitations in TbMnO3 and MnWO4. The spin-wave model for TbMnO3 obtained in this thesis is able to correctly describe the dispersion of this component. The double tungstate NaFe(WO4)2 is isostructural to the multiferroic MnWO4 and develops a complex magnetic phase diagram. By the use of neutron diffraction techniques, the zero-field structure and high-field structures in magnetic field applied along the b-axis could be determined. The data reveal a direct transition into an incommensurate spin-spiral structure. The value of the incommensurability is driven by anharmonic modulations and shows strong hysteresis effects. The static and dynamic properties in the magnetoelectric spin-glass phase of Ni0.42Mn0.58TiO3 were studied in detail. The spin-glass phase is composed of short-ranged MnTiO3 and NiTiO3-type order. The antiferromagnetic domains could be controlled by crossed magnetic and electric fields, which was visualized using spherical neutron polarimetry. A comprehensive dataset of the magnetic excitations in the spin-glass phase was collected. The dataset revealed correlations in the hexagonal plane which are only weakly coupled along the c-axis. The excitation spectra could be simulated by taking into account the MnTiO3-type order.

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.