Refining hydrogen positions in {\alpha}-FeOOH through combined neutron diffraction and computational techniques

This study combines neutron diffraction and first-principles calculations to precisely determine the hydrogen positions and confirm the $b$-axis antiferromagnetic spin arrangement in $\alpha$-FeOOH, demonstrating the effectiveness of this dual approach for undeuterated compounds to aid in understanding CO$_2$ reduction mechanisms.

The Gist

The Big Picture: Finding the "Invisible Ghost" in Rust

Imagine you are trying to solve a 3D puzzle of a complex machine (in this case, a mineral called Goethite, which is basically the scientific name for the red rust you see on old bikes or nails). You know exactly where the heavy metal parts (Iron) and the oxygen parts are. But there is a tiny, invisible piece of the puzzle: Hydrogen.

Hydrogen is the "ghost" of the atomic world. It's so small and light that standard tools (like X-rays, which are used to take pictures of atoms) just bounce right off it. It's like trying to find a specific firefly in a stadium full of bright spotlights; the firefly is there, but you can't see it against the glare.

This paper is about a team of scientists who finally managed to pinpoint exactly where these "ghosts" (hydrogen atoms) are hiding inside the rust, and they did it by using two different superpowers at the same time.

---

The Two Superpowers

To find the invisible hydrogen, the scientists used a "tag-team" approach:

1. The Neutron Flashlight (Neutron Diffraction):

   • The Problem: X-rays see electrons. Hydrogen has very few electrons, so it's invisible to them.
   • The Solution: Neutrons are different. They don't care about electrons; they bounce off the nucleus (the center) of the atom. Hydrogen has a very strong nucleus for its size.
   • The Analogy: If X-rays are like a spotlight that misses a ghost, Neutrons are like a thermal camera that sees the ghost's body heat. The scientists shot a beam of neutrons at the rust powder. When the neutrons bounced off, they created a pattern that revealed where the hydrogen was.
   • The Catch: Hydrogen is tricky. It has a "negative" scattering length (a fancy way of saying it messes with the math in a weird way), which can sometimes make the picture blurry.

2. The Crystal Ball (First-Principles Calculations):

   • The Solution: Since the neutron picture wasn't 100% perfect on its own, they used a supercomputer to run a simulation. They built a virtual model of the rust and asked the computer: "If we arrange the atoms this way, what would the energy look like?"
   • The Analogy: This is like having a master architect who knows the laws of physics perfectly. Even if the photo is blurry, the architect can say, "Based on how gravity and structure work, this beam must be here."

The Discovery: A Perfect Match

When the scientists compared the "Neutron Flashlight" photo with the "Crystal Ball" simulation, they matched up perfectly.

• The Magnetic Spin: They also figured out how the tiny magnets inside the iron atoms are arranged. Imagine the iron atoms are like tiny compasses. The scientists found that these compasses are all pointing in opposite directions (North-South, South-North) along a specific line. This is called antiferromagnetism. It's like a crowd of people doing a "wave" where everyone stands up and sits down in perfect alternation.
• The Hydrogen Location: They found that the hydrogen atoms are attached to oxygen atoms, forming "OH" groups. Crucially, they mapped out exactly how these groups are angled and how close they are to each other.

Why Does This Matter? (The "Why Should I Care?" Factor)

You might be thinking, "So, we found where the hydrogen is in rust. Big deal."

Here is the deal: This rust is actually a superhero for the environment.

Scientists have discovered that this specific type of rust (Goethite) is a great catalyst for turning Carbon Dioxide (CO2)—the gas that causes climate change—into useful things like fuel or formic acid.

• The Analogy: Think of the rust surface as a factory floor. To turn CO2 into fuel, you need to deliver "protons" (which are just hydrogen nuclei) to the CO2 molecule.
• The Problem: If you don't know exactly where the hydrogen "delivery trucks" (the OH groups) are parked on the factory floor, you can't design an efficient delivery route.
• The Solution: By knowing the exact position and angle of the hydrogen atoms, scientists can now build better computer models to figure out how to make this rust work even better at cleaning up our atmosphere.

The "No-Deuterium" Trick

Usually, to make hydrogen easier to see in these experiments, scientists replace normal hydrogen with Deuterium (a heavier version of hydrogen, like swapping a ping-pong ball for a golf ball). This makes it easier to spot.

However, this team did something special: They didn't use Deuterium. They used normal hydrogen.

• Why is this cool? It proves that with the right combination of Neutron data and Supercomputer simulations, you don't need to change the chemistry of the sample to get a perfect picture. It's like solving a mystery using only a blurry photo and a brilliant detective, without needing to paint the suspect's face red to make them stand out.

The Takeaway

This paper is a victory for teamwork. By combining the physical "eyes" of neutron scattering with the "brain" of computer simulations, the scientists successfully mapped the invisible hydrogen atoms in rust. This map is now a crucial blueprint for developing new technologies that can help us fight climate change by turning pollution into fuel.

Technical Summary

1. Problem Statement

• Scientific Context: Goethite ($\alpha$-FeOOH) is a primary component of iron rust and a promising catalyst for the photochemical reduction of $CO_2$ to formic acid ($HCOOH$). The catalytic mechanism relies heavily on surface hydroxyl ($OH^-$) groups involved in proton-coupled electron transfer.
• The Challenge: Accurate knowledge of bulk hydrogen positions and hydrogen-bond topology is essential for modeling surface terminations and reaction mechanisms. However, determining hydrogen positions is difficult:
  • X-ray Diffraction: Fails due to hydrogen's low electron density.
  • Neutron Diffraction: While sensitive to hydrogen, standard Rietveld analysis on undeuterated (light hydrogen) samples can suffer from ambiguity because hydrogen has a negative coherent scattering length ($b_H = -3.739$ fm), which can cancel out contributions from other elements with positive scattering lengths, complicating intensity interpretation.
• Objective: To precisely determine the hydrogen positions and magnetic structure of $\alpha$-FeOOH using a powder sample without deuteration, by combining neutron diffraction with first-principles calculations.

2. Methodology

The study employed a dual approach combining experimental neutron diffraction with Density Functional Theory (DFT) calculations.

Experimental Approach:

• Sample: Commercial $\alpha$-FeOOH powder (99%+).
• Technique: Neutron powder diffraction using the HERMES diffractometer at JRR-3 (Japan Research Reactor).
• Conditions: Data collected at wavelengths $\lambda = 1.34239(7)$ Å over a temperature range of 60–485 K.
• Analysis:
  • Magnetic Structure: Analyzed using group theory (irreducible representations and magnetic space groups) under the space group $Pnma$ with magnetic wavevector $q_m = (0, 0, 0)$.
  • Refinement: Rietveld refinement performed using FullProf Suite with a pseudo-Voigt profile.
  • Temperature Dependence: Integrated intensities of magnetic peaks (200, 101) were used to determine the Néel temperature ($T_N$).

Computational Approach:

• Code: VASP (Projector Augmented-Wave method).
• Functional: PBEsol for exchange-correlation; DFT+U was applied to Fe 3d orbitals ($U = 5.3$ eV) to account for strong correlation.
• Configurations: Spin-polarized calculations for both ferromagnetic and various antiferromagnetic (AFM) configurations in conventional unit cells and supercells.
• Validation: Calculated lattice constants, atomic positions, and magnetization densities were compared against experimental data.

3. Key Contributions

• Group-Theoretical Magnetic Refinement: The authors systematically analyzed all 8 possible magnetic space groups (MSGs) allowed for the $Pnma$ structure with $q_m = (0, 0, 0)$. They identified the correct magnetic structure by comparing reliability factors ($R_{mag}$) for each MSG.
• Validation of Undeuterated Samples: The study demonstrates that accurate hydrogen positioning is achievable without deuteration or polarized neutrons, provided that neutron diffraction data is cross-validated with first-principles calculations.
• Comprehensive Structural Model: The paper provides a fully refined structural model including precise fractional coordinates for hydrogen, magnetic moments, and hydrogen-bond geometry across a wide temperature range.

4. Key Results

• Magnetic Structure:

  • Néel Temperature ($T_N$): Determined to be 382(2) K, consistent with previous literature.
  • Magnetic Ordering: The magnetic moments align along the b-axis in an antiferromagnetic arrangement.
  • Symmetry: The best fit corresponds to the magnetic space group $Pnma'$ (#62.445) with a reliability factor $R_{mag} = 8.53\%$. The magnetic unit cell matches the crystallographic unit cell.
  • Moment Size: The refined magnetic moment is 3.64(4) $\mu_B$ per $Fe^{3+}$ site.
  • Consistency: This result aligns with previous studies (when accounting for axis transformation between $Pbnm$ and $Pnma$ settings) and DFT predictions showing AFM coupling along [100] and FM coupling along [010].

• Hydrogen Positions and Geometry:

  • O–H Bond Length: Refined experimentally as 1.008(7) Å at 64.8 K. DFT calculated 1.02253 Å (within 1.5% error).
  • Hydrogen Bonding: The study visualized the evolution of the hydrogen bond network ($O2-H \cdots O1$).
    • Experimental (64.8 K): $H \cdots O1 = 1.781(7)$ Å, $O1 \cdots O2 = 2.734(4)$ Å, angle $\angle O2-H \cdots O1 = 156.4(6)^\circ$.
    • DFT: $H \cdots O1 = 1.586$ Å, $O1 \cdots O2 = 2.704$ Å, angle $164.3^\circ$.
  • Electronic Structure: DFT calculated a band gap of 2.16 eV. The valence band top and conduction band bottom are dominated by Fe(3d) and O(2p) orbitals. OH molecular orbitals are located in deeper valence bands compared to single oxide ions.

• Temperature Effects: No abrupt changes or anomalies were observed in the OH distance or hydrogen-bond geometry across the measured temperature range, indicating no proton order-disorder transition or symmetry change in the average structure.

5. Significance

• Methodological Advancement: The paper establishes a robust protocol for determining hydrogen positions in inorganic materials using undeuterated samples. This simplifies experimental procedures and broadens the applicability of neutron diffraction, as deuteration is often difficult or impossible for certain compounds.
• Catalytic Mechanism Insight: By providing a precise bulk reference for hydrogen positions and hydrogen-bond topology, the study offers a critical foundation for modeling surface reactions. This is vital for understanding the mechanism of $CO_2$ reduction and water oxidation on $\alpha$-FeOOH, where proton transfer steps are rate-limiting.
• Theoretical-Experimental Synergy: The high degree of agreement between the group-theoretical refinement of neutron data and DFT calculations validates the use of complementary techniques to resolve structural ambiguities that neither method could solve alone.
• Future Applications: The refined structural data serves as a consistent starting point for investigating surface defects, proton mobility, and in situ catalytic processes under $CO_2$ atmospheres.