Showing 1–50 of 54 results for author: Patrick, C E

Origin of Incommensurate Magnetic Order in Rare-Earth Magnetic Weyl Semimetals

Authors: Juba Bouaziz, Gustav Bihlmayer, Christopher E. Patrick, Julie B. Staunton, Stefan Blügel

Abstract: We investigate rare-earth magnetic Weyl semimetals through first-principles simulations, analyzing the connection between incommensurate magnetic order and the presence of Weyl nodes in the electronic band structure. Focusing on PrAlSi, NdAlSi, and SmAlSi, we demonstrate that the reported helical ordering does not originate from the nesting of topological features at the Fermi Surface or the Dzyal… ▽ More We investigate rare-earth magnetic Weyl semimetals through first-principles simulations, analyzing the connection between incommensurate magnetic order and the presence of Weyl nodes in the electronic band structure. Focusing on PrAlSi, NdAlSi, and SmAlSi, we demonstrate that the reported helical ordering does not originate from the nesting of topological features at the Fermi Surface or the Dzyaloshinskii-Moriya interaction. Instead, the helical order arises from frustrated isotropic short-range superexchange between the 4f moments facilitated by pd-hybridization with the main group elements. Employing a spin Hamiltonian with isotropic exchange and single-ion anisotropy we replicate the experimentally observed helical modulation. △ Less

Submitted 13 November, 2023; originally announced November 2023.

Comments: 7 pages, 3 figures

Theory of defect-induced crystal field perturbations in rare earth magnets

Authors: Christopher E. Patrick, Yixuan Huang, Laura H. Lewis, Julie B. Staunton

Abstract: We present a theory describing the single-ion anisotropy of rare earth (RE) magnets in the presence of point defects. Taking the RE-lean 1:12 magnet class as a prototype, we use first-principles calculations to show how the introduction of Ti substitutions into SmFe$_{12}$ perturbs the crystal field, generating new coefficients due to the lower symmetry of the RE environment. We then demonstrate t… ▽ More We present a theory describing the single-ion anisotropy of rare earth (RE) magnets in the presence of point defects. Taking the RE-lean 1:12 magnet class as a prototype, we use first-principles calculations to show how the introduction of Ti substitutions into SmFe$_{12}$ perturbs the crystal field, generating new coefficients due to the lower symmetry of the RE environment. We then demonstrate that these perturbations can be described extremely efficiently using a screened point charge model. We provide analytical expressions for the anisotropy energy which can be straightforwardly implemented in atomistic spin dynamics simulations, meaning that such simulations can be carried out for an arbitrary arrangement of point defects. The significant crystal field perturbations calculated here demonstrate that a sample which is single-phase from a structural point of view can nonetheless have a dramatically varying anisotropy profile at the atomistic level if there is compositional disorder, which may influence localized magnetic objects like domain walls or skyrmions. △ Less

Submitted 22 August, 2023; originally announced August 2023.

Comments: 14 pages, 3 figures

Revisiting Néel 60 years on: the magnetic anisotropy of $\mathrm{L}1_0$ FeNi (tetrataenite)

Authors: Christopher D. Woodgate, Christopher E. Patrick, Laura H. Lewis, Julie B. Staunton

Abstract: The magnetocrystalline anisotropy energy of atomically ordered $\mathrm{L}1_0$ FeNi (the meteoritic mineral tetrataenite) is studied within a first-principles electronic structure framework. Two compositions are examined: equiatomic Fe$_{0.5}$Ni$_{0.5}$ and an Fe-rich composition, Fe$_{0.56}$Ni$_{0.44}$. It is confirmed that, for the single crystals modelled in this work, the leading-order anisotr… ▽ More The magnetocrystalline anisotropy energy of atomically ordered $\mathrm{L}1_0$ FeNi (the meteoritic mineral tetrataenite) is studied within a first-principles electronic structure framework. Two compositions are examined: equiatomic Fe$_{0.5}$Ni$_{0.5}$ and an Fe-rich composition, Fe$_{0.56}$Ni$_{0.44}$. It is confirmed that, for the single crystals modelled in this work, the leading-order anisotropy coefficient $K_1$ dominates the higher-order coefficients $K_2$ and $K_3$. To enable comparison with experiment, the effects of both imperfect atomic long-range order and finite temperature are included. While our computational results initially appear to undershoot the measured experimental values for this system, careful scrutiny of the original analysis due to Néel et al. [J. Appl. Phys. 35, 873 (1964)] suggests that our computed value of $K_1$ is, in fact, consistent with experimental values, and that the noted discrepancy has its origins in the nanoscale polycrystalline, multivariant nature of experimental samples, that yields much larger values of $K_2$ and $K_3$ than expected a priori. These results provide fresh insight into the existing discrepancies in the literature regarding the value of tetrataenite's uniaxial magnetocrystalline anisotropy in both natural and synthetic samples. △ Less

Submitted 12 October, 2023; v1 submitted 28 July, 2023; originally announced July 2023.

Comments: 9 pages, 3 figures, 2 tables

Journal ref: J. Appl. Phys. 134, 163905 (2023)

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.

Quantification of electronic and magnetoelastic mechanisms of first-order magnetic phase transitions from first principles: application to caloric effects in La(Fe$_x$Si$_{1-x}$)$_{13}$

Authors: Eduardo Mendive Tapia, Christopher E. Patrick, Tilmann Hickel, Jörg Neugebauer, Julie B. Staunton

Abstract: La(Fe$_x$Si$_{1-x}$)$_{13}$ and derived quaternary compounds are well-known for their giant, tunable, magneto- and barocaloric responses around a first-order paramagnetic-ferromagnetic transition near room temperature with low hysteresis. Remarkably, such a transition shows a large spontaneous volume change together with itinerant electron metamagnetic features. While magnetovolume effects are wel… ▽ More La(Fe$_x$Si$_{1-x}$)$_{13}$ and derived quaternary compounds are well-known for their giant, tunable, magneto- and barocaloric responses around a first-order paramagnetic-ferromagnetic transition near room temperature with low hysteresis. Remarkably, such a transition shows a large spontaneous volume change together with itinerant electron metamagnetic features. While magnetovolume effects are well-established mechanisms driving first-order transitions, purely electronic sources have a long, subtle history and remain poorly understood. Here we apply a disordered local moment picture to quantify electronic and magnetoelastic effects at finite temperature in La(Fe$_x$Si$_{1-x}$)$_{13}$ from first-principles. We obtain results in very good agreement with experiment and demonstrate that the magnetoelastic coupling, rather than purely electronic mechanisms, drives the first-order character and causes at the same time a huge electronic entropy contribution to the caloric response. △ Less

Submitted 24 October, 2023; v1 submitted 13 February, 2023; originally announced February 2023.

Comments: 20 pages, 6 figures

Journal ref: Journal of Physics: Energy, 5, 034004 (2023)

Crucial role of Fe in determining the hard magnetic properties of Nd$_2$Fe$_{14}$B

Authors: Juba Bouaziz, Christopher E. Patrick, Julie B. Staunton

Abstract: Nd$_2$Fe$_{14}$B's unsurpassed, hard magnetic properties for a wide range of temperatures result from a combination of a large volume magnetization from Fe and a strong single-ion anisotropy from Nd. Here, using finite temperature first-principles calculations, we focus on the other crucial roles played by the Fe atoms in maintaining the magnetic order on the Nd sublattices, and hence the large ma… ▽ More Nd$_2$Fe$_{14}$B's unsurpassed, hard magnetic properties for a wide range of temperatures result from a combination of a large volume magnetization from Fe and a strong single-ion anisotropy from Nd. Here, using finite temperature first-principles calculations, we focus on the other crucial roles played by the Fe atoms in maintaining the magnetic order on the Nd sublattices, and hence the large magnetic anisotropy, and directly generating significant uniaxial anisotropy at high temperatures. We identify effective spins for atomistic modelling from the material's interacting electrons and {quantify pairwise and higher order, non-pairwise magnetic interactions among them. We find the Nd spins couple most strongly to spins on sites belonging to two specific Fe sublattices, 8$j_1$, 8$j_2$. Moreover the Fe 8$j_1$ sublattice also provides the electronic origin of the unusual, nonmonotonic temperature dependence of the anisotropy of Y$_2$Fe$_{14}$B.} Our work provides atomic-level resolution of the properties of this fascinating magnetic material. △ Less

Submitted 7 January, 2023; originally announced January 2023.

Journal ref: Phys. Rev. B 107, L020401 (2023)

Spin orientation and magnetostriction of Tb$_{1-x}$Dy$_{x}$Fe$_{2}$ from first principles

Authors: Christopher E. Patrick, George A. Marchant, Julie B. Staunton

Abstract: The optimal amount of dysprosium in the highly magnetostrictive rare-earth compounds Tb$_{1-x}$Dy$_x$Fe$_2$ for room temperature applications has long been known to be $x$=0.73 (Terfenol-D). Here, we derive this value from first principles by calculating the easy magnetization direction and magnetostriction as a function of composition and temperature. We use crystal field coefficients obtained wi… ▽ More The optimal amount of dysprosium in the highly magnetostrictive rare-earth compounds Tb$_{1-x}$Dy$_x$Fe$_2$ for room temperature applications has long been known to be $x$=0.73 (Terfenol-D). Here, we derive this value from first principles by calculating the easy magnetization direction and magnetostriction as a function of composition and temperature. We use crystal field coefficients obtained within density-functional theory to construct phenomenological anisotropy and magnetoelastic constants. The temperature dependence of these constants is obtained from disordered local moment calculations of the rare earth magnetic order parameter. Our calculations find the critical Dy concentration required to switch the magnetization direction at room temperature to be $x_c$=0.78, with magnetostrictions $λ_{111}$=2700 and $λ_{100}$=-430~ppm, close to the Terfenol-D values. △ Less

Submitted 30 July, 2020; originally announced July 2020.

Comments: 27 pages, 4 figures

Journal ref: Phys. Rev. Applied 14, 014091 (2020)

Temperature-dependent magnetocrystalline anisotropy of rare earth/transition metal permanent magnets from first principles: The light RCo$_5$ (R=Y, La-Gd) intermetallics

Authors: Christopher E. Patrick, Julie B. Staunton

Abstract: Computational design of more efficient rare earth/transition metal (RE-TM) permanent magnets requires accurately calculating the magnetocrystalline anisotropy (MCA) at finite temperature, since this property places an upper bound on the coercivity. Here, we present a first-principles methodology to calculate the MCA of RE-TM magnets which fully accounts for the effects of temperature on the underl… ▽ More Computational design of more efficient rare earth/transition metal (RE-TM) permanent magnets requires accurately calculating the magnetocrystalline anisotropy (MCA) at finite temperature, since this property places an upper bound on the coercivity. Here, we present a first-principles methodology to calculate the MCA of RE-TM magnets which fully accounts for the effects of temperature on the underlying electrons. The itinerant electron TM magnetism is described within the disordered local moment picture, and the localized RE-4f magnetism is described within crystal field theory. We use our model, which is free of adjustable parameters, to calculate the MCA of the RCo$_5$ (R=Y, La-Gd) magnet family for temperatures 0--600 K. We correctly find a huge uniaxial anisotropy for SmCo$_5$ (21.3 MJm$^{-3}$ at 300 K) and two finite temperature spin reorientation transitions for NdCo$_5$. The calculations also demonstrate dramatic valency effects in CeCo$_5$ and PrCo$_5$. Our calculations provide quantitative, first-principles insight into several decades of RE-TM experimental studies. △ Less

Submitted 16 October, 2019; originally announced October 2019.

Comments: 8 pages, 4 figures

Journal ref: Phys. Rev. Materials 3, 101401(R) (2019)

Accurate ground state- and quasiparticle energies: beyond the RPA and GW methods with adiabatic exchange-correlation kernels

Authors: Thomas Olsen, Christopher E. Patrick, Jefferson E. Bates, Adrienn Ruzsinszky, Kristian S. Thygesen

Abstract: We review the theory and application of adiabatic exchange-correlation (xc-) kernels for ab initio calculations of ground state energies and quasiparticle excitations within the frameworks of the adiabatic connection fluctuation dissipation theorem and Hedin's equations, respectively. Various different xc-kernels, which are all rooted in the homogeneous electron gas, are introduced but hereafter w… ▽ More We review the theory and application of adiabatic exchange-correlation (xc-) kernels for ab initio calculations of ground state energies and quasiparticle excitations within the frameworks of the adiabatic connection fluctuation dissipation theorem and Hedin's equations, respectively. Various different xc-kernels, which are all rooted in the homogeneous electron gas, are introduced but hereafter we focus on the specific class of renormalized adiabatic kernels, in particular the rALDA and rAPBE. The kernels drastically improve the description of short-range correlations as compared to the random phase approximation (RPA), resulting in significantly better correlation energies. This effect greatly reduces the reliance on error cancellations, which is essential in RPA, and systematically improves covalent bond energies while preserving the good performance of the RPA for dispersive interactions. For quasiparticle energies, the xc-kernels account for vertex corrections that are missing in the GW self-energy. In this context, we show that the short-range correlations mainly correct the absolute band positions while the band gap is less affected in agreement with the known good performance of GW for the latter. The renormalized xc-kernels offer a rigorous extension of the RPA and GW methods with clear improvements in terms of accuracy at little extra computational cost. △ Less

Submitted 18 June, 2019; originally announced June 2019.

Comments: Review article. 53 pages-

Journal ref: npj Computational Materials 5, 106 (2019)

Crystal field coefficients for yttrium analogues of rare-earth/transition-metal magnets using density-functional theory in the projector-augmented wave formalism

Authors: Christopher E. Patrick, Julie B. Staunton

Abstract: We present a method of calculating crystal field coefficients of rare-earth/transition-metal (RE-TM) magnets within density-functional theory (DFT). The principal idea of the method is to calculate the crystal field potential of the yttrium analogue ("Y-analogue") of the RE-TM magnet, i.e. the material where the lanthanide elements have been substituted with yttrium. The advantage of dealing with… ▽ More We present a method of calculating crystal field coefficients of rare-earth/transition-metal (RE-TM) magnets within density-functional theory (DFT). The principal idea of the method is to calculate the crystal field potential of the yttrium analogue ("Y-analogue") of the RE-TM magnet, i.e. the material where the lanthanide elements have been substituted with yttrium. The advantage of dealing with Y-analogues is that the methodological and conceptual difficulties associated with treating the highly-localized 4f electrons in DFT are avoided, whilst the nominal valence electronic structure principally responsible for the crystal field is preserved. In order to correctly describe the crystal field potential in the core region of the atoms we use the projector-augmented wave formalism of DFT, which allows the reconstruction of the full charge density and electrostatic potential. The Y-analogue crystal field potentials are combined with radial 4f charge densities obtained in self-interaction-corrected calculations on the lanthanides to obtain crystal field coefficients. We demonstrate our method on a test set of 10 materials comprising 9 RE-TM magnets and elemental Tb. We show that the calculated easy directions of magnetization agree with experimental observations, including a correct description of the anisotropy within the basal plane of Tb and NdCo$_5$. We further show that the Y-analogue calculations generally agree quantitatively with previous calculations using the open-core approximation to treat the 4f electrons, and argue that our simple approach may be useful for large-scale computational screening of new magnetic materials. △ Less

Submitted 15 April, 2019; originally announced April 2019.

Comments: 22 pages, 2 figures Accepted for publication in Journal of Physics: Condensed Matter

Journal ref: J. Phys.: Condens. Matter 31, 305901 (2019)

Structural and magnetic properties of GdCo$_{5-x}$Ni$_x$

Authors: Amy L. Tedstone, Christopher E. Patrick, Santosh Kumar, Rachel S. Edwards, Martin R. Lees, Geetha Balakrishnan, Julie B. Staunton

Abstract: GdCo$_5$ may be considered as two sublattices - one of Gd and one of Co - whose magnetizations are in antiparallel alignment, forming a ferrimagnet. Substitution of nickel in the cobalt sublattice of GdCo$_5$ has been investigated to gain insight into how the magnetic properties of this prototype rare-earth/transition-metal magnet are affected by changes in the transition metal sublattice. Polycry… ▽ More GdCo$_5$ may be considered as two sublattices - one of Gd and one of Co - whose magnetizations are in antiparallel alignment, forming a ferrimagnet. Substitution of nickel in the cobalt sublattice of GdCo$_5$ has been investigated to gain insight into how the magnetic properties of this prototype rare-earth/transition-metal magnet are affected by changes in the transition metal sublattice. Polycrystalline samples of GdCo$_{5-x}$Ni$_x$ for 0 $ \leq x \leq $ 5 were synthesized by arc melting. Structural characterization was carried out by powder x-ray diffraction and optical and scanning electron microscope imaging of metallographic slides, the latter revealing a low concentration of Gd$_2$(Co, Ni)$_7$ lamellae for $x \leq 2.5$. Compensation - i.e. the cancellation of the opposing Gd and transition metal moments is observed for $1 \leq x \leq 3$ at a temperature which increases with Ni content; for larger $x$, no compensation is observed below 360 K. A peak in the coercivity is seen at $x \approx 1$ at 10K coinciding with a minimum in the saturation magnetization. Density-functional theory calculations within the disordered local moment picture reproduce the dependence of the magnetization on Ni content and temperature. The calculations also show a peak in the magnetocrystalline anisotropy at similar Ni concentrations to the experimentally observed coercivity maximum. △ Less

Submitted 27 March, 2019; originally announced March 2019.

Comments: 12 pages, 12 figures, 2 tables

Journal ref: Phys. Rev. Materials 3, 034409 (2019)

First-principles calculations of the magnetocrystalline anisotropy of the prototype 2:17 cell boundary phase Y(Co$_{1-x-y}$Fe$_x$Cu$_y$)$_5$

Authors: Christopher E. Patrick, Munehisa Matsumoto, Julie B. Staunton

Abstract: We present a computational study of the compound Y(Co$_{1-x-y}$Fe$_x$Cu$_y$)$_5$ for 0 $\leq x,y \leq 0.2$. This compound was chosen as a prototype for investigating the cell boundary phase believed to play a key role in establishing the high coercivity of commercial Sm-Co 2:17 magnets. Using density-functional theory, we have calculated the magnetization and magnetocrystalline anisotropy at zero… ▽ More We present a computational study of the compound Y(Co$_{1-x-y}$Fe$_x$Cu$_y$)$_5$ for 0 $\leq x,y \leq 0.2$. This compound was chosen as a prototype for investigating the cell boundary phase believed to play a key role in establishing the high coercivity of commercial Sm-Co 2:17 magnets. Using density-functional theory, we have calculated the magnetization and magnetocrystalline anisotropy at zero temperature for a range of compositions, modeling the doped compounds within the coherent potential approximation. We have also performed finite temperature calculations for YCo$_5$, Y(Co$_{0.838}$Cu$_{0.162}$)$_5$ and Y(Co$_{0.838}$Fe$_{0.081}$Cu$_{0.081}$)$_5$ within the disordered local moment picture. Our calculations find that substituting Co with small amounts of either Fe or Cu boosts the magnetocrystalline anisotropy $K$, but the change in $K$ depends strongly on the location of the dopants. Furthermore, the calculations do not show a particularly large difference between the magnetic properties of Cu-rich Y(Co$_{0.838}$Cu$_{0.162}$)$_5$ and equal Fe-Cu Y(Co$_{0.838}$Fe$_{0.081}$Cu$_{0.081}$)$_5$, despite these two compositions showing different coercivity behavior when found in the cell boundary phase of 2:17 magnets. Our study lays the groundwork for studying the rare earth contribution to the anisotropy of Sm(Co$_{1-x-y}$Fe$_x$Cu$_y$)$_5$, and also shows how a small amount of transition metal substitution can boost the anisotropy field of YCo$_5$. △ Less

Submitted 14 January, 2019; originally announced January 2019.

Comments: 23 pages, 7 figures

Measurement of Quasielastic-Like Neutrino Scattering at $\left< E_ν\right> \sim 3.5$~ GeV on a Hydrocarbon Target

Authors: D. Ruterbories, K. Hurtado, J. Osta, F. Akbar, L. Aliaga, D. A. Andrade, M. V. Ascencio, A. Bashyal, A. Bercellie, M. Betancourt, A. Bodek, H. Budd, G. Caceres, T. Cai, M. F. Carneiro, J. Chaves, D. Coplowe, H. da Motta, S. A. Dytman, G. A. Díaz, J. Felix, L. Fields, A. Filkins, R. Fine, A. M. Gago , et al. (48 additional authors not shown)

Abstract: MINERvA presents a new analysis of neutrino induced quasielastic-like interactions in a hydrocarbon tracking target. We report a double-differential cross section using the muon transverse and longitudinal momentum. In addition, differential cross sections as a function of the square of the four-momentum transferred and the neutrino energy are calculated using a quasielastic hypothesis. Finally, a… ▽ More MINERvA presents a new analysis of neutrino induced quasielastic-like interactions in a hydrocarbon tracking target. We report a double-differential cross section using the muon transverse and longitudinal momentum. In addition, differential cross sections as a function of the square of the four-momentum transferred and the neutrino energy are calculated using a quasielastic hypothesis. Finally, an analysis of energy deposited near the interaction vertex is presented. These results are compared to modified GENIE predictions as well as a NuWro prediction. All results use a dataset produced by $3.34\times10^{20}$ protons on target creating a neutrino beam with a peak energy of approximately 3.5 GeV △ Less

Submitted 7 November, 2018; originally announced November 2018.

Comments: 38 pages, 47 figures

Report number: FERMILAB-PUB-18-618-ND

Journal ref: Phys. Rev. D 99, 012004 (2019)

Field-induced canting of magnetic moments in GdCo5 at finite temperature: first-principles calculations and high-field measurements

Authors: Christopher E. Patrick, Santosh Kumar, Kathrin Götze, Matthew J. Pearce, John Singleton, George Rowlands, Geetha Balakrishnan, Martin R. Lees, Paul A. Goddard, Julie B. Staunton

Abstract: We present calculations and experimental measurements of the temperature-dependent magnetization of a single crystal of GdCo$_5$ in magnetic fields of order 60 T. At zero temperature the calculations, based on density-functional theory in the disordered-local-moment picture, predict a field-induced transition from an antiferromagnetic to a canted alignment of Gd and Co moments at 46.1 T. At higher… ▽ More We present calculations and experimental measurements of the temperature-dependent magnetization of a single crystal of GdCo$_5$ in magnetic fields of order 60 T. At zero temperature the calculations, based on density-functional theory in the disordered-local-moment picture, predict a field-induced transition from an antiferromagnetic to a canted alignment of Gd and Co moments at 46.1 T. At higher temperatures the calculations find this critical field to increase along with the zero-field magnetization. The experimental measurements observe this transition to occur between 44-48 T at 1.4 K. Up to temperatures of at least 100 K, the experiments continue to observe the transition; however, at variance with the calculations, no strong temperature dependence of the critical field is apparent. We assign this difference to the inaccurate description of the zero-field magnetization of the calculations at low temperatures, due to the use of classical statistical mechanics. Correcting for this effect, we recover a consistent description of the high-field magnetization of GdCo$_5$ from theory and experiment. △ Less

Submitted 25 July, 2018; originally announced July 2018.

Comments: 11 pages, 4 figures. Data presented in this paper are available at http://wrap.warwick.ac.uk/103805

Journal ref: J. Phys.: Condens. Matter 30, 32LT01 (2018)

Rare-earth/transition-metal magnets at finite temperature: Self-interaction-corrected relativistic density functional theory in the disordered local moment picture

Authors: Christopher E. Patrick, Julie B. Staunton

Abstract: Atomic-scale computational modeling of technologically relevant permanent magnetic materials faces two key challenges. First, a material's magnetic properties depend sensitively on temperature, so the calculations must account for thermally induced magnetic disorder. Second, the most widely-used permanent magnets are based on rare-earth elements, whose highly localized 4$f$ electrons are poorly de… ▽ More Atomic-scale computational modeling of technologically relevant permanent magnetic materials faces two key challenges. First, a material's magnetic properties depend sensitively on temperature, so the calculations must account for thermally induced magnetic disorder. Second, the most widely-used permanent magnets are based on rare-earth elements, whose highly localized 4$f$ electrons are poorly described by standard electronic structure methods. Here, we take two established theories, the disordered local moment picture of thermally induced magnetic disorder and self-interaction-corrected density functional theory, and devise a computational framework to overcome these challenges. Using the new approach, we calculate magnetic moments and Curie temperatures of the rare-earth cobalt (RECo$_5$) family for RE=Y--Lu. The calculations correctly reproduce the experimentally measured trends across the series and confirm that, apart from the hypothetical compound EuCo$_5$, SmCo$_5$ has the strongest magnetic properties at high temperature. An order-parameter analysis demonstrates that varying the RE has a surprisingly strong effect on the Co--Co magnetic interactions determining the Curie temperature, even when the lattice parameters are kept fixed. We propose the origin of this behavior is a small contribution to the density from $f$-character electrons located close to the Fermi level. △ Less

Submitted 14 June, 2018; originally announced June 2018.

Comments: 17 pages, 11 figures

Journal ref: Phys. Rev. B 97, 224415 (2018)

Measurement of final-state correlations in neutrino muon-proton mesonless production on hydrocarbon at $\langle E_ν\rangle=3$ GeV

Authors: X. -G. Lu, M. Betancourt, T. Walton, F. Akbar, L. Aliaga, O. Altinok, D. A. Andrade, M. Ascencio, L. Bellantoni, A. Bercellie, A. Bodek, A. Bravar, H. Budd, T. Cai, M. F. Carneiro, J. Chaves, D. Coplowe, H. da Motta, S. A. Dytman, G. A. Diaz, J. Felix, L. Fields, R. Fine, A. M. Gago, R. Galindo , et al. (47 additional authors not shown)

Abstract: Final-state kinematic imbalances are measured in mesonless production of $ν_μ+ A \to μ^- + p + X$ in the MINERvA tracker. Initial- and final-state nuclear effects are probed using the direction of the $μ^-$-p transverse momentum imbalance and the initial-state momentum of the struck neutron. Differential cross sections are compared to predictions based on current approaches to medium modeling. The… ▽ More Final-state kinematic imbalances are measured in mesonless production of $ν_μ+ A \to μ^- + p + X$ in the MINERvA tracker. Initial- and final-state nuclear effects are probed using the direction of the $μ^-$-p transverse momentum imbalance and the initial-state momentum of the struck neutron. Differential cross sections are compared to predictions based on current approaches to medium modeling. These models under-predict the cross section at intermediate intranuclear momentum transfers that generally exceed the Fermi momenta. As neutrino interaction models need to correctly incorporate the effect of the nucleus in order to predict neutrino energy resolution in oscillation experiments, this result points to a region of phase space where additional cross section strength is needed in current models, and demonstrates a new technique that would be suitable for use in fine grained liquid argon detectors where the effect of the nucleus may be even larger. △ Less

Submitted 18 September, 2018; v1 submitted 14 May, 2018; originally announced May 2018.

Comments: 6 pages, 6 figures, accepted by PRL, data release package updated

Report number: FERMILAB-PUB-18-179-ND

Journal ref: Phys. Rev. Lett. 121, 022504 (2018)

Antineutrino Charged-Current reactions on Hydrocarbon with Low Momentum Transfer

Authors: R. Gran, M. Betancourt, M. Elkins, P. A. Rodrigues, F. Akbar, L. Aliaga, D. A. Andrade, A. Bashyal, L. Bellantoni, A. Bercellie, A. Bodek, A. Bravar, H. Budd, G. F. R. Caceres Vera, T. Cai, M. F. Carneiro, D. Coplowe, H. da Motta, S. A. Dytman, G. A. Diaz, J. Felix, L. Fields, R. Fine, H. Gallagher, A. Ghosh , et al. (46 additional authors not shown)

Abstract: We report on multinucleon effects in low momentum transfer ($< 0.8$ GeV/c) anti-neutrino interactions on plastic (CH) scintillator. These data are from the 2010-2011 antineutrino phase of the MINERvA experiment at Fermilab. The hadronic energy spectrum of this inclusive sample is well described when a screening effect at low energy transfer and a two-nucleon knockout process are added to a relativ… ▽ More We report on multinucleon effects in low momentum transfer ($< 0.8$ GeV/c) anti-neutrino interactions on plastic (CH) scintillator. These data are from the 2010-2011 antineutrino phase of the MINERvA experiment at Fermilab. The hadronic energy spectrum of this inclusive sample is well described when a screening effect at low energy transfer and a two-nucleon knockout process are added to a relativistic Fermi gas model of quasielastic, $Δ$ resonance, and higher resonance processes. In this analysis, model elements introduced to describe previously published neutrino results have quantitatively similar benefits for this antineutrino sample. We present the results as a double-differential cross section to accelerate investigation of alternate models for antineutrino scattering off nuclei. △ Less

Submitted 15 January, 2019; v1 submitted 25 March, 2018; originally announced March 2018.

Comments: 7 pages, 4 figures. Updated to synchronize with reviewers and proofs recommendations and add supplemental cross section data file

Report number: FERMILAB-PUB-18-083-ND

Journal ref: Phys. Rev. Lett. 120, 221805 (2018)

Calculating the Magnetic Anisotropy of Rare-Earth-Transition-Metal Ferrimagnets

Authors: Christopher E. Patrick, Santosh Kumar, Geetha Balakrishnan, Rachel S. Edwards, Martin R. Lees, Leon Petit, Julie B. Staunton

Abstract: Magnetocrystalline anisotropy, the microscopic origin of permanent magnetism, is often explained in terms of ferromagnets. However, the best performing permanent magnets based on rare earths and transition metals (RE-TM) are in fact ferrimagnets, consisting of a number of magnetic sublattices. Here we show how a naive calculation of the magnetocrystalline anisotropy of the classic RE-TM ferrimagne… ▽ More Magnetocrystalline anisotropy, the microscopic origin of permanent magnetism, is often explained in terms of ferromagnets. However, the best performing permanent magnets based on rare earths and transition metals (RE-TM) are in fact ferrimagnets, consisting of a number of magnetic sublattices. Here we show how a naive calculation of the magnetocrystalline anisotropy of the classic RE-TM ferrimagnet GdCo$_5$ gives numbers which are too large at 0 K and exhibit the wrong temperature dependence. We solve this problem by introducing a first-principles approach to calculate temperature-dependent magnetization vs. field (FPMVB) curves, mirroring the experiments actually used to determine the anisotropy. We pair our calculations with measurements on a recently-grown single crystal of GdCo$_5$, and find excellent agreement. The FPMVB approach demonstrates a new level of sophistication in the use of first-principles calculations to understand RE-TM magnets. △ Less

Submitted 1 March, 2018; originally announced March 2018.

Comments: 13 pages, 3 figures; Supplemental Information included in submission

Journal ref: Phys. Rev. Lett. 120, 097202 (2018)

Measurement of the muon anti-neutrino double-differential cross section for quasi-elastic scattering on hydrocarbon at~$E_ν\sim 3.5$ GeV

Authors: C. E. Patrick, L. Aliaga, A. Bashyal, L. Bellantoni, A. Bercellie, M. Betancourt, A. Bodek, A. Bravar, H. Budd, G. F. R. Caceres, M. F. Carneiro, E. Chavarria, H. da Motta, S. A. Dytman, G. A. Diaz, J. Felix, L. Fields, R. Fine, A. M. Gago, R. Galindo, H. Gallager, A. Ghosh, R. Gran, J. Y. Han, D. A. Harris , et al. (42 additional authors not shown)

Abstract: We present double-differential measurements of anti-neutrino quasi-elastic scattering in the MINERvA detector. This study improves on a previous single differential measurement by using updated reconstruction algorithms and interaction models, and provides a complete description of observed muon kinematics in the form of a double-differential cross section with respect to muon transverse and longi… ▽ More We present double-differential measurements of anti-neutrino quasi-elastic scattering in the MINERvA detector. This study improves on a previous single differential measurement by using updated reconstruction algorithms and interaction models, and provides a complete description of observed muon kinematics in the form of a double-differential cross section with respect to muon transverse and longitudinal momentum. We include in our signal definition zero-meson final states arising from multi-nucleon interactions and from resonant pion production followed by pion absorption in the primary nucleus. We find that model agreement is considerably improved by a model tuned to MINERvA inclusive neutrino scattering data that incorporates nuclear effects such as weak nuclear screening and two-particle, two-hole enhancements. △ Less

Submitted 9 May, 2018; v1 submitted 3 January, 2018; originally announced January 2018.

Comments: 47 pages, 31 figures

Journal ref: Phys. Rev. D 97, 052002 (2018)

Simple vertex correction improves GW band energies of bulk and two-dimensional crystals

Authors: Per S. Schmidt, Christopher E. Patrick, Kristian S. Thygesen

Abstract: The GW self-energy method has long been recognized as the gold standard for quasiparticle (QP) calculations of solids in spite of the fact that the neglect of vertex corrections and the use of a DFT starting point lacks rigorous justification. In this work we remedy this situation by including a simple vertex correction that is consistent with an LDA starting point. We analyse the effect of the se… ▽ More The GW self-energy method has long been recognized as the gold standard for quasiparticle (QP) calculations of solids in spite of the fact that the neglect of vertex corrections and the use of a DFT starting point lacks rigorous justification. In this work we remedy this situation by including a simple vertex correction that is consistent with an LDA starting point. We analyse the effect of the self-energy by splitting it into a short-range and long-range term which are shown to govern respectively the center and size of the band gap. The vertex mainly improves the short-range correlations and therefore has a small effect on the band gap, while it shifts the band gap center up in energy by around 0.5 eV in good agreement with experiments. Our analysis also explains how the relative importance of short- and long-range interactions in structures of different dimensionality is reflected in their QP energies. Inclusion of the vertex comes at practically no extra computational cost and even improves the basis set convergence compared to GW. The method thus provides an efficient and rigorous improvement over the GW approximation and sets a new standard for quasiparticle calculations of solids. △ Less

Submitted 8 November, 2017; originally announced November 2017.

Journal ref: Phys. Rev. B 96, 205206 (2017)

Measurement of $ν_μ$ charged-current single $π^{0}$ production on hydrocarbon in the few-GeV region using MINERvA

Authors: O. Altinok, T. Le, L. Aliaga, L. Bellantoni, A. Bercellie, M. Betancourt, A. Bodek, A. Bravar, H. Budd, G. F. R. Caceres Vera, T. Cai, M. F. Carneiro, H. da Motta, S. A. Dytman, G. A. Díaz, J. Felix, L. Fields, R. Fine, A. M. Gago, R. Galindo, H. Gallagher, A. Ghosh, R. Gran, J. Y. Han, D. A. Harris , et al. (35 additional authors not shown)

Abstract: The semi-exclusive channel $ν_μ+\textrm{CH}\rightarrowμ^{-}π^{0}+\textrm{nucleon(s)}$ is analyzed using MINERvA exposed to the low-energy NuMI $ν_μ$ beam with spectral peak at $E_ν \simeq 3$ GeV. Differential cross sections for muon momentum and production angle, $π^{0}$ kinetic energy and production angle, and for squared four-momentum transfer are reported, and the cross section $σ(E_ν)$ is obta… ▽ More The semi-exclusive channel $ν_μ+\textrm{CH}\rightarrowμ^{-}π^{0}+\textrm{nucleon(s)}$ is analyzed using MINERvA exposed to the low-energy NuMI $ν_μ$ beam with spectral peak at $E_ν \simeq 3$ GeV. Differential cross sections for muon momentum and production angle, $π^{0}$ kinetic energy and production angle, and for squared four-momentum transfer are reported, and the cross section $σ(E_ν)$ is obtained over the range 1.5 GeV $\leq E_ν <$ 20 GeV. Results are compared to GENIE and NuWro predictions and to published MINERvA cross sections for $ν_μ\textrm{-CC}(π^{+})$ and $\barν_μ\textrm{-CC}(π^{0})$. Disagreements between data and simulation are observed at very low and relatively high values for muon angle and for $Q^2$ that may reflect shortfalls in modeling of interactions on carbon. For $π^{0}$ kinematic distributions however, the data are consistent with the simulation and provide support for generator treatments of pion intranuclear scattering. Using signal-event subsamples that have reconstructed protons as well as $π^{0}$ mesons, the $pπ^{0}$ invariant mass distribution is obtained, and the decay polar and azimuthal angle distributions in the rest frame of the $pπ^{0}$ system are measured in the region of $Δ(1232)^+$ production, $W < 1.4$ GeV. △ Less

Submitted 1 October, 2019; v1 submitted 11 August, 2017; originally announced August 2017.

Comments: 19 pages, 22 figures

Report number: FERMILAB-PUB-17-220-ND

Journal ref: Phys. Rev. D 96, 072003 (2017)

Rare-earth/transition-metal magnetic interactions in pristine and (Ni,Fe)-doped YCo5 and GdCo5

Authors: Christopher E. Patrick, Santosh Kumar, Geetha Balakrishnan, Rachel S. Edwards, Martin R. Lees, Eduardo Mendive-Tapia, Leon Petit, Julie B. Staunton

Abstract: We present an investigation into the intrinsic magnetic properties of the compounds YCo5 and GdCo5, members of the RETM5 class of permanent magnets (RE = rare earth, TM = transition metal). Focusing on Y and Gd provides direct insight into both the TM magnetization and RE-TM interactions without the complication of strong crystal field effects. We synthesize single crystals of YCo5 and GdCo5 using… ▽ More We present an investigation into the intrinsic magnetic properties of the compounds YCo5 and GdCo5, members of the RETM5 class of permanent magnets (RE = rare earth, TM = transition metal). Focusing on Y and Gd provides direct insight into both the TM magnetization and RE-TM interactions without the complication of strong crystal field effects. We synthesize single crystals of YCo5 and GdCo5 using the optical floating zone technique and measure the magnetization from liquid helium temperatures up to 800 K. These measurements are interpreted through calculations based on a Green's function formulation of density-functional theory, treating the thermal disorder of the local magnetic moments within the coherent potential approximation. The rise in the magnetization of GdCo5 with temperature is shown to arise from a faster disordering of the Gd magnetic moments compared to the antiferromagnetically aligned Co sublattice. We use the calculations to analyze the different Curie temperatures of the compounds and also compare the molecular (Weiss) fields at the RE site with previously published neutron scattering experiments. To gain further insight into the RE-TM interactions, we perform substitutional doping on the TM site, studying the compounds RECo4.5Ni0.5, RECo4Ni, and RECo4.5Fe0.5. Both our calculations and experiments on powdered samples find an increased/decreased magnetization with Fe/Ni doping, respectively. The calculations further reveal a pronounced dependence on the location of the dopant atoms of both the Curie temperatures and the Weiss field at the RE site. △ Less

Submitted 1 August, 2017; originally announced August 2017.

Comments: 13 figures, 1 table

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

Direct Measurement of Nuclear Dependence of Charged Current Quasielastic-like Neutrino Interactions using MINERvA

Authors: M. Betancourt, A. Ghosh, T. Walton, O. Altinok, L. Bellantoni, A. Bercellie, A. Bodek, A. Bravar, T. Cai, D. A. Martinez Caicedo, M. F. Carneiro, S. A. Dytman, G. A. Dıaz, J. Felix, L. Fields, R. Fine, R. Galindo, H. Gallagher, A. Ghosh, T. Golan, R. Gran, D. A. Harris, A. Higuera, K. Hurtado, M. Kiveni , et al. (32 additional authors not shown)

Abstract: Charged-current $ν_μ$ interactions on carbon, iron, and lead with a final state hadronic system of one or more protons with zero mesons are used to investigate the influence of the nuclear environment on quasielastic-like interactions. The transfered four-momentum squared to the target nucleus, $Q^2$, is reconstructed based on the kinematics of the leading proton, and differential cross sections v… ▽ More Charged-current $ν_μ$ interactions on carbon, iron, and lead with a final state hadronic system of one or more protons with zero mesons are used to investigate the influence of the nuclear environment on quasielastic-like interactions. The transfered four-momentum squared to the target nucleus, $Q^2$, is reconstructed based on the kinematics of the leading proton, and differential cross sections versus $Q^2$ and the cross-section ratios of iron, lead and carbon to scintillator are measured for the first time in a single experiment. The measurements show a dependence on atomic number. While the quasielastic-like scattering on carbon is compatible with predictions, the trends exhibited by scattering on iron and lead favor a prediction with intranuclear rescattering of hadrons accounted for by a conventional particle cascade treatment. These measurements help discriminate between different models of both initial state nucleons and final state interactions used in the neutrino oscillation experiments. △ Less

Submitted 10 May, 2017; originally announced May 2017.

Report number: FERMILAB-PUB-17-154-ND

Journal ref: Phys. Rev. Lett. 119, 082001 (2017)

Measurement of the antineutrino to neutrino charged-current interaction cross section ratio in MINERvA

Authors: Minerva Collaboration, L. Ren, L. Aliaga, O. Altinok, L. Bellantoni, A. Bercellie, M. Betancourt, A. Bodek, A. Bravar, H. Budd, T. Cai, M. F. Carneiro, H. da Motta, S. A. Dytman, G. A. Diaz, B. Eberly, E. Endress, J. Felix, L. Fields, R. Fine, A. M. Gago, R. Galindo, H. Gallagher, A. Ghosh, T. Golan , et al. (44 additional authors not shown)

Abstract: We present measurements of the neutrino and antineutrino total charged-current cross sections on carbon and their ratio using the MINERvA scintillator-tracker. The measurements span the energy range 2-22 GeV and were performed using forward and reversed horn focusing modes of the Fermilab low-energy NuMI beam to obtain large neutrino and antineutrino samples. The flux is obtained using a sub-sampl… ▽ More We present measurements of the neutrino and antineutrino total charged-current cross sections on carbon and their ratio using the MINERvA scintillator-tracker. The measurements span the energy range 2-22 GeV and were performed using forward and reversed horn focusing modes of the Fermilab low-energy NuMI beam to obtain large neutrino and antineutrino samples. The flux is obtained using a sub-sample of charged-current events at low hadronic energy transfer along with precise higher energy external neutrino cross section data overlapping with our energy range between 12-22 GeV. We also report on the antineutrino-neutrino cross section ratio, Rcc, which does not rely on external normalization information. Our ratio measurement, obtained within the same experiment using the same technique, benefits from the cancellation of common sample systematic uncertainties and reaches a precision of 5% at low energy. Our results for the antineutrino-nucleus scattering cross section and for Rcc are the most precise to date in the energy range $E_ν <$ 6GeV. △ Less

Submitted 1 January, 2018; v1 submitted 17 January, 2017; originally announced January 2017.

Journal ref: Phys. Rev. D 95, 072009 (2017)

Measurement of neutral-current $K^+$ production by neutrinos using MINERvA

Authors: C. M. Marshall, L. Aliaga, O. Altinok, L. Bellantoni, A. Bercellie, M. Betancourt, A. Bodek, A. Bravar, T. Cai, M. F. Carneiro, H. da Motta, S. A. Dytman, G. A. Dìaz, M. Dunkman, B. Eberly, E. Endress, J. Felix, L. Fields, R. Fine, A. M. Gago, R. Galindo, H. Gallagher, A. Ghosh, T. Golan, R. Gran , et al. (42 additional authors not shown)

Abstract: Neutral-current production of $K^{+}$ by atmospheric neutrinos is a background in searches for the proton decay $p \rightarrow K^{+} \barν$. Reactions such as $νp \rightarrow νK^{+} Λ$ are indistinguishable from proton decays when the decay products of the $Λ$ are below detection threshold. Events with $K^{+}$ are identified in MINERvA by reconstructing the timing signature of a $K^{+}$ decay at r… ▽ More Neutral-current production of $K^{+}$ by atmospheric neutrinos is a background in searches for the proton decay $p \rightarrow K^{+} \barν$. Reactions such as $νp \rightarrow νK^{+} Λ$ are indistinguishable from proton decays when the decay products of the $Λ$ are below detection threshold. Events with $K^{+}$ are identified in MINERvA by reconstructing the timing signature of a $K^{+}$ decay at rest. A sample of 201 neutrino-induced neutral-current $K^{+}$ events is used to measure differential cross sections with respect to the $K^{+}$ kinetic energy, and the non-$K^{+}$ hadronic visible energy. An excess of events at low hadronic visible energy is observed relative to the prediction of the NEUT event generator. Good agreement is observed with the cross section prediction of the GENIE generator. A search for photons from $π^{0}$ decay, which would veto a neutral-current $K^{+}$ event in a proton decay search, is performed, and a 2$σ$ deficit of detached photons is observed relative to the GENIE prediction. △ Less

Submitted 2 June, 2017; v1 submitted 7 November, 2016; originally announced November 2016.

Journal ref: Phys. Rev. Lett. 119, 011802 (2017)

Measurements of the Inclusive Neutrino and Antineutrino Charged Current Cross Sections in MINERvA Using the Low-$ν$ Flux Method

Authors: MINERvA Collaboration, J. DeVan, L. Ren, L. Aliaga, O. Altinok, L. Bellantoni, A. Bercellie, M. Betancourt, A. Bodek, H. Budd, T. Cai, M. F. Carneiro, H. da Motta, J. Devan, S. A. Dytman, G. A. Díaz, B. Eberly, E. Endress, J. Felix, L. Fields, R. Fine, A. M. Gago, R. Galindo, H. Gallagher, A. Ghosh , et al. (43 additional authors not shown)

Abstract: The total cross sections are important ingredients for the current and future neutrino oscillation experiments. We present measurements of the total charged-current neutrino and antineutrino cross sections on scintillator (CH) in the NuMI low-energy beamline using an {\em in situ} prediction of the shape of the flux as a function of neutrino energy from 2--50 GeV. This flux prediction takes advant… ▽ More The total cross sections are important ingredients for the current and future neutrino oscillation experiments. We present measurements of the total charged-current neutrino and antineutrino cross sections on scintillator (CH) in the NuMI low-energy beamline using an {\em in situ} prediction of the shape of the flux as a function of neutrino energy from 2--50 GeV. This flux prediction takes advantage of the fact that neutrino and antineutrino interactions with low nuclear recoil energy ($ν$) have a nearly constant cross section as a function of incident neutrino energy. This measurement is the lowest energy application of the low-$ν$ flux technique, the first time it has been used in the NuMI antineutrino beam configuration, and demonstrates that the technique is applicable to future neutrino beams operating at multi-GeV energies. The cross section measurements presented are the most precise measurements to date below 5 GeV. △ Less

Submitted 15 October, 2016; originally announced October 2016.

Comments: 8 pages, 6 figures, and 12 supplemental tables of fluxes, cross sections, breakdowns of their uncertainties, and their covariance matrices

Report number: FERMILAB-PUB-16-368-ND

Journal ref: Phys. Rev. D 94, 112007 (2016)

Frustrated Magnetism and Caloric Effects in Mn-antiperovskite Nitrides: $Ab~Initio$ Theory

Authors: J. Zemen, E. Mendive Tapia, Z. Gercsi, R. Banerjee, C. E. Patrick, J. B. Staunton, K. G. Sandeman

Abstract: We model changes of magnetic ordering in Mn-antiperovskite nitrides driven by biaxial lattice strain at zero and at finite temperature. We employ a non-collinear spin-polarised density functional theory to compare the response of the geometrically frustrated exchange interactions to a tetragonal symmetry breaking (the so called piezomagnetic effect) across a range of Mn$_3$AN (A = Rh, Pd, Ag, Co,… ▽ More We model changes of magnetic ordering in Mn-antiperovskite nitrides driven by biaxial lattice strain at zero and at finite temperature. We employ a non-collinear spin-polarised density functional theory to compare the response of the geometrically frustrated exchange interactions to a tetragonal symmetry breaking (the so called piezomagnetic effect) across a range of Mn$_3$AN (A = Rh, Pd, Ag, Co, Ni, Zn, Ga, In, Sn) at zero temperature. Building on the robustness of the effect we focus on Mn$_3$GaN and extend our study to finite temperature using the disordered local moment (DLM) first-principles electronic structure theory to model the interplay between the ordering of Mn magnetic moments and itinerant electron states. We discover a rich temperature-strain magnetic phase diagram with two previously unreported phases stabilised by strains larger than 0.75\% and with transition temperatures strongly dependent on strain. We propose an elastocaloric cooling cycle crossing two of the available phase transitions to achieve simultaneously a large isothermal entropy change (due to the first order transition) and a large adiabatic temperature change (due to the second order transition). △ Less

Submitted 19 April, 2017; v1 submitted 12 September, 2016; originally announced September 2016.

Comments: 12 pages, 10 figures

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

Quantitative analysis of valence photoemission spectra and quasiparticle excitations at chromophore-semiconductor interfaces

Authors: Christopher E. Patrick, Feliciano Giustino

Abstract: Investigating quasiparticle excitations of molecules on surfaces through photoemission spectroscopy forms a major part of nanotechnology research. Resolving spectral features at these interfaces requires a comprehensive theory of electron removal and addition processes in molecules and solids which captures the complex interplay of image charges, thermal effects and configurational disorder. We he… ▽ More Investigating quasiparticle excitations of molecules on surfaces through photoemission spectroscopy forms a major part of nanotechnology research. Resolving spectral features at these interfaces requires a comprehensive theory of electron removal and addition processes in molecules and solids which captures the complex interplay of image charges, thermal effects and configurational disorder. We here develop such a theory and calculate the quasiparticle energy-level alignment and the valence photoemission spectrum for the prototype biomimetic solar cell interface between anatase TiO2 and the N3 chromophore. By directly matching our calculated photoemission spectrum to experimental data we clarify the atomistic origin of the chromophore peak at low binding energy. This case study sets a new standard in the interpretation of photoemission spectroscopy at complex chromophore/semiconductor interfaces. △ Less

Submitted 27 July, 2016; originally announced July 2016.

Comments: 4 figures. Article published by in Phys. Rev. Lett. by the American Physical Society (2012)

Journal ref: Phys. Rev. Lett. 109, 116801 (2012)

Structure of a Water Monolayer on the Anatase TiO$_2$(101) Surface

Authors: Christopher E. Patrick, Feliciano Giustino

Abstract: Titanium dioxide (TiO$_2$) plays a central role in the study of artificial photosynthesis, owing to its ability to perform photocatalytic water splitting. Despite over four decades of intense research efforts in this area, there is still some debate over the nature of the first water monolayer on the technologically-relevant anatase TiO$_2$ (101) surface. In this work we use first-principles calcu… ▽ More Titanium dioxide (TiO$_2$) plays a central role in the study of artificial photosynthesis, owing to its ability to perform photocatalytic water splitting. Despite over four decades of intense research efforts in this area, there is still some debate over the nature of the first water monolayer on the technologically-relevant anatase TiO$_2$ (101) surface. In this work we use first-principles calculations to reverse-engineer the experimental high-resolution X-ray photoelectron spectra measured for this surface in [Walle et al., J. Phys. Chem. C 115, 9545 (2011)], and find evidence supporting the existence of a mix of dissociated and molecular water in the first monolayer. Using both semilocal and hybrid functional calculations we revise the current understanding of the adsorption energetics by showing that the energetic cost of water dissociation is reduced via the formation of a hydrogen-bonded hydroxyl-water complex. We also show that such a complex can provide an explanation of an unusual superstructure observed in high-resolution scanning tunneling microscopy experiments. △ Less

Submitted 27 July, 2016; originally announced July 2016.

Comments: 11 pages, 9 figures. Published in Phys. Rev. Appl. by the American Physical Society (2014)

Journal ref: Phys. Rev. Applied 2, 014001 (2014)

Unified theory of electron-phonon renormalization and phonon-assisted optical absorption

Authors: Christopher E. Patrick, Feliciano Giustino

Abstract: We present a theory of electronic excitation energies and optical absorption spectra which incorporates energy-level renormalization and phonon-assisted optical absorption within a unified framework. Using time-independent perturbation theory we show how the standard approaches for studying vibronic effects in molecules and those for addressing electron-phonon interactions in solids correspond to… ▽ More We present a theory of electronic excitation energies and optical absorption spectra which incorporates energy-level renormalization and phonon-assisted optical absorption within a unified framework. Using time-independent perturbation theory we show how the standard approaches for studying vibronic effects in molecules and those for addressing electron-phonon interactions in solids correspond to slightly different choices for the non-interacting Hamiltonian. Our present approach naturally leads to the Allen-Heine theory of temperature-dependent energy levels, the Franck-Condon principle, the Herzberg-Teller effect and to phonon-assisted optical absorption in indirect band gap materials. In addition, our theory predicts sub-gap phonon-assisted optical absorption in direct gap materials, as well as an exponential edge which we tentatively assign to the Urbach tail. We also consider a semiclassical approach to the calculation of optical absorption spectra which simultaneously captures energy-level renormalization and phonon-assisted transitions and is especially suited to first-principles electronic structure calculations. We demonstrate this approach by calculating the phonon-assisted optical absorption spectrum of bulk silicon. △ Less

Submitted 27 July, 2016; originally announced July 2016.

Comments: This is an author-created, un-copyedited version of an article published in J. Phys.: Condens. Matter. IOP Publishing Ltd is not responsible for any errors or omissions in this version of the manuscript or any version derived from it. The Version of Record is available online at http://dx.doi.org/10.1088/0953-8984/26/36/365503

Journal ref: J. Phys.: Condens. Matter 26, 365503 (2014)

Neutrino Flux Predictions for the NuMI Beam

Authors: MINERvA Collaboration, L. Aliaga, M. Kordosky, T. Golan, O. Altinok, L. Bellantoni, A. Bercellie, M. Betancourt, A. Bravar, H. Budd, M. F. Carneiro, G. A. Diaz, E. Endress, J. Felix, L. Fields, R. Fine, A. M. Gago, R. Galindo, H. Gallagher, R. Gran, D. A. Harris, A. Higuera, K. Hurtado, M. Kiveni, J. Kleykamp , et al. (36 additional authors not shown)

Abstract: Knowledge of the neutrino flux produced by the Neutrinos at the Main Injector (NuMI) beamline is essential to the neutrino oscillation and neutrino interaction measurements of the MINERvA, MINOS+, NOvA and MicroBooNE experiments at Fermi National Accelerator Laboratory. We have produced a flux prediction which uses all available and relevant hadron production data, incorporating measurements of pa… ▽ More Knowledge of the neutrino flux produced by the Neutrinos at the Main Injector (NuMI) beamline is essential to the neutrino oscillation and neutrino interaction measurements of the MINERvA, MINOS+, NOvA and MicroBooNE experiments at Fermi National Accelerator Laboratory. We have produced a flux prediction which uses all available and relevant hadron production data, incorporating measurements of particle production off of thin targets as well as measurements of particle yields from a spare NuMI target exposed to a 120 GeV proton beam. The result is the most precise flux prediction achieved for a neutrino beam in the one to tens of GeV energy region. We have also compared the prediction to in situ measurements of the neutrino flux and find good agreement. △ Less

Submitted 11 July, 2016; v1 submitted 3 July, 2016; originally announced July 2016.

Comments: v2 includes supplemental material consisting of the flux prediction and uncertainties in ascii and root format, a program to read the ascii datafiles, and a short guide with additional details. v3 fixes a couple of typographical errors on the units for some quantities in the beam focusing section

Report number: Fermilab PUB-16-091-ND

Journal ref: Phys. Rev. D 94, 092005 (2016)

First evidence of coherent $K^{+}$ meson production in neutrino-nucleus scattering

Authors: Z. Wang, C. M. Marshall, L. Aliaga, O. Altinok, L. Bellantoni, A. Bercellie, M. Betancourt, A. Bodek, A. Bravar, H. Budd, T. Cai, M. F. Carneiro, H. da Motta, S. A. Dytman, G. A. Dìaz, B. Eberly, E. Endress, J. Felix, L. Fields, R. Fine, R. Galindo, H. Gallagher, A. Ghosh, T. Golan, R. Gran , et al. (46 additional authors not shown)

Abstract: Neutrino-induced charged-current coherent kaon production, $ν_μA\rightarrowμ^{-}K^{+}A$, is a rare, inelastic electroweak process that brings a $K^+$ on shell and leaves the target nucleus intact in its ground state. This process is significantly lower in rate than neutrino-induced charged-current coherent pion production, because of Cabibbo suppression and a kinematic suppression due to the large… ▽ More Neutrino-induced charged-current coherent kaon production, $ν_μA\rightarrowμ^{-}K^{+}A$, is a rare, inelastic electroweak process that brings a $K^+$ on shell and leaves the target nucleus intact in its ground state. This process is significantly lower in rate than neutrino-induced charged-current coherent pion production, because of Cabibbo suppression and a kinematic suppression due to the larger kaon mass. We search for such events in the scintillator tracker of MINERvA by observing the final state $K^+$, $μ^-$ and no other detector activity, and by using the kinematics of the final state particles to reconstruct the small momentum transfer to the nucleus, which is a model-independent characteristic of coherent scattering. We find the first experimental evidence for the process at $3σ$ significance. △ Less

Submitted 12 July, 2016; v1 submitted 28 June, 2016; originally announced June 2016.

Comments: added ancillary file with information about the six kaon candidates

Report number: FERMILAB-PUB-16-282-ND

Journal ref: Phys. Rev. Lett. 117, 061802 (2016)

Cross sections for neutrino and antineutrino induced pion production on hydrocarbon in the few-GeV region using MINERvA

Authors: C. L. McGivern, T. Le, B. Eberly, L. Aliaga, O. Altinok, L. Bellantoni, A. Bercellie, M. Betancourt, A. Bodek, A. Bravar, H. Budd, T. Cai, M. F. Carneiro, M. E. Christy, H. da Motta, S. A. Dytman, G. A. Diaz, E. Endress, J. Felix, L. Fields, R. Fine, R. Galindo, H. Gallagher, T. Golan, R. Gran , et al. (45 additional authors not shown)

Abstract: Separate samples of charged-current pion production events representing two semi-inclusive channels $ν_μ$-CC($π^{+}$) and $\barν_μ$-CC($π^{0}$) have been obtained using neutrino and antineutrino exposures of the MINERvA detector. Distributions in kinematic variables based upon $μ^{\pm}$-track reconstructions are analyzed and compared for the two samples. The differential cross sections for muon pr… ▽ More Separate samples of charged-current pion production events representing two semi-inclusive channels $ν_μ$-CC($π^{+}$) and $\barν_μ$-CC($π^{0}$) have been obtained using neutrino and antineutrino exposures of the MINERvA detector. Distributions in kinematic variables based upon $μ^{\pm}$-track reconstructions are analyzed and compared for the two samples. The differential cross sections for muon production angle, muon momentum, and four-momentum transfer $Q^2$, are reported, and cross sections versus neutrino energy are obtained. Comparisons with predictions of current neutrino event generators are used to clarify the role of the $Δ(1232)$ and higher-mass baryon resonances in CC pion production and to show the importance of pion final-state interactions. For the $ν_μ$-CC($π^{+}$) ($\barν_μ$-CC($π^{0}$)) sample, the absolute data rate is observed to lie below (above) the predictions of some of the event generators by amounts that are typically 1-to-2 $σ$. However the generators are able to reproduce the shapes of the differential cross sections for all kinematic variables of either data set. △ Less

Submitted 8 October, 2018; v1 submitted 22 June, 2016; originally announced June 2016.

Comments: 17 pages, 16 figures; supplementary materials in Supplement file (TeX and pdf)

Journal ref: Phys. Rev. D 94, 052005 (2016)

Measurement of $K^{+}$ production in charged-current $ν_μ$ interactions

Authors: C. M. Marshall, L. Aliaga, O. Altinok, L. Bellantoni, A. Bercellie, M. Betancourt, A. Bodek, A. Bravar, H. Budd, T. Cai, M. F. Carneiro, J. Chvojka, H. da Motta, J. Devan, S. A. Dytman, G. A. Díaz, B. Eberly, E. Endress, J. Felix, L. Fields, A. Filkins, R. Fine, A. M. Gago, R. Galindo, H. Gallagher , et al. (57 additional authors not shown)

Abstract: Production of K^{+} mesons in charged-current ν_μ interactions on plastic scintillator (CH) is measured using MINERvA exposed to the low-energy NuMI beam at Fermilab. Timing information is used to isolate a sample of 885 charged-current events containing a stopping K^{+} which decays at rest. The differential cross section in K^{+} kinetic energy, dσ/dT_{K}, is observed to be relatively flat betwe… ▽ More Production of K^{+} mesons in charged-current ν_μ interactions on plastic scintillator (CH) is measured using MINERvA exposed to the low-energy NuMI beam at Fermilab. Timing information is used to isolate a sample of 885 charged-current events containing a stopping K^{+} which decays at rest. The differential cross section in K^{+} kinetic energy, dσ/dT_{K}, is observed to be relatively flat between 0 and 500 MeV. Its shape is in good agreement with the prediction by the \textsc{genie} neutrino event generator when final-state interactions are included, however the data rate is lower than the prediction by 15\%. △ Less

Submitted 25 July, 2016; v1 submitted 13 April, 2016; originally announced April 2016.

Comments: added ancillary files with cross-section, statistical uncertainty covariance matrix and systematic uncertainty covariance matrix decomposed into flux and non-flux components

Report number: FERMILAB PUB-16-109-ND

Journal ref: Phys. Rev. D 94, 012002 (2016)

Evidence for neutral-current diffractive neutral pion production from hydrogen in neutrino interactions on hydrocarbon

Authors: MINERvA Collaboration, J. Wolcott, L. Aliaga, O. Altinok, A. Bercellie, M. Betancourt, A. Bodek, A. Bravar, H. Budd, T. Cai, M. F. Carneiro, J. Chvojka, H. da Motta, J. Devan, S. A. Dytman, G. A. Díaz, B. Eberly, E. Endress, J. Felix, L. Fields, R. Fine, R. Galindo, H. Gallagher, T. Golan, R. Gran , et al. (46 additional authors not shown)

Abstract: The MINERvA experiment observes an excess of events containing electromagnetic showers relative to the expectation from Monte Carlo simulations in neutral-current neutrino interactions with mean beam energy of 4.5 GeV on a hydrocarbon target. The excess is characterized and found to be consistent with neutral-current neutral pion production with a broad energy distribution peaking at 7 GeV and a t… ▽ More The MINERvA experiment observes an excess of events containing electromagnetic showers relative to the expectation from Monte Carlo simulations in neutral-current neutrino interactions with mean beam energy of 4.5 GeV on a hydrocarbon target. The excess is characterized and found to be consistent with neutral-current neutral pion production with a broad energy distribution peaking at 7 GeV and a total cross section of 0.26 +- 0.02 (stat) +- 0.08 (sys) x 10^{-39} cm^{2}. The angular distribution, electromagnetic shower energy, and spatial distribution of the energy depositions of the excess are consistent with expectations from neutrino neutral-current diffractive neutral pion production from hydrogen in the hydrocarbon target. These data comprise the first direct experimental observation and constraint for a reaction that poses an important background process in neutrino oscillation experiments searching for muon neutrino to electron neutrino oscillations. △ Less

Submitted 28 July, 2016; v1 submitted 6 April, 2016; originally announced April 2016.

Comments: 15 pages, 7 figures; accepted by Phys. Rev. Lett

Report number: FERMILAB-PUB-16-108-ND

Journal ref: Phys. Rev. Lett. 117, 111801 (2016)

Hubbard-$U$ corrected Hamiltonians for non-self-consistent random-phase approximation total-energy calculations: A study of ZnS, TiO$_2$, and NiO

Authors: Christopher E. Patrick, Kristian S. Thygesen

Abstract: In non-self-consistent calculations of the total energy within the random-phase approximation (RPA) for electronic correlation, it is necessary to choose a single-particle Hamiltonian whose solutions are used to construct the electronic density and non-interacting response function. Here we investigate the effect of including a Hubbard-$U$ term in this single-particle Hamiltonian, to better descri… ▽ More In non-self-consistent calculations of the total energy within the random-phase approximation (RPA) for electronic correlation, it is necessary to choose a single-particle Hamiltonian whose solutions are used to construct the electronic density and non-interacting response function. Here we investigate the effect of including a Hubbard-$U$ term in this single-particle Hamiltonian, to better describe the on-site correlation of 3$d$ electrons in the transition metal compounds ZnS, TiO$_2$ and NiO. We find that the RPA lattice constants are essentially independent of $U$, despite large changes in the underlying electronic structure. We further demonstrate that the non-self-consistent RPA total energies of these materials have minima at nonzero $U$. Our RPA calculations find the rutile phase of TiO$_2$ to be more stable than anatase independent of $U$, a result which is consistent with experiments and qualitatively different to that found from calculations employing $U$-corrected (semi)local functionals. However we also find that the +$U$ term cannot be used to correct the RPA's poor description of the heat of formation of NiO. △ Less

Submitted 27 January, 2016; originally announced January 2016.

Comments: 31 pages, 10 figures

Journal ref: Phys. Rev. B 93, 035133 (2016)

Measurement of Partonic Nuclear Effects in Deep-Inelastic Neutrino Scattering using MINERvA

Authors: MINERvA Collaboration, J. Mousseau, M. Wospakrik, L. Aliaga, O. Altinok, L. Bellantoni, A. Bercellie, M. Betancourt, A. Bodek, A. Bravar, H. Budd, T. Cai, M. F. Carneiro, M. E. Christy, J. Chvojka, H. da Motta, J. Devan, S. A. Dytman, G. A. Díaz, B. Eberly, J. Felix, L. Fields, R. Fine, A. M. Gago, R. Galindo , et al. (50 additional authors not shown)

Abstract: The MINERvA collaboration reports a novel study of neutrino-nucleus charged-current deep inelastic scattering (DIS) using the same neutrino beam incident on targets of polystyrene, graphite, iron, and lead. Results are presented as ratios of C, Fe, and Pb to CH. The ratios of total DIS cross sections as a function of neutrino energy and flux-integrated differential cross sections as a function of… ▽ More The MINERvA collaboration reports a novel study of neutrino-nucleus charged-current deep inelastic scattering (DIS) using the same neutrino beam incident on targets of polystyrene, graphite, iron, and lead. Results are presented as ratios of C, Fe, and Pb to CH. The ratios of total DIS cross sections as a function of neutrino energy and flux-integrated differential cross sections as a function of the Bjorken scaling variable x are presented in the neutrino-energy range of 5 - 50 GeV. Good agreement is found between the data and predicted ratios, based on charged-lepton nucleus scattering, at medium x and low neutrino energies. However, the data rate appears depleted in the vicinity of the nuclear shadowing region, x < 0.1. This apparent deficit, reflected in the DIS cross-section ratio at high neutrino energy , is consistent with previous MINERvA observations and with the predicted onset of nuclear shadowing with the the axial-vector current in neutrino scattering. △ Less

Submitted 30 September, 2016; v1 submitted 23 January, 2016; originally announced January 2016.

Report number: Fermilab-Pub-15-556-ND

Journal ref: Phys. Rev. D 93, 071101 (2016)

Measurement of Neutrino Flux from Neutrino-Electron Elastic Scattering

Authors: MINERvA Collaboration, J. Park, L. Aliaga, O. Altinok, L. Bellantoni, A. Bercellie, M. Betancourt, A. Bodek, A. Bravar, H. Budd, T. Cai, M. F. Carneiro, M. E. Christy, J. Chvojka, H. da Motta, S. A. Dytman, G. A. Diaz, B. Eberly, J. Felix, L. Fields, R. Fine, A. M. Gago, R. Galindo, A. Ghosh, T. Golan , et al. (44 additional authors not shown)

Abstract: Muon-neutrino elastic scattering on electrons is an observable neutrino process whose cross section is precisely known. Consequently a measurement of this process in an accelerator-based $ν_μ$ beam can improve the knowledge of the absolute neutrino flux impinging upon the detector; typically this knowledge is limited to $\sim$ 10% due to uncertainties in hadron production and focusing. We have iso… ▽ More Muon-neutrino elastic scattering on electrons is an observable neutrino process whose cross section is precisely known. Consequently a measurement of this process in an accelerator-based $ν_μ$ beam can improve the knowledge of the absolute neutrino flux impinging upon the detector; typically this knowledge is limited to $\sim$ 10% due to uncertainties in hadron production and focusing. We have isolated a sample of 135 $\pm$ 17 neutrino-electron elastic scattering candidates in the segmented scintillator detector of MINERvA, after subtracting backgrounds and correcting for efficiency. We show how this sample can be used to reduce the total uncertainty on the NuMI $ν_μ$ flux from 9% to 6%. Our measurement provides a flux constraint that is useful to other experiments using the NuMI beam, and this technique is applicable to future neutrino beams operating at multi-GeV energies. △ Less

Submitted 15 June, 2016; v1 submitted 23 December, 2015; originally announced December 2015.

Comments: 11 pages, 11 figures

Report number: FERMILAB-PUB-15-575-ND

Journal ref: Phys. Rev. D 93, 112007 (2016)

Identification of nuclear effects in neutrino-carbon interactions at low three-momentum transfer

Authors: P. A. Rodrigues, J. Demgen, E. Miltenberger, L. Aliaga, O. Altinok, L. Bellantoni, A. Bercellie, M. Betancourt, A. Bodek, A. Bravar, H. Budd, T. Kai, M. F. Carneiro, J. Chvojka, J. Devan, S. A. Dytman, G. A. Diaz, B. Eberly, M. Elkins, J. Felix, L. Fields, R. Fine, A. M. Gago, R. Galindo, H. Gallagher , et al. (51 additional authors not shown)

Abstract: Two different nuclear-medium effects are isolated using a low three-momentum transfer subsample of neutrino-carbon scattering data from the MINERvA neutrino experiment. The observed hadronic energy in charged-current $ν_μ$ interactions is combined with muon kinematics to permit separation of the quasielastic and $Δ$(1232) resonance processes. First, we observe a small cross section at very low ene… ▽ More Two different nuclear-medium effects are isolated using a low three-momentum transfer subsample of neutrino-carbon scattering data from the MINERvA neutrino experiment. The observed hadronic energy in charged-current $ν_μ$ interactions is combined with muon kinematics to permit separation of the quasielastic and $Δ$(1232) resonance processes. First, we observe a small cross section at very low energy transfer that matches the expected screening effect of long-range nucleon correlations. Second, additions to the event rate in the kinematic region between the quasielastic and $Δ$ resonance processes are needed to describe the data. The data in this kinematic region also has an enhanced population of multi-proton final states. Contributions predicted for scattering from a nucleon pair have both properties; the model tested in this analysis is a significant improvement but does not fully describe the data. We present the results as a double-differential cross section to enable further investigation of nuclear models. Improved description of the effects of the nuclear environment are required by current and future neutrino oscillation experiments. △ Less

Submitted 20 September, 2019; v1 submitted 18 November, 2015; originally announced November 2015.

Comments: 7 pages, 3 figures, supplemental material. v2 includes updates from the review process. v3 includes ancillary files. v4 fixes incorrect exponent in Fig. 2 and supplemental files, synchronized with the erratum in the journal version

Journal ref: Phys. Rev. Lett. 116, 071802 (2016)

Stochastic approach to phonon-assisted optical absorption

Authors: Marios Zacharias, Christopher E. Patrick, Feliciano Giustino

Abstract: We develop a first-principles theory of phonon-assisted optical absorption in semiconductors and insulators which incorporates the temperature dependence of the electronic structure. We show that the Hall-Bardeen-Blatt theory of indirect optical absorption and the Allen-Heine theory of temperature-dependent band structures can be derived from the present formalism by retaining only one-phonon proc… ▽ More We develop a first-principles theory of phonon-assisted optical absorption in semiconductors and insulators which incorporates the temperature dependence of the electronic structure. We show that the Hall-Bardeen-Blatt theory of indirect optical absorption and the Allen-Heine theory of temperature-dependent band structures can be derived from the present formalism by retaining only one-phonon processes. We demonstrate this method by calculating the optical absorption coefficient of silicon using an importance sampling Monte Carlo scheme, and we obtain temperature-dependent lineshapes and band gaps in good agreement with experiment. The present approach opens the way to predictive calculations of the optical properties of solids at finite temperature. △ Less

Submitted 21 October, 2015; originally announced October 2015.

Comments: 11 pages, 9 figures

Journal ref: Phys. Rev. Lett. 115, 177401, 2015

Anharmonic stabilization and band gap renormalization in the perovskite CsSnI$_3$

Authors: Christopher E. Patrick, Karsten W. Jacobsen, Kristian S. Thygesen

Abstract: Amongst the X(Sn,Pb)Y$_3$ perovskites currently under scrutiny for their photovoltaic applications, the cubic B-$α$ phase of CsSnI$_3$ is arguably the best characterized experimentally. Yet, according to the standard harmonic theory of phonons, this deceptively simple phase should not exist at all due to rotational instabilities of the SnI$_6$ octahedra. Here, employing self-consistent phonon theo… ▽ More Amongst the X(Sn,Pb)Y$_3$ perovskites currently under scrutiny for their photovoltaic applications, the cubic B-$α$ phase of CsSnI$_3$ is arguably the best characterized experimentally. Yet, according to the standard harmonic theory of phonons, this deceptively simple phase should not exist at all due to rotational instabilities of the SnI$_6$ octahedra. Here, employing self-consistent phonon theory we show that these soft modes are stabilized at experimental conditions through anharmonic phonon-phonon interactions between the Cs ions and their iodine cages. We further calculate the renormalization of the electronic energies due to vibrations and find an unusual opening of the band gap, estimated as 0.24 and 0.11 eV at 500 and 300 K, which we attribute to the stretching of Sn-I bonds. Our work demonstrates the important role of temperature in accurately describing these materials. △ Less

Submitted 21 October, 2015; originally announced October 2015.

Comments: 12 pages, 4 figures

Journal ref: Phys. Rev. B 92, 201205 (2015)

Measurement of electron neutrino quasielastic and quasielastic-like scattering on hydrocarbon at $\langle E_ν \rangle $ = 3.6 GeV

Authors: MINERvA Collaboration, J. Wolcott, O. Altinok, L. Bellantoni, A. Bercellie, M. Betancourt, A. Bodek, A. Bravar, H. Budd, M. F. Carneiro, J. Chvojka, H. da Motta, J. Devan, S. A. Dytman, G. A. Diaz, B. Eberly, J. Felix, L. Fields, R. Fine, R. Galindo, H. Gallagher, A. Ghosh, T. Golan, R. Gran, D. A. Harris , et al. (45 additional authors not shown)

Abstract: The first direct measurement of electron-neutrino quasielastic and quasielastic-like scattering on hydrocarbon in the few-GeV region of incident neutrino energy has been carried out using the MINERvA detector in the NuMI beam at Fermilab. The flux-integrated differential cross sections in electron production angle, electron energy and $Q^{2}$ are presented. The ratio of the quasielastic, flux-inte… ▽ More The first direct measurement of electron-neutrino quasielastic and quasielastic-like scattering on hydrocarbon in the few-GeV region of incident neutrino energy has been carried out using the MINERvA detector in the NuMI beam at Fermilab. The flux-integrated differential cross sections in electron production angle, electron energy and $Q^{2}$ are presented. The ratio of the quasielastic, flux-integrated differential cross section in $Q^{2}$ for $ν_{e}$ with that of similarly-selected $ν_μ$-induced events from the same exposure is used to probe assumptions that underpin conventional treatments of charged-current $ν_{e}$ interactions used by long-baseline neutrino oscillation experiments. The data are found to be consistent with lepton universality and are well-described by the predictions of the neutrino event generator GENIE. △ Less

Submitted 9 November, 2016; v1 submitted 18 September, 2015; originally announced September 2015.

Comments: includes ancillary files with flux, measured cross sections, and covariance matrices (v9 includes tables of nu_e + antinu_e / nu_mu ratio and associated covariances)

Report number: FERMILAB-PUB-15-399-ND

Journal ref: Phys. Rev. Lett. 116, 081802 (2016)

Single neutral pion production by charged-current $\barν_μ$ interactions on hydrocarbon at $\langle E_ν\rangle = $ 3.6 GeV

Authors: T. Le, J. L. Palomino, L. Aliaga, O. Altinok, A. Bercellie, A. Bodek, A. Bravar, W. K. Brooks, A. Butkevich, D. A. Martinez Caicedo, M. F. Carneiro, M. E. Christy, J. Chvojka, H. da Motta, J. Devan, S. A. Dytman, G. A. Díaz, B. Eberly, J. Felix, L. Fields, R. Fine, A. M. Gago, H. Gallagher, R. Gran, D. A. Harris , et al. (40 additional authors not shown)

Abstract: Single neutral pion production via muon antineutrino charged-current interactions in plastic scintillator (CH) is studied using the \minerva detector exposed to the NuMI low-energy, wideband antineutrino beam at Fermilab. Measurement of this process constrains models of neutral pion production in nuclei, which is important because the neutral-current analog is a background for $\barν_e$ appearance… ▽ More Single neutral pion production via muon antineutrino charged-current interactions in plastic scintillator (CH) is studied using the \minerva detector exposed to the NuMI low-energy, wideband antineutrino beam at Fermilab. Measurement of this process constrains models of neutral pion production in nuclei, which is important because the neutral-current analog is a background for $\barν_e$ appearance oscillation experiments. The differential cross sections for $π^0$ momentum and production angle, for events with a single observed $π^0$ and no charged pions, are presented and compared to model predictions. These results comprise the first measurement of the $π^0$ kinematics for this process. △ Less

Submitted 24 August, 2015; v1 submitted 6 March, 2015; originally announced March 2015.

Comments: 6 pages, 5 figures, submitted to Physics Letters B

Report number: FERMILAB-PUB-15-075-ND

Journal ref: Phys.Lett. B749 (2015) 130-136

MINERvA neutrino detector response measured with test beam data

Authors: MINERvA Collaboration, L. Aliaga, O. Altinok, C. Araujo Del Castillo, L. Bagby, L. Bellantoni, W. F. Bergan, A. Bodek, R. Bradford, A. Bravar, H. Budd, A. Butkevich, D. A. Martinez Caicedo, M. F. Carneiro, M. E. Christy, J. Chvojka, H. da Motta, J. Devan, G. A. Diaz, S. A. Dytman, B. Eberly, J. Felix, L. Fields, R. Fine, R. Flight , et al. (63 additional authors not shown)

Abstract: The MINERvA collaboration operated a scaled-down replica of the solid scintillator tracking and sampling calorimeter regions of the MINERvA detector in a hadron test beam at the Fermilab Test Beam Facility. This article reports measurements with samples of protons, pions, and electrons from 0.35 to 2.0 GeV/c momentum. The calorimetric response to protons, pions, and electrons are obtained from the… ▽ More The MINERvA collaboration operated a scaled-down replica of the solid scintillator tracking and sampling calorimeter regions of the MINERvA detector in a hadron test beam at the Fermilab Test Beam Facility. This article reports measurements with samples of protons, pions, and electrons from 0.35 to 2.0 GeV/c momentum. The calorimetric response to protons, pions, and electrons are obtained from these data. A measurement of the parameter in Birks' law and an estimate of the tracking efficiency are extracted from the proton sample. Overall the data are well described by a Geant4-based Monte Carlo simulation of the detector and particle interactions with agreements better than 4%, though some features of the data are not precisely modeled. These measurements are used to tune the MINERvA detector simulation and evaluate systematic uncertainties in support of the MINERvA neutrino cross section measurement program. △ Less

Submitted 7 April, 2015; v1 submitted 26 January, 2015; originally announced January 2015.

Comments: as accepted by NIM A

Report number: FERMILAB-PUB-15-018-ND

Adiabatic-connection fluctuation-dissipation DFT for the structural properties of solids-the renormalized ALDA and electron gas kernels

Authors: Christopher E. Patrick, Kristian S. Thygesen

Abstract: We present calculations of the correlation energies of crystalline solids and isolated systems within the adiabatic-connection fluctuation-dissipation formulation of density-functional theory. We perform a quantitative comparison of a set of model exchange-correlation kernels originally derived for the homogeneous electron gas (HEG), including the recently-introduced renormalized adiabatic local-d… ▽ More We present calculations of the correlation energies of crystalline solids and isolated systems within the adiabatic-connection fluctuation-dissipation formulation of density-functional theory. We perform a quantitative comparison of a set of model exchange-correlation kernels originally derived for the homogeneous electron gas (HEG), including the recently-introduced renormalized adiabatic local-density approximation (rALDA) and also kernels which (a) satisfy known exact limits of the HEG, (b) carry a frequency dependence or (c) display a 1/$k^2$ divergence for small wavevectors. After generalizing the kernels to inhomogeneous systems through a reciprocal-space averaging procedure, we calculate the lattice constants and bulk moduli of a test set of 10 solids consisting of tetrahedrally-bonded semiconductors (C, Si, SiC), ionic compounds (MgO, LiCl, LiF) and metals (Al, Na, Cu, Pd). We also consider the atomization energy of the H$_2$ molecule. We compare the results calculated with different kernels to those obtained from the random-phase approximation (RPA) and to experimental measurements. We demonstrate that the model kernels correct the RPA's tendency to overestimate the magnitude of the correlation energy whilst maintaining a high-accuracy description of structural properties. △ Less

Submitted 5 January, 2015; originally announced January 2015.

Comments: 41 pages, 7 figures

Journal ref: J. Chem. Phys. 143, 102802 (2015)

Measurement of muon plus proton final states in $ν_μ$ Interactions on Hydrocarbon at $\langle$$E_ν$$\rangle$ = 4.2 GeV

Authors: T. Walton, M. Betancourt, L. Aliaga, O. Altinok, A. Bodek, A. Bravar, H. Budd, M. J. Bustamante, A. Butkevich, D. A. Martinez Caicedo, M. F. Carneiro, C. M. Castromonte, M. E. Christy, J. Chvojka, H. da Motta, M. Datta, J. Devan, S. A. Dytman, G. A. Díaz, B. Eberly, J. Felix, L. Fields, R. Fine, G. A. Fiorentini, A. M. Gago , et al. (51 additional authors not shown)

Abstract: A study of charged-current muon neutrino scattering on hydrocarbon in which the final state includes a muon and a proton and no pions is presented. Although this signature has the topology of neutrino quasielastic scattering from neutrons, the event sample contains contributions from both quasielastic and inelastic processes where pions are absorbed in the nucleus. The analysis accepts events with… ▽ More A study of charged-current muon neutrino scattering on hydrocarbon in which the final state includes a muon and a proton and no pions is presented. Although this signature has the topology of neutrino quasielastic scattering from neutrons, the event sample contains contributions from both quasielastic and inelastic processes where pions are absorbed in the nucleus. The analysis accepts events with muon production angles up to 70$^{\circ}$ and proton kinetic energies greater than 110 MeV. The extracted cross section, when based completely on hadronic kinematics, is well-described by a simple relativistic Fermi gas nuclear model including the neutrino event generator modeling for inelastic processes and particle transportation through the nucleus. This is in contrast to the quasielastic cross section based on muon kinematics, which is best described by an extended model that incorporates multi-nucleon correlations. This measurement guides the formulation of a complete description of neutrino-nucleus interactions that encompasses the hadronic as well as the leptonic aspects of this process. △ Less

Submitted 6 April, 2015; v1 submitted 16 September, 2014; originally announced September 2014.

Comments: total 10 pages, 5 figures, 2 pages of supplemental materials. Published in Physical Review D (Rapid Communication) April 1, 2015

Report number: FERMILAB-PUB-14-330-E

Journal ref: Phys. Rev. D 91, 071301 (2015)

Measurement of Coherent Production of $π^\pm$ in Neutrino and Anti-Neutrino Beams on Carbon from $E_ν$ of $1.5$ to $20$ GeV

Authors: A. Higuera, A. Mislivec, L. Aliaga, O. Altinok, A. Bercellie, M. Betancourt, A. Bodek, A. Bravar, W. K. Brooks, H. Budd, A. Butkevich, M. F. Carneiro, C. M. Castromonte, M. E. Christy, J. Chvojka, H. da Motta, J. Devan, S. A. Dytman, G. A. Díaz, B. Eberly, J. Felix, L. Fields, R. Fine, G. A. Fiorentini, H. Gallagher , et al. (47 additional authors not shown)

Abstract: Neutrino-induced coherent charged pion production on nuclei, $\stackrel{(-)}ν_μA\toμ^\pmπ^\mp A$ is a rare, inelastic interaction in which a small squared four-momentum $| t|$ is transferred to the recoil nucleus leaving it intact in the reaction. In the scintillator tracker of MINERvA, we remove events with evidence of particles from nuclear breakup and reconstruct $| t|$ from the final state pio… ▽ More Neutrino-induced coherent charged pion production on nuclei, $\stackrel{(-)}ν_μA\toμ^\pmπ^\mp A$ is a rare, inelastic interaction in which a small squared four-momentum $| t|$ is transferred to the recoil nucleus leaving it intact in the reaction. In the scintillator tracker of MINERvA, we remove events with evidence of particles from nuclear breakup and reconstruct $| t|$ from the final state pion and muon. We select low $| t|$ events to isolate a sample rich in coherent candidates. By selecting low $| t|$ events we produce a model-independent measurement of the differential cross section for coherent scattering of neutrinos and anti-neutrinos on carbon. We find poor agreement with the predicted kinematics in neutrino generators used by current oscillation experiments. △ Less

Submitted 25 November, 2014; v1 submitted 12 September, 2014; originally announced September 2014.

Comments: 13 pages, 4 figures. v3->v4 corrections to several inline formulas and fixes to citations. v2->v3 revision includes minor corrections to text and tables, the final ArgoNeuT result in Figure 3, and ancillary files with the tables in CSV format. v1->v2 revision added comparisons to NEUT generator and updated estimate of the diffractive-like process on protons

Report number: FERMILAB-PUB-14-331-E

Journal ref: Phys. Rev. Lett. 113, 261802 (2014)

Charged Pion Production in $ν_μ$ Interactions on Hydrocarbon at $\langle E_ν\rangle$= 4.0 GeV

Authors: B. Eberly, L. Aliaga, O. Altinok, M. G. Barrios Sazo, L. Bellantoni, M. Betancourt, A. Bodek, A. Bravar, H. Budd, M. J. Bustamante, A. Butkevich, D. A. Martinez Caicedo, M. F. Carneiro, M. E. Christy, J. Chvojka, H. da Motta, M. Datta, J. Devan, S. A. Dytman, G. A. Díaz, J. Felix, L. Fields, R. Fine, G. A. Fiorentini, A. M. Gago , et al. (56 additional authors not shown)

Abstract: Charged pion production via charged current $ν_μ$ interactions on plastic (CH) is studied using the MINERvA detector exposed to the NuMI wideband neutrino beam at Fermilab. Events with hadronic invariant mass W $<$ 1.4 GeV are selected to isolate single pion production, which is expected to occur primarily through the $Δ(1232)$ resonance. Cross sections as functions of pion production angle and ki… ▽ More Charged pion production via charged current $ν_μ$ interactions on plastic (CH) is studied using the MINERvA detector exposed to the NuMI wideband neutrino beam at Fermilab. Events with hadronic invariant mass W $<$ 1.4 GeV are selected to isolate single pion production, which is expected to occur primarily through the $Δ(1232)$ resonance. Cross sections as functions of pion production angle and kinetic energy are reported and compared to predictions from different theoretical calculations and generator-based models, for neutrinos ranging in energy from 1.5 GeV to 10 GeV. The data are best described by calculations which include significant contributions from pion intranuclear rescattering. These measurements constrain the primary interaction rate and the role of final state interactions in pion production, both of which need to be well understood by neutrino oscillation experiments. △ Less

Submitted 31 July, 2015; v1 submitted 24 June, 2014; originally announced June 2014.

Comments: 23 pages, 13 figures. Submitted to Physical Review D

Report number: FERMILAB-PUB-14-193-E

Journal ref: Phys. Rev. D 92, 092008 (2015)

Measurement of Ratios of $ν_μ$ Charged-Current Cross Sections on C, Fe, and Pb to CH at Neutrino Energies 2-20 GeV

Authors: B. G. Tice, M. Datta, J. Mousseau, L. Aliaga, O. Altinok, M. G. Barrios Sazo, M. Betancourt, A. Bodek, A. Bravar, W. K. Brooks, H. Budd, M. J. Bustamante, A. Butkevich, D. A. Martinez Caicedo, C. M. Castromonte, M. E. Christy, J. Chvojka, H. da Motta, J. Devan, S. A. Dytman, G. A. Díaz, B. Eberly, J. Felix, L. Fields, G. A. Fiorentini , et al. (52 additional authors not shown)

Abstract: We present measurements of $ν_μ$ charged-current cross section ratios on carbon, iron, and lead relative to a scintillator (CH) using the fine-grained MINERvA detector exposed to the NuMI neutrino beam at Fermilab. The measurements utilize events of energies $2<E_ν<20~GeV$, with $\left< E_ν\right>=8~GeV$, which have a reconstructed $μ^{-}$ scattering angle less than $17^\circ$ to extract ratios of… ▽ More We present measurements of $ν_μ$ charged-current cross section ratios on carbon, iron, and lead relative to a scintillator (CH) using the fine-grained MINERvA detector exposed to the NuMI neutrino beam at Fermilab. The measurements utilize events of energies $2<E_ν<20~GeV$, with $\left< E_ν\right>=8~GeV$, which have a reconstructed $μ^{-}$ scattering angle less than $17^\circ$ to extract ratios of inclusive total cross sections as a function of neutrino energy $E_ν$ and flux-integrated differential cross sections with respect to the Bjorken scaling variable $x$. These results provide the first high-statistics direct measurements of nuclear effects in neutrino scattering using different targets in the same neutrino beam. Measured cross section ratios exhibit a relative depletion at low $x$ and enhancement at large $x$. Both become more pronounced as the nucleon number of the target nucleus increases. The data are not reproduced by GENIE, a conventional neutrino-nucleus scattering simulation, or by the alternative models for the nuclear dependence of inelastic scattering that are considered. △ Less

Submitted 10 March, 2015; v1 submitted 9 March, 2014; originally announced March 2014.

Comments: 6 pages, 2 figures, 7 pages of supplemental material. Published in Physical Review Letters June 13, 2014

Report number: FERMILAB-PUB-14-025-E

Journal ref: Phys. Rev. Lett. 112, 231801 (2014)

Measurement of Muon Neutrino Quasi-Elastic Scattering on a Hydrocarbon Target at E_ν ~ 3.5 GeV

Authors: The MINERvA collaboration, G. A. Fiorentini, D. W. Schmitz, P. A. Rodrigues, L. Aliaga, O. Altinok, B. Baldin, A. Baumbaugh, A. Bodek, D. Boehnlein, S. Boyd, R. Bradford, W. K. Brooks, H. Budd, A. Butkevich, D. A. Martinez Caicedo, C. M. Castromonte, M. E. Christy, H. Chung, J. Chvojka, M. Clark, H. da Motta, D. S. Damiani, I. Danko, M. Datta , et al. (93 additional authors not shown)

Abstract: We report a study of muon neutrino charged-current quasi-elastic events in the segmented scintillator inner tracker of the MINERvA experiment running in the NuMI neutrino beam at Fermilab. The events were selected by requiring a μ^- and low calorimetric recoil energy separated from the interaction vertex. We measure the flux-averaged differential cross-section, dσ/dQ^2, and study the low energy pa… ▽ More We report a study of muon neutrino charged-current quasi-elastic events in the segmented scintillator inner tracker of the MINERvA experiment running in the NuMI neutrino beam at Fermilab. The events were selected by requiring a μ^- and low calorimetric recoil energy separated from the interaction vertex. We measure the flux-averaged differential cross-section, dσ/dQ^2, and study the low energy particle content of the final state. Deviations are found between the measured dσ/dQ^2 and the expectations of a model of independent nucleons in a relativistic Fermi gas. We also observe an excess of energy near the vertex consistent with multiple protons in the final state. △ Less

Submitted 30 March, 2014; v1 submitted 9 May, 2013; originally announced May 2013.

Comments: 7 pages, 5 figures. Added correlation between neutrino and anti-neutrino result in ancillary files (CSV)

Report number: FERMILAB-PUB-13-146-E

Journal ref: Phys. Rev. Lett. 111, 022502 (2013)