Crystalline Current: Structure & Transport

• Crystalline current is a multifaceted transport phenomenon governed by crystal symmetry and order that induces charge, spin, ionic, and optical currents.
• The phenomenon elucidates how microscopic features like chirality, phase transitions, and topological constraints define key responses such as Edelstein magnetizations and vortex-antivortex lattice formation.
• Emerging theoretical models and experimental techniques are refining our understanding of crystalline current, offering new routes for high-efficiency energy, information, and quantum device applications.

In the literature surveyed here, “crystalline current” appears not as a single narrowly standardized term but as a family of transport phenomena in which the existence, direction, magnitude, or dynamical role of current is fixed by crystalline order, crystal symmetry, or the motion of a crystalline state. In this sense it includes longitudinal orbital and spin magnetizations induced by current in helical crystals, rigorously defined current per unit volume in perfect periodic lattices, topologically constrained current of a sliding electron crystal, crystallinity-dependent spin-current generation in metallic thin films, ionic conduction in plastic crystals, bulk photovoltaic and photon-drag currents, magnonic energy current carried by crystalline-symmetry-protected surface states, and current-induced formation of magnetic vortex crystals (Yoda et al., 2015, Cancès et al., 2020, Soejima et al., 24 Jul 2025, Kumar et al., 2018, Geirhos et al., 2015, Resta, 2024, Kondo et al., 2020, Gaididei et al., 2012).

1. Chirality, phase, and symmetry as controllers of current response

A particularly direct realization of crystalline current is the helical crystal studied in “Current-induced Orbital and Spin Magnetizations in Crystals with Helical Structure” (Yoda et al., 2015). There, a current driven along the helical axis induces both an orbital magnetization and, when spin–orbit coupling is included, a spin magnetization along that same axis. The sign reverses between right-handed and left-handed helices, making the response the condensed-matter analogue of reversing the winding of a solenoid. In the spinless model, the induced orbital magnetization is a metallic Fermi-surface effect controlled by the orbital magnetic moment and is strongly enhanced near Dirac points on the $K$-$H$ and $K'$-$H'$ lines. In the spinful extension, the relevant spin texture is radial rather than Rashba-like tangential, so the induced spin magnetization is parallel to the helical-axis current rather than perpendicular to it. The paper emphasizes that this longitudinal Edelstein-like response is symmetry-allowed specifically in chiral crystals without inversion and mirror symmetries, and highlights Te and Se as promising materials (Yoda et al., 2015).

Crystalline phase can also control how efficiently charge current is converted into spin current. In Py/Ta bilayers, the Ta crystalline phase is tuned by sputtering growth rate between predominantly $\alpha$-Ta, mixed $(\alpha+\beta)$-Ta, and predominantly $\beta$-Ta. The mixed-phase films exhibit the largest spin mixing conductance, the largest reported spin Hall angle in that data set, $\theta_{\rm SH}=-0.16\pm0.01$, and a spin Hall conductivity $$\sigma_{\rm SH}=-2439~(\hbar/e)~\Omega^{-1}\mathrm{cm}^{-1}$$ at $\rho=68~\mu\Omega\text{-cm}$, substantially exceeding the cited intrinsic benchmarks for pure $H$0- and $H$1-Ta (Kumar et al., 2018). This directly contradicts the simplified expectation that either elemental identity alone or the most resistive $H$2-Ta phase necessarily optimizes spin-current generation. A plausible implication is that “crystalline current” in spin-orbit torque materials is often microstructure-dependent rather than phase-purity-limited.

2. Rigorous and topological definitions of current in crystals

In perfect periodic crystals, current can be defined rigorously as an infinite-volume observable. “Coherent electronic transport in periodic crystals” formulates the current per unit volume as

$H$3

with $H$4 the trace per unit volume and $H$5 the time-evolved Fermi projector after a uniform electric field is switched on (Cancès et al., 2020). After a gauge transform, the field enters as a uniform drift in crystal momentum, $H$6. Within that single framework the paper derives, depending on how the Fermi level intersects the Bloch spectrum, the Berry-curvature/TKNN Hall response of insulators, ballistic growth $H$7 in clean metals, Bloch oscillations at longer times, and the finite universal conductivity of graphene-type semimetals (Cancès et al., 2020). The resulting taxonomy makes the meaning of “current in a crystal” explicitly spectral: gapped projector, Fermi surface, and isolated Dirac points produce qualitatively different asymptotics.

A more recent topological generalization concerns the current of a moving crystal itself. “Topological constraint on crystalline current” defines a sliding electron crystal by adiabatically dragging a pinned crystalline insulator and obtains the exact formula

$H$8

where $H$9 is electron density, $K'$0 magnetic flux density, $K'$1 the many-body Chern number, and $K'$2 the sliding velocity (Soejima et al., 24 Jul 2025). The topological correction $K'$3 means the current is not generically the naive convective current $K'$4. When $K'$5, a “full Hall crystal” slides with zero current. The same quantity fixes the Lorentz/Berry term

$K'$6

and with it the phonon counting rule: one gapless phonon if $K'$7, two if $K'$8 (Soejima et al., 24 Jul 2025). This sharpens a common misconception: motion of a charge-ordered state need not imply transport of the full electronic density.

3. Optical and photovoltaic crystalline currents

Optically generated crystalline currents occupy a parallel but conceptually related domain. “Geometrical theory of the shift current in presence of disorder and interaction” recasts shift current as a polarization-change problem: under illumination, population is transferred from the many-electron ground state to excited states with different polarization, so the dc current is the time derivative of the induced polarization (Resta, 2024). The central geometric quantity is the many-body shift vector

$K'$9

which yields the exact many-body formula

$H'$0

In the crystalline independent-electron limit this reduces to the standard Brillouin-zone formula with the band shift vector $H'$1, while in disordered systems it remains computable in a supercell formulation (Resta, 2024). The paper therefore places shift current within a general geometric framework that survives beyond ideal Bloch crystals.

“Jerk current: A novel bulk photovoltaic effect” introduces a distinct third-order dc photocurrent in biased crystalline insulators (Fregoso et al., 2018). In the clean, unsaturated limit it satisfies

$H'$2

so the current grows as $H'$3. The paper identifies two contributions: acceleration of optically injected carriers by the static field and a bias-induced modification of the injection rate itself. The latter is not captured by standard hydrodynamic descriptions, and the total jerk current can have a component perpendicular to the static field (Fregoso et al., 2018). A plausible implication is that nonlinear crystalline current cannot, in general, be reduced to “injection plus drift.”

The relation between symmetry and optical current becomes especially subtle in inversion-symmetric crystals. “Geometric photon-drag effect and nonlinear shift current in centrosymmetric crystals” argues that finite photon momentum activates a shift current that is otherwise absent at $H'$4. Expanding in photon momentum gives

$H'$5

with $H'$6 the “shift current dipole,” a geometric object built from interband transitions (Shi et al., 2020). Under suitable mirror or composite $H'$7 symmetry, the response in centrosymmetric crystals is purely transverse. By contrast, experiments on Bi$H'$8SiO$H'$9 show that Hall-based separation of shift and ballistic currents can be complicated by a magnetically induced bulk photocurrent $\alpha$0, which may be a significant fraction of the Hall-direction current under specific conditions (Burger et al., 2020). The controversy is not over whether shift current exists in such measurements, but over whether its quantitative separation from other bulk channels is unique when magnetic-field-induced photocurrent is present.

4. Nonperturbative and Floquet regimes in artificial crystals

Large-period artificial crystals make another regime accessible: current that is nonlinear to all orders in a static field but still describable within a single isolated miniband. In “Roses in the Nonperturbative Current Response of Artificial Crystals,” the relevant control parameter is the Bloch frequency

$\alpha$1

and the nonperturbative regime is $\alpha$2 (Beule et al., 2023). Within the relaxation-time approximation the steady-state distribution is resummed exactly in lattice-vector space, leading to all-orders expressions for the Bloch and geometric contributions to current. In trigonal systems these responses become symmetry-constrained rose curves as functions of field angle. For the first coordination shell, the transverse Bloch current carries a $\alpha$3 structure, yielding twelve petals in $\alpha$4, while the geometric current carries a $\alpha$5 structure, yielding six petals (Beule et al., 2023). Periodically buckled graphene is first-shell-dominated, whereas twisted double bilayer graphene receives significant multi-shell contributions. The result is a literal crystallization of the angular current pattern by the real-space harmonic content of miniband dispersion and Berry curvature.

Periodic driving can also imprint crystalline structure in time rather than angle. “Emergence and Dynamical Stability of Charge Time-Crystal in a Current-Carrying Quantum Dot Simulator” studies a driven open spinless Hubbard ladder implementable as a quantum-dot array and shows that a discrete time-crystal response can be read out directly in the drain current

$\alpha$6

The underlying mechanism is a hidden SU(2) symmetry of the ladder, a Floquet dynamical symmetry, and a special rung-wise dephasing channel (Sarkar et al., 2022). The current spectrum develops a sharp peak at a subharmonic frequency equivalent modulo $\alpha$7 to the DTC frequency, and the system can be tuned out of and back into the time-crystal phase by changing $\alpha$8 (Sarkar et al., 2022). This extends the notion of crystalline current from spatial symmetry to discrete time-translation symmetry breaking in transport itself.

5. Ionic, energy, and magnonic transport in crystalline media

Crystalline current is not restricted to electrons. In “Conductivity enhancement in plastic-crystalline solid-state electrolytes,” the material $\alpha$9 remains translationally crystalline while retaining orientational molecular dynamics (Geirhos et al., 2015). Replacing some succinonitrile by larger glutaronitrile enhances Li$(\alpha+\beta)$0 conductivity by more than one decade already at $(\alpha+\beta)$1 and by up to three decades at higher $(\alpha+\beta)$2. The conductivity crosses over from approximately Arrhenius behavior in pure SN to VFT-like behavior in GN-containing mixtures, and for $(\alpha+\beta)$3 the reorientational relaxation time $(\alpha+\beta)$4 and resistivity $(\alpha+\beta)$5 can be scaled almost perfectly onto one another (Geirhos et al., 2015). The interpretation is an optimized “revolving door” mechanism: molecular reorientation transiently opens pathways for ion hopping. This is a direct example of substantial current in a translationally ordered crystal enabled by orientational disorder.

An analogous rigor exists for lattice energy flow. “On the energy current for harmonic crystals” studies infinite harmonic crystals with random initial data asymptotically close to different translation-invariant states in different spatial sectors (Dudnikova, 2017). The paper proves convergence to a limiting Gaussian measure and derives explicit formulas for the mean energy current density. For Gibbs asymptotic states at different temperatures and under symmetry conditions $(\alpha+\beta)$6, the full-space stationary current for $(\alpha+\beta)$7 becomes

$(\alpha+\beta)$8

with $(\alpha+\beta)$9 (Dudnikova, 2017). The current is constant in space in full space, ballistic in character, and carried by phonon modes weighted by $\beta$0. In the half-space problem the current vanishes at the boundary and approaches the bulk limit far from the wall (Dudnikova, 2017).

A further bosonic realization arises in “Dirac surface states in magnonic analogs of topological crystalline insulators” (Kondo et al., 2020). There the protecting symmetry is

$\beta$1

the combination of time reversal and half translation, with $\beta$2 on the $\beta$3 plane. This yields a magnonic $\beta$4-type topology and a single Dirac cone on an $\beta$5-preserving surface. Under a homogeneous electric field, magnons couple through the Aharonov–Casher effect, shifting the two opposite surface Dirac cones in opposite directions in energy; with surface spin-momentum locking this produces an energy current quantified by

$\beta$6

Monoclinic CrI$\beta$7 is proposed as a candidate realization (Kondo et al., 2020). This suggests that crystalline current can be symmetry-protected surface energy flow in a neutral bosonic topological medium.

6. Current as an agent of crystalline order and high-current functionality

In some systems current does not merely respond to crystalline order but creates it. “Magnetic vortex-antivortex crystals generated by spin-polarized current” considers a soft ferromagnetic thin film driven by a dc spin-polarized current normal to the film plane (Gaididei et al., 2012). Above the critical current

$\beta$8

the film reaches a saturated out-of-plane state; slightly below $\beta$9, a finite-$\theta_{\rm SH}=-0.16\pm0.01$0 instability selects four symmetry-related modes that lock into a square vortex-antivortex superlattice (Gaididei et al., 2012). The weakly nonlinear solution has $\theta_{\rm SH}=-0.16\pm0.01$1 symmetry, alternating topological density, and a lattice spacing controlled by the selected wave number $\theta_{\rm SH}=-0.16\pm0.01$2. As the current is reduced further, the crystal melts into a fluid-like vortex-antivortex state. Here the phrase “crystalline current” is nearly literal: current is the organizing principle of the emergent lattice.

At the opposite end of the materials spectrum, current can expose unusually favorable crystalline functionality even in structurally clean samples. “Unusually high critical current of clean P-doped BaFe$\theta_{\rm SH}=-0.16\pm0.01$3As$\theta_{\rm SH}=-0.16\pm0.01$4 single crystalline thin film” reports a microstructurally clean, epitaxial, single-crystalline film with self-field $\theta_{\rm SH}=-0.16\pm0.01$5 at 4.2 K, $\theta_{\rm SH}=-0.16\pm0.01$6 at $\theta_{\rm SH}=-0.16\pm0.01$7 for $\theta_{\rm SH}=-0.16\pm0.01$8, and $\theta_{\rm SH}=-0.16\pm0.01$9 at $$\sigma_{\rm SH}=-2439~(\hbar/e)~\Omega^{-1}\mathrm{cm}^{-1}$$0 for $$\sigma_{\rm SH}=-2439~(\hbar/e)~\Omega^{-1}\mathrm{cm}^{-1}$$1, despite the absence of obvious structural pinning centers in XRD or TEM (Kurth et al., 2015). The interpretation given is enhanced vortex core energy near optimal $$\sigma_{\rm SH}=-2439~(\hbar/e)~\Omega^{-1}\mathrm{cm}^{-1}$$2, driven by in-plane tensile strain and P doping. This resists the common intuition that exceptionally high critical current in a crystal must imply a dense visible defect landscape.

Taken together, these results suggest that “crystalline current” is best understood as a structural principle rather than a single transport coefficient. In some contexts the crystal constrains current through symmetry, topology, or band geometry; in others current reveals, stabilizes, or even generates crystalline order. Across charge, spin, orbital, ionic, phononic, magnonic, and superconducting settings, the unifying feature is that transport is inseparable from the specific way order is organized in a crystal.