Strain already writes a three-hundred-tesla field into graphene, and no magnet did it

Scan a strained graphene nanobubble with a tunnelling microscope and the electrons inside it behave as though they were sitting in a magnetic field of more than 300 tesla — a field with no coil, no current, and no magnet anywhere near the sample. Nalini Levy and collaborators reported the effect in 2010, reading pseudo-magnetic fields in excess of 300 tesla out of bubbles only a few nanometres across, with theory placing the peak curvature regions at 200 to 400 tesla, against a strongest sustained laboratory magnet at the time of roughly 85 tesla [1]. The comparison is not close. The strained sheet of carbon beats the best magnet physics can build by a factor of three to five, and it does so because of a fact about graphene’s electrons rather than any exotic material: strain, in a material whose low-energy electrons obey a Dirac equation, enters that equation exactly the way a gauge field would. Bend the lattice and the electrons cannot tell the difference between the lattice distortion and an honest magnetic field.

A new paper, “A Strained Crystal Is a Designer Spacetime,” the fourth entry in The Crossed Fields Papers, takes that fact and asks how far it goes. The 300-tesla record is not a laboratory curiosity that happened once: architected nanopillar substrates push strain-free monolayer graphene to pseudo-fields near 800 tesla at room temperature over a scalable area [14], and heteroepitaxial silicon-carbide nanoprisms produce a uniform 41-tesla field, wafer-scale, with its pseudo-Landau levels resolved by photoemission at room temperature rather than deep cryogenic cooling [13]. When the object being distorted is the electron’s own dispersion relation, geometry is not a metaphor for the physics — it is the largest effective field condensed-matter physics has ever produced, and it is designed rather than discovered. The paper’s project is to take that same strain-to-geometry mapping into three dimensions, where a different class of material and a different feature of the electron’s motion — not its magnetic response but the tilt of its energy cone — turns an engineered strain profile into something closer to a curved spacetime than a magnetic field.

Read the full paper (PDF)

ADVERTISEMENT

A tilted Dirac cone is a curved spacetime for the electrons that live in it

Graphene’s electrons obey a two-dimensional Dirac equation with an upright cone: plot energy against momentum near the point where the bands touch, and the result is a symmetric cone, tip down, identical in every direction. A three-dimensional Weyl semimetal has the same kind of cone, but a strain field acting on it does two things at once rather than one. It shifts the cone’s tip sideways in momentum space — a chiral gauge field, the direct three-dimensional cousin of graphene’s pseudo-magnetic field — and it also tilts the cone, leaning it over in a particular direction the way a Dirac cone in flat graphene never leans. A cone that leans far enough stops looking like a cone at all: instead of a closed surface of constant energy around the tip, the surface opens into a hyperboloid, and the material has crossed from what the field calls a type-I Weyl phase into a type-II one, a transition Alexey Soluyanov and collaborators predicted theoretically for the layered semimetal tungsten ditelluride in 2015 [2].

The tilt is where the geometry stops being a figure of speech. Once the tilt varies smoothly from place to place across a crystal, an electron propagating through it is, at leading order in the standard low-energy expansion, equivalent to a particle moving through a curved effective spacetime of the same mathematical form — a Painlevé–Gullstrand metric — that describes an object falling into a black hole, with the tilt itself playing the role of the metric’s off-diagonal, “flow” components. The paper writes the whole statement as a single Hamiltonian, fixing the Fermi velocity vFv_F↗, the Pauli matrices σ\boldsymbol{\sigma}↗ acting on the electron’s internal degree of freedom, the momentum-space shift b\mathbf{b}↗, and the tilt vector W(x)W(x)↗:

H=vF σ⋅(k−b)+vF W(x) k⋅1 H = v_F\,\boldsymbol{\sigma}\cdot(\mathbf{k}-\mathbf{b}) + v_F\,W(x)\,\mathbf{k}\cdot\mathbb{1} ↗

The first term is ordinary tilted-free Weyl physics, shifted by the strain-induced gauge field. The second term is the one that matters here: it multiplies the momentum by the tilt and by the identity matrix, so it does not distinguish between the electron’s two internal states the way the first term does — it simply drags the whole cone in one direction, by an amount that depends on position through W(x)W(x)↗. Where ∣W∣|W|↗ crosses the value 1, the type-I/type-II boundary, the cone tips over far enough that one branch of carriers can only move one way past that line, in exactly the sense that nothing can move outward past an event horizon. That correspondence between an over-tilted Weyl cone and a horizon is not this paper’s invention — it descends from Grigory Volovik’s 2016 observation that a type-II Weyl semimetal could in principle host a Hawking-radiating horizon at room temperature if the tilted region were made small enough, stated as a symbolic relation between the Hawking temperature and the velocity gradient at the horizon, with no material, no computed number, and no fabrication route attached [15]. What the new paper adds is everything Volovik’s ancestor claim left as a symbol: a specific material, a specific device geometry, and a computed temperature.

Bend a near-critical crystal until the tilt crosses one, and the crossing is a horizon

Turning that mathematics into a device means choosing a crystal whose tilt sits close to the critical value already, at rest, so that an achievable amount of strain can push it the rest of the way across. Push too little strain onto a crystal whose tilt starts far below 1 and the crossing never arrives before the crystal’s elastic limit does; start from a crystal whose tilt is already far above 1 and the crossing needs a strain the lattice cannot survive. The paper’s candidate is the tungsten ditelluride and molybdenum ditelluride family, sitting close enough to the type-I/type-II boundary that density-functional calculations already show strain can drive the transition, with the tungsten-doped compound Mo0.91_{0.91}↗W0.09_{0.09}↗Te2_2↗ as the concrete lead, because it is the one member of that family on which strain has actually been delivered and measured in the laboratory [3].

The device geometry is a bending-beam flexure cell: clamp a thin slab of the crystal at its ends, bend it to a controlled radius of curvature, and the resulting surface strain varies smoothly along the beam, sweeping the local tilt through the critical value across a channel roughly a hundred nanometres to a micrometre long. The line where the tilt crosses 1 is the horizon, with an analog surface gravity set by how fast the tilt changes along the channel and an analog Hawking temperature that follows from it in the same way a real event horizon’s surface gravity fixes its radiation temperature. At the paper’s chosen design point — a Fermi velocity of 2×1052\times10^5↗ metres per second, a strain of three-tenths of one percent delivered over a hundred-nanometre channel — the analog Hawking temperature comes out around a tenth of a kelvin, comfortably inside dilution-refrigerator range, with the design envelope spanning roughly 0.04 to 0.4 kelvin depending on how the channel length and strain sensitivity are chosen. None of the mechanical pieces of that design are new on their own: the strain itself can be delivered by a piezoelectric bending-beam cell of the kind Clifford Hicks and Andrew Mackenzie’s group built for tuning correlated-electron materials, reaching continuously tunable, reversible strain up to about 0.23 percent at cryogenic temperature [7], and ohmic contacts to a transition-metal ditelluride channel are available by laser phase-patterning a metallic region directly into the semiconducting crystal, a technique that raised a MoTe2_2↗ transistor’s mobility roughly fifty-fold while holding a million-to-one on/off ratio [8].

ADVERTISEMENT
A wedge wire bonder's capillary needle mid-arc over a chip mounted in a sample puck, one bond pad still unconnected beside a completed gold wire loop
Figure 1. Getting a strained-crystal device from a bent flake to a working circuit still runs through ordinary microfabrication rather than anything exotic — the same contact technology that raised a MoTe2 transistor's mobility about fifty-fold once its electrodes were laser phase-patterned into a metallic, ohmic region rather than left as a bare semiconductor junction [@cho-et-al-2015].

The number nobody has measured yet: how hard strain actually tilts a real cone

Here the paper states its first honest limitation before drawing another line of the device. Every strain-to-tilt sensitivity that goes into the temperature above is computed, not measured — no experiment between 2016 and 2026 reports a tilt parameter under strain in any real type-II material. The most recent attempt is the same Mo0.91_{0.91}↗W0.09_{0.09}↗Te2_2↗ operando-strain study the device borrows its material choice from: Lanz and collaborators delivered tensile strain up to 0.34 percent along the crystal’s a-axis, reversible below about 0.17 percent, and measured a roughly 30-percent increase in the density of states at the Fermi level using strain-resolved photoemission — a careful, state-of-the-art measurement of how the bands respond to strain, and one that stops short, by the authors’ own account, of reporting a tilt parameter or locating a type-I/type-II crossing [3]. The gap is not a matter of a less careful group choosing not to look; it is that nobody has yet built the measurement that isolates the tilt itself from the other things strain does to a band structure.

Because that sensitivity is unmeasured, the paper computes it instead, anchoring the calculation to the density-functional strain sensitivities already published for this material family and to a Grüneisen-type estimate relating the tilt to the same electron-lattice coupling that fixes graphene’s own strain response, and reports the result as a band spanning roughly a factor of six rather than a single number. That spread propagates directly into the temperature: sweeping the computed sensitivity across its conservative-to-optimistic range moves the analog Hawking temperature by roughly the same factor of six the sensitivity itself carries. The paper does not treat that spread as a flaw to explain away. It treats measuring ∂W/∂strain\partial W/\partial\text{strain}↗ directly, in an operando-strain photoemission sweep of exactly the kind Lanz and collaborators already built most of the apparatus for, as the device’s first and most important milestone — the single experiment that would collapse today’s sixfold uncertainty into a real number.

A dilution-refrigerator sample puck carrying a flexure bending-beam strain cell, suspended a hair above its mixing-chamber socket with one cold-finger clip not yet closed
Figure 2. This is the part of the design that has to work at cryogenic temperature and still bend a crystal on demand: a piezoelectric bending-beam cell in the Hicks–Mackenzie lineage delivers continuously tunable, reversible strain up to about 0.23 percent at cryogenic temperature, which is the strain-delivery method the device's tilt-sweep arithmetic assumes [@hicks-et-al-2014].

A horizon is a race against disorder, and today no material wins it

A tilted cone crossing 1 on a band diagram is not automatically a horizon an electron experiences. For the crossing to be a physical event rather than a line on a chart, a carrier has to travel across it ballistically — without being scattered by disorder — for at least as long as the horizon takes to act on it. Turning that single requirement into arithmetic sets a minimum mobility: too dirty a crystal, and disorder buries the horizon under ordinary resistive scattering before any carrier gets close enough to feel the tipped cone. At the paper’s design point, a tenth-of-a-kelvin horizon demands a mobility on the order of 7×105 cm2/V⋅s7\times10^5\,\text{cm}^2/\text{V·s}↗, with carriers needing to cross roughly fourteen micrometres without scattering to register it.

Measured against real crystals, that bar splits the shortlist cleanly in two, and neither half clears it whole. The near-critical type-II hosts that actually carry the tilt — tungsten ditelluride, molybdenum ditelluride, and the strain-tunable tungsten-doped compound — sit at measured mobilities in the low thousands to low tens of thousands of square centimetres per volt-second, with tungsten ditelluride’s magnetotransport mobility around 5×103 cm2/V⋅s5\times10^3\,\text{cm}^2/\text{V·s}↗ [12], roughly two orders of magnitude below what the design point needs. The crystals that clear the mobility bar with room to spare sit on the other side of the same divide: the Dirac semimetal cadmium arsenide reaches mobilities above 107 cm2/V⋅s10^7\,\text{cm}^2/\text{V·s}↗ below four kelvin [5], and the Weyl semimetal niobium phosphide reaches 5×106 cm2/V⋅s5\times10^6\,\text{cm}^2/\text{V·s}↗ [6] — both comfortably ultraclean, and both type-I, with cones that sit upright and carry no tilt to bend into a horizon at all. Even suspended graphene, the two-dimensional mobility champion at roughly 2×105 cm2/V⋅s2\times10^5\,\text{cm}^2/\text{V·s}↗ [10], fails on both counts at once: short of the threshold, and a two-dimensional Dirac system with no cone tilt in the first place. The property that lets a crystal carry the geometry and the property that lets a carrier survive long enough to see it currently live in different crystals, and no material shortlist entry has both.

That is not, on the paper’s own reading, a wall that closes the project. Because the required mobility falls as the horizon warms, a design aimed at a few kelvin rather than a tenth of a kelvin needs roughly two orders of magnitude less mobility than the sub-kelvin design does — pulling the target down toward what today’s type-II crystal growth can plausibly deliver rather than what it delivers today. The honest first target, on this accounting, is not a sub-kelvin horizon in today’s tungsten ditelluride, but a few-kelvin horizon in a demonstrably cleaner type-II crystal — a growth problem stated as a number, rather than a physics problem stated as an obstacle.

A sputter and thermal-evaporator deposition chamber with its viewport door swung open and an internal source shutter caught half open above a mounted substrate
Figure 3. The suspended-membrane route this device follows — exfoliation over a trench, clamped at the contacts, released to bend — depends on exactly this kind of deposition to lay down clean contact metal without contaminating the released crystal underneath.

What “analog” buys here, and what it forbids

None of this is a claim about gravity, and the paper is explicit about exactly where the line sits. A companion paper in the same series, “What a Bathtub Can and Cannot Know About a Black Hole,” draws the general version of that line: an observable in an analog-gravity system is kinematic when it is fixed by the wave equation on an effective metric, the field’s quantum state, and the dispersion relation alone, and dynamic when its black-hole counterpart depends on the Einstein field equations, their matter sources, or the area-entropy law that ties a horizon’s size to its entropy [16]. This device sits entirely on the kinematic side of that line, by construction. The tilted-cone Hamiltonian above is a statement about the electrons’ dispersion and nothing else, so everything it can produce — an emission spectrum, a transmission probability, a dissipation pattern — is a genuine analog of the corresponding black-hole quantity in exactly the narrow sense the partition allows. What it cannot produce is anything that depends on spacetime actually responding to the electrons: there is no Bekenstein-Hawking entropy here, because that quantity comes from the area law and the field equations, not from a dispersion relation; there is no evaporation lifetime or Page curve, because those are histories of a dynamical spacetime; and the device’s own “back-reaction” is the crystal’s ordinary elasticity and electronic screening, not spacetime evolving under Einstein’s equations. The paper engineers a designer band-geometry for its electrons. It does not, at any point, claim to have built a piece of gravity.

ADVERTISEMENT

Three ways to see it, at three very different prices

Assume, for the sake of the argument, that the mobility bar above is eventually cleared in a cleaner type-II crystal at a warmer design point. Three measurements would then distinguish a genuine horizon from an ordinary strained resistor, each kinematic in the sense just fixed, and each computed by the paper at magnitudes rather than merely asserted to exist. The most direct, and the most demanding, is shot noise across the crossing: charge quantization makes current noisy in a way that carries a temperature, and the horizon should add an excess noise temperature set by the analog Hawking temperature itself rather than by the cryostat’s bath, visible once the bath is cooled to roughly that same temperature or below — a dilution-refrigerator measurement, not merely a cold one.

The cheapest of the three needs no such cooling at all. Christophe De Beule and collaborators showed that a slowly varying tilt horizon transmits carriers through Hawking-like thermal channels, and that the effect shows up as a peak in the ordinary two-terminal differential conductance, of order 2e2/h2e^2/h↗ per conducting mode, whose position tracks the strain-set tilt slope directly [4]. Because the peak is a property of how carriers transmit through the tipped cone rather than a property of thermal equilibrium, it is present as soon as the horizon exists, with no requirement to cool below the analog Hawking temperature — the signature a first, warmer device should chase.

A tungsten probe tip hovering just above a contact pad on a bent flake device mounted in a probe station, the chuck's vacuum line still unclipped beside it
Figure 4. This is the measurement the whole design is built to reach: a two-terminal conductance peak of order $2e^2/h$ per mode across the strained channel, tracking the tilt slope rather than sitting still under disorder — the signature De Beule and collaborators showed does not require cooling the device below its own analog Hawking temperature to see [@de-beule-2021].

The third sits exactly at its own instrument’s limit. A scanning SQUID-on-tip thermometer, a probe with a sub-50-nanometre diameter capable of imaging where power is actually dissipated in a device, has already resolved individual atomic-defect dissipation in graphene and reports a minimum detectable continuous power of about 40 femtowatts [9]. The horizon’s computed dissipation at the paper’s design point comes out to roughly half that floor — a signal the instrument can only just miss rather than comfortably clear, which the paper reports as the honest number rather than rounding it toward either a clean pass or a clean failure. No one of the three signatures is decisive by itself: a conductance peak can come from ordinary disorder, an excess noise temperature can come from stray heating, and a dissipation hot spot can come from a leaky contact. Together, tracking one strain-set channel at once, they are the argument the paper is actually proposing to run.

The strain-to-geometry mapping this paper works from belongs to the graphene community that built it over two decades, and the 300-tesla record that opens this story belongs to Levy and collaborators [1]; the type-II Weyl phase belongs to Soluyanov and collaborators [2]; the idea that an over-tilted cone reproduces a black-hole horizon traces back through Grigory Volovik [15] to William Unruh’s founding 1981 proposal that a sonic horizon in a moving fluid should radiate the same way a black hole does [11]; and the conductance-peak signature belongs to De Beule and collaborators [4]. What the new paper adds to that inheritance is the engineering: a quantitative strain-to-geometry table across named materials, a concrete device geometry with a temperature computed from real material parameters rather than left as a symbol, a disorder threshold stated as a mobility number rather than a caveat, and three signature magnitudes computed at one design point rather than asserted in the abstract. A verified audit of eighteen prior horizon proposals spanning 1998 to 2023 found no earlier paper combining even two of those four pieces. Strain already writes a designer field into graphene strong enough to embarrass every magnet on Earth; this paper’s contribution is showing, honestly and with the failure modes stated in the same breath as the design, what it would take to make the three-dimensional version of that field show an electron its own horizon.

Read the full paper (PDF)