Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Spatial dependence”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,027 records · Page 57

The dependence of aluminum lattice orientation on semiconductor lattice parameter in planar InAs/Al hybrid heterostructures

We present a detailed study of 7-nm thick Al layers epitaxially grown on In x Al 1-x As or Al 1-x Ga x Sb y As 1-y (0 0 1) semiconductors using high-resolution transmission electron microscopy (HRTEM) and synchrotron X-ray diffraction (XRD). These techniques provide high spatial and high reciprocal space resolution information about Al lattice orientations and interfacial structure, highlighting the existence of a strong correlation between the surface lattice parameter of the semiconductor and the normal orientation of the aluminum film. The normal orientation of Al layers is found to be [1 1 1] when the in-plane semiconductor surface lattice parameter, a S , is smaller than 5.98 Å, and is [1 1 0] when aS is larger than this value. In the transition between these lattice parameters where one orientation is preferred, both orientations may be observed. Concomitant with the change in normal orientation is an inversion in the strain state of the aluminum film. When the normal orientation is [1 1 1], the Al film is compressively strained; while for [1 1 0] normal orientations, the Al in-plane strain is tensile. An energy balance model, accounting for surface, interfacial, side-face, and strain energy, and the registry of lattice planes across the semiconductor/aluminum interface provide explanations for the observed normal orientation and strain inversion.

36 MATERIALS SCIENCE↗

Windowed least-squares model reduction for dynamical systems

Here we propose a windowed least-squares (WLS) approach for model reduction of dynamical systems. The proposed approach sequentially minimizes the time-continuous full-order-model residual within a low-dimensional space–time trial subspace over time windows. The approach comprises a generalization of existing model reduction approaches, as particular instances of the methodology recover Galerkin, least-squares Petrov–Galerkin (LSPG), and space–time LSPG projection. In addition, the approach addresses key deficiencies in existing model reduction techniques, e.g., the dependence of LSPG and space–time LSPG projection on the time discretization and the exponential growth in time exhibited by a posteriori error bounds for both Galerkin and LSPG projection. We consider two types of space–time trial subspaces within the proposed approach: one that reduces only the spatial dimension of the full-order model, and one that reduces both the spatial and temporal dimensions of the full-order model. For each type of trial subspace, we consider two different solution techniques: direct (i.e., discretize then optimize) and indirect (i.e., optimize then discretize). Numerical experiments conducted using trial subspaces characterized by spatial dimension reduction demonstrate that the WLS approach can yield more accurate solutions with lower space–time residuals than Galerkin and LSPG projection.

97 MATHEMATICS AND COMPUTING↗

An Atomistic Study of Reactivity in Solid-State Electrolyte Interphase Formation for Li/Li7P3S11

Lithium metal batteries offer superior volumetric and gravimetric specific capacities compared to those based on traditional graphite anodes. Although advancements in solid-state electrolytes address safety concerns, challenges remain, particularly regarding interphase formation in lithium metal anodes. This work presents a computational framework based on high-throughput first-principles density functional theory and machine-learning interatomic potentials (MLIPs) including automated iterative, active learning to enable robust computational exploration of interphase formation between lithium metal anodes and an inorganic solid-state electrolyte. As a demonstration, we apply the framework to a Li/Li7P3S11 interface and find that it accurately identifies the experimentally observed, thermodynamically stable interphase products as well as their overall spatial arrangement within a heterogeneous, amorphous layered structure, with Li2S domains of nanocrystallinity. Our simulations show two stages, a fast and slow diffusion reaction regime, that corroborate the relative phase formation rate of Li x P, Li2S, and Li3P. Using the Onsager transport theory, we capture time-dependent ionic diffusion within the reacting interface, including cross-correlation effects. We found that cross-correlation effects between Li-P and P-S ionic motion significantly influence P-ion diffusion, making it highly sensitive to the local environment and potentially leading to "kinetic trapping" of Li-P phases. The passivation of the interface is shown as the ionic fluxes all approach zero, effectively halting interphase growth.

Diffusion↗

ZeoNet: 3D convolutional neural networks for predicting adsorption in nanoporous zeolites

Zeolites are one of the most widely used materials in the chemical industry due to their nanometer-sized pores that can adsorb and react upon molecules selectively. With hundreds of known framework topologies and hundreds of thousands of computationally predicted structures, the ability to rapidly predict zeolite performance allows researchers to prioritize their efforts on the most promising structures for a given application. Although the accuracy of forcefield-based atomistic simulations has advanced significantly in the past two decades, these simulations can be computationally expensive, especially for long-chain, complex molecules. Here, we present ZeoNet, a representation learning framework using convolutional neural networks (ConvNets) and 3D volumetric representations for predicting adsorption in zeolites. ZeoNet was trained on the task of predicting Henry's constants for adsorption, k H , of n-octadecane in more than 330 000 known and predicted zeolite materials. Employing a 3D grid based on the distances to solvent-accessible surfaces, a volumetric representation that can be generated efficiently, the best-performing ZeoNet achieved a correlation coefficient r 2 = 0.977 and a mean-squared error MSE = 3.8 in ln k H , which corresponds to an error of 9.3 kJ mol -1 in adsorption free energy. In comparison, a model based on hand-designed geometric features has values of r 2 = 0.783 and MSE = 35.7. ZeoNet is also relatively efficient and can process ≈8 structures per second on an Nvidia RTX 2080TI GPU, orders of magnitude faster than forcefield-based simulations. A systematic analysis was conducted to investigate how the choice of ConvNet architectures, the linear dimension (L) and spatial resolution (Δd) of the distance grids, batch size, optimizer, and learning rate impact the model performance. We found that ConvNets based on the ResNet architecture offer the best tradeoff between expressiveness and efficiency. The performance for all models reaches a plateau at L = 30–45 Å and depends less sensitively on grid resolution, with a small benefit around Δd = 0.30–0.45 Å. Finally, saliency maps were visualized to identify which regions of the materials contributed the most to model predictions. It was found, interestingly, that the predictions are driven primarily by the accessible pore volume rather than the region occupied by the framework atoms.

36 MATERIALS SCIENCE↗

Topological Defect Mediated Helical Phase Reorientation by Uniaxial Stress

Strain engineering enables precise, energy-efficient control of nanoscale magnetism. However, unlike well-studied strain-dislocation interactions in mechanical deformation, the spatial evolution of strain-induced spin rearrangement remains poorly understood. Using Lorentz transmission electron microscopy, we manipulate and observe helical domain reorientation under quantitatively applied uniaxial tensile stress. Our findings reveal striking similarity to plastic deformation in metals, where the critical stress for propagation vector ( Q ) reorientation depends on its angle with the stress direction. Magnetic defects mediate reorientation via “break-and-reconnect” or “dislocation gliding–annihilation” processes. Simulations confirm that strain-induced anisotropic Dzyaloshinskii-Moriya interaction may play a key role. These insights advance strain-driven magnetism and offer a promising route for energy-efficient magnetic nanophase control in next-generation information technology. Published by the American Physical Society 2025

36 MATERIALS SCIENCE↗

Review of the November 11, 1905 Earthquake: Seismic Intensity Data, Magnitude, and Location Report

The macroseismic location and magnitude of the Shoshone 1905 earthquake has been based on the application of objective and quantitative approaches that produce estimates of earthquake source location and magnitude, and associated uncertainties. The performance of all three approaches applied depends critically on the characteristics of the seismic intensity field. Regional intensity attenuation relations were derived based on a calibration set of earthquakes and then, the procedure was validated using three independent very recent events. The macroseismic epicenter location was determined by either a grid search procedure, the spatial distribution of the sites experiencing the largest intensities, or a combination of both. The macroseismic magnitude estimates are based on the derived regional intensity attenuation relations. There is consistency in the location and magnitude of the November 1905 earthquake as obtained by the three methods. The overlapping region of the 68% bootstrap resampling solutions suggests a reliable epicenter location within 41.550N ? 42.060N and 113.478W ? 112.320W. This region encloses the original Oaks (1992) epicenter location. Macroseismic magnitude estimates range between 5.87 and 6.05.

58 GEOSCIENCES↗

Science Area 1: Standard Award: Model-Data Fusion to Examine Multiscale Dynamical Controls on Snow Cover and Critical Zone Moisture Inputs (Final Report)

In many mountain watersheds of the world, seasonal snowpacks play an important role as natural reservoirs of water. Seasonal snowpacks accumulate water during cold, wet winter months that subsequently melts. Downstream communities depend on water from melting seasonal snowpacks to support agricultural, industrial, and municipal water needs. Rapidly melting snowpacks can also present a flooding hazard, particularly if snowpacks melt at rates faster than anticipated and where adequate reservoir capacity is unavailable to buffer river flows associated with melt. The spatial and temporal dynamics of snow accumulation and melt also play an important role in supporting upland ecosystems in mountain landscapes. Snowmelt provides soil moisture that enable terrestrial ecosystem productivity and exert control on soil microorganisms that play important roles in global carbon cycles. Climate warming is gradually decreasing the amount of precipitation in mountain watersheds arriving as snow, presenting potentially profound disruptions to mountain ecosystems, as well as downstream delivery of water. The overarching goal of this project was to understand how interactions between the near-surface atmosphere and surface topography control the input, accumulation, retention, and release of water from mountain snowpacks. Over a 5-year period, this project pursued an approach combining high-resolution regional climate modeling, satellite and airborne remote sensing data, and ground-based observations to develop and analyze benchmark datasets to address overarching science questions and hypotheses. Key products include a continuous, long-term, high spatiotemporal resolution (1 km/1 hr) dataset characterizing key climate variables in the Upper Colorado River Basin. The dataset included historical estimates of precipitation, temperature, humidity, solar and longwave radiation, and wind speeds. Additionally, the project developed a 20+ year long, 30 m spatial, daily temporal multi-sensor dataset characterizing snow presence/absence in the East/Taylor River watersheds in the Central Rocky Mountains of Colorado. The project supported training of 1 postdoctoral scholar, 1 Ph.D. student, and 1 M.S. student.

54 ENVIRONMENTAL SCIENCES↗

Optimization of scintillated imaging detectors for MeV X-ray radiography

Radiographic imaging at X-ray energies in the MeV range poses a complex problem. Detection of high energy X-rays depends heavily on the ability of the detection material, in this work the scintillator, to attenuate and ultimately stop the X-ray. Attenuation is more readily achieved by increasing material thickness, however thick scintillators are plagued with relatively high internal scattering, leading to an overall decrease in spatial resolution and image quality.

42 ENGINEERING↗

The effects of storms and a transient sandy veneer on the interannual planform evolution of a low-relief coastal cliff and shore platform at Sargent Beach, Texas, USA

Coastal cliff erosion is alongshore-variable and episodic, with retreat rates that depend upon sediment as either tools of abrasion or protective cover. However, the feedbacks between coastal cliff planform morphology, retreat rate, and sediment cover are poorly quantified. This study investigates Sargent Beach, Texas, USA, at the annual to interannual scale to explore (1) the relationship between temporal and spatial variability in cliff retreat rate, roughness, and sinuosity and (2) the response of retreat rate and roughness to changes in sand and shell hash cover of the underlying mud substrate as well as the impact of major storms using field measurements of sediment cover, erosion, and aerial images to measure shore platform morphology and retreat. A storm event in 2009 increased the planform roughness and sinuosity of the coastal cliff at Sargent Beach. Following the storm, aerial-image-derived shorelines with annual resolution show a decrease in average alongshore erosion rates from 12 to 4 m yr –1 , coincident with a decrease in shoreline roughness and sinuosity (smoothing). Like the previous storm, a storm event in 2017 increased the planform roughness and sinuosity of the cliff. Over shorter timescales, monthly retreat of the sea cliff occurred only when the platform was sparsely covered with sediment cover on the shore platform, indicating that the tools and cover effects can significantly affect short-term erosion rates. The timescale to return to a smooth shoreline following a storm or roughening event, given a steady-state erosion rate, is approximately 24 years, with the long-term rate suggesting a maximum of ~107 years until Sargent Beach breaches, compromising the Gulf Intracoastal Waterway (GIWW) under current conditions and assuming no future storms or intervention. The observed retreat rate varies, both spatially and temporally, with cliff face morphology, demonstrating the importance of multi-scale measurements and analysis for interpretation of coastal processes and patterns of cliff retreat.

58 GEOSCIENCES↗

Cadmium alters whole animal ionome and promotes the re-distribution of iron in intestinal cells of Caenorhabditis elegans

The chronic exposure of humans to the toxic metal cadmium (Cd), either occupational or from food and air, causes various diseases, including neurodegenerative conditions, dysfunction of vital organs, and cancer. While the toxicology of Cd and its effect on the homeostasis of biologically relevant elements is increasingly recognized, the spatial distribution of Cd and other elements in Cd toxicity-caused diseases is still poorly understood. Here, we use Caenorhabditis elegans as a non-mammalian multicellular model system to determine the distribution of Cd at the tissue and cellular resolution and its effect on the internal levels and the distribution of biologically relevant elements. Using inductively coupled plasma-mass spectrophotometry (ICP-MS), we show that exposure of worms to Cd not only led to its internal accumulation but also significantly altered the C. elegans ionome. Specifically, Cd treatment was associated with increased levels of toxic elements such as arsenic (As) and rubidium (Rb) and a decreased accumulation of essential elements such as zinc (Zn), copper (Cu), manganese (Mn), calcium (Ca), cobalt (Co) and, depending on the Cd-concentration used in the assay, iron (Fe). We regarded these changes as an ionomic signature of Cd toxicity in C. elegans. We also show that supplementing nematode growth medium with Zn but not Cu, rescues Cd toxicity and that mutant worms lacking Zn transporters CDF-1 or SUR-7, or both are more sensitive to Cd toxicity. Finally, using synchrotron X-Ray fluorescence Microscopy (XRF), we showed that Cd significantly alters the spatial distribution of mineral elements. The effect of Cd on the distribution of Fe was particularly striking: while Fe was evenly distributed in intestinal cells of worms grown without Cd, in the presence of Cd, Fe, and Cd co-localized in punctum-like structures in the intestinal cells. Together, this study advances our understanding of the effect of Cd on the accumulation and distribution of biologically relevant elements. Considering that C. elegans possesses the principal tissues and cell types as humans, our data may have important implications for future therapeutic developments aiming to alleviate Cd-related pathologies in humans.

59 BASIC BIOLOGICAL SCIENCES↗

What's New in the NSRDB

The National Solar Radiation Database (NSRDB) provides solar resource data across the globe at a high temporal and spatial resolution. This data is primarily used in solar energy modeling. The NSRDB is updated annually for the United States and North, Central and South America and the data is currently available from 1998-2021. In 2022 the NSRDB was updated using the latest version of the underlying Physical Solar Model (PSM). This update includes improved surface albedo and gap-filling of cloud properties. The inclusion of these updates reduced the uncertainty in the data compared to previous versions of the NSRDB. The Himawari and Meteosat Indian Ocean Data Coverage (IODC) satellites were added to the Geostationary Operational Environmental Satellite (GOES) and made our coverage global. While standard data from the GOES continues to be served at an hourly 4km x 4km resolution, full resolution data has also been made available to the user. The NSRDB now contains over 200Tb of data with nearly 40Tb being added annually. We provide significant flexibility for data download depending on the amount of data required by the users. In this paper we provide an update on the current status on the NSRDB.

photovoltaic systems↗

The balance of orbital overlap and orbital energy in the activation of methane by actinide cations: insights from inductively coupled plasma tandem mass spectrometry

The actinides present a unique challenge to chemical theory. The classical view of covalent bonding is driven by the extent of spatial overlap of valence orbitals. Modern theory has expanded assessments of covalency to include considerations of orbital energy degeneracy to assess orbital energy mixing between metal and ligand valence orbitals. Actinide–ligand (An–L) bonding has more recently been described as a balance between orbital overlap and orbital energy mixing, where 5f and L valence orbital overlap decreases while energy mixing between An 5f and L valence orbitals increases across the series. To test these existing views, we employed inductively coupled plasma tandem mass spectrometry to examine the kinetic energy dependences of reactions of actinide cations, Th + –Am + , with methane. Further, this is the first experimental report of the energy dependences of methane activation reactions involving the cations of Pa, Np, Pu, and Am and the first experimental determination of transuranic An + –D, An + –CD 2 , An + –CD 3 , and An + –CD bond dissociation energies. The correlation of the measured An + –CD 2 bond energies with E p (6d 2 ) indicates that An + 6d orbitals are the dominant contributors in the An + –CD 2 bonds. Close examination of the relative reactivities of An + offers additional support that the balance of classical and modern views of molecular bonding may lie between Np + and Pu + and that the increased reactivity of Th + –Np + may be attributed to the increased spatial extension of the 5f orbitals whereas covalent An + bond formation may be more driven by the decreasing energies of the 5f orbitals across the actinide series.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Chemical and Structural Insights into Solid Electrolyte Interphase Evolution for Sodium Metal Electrodes

The solid electrolyte interphase (SEI) critically governs the reversibility of sodium metal batteries, through dynamically mediating ion transport and interfacial reactions. However, its kinetic evolution under operating conditions, and how it influences interfacial stability, remains poorly understood. Here, in this study, we reveal that the SEI undergoes coupled chemical and mechanical changes during sodium plating and stripping, leading to spatial and temporal heterogeneity that drives interfacial degradation. Synchrotron operando grazing-incidence wide-angle X-ray scattering and soft X-ray absorption spectroscopy capture the sequential formation and dissolution of inorganic SEI phases (NaF, NaH, NaOH, Na 2 PO 3 F), accompanied by depth-dependent alterations in organic SEI components. Mesoscale modeling connects this evolving SEI heterogeneity to localized current density fluctuations and stress accumulation at the Na interface, identifying pathways to electrically isolated sodium formation. These findings show that SEI instability fundamentally limits reversibility in sodium metal batteries, and that controlling SEI chemistry–mechanics coupling is essential to achieving its durability.

36 MATERIALS SCIENCE↗

The trailing edge is everywhere: tree rings reveal the transient risk of extinction hidden inside climate envelope forecasts

Given the importance of climate in shaping species’ geographic distributions, climate change poses an existential threat to biodiversity. Climate envelope modeling, the predominant approach used to quantify this threat, presumes that individuals in populations respond to climate variability and change according to species-level responses inferred from spatial occurrence data—such that individuals at the cool edge of a species’ distribution should benefit from warming (the “leading edge”), whereas individuals at the warm edge should suffer (the “trailing edge”). Using 1,558 tree-ring time series of an aridland pine (Pinus edulis) collected at 977 locations across the species’ distribution, we found that trees everywhere grow less in warmer-than-average and drier-than-average years. Ubiquitous negative temperature sensitivity indicates that individuals across the entire distribution should suffer with warming—the entire distribution is a trailing edge. Species-level responses to spatial climate variation are opposite in sign to individual-scale responses to time-varying climate for approximately half the species’ distribution with respect to temperature and the majority of the species’ distribution with respect to precipitation. These findings, added to evidence from the literature for scale-dependent climate responses in hundreds of species, suggest that correlative, equilibrium-based range forecasts may fail to accurately represent how individuals in populations will be impacted by changing climate. A scale-dependent view of the impact of climate change on biodiversity highlights the transient risk of extinction hidden inside climate envelope forecasts and the importance of evolution in rescuing species from extinction whenever local climate variability and change exceeds individual-scale climate tolerances.

54 ENVIRONMENTAL SCIENCES↗

Stress-Drop and Source Scaling of the 2019 Ridgecrest, California, Earthquake Sequence

Stress drop, while difficult to measure reliably and at scale, is a key source parameter for understanding the earthquake rupture process and its relationship to strong ground motion. Here, we use a P-wave spectral decomposition approach, designed for large and densely sampled datasets, to measure earthquake stress drop in the region surrounding the 2019 Ridgecrest, California, earthquake sequence. With more than 11,000 measurements of earthquake stress drop in the 20-yr time period from 2000 through 2019, this dataset provides an opportunity to understand how coseismic stress changes and how other geophysical factors relate to the distribution of stress drop and its evolution in space and time. We observe a mild but persistent deviation from self-similar scaling, with larger events having systematically higher stress drops, though this trend depends on the assumption of an omega-square source spectral model. Earthquake stress drop increases with hypocentral depth in this study region, and the Ridgecrest aftershocks tend to have higher stress drops than the pre-event seismicity. This is in part due to their deeper hypocenters. Coherent spatial patterns of stress drop in the aftershock sequence correlate with the slip distribution of the M 7.1 mainshock, whose northwest rupture tip terminated in a long-lived zone of enervated stress drop. Although physical interpretation of these results is complicated by the trade-offs between the timing, depth, and location of these earthquakes, the observations provide new insight into the physics of the earthquake source in an area of renewed seismic activity in southern California.

58 GEOSCIENCES↗

Tree rings reveal the transient risk of extinction hidden inside climate envelope forecasts

Given the importance of climate in shaping species’ geographic distributions, climate change poses an existential threat to biodiversity. Climate envelope modeling, the predominant approach used to quantify this threat, presumes that individuals in populations respond to climate variability and change according to species-level responses inferred from spatial occurrence data—such that individuals at the cool edge of a species’ distribution should benefit from warming (the “leading edge”), whereas individuals at the warm edge should suffer (the “trailing edge”). Using 1,558 tree-ring time series of an aridland pine (Pinus edulis) collected at 977 locations across the species’ distribution, we found that trees everywhere grow less in warmer-than-average and drier-than-average years. Ubiquitous negative temperature sensitivity indicates that individuals across the entire distribution should suffer with warming—the entire distribution is a trailing edge. Species-level responses to spatial climate variation are opposite in sign to individual-scale responses to time-varying climate for approximately half the species’ distribution with respect to temperature and the majority of the species’ distribution with respect to precipitation. These findings, added to evidence from the literature for scale-dependent climate responses in hundreds of species, suggest that correlative, equilibrium-based range forecasts may fail to accurately represent how individuals in populations will be impacted by changing climate. A scale-dependent view of the impact of climate change on biodiversity highlights the transient risk of extinction hidden inside climate envelope forecasts and the importance of evolution in rescuing species from extinction whenever local climate variability and change exceeds individual-scale climate tolerances.

54 ENVIRONMENTAL SCIENCES↗

Modeling detonation with CartaBlanca simulations

The accurate modeling of high explosive (HE) detonation and the con- sequent large solid deformation, failure, plastic flow, porosity growth, and shock wave propagation is important because simulations can capture spatial and temporal features that experimental diagnostics cannot capture. However, the simulation of the explosive event poses challenges to a computational scientist. These include the accurate modeling of ductile damage, crack formation, plastic deformation, as well as physical and nu- merical instabilities. The material response can be history-dependent and subject to large material deformation. Our research simulates the impact of a high explosive (Detasheet) onto a tantulum metal plate. We performed the simulations using CartaBlanca at different mesh resolutions. We decided the study would be impactful if we perform the simulations with the Material Point Method. Differences were observed at the different mesh resolutions in velocity and nodal stress magnitude, so increased mesh resolutions may be required. In addition, the Discontinuous Galerkin method would be needed to account for the large gas expansion.

97 MATHEMATICS AND COMPUTING↗

Search for Q ~0 Order near a Forbidden Bragg Position in Bi 2.1 Sr 1.9 CaCu 2 O 8+x with Resonant Soft X-ray Scattering

Identifying what broken symmetries are present in the cuprates has become a major area of research. Many authors have reported evidence for so-called “Q ~ 0” order that involves broken inversion, mirror, chiral, or time-reversal symmetry that is uniform in space. Not all these observations are well understood and new experimental probes are needed. Here we use resonant soft x-ray scattering (RSXS) to search for Q ~ 0 order in Bi 2.1 Sr 1.9 CaCu 2 O 8+x (Bi-2212) by measuring the region of a forbidden Bragg peak, (0,0,3), which is normally extinguished by symmetry but may become allowed on resonance if valence band order is present. Using circularly polarized light, we found that this reflection becomes allowed on the CuL 3 resonance for temperatures T c < T < T*, though remains absent in linear polarization and at other temperatures. This observation suggests the existence of spatially uniform valence band order near the pseudogap temperature. In addition, we observed periodic oscillations in the specular reflectivity from the sample surface that resemble thin film interference fringes, though no known film is present. These fringes are highly resonant, appear in all polarizations, and exhibit a period that depends on the location where the beam strikes the sample surface. We speculate that these fringes arise from interaction between some intrinsic valence band instability and extrinsic structural surface morphologies of the material. Our study supports the existence of some kind of Q ~ 0 broken symmetry state in Bi-2212 at intermediate temperatures. Further studies using a microfocused beam, which could disentangle microscopic effects from macroscopic heterogeneities, are needed to ascertain the universality of the effect.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗