Can model Hamiltonians describe the electron–electron interaction in π-conjugated systems?: PAH and graphene

Model Hamiltonians have been, and still are, a valuable tool for investigating the electronic structure of systems for which mean field theories work poorly. This review will concentrate on the application of Pariser–Parr–Pople (PPP) and Hubbard Hamiltonians to investigate some relevant properties of polycyclic aromatic hydrocarbons (PAH) and graphene. When presenting these two Hamiltonians we will resort to second quantisation which, although not the way chosen in its original proposal of the former, is much clearer. We will not attempt to be comprehensive, but rather our objective will be to try to provide the reader with information on what kinds of problems they will encounter and what tools they will need to solve them. One of the key issues concerning model Hamiltonians that will be treated in detail is the choice of model parameters. Although model Hamiltonians reduce the complexity of the original Hamiltonian, they cannot be solved in most cases exactly. So, we shall first consider the Hartree–Fock approximation, still the only tool for handling large systems, besides density functional theory (DFT) approaches. We proceed by discussing to what extent one may exactly solve model Hamiltonians and the Lanczos approach. We shall describe the configuration interaction (CI) method, a common technology in quantum chemistry but one rarely used to solve model Hamiltonians. In particular, we propose a variant of the Lanczos method, inspired by CI, that has the novelty of using as the seed of the Lanczos process a mean field (Hartree–Fock) determinant (the method will be named LCI). Two questions of interest related to model Hamiltonians will be discussed: (i) when including long-range interactions, how crucial is including in the Hamiltonian the electronic charge that compensates ion charges? (ii) Is it possible to reduce a Hamiltonian incorporating Coulomb interactions (PPP) to an 'effective' Hamiltonian including only on-site interactions (Hubbard)? The performance of CI will be checked on small molecules. The electronic structure of azulene and fused azulene will be used to illustrate several aspects of the method. As regards graphene, several questions will be considered: (i) paramagnetic versus antiferromagnetic solutions, (ii) forbidden gap versus dot size, (iii) graphene nano-ribbons, and (iv) optical properties.

Noncollinear magnetic phases and edge states in graphene quantum Hall bars

Application of a perpendicular magnetic field to charge neutral graphene is expected to result in a variety of broken symmetry phases, including antiferromagnetic, canted, and ferromagnetic. All these phases open a gap in bulk but have very different edge states and noncollinear spin order, recently confirmed experimentally. Here we provide an integrated description of both edge and bulk for the various magnetic phases of graphene Hall bars making use of a noncollinear mean field Hubbard model. Our calculations show that, at the edges, the three types of magnetic order are either enhanced (zigzag) or suppressed (armchair). Interestingly, we find that preformed local moments in zigzag edges interact with the quantum spin Hall like edge states of the ferromagnetic phase and can induce backscattering.

Majorana Zero Modes in Graphene

A clear demonstration of topological superconductivity (TS) and Majorana zero modes remains one of the major pending goals in the field of topological materials. One common strategy to generate TS is through the coupling of an s-wave superconductor to a helical half-metallic system. Numerous proposals for the latter have been put forward in the literature, most of them based on semiconductors or topological insulators with strong spin-orbit coupling. Here, we demonstrate an alternative approach for the creation of TS in graphene-superconductor junctions without the need for spin-orbit coupling. Our prediction stems from the helicity of graphene’s zero-Landau-level edge states in the presence of interactions and from the possibility, experimentally demonstrated, of tuning their magnetic properties with in-plane magnetic fields. We show how canted antiferromagnetic ordering in the graphene bulk close to neutrality induces TS along the junction and gives rise to isolated, topologically protected Majorana bound states at either end. We also discuss possible strategies to detect their presence in graphene Josephson junctions through Fraunhofer pattern anomalies and Andreev spectroscopy. The latter, in particular, exhibits strong unambiguous signatures of the presence of the Majorana states in the form of universal zero-bias anomalies. Remarkable progress has recently been reported in the fabrication of the proposed type of junctions, which offers a promising outlook for Majorana physics in graphene systems.

Strain-induced structural instability in FeRh

We perform density functional calculations to investigate the structure of the intermetallic alloy FeRh under epitaxial strain. Bulk FeRh exhibits a metamagnetic transition from a low-temperature antiferromagnetic (AFM) phase to a ferromagnetic phase at 350 K, and its strain dependence is of interest for tuning the transition temperature to the room-temperature operating conditions of typical memory devices. We find an unusually strong dependence of the structural energetics on the choice of exchange-correlation functional, with the usual local density approximation yielding the wrong ground-state structure, and generalized gradient (GGA) extensions being in better agreement with the bulk experimental structure. Using the GGA we show the existence of a metastable face-centered-cubic-like AFM structure that is reached from the ground-state body-centered-cubic-like AFM structure by following the epitaxial Bain path. We show that the behavior is well described using nonlinear elasticity theory, which captures the softening and eventual sign change of the orthorhombic shear modulus under compressive strain, consistent with this structural instability. Finally, we predict the existence of an additional unit-cell-doubling lattice instability, which should be observable at low temperature.

Vitrinite reflectance of ODP Leg 139 samples (Table 1)

Isolated kerogens from four sites in the Middle Valley hydrothermal region of the Juan de Fuca Ridge (Ocean Drilling Project (ODP) Leg 139) were analysed by electron paramagnetic resonance (EPR) spectroscopy. Measurements of peak width, spin density and power saturation for site 857 kerogens, which increased regularly in maturity downhole, show correlation with vitrinite reflectance values from 0.61 to 2.5%, indicating the start of the oil window at depths from 200 to 400 m. Spin density increases to 1.56 * 10**17 spins per gram and peak width decreases to 3.45 G (gauss) with increasing depth. The tendency to power saturate also decreases with increasing maturity and increasing vitrinite reflectance within the oil window. These trends are consistent with a model in which exchange processes are occurring and cause changes in the EPR behavior of samples from this site. Sediments from other Middle Valley sites, 855, 856 and 858 contain large quantities of pyrite with Mn2+ impurities which interact with the carbon radical to distort the EPR measurements.

Mode Localization in a Nearly Cyclic Symmetric System

Thesis (Ph.D.)--University of Washington, 2016-06

Optimal measurements for relative quantum information

We provide optimal measurement schemes for estimating relative parameters of the quantum state of a pair of spin systems. We prove that the optimal measurements are joint measurements on the pair of systems, meaning that they cannot be achieved by local operations and classical communication. We also demonstrate that in the limit where one of the spins becomes macroscopic, our results reproduce those that are obtained by treating that spin as a classical reference direction.

Entanglement sharing and decoherence in the spin-bath

The monogamous nature of entanglement has been illustrated by the derivation of entanglement-sharing inequalities-bounds on the amount of entanglement that can be shared among the various parts of a multipartite system. Motivated by recent studies of decoherence, we demonstrate an interesting manifestation of this phenomena that arises in system-environment models where there exists interactions between the modes or subsystems of the environment. We investigate this phenomenon in the spin-bath environment, constructing an entanglement-sharing inequality bounding the entanglement between a central spin and the environment in terms of the pairwise entanglement between individual bath spins. The relation of this result to decoherence will be illustrated using simplified system-bath models of decoherence.

Quantum entanglement and fixed-point bifurcations

How does the classical phase-space structure for a composite system relate to the entanglement characteristics of the corresponding quantum system? We demonstrate how the entanglement in nonlinear bipartite systems can be associated with a fixed-point bifurcation in the classical dynamics. Using the example of coupled giant spins we show that when a fixed point undergoes a supercritical pitchfork bifurcation, the corresponding quantum state-the ground state-achieves its maximum amount of entanglement near the critical point. We conjecture that this will be a generic feature of systems whose classical limit exhibits such a bifurcation.

Global control and fast solid-state donor electron spin quantum computing

We propose a scheme for quantum information processing based on donor electron spins in semiconductors, with an architecture complementary to the original Kane proposal. We show that a naive implementation of electron spin qubits provides only modest improvement over the Kane scheme, however through the introduction of global gate control we are able to take full advantage of the fast electron evolution timescales. We estimate that the latent clock speed is 100-1000 times that of the nuclear spin quantum computer with the ratio T-2/T-ops approaching the 10(6) level.

Spin exchange and superconductivity in a t-J(')-V model for two-dimensional quarter-filled systems

The effect of antiferromagnetic spin fluctuations on two-dimensional quarter-filled systems is studied theoretically. An effective t-J(')-V model on a square lattice which accounts for checkerboard charge fluctuations and next-nearest-neighbor antiferromagnetic spin fluctuations is considered. From calculations based on large-N theory on this model it is found that the exchange interaction J(') increases the attraction between electrons in the d(xy) channel only, so that both charge and spin fluctuations work cooperatively to produce d(xy) pairing.

Measurement-based teleportation along quantum spin chains

We examine the teleportation of an unknown spin-1/2 quantum state along a quantum spin chain with an even number of sites. Our protocol, using a sequence of Bell measurements, may be viewed as an iterated version of the 2-qubit protocol of C. H. Bennett et al. [Phys. Rev. Lett. 70, 1895 (1993)]. A decomposition of the Hilbert space of the spin chain into 4 vector spaces, called Bell subspaces, is given. It is established that any state from a Bell subspace may be used as a channel to perform unit fidelity teleportation. The space of all spin-0 many-body states, which includes the ground states of many known antiferromagnetic systems, belongs to a common Bell subspace. A channel-dependent teleportation parameter O is introduced, and a bound on the teleportation fidelity is given in terms of O.

On the chaotic instability of a nonsliding liquid-filled top with a small spheroidal base via Melnikov-Holmes-Marsden integrals

Chaotic orientations of a top containing a fluid filled cavity are investigated analytically and numerically under small perturbations. The top spins and rolls in nonsliding contact with a rough horizontal plane and the fluid in the ellipsoidal shaped cavity is considered to be ideal and describable by finite degrees of freedom. A Hamiltonian structure is established to facilitate the application of Melnikov-Holmes-Marsden (MHM) integrals. In particular, chaotic motion of the liquid-filled top is identified to be arisen from the transversal intersections between the stable and unstable manifolds of an approximated, disturbed flow of the liquid-filled top via the MHM integrals. The developed analytical criteria are crosschecked with numerical simulations via the 4th Runge-Kutta algorithms with adaptive time steps.

Effects of cation substitutions on the physical properties of M2M1T2O6 pyroxenes

In this PhD study, the effects of the cation substitutions on the physical properties of pyroxenes have been discussed. The results of this work extend the knowledge on pyroxenes with different chemical compositions. These properties might be used in the development of ceramic pigments, advanced materials and for the mineralogical phase identification. First of all, the crystallographic differences between Ge and Si pyroxenes have been examined. The structure of C2/c Ca rich Ge clinopyroxenes is very close to the low pressure C2/c structural configuration found in Ca-rich Si-pyroxenes. The shear of the unit cell is very similar, and the difference between a Ge end member and the corresponding Si-rich one is less than 1°. Instead, a remarkable difference exists between Ca-poor Si and Ge clinopyroxenes. First, Ca-poor Ge pyroxenes do not display a P21/c symmetry, but retain the C2/c symmetry; second, the observed C2/c structure shows, at room pressure, the configuration with highly kinked tetrahedral chains characteristic of the high pressure C2/c symmetry of Si Ca-poor pyroxenes. In orthopyroxenes, with Pbca symmetry, Ge-pyroxenes have volume larger than Si-pyroxenes. Samples along the system CaCoGe2O6 - CoCoGe2O6 have been synthesized at three different temperatures: 1050 °C, 1200 °C and 1250 °C. The aim of these solid state syntheses was to obtain a solid solution at ambient pressure, since the analogues Si-system needs high pressure. Unfortunately, very limited solution occurs because the structure forms of the two end member (high temperature for CaCoGe2O6 and high pressure CoCoGe2O6) are incompatible. The phase diagram of this system has been sketched and compared to that of Si. The cobalt end member (CoCoGe2O6) is stable at ambient pressure in two symmetries: at 1050 °C C2/c and 1200 °C Pbca. The impurity phase formed during these experiments is cobalt spinel. Raman spectroscopy has been used to investigate the vibrational properties of Ca-pyroxenes CaCoGe2O6, CaMgGe2O6, CaMgSi2O6 and CaCoSi2O6. A comparison between silicate and germanate pyroxenes shows significant changes in peak positions of the corresponding modes caused mainly by the difference of the Ge-Si atomic weight along with the distortion and compression of the coordination polyhedra. Red shift in Raman spectra of germanates has been calculated by a rough scale factor calculated by a simple harmonic oscillator model, considering the different bond lengths for 4-coordinated Si ~ 1.60- 1.65 Å vs Ge–O distance ~1.70 - 1.80 Å. The Raman spectra of CaMgGe2O6 and CaCoGe2O6 have been classified, in analogy with silicate (Wang et al., 2001) counterparts, in different ranges: - R1 (880-640 cm-1): strong T-O stretching modes of Ge and non-bridging O1 and O2 atoms within the GeO4 tetrahedron; - R2 (640-480 cm-1): stretching/bending modes of Ge-Obr-Ge bonds (chain stretching and chain bending); - R4 (480-360 cm-1): O-Ge-O vibrations; - R3 (360-240 cm-1): motions of the cations in M2 and M1 sites correlated with tetrahedral chain motion and tilting tetrahedra; - R5 (below 240 cm-1): lattice modes. The largest shift with respect to CaMgSi2O6 - CaCoSi2O6 is shown by the T-O stretching and chain modes. High-pressure Raman spectroscopy (up to about 8 GPa) on the same samples of Ca-pyroxenes using an ETH-type diamond anvil cell shows no phase transition within the P-ranges investigated, as all the peak positions vary linearly as a function of pressure. Our data confirm previous experimental findings on Si-diopside (Chopelas and Serghiou, 2000). In the investigated samples, all the Raman peaks shift upon compression, but the major changes in wavenumber with pressure are attributed to the chain bending (Ge-Obr-Ge bonds) and tetrahedra stretching modes (Ge-Onbr). Upon compression, the kinking angle, the bond lengths and T-T distances between tetrahedra decrease and consequently the wavenumber of the bending chain mode and tetrahedra stretching mode increases. Ge-pyroxenes show the higher P-induced peak-position shifts, being more compressible than corresponding silicates. The vibrational properties of CaM2+Ge2O6 (M2+ =Mg, Mn, Fe, Co, Ni, Zn) are reported for the first time. The wavenumber of Ge-Obr-Ge bending modes decreases linearly with increasing ionic radius of the M1 cation. No simple correlation has been found with M1 atomic mass or size or crystallographic parameters for the peak at ~850 cm-1 and in the low wavenumber regions. The magnetic properties of the system CaCoSi2O6 - CoCoSi2O6 have been investigated by magnetometry. The join is always characterized by 1 a.p.f.u. of cobalt in M1 site and this causes a pure collinear antiferromagnetic behaviour of the intra-chain superexchange interaction involving Co ions detected in all the measurements, while the magnetic order developed by the cobalt ions in M2 site (intra-chain) is affected by weak ferromagnetism, due to the non-collinearity of their antiferromagnetic interaction. In magnetically ordered systems, this non-collinearity effect promotes a spin canting of anti-parallel aligned magnetic moments and thus is a source of weak ferromagnetic behaviour in an antiferromagnetic. The weak ferromagnetism can be observed only for the samples with Co content higher than 0.5 a.p.f.u. in M2, when the concentration is sufficiently high to create a long range order along the M2 chain which is magnetically independent of M1 chain. The ferromagnetism was detected both in the M(T) at 10 Oe and M(H).

Magnetic properties and Spin dynamics in magnetic Molecular Rings

Molecular nanomagnets are spin clusters whose topology and magnetic interactions can be modulated at the level of the chemical synthesis. They are formed by a small number of transition metal ions coupled by the Heisenberg's exchange interactions. Each cluster is magnetically isolated from its neighbors by organic ligands, making each unit not interacting with the others. Therefore, we can investigate the magnetic properties of an isolated molecular nanomagnet by bulk measurements. The present thesis has been mostly devoted to the experimental investigation of the magnetic properties and spin dynamics of different classes of antiferromagnetic (AF) molecular rings. This study has been exploiting various techniques of investigations, such as Nuclear Magnetic Resonance (NMR), muon spin relaxation (muSR) and SQUiD magnetometry. We investigate the magnetic properties and the phonon-induced relaxation dynamics of the first regular Cr9 antiferromagnetic (AF) ring, which represents a prototype frustrated AF ring. The magnetically-open AF rings like Cr8Cd are model systems for the study of the microscopic magnetic behaviour of finite AF Heisenberg chains. In this type of system the different magnetic behaviour depends length and on the parity of the chain (odd or even). In order to study the local spin densities on the Cr sites, the Cr-NMR spectra was collected at low temperature. The experimental result confirm the theoretical predictions for the spin configuration. Finally, the study of Dy6, the first rare-earth based ring that has been ever synthesized, has been performed by AC-SQuID and muSR measurements. We found that the dynamics is characterized by more than one characteristic correlation time, whose values depend strongly on the applied field.