Showing 1–50 of 55 results for author: Nilsen, J

Kitaev Interactions Through an Extended Superexchange Pathway in the jeff = 1/2 Ru3+ Honeycomb Magnet, RuP3SiO11

Authors: Aly H. Abdeldaim, Hlynur Gretarsson, Sarah J. Day, M. Duc Le, Gavin B. G. Stenning, Pascal Manuel, Robin S. Perry, Alexander A. Tsirlin, Gøran J. Nilsen, Lucy Clark

Abstract: Magnetic materials are composed of the simple building blocks of magnetic moments on a crystal lattice that interact via magnetic exchange. Yet from this simplicity emerges a remarkable diversity of magnetic states. Some reveal the deep quantum mechanical origins of magnetism, for example, quantum spin liquid (QSL) states in which magnetic moments remain disordered at low temperatures despite bein… ▽ More Magnetic materials are composed of the simple building blocks of magnetic moments on a crystal lattice that interact via magnetic exchange. Yet from this simplicity emerges a remarkable diversity of magnetic states. Some reveal the deep quantum mechanical origins of magnetism, for example, quantum spin liquid (QSL) states in which magnetic moments remain disordered at low temperatures despite being strongly correlated through quantum entanglement. A promising theoretical model of a QSL is the Kitaev model, composed of unusual bond-dependent exchange interactions, but experimentally, this model is challenging to realise. Here we show that the material requirements for the Kitaev QSL survive an extended pseudo-edge-sharing superexchange pathway of Ru3+ octahedra within the honeycomb layers of the inorganic framework solid, RuP3SiO11. We confirm the requisite jeff = 1/2 state of Ru3+ in RuP3SiO11 and resolve the hierarchy of exchange interactions that provide experimental access to an unexplored region of the Kitaev model. △ Less

Submitted 14 October, 2024; v1 submitted 28 March, 2024; originally announced March 2024.

Comments: 28 Pages, 5 figures, 1 table

Quasi van der Waals Epitaxial Growth of GaAsSb Nanowires on Graphitic Substrate for Photonic Applications

Authors: Dingding Ren, Tron A. Nilsen, Julie S. Nilsen, Lyubomir Ahtapodov, Anjan Mukherjee, Yang Li, Antonius T. J. van Helvoort, Helge Weman, Bjørn-Ove Fimland

Abstract: III-V semiconductor nanowires are considered promising building blocks for advanced photonic devices. One of the key advantages is that the lattice mismatch can easily be accommodated in 1D structures, resulting in superior heteroepitaxial quality compared to thin films. However, few reports break the limitation of using bulk crystalline materials as substrates for epitaxial growth of high-quality… ▽ More III-V semiconductor nanowires are considered promising building blocks for advanced photonic devices. One of the key advantages is that the lattice mismatch can easily be accommodated in 1D structures, resulting in superior heteroepitaxial quality compared to thin films. However, few reports break the limitation of using bulk crystalline materials as substrates for epitaxial growth of high-quality photonic 1D components, making monolithic integration of III-V components on arbitrary substrates challenging. In this work, we show that the growth of self-catalyzed GaAsSb nanowires on graphitic substrates can be promoted by creating step edges of monolayer thickness on kish graphite before the growth. By further alternating the deposition sequence of the group-III element Al and the group-V elements As and Sb, it was found that triangular crystallites form when Al is deposited first. This indicates that the surface binding energy between the graphitic surface and the III-V nucleus profoundly influences the epitaxial growth of III-V materials on graphitic surfaces. Using the optimized growth recipe with an AlAsSb buffer nuclei, vertical [111]-oriented GaAsSb/GaAs nanowires with GaAsSb-based multiple axial superlattices were grown on exfoliated graphite, which was attached to a (001) AlAs/GaAs distributed Bragg reflector (DBR) using the simple Scotch tape method. Fabry-Pérot resonance modes were observed under optical excitation at room temperature, indicating a successful monolithic integration with optical feedback from the DBR system. These results demonstrate the great potential for flexible integration of high-efficiency III-V nanowire photonic devices on arbitrary photonic platforms using a 2D material buffer layer, e.g., graphene, without breaking the orientation registry. △ Less

Submitted 9 February, 2024; originally announced February 2024.

Inelastic neutron scattering investigation of the crystal field excitations of NdCo$_5$

Authors: F. de Almeida Passos, G. J. Nilsen, C. E. Patrick, M. D. Le, G. Balakrishnan, Santosh Kumar, A. Thamizhavel, D. R. Cornejo, J. Larrea Jiménez

Abstract: We present an inelastic neutron scattering study of the crystal electric field levels in the intermetallic ferrimagnets RECo$_{5}$ (RE = Nd and Y). In NdCo$_{5}$, measurements at $5~$K reveal two levels at approximately 28.9 and 52.9 meV. Crystal field calculations including the exchange field $B_{\textrm{exc}}$ from the Co sites account for both of these, as well as the spectrum at temperatures a… ▽ More We present an inelastic neutron scattering study of the crystal electric field levels in the intermetallic ferrimagnets RECo$_{5}$ (RE = Nd and Y). In NdCo$_{5}$, measurements at $5~$K reveal two levels at approximately 28.9 and 52.9 meV. Crystal field calculations including the exchange field $B_{\textrm{exc}}$ from the Co sites account for both of these, as well as the spectrum at temperatures above the spin-reorientation transition at $\sim 280$~K. In particular, it is found that both a large hexagonal crystal field parameter $A_{6}^6\langle r^6 \rangle$ and $B_{\textrm{exc}}$ are required to reproduce the data, with the latter having a much larger value than that deduced from previous computational and experimental studies. Our study sheds light on the delicate interplay of terms in the rare-earth Hamiltonian of RECo$_5$ systems, and is therefore expected to stimulate further experimental and computational work on the broader family of rare-earth permanent magnets. △ Less

Submitted 1 June, 2023; originally announced June 2023.

Dynamics of the antiferromagnetic spin ice phase in pyrochlore spinels

Authors: Attila Szabó, Gøran J. Nilsen

Abstract: Motivated by the classical spin-nematic state observed in the breathing pyrochlore spinel LiGa$_{0.95}$In$_{0.05}$Cr$_4$O$_8$, we theoretically discuss spin dynamics in models of spin-lattice coupling in these materials. Semiclassical dynamical simulations successfully recover the key features of inelastic neutron-scattering experiments on LiGa$_{0.95}$In$_{0.05}$Cr$_4$O$_8$: a broad finite-energy… ▽ More Motivated by the classical spin-nematic state observed in the breathing pyrochlore spinel LiGa$_{0.95}$In$_{0.05}$Cr$_4$O$_8$, we theoretically discuss spin dynamics in models of spin-lattice coupling in these materials. Semiclassical dynamical simulations successfully recover the key features of inelastic neutron-scattering experiments on LiGa$_{0.95}$In$_{0.05}$Cr$_4$O$_8$: a broad finite-energy peak alongside a continuum of scattering near the (200) wave vector that extends from the elastic line to high energies. To interpret this result, we develop a small-fluctuation theory for the spin-ice-like nematic ground states, analogous to linear spin-wave theory for conventionally ordered magnets, which reproduces the numerical simulation results quantitatively. In particular, we find that the inelastic peak is well explained by collective modes confined to ferromagnetic loops of the underlying nematic order. In addition, we find a sharp, linearly dispersing mode in the dynamical structure factor, which originates in long-wavelength fluctuations of the nematic director: We believe identifying this mode will be an interesting target for future experiments on these materials. We also outline potential future applications of our methods to both pyrochlore spinels and other spin-nematic systems. △ Less

Submitted 25 May, 2023; originally announced May 2023.

Comments: 14 pages, 11 figures

Inelastic Neutron Scattering Study of the Spin Dynamics in the Breathing Pyrochlore System LiGa0.95In0.05Cr4O8

Authors: Yu Tanaka, Rafał Wawrzyńczak, Manh Duc Le, Tatiana Guidi, Yoshihiko Okamoto, Takeshi Yajima, Zenji Hiroi, Masashi Takigawa, Gøran J. Nilsen

Abstract: The A-site ordered chromate spinels LiGa1-xInxCr4O8 host a network of size-alternating spin-3/2 Cr3+ tetrahedra known as a 'breathing' pyrochlore lattice. For the x=0.05 composition, the complex magneto-structural ordering observed in the parent x=0 material is replaced by a single transition at Tf=11 K, ascribed to the collinear nematic order caused by strong spin-lattice coupling. We present her… ▽ More The A-site ordered chromate spinels LiGa1-xInxCr4O8 host a network of size-alternating spin-3/2 Cr3+ tetrahedra known as a 'breathing' pyrochlore lattice. For the x=0.05 composition, the complex magneto-structural ordering observed in the parent x=0 material is replaced by a single transition at Tf=11 K, ascribed to the collinear nematic order caused by strong spin-lattice coupling. We present here an inelastic neutron scattering study of the spin dynamics in this composition. Above Tf , the dynamical scattering function S(Q,E) is ungapped and quasi-elastic, similar to undoped LiGaCr4O8. Below Tf , the spectral weight splits between a broad inelastic feature at 5.8 meV and toward the elastic line. The former feature can be ascribed to spin precessions within antiferromagnetic loops, lifted to finite energy by the effective biquadratic spin-lattice term in the spin Hamiltonian. △ Less

Submitted 12 January, 2023; originally announced January 2023.

Journal ref: J. Phys. Soc. Jpn. 87, 073710 (2018)

One-dimensional magnetism in synthetic Pauflerite, $β$-VOSO$_4$

Authors: Diana Lucia Quintero-Castro, Gøran J. Nilsen, Katrin Meier-Kirchner, Angelica Benitez-Castro, Gerrit Guenther, Toshiro Sakakibara, Masashi Tokunaga, Chidozie Agu, Ipsita Mandal, Alexander A. Tsirlin

Abstract: We have synthesized single-crystal samples of $β$-VOSO$_4$ and fully characterized their magnetic properties. Our magnetic susceptibility, high field magnetization and powder inelastic neutron scattering results are in excellent agreement with theoretical expressions for a one-dimensional spin-1/2 Heisenberg chain with an exchange parameter of $3.83(2)$ meV. Ab initio calculations identify the sup… ▽ More We have synthesized single-crystal samples of $β$-VOSO$_4$ and fully characterized their magnetic properties. Our magnetic susceptibility, high field magnetization and powder inelastic neutron scattering results are in excellent agreement with theoretical expressions for a one-dimensional spin-1/2 Heisenberg chain with an exchange parameter of $3.83(2)$ meV. Ab initio calculations identify the superexchange pathway, revealing that the spin-chain does not run along the expected crystallographic chain $a$ direction but instead between V$^{4+}$O$_{6}$ octahedra that are linked via SO$_{4}$ tetrahedra along the $b$ axis. We do not detect any phase transition to a long-range magnetic order within our experimental conditions, indicating $β$-VOSO$_4$ is very close to an ideal one-dimensional magnetic system. △ Less

Submitted 29 April, 2023; v1 submitted 16 December, 2022; originally announced December 2022.

Journal ref: Phys. Rev. Materials 7, 045003 (2023)

Magnetic, electronic, and structural investigation of the strongly correlated Y$_{1-x}$Sm$_x$Co$_5$ system

Authors: Fernando de Almeida Passos, Daniel Cornejo, Márcia Fantini, Gøran J. Nilsen, Valentina Martelli, Julio Larrea Jiménez

Abstract: SmCo$_5$ and YCo$_5$ are isostructural compounds both showing large magnetocrystalline anisotropy, where the former originates mainly from the crystal-electric field and magnetic interactions. We investigate the contribution of both interactions by partially substituting Y by Sm in as-cast polycrystalline Y$_{1-x}$Sm$_{x}$Co$_5$ (with $x$=0, 0.1, 0.2, 0.3, and 0.4) and measuring their structural,… ▽ More SmCo$_5$ and YCo$_5$ are isostructural compounds both showing large magnetocrystalline anisotropy, where the former originates mainly from the crystal-electric field and magnetic interactions. We investigate the contribution of both interactions by partially substituting Y by Sm in as-cast polycrystalline Y$_{1-x}$Sm$_{x}$Co$_5$ (with $x$=0, 0.1, 0.2, 0.3, and 0.4) and measuring their structural, magnetic, and electrical properties through X-ray diffraction, magnetization, and electrical transport measurements. Our results suggest an interplay between microstructure strain in as-cast samples and the electronic and magnetic interactions. △ Less

Submitted 17 August, 2022; originally announced August 2022.

Journal ref: SciPost Phys. Proc. 11, 021 (2023)

One-Dimensional Quantum Magnetism in the S = 1/2 Mo(V) system, KMoOP2O7

Authors: Aly H. Abdeldaim, Alexander A. Tsirlin, Jacques Ollivier, Clemens Ritter, Dominic Fortes, Robin S. Perry, Lucy Clark, Gøran J. Nilsen

Abstract: We present a comprehensive experimental and ab-initio study of the $S=1/2$ Mo$^{5+}$ system, KMoOP$_2$O$_7$, and show that it realizes the $S = 1/2$ Heisenberg chain antiferromagnet model. Powder neutron diffraction reveals that KMoOP$_2$O$_7$ forms a magnetic network comprised of pairs of Mo$^{5+}$ chains within its monoclinic $P2_1/n$ structure. Antiferromagnetic interactions within the Mo… ▽ More We present a comprehensive experimental and ab-initio study of the $S=1/2$ Mo$^{5+}$ system, KMoOP$_2$O$_7$, and show that it realizes the $S = 1/2$ Heisenberg chain antiferromagnet model. Powder neutron diffraction reveals that KMoOP$_2$O$_7$ forms a magnetic network comprised of pairs of Mo$^{5+}$ chains within its monoclinic $P2_1/n$ structure. Antiferromagnetic interactions within the Mo$^{5+}$ chains are identified through magnetometry measurements and confirmed by analysis of the magnetic specific heat. The latter reveals a broad feature centred on $T_\textrm{N} = 0.54$ K, which we ascribe to the onset of long-range antiferromagnetic order. No magnetic Bragg scattering is observed in powder neutron diffraction data collected at 0.05 K, however, which is consistent with a strongly suppressed ordered moment with an upper limit $μ_\textrm{ord} < 0.15 μ_\textrm{B}$. The one-dimensional character of the magnetic correlations in KMoOP$_2$O$_7$ is verified through analysis of inelastic neutron scattering data, resulting in a model with $J_\textrm{1} \approx 34$ K and $J_\textrm{2} \approx -2$ K for the intrachain and interchain exchange interactions, respectively. The origin of these experimental findings are addressed through density-functional theory calculations. △ Less

Submitted 14 November, 2022; v1 submitted 4 July, 2022; originally announced July 2022.

Comments: 10 pages, 7 figures

Shock Hugoniot of diamond from 3 to 80 TPa

Authors: Damian C. Swift, Andrea L. Kritcher, Amy Lazicki, James A. Hawreliak, Tilo Doeppner, Heather D. Whitley, Joseph Nilsen, Benjamin Bachmann, Michael MacDonald, Brian Maddox, Natalie Kostinski, Gilbert W. Collins, Siegfried Glenzer, Stephen D. Rothman, Dominik Kraus, Roger W. Falcone

Abstract: The principal Hugoniot of carbon, initially diamond, was measured from 3 to 80 TPa (30 to 800 million atmospheres), the highest pressure ever achieved, using radiography of spherically-converging shocks. The shocks were generated by ablation of a plastic coating by soft x-rays in a laser-heated hohlraum at the National Ignition Facility (NIF). Experiments were performed with low and high drive pow… ▽ More The principal Hugoniot of carbon, initially diamond, was measured from 3 to 80 TPa (30 to 800 million atmospheres), the highest pressure ever achieved, using radiography of spherically-converging shocks. The shocks were generated by ablation of a plastic coating by soft x-rays in a laser-heated hohlraum at the National Ignition Facility (NIF). Experiments were performed with low and high drive powers, spanning different but overlapping pressure ranges. The radius-time history of the shock, and the profile of mass density behind, were determined by profile-matching from a time-resolved x-ray radiograph across the diameter of the sphere. Above ~50 TPa, the heating induced by the shock was great enough to ionize a significant fraction of K-shell electrons, reducing the opacity to the 10.2 keV probe x-rays. The opacity and mass density were deduced simultaneously using the constraint that the total mass of the sample was constant. The Hugoniot and opacity were consistent with density functional theory calculations of the electronic states and equation of state (EOS), and varied significantly from theoretical Hugoniots based on Thomas-Fermi theory. Theoretical models used to predict the compressibility of diamond ablator experiments at the NIF, producing the highest neutron yields so far from inertial confinement fusion experiments, are qualitatively consistent with our EOS measurements but appear to overpredict the compressibility slightly. These measurements help to evaluate theoretical techniques and constrain wide-range EOS models applicable to white dwarf stars, which are the ultimate evolutionary form of at least 97% of stars in the galaxy. △ Less

Submitted 30 March, 2022; originally announced March 2022.

Comments: arXiv admin note: text overlap with arXiv:2203.08891

Report number: LLNL-JRNL-815773

Improved analysis of converging shock radiographs

Authors: Damian C. Swift, Andrea L. Kritcher, Amy Lazicki, Natalie Kostinski, Brian R. Maddox, Madison E. Martin, Tilo Doeppner, Joseph Nilsen, Heather D. Whitley

Abstract: We previously reported an experimental platform to induce a spherically-convergent shock in a sample using laser-driven ablation, probed with time-resolved x-ray radiography, and an analysis method to deduce states along the principal shock Hugoniot simultaneously with the x-ray opacity. We have now developed a modified method of analysis that is numerically better-conditioned and faster, and usua… ▽ More We previously reported an experimental platform to induce a spherically-convergent shock in a sample using laser-driven ablation, probed with time-resolved x-ray radiography, and an analysis method to deduce states along the principal shock Hugoniot simultaneously with the x-ray opacity. We have now developed a modified method of analysis that is numerically better-conditioned and faster, and usually provides a better representation of the radiograph with correspondingly lower uncertainties. The previous approach was based on optimizing parameters in a model of the density distribution as a function of radius and time, warped to follow loci such as the shock and the outside of the sample. The converging shock configuration can be described more efficiently in terms of the shocked density and sound speed, expressed as functions of the shock speed Studies of the Hugoniot from various theoretical equations of state (EOS) indicate that, in the typical range of states explored by these experiments, these functions can be described by low-order polynomials. Similarly, few-parameter functions were found suitable for representing the variation of x-ray opacity with shock pressure. This approach was found to perform better in most cases than an alternative method based on parameterization of the EOS. △ Less

Submitted 16 March, 2022; originally announced March 2022.

Report number: LLNL-JRNL-832834

Magnetic frustration in a van der Waals metal CeSiI

Authors: Ryutaro Okuma, Clemens Ritter, Gøran J. Nilsen, Yoshinori Okada

Abstract: The realization of magnetic frustration in a metallic van der Waals (vdW) coupled material has been sought as a promising platform to explore novel phenomena both in bulk matter and in exfoliated devices. However, a suitable material platform has been lacking so far. Here, we demonstrate that CeSiI hosts itinerant electrons coexisting with exotic magnetism. In CeSiI, the magnetic cerium atoms form… ▽ More The realization of magnetic frustration in a metallic van der Waals (vdW) coupled material has been sought as a promising platform to explore novel phenomena both in bulk matter and in exfoliated devices. However, a suitable material platform has been lacking so far. Here, we demonstrate that CeSiI hosts itinerant electrons coexisting with exotic magnetism. In CeSiI, the magnetic cerium atoms form a triangular bilayer structure sandwiched by van der Waals stacked iodine layers. From resistivity and magnetometry measurements, we confirm the coexistence of itinerant electrons with magnetism with dominant antiferromagnetic exchange between the strongly Ising-like Ce moments below 7 K. Neutron diffraction directly confirms magnetic order with an incommensurate propagation vector k ~ (0.28, 0, 0.19) at 1.6 K, which points to the importance of further neighbor magnetic interactions in this system. The presence of a two-step magnetic-field-induced phase transition along c axis further suggests magnetic frustration in the ground state. Our findings provide a novel material platform hosting a coexistence of itinerant electron and frustrated magnetism in a vdW system, where exotic phenomena arising from rich interplay between spin, charge and lattice in low dimension can be explored. △ Less

Submitted 29 July, 2021; originally announced July 2021.

Comments: 15 pages, 9 figures

Journal ref: Phys. Rev. Materials 5, L121401 (2021)

Diffraction by multipoles in a 5d2 rhenium double perovskite

Authors: S. W. Lovesey, D. D. Khalyavin, G. van der Laan, G. J. Nilsen

Abstract: A recent polarized neutron diffraction experiment on the 5d2 rhenium double perovskite Ba2YReO6 held at a low temperature uncovered weak magnetic diffraction peaks. Data analysis inferred a significantly reduced Re dipole moment, and long-range order compatible with an antiferromagnet, non-collinear motif. To interpret the experimental findings, we present a model wavefunction for Re ions derived… ▽ More A recent polarized neutron diffraction experiment on the 5d2 rhenium double perovskite Ba2YReO6 held at a low temperature uncovered weak magnetic diffraction peaks. Data analysis inferred a significantly reduced Re dipole moment, and long-range order compatible with an antiferromagnet, non-collinear motif. To interpret the experimental findings, we present a model wavefunction for Re ions derived from the crystal field potential, Coulomb interaction, and spin-orbit coupling that fully respects the symmetry of the low-temperature ordered state. It is used to calculate in analytic form all multipole moments visible in neutron and resonance enhanced x-ray diffraction. A minimal model consistent with available neutron diffraction data predicts significant multipolar moments up to the hexadecapole, and, in particular, a dominant charge-like quadrupole moment. Calculated diffraction patterns embrace single crystal x-ray diffraction at the Re L-edge, and renewed neutron diffraction, to probe the presumed underlying multipolar order. △ Less

Submitted 8 November, 2020; originally announced November 2020.

Journal ref: Phys. Rev. B 103, 104429 (2021)

Realising square and diamond lattice $S=1/2$ Heisenberg antiferromagnet models in the $α$ and $β$ phases of the coordination framework, KTi(C$_2$O$_4$)$_2\cdot$\textit{x}H$_2$O

Authors: Aly H. Abdeldaim, Teng Li, Lewis Farrar, Alexander A. Tsirlin, Wenjiao Yao, Alexandra S. Gibbs, Pascal Manuel, Philip Lightfoot, Gøran J. Nilsen, Lucy Clark

Abstract: We report the crystal structures and magnetic properties of two psuedo-polymorphs of the $S=1/2$ Ti$^{3+}$ coordination framework, KTi(C$_2$O$_4$)$_2\cdot$xH$_2$O. Single-crystal X-ray and powder neutron diffraction measurements on $α$-KTi(C$_2$O$_4$)$_2\cdot$xH$_2$O confirm its structure in the tetragonal $I4/mcm$ space group with a square planar arrangement of Ti$^{3+}$ ions. Magnetometry and sp… ▽ More We report the crystal structures and magnetic properties of two psuedo-polymorphs of the $S=1/2$ Ti$^{3+}$ coordination framework, KTi(C$_2$O$_4$)$_2\cdot$xH$_2$O. Single-crystal X-ray and powder neutron diffraction measurements on $α$-KTi(C$_2$O$_4$)$_2\cdot$xH$_2$O confirm its structure in the tetragonal $I4/mcm$ space group with a square planar arrangement of Ti$^{3+}$ ions. Magnetometry and specific heat measurements reveal weak antiferromagnetic interactions, with $J_1\approx7$ K and $J_2/J_1=0.11$ indicating a slight frustration of nearest- and next-nearest-neighbor interactions. Below $1.8$ K, $α$ undergoes a transition to G-type antiferromagnetic order with magnetic moments aligned along the $c$ axis of the tetragonal structure. The estimated ordered moment of Ti$^{3+}$ in $α$ is suppressed from its spin-only value to $0.62(3)~μ_B$, thus verifying the two-dimensional nature of the magnetic interactions within the system. $β$-KTi(C$_2$O$_4$)$_2\cdot$2H$_2$O, on the other hand, realises a three-dimensional diamond-like magnetic network of Ti$^{3+}$ moments within a hexagonal $P6_222$ structure. An antiferromagnetic exchange coupling of $J\approx54$ K -- an order of magnitude larger than in $α$ -- is extracted from magnetometry and specific heat data. $β$ undergoes Néel ordering at $T_N=28$ K, with the magnetic moments aligned within the $ab$ plane and a slightly reduced ordered moment of $0.79~μ_B$ per Ti$^{3+}$. Through density-functional theory calculations, we address the origin of the large difference in the exchange parameters between the $α$ and $β$ psuedo-polymorphs. Given their observed magnetic behaviors, we propose $α$-KTi(C$_2$O$_4$)$_2\cdot$xH$_2$O and $β$-KTi(C$_2$O$_4$)$_2\cdot$2H$_2$O as close to ideal model $S=1/2$ Heisenberg square and diamond lattice antiferromagnets, respectively. △ Less

Submitted 7 October, 2020; v1 submitted 7 September, 2020; originally announced September 2020.

Comments: 12 pages, 7 figures

Benchmarking boron carbide equation of state using computation and experiment

Authors: Shuai Zhang, Michelle C. Marshall, Lin H. Yang, Philip A. Sterne, Burkhard Militzer, Markus Daene, James A. Gaffney, Andrew Shamp, Tadashi Ogitsu, Kyle Caspersen, Amy E. Lazicki, David Erskine, Richard A. London, Peter M. Celliers, Joseph Nilsen, Heather D. Whitley

Abstract: Boron carbide (B$_4$C) is of both fundamental scientific and practical interest in inertial confinement fusion (ICF) and high energy density physics experiments. We report the results of a comprehensive computational study of the equation of state (EOS) of B$_4$C in the liquid, warm dense matter, and plasma phases. Our calculations are cross-validated by comparisons with Hugoniot measurements up t… ▽ More Boron carbide (B$_4$C) is of both fundamental scientific and practical interest in inertial confinement fusion (ICF) and high energy density physics experiments. We report the results of a comprehensive computational study of the equation of state (EOS) of B$_4$C in the liquid, warm dense matter, and plasma phases. Our calculations are cross-validated by comparisons with Hugoniot measurements up to 61 megabar from planar shock experiments performed at the National Ignition Facility (NIF). Our computational methods include path integral Monte Carlo, activity expansion, as well as all-electron Green's function Korringa-Kohn-Rostoker and molecular dynamics that are both based on density functional theory. We calculate the pressure-internal energy EOS of B$_4$C over a broad range of temperatures ($\sim$6$\times$10$^3$--5$\times$10$^8$ K) and densities (0.025--50 g/cm$^{3}$). We assess that the largest discrepancies between theoretical predictions are $\lesssim$5% near the compression maximum at 1--2$\times10^6$ K. This is the warm-dense state in which the K shell significantly ionizes and has posed grand challenges to theory and experiment. By comparing with different EOS models, we find a Purgatorio model (LEOS 2122) that agrees with our calculations. The maximum discrepancies in pressure between our first-principles predictions and LEOS 2122 are $\sim$18% and occur at temperatures between 6$\times$10$^3$--2$\times$10$^5$ K, which we believe originate from differences in the ion thermal term and the cold curve that are modeled in LEOS 2122 in comparison with our first-principles calculations. In addition, we have developed three new equation of state models and applied them to 1D hydrodynamic simulations of a polar direct-drive NIF implosion, demonstrating that these new models are now available for future ICF design studies. △ Less

Submitted 3 August, 2020; originally announced August 2020.

Comments: 20 pages, 12 figures, 2 tables

Report number: LLNL-JRNL-812984

Journal ref: Phys. Rev. E 102, 053203 (2020)

Comparison of ablators for the polar direct drive exploding pusher platform

Authors: Heather D. Whitley, G. Elijah Kemp, Charles Yeamans, Zachary Walters, Brent E. Blue, Warren Garbett, Marilyn Schneider, R. Stephen Craxton, Emma M. Garcia, Patrick W. McKenty, Maria Gatu-Johnson, Kyle Caspersen, John I. Castor, Markus Däne, C. Leland Ellison, James Gaffney, Frank R. Graziani, John Klepeis, Natalie Kostinski, Andrea Kritcher, Brandon Lahmann, Amy E. Lazicki, Hai P. Le, Richard A. London, Brian Maddox , et al. (14 additional authors not shown)

Abstract: We examine the performance of pure boron, boron carbide, high density carbon, and boron nitride ablators in the polar direct drive exploding pusher (PDXP) platform. The platform uses the polar direct drive configuration at the National Ignition Facility to drive high ion temperatures in a room temperature capsule and has potential applications for plasma physics studies and as a neutron source. Th… ▽ More We examine the performance of pure boron, boron carbide, high density carbon, and boron nitride ablators in the polar direct drive exploding pusher (PDXP) platform. The platform uses the polar direct drive configuration at the National Ignition Facility to drive high ion temperatures in a room temperature capsule and has potential applications for plasma physics studies and as a neutron source. The higher tensile strength of these materials compared to plastic enables a thinner ablator to support higher gas pressures, which could help optimize its performance for plasma physics experiments, while ablators containing boron enable the possiblity of collecting addtional data to constrain models of the platform. Applying recently developed and experimentally validated equation of state models for the boron materials, we examine the performance of these materials as ablators in 2D simulations, with particular focus on changes to the ablator and gas areal density, as well as the predicted symmetry of the inherently 2D implosion. △ Less

Submitted 30 December, 2020; v1 submitted 28 June, 2020; originally announced June 2020.

Report number: LLNL-JRNL-803851

Simultaneous compression and opacity data from time-series radiography with a Lagrangian marker

Authors: Damian C. Swift, Andrea L. Kritcher, James A. Hawreliak, James Gaffney, Amy Lazicki, Andrew MacPhee, Benjamin Bachmann, Tilo Doeppner, Joseph Nilsen, Heather D. Whitley, Gilbert W. Collins, Siegfried Glenzer, Stephen D. Rothman, Dominik Kraus, Roger W. Falcone

Abstract: Time-resolved radiography can be used to obtain absolute shock Hugoniot states by simultaneously measuring at least two mechanical parameters of the shock, and this technique is particularly suitable for one-dimensional converging shocks where a single experiment probes a range of pressures as the converging shock strengthens. However, at sufficiently high pressures, the shocked material becomes h… ▽ More Time-resolved radiography can be used to obtain absolute shock Hugoniot states by simultaneously measuring at least two mechanical parameters of the shock, and this technique is particularly suitable for one-dimensional converging shocks where a single experiment probes a range of pressures as the converging shock strengthens. However, at sufficiently high pressures, the shocked material becomes hot enough that the x-ray opacity falls significantly. If the system includes a Lagrangian marker, such that the mass within the marker is known, this additional information can be used to constrain the opacity as well as the Hugoniot state. In the limit that the opacity changes only on shock heating, and not significantly on subsequent isentropic compression, the opacity of shocked material can be determined uniquely. More generally, it is necessary to assume the form of the variation of opacity with isentropic compression, or to introduce multiple marker layers. Alternatively, assuming either the equation of state or the opacity, the presence of a marker layer in such experiments enables the non-assumed property to be deduced more accurately than from the radiographic density reconstruction alone. An example analysis is shown for measurements of a converging shock wave in polystyrene, at the National Ignition Facility. △ Less

Submitted 15 October, 2020; v1 submitted 23 June, 2020; originally announced June 2020.

Report number: LLNL-JRNL-805296

Journal ref: Rev. Sci. Instrum. 92, 063514 (2021)

Zero-energy excitation in the classical kagome antiferromagnet NaBa$_{2}$Mn$_{3}$F$_{11}$

Authors: Shohei Hayashida, Hajime Ishikawa, Yoshihiko Okamoto, Tsuyoshi Okubo, Zenji Hiroi, Gøran J. Nilsen, Hannu Mutka, Takatsugu Masuda

Abstract: We performed inelastic neutron scattering measurements on a polycrystalline sample of a classical kagome antiferromagnet NaBa$_{2}$Mn$_{3}$F$_{11}$ to investigate the possibility of a dispersionless zero-energy excitation associated with rotation of spins along the chains. The observed spectra indeed exhibit such an excitation with strong intensity at low energy, as well as dispersive excitations… ▽ More We performed inelastic neutron scattering measurements on a polycrystalline sample of a classical kagome antiferromagnet NaBa$_{2}$Mn$_{3}$F$_{11}$ to investigate the possibility of a dispersionless zero-energy excitation associated with rotation of spins along the chains. The observed spectra indeed exhibit such an excitation with strong intensity at low energy, as well as dispersive excitations with weak intensity at high energy. Combining the measurements with calculations from linear spin-wave theory reveals that NaBa$_{2}$Mn$_{3}$F$_{11}$ is a good realization of the classical kagome antiferromagnet which exhibits a dispersionless mode lifted by the magnetic dipole-dipole interaction. △ Less

Submitted 4 June, 2020; v1 submitted 2 June, 2020; originally announced June 2020.

Journal ref: Phys. Rev. B 101, 214409 (2020)

Large Easy Axis Anisotropy in the One-Dimensional Magnet BaMo(PO$_4$)$_2$

Authors: Aly H. Abdeldaim, Danis I. Badrtdinov, Alexandra S. Gibbs, Pascal Manuel, Helen C. Walker, Manh Duc Le, Chien Hung Wu, Dariusz Wardecki, Sten-Gunnar Eriksson, Yaroslav O. Kvashnin, Alexander A. Tsirlin, Gøran J. Nilsen

Abstract: We present an extensive experimental and theoretical study on the low-temperature magnetic properties of the monoclinic anhydrous alum compound BaMo(PO$_4$)$_2$. The magnetic susceptibility reveals strong antiferromagnetic interactions $θ_{CW} = -167$~K and long-range magnetic order at $T_N=22$~K, in agreement with a recent report. Powder neutron diffraction furthermore shows that the order is col… ▽ More We present an extensive experimental and theoretical study on the low-temperature magnetic properties of the monoclinic anhydrous alum compound BaMo(PO$_4$)$_2$. The magnetic susceptibility reveals strong antiferromagnetic interactions $θ_{CW} = -167$~K and long-range magnetic order at $T_N=22$~K, in agreement with a recent report. Powder neutron diffraction furthermore shows that the order is collinear, with the moments near the $ac$ plane. Neutron spectroscopy reveals a large excitation gap $Δ= 15$~meV in the low-temperature ordered phase, suggesting a much larger easy-axis spin anisotropy than anticipated. However, the large anisotropy justifies the relatively high ordered moment, Néel temperature, and collinear order observed experimentally, and is furthermore reproduced in a first principles calculations using a new computational scheme. We therefore propose BaMo(PO$_4$)$_2$ to host $S=1$ antiferromagnetic chains with large easy-axis anisotropy, which has been theoretically predicted to realize novel excitation continua. △ Less

Submitted 9 December, 2019; originally announced December 2019.

Comments: 11 pages, 7 figures

Temperature dependence of the $(π,0)$ anomaly in the excitation spectrum of the 2D quantum Heisenberg antiferromagnet

Authors: W. Wan, N. B. Christensen, A. W. Sandvik, P. Tregenna-Piggott, G. J. Nilsen, M. Mourigal, T. G. Perring, C. D. Frost, D. F. McMorrow, H. M. Rønnow

Abstract: It is well established that in the low-temperature limit, the two-dimensional quantum Heisenberg antiferromagnet on a square lattice (2DQHAFSL) exhibits an anomaly in its spectrum at short-wavelengths on the zone-boundary. In the vicinity of the $(π,0)$ point the pole in the one-magnon response exhibits a downward dispersion, is heavily damped and attenuated, giving way to an isotropic continuum o… ▽ More It is well established that in the low-temperature limit, the two-dimensional quantum Heisenberg antiferromagnet on a square lattice (2DQHAFSL) exhibits an anomaly in its spectrum at short-wavelengths on the zone-boundary. In the vicinity of the $(π,0)$ point the pole in the one-magnon response exhibits a downward dispersion, is heavily damped and attenuated, giving way to an isotropic continuum of excitations extending to high energies. The origin of the anomaly and the presence of the continuum are of current theoretical interest, with suggestions focused around the idea that the latter evidences the existence of spinons in a two-dimensional system. Here we present the results of neutron inelastic scattering experiments and Quantum Monte Carlo calculations on the metallo-organic compound Cu(DCOO)$_2\cdot 4$D$_2$O (CFTD), an excellent physical realisation of the 2DQHAFSL, designed to investigate how the anomaly at $(π,0)$ evolves up to finite temperatures $T/J\sim2/3$. Our data reveal that on warming the anomaly survives the loss of long-range, three-dimensional order, and is thus a robust feature of the two-dimensional system. With further increase of temperature the zone-boundary response gradually softens and broadens, washing out the $(π,0)$ anomaly. This is confirmed by a comparison of our data with the results of finite-temperature Quantum Monte Carlo simulations where the two are found to be in good accord. At lower energies, in the vicinity of the antiferromagnetic zone centre, there was no significant softening of the magnetic excitations over the range of temperatures investigated. △ Less

Submitted 3 December, 2019; originally announced December 2019.

Comments: Dedicated to the life and work of Professor Roger Cowley. 22 pages, 8 figures

A quantum liquid of magnetic octupoles on the pyrochlore lattice

Authors: Romain Sibille, Nicolas Gauthier, Elsa Lhotel, Victor Porée, Vladimir Pomjakushin, Russell A. Ewings, Toby G. Perring, Jacques Ollivier, Andrew Wildes, Clemens Ritter, Thomas C. Hansen, David A. Keen, Gøran J. Nilsen, Lukas Keller, Sylvain Petit, Tom Fennell

Abstract: Spin liquids are highly correlated yet disordered states formed by the entanglement of magnetic dipoles$^1$. Theories typically define such states using gauge fields and deconfined quasiparticle excitations that emerge from a simple rule governing the local ground state of a frustrated magnet. For example, the '2-in-2-out' ice rule for dipole moments on a tetrahedron can lead to a quantum spin ice… ▽ More Spin liquids are highly correlated yet disordered states formed by the entanglement of magnetic dipoles$^1$. Theories typically define such states using gauge fields and deconfined quasiparticle excitations that emerge from a simple rule governing the local ground state of a frustrated magnet. For example, the '2-in-2-out' ice rule for dipole moments on a tetrahedron can lead to a quantum spin ice in rare-earth pyrochlores - a state described by a lattice gauge theory of quantum electrodynamics$^{2-4}$. However, f-electron ions often carry multipole degrees of freedom of higher rank than dipoles, leading to intriguing behaviours and 'hidden' orders$^{5-6}$. Here we show that the correlated ground state of a Ce$^{3+}$-based pyrochlore, Ce$_2$Sn$_2$O$_7$, is a quantum liquid of magnetic octupoles. Our neutron scattering results are consistent with the formation of a fluid-like state of matter, but the intensity distribution is weighted to larger scattering vectors, which indicates that the correlated degrees of freedom have a more complex magnetization density than that typical of magnetic dipoles in a spin liquid. The temperature evolution of the bulk properties in the correlated regime below 1 Kelvin is well reproduced using a model of dipole-octupole doublets on a pyrochlore lattice$^{7-8}$. The nature and strength of the octupole-octupole couplings, together with the existence of a continuum of excitations attributed to spinons, provides further evidence for a quantum ice of octupoles governed by a '2-plus-2-minus' rule. Our work identifies Ce$_2$Sn$_2$O$_7$ as a unique example of a material where frustrated multipoles form a 'hidden' topological order, thus generalizing observations on quantum spin liquids to multipolar phases that can support novel types of emergent fields and excitations. △ Less

Submitted 2 December, 2019; originally announced December 2019.

Journal ref: Nature Physics (2020)

Magnetic order and single-ion anisotropy in Tb$_3$Ga$_5$O$_{12}$

Authors: R. Wawrzyńczak, B. Tomasello, P. Manuel, D. Khalyavin, M. D. Le, T. Guidi, A. Cervellino, T. Ziman, M. Boehm, G. J. Nilsen, T. Fennell

Abstract: Terbium gallium garnet (TGG), Tb$_3$Ga$_5$O$_{12}$, is well known for its applications in laser optics, but also exhibits complex low-temperature magnetism that is not yet fully understood. Its low-temperature magnetic order is determined by means of time-of-flight neutron powder diffraction. It is found to be a multiaxial antiferromagnet with magnetic Tb$^{3+}$ ions forming six sublattices of mag… ▽ More Terbium gallium garnet (TGG), Tb$_3$Ga$_5$O$_{12}$, is well known for its applications in laser optics, but also exhibits complex low-temperature magnetism that is not yet fully understood. Its low-temperature magnetic order is determined by means of time-of-flight neutron powder diffraction. It is found to be a multiaxial antiferromagnet with magnetic Tb$^{3+}$ ions forming six sublattices of magnetic moments aligned parallel and anti-parallel to the $\langle100\rangle$ crystallographic directions of the cubic unit cell. The structure displays strong easy-axis anisotropy with respect to a two-fold axis of symmetry in the local orthorhombic environment of the Tb$^{3+}$ sites. The crystal-field splitting within the single-ion ground-state manifold is investigated by inelastic neutron scattering on powder samples. A strong temperature dependence of the quasidoublet ground-state is observed and revised parameters of the crystal-field Hamiltonian are given. The results of bulk magnetic susceptibility and magnetisation measurements are in good agreement with values based on the crystal-field model down to 20~K, where the onset of magnetic correlations is observed. △ Less

Submitted 18 October, 2019; originally announced October 2019.

Comments: 15 pages, 12 figures

Journal ref: Phys. Rev. B 100, 094442 (2019)

Putative spin-nematic phase in BaCdVO(PO$_{4}$)$_{2}$

Authors: M. Skoulatos, F. Rucker, G. J. Nilsen, A. Bertin, E. Pomjakushina, J. Ollivier, A. Schneidewind, R. Georgii, O. Zaharko, L. Keller, Ch. Rüegg, C. Pfleiderer, B. Schmidt, N. Shannon, A. Kriele, A. Senyshyn, A. Smerald

Abstract: We report neutron scattering and AC magnetic susceptibility measurements of the 2D spin-1/2 frustrated magnet BaCdVO(PO$_{4}$)$_{2}$. At temperatures well below $T_{\sf N}\approx 1K$, we show that only 34 % of the spin moment orders in an up-up-down-down strip structure. Dominant magnetic diffuse scattering and comparison to published $μ$sr measurements indicates that the remaining 66 % is fluctua… ▽ More We report neutron scattering and AC magnetic susceptibility measurements of the 2D spin-1/2 frustrated magnet BaCdVO(PO$_{4}$)$_{2}$. At temperatures well below $T_{\sf N}\approx 1K$, we show that only 34 % of the spin moment orders in an up-up-down-down strip structure. Dominant magnetic diffuse scattering and comparison to published $μ$sr measurements indicates that the remaining 66 % is fluctuating. This demonstrates the presence of strong frustration, associated with competing ferromagnetic and antiferromagnetic interactions, and points to a subtle ordering mechanism driven by magnon interactions. On applying magnetic field, we find that at $T=0.1$ K the magnetic order vanishes at 3.8 T, whereas magnetic saturation is reached only above 4.5 T. We argue that the putative high-field phase is a realization of the long-sought bond-spin-nematic state. △ Less

Submitted 29 March, 2019; originally announced March 2019.

Journal ref: Phys. Rev. B 100, 014405 (2019)

Role of defects in determining the magnetic ground state of ytterbium titanate

Authors: D. F. Bowman, E. Cemal, T. Lehner, A. R. Wildes, L. Mangin-Thro, G. J. Nilsen, M. J. Gutmann, D. J. Voneshen, D. Prabhakaran, A. T. Boothroyd, D. G. Porter, C. Castelnovo, K. Refson, J. P. Goff

Abstract: Pyrochlore systems are ideally suited to the exploration of geometrical frustration in three dimensions, and their rich phenomenology encompasses topological order and fractional excitations. Classical spin ices provide the first context in which it is possible to control emergent magnetic monopoles, and anisotropic exchange leads to even richer behaviour associated with large quantum fluctuations… ▽ More Pyrochlore systems are ideally suited to the exploration of geometrical frustration in three dimensions, and their rich phenomenology encompasses topological order and fractional excitations. Classical spin ices provide the first context in which it is possible to control emergent magnetic monopoles, and anisotropic exchange leads to even richer behaviour associated with large quantum fluctuations. Whether the magnetic ground state of Yb2Ti2O7 is a quantum spin liquid or a ferromagnetic phase induced by a Higgs transition appears to be sample dependent. Here we have determined the role of structural defects on the magnetic ground state via the diffuse scattering of neutrons. We find that oxygen vacancies stabilise the spin liquid phase and the stuffing of Ti sites by Yb suppresses it. Samples in which the oxygen vacancies have been eliminated by annealing in oxygen exhibit a transition to a ferromagnetic phase, and this is the true magnetic ground state. △ Less

Submitted 19 February, 2019; originally announced February 2019.

Comments: 15 pages, 1 table, 4 figures, journal article

Journal ref: Nature Communications 10:637 (2019)

Equation of state of warm-dense boron nitride combining computation, modeling, and experiment

Authors: Shuai Zhang, Amy Lazicki, Burkhard Militzer, Lin H. Yang, Kyle Caspersen, Jim A. Gaffney, Markus W. Däne, John E. Pask, Walter R. Johnson, Abhiraj Sharma, Phanish Suryanarayana, Duane D. Johnson, Andrey V. Smirnov, Philip A. Sterne, David Erskine, Richard A. London, Federica Coppari, Damian Swift, Joseph Nilsen, Art J. Nelson, Heather D. Whitley

Abstract: The equation of state (EOS) of materials at warm dense conditions poses significant challenges to both theory and experiment. We report a combined computational, modeling, and experimental investigation leveraging new theoretical and experimental capabilities to investigate warm-dense boron nitride (BN). The simulation methodologies include path integral Monte Carlo (PIMC), several density functio… ▽ More The equation of state (EOS) of materials at warm dense conditions poses significant challenges to both theory and experiment. We report a combined computational, modeling, and experimental investigation leveraging new theoretical and experimental capabilities to investigate warm-dense boron nitride (BN). The simulation methodologies include path integral Monte Carlo (PIMC), several density functional theory (DFT) molecular dynamics methods [plane-wave pseudopotential, Fermi operator expansion (FOE), and spectral quadrature (SQ)], activity expansion (ACTEX), and all-electron Green's function Korringa-Kohn-Rostoker (MECCA), and compute the pressure and internal energy of BN over a broad range of densities ($ρ$) and temperatures ($T$). Our experiments were conducted at the Omega laser facility and measured the Hugoniot of BN to unprecedented pressures (12--30 Mbar). The EOSs computed using different methods cross validate one another, and the experimental Hugoniot are in good agreement with our theoretical predictions. We assess that the largest discrepancies between theoretical predictions are $<$4% in pressure and $<$3% in energy and occur at $10^6$ K. We find remarkable consistency between the EOS from DFT calculations performed on different platforms and using different exchange-correlation functionals and those from PIMC using free-particle nodes. This provides strong evidence for the accuracy of both PIMC and DFT in the warm-dense regime. Moreover, SQ and FOE data have significantly smaller error bars than PIMC, and so represent significant advances for efficient computation at high $T$. We also construct tabular EOS models and clarify the ionic and electronic structure of BN over a broad $T-ρ$ range and quantify their roles in the EOS. The tabular models may be utilized for future simulations of laser-driven experiments that include BN as a candidate ablator material. △ Less

Submitted 2 February, 2019; originally announced February 2019.

Comments: 19 pages, 14 figures, 4 tables

Journal ref: Phys. Rev. B 99, 165103 (2019)

Nuclear and Magnetic Structures of the Frustrated Quantum Antiferromagnet Barlowite, Cu$_{4}$(OH)$_{6}$FBr

Authors: K. Tustain, G. J. Nilsen, C. Ritter, I. da Silva, L. Clark

Abstract: Barlowite, Cu$_{4}$(OH)$_{6}$FBr, has attracted much attention as the parent compound of a new series of quantum spin liquid candidates, Zn$_{x}$Cu$_{4-x}$(OH)$_{6}$FBr. While it is known to undergo a magnetic phase transition to a long-range ordered state at $T_{N} = 15$ K, there is still no consensus over either its nuclear or magnetic structures. Here, we use comprehensive powder neutron diffra… ▽ More Barlowite, Cu$_{4}$(OH)$_{6}$FBr, has attracted much attention as the parent compound of a new series of quantum spin liquid candidates, Zn$_{x}$Cu$_{4-x}$(OH)$_{6}$FBr. While it is known to undergo a magnetic phase transition to a long-range ordered state at $T_{N} = 15$ K, there is still no consensus over either its nuclear or magnetic structures. Here, we use comprehensive powder neutron diffraction studies on deuterated samples of barlowite to demonstrate that the only space group consistent with the observed nuclear and magnetic diffraction at low-temperatures is the orthorhombic $Pnma$ space group. We furthermore conclude that the magnetic intensity at $T < T_{N}$ is correctly described by the $Pn^\prime m^\prime a$ magnetic space group, which crucially allows the ferromagnetic component observed in previous single-crystal and powder magnetisation measurements. As such, the magnetic structure of barlowite resembles that of the related material clinoatacamite, Cu$_{4}$(OH)$_{6}$Cl$_{2}$, the parent compound of the well-known quantum spin liquid candidate hebertsmithite, ZnCu$_{3}$(OH)$_{6}$Cl$_{2}$. △ Less

Submitted 3 October, 2018; originally announced October 2018.

Comments: 11 pages, 4 figures, supplemental material

Journal ref: Phys. Rev. Materials 2, 111405 (2018)

Opacity of Shock-Heated Boron Plasmas

Authors: W. R. Johnson, J. Nilsen

Abstract: Standard measures of opacity, the imaginary part of the atomic scattering factor $f_2$ and the x-ray mass attenuation coefficient $μ/ρ$, are evaluated in shock-heated boron, boron carbide and boron nitride plasmas. The Hugoniot equation, relating the temperature $T$ behind a shock wave to the compression ratio $ρ/ρ_0$ across the shock front, is used in connection with the plasma equation of state… ▽ More Standard measures of opacity, the imaginary part of the atomic scattering factor $f_2$ and the x-ray mass attenuation coefficient $μ/ρ$, are evaluated in shock-heated boron, boron carbide and boron nitride plasmas. The Hugoniot equation, relating the temperature $T$ behind a shock wave to the compression ratio $ρ/ρ_0$ across the shock front, is used in connection with the plasma equation of state to determine the pressure $p$, effective plasma charge $Z^\ast$ and the K-shell occupation in terms of $ρ/ρ_0$. Solutions of the Hugoniot equation (determined within the framework of the generalized Thomas-Fermi theory) reveal that the K-shell occupation in low-Z ions decreases rapidly from 2 to 0.1 as the temperature increases from 20eV to 500eV; a temperature range in which the shock compression ratio is near 4. The average-atom model (a quantum mechanical version of the generalized Thomas-Fermi theory) is used to determine K-shell and continuum wave functions and the photoionization cross section for x-rays in the energy range $ω=1$eV to 10keV, where the opacity is dominated by the atomic photoionization process. For an uncompressed boron plasma at $T=10$eV, where the K-shell is filled, the average-atom cross section, the atomic scattering factor and the mass attenuation coefficient are all shown to agree closely with previous (cold matter) tabulations. For shock-compressed plasmas, the dependence of the $μ/ρ$ on temperature can be approximated by scaling previously tabulated cold-matter values by the relative K-shell occupation, however, there is a relatively small residual dependence arising from the photoionization cross section. Attenuation coefficients $μ$ for a 9 keV x-ray are given as functions of $T$ along the Hugoniot for B, C, B$_4$C and BN plasmas. △ Less

Submitted 15 March, 2019; v1 submitted 29 September, 2018; originally announced October 2018.

Comments: 15 pages, 8 figures, prepared as a HEDP preprint to accompany a presentation at the RPHDM18 conference

Report number: LLNL-JRNL-758940-DRAFT

Theoretical and experimental investigation of the equation of state of boron plasmas

Authors: Shuai Zhang, Burkhard Militzer, Michelle C. Gregor, Kyle Caspersen, Lin H. Yang, Tadashi Ogitsu, Damian Swift, Amy Lazicki, D. Erskine, Richard A. London, P. M. Celliers, Joseph Nilsen, Philip A. Sterne, Heather D. Whitley

Abstract: We report a theoretical equation of state (EOS) table for boron across a wide range of temperatures (5.1$\times$10$^4$-5.2$\times$10$^8$ K) and densities (0.25-49 g/cm$^3$), and experimental shock Hugoniot data at unprecedented high pressures (5608$\pm$118 GPa). The calculations are performed with full, first-principles methods combining path integral Monte Carlo (PIMC) at high temperatures and de… ▽ More We report a theoretical equation of state (EOS) table for boron across a wide range of temperatures (5.1$\times$10$^4$-5.2$\times$10$^8$ K) and densities (0.25-49 g/cm$^3$), and experimental shock Hugoniot data at unprecedented high pressures (5608$\pm$118 GPa). The calculations are performed with full, first-principles methods combining path integral Monte Carlo (PIMC) at high temperatures and density functional theory molecular dynamics (DFT-MD) methods at lower temperatures. PIMC and DFT-MD cross-validate each other by providing coherent EOS (difference $<$1.5 Hartree/boron in energy and $<$5% in pressure) at 5.1$\times$10$^5$ K. The Hugoniot measurement is conducted at the National Ignition Facility using a planar shock platform. The pressure-density relation found in our shock experiment is on top of the shock Hugoniot profile predicted with our first-principles EOS and a semi-empirical EOS table (LEOS 50). We investigate the self diffusivity and the effect of thermal and pressure-driven ionization on the EOS and shock compression behavior in high pressure and temperature conditions We study the performance sensitivity of a polar direct-drive exploding pusher platform to pressure variations based on comparison of the first-principles calculations with LEOS 50 via 1D hydrodynamic simulations. The results are valuable for future theoretical and experimental studies and engineering design in high energy density research. (LLNL-JRNL-748227) △ Less

Submitted 30 April, 2018; originally announced April 2018.

Comments: 12 pages, 9 figures, 2 tables

Journal ref: Phys. Rev. E 98, 023205 (2018)

X-ray line coincidence photopumping in a solar flare

Authors: F. P. Keenan, K. Poppenhaeger, M. Mathioudakis, S. J. Rose, J. Flowerdew, D. Hynes, D. J. Christian, J. Nilsen, W. R. Johnson

Abstract: Line coincidence photopumping is a process where the electrons of an atomic or molecular species are radiatively excited through the absorption of line emission from another species at a coincident wavelength. There are many instances of line coincidence photopumping in astrophysical sources at optical and ultraviolet wavelengths, with the most famous example being Bowen fluorescence (pumping of O… ▽ More Line coincidence photopumping is a process where the electrons of an atomic or molecular species are radiatively excited through the absorption of line emission from another species at a coincident wavelength. There are many instances of line coincidence photopumping in astrophysical sources at optical and ultraviolet wavelengths, with the most famous example being Bowen fluorescence (pumping of O III 303.80 A by He II), but none to our knowledge in X-rays. However, here we report on a scheme where a He-like line of Ne IX at 11.000 A is photopumped by He-like Na X at 11.003 A, which predicts significant intensity enhancement in the Ne IX 82.76 A transition under physical conditions found in solar flare plasmas. A comparison of our theoretical models with published X-ray observations of a solar flare obtained during a rocket flight provides evidence for line enhancement, with the measured degree of enhancement being consistent with that expected from theory, a truly surprising result. Observations of this enhancement during flares on stars other than the Sun would provide a powerful new diagnostic tool for determining the sizes of flare loops in these distant, spatially-unresolved, astronomical sources. △ Less

Submitted 21 November, 2017; originally announced November 2017.

Comments: 5 pages, 2 figures, MNRAS in press

Single-mode Near-infrared Lasing in a GaAsSb/GaAs Nanowire Superlattice at Room Temperature

Authors: Dingding Ren, Lyubomir Ahtapodov, Julie S. Nilsen, Jianfeng Yang, Anders Gustafsson, Junghwan Huh, Gavin J. Conibeer, Antonius T. J. van Helvoort, Bjørn-Ove Fimland, Helge Weman

Abstract: Semiconductor nanowire lasers can produce guided coherent light emission with miniaturized geometry, bringing about new possibility for a variety of applications including nanophotonic circuits, optical sensing, and on-chip and chip-to-chip optical communications. Here, we report on the realization of single-mode room-temperature lasing from 890 nm to 990 nm utilizing a novel design of single nano… ▽ More Semiconductor nanowire lasers can produce guided coherent light emission with miniaturized geometry, bringing about new possibility for a variety of applications including nanophotonic circuits, optical sensing, and on-chip and chip-to-chip optical communications. Here, we report on the realization of single-mode room-temperature lasing from 890 nm to 990 nm utilizing a novel design of single nanowires with GaAsSb-based multiple superlattices as gain medium under optical pumping. The wavelength tunability with comprehensively enhanced lasing performance is shown to result from the unique nanowire structure with efficient gain materials, which delivers a lasing quality factor as high as 1250, a reduced lasing threshold ~ 6 kW cm-2 and a high characteristic temperature ~ 129 K. These results present a major advancement for the design and synthesis of nanowire laser structures, which can pave the way towards future nanoscale integrated optoelectronic systems with stunning performance. △ Less

Submitted 23 August, 2017; originally announced August 2017.

Real-space investigation of short-range magnetic correlations in fluoride pyrochlores NaCaCo$_2$F$_7$ and NaSrCo$_2$F$_7$ with magnetic pair distribution function analysis

Authors: Benjamin A. Frandsen, Kate A. Ross, Jason W. Krizan, Gøran J. Nilsen, Andrew R. Wildes, Robert J. Cava, Robert J. Birgeneau, Simon J. L. Billinge

Abstract: We present time-of-flight neutron total scattering and polarized neutron scattering measurements of the magnetically frustrated compounds NaCaCo$_2$F$_7$ and NaSrCo$_2$F$_7$, which belong to a class of recently discovered pyrochlore compounds based on transition metals and fluorine. The magnetic pair distribution function (mPDF) technique is used to analyze and model the total scattering data in r… ▽ More We present time-of-flight neutron total scattering and polarized neutron scattering measurements of the magnetically frustrated compounds NaCaCo$_2$F$_7$ and NaSrCo$_2$F$_7$, which belong to a class of recently discovered pyrochlore compounds based on transition metals and fluorine. The magnetic pair distribution function (mPDF) technique is used to analyze and model the total scattering data in real space. We find that a previously-proposed model of short-range XY-like correlations with a length scale of 10-15 Å, combined with nearest-neighbor collinear antiferromagnetic correlations, accurately describes the mPDF data at low temperature, confirming the magnetic ground state in these materials. This model is further verified by the polarized neutron scattering data. From an analysis of the temperature dependence of the mPDF and polarized neutron scattering data, we find that short-range correlations persist on the nearest-neighbor length scale up to 200 K, approximately two orders of magnitude higher than the spin freezing temperatures of these compounds. These results highlight the opportunity presented by these new pyrochlore compounds to study the effects of geometric frustration at relatively high temperatures, while also advancing the mPDF technique and providing a novel opportunity to investigate a genuinely short-range-ordered magnetic ground state directly in real space. △ Less

Submitted 30 December, 2017; v1 submitted 18 August, 2017; originally announced August 2017.

Journal ref: Phys. Rev. Materials 1, 074412 (2017)

Classical Spin Nematic Transition in LiGa$_{0.95}$In$_{0.05}$Cr$_4$O$_8$

Authors: R. Wawrzyńczak, Y. Tanaka, M. Yoshida, Y. Okamoto, P. Manuel, N. Casati, Z. Hiroi, M. Takigawa, G. J. Nilsen

Abstract: We present the results of a combined $^7$Li NMR and diffraction study on LiGa$_{0.95}$In$_{0.05}$Cr$_4$O$_8$, a member of LiGa$_{1-x}$In$_{x}$Cr$_4$O$_8$ "breathing" pyrochlore family. Via specific heat and NMR measurements, we find that the complex sequence of first-order transitions observed for $x=0$ is replaced by a single, apparently second-order transition at $T_f=11$ K. Neutron and X-ray di… ▽ More We present the results of a combined $^7$Li NMR and diffraction study on LiGa$_{0.95}$In$_{0.05}$Cr$_4$O$_8$, a member of LiGa$_{1-x}$In$_{x}$Cr$_4$O$_8$ "breathing" pyrochlore family. Via specific heat and NMR measurements, we find that the complex sequence of first-order transitions observed for $x=0$ is replaced by a single, apparently second-order transition at $T_f=11$ K. Neutron and X-ray diffraction rule out both structural symmetry lowering and magnetic long-range order as the origin of this transition. Instead, reverse Monte Carlo fitting of the magnetic diffuse scattering indicates that the low temperature phase may be described as a collinear spin nematic state, characterized by a quadrupolar order parameter. This state also shows signs of short range order between collinear spin arrangements on tetrahedra, revealed by mapping the reverse Monte Carlo spin configurations onto a three-state color model describing the manifold of nematic states. △ Less

Submitted 26 June, 2017; originally announced June 2017.

Comments: 11 pages, 10 figures

Phase diagram of multiferroic KCu$_3$As$_2$O$_7$(OD)$_3$

Authors: Gøran J. Nilsen, Virginie Simonet, Claire V. Colin, Ryutaro Okuma, Yoshihiko Okamoto, Masashi Tokunaga, Thomas C. Hansen, Dmitry D. Khalyavin, Zenji Hiroi

Abstract: The layered compound KCu$_3$As$_2$O$_7$(OD)$_3$, comprising distorted kagome planes of $S=1/2$ Cu$^{2+}$ ions, is a recent addition to the family of type-II multiferroics. Previous zero field neutron diffraction work has found two helically ordered regimes in \kns, each showing a distinct coupling between the magnetic and ferroelectric order parameters. Here, we extend this work to magnetic fields… ▽ More The layered compound KCu$_3$As$_2$O$_7$(OD)$_3$, comprising distorted kagome planes of $S=1/2$ Cu$^{2+}$ ions, is a recent addition to the family of type-II multiferroics. Previous zero field neutron diffraction work has found two helically ordered regimes in \kns, each showing a distinct coupling between the magnetic and ferroelectric order parameters. Here, we extend this work to magnetic fields up to $20$~T using neutron powder diffraction, capacitance, polarization, and high-field magnetization measurements, hence determining the $H-T$ phase diagram. We find metamagnetic transitions in both low temperatures phases around $μ_0 H_c \sim 3.7$~T, which neutron powder diffraction reveals to correspond to a rotation of the helix plane away from the easy plane, as well as a small change in the propagation vector. Furthermore, we show that the sign of the ferroelectric polarization is reversible in a magnetic field, although no change is observed (or expected on the basis of the magnetic structure) due to the transition at $3.7$~T. We finally justify the temperature dependence of the polarization in both zero-field ordered phases by a symmetry analysis of the free energy expansion. △ Less

Submitted 10 February, 2017; originally announced February 2017.

Journal ref: Phys. Rev. B 95, 214415 (2017)

From 3D static order to 2D dynamic correlations in the tetragonal double perovskite Sr2CuWO6

Authors: Oliver J. Burrows, Gøran J. Nilsen, Emmanuelle Suard, Mark Telling, J. Ross Stewart, Mark A. de Vries

Abstract: \SCWO is a tetragonal ($I4/m$) double perovskite with diamagnetic W$^{6+}$ (5d$^0$) and (Jahn-Teller-active) Cu$^{2+}$ (3$d^9$) in a rocksalt-ordered arrangement on the B sites, modelling an anisotropic fcc-lattice antiferromagnet with $S=1/2$. We have studied the magnetic structure, dynamics and thermodynamic properties of this material above and below the \Neel{}-ordering transition at… ▽ More \SCWO is a tetragonal ($I4/m$) double perovskite with diamagnetic W$^{6+}$ (5d$^0$) and (Jahn-Teller-active) Cu$^{2+}$ (3$d^9$) in a rocksalt-ordered arrangement on the B sites, modelling an anisotropic fcc-lattice antiferromagnet with $S=1/2$. We have studied the magnetic structure, dynamics and thermodynamic properties of this material above and below the \Neel{}-ordering transition at $T_{\textrm{N}} = 24$~K. Muon spin relaxation spectroscopy, neutron diffraction and neutron spectroscopy experiments show that the transition at $T_{\textrm{N}}$ on heating is from static 3D type-2 antiferromagnetism to 2D dynamic correlations that remain detectable up to 100~K. Above $T_{\textrm{N}}$ the muon relaxation is described by a compressed exponential decay, revealing a strong correlation between the spin fluctuation rate and the 2D correlation length. The low-temperature muon data and the entropy release around $T_{\textrm{N}}$ as obtained from the heat capacity data point to small magnetic domains, probably of the order of 60~Å. We note that the material behaves more like a quasi-2D square lattice antiferromagnet than could be expected considering its structure only. A possible explanation might lie in the unusual strong spin-orbit-coupled magnetism of the W 5$d$ levels, as the dominant exchange interaction is via the Cu-O-W-O-Cu pathway. △ Less

Submitted 12 February, 2016; originally announced February 2016.

Comments: Just over 9 pages, 12 figures. Submitted

Average-atom treatment of relaxation time in X-ray Thomson scattering from warm-dense matter

Authors: Walter R. Johnson, Joseph Nilsen

Abstract: The influence of finite relaxation times on Thomson scattering from warm-dense plasmas is examined within the framework of the average-atom approximation. Presently most calculations use the collision-free Lindhard dielectric function to evaluate the free-electron contribution to the Thomson cross section. In this work, we use the Mermin dielectric function, which includes relaxation time explicit… ▽ More The influence of finite relaxation times on Thomson scattering from warm-dense plasmas is examined within the framework of the average-atom approximation. Presently most calculations use the collision-free Lindhard dielectric function to evaluate the free-electron contribution to the Thomson cross section. In this work, we use the Mermin dielectric function, which includes relaxation time explicitly. The relaxation time is evaluated by treating the average atom as an impurity in a uniform electron gas and depends critically on the transport cross section. The calculated relaxation rates agree well with values inferred from the Ziman formula for the static conductivity and also with rates inferred from a fit to the frequency-dependent conductivity. Transport cross sections determined by the phase-shift analysis in the average-atom potential are compared with those evaluated in the commonly used Born approximation. The Born approximation converges to the exact cross sections at high energies; however, differences that occur at low energies lead to corresponding differences in relaxation rates. The relative importance of including relaxation time when modeling X-ray Thomson scattering spectra is examined by comparing calculations of the free-electron dynamic structure function for Thomson scattering using Lindhard and Mermin dielectric functions. Applications are given to warm-dense Be plasmas, with temperatures ranging from 2 to 32~eV and densities ranging from 2 to 64~g/cc. △ Less

Submitted 1 April, 2016; v1 submitted 30 December, 2015; originally announced December 2015.

Comments: 7 pages, 6 figures, 2 tables

Journal ref: Phys. Rev. E 93, 033205 (2016)

Magnetic Phase Diagram of the Breathing Pyrochlore Antiferromagnet LiGa1-xInxCr4O8

Authors: Yoshihiko Okamoto, Gøran J. Nilsen, Taishi Nakazono, Zenji Hiroi

Abstract: The spinel oxides LiGaCr4O8 and LiInCr4O8 contain size-alternating pyrochlore lattices of spin-3/2 Cr3+ tetrahedra with different magnitudes of alternation. We show here that the solid solutions LiGa1-xInxCr4O8 between these two 'breathing' pyrochlore compounds display (i) rapid suppression of magnetic and structural transitions upon doping the end members, (ii) spin-glass-like freezing above 2 K… ▽ More The spinel oxides LiGaCr4O8 and LiInCr4O8 contain size-alternating pyrochlore lattices of spin-3/2 Cr3+ tetrahedra with different magnitudes of alternation. We show here that the solid solutions LiGa1-xInxCr4O8 between these two 'breathing' pyrochlore compounds display (i) rapid suppression of magnetic and structural transitions upon doping the end members, (ii) spin-glass-like freezing above 2 K in the range 0.1 $\lesssim$ x $\lesssim$ 0.6, and (iii) apparent spin-gap behavior for x $\gtrsim$ 0.7. Furthermore, no transitions are observed above 2 K at x ~ 0.9, where magnetic susceptibility remains finite at 2 K and magnetic heat capacity shows a quadratic temperature dependence at 1-5 K. Our work shows that breathing pyrochlore compounds provide a unique opportunity for studying both geometrical frustration and bond alternation. △ Less

Submitted 30 March, 2015; originally announced March 2015.

Comments: 5 pages, 4 figures

Journal ref: Journal of the Physical Society of Japan, 84, 043707 (2015)

Fractional excitations in the square-lattice quantum antiferromagnet

Authors: B. Dalla Piazza, M. Mourigal, N. B. Christensen, G. J. Nilsen, P. Tregenna-Piggott, T. G. Perring, M. Enderle, D. F. McMorrow, D. A. Ivanov, H. M. Rønnow

Abstract: The square-lattice quantum Heisenberg antiferromagnet displays a pronounced anomaly of unknown origin in its magnetic excitation spectrum. The anomaly manifests itself only for short wavelength excitations propagating along the direction connecting nearest neighbors. Using polarized neutron spectroscopy, we have fully characterized the magnetic fluctuations in the model metal-organic compound CFTD… ▽ More The square-lattice quantum Heisenberg antiferromagnet displays a pronounced anomaly of unknown origin in its magnetic excitation spectrum. The anomaly manifests itself only for short wavelength excitations propagating along the direction connecting nearest neighbors. Using polarized neutron spectroscopy, we have fully characterized the magnetic fluctuations in the model metal-organic compound CFTD, revealing an isotropic continuum at the anomaly indicative of fractional excitations. A theoretical framework based on the Gutzwiller projection method is developed to explain the origin of the continuum at the anomaly. This indicates that the anomaly arises from deconfined fractional spin-1/2 quasiparticle pairs, the 2D analog of 1D spinons. Away from the anomaly the conventional spin-wave spectrum is recovered as pairs of fractional quasiparticles bind to form spin-1 magnons. Our results therefore establish the existence of fractional quasiparticles in the simplest model two dimensional antiferromagnet even in the absence of frustration. △ Less

Submitted 8 January, 2015; originally announced January 2015.

Comments: 41 pages, 10 figures

Journal ref: Nature Physics 11 (2015) 62-68

From Spin Glass to Quantum Spin Liquid Ground States in Molybdate Pyrochlores

Authors: L. Clark, G. J. Nilsen, E. Kermarrec, G. Ehlers, K. S. Knight, A. Harrison, J. P. Attfield, B. D. Gaulin

Abstract: We present new magnetic heat capacity and neutron scattering results for two magnetically frustrated molybdate pyrochlores: $S=1$ oxide Lu$_2$Mo$_2$O$_7$ and $S={\frac{1}{2}}$ oxynitride Lu$_2$Mo$_2$O$_5$N$_2$. Lu$_2$Mo$_2$O$_7$ undergoes a transition to an unconventional spin glass ground state at $T_f {\sim} 16$ K. However, the preparation of the corresponding oxynitride tunes the nature of the… ▽ More We present new magnetic heat capacity and neutron scattering results for two magnetically frustrated molybdate pyrochlores: $S=1$ oxide Lu$_2$Mo$_2$O$_7$ and $S={\frac{1}{2}}$ oxynitride Lu$_2$Mo$_2$O$_5$N$_2$. Lu$_2$Mo$_2$O$_7$ undergoes a transition to an unconventional spin glass ground state at $T_f {\sim} 16$ K. However, the preparation of the corresponding oxynitride tunes the nature of the ground state from spin glass to quantum spin liquid. The comparison of the static and dynamic spin correlations within the oxide and oxynitride phases presented here reveals the crucial role played by quantum fluctuations in the selection of a ground state. Furthermore, we estimate an upper limit for a gap in the spin excitation spectrum of the quantum spin liquid state of the oxynitride of $Δ {\sim} 0.05$ meV or ${\fracΔ{|θ|}}\sim0.004$, in units of its antiferromagnetic Weiss constant $θ {\sim}-121$ K. △ Less

Submitted 13 May, 2014; originally announced May 2014.

Comments: 8 pages, 9 figures

Helical order and multiferroicity in the $S=1/2$ quasi-kagome system KCu$_3$As$_2$O$_7$(OD)$_3$

Authors: Gøran J. Nilsen, Yoshihiko Okamoto, Hajime Ishikawa, Virginie Simonet, Claire V. Colin, Andres Cano, Laurent Chapon, Thomas Hansen, Hannu Mutka, Zenji Hiroi

Abstract: Several Cu$^{2+}$ hydroxide minerals have recently been identified as candidate realizations of the $S=1/2$ kagome Heisenberg model. In this context, we have studied the distorted system KCu$_3$As$_2$O$_7$(OD)$_3$ using neutron scattering and bulk measurements. Although the distortion favors magnetic order over a spin liquid ground state, refinement of the magnetic diffraction pattern below… ▽ More Several Cu$^{2+}$ hydroxide minerals have recently been identified as candidate realizations of the $S=1/2$ kagome Heisenberg model. In this context, we have studied the distorted system KCu$_3$As$_2$O$_7$(OD)$_3$ using neutron scattering and bulk measurements. Although the distortion favors magnetic order over a spin liquid ground state, refinement of the magnetic diffraction pattern below $T_{N1}=7.05(5)$ K yields a complex helical structure with $\mathbf{k}=(0.77~0~0.11)$. This structure, as well as the spin excitation spectrum, are well described by a classical Heisenberg model with ferromagnetic nearest neighbor couplings. Multiferroicity is observed below $T_{N1}$, with an unusual crossover between improper and pseudo-proper behavior occurring at $T_{N2}=5.5$ K. The polarization at $T=2$ K is $P=1.5 μ$Cm$^{-2}$. The properties of KCu$_3$As$_2$O$_7$(OD)$_3$ highlight the variety of physics which arise from the interplay of spin and orbital degrees of freedom in Cu$^{2+}$ kagome systems. △ Less

Submitted 11 April, 2014; originally announced April 2014.

Comments: Accepted for publication in Phys. Rev. B

X-ray Thomson scattering for partially ionized plasmas including the effect of bound levels

Authors: J. Nilsen, W. R. Johnson, K. T. Cheng

Abstract: X-ray Thomson scattering is being developed as a method to measure the temperature, electron density, and ionization state of high energy density plasmas such as those used in inertial confinement fusion. Most experiments are currently done at large laser facilities that can create bright X-ray sources, however the advent of the X-ray free electron laser (X-FEL) provides a new bright source to use… ▽ More X-ray Thomson scattering is being developed as a method to measure the temperature, electron density, and ionization state of high energy density plasmas such as those used in inertial confinement fusion. Most experiments are currently done at large laser facilities that can create bright X-ray sources, however the advent of the X-ray free electron laser (X-FEL) provides a new bright source to use in these experiments. One challenge with X-ray Thomson scattering experiments is understanding how to model the scattering for partially ionized plasmas in order to include the contributions of the bound electrons in the scattered intensity. In this work we take the existing models of Thomson scattering that include elastic ion-ion scattering and the electron-electron plasmon scattering and add the contribution of the bound electrons in the partially ionized plasmas. We validated our model by analyzing existing beryllium experimental data. We then consider several higher Z materials such as Cr and predict the existence of additional peaks in the scattering spectrum that requires new computational tools to understand. We also show examples of experiments in CH and Al that have bound contributions that change the shape of the scattered spectra. △ Less

Submitted 24 September, 2013; originally announced September 2013.

Comments: SPIE 2013 Optics and Photonics, San Diego, CA, United States August 25, 2013 through August 29, 2013. arXiv admin note: substantial text overlap with arXiv:1212.5972, arXiv:1207.5075

Report number: LLNL-PROC-643794

Electronic structure of KTi(SO4)2.H2O - a S=1/2 frustrated chain antiferromagnet

Authors: Deepa Kasinathan, K. Koepernik, O. Janson, G. J. Nilsen, J. O. Piatek, H. M. Ronnow, H. Rosner

Abstract: The compound KTi(SO4)2.H2O was recently reported as a quasi one-dimensional spin 1/2 compound with competing antiferromagnetic nearest neighbor exchange J1 and next-nearest neighbor exchange J2 along the chain with a frustration ratio alpha = J2/J1 ~ 0.29 [Chem. Mater. vol. 20, pg. 8 (2008)]. Here, we report a microscopically based magnetic model for this compound derived from density functional e… ▽ More The compound KTi(SO4)2.H2O was recently reported as a quasi one-dimensional spin 1/2 compound with competing antiferromagnetic nearest neighbor exchange J1 and next-nearest neighbor exchange J2 along the chain with a frustration ratio alpha = J2/J1 ~ 0.29 [Chem. Mater. vol. 20, pg. 8 (2008)]. Here, we report a microscopically based magnetic model for this compound derived from density functional electronic structure calculations along with respective tight-binding models. Our calculations confirm the quasi one-dimensional nature of the system with antiferromagnetic J1 and J2, but suggest a significantly larger frustration ratio alpha ~ 1.1 +- 0.2. Based on transfer matrix renormalization group calculations we found that, due to an intrinsic symmetry of the J1-J2 model, our larger frustration ratio alpha is also consistent with the previous thermodynamic data. To resolve this issue, we propose performing high-field magnetization measurements and low temperature susceptibility measurements which should allow to precisely identify the frustration ratio alpha. △ Less

Submitted 11 April, 2013; originally announced April 2013.

Comments: 11 pages, 8 figures

Breathing Pyrochlore Lattice Realized in A-Site Ordered Spinel Oxides LiGaCr4O8 and LiInCr4O8

Authors: Yoshihiko Okamoto, Gøran J. Nilsen, J. Paul Attfield, Zenji Hiroi

Abstract: A unique type of frustrated lattice is found in two A-site ordered spinel oxides, LiGaCr4O8 and LiInCr4O8. Because of the large size mismatch between Li+ and Ga3+/In3+ ions at the A site, the pyrochlore lattice, made up of Cr3+ ions carrying spin 3/2, becomes an alternating array of small and large tetrahedra, i.e., a "breathing" pyrochlore lattice. We introduce a parameter, the breathing factor B… ▽ More A unique type of frustrated lattice is found in two A-site ordered spinel oxides, LiGaCr4O8 and LiInCr4O8. Because of the large size mismatch between Li+ and Ga3+/In3+ ions at the A site, the pyrochlore lattice, made up of Cr3+ ions carrying spin 3/2, becomes an alternating array of small and large tetrahedra, i.e., a "breathing" pyrochlore lattice. We introduce a parameter, the breathing factor Bf, which quantifies the degree of frustration in the pyrochlore lattice: Bf is defined as J'/J, where J' and J are nearest-neighbor magnetic interactions in the large and small tetrahedra, respectively. LiGaCr4O8 with Bf ~ 0.6 shows magnetic susceptibility similar to that of conventional Cr spinel oxides such as ZnCr2O4. In contrast, LiInCr4O8 with a small Bf ~ 0.1 exhibits a spin-gap behavior in its magnetic susceptibility, suggesting a proximity to an exotic singlet ground state. Magnetic long-range order occurs at 13.8 and 15.9 K for LiGaCr4O8 and LiInCr4O8, respectively, in both cases likely owing to the coupling to structural distortions. △ Less

Submitted 29 January, 2013; originally announced January 2013.

Comments: 5 pages, 5 figures, accepted for publication in Phys. Rev. Lett

Resonant Bound-Free Contributions to Thomson Scattering of X-rays by Warm Dense Matter

Authors: W. R. Johnson, J. Nilsen, K. T. Cheng

Abstract: Recent calculations [Nilsen et al. arXiv:1212.5972] predict that contributions to the scattered photon spectrum from 3s and 3p bound states in chromium (Z=24) at metallic density and T=12 eV resonate below the respective bound-state thresholds. These resonances are shown to be closely related to continuum lowering, where 3d bound states in the free atom dissolve into a resonant l=2 partial wave in… ▽ More Recent calculations [Nilsen et al. arXiv:1212.5972] predict that contributions to the scattered photon spectrum from 3s and 3p bound states in chromium (Z=24) at metallic density and T=12 eV resonate below the respective bound-state thresholds. These resonances are shown to be closely related to continuum lowering, where 3d bound states in the free atom dissolve into a resonant l=2 partial wave in the continuum. The resulting d-state resonance dominates contributions to the bound-free dynamic structure function, leading to the predicted resonances in the scattered X-ray spectrum. Similar resonant features are shown to occur in all elements in the periodic table between Ca and Mn (20 <= Z <= 25). △ Less

Submitted 15 January, 2013; originally announced January 2013.

Comments: 6 pages, 3 figures, 1 table

Journal ref: High Energy Density Physics, 9 407-409 (2013)

The effect of bound states on X-ray Thomson scattering for partially ionized plasmas

Authors: J. Nilsen, W. R. Johnson, K. T. Cheng

Abstract: X-ray Thomson scattering is being developed as a method to measure the temperature, electron density, and ionization state of high energy density plasmas such as those used in inertial confinement fusion. X-ray laser sources have always been of interest because of the need to have a bright monochromatic x-ray source to overcome plasma emission and eliminate other lines in the background that compl… ▽ More X-ray Thomson scattering is being developed as a method to measure the temperature, electron density, and ionization state of high energy density plasmas such as those used in inertial confinement fusion. X-ray laser sources have always been of interest because of the need to have a bright monochromatic x-ray source to overcome plasma emission and eliminate other lines in the background that complicate the analysis. With the advent of the xray free electron laser (X-FEL) at the SLAC Linac Coherent Light Source (LCLS) and other facilities coming online worldwide, we now have such a source available in the keV regime. Most Thomson scattering codes used to model experimental data greatly simplify or neglect the contributions of the bound electrons to the scattered intensity. In this work we take the existing models of Thomson scattering that include elastic ion-ion scattering and inelastic electron-electron scattering and add the contribution of bound electrons in the partially ionized plasmas. To date, most experiments have studied hydrogen or beryllium plasmas. We first analyze existing experimental data for beryllium to validate the code. We then consider several higher Z materials such as Cr and predict the existence of additional peaks in the scattering spectrum that requires new computational tools to understand. For a Sn plasma, we show that bound contributions change the shape of the scattered spectrum in a way that would change the plasma temperature and density inferred from experiment. △ Less

Submitted 15 January, 2013; v1 submitted 24 December, 2012; originally announced December 2012.

Comments: Presented at the workshop "radiative Properties of Warm Dense Matter," Santa Barbara, CA, Nov. 6-10, 2012

Report number: LLNL-JRNL-608821

Journal ref: High Energy Density Physics, 9 388-391 (2013)

Average-Atom Model for X-ray Scattering from Warm Dense Matter

Authors: W. R. Johnson, J. Nilsen, K. T. Cheng

Abstract: A scheme for analyzing Thomson scattering of x-rays by warm dense matter, based on the average-atom model, is developed. Emphasis is given to x-ray scattering by bound electrons. Contributions to the scattered x-ray spectrum from elastic scattering by electrons moving with the ions and from inelastic scattering by free and bound electrons are evaluated using parameters (chemical potential, average… ▽ More A scheme for analyzing Thomson scattering of x-rays by warm dense matter, based on the average-atom model, is developed. Emphasis is given to x-ray scattering by bound electrons. Contributions to the scattered x-ray spectrum from elastic scattering by electrons moving with the ions and from inelastic scattering by free and bound electrons are evaluated using parameters (chemical potential, average ionic charge, free electron density, bound and continuum wave functions, and occupation numbers) taken from the average-atom model. The resulting scheme provides a relatively simple diagnostic for use in connection with x-ray scattering measurements. Applications are given to dense hydrogen, beryllium, aluminum, titanium, and tin plasmas. At high momentum transfer, contributions from inelastic scattering by bound electrons are dominant features of the scattered x-ray spectrum for aluminum, titanium, and tin. △ Less

Submitted 1 November, 2012; originally announced November 2012.

Comments: 22 pages, 10 figures Presentation at Workshop IV: Computational Challenges in Warm Dense Matter at IPAM (UCLA) May 21 - 25, 2012

Report number: LLNL-PROC-597272

A novel crystal polymorph of volborthite, Cu3V2O7(OH)2 2H2O

Authors: Hajime Ishikawa, Jun-ichi Yamaura, Yoshihiko Okamoto, Hiroyuki Yoshida, Gøran J. Nilsen, Zenji Hiroi

Abstract: A new polymorph of volborthite [tricopper(II) divanadium(V) heptaoxide dihydroxide dihydrate], Cu3V2O7(OH)2 2H2O, has been discovered in a single crystal prepared by hydrothermal synthesis. X-ray analysis reveals that the monoclinic structure has the space group C2/c at room temperature, which is different from that of the previously reported C2/m structure. Both structures have Cu3O6(OH)2 layers… ▽ More A new polymorph of volborthite [tricopper(II) divanadium(V) heptaoxide dihydroxide dihydrate], Cu3V2O7(OH)2 2H2O, has been discovered in a single crystal prepared by hydrothermal synthesis. X-ray analysis reveals that the monoclinic structure has the space group C2/c at room temperature, which is different from that of the previously reported C2/m structure. Both structures have Cu3O6(OH)2 layers composed of edge-sharing CuO4(OH)2 octahedra, with V2O7 pillars and water molecules between the layers. The Cu atoms occupy two and three independent crystallographic sites in the C2/m and C2/c structures, respectively, likely giving rise to different magnetic interactions between CuII spins in the kagome lattices embedded in the Cu3O6(OH)2 layers. △ Less

Submitted 14 September, 2012; originally announced September 2012.

Comments: 7pages, 2figures

Journal ref: Acta Cryst. C68, i41-i44 (2012)

Using the X-FEL to understand X-ray Thomson scattering for partially ionized plasmas

Authors: J. Nilsen, W. R. Johnson, K. T. Cheng

Abstract: For the last decade numerous researchers have been trying to develop experimental techniques to use X-ray Thomson scattering as a method to measure the temperature, electron density, and ionization state of high energy density plasmas such as those used in inertial confinement fusion. With the advent of the X-ray free electron laser (X-FEL) at the SLAC Linac Coherent Light Source (LCLS) we now hav… ▽ More For the last decade numerous researchers have been trying to develop experimental techniques to use X-ray Thomson scattering as a method to measure the temperature, electron density, and ionization state of high energy density plasmas such as those used in inertial confinement fusion. With the advent of the X-ray free electron laser (X-FEL) at the SLAC Linac Coherent Light Source (LCLS) we now have such a source available in the keV regime. One challenge with X-ray Thomson scattering experiments is understanding how to model the scattering for partially ionized plasmas. Most Thomson scattering codes used to model experimental data greatly simplify or neglect the contributions of the bound electrons to the scattered intensity. In this work we take the existing models of Thomson scattering that include elastic ion-ion scattering and the electron-electron plasmon scattering and add the contribution of the bound electrons in the partially ionized plasmas. Except for hydrogen plasmas almost every plasma that is studied today has bound electrons and it is important to understand their contribution to the Thomson scattering, especially as new X-ray sources such as the X-FEL will allow us to study much higher Z plasmas. Currently most experiments have looked at hydrogen or beryllium. We will first look at the bound electron contributions to beryllium by analysing existing experimental data. We then consider several higher Z materials such as Cr and predict the existence of additional peaks in the scattering spectrum that requires new computational tools to understand. For a Sn plasma we show that the bound contributions changes the shape of the scattered spectrum in a way that would change the plasma temperature and density inferred by the experiment. △ Less

Submitted 1 October, 2012; v1 submitted 20 July, 2012; originally announced July 2012.

Comments: 13th International Conference on X-ray Lasers Paris, France June 10, 2012 through June 15, 2012

Report number: LLNL-PROC-564720

Thomson scattering in the average-atom approximation

Authors: W. R. Johnson, J. Nilsen, K. T. Cheng

Abstract: The average-atom model is applied to study Thomson scattering of x-rays from warm-dense matter with emphasis on scattering by bound electrons. Parameters needed to evaluate the dynamic structure function (chemical potential, average ionic charge, free electron density, bound and continuum wave-functions and occupation numbers) are obtained from the average-atom model. The resulting analysis provid… ▽ More The average-atom model is applied to study Thomson scattering of x-rays from warm-dense matter with emphasis on scattering by bound electrons. Parameters needed to evaluate the dynamic structure function (chemical potential, average ionic charge, free electron density, bound and continuum wave-functions and occupation numbers) are obtained from the average-atom model. The resulting analysis provides a relatively simple diagnostic for use in connection with x-ray scattering measurements. Applications are given to dense hydrogen, beryllium, aluminum and titanium plasmas. In the case of titanium, bound states are predicted to modify the spectrum significantly. △ Less

Submitted 1 October, 2012; v1 submitted 1 July, 2012; originally announced July 2012.

Comments: 9 pages, 9 figures

Journal ref: Physical Review E 86, 036410 (2012)

Distorted Kagome Lattice Generated by a Unique Orbital Arrangement in the Copper Mineral KCu3As2O7(OH)3

Authors: Yoshihiko Okamoto, Hajime Ishikawa, Gøran J. Nilsen, Zenji Hiroi

Abstract: We study polycrystalline samples of KCu3As2O7(OH)3, a new candidate spin-1/2 kagome antiferromagnet, by magnetic susceptibility and heat capacity measurements above 2 K. The unique arrangement of the 3z2 - r2 and x2 - y2 orbitals on the Cu2+ kagome net is noted and compared to the orbital patterns found in other kagome minerals. It is suggested that this orbital arrangement gives rise to one antif… ▽ More We study polycrystalline samples of KCu3As2O7(OH)3, a new candidate spin-1/2 kagome antiferromagnet, by magnetic susceptibility and heat capacity measurements above 2 K. The unique arrangement of the 3z2 - r2 and x2 - y2 orbitals on the Cu2+ kagome net is noted and compared to the orbital patterns found in other kagome minerals. It is suggested that this orbital arrangement gives rise to one antiferromagnetic and two ferromagnetic interactions on isosceles triangles forming a highly distorted kagome lattice. KCu3As2O7(OH)3 is found to show an antiferromagnetic long-range order at T_N = 7.2 K. Remarkably, a spin entropy is more gradually released upon cooling below T_N compared with a conventional magnetic long-range order, which may originate from the geometrical frustration still present in this highly distorted kagome lattice. △ Less

Submitted 29 January, 2013; v1 submitted 13 February, 2012; originally announced February 2012.

Comments: 5 pages, 4 figures

Journal ref: J. Phys. Soc. Jpn. 81, 033707 (2012)

Performance and Stability of the Chelonia Storage Cloud

Authors: Jon K. Nilsen, Salman Toor, Zsombor Nagy, Bjarte Mohn, Alex L. Read

Abstract: In this paper we present the Chelonia storage cloud middleware. It was designed to fill the requirements gap between those of large, sophisticated scientific collaborations which have adopted the grid paradigm for their distributed storage needs, and of corporate business communities which are gravitating towards the cloud paradigm. The similarities to and differences between Chelonia and severa… ▽ More In this paper we present the Chelonia storage cloud middleware. It was designed to fill the requirements gap between those of large, sophisticated scientific collaborations which have adopted the grid paradigm for their distributed storage needs, and of corporate business communities which are gravitating towards the cloud paradigm. The similarities to and differences between Chelonia and several well-known grid- and cloud-based storage solutions are commented. The design of Chelonia has been chosen to optimize high reliability and scalability of an integrated system of heterogeneous, geographically dispersed storage sites and the ability to easily expand the system dynamically. The architecture and implementation in term of web-services running inside the Advanced Resource Connector Hosting Environment Dameon (ARC HED) are described. We present results of tests in both local-area and wide-area networks that demonstrate the fault-tolerance, stability and scalability of Chelonia. △ Less

Submitted 3 February, 2010; originally announced February 2010.

Comments: 29 pages, 10 figures, submitted to Future Generation Computing Systems

Simplifying Parallelization of Scientific Codes by a Function-Centric Approach in Python

Authors: Jon K. Nilsen, Xing Cai, Bjorn Hoyland, Hans Petter Langtangen

Abstract: The purpose of this paper is to show how existing scientific software can be parallelized using a separate thin layer of Python code where all parallel communication is implemented. We provide specific examples on such layers of code, and these examples may act as templates for parallelizing a wide set of serial scientific codes. The use of Python for parallelization is motivated by the fact tha… ▽ More The purpose of this paper is to show how existing scientific software can be parallelized using a separate thin layer of Python code where all parallel communication is implemented. We provide specific examples on such layers of code, and these examples may act as templates for parallelizing a wide set of serial scientific codes. The use of Python for parallelization is motivated by the fact that the language is well suited for reusing existing serial codes programmed in other languages. The extreme flexibility of Python with regard to handling functions makes it very easy to wrap up decomposed computational tasks of a serial scientific application as Python functions. Many parallelization-specific components can be implemented as generic Python functions, which may take as input those functions that perform concrete computational tasks. The overall programming effort needed by this parallelization approach is rather limited, and the resulting parallel Python scripts have a compact and clean structure. The usefulness of the parallelization approach is exemplified by three different classes of applications in natural and social sciences. △ Less

Submitted 3 February, 2010; originally announced February 2010.

Comments: 29 pages, submitted to Computational Science and Discovery