Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “mean current vector”

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 73 records · Page 4

Field-aligned currents in the Jovian magnetosphere during the Ulysses flyby

The data recorded by the dual vector helium/fluxgate magnetometer flown onboard the Ulysses spacecraft during the flyby of Jupiter in February 1992 are analyzed with the aim of identifying the presence of field-aligned current signatures. Field-aligned current flow is expected wherever stress is being transmitted electromagnetically along the magnetic field direction. Sources of such currents at Jupiter are departures within the magnetosphere from corotation, momentum transfer from the solar wind, or centrifugally driven magnetospheric outflow. It is pointed out that the azimuthal field component provides a simple first-order means of monitoring the presence of currents, the currents occurring in regions where the azimuthal component changes significantly. The data from both inbound and outbound passes show evidence of 'leading' and 'lagging' azimuthal field signatures where the field bends out of the meridian, and which are signatures symptomatic of current systems associated with departures from corotation. On the outbound pass, the most intense signatures are found where the field switches from a configuration symptomatic of the field lagging corotation to a configuration representing the field leading. The latter configuration also corresponds to a tail-like displacement of the field and, indeed, the magnetometer data alone cannot distinguish the source of the current system, which could be due to solar wind magnetosphere coupling or which may arise from internal stress imbalance.

Dougherty, M. K.↗

Developing New Beam Constraints to Improve Sterile Neutrino Searches at NOvA

NOvA is a long-baseline neutrino oscillation experiment utilizing two functionally identical liquid scintillator detectors located 14 mrad off-axis from the NuMI beam at Fermilab. The position of the Near Detector (ND), 1 km from the neutrino production target, means it sees a range of off-axis angles (12.6-17.2 mrad), resulting in a noticeable energy gradient across the detector face. In addition, the neutrino production points vary widely along the decay pipe, which leads to a distribution of true neutrino baselines (L), and limits the sensitivity of our current sterile neutrino searches. This work presents a new "baseline reconstruction" technique designed to exploit our off-axis position to triangulate the precise origin of neutrinos within the NuMI beamline. By projecting the total momentum vector of Charged Current (CC) and Neutral Current (NC) interactions back to the central beam axis, we estimate the individual baseline for each event. We describe the technique, evaluate its performance, and demonstrate how this can improve the sensitivity of future searches for sterile neutrinos through reduced flux systematic uncertainties and improved L/E resolution.

Burns, Jessica [Cincinnati U.]↗

The use of the three dimensional coherency matrix in analyzing the polarization properties of plane waves

A technique for analyzing the polarization properties of plane waves is developed which offers a number of advantages over methods currently used in the analysis of both ground and satellite observations of waves. This technique reduces the computations required to find the wave normal vector, is less sensitive to common noise sources, and is amenable to analog implementation. This technique is applied specifically to the analysis of a proton whistler, but may also be used in most studies of ULF, ELF, and VLF magnetic wave phenomena.

Means, J. D.↗

Field Encapsulation Library The FEL 2.2 User Guide

This document describes version 2.2 of the Field Encapsulation Library (FEL), a library of mesh and field classes. FEL is a library for programmers - it is a "building block" enabling the rapid development of applications by a user. Since FEL is a library intended for code development, it is essential that enough technical detail be provided so that one can make full use of the code. Providing such detail requires some assumptions with respect to the reader's familiarity with the library implementation language, C++, particularly C++ with templates. We have done our best to make the explanations accessible to those who may not be completely C++ literate. Nevertheless, familiarity with the language will certainly help one's understanding of how and why things work the way they do. One consolation is that the level of understanding essential for using the library is significantly less than the level that one should have in order to modify or extend the library. One more remark on C++ templates: Templates have been a source of both joy and frustration for us. The frustration stems from the lack of mature or complete implementations that one has to work with. Template problems rear their ugly head particularly when porting. When porting C code, successfully compiling to a set of object files typically means that one is almost done. With templated C++ and the current state of the compilers and linkers, generating the object files is often only the beginning of the fun. On the other hand, templates are quite powerful. Used judiciously, templates enable more succinct designs and more efficient code. Templates also help with code maintenance. Designers can avoid creating objects that are the same in many respects, but not exactly the same. For example, FEL fields are templated by node type, thus the code for scalar fields and vector fields is shared. Furthermore, node type templating allows the library user to instantiate fields with data types not provided by the FEL authors. This type of flexibility would be difficult to offer without the support of the language. For users who may be having template-related problems, we offer the consolation that support for C++ templates is destined to improve with time. Efforts such as the Standard Template Library (STL) will inevitably drive vendors to provide more thorough, optimized tools for template code development. Furthermore, the benefits will become harder to resist for those who currently subscribe to the least-common-denominator "code it all in C" strategy. May FEL bring you both increased productivity and aesthetic satisfaction.

Moran, Patrick J.↗

Nucleon Isovector Axial Form Factors

We present results for the isovector axial vector form factors obtained using thirteen 2 + 1 + 1-flavor highly improved staggered quark (HISQ) ensembles generated by the MILC collaboration. The calculation of nucleon two- and three-point correlation functions has been done using Wilson-clover fermions. In the analysis of these data, we quantify the sensitivity of the results to strategies used for removing excited state contamination and invoke the partially conserved axial current relation between the form factors to choose between them. Our data driven analysis includes removing contributions from multihadron $Nπ$ states that make significant contributions. Our final results are $g_A$ = 1.292(53) stat (24) sys for the axial charge; $g_S$ = 1.085(50) stat (103) sys and $g_T$ = 0.991(21) stat (10) sys for the scalar and tensor charges; $\langle{r^2_A}\rangle$ = 0.439(56) sta t(34) sys fm 2 for the mean squared axial charge radius, $g^*_P$ = 9.03(47) stat (42) sys for the induced pseudoscalar charge; and $g_{πNN}$ = 14.14(81) stat (85) sys for the pion-nucleon coupling. We also provide a parametrization of the axial form factor $G_A(Q^2)$ over the range 0 ≤ $Q^2$ ≤ 1 GeV 2 for use in phenomenology and a comparison with other lattice determinations. We find that the various lattice data agree within 10% but are significantly different from the extraction of $G_A(Q^2)$ from the $ν$-deuterium scattering data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Observed currents on the earth's high-latitude magnetopause

A survey of electrical currents of the earth's magnetosphere, principally at high latitudes, as inferred from magnetic vector data acquired by the Hawkeye 1 satellite, is reported. A total of 536 candidate crossings of the magnetopause were examined. A reduced data set of 139 selected cases was analyzed in detail though solar wind dynamic pressure data were available for only 117 of these cases. Inferred values of the lineal current densities on the magnetopause are in the range 5.5 to 157.5 mA/m over a wide range of solar wind dynamic pressure from 1.17 to 16.1 nPa. The apparent normal thickness of the magnetopause current sheet ranges from 30 to 850 km with mean and median values of 185 and 158 km, respectively. It is argued that the radial rate of motion of the magnetopause is of the order of 2 km/s and hence that its true thickness is of similar magnitude. The relationship of these results to models of the geomagnetic field and to other related work is discussed.

Van Allen, J. A.↗

Nucleon isovector axial form factors

We present results for the isovector axial vector form factors obtained using thirteen 2 + 1 + 1 -flavor highly improved staggered quark (HISQ) ensembles generated by the MILC collaboration. The calculation of nucleon two- and three-point correlation functions has been done using Wilson-clover fermions. In the analysis of these data, we quantify the sensitivity of the results to strategies used for removing excited state contamination and invoke the partially conserved axial current relation between the form factors to choose between them. Our data driven analysis includes removing contributions from multihadron N π states that make significant contributions. Our final results are g A = 1.292 ( 53 ) stat ( 24 ) sys for the axial charge; g S = 1.085 ( 50 ) stat ( 103 ) sys and g T = 0.991 ( 21 ) stat ( 10 ) sys for the scalar and tensor charges; ⟨ r A 2 ⟩ = 0.439 ( 56 ) stat ( 34 ) sys fm 2 for the mean squared axial charge radius, g P * = 9.03 ( 47 ) stat ( 42 ) sys for the induced pseudoscalar charge; and g π N N = 14.14 ( 81 ) stat ( 85 ) sys for the pion-nucleon coupling. We also provide a parametrization of the axial form factor G A ( Q 2 ) over the range 0 ≤ Q 2 ≤ 1 GeV 2 for use in phenomenology and a comparison with other lattice determinations. We find that the various lattice data agree within 10% but are significantly different from the extraction of G A ( Q 2 ) from the ν -deuterium scattering data. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Post-hazard condition assessment of nuclear piping-equipment systems: Novel approach to feature extraction and deep learning

Over the past decade, the use of artificial intelligence techniques in the field of health-monitoring has gained significant interest, especially for structures such as building and bridges. However, applications to industrial systems such as equipment-piping systems in nuclear plants have not been explored. Here, in this paper, it is shown that the existing techniques developed for buildings and bridges cannot be extended directly to equipment-piping systems as the response of such systems is governed by multiple localized modes unlike that in buildings and bridges. This paper proposes a new approach that consists of three key aspects: (i) a novel vector of degradation-sensitive features extracted from measured data, (ii) using a deep Artificial Neural Network (ANN) for diagnosis of degradation location and degradation severity, and (iii) consideration of uncertainty in degradation severity when training the ANN. Degradation in piping-equipment systems can occur due to flow-accelerated erosion and corrosion. These locations can potentially exhibit damage such as localized yielding or initiation of cracking due to an external event such as an earthquake. Moreover, such locations can at times go undetected by current inspection techniques. Therefore, a robust framework is needed for detection of degradation after a seismic event. This manuscript proposes a proof-of-concept framework, which utilizes data collected from sensors to generate a deep ANN database for predicting degraded locations and severity in a piping-equipment system. Degradation severity is classified as minor, moderate, and severe. In the suggested methodology, a novel vector of degradation-sensitive features is extracted from the sensor data to train the ANN. A simple piping-equipment system is selected to demonstrate feature extraction as a means to simplify pattern recognition, explore the design and parameters of an ANN, and develop a sensor placement strategy. The effectiveness of the proposed framework is demonstrated on a realistic primary safety system of a two-loop nuclear reactor. It is shown that the proposed post-hazard condition assessment framework is able to detect degraded locations along with the severity levels, including minor degradation, with considerably higher accuracy.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Machine Learning with Gradient-Based Optimization of Nuclear Waste Vitrification with Uncertainties and Constraints

Gekko is an optimization suite in Python that solves optimization problems involving mixed-integer, nonlinear, and differential equations. The purpose of this study is to integrate common Machine Learning (ML) algorithms such as Gaussian Process Regression (GPR), support vector regression (SVR), and artificial neural network (ANN) models into Gekko to solve data based optimization problems. Uncertainty quantification (UQ) is used alongside ML for better decision making. These methods include ensemble methods, model-specific methods, conformal predictions, and the delta method. An optimization problem involving nuclear waste vitrification is presented to demonstrate the benefit of ML in this field. ML models are compared against the current partial quadratic mixture (PQM) model in an optimization problem in Gekko. GPR with conformal uncertainty was chosen as the best substitute model as it had a lower mean squared error of 0.0025 compared to 0.018 and more confidently predicted a higher waste loading of 37.5 wt% compared to 34 wt%. The example problem shows that these tools can be used in similar industry settings where easier use and better performance is needed over classical approaches. Future works with these tools include expanding them with other regression models and UQ methods, and exploration into other optimization problems or dynamic control.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Subphotospheric current systems and flares

Subphotospheric current systems inferred from recent vector magnetograph observations imply the existence of electric currents penetrating the photosphere and thus flowing deep in the solar convection zone. These currents presumably originate in an internal dynamo that supplies the observed photospheric magnetic fields through the buoyant motions of the initially deeply-buried flux tubes. The coronal fields resulting from this process therefore must carry slowly-varying currents driven by emfs remote from the surface. These currents may then drive solar-flare energy release. This paper discusses the consequences of such a deep origin of the coronal parallel currents. Simple estimates for a large active region suggest a mean current-closure depth of not less than 10,000 km, with a subphotospheric inductance of not less than 100 H and a subphotospheric stored energy of not less than 10 to the 33rd ergs.

Hudson, Hugh S.↗

Effects of Nose Bluntness on Hypersonic Boundary-Layer Receptivity and Stability Over Cones

The receptivity to freestream acoustic disturbances and the stability properties of hypersonic boundary layers are numerically investigated for boundary-layer flows over a 5 straight cone at a freestream Mach number of 6.0. To compute the shock and the interaction of the shock with the instability waves, the Navier-Stokes equations in axisymmetric coordinates were solved. In the governing equations, inviscid and viscous flux vectors are discretized using a fifth-order accurate weighted-essentially-non-oscillatory scheme. A third-order accurate total-variation-diminishing Runge-Kutta scheme is employed for time integration. After the mean flow field is computed, disturbances are introduced at the upstream end of the computational domain. The appearance of instability waves near the nose region and the receptivity of the boundary layer with respect to slow mode acoustic waves are investigated. Computations confirm the stabilizing effect of nose bluntness and the role of the entropy layer in the delay of boundary-layer transition. The current solutions, compared with experimental observations and other computational results, exhibit good agreement.

Kara, Kursat↗

Mass Inferencing Model Creation and Deployment to the RASSOR Lunar Excavation Robot

The Regolith Advanced Surface Systems Operations Robot (RASSOR) Excavator is a teleoperated mobile robotic platform with a unique space regolith excavation capability. The Intelligent Capabilities Enhanced RASSOR research project developed functionality for inferencing regolith mass ingested during RASSOR operation, enhancing RASSOR’s ability to successfully complete ISRU missions. To teleoperate or run autonomously, it is crucial for the quantity of regolith mass ingested by RASSOR to be available as a system state for efficient operation. For example, during autonomous operation, RASSOR should navigate and move to a processing plant to offload the collected regolith when the drums are full; without knowledge of how much mass is in the drums, this type of high-level planning is not possible. Four distinct modeling approaches were employed in developing a mass inferencing approach that could work on RASSOR. All take in system states, such as arm/drum positions, velocities, currents, voltages, and robot pose, and output a mass prediction for each set of the robot’s bucket drums.1) A neural network model that takes a vector of normalized system states; 2) A model that uses the integrated power consumption of an arm-raise (normalized by velocity); 3) A model that uses average drum current over a variable length interval of the drum disengaged from the surface; and 4) A real-time estimation model that aggregates excavation drum current. The developed models run in real time, outputting predictions for the front and rear drums, timestamp of the last prediction, and total mass in RASSOR’s drums. Further testing is required to validate the arm-raise model (2), though initial tests indicate reasonable performance (<10% mean error) on the hardware. The linear fit of average drum-current model (3) had a front value of r^2=0.99 and a rear value of r^2=0.98 on the validation dataset. This model currently has the best performance on unseen data. The real time model (4) is still in development, though initial results on a small subset of the training data show that it has high accuracy in predicting the increase in mass during excavation. Though work remains to be done with deploying a high-fidelity model to the physical system that makes predictions with error below the desired threshold, the modular architecture for model development allows quick adjustment of parameters to increase model fidelity. This architecture can also be adapted to use lunar excavation data to create models that are reflective of RASSOR’s dynamics when operating on the lunar surface. The results are promising as it has been shown that models can be developed that accurately estimate excavated regolith mass.

rassor↗

New models of Saturn's magnetic field using Pioneer 11 Vector Helium Magnetometer data

In a reanalysis of the Vector Helium Magnetometer data taken by Pioneer 11 during its Saturn encounter in 1979, using improvements in the data set and in the procedures, studies are made of a variety of models. The best is the P(11)84 model, an axisymmetric spherical harmonic model of Saturn's magnetic field within 8 Saturn radii of the planet. The appropriately weighted root mean square average of the difference between the observed and the modeled field is 1.13 percent. For the Voyager-based Z3 model of Connerney, Acuna, and Ness, this average difference from the Pioneer 11 data is 1.81 percent. The external source currents in the magnetopause, tail, bow shock, and perhaps ring currents vary with time and can only be crudely modeled. An algebraic formula is derived for calculating the L shells on which energetic charged particles drift in axisymmetric fields.

Davis, L., Jr.↗

S-band Smallsat InSAR Constellation for Surface Deformation Science

Surface deformation studies using repeat-pass interferometric SAR have evolved into a powerful tool for geophysicists studying earthquake fault zones, volcanoes, ice sheet motion, and subterranean aquifers. Longer wavelengths (S-Band and L-Band) are preferred because they do not decorrelate as quickly as shorter wavelengths. Rapid revisit (1-3 days) is preferred because it allows the study of these phenomena at the timescales at which they commonly occur. Global access on such timescales is also required. Vector surface deformation measurements, taken from more than one direction, are a desired feature. This paper describes the conceptual architecture of a longer wave length, Smallsat SAR constellation of up to 12 satellites for rapid revisit surface deformation studies. The key to making such a constellation affordable is to lower launch costs, spacecraft costs, and instrument (SAR) costs. The first two objectives can be achieved using an ESPA-ring class, or Smallsat spacecraft. The third objective requires a SAR instrument sized to fit the mass and volume constraints imposed by such a spacecraft. Current state-of-the-art in miniaturization of electronics means that the radar transmit, receive and data handling functions can easily be implemented in a compact, low mass solution. The most significant challenge in designing a SAR to fit the Smallsat paradigm is in the dimensions of the antenna. The antenna sizing problem is addressed by adopting a smaller antenna than allowed by conventional SAR design rules. The baseline antenna design is simple, requiring no electronic beam-steering or beam-forming capability. Both reflectarray and microstrip patch antenna solutions are considered. The antenna structure is dual-purpose, to limit the overall system mass, with solar panels on the backplane providing power for the radar and spacecraft. The proposed solution easily accommodates radar squint angles of +/-30 degrees for repeat-pass interferometry measurements from multiple direct

Freeman, Anthony↗

Data-Driven Transition Path Analysis Yields a Statistical Understanding of Sudden Stratospheric Warming Events in an Idealized Model

Abstract Atmospheric regime transitions are highly impactful as drivers of extreme weather events, but pose two formidable modeling challenges: predicting the next event (weather forecasting) and characterizing the statistics of events of a given severity (the risk climatology). Each event has a different duration and spatial structure, making it hard to define an objective “average event.” We argue here that transition path theory (TPT), a stochastic process framework, is an appropriate tool for the task. We demonstrate TPT’s capacities on a wave–mean flow model of sudden stratospheric warmings (SSWs) developed by Holton and Mass, which is idealized enough for transparent TPT analysis but complex enough to demonstrate computational scalability. Whereas a recent article (Finkel et al. 2021) studied near-term SSW predictability, the present article uses TPT to link predictability to long-term SSW frequency. This requires not only forecasting forward in time from an initial condition, but also backward in time to assess the probability of the initial conditions themselves. TPT enables one to condition the dynamics on the regime transition occurring, and thus visualize its physical drivers with a vector field called the reactive current . The reactive current shows that before an SSW, dissipation and stochastic forcing drive a slow decay of vortex strength at lower altitudes. The response of upper-level winds is late and sudden, occurring only after the transition is almost complete from a probabilistic point of view. This case study demonstrates that TPT quantities, visualized in a space of physically meaningful variables, can help one understand the dynamics of regime transitions.

Meteorology & Atmospheric Sciences↗

Lunar Reconnaissance Orbiter (LRO) Guidance, Navigation and Control (GN&C) Overview

The National Aeronautics and Space Administration s (NASA) Lunar Reconnaissance Orbiter (LRO) launched on June 18, 2009 from the Cape Canaveral Air Force Station aboard an Atlas V launch vehicle and into a direct insertion trajectory to the oon. LRO, which was designed, built, and operated by the NASA Goddard Space Flight Center in Greenbelt, MD, is gathering crucial data on the lunar environment that will help astronauts prepare for long-duration lunar expeditions. The mission has a nominal life of 1 year as its seven instruments find safe landing sites, locate potential resources, characterize the radiation environment, and test new technology. To date, LRO has been operating well within the bounds of its requirements and has been collecting excellent science data images taken from the LRO Camera Narrow Angle Camera of the Apollo landing sites appeared on cable news networks. A significant amount of information on LRO s science instruments is provided at the LRO mission webpage. LRO s Guidance, Navigation and Control (GN&C) subsystem is made up of an onboard attitude control system (ACS) and a hardware suite of sensors and actuators. The LRO onboard ACS is a collection of algorithms based on high level and derived requirements, and reflect the science and operational events throughout the mission lifetime. The primary control mode is the Observing mode, which maintains the lunar pointing orientation and any offset pointing from this baseline. It is within this mode that all science instrument calibrations, slews and science data is collected. Because of a high accuracy requirement for knowledge and pointing, the Observing mode makes use of star tracker (ST) measurement data to determine an instantaneous attitude pointing. But even the star trackers alone do not meet the tight requirements, so a six-state Kalman Filter is employed to improve the noisy measurement data. The Observing mode obtains its rate information from an inertial reference unit (IRU) and in the event of an IRU failure, the rate data is be derived from the star tracker, but with degraded pointing performance. The Delta-V control mode responsibility is to maintain attitude pointing during the cruise trajectory, insertion burns and lunar orbit maintenance by adjustments made to the spacecraft s velocity magnitude and vector direction. The ACS also provides for a thruster based system momentum management algorithm (known as Delta-H) to maintain the system and wheel momentum to within acceptable levels. In the event an anomaly causes the LRO spacecraft to lose the ability to maintain its current attitude pointing, a Sun Safe mode is included in the ACS for the purpose of providing a known power and thermally safe coarse inertial sun attitude for an indefinite period of time, within the manageable limits of the reaction wheels. The Sun Safe mode is also the initial spacecraft control mode off of the launch vehicle and provides for a means to null tip-off rates immediately after separation. The nominal configuration is to use the IRU for rate information in the controller. In the event of a gyro failure a gyroless control mode was developed that computes rate information from the CSS data.

Garrick, Joseph↗

Algorithmic Classification of Raman Spectra Biosignatures: Improving Life Detection Confidence

“Agnostic” biosignatures – indicators of life (or the absence of life), independent of a particular biochemistry – are increasingly considered a high standard for life detection. The Ladder of Life Detection (2018) called for investigating how combinations of independent and different potential biosignatures affect confidence. To address this gap, statistical classification of elemental abundances, isotopic fractionation, and reflectance spectroscopy (VNIR) has been implemented. Raman spectroscopy, highly desirable due to its wide availability, has the potential to improve this predictive power. This work implemented biosignature classification algorithms on Raman data alone, in preparation for combination with the other data types. Raman spectroscopy data was collected from published databases and papers as part of a manually curated dataset of “indicative” and “non-indicative of life” samples. These currently include 61 non-indicative samples (meteorites, magnetite); 3 indicative living samples (bacteria); 20 indicative non-living samples (chalk, bone); and 12 indicative mixed (with non-indicative material) samples (soil, microbial mats). Laboratory work is ongoing to characterize additional samples, particularly a greater breadth of mixed systems. Spectra were interpolated, filtered with the Savitzsky-Golay filter, and de-noised. For a preliminary examination, agnostic features were manually extracted including mean intensity, number of peaks, and mean peak width. Different peak prominences and filtering polynomials were used to refine features. Classification algorithms were implemented: k-nearest neighbors (KNN), logistic regression (LR), linear support vector machines (SVM), random forest (RF), Gaussian naïve bayes (GNB). Lastly, Monte Carlo simulations on 1,000 50%-train-test-splits were used to validate classification performance and feature significance. The preliminary feature set achieved its highest AUC of 0.52 with LR, with no strongly discriminatory features. Work to improve feature extraction, such as through deep learning with back propagation, is planned. In future work, the Raman data will be combined with the other data types, and potentially new data types such as enantiomeric excess. This project was partially supported through the NASA Ames Project EXcellence (APEX) incubator program.

Astrobiology↗

Light-driven nanoscale vectorial currents

Controlled charge flows are fundamental to many areas of science and technology, serving as carriers of energy and information, as probes of material properties and dynamics and as a means of revealing or even inducing broken symmetries. Emerging methods for light-based current control offer particularly promising routes beyond the speed and adaptability limitations of conventional voltage-driven systems. However, optical generation and manipulation of currents at nanometre spatial scales remains a basic challenge and a crucial step towards scalable optoelectronic systems for microelectronics and information science. Here we introduce vectorial optoelectronic metasurfaces in which ultrafast light pulses induce local directional charge flows around symmetry-broken plasmonic nanostructures, with tunable responses and arbitrary patterning down to subdiffractive nanometre scales. Local symmetries and vectorial currents are revealed by polarization-dependent and wavelength-sensitive electrical readout and terahertz (THz) emission, whereas spatially tailored global currents are demonstrated in the direct generation of elusive broadband THz vector beams 17 . We show that, in graphene, a detailed interplay between electrodynamic, thermodynamic and hydrodynamic degrees of freedom gives rise to rapidly evolving nanoscale driving forces and charge flows under the extremely spatially and temporally localized excitation. These results set the stage for versatile patterning and optical control over nanoscale currents in materials diagnostics, THz spectroscopies, nanomagnetism and ultrafast information processing.

36 MATERIALS SCIENCE↗