Frustration and anomalous behavior in the Bell-Lavis model of liquid water

We have reconsidered the Bell-Lavis model of liquid water and investigated its relation to its isotropic version, the antiferromagnetic Blume-Emery-Griffiths model on the triangular lattice. Our study was carried out by means of an exact solution on the sequential Husimi cactus. We show that the ground states of both models share the same topology and that fluid phases (gas and low- and high-density liquids) can be mapped onto magnetic phases (paramagnetic, antiferromagnetic, and dense paramagnetic, respectively). Both models present liquid-liquid coexistence and several thermodynamic anomalies. This result suggests that anisotropy introduced through orientational variables play no specific role in producing the density anomaly, in agreement with a similar conclusion discussed previously following results for continuous soft core,models. We propose that the presence of liquid anomalies may be related to energetic frustration, a feature common to both models.

Effect of long cyclic exchanges on the magnetic properties of bcc (3)He

Using path-integral Monte Carlo calculations, we have calculated ring exchange frequencies in the bcc phase of solid (3)He for densities from melting to the highest stable density. We evaluate 42 different exchange frequencies from two atoms up to eight atoms and find their Gruneisen exponents. Using a fit to these frequencies, we calculate the contribution to the Curie-Weiss temperature, Theta(CW), and upper critical magnetic field, B(c2), for even longer exchanges using a lattice Monte Carlo procedure. We find that contributions from seven-and eight-particle exchanges make a significant contribution to Theta(CW) and B(c2) at melting density. Comparison with experimental data is given.

A 2D numerical model for simulating the physics of fault systems

Simulations provide a powerful means to help gain the understanding of crustal fault system physics required to progress towards the goal of earthquake forecasting. Cellular Automata are efficient enough to probe system dynamics but their simplifications render interpretations questionable. In contrast, sophisticated elasto-dynamic models yield more convincing results but are too computationally demanding to explore phase space. To help bridge this gap, we develop a simple 2D elastodynamic model of parallel fault systems. The model is discretised onto a triangular lattice and faults are specified as split nodes along horizontal rows in the lattice. A simple numerical approach is presented for calculating the forces at medium and split nodes such that general nonlinear frictional constitutive relations can be modeled along faults. Single and multi-fault simulation examples are presented using a nonlinear frictional relation that is slip and slip-rate dependent in order to illustrate the model.

Fractal-to-Euclidean crossover of the isotropy restoration feature in a family of fractal resistor networks

Fractal with microscopic anisotropy shows a unique type of macroscopic isotropy restoration phenomenon that is absent in Euclidean space [M. T. Barlow et al., Phys. Rev. Lett. 75, 3042]. In this paper the isotropy restoration feature is considered for a family of two-dimensional Sierpinski gasket type fractal resistor networks. A parameter xi is introduced to describe this phenomenon. Our numerical results show that xi satisfies the scaling law xi similar to l(-alpha), where l is the system size and alpha is an exponent independent of the degree of microscopic anisotropy, characterizing the isotropy restoration feature of the fractal systems. By changing the underlying fractal structure towards the Euclidean triangular lattice through increasing the side length b of the gasket generators, the fractal-to-Euclidean crossover behavior of the isotropy restoration feature is discussed.

Compressible Sherrington-Kirkpatrick spin-glass model

We introduce a Sherrington-Kirkpatrick spin-glass model with the addition of elastic degrees of freedom. The problem is formulated in terms of an effective four-spin Hamiltonian in the pressure ensemble, which can be treated by the replica method. In the replica-symmetric approximation, we analyze the pressure-temperature phase diagram, and obtain expressions for the critical boundaries between the disordered and the ordered (spin-glass and ferromagnetic) phases. The second-order para-ferromagnetic border ends at a tricritical point, beyond which the transition becomes discontinuous. We use these results to make contact with the temperature-concentration phase diagrams of mixtures of hydrogen-bonded crystals.

DFT study of the electronic, vibrational, and optical properties of SnO(2)

We report results on the electronic, vibrational, and optical properties of SnO(2) obtained using first-principles calculations performed within the density functional theory. All the calculated phonon frequencies, real and imaginary parts of complex dielectric function, the energy-loss spectrum, the refractive index, the extinction, and the absorption coefficients show good agreement with experimental results. Based on our calculations, the SnO(2) electron and hole effective masses were found to be strongly anisotropic. The lattice contribution to the low-frequency region of the SnO(2) dielectric function arising from optical phonons was also determined resulting the values of E > (1aSyen) (latt) (0) = 14.6 and E > (1ayen) (latt) (0) = 10.7 for directions perpendicular and parallel to the tetragonal c-axis, respectively. This is in excellent agreement with the available experimental data. After adding the electronic contribution to the lattice contribution, a total average value of E >(1)(0) = 18.2 is predicted for the static permittivity constant of SnO(2).

Propagação de danos no modelo de ising em redes de bravais e em fractais

In this work we have studied, by Monte Carlo computer simulation, several properties that characterize the damage spreading in the Ising model, defined in Bravais lattices (the square and the triangular lattices) and in the Sierpinski Gasket. First, we investigated the antiferromagnetic model in the triangular lattice with uniform magnetic field, by Glauber dynamics; The chaotic-frozen critical frontier that we obtained coincides , within error bars, with the paramegnetic-ferromagnetic frontier of the static transition. Using heat-bath dynamics, we have studied the ferromagnetic model in the Sierpinski Gasket: We have shown that there are two times that characterize the relaxation of the damage: One of them satisfy the generalized scaling theory proposed by Henley (critical exponent z~A/T for low temperatures). On the other hand, the other time does not obey any of the known scaling theories. Finally, we have used methods of time series analysis to study in Glauber dynamics, the damage in the ferromagnetic Ising model on a square lattice. We have obtained a Hurst exponent with value 0.5 in high temperatures and that grows to 1, close to the temperature TD, that separates the chaotic and the frozen phases

Existence and stability of new nanoreactors: Highly swollen hexagonal liquid crystals

We report the preparation of direct hexagonal liquid crystals, constituted of oil-swollen cylinders arranged on a triangular lattice in water. The volume ratio of oil over water, rho, can be as large as 3.8. From the lattice parameter measured by small-angle X-ray scattering, we show that all the oil is indeed incorporated into the cylinders, thus allowing the diameter of the cylinders to be controlled over one decade range, provided that the ionic strength of the aqueous medium and rho are varied concomitantly. These hexagonal swollen liquid crystals (SLCs) have been first reported with sodium dodecyl sulfate as anionic surfactant, cyclohexane as solvent, 1-pentanol as co-surfactant, and sodium chloride as salt (Ramos, L.; Fabre, P. Langmuir 1997, 13, 13). The stability of these liquid crystals is investigated when the pH of the aqueous medium or the chemical nature of the components (salt and surfactant) is changed. We demonstrate that the range of stability is quite extended, rendering swollen hexagonal phases potentially useful for the fabrication of nanomaterials. As illustrations, we finally show that gelation of inorganic particles in the continuous aqueous medium of a SLC and polymerization within the oil-swollen cylinders of a SLC can be conducted without disrupting the hexagonal order of the system.

Temperature dependence of structural and electronic properties of the spin-liquid candidate κ-(BEDT-TTF) 2Cu 2(CN) 3

We investigate the effect that the temperature dependence of the crystal structure of a two-dimensional organic charge-transfer salt has on the low-energy Hamiltonian representation of the electronic structure. For that, we determine the crystal structure of κ-(BEDT-TTF) 2Cu 2(CN) 3 for a series of temperatures between T=5 and 300 K by single crystal X-ray diffraction and analyze the evolution of the electronic structure with temperature by using density functional theory and tight binding methods. We find a considerable temperature dependence of the corresponding triangular lattice Hubbard Hamiltonian parameters. We conclude that even in the absence of a change of symmetry, the temperature dependence of quantities like frustration and interaction strength can be significant and should be taken into account. © 2012 American Physical Society.

Computer simulations of two-dimensional colloidal crystals in confinement

Monte Carlo simulations are used to study the effect of confinement on a crystal of point particles interacting with an inverse power law potential in d=2 dimensions. This system can describe colloidal particles at the air-water interface, a model system for experimental study of two-dimensional melting. It is shown that the state of the system (a strip of width D) depends very sensitively on the precise boundary conditions at the two ``walls'' providing the confinement. If one uses a corrugated boundary commensurate with the order of the bulk triangular crystalline structure, both orientational order and positional order is enhanced, and such surface-induced order persists near the boundaries also at temperatures where the system in the bulk is in its fluid state. However, using smooth repulsive boundaries as walls providing the confinement, only the orientational order is enhanced, but positional (quasi-) long range order is destroyed: The mean-square displacement of two particles n lattice parameters apart in the y-direction along the walls then crosses over from the logarithmic increase (characteristic for $d=2$) to a linear increase (characteristic for d=1). The strip then exhibits a vanishing shear modulus. These results are interpreted in terms of a phenomenological harmonic theory. Also the effect of incommensurability of the strip width D with the triangular lattice structure is discussed, and a comparison with surface effects on phase transitions in simple Ising- and XY-models is made

Optical hysteresis of ridge waveguide magnetic garnet films and cavities in two-dimensional magneto-photonic crystal slabs

Magnetic iron garnets as well as magnetic photonic crystals are of great interests in magneto-optic applications such as isolators, current captors, circulators, TE-TM mode conversion, wavelength accordable filters, optical sensors and switches, all of which provide a promising platform for future integrated optical circuits. In the present work, two topics are studied based on magnetic iron garnet films. In the first part, the characteristics of the magnetization are investigated for ridge waveguides fabricated on (100) oriented iron garnet thin films. The magnetic response in magneto-optic waveguides patterned on epitaxial magnetic garnet films depends on the crystallographic orientation of the waveguides and the magnetic anisotropy of the material. These can be studied by polarization rotation hysteresis loops, which are related to the component of magnetization parallel to the light propagation direction and the linear birefringence. Polarization rotation hysteresis loops for low birefringence waveguides with different orientations are experimentally investigated. Asymmetric stepped curves are obtained from waveguides along, due to the large magnetocrystalline anisotropy in the plane. A model based on the free energy density is developed to demonstrate the motion of the magnetization and can be used in the design of magneto-optic devices. The second part of this thesis focuses on the design and fabrication of high-Q cavities in two-dimensional magneto-photonic crystal slabs. The device consists of a layer of silicon and a layer of iron garnet thin film. Triangular lattice elliptical air holes are patterned in the slab. The fundamental TM band gap overlaps with the first-order TE band gap from 0374~0.431(a/λ) showing that both TE and TM polarization light can be confined in the photonic crystals. A nanocavity is designed to obtain both TE and TM defect modes in the band gaps. Additional work is needed to overlap the TE and TM defect modes and obtain a high-Q cavity so as to develop miniaturized Faraday rotators.

Anomalies in stiffness and damping of a 2D discrete viscoelastic system due to negative stiffness components

The recent development of using negative stiffness inclusions to achieve extreme overall stiffness and mechanical damping of composite materials reveals a new avenue for constructing high performance materials. One of the negative stiffness sources can be obtained from phase transforming materials in the vicinity of their phase transition, as suggested by the Landau theory. To understand the underlying mechanism from a microscopic viewpoint, we theoretically analyze a 2D, nested triangular lattice cell with pre-chosen elements containing negative stiffness to demonstrate anomalies in overall stiffness and damping. Combining with current knowledge from continuum models, based on the composite theory, such as the Voigt, Reuss, and Hashin-Shtrikman model, we further explore the stability of the system with Lyapunov's indirect stability theorem. The evolution of the microstructure in terms of the discrete system is discussed. A potential application of the results presented here is to develop special thin films with unusual in-plane mechanical properties. © 2006 Elsevier B.V. All rights reserved.

The correlation of structure and electronic properties near the surface of transition metal oxides

The strong couplings between different degrees of freedom are believed to be responsible for novel and complex phenomena discovered in transition metal oxides (TMOs). The physical complexity is directly responsible for their tunability. Creating surfaces/interfaces add an additional ' man-made' twist, approaching the quantum phenomena of correlated materials. ^ The dissertation focused on the structural and electronic properties in proximity of surface of three prototype TMO compounds by using three complementary techniques: scanning tunneling microscopy, angle-resolved photoelectron spectroscopy and low energy electron diffraction, particularly emphasized the effects of broken symmetry and imperfections like defects on the coupling between charge and lattice degrees of freedom. ^ Ca1.5Sr0.5RuO4 is a layered ruthenate with square lattice and at the boundary of magnetic/orbital instability in Ca2-xSrxRuO4. That the substitution of Sr 2+ with Ca2+ causing RuO6 rotation narrows the dxy band width and changes the Fermi surface topology. Particularly, the γ(dxy) Fermi surface sheet exhibited hole-like in Ca1.5Sr0.5RuO4 in contrast to electron-like in Sr2RuO4, showing a strong charge-lattice coupling. ^ Na0.75CoO2 is a layered cobaltite with triangular lattice exhibiting extraordinary thermoelectric properties. The well-ordered CoO2-terminated surface with random Na distribution was observed. However, lattice constants of the surface are smaller than that in bulk. The surface density of states (DOS) showed strong temperature dependence. Especially, an unusual shift of the minimum DOS occurs below 230 K, clearly indicating a local charging effect on the surface. ^ Cd2Re2O7 is the first known pyrochlore oxide superconductor (Tc ∼ 1K). It exhibited an unusual second-order phase transition occurring at TS1 = 200 K and a controversial first-order transition at TS2 = 120 K. While bulk properties display large anomalies at TS1 but rather subtle and sample-dependent changes at TS2, the surface DOS near the EF show no change at T s1 but a substantial increase below TS2---a complete reversal as the signature for the transitions. We argued that crystal imperfections, mainly defects, which were considerably enhanced at the surface, resulted in the transition at TS2. ^

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.

Study of charge, spin, and structural orderings in quantum materials using nuclei and muons as local probes

Quantum Materials are many body systems displaying emergent phenomena caused by quantum collective behaviour, such as superconductivity, charge density wave, fractional hall effect, and exotic magnetism. Among quantum materials, two families have recently attracted attention: kagome metals and Kitaev materials. Kagome metals have a unique crystal structure made up of triangular lattice layers that are used to form the kagome layer. Due to superconductivity, magnetism, and charge ordering states such as the Charge Density Wave (CDW), unexpected physical phenomena such as the massive Anomalous Hall Effect (AHE) and possible Majorana fermions develop in these materials. Kitaev materials are a type of quantum material with a unique spin model named after Alexei Kitaev. They include fractional fluctuations of Majorana fermions and non-topological abelian anyons, both of which might be used in quantum computing. Furthermore, they provide a realistic framework for the development of quantum spin liquid (QSL), in which quantum fluctuations produce long-range entanglements between electronic states despite the lack of classical magnetic ordering. In my research, I performed several nuclear magnetic resonance (NMR), nuclear quadrupole resonance (NQR), and muon spin spectroscopy (µSR) experiments to explain and unravel novel phases of matter within these unusual families of materials. NMR has been found to be an excellent tool for studying these materials’ local electronic structures and magnetic properties. I could use NMR to determine, for the first time, the structure of a novel kagome superconductor, RbV3Sb5, below the CDW transition, and to highlight the role of chemical doping in the CDW phase of AV3Sb5 superconductors. µSR has been used to investigate the effect of doping on kagome material samples in order to study the presence and behaviour of an anomalous phase developing at low temperatures and possibly related to time-reversal symmetry breaking.